
книги из ГПНТБ / Кудрин, Л. П. Статистическая физика плазмы
.pdf2 |
|
|
lim J /(r) - g C«.k<Pn, k dr = |
0, |
(28.9) |
где |
|
|
'/( r ) 4>;ik (r) dr. |
|
|
Пусть F ( r) — произвольная непрерывная |
функция, |
опреде |
ленная в объеме V. Обозначим С„,к коэффициенты Фурье этой функции по системе (28.6). Тогда имеет место очевидное ра венство
Сп , к = |
VJ ■F (г) Фп, к (r) dr |
=Vj F (г) Ф„, к (r) dr |
= |
к* (28.10) |
|
|
|||||
С учетом выражений (28.8) |
и (28.9) |
|
|
||
Нт Г|/г(г)— |
V |
с п, к Фп и (г) |2 cfr - |
|
|
|
|
V |
п.к |
|
|
|
= y j lini П ^ |
— yjCn.k4P„ к (г) |2 dr = |
0. |
|
||
« |
vn |
1г |
|
|
|
Этим доказана полнота системы функций {cfn,k} в объеме V. Таким образом, система функций (28.5) образует ортонормироваипый базис в пространстве состояний одной частицы и мо жет быть использована для изучения объемных свойств макро скопической системы наряду с обычным одпочастичным бази сом в виде плоских волн.
До сих пор не использовался явный вид функций (28.5) внутри ячейки. Поэтому описанная процедура (с небольшими изменениями) пригодна для построения ячеистого базиса на произвольной полной ортонормированпой системе функций, определенных внутри ячейки.
Вторичное квантование с помощью ячеистого базиса. В представлении вторичного квантования гамильтониан систе
мы фермионов, заключенных в объеме |
V, |
имеет вид |
||||
|
Я = J d r $+(!■) \^£- + U(r) |
ЧДг) |
||||
+ ~ |
f jd rd r'¥ + (r)f+ (r')« (r — r')V (r') Y (г). (28.11) |
|||||
Здесь TF (г)— полевой оператор, характеризующий уничтоже- |
||||||
ние электрона |
в точке г. |
Операторы |
Ч/+ |
/ч |
||
и XF подчиняются |
||||||
обычным соотношениям аптикоммутации: |
|
|
|
|||
(?(г), ¥+ (г')} = б (г — г'); |
{Ф(Г), |
Y (!-')} = |
||||
|
= {?+ |
(г) Y+ (г')} |
- |
о. |
|
( 2 8 . 1 2 ) |
29#
Введем операторы |
а+ к |
и о„,к. описывающие соответствен |
|
но возникновение и поглощение частицы в |
состоянии фп,к (г). |
||
Тогда |
|
|
|
¥ |
(Г) = |
2 а п . к Ф „ , к ( Г)- |
(28.13) |
|
|
п .к |
|
Операторы а и а+ подчиняются соотношениям антикоммутации:
|
|
{#п, к» @п', к'} |
==: ®п, п'^к , к'» |
|
|
|
|||
|
/N |
/N |
/N |
/S |
|
|
|
|
|
|
(“ n . k f l n ' . k ' l |
= № к ап'.к'}= 0 - |
|
|
(28.14) |
||||
В новом представлении получим вместо выражения |
(28.11) |
||||||||
-f 4 -Sn, |
|
|
п2 |
|
\ ^ |
I-'- |
|
|
|
|
|
2М |
’ U | к3 / а+ t ап . + |
|
|||||
n |
kt , к |
I |
- */ п П, k t п, к* 1 |
|
|||||
|
|
|
|
|
|
|
|||
n' ki, к2,Sк3, к4< k i k 2 l « l k 3 k 4 > n , n ' a + k i a + i k 2 a n. |
U i a n |
к>( |
|||||||
|
|
|
|
|
|
|
|
|
(28.15) |
где |
|
|
|
|
|
|
|
|
|
<ki |
£ - + и |
кз) п = I |
dr(p"’ к‘ (Г) |
2М |
+ г/(г) |
« V к, (г). |
|||
ж |
|
|
|
|
|||||
< k lf к2 | и | к3, к4>„, |
= | drdr' |
ki (г) ф^, kj (г') X |
X и (г, r')<Pn.ik4(r')«Pn, ki(r).
Гамильтониан (28.15) можно представить в виде суммы «адди тивного» гамильтониана и члена, учитывающего взаимодейст вие между ячейками:
|
|
|
Я = |
Яад + |
W, |
|
|
(28.16) |
где |
|
|
|
|
|
|
|
|
|
|
|
* = 2 ЯП; |
|
|
(28.17) |
||
Я„ = |
\ |
ki |
2М |
+ и |
к, / а + . |
а |
. 4- |
|
|
|
|
J/ n |
к, |
|
П, к, 1 |
||
|
к,, к, |
|
|
|
|
|
|
|
+ -2“ S |
< к 1( к2 |
[ и | к3> |
k 4 > n . n ' f l + k i a + k2a n k t a, к,’ |
|||||
ki • ka, ki» к« |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( 2 8 .1 8 ) |
II* 291
n / п' к ,, к., k3, к4
(28.19)
При этом номер ячейки п выступает в задаче как дополни тельное квантовое число. Выражения (28.16) —(28.19) являются точными в том смысле, что гамильтониан в формуле (28.16) в точности эквивалентен исходному гамильтониану (28.11), так что.введение ячеистого базиса само по себе не вносит никаких приближений и не связано ни с какими моделями. В частности,
неаддитивная часть взаимодействия W дает вклад в энергию, который в случае однородной системы пропорционален числу частиц. В то же время введение ячеистого базиса позволяет устранить трудности, возникающие при попытках распростра нить методы оценок, приведенные в § 27, на системы с макро скопически большим числом степеней свободы.
Собственные функции в представлении ячеистого базиса. Оценки энергии в аддитивном приближении. Физически очевид но, что при достаточно большом размере ячеек неаддитивный
член W гамильтониана (28.16) мал по сравнению с # ад и его можно учитывать по теории возмущений. Отметим, что ячей ка вовсе не должна содержать макроскопически большое число частиц, как этого требует обычная статистическая механика, которая полностью пренебрегает взаимодействием между «под системами». Поскольку для применимости теории возмущений
требуется лишь относительная малость W, то ограничения на размер ячейки могут оказаться значительно менее жесткими, чем этого требует статистическая механика. В частности, рас четы для предельного случая г8^>1 подтверждают это.
Энергия системы в нулевом приближении по W равна наинизшему собственному значению аддитивного гамильтониана
(28.20)
П
Поскольку Н„ с различными п коммутируют, то собствен
ный вектор гамильтониана Яад можно представить в виде пря мого произведения собственных векторов гамильтонианов яче-
П
где |Ч ^П> определяется из уравнения Шредингера, описываю щего движение частиц в n-й ячейке:
н п I Ч'п, k > = en>k I 4 V k > . |
( 2 8 .2 2 ) |
292
В этом уравнении п характеризует номер ячейки, а к — полный набор квантовых чисел, определяющих собственную функцию.
Собственные функции полного гамильтониана (28.16) можно представить в виде
I |
Y > |
= 2 C ( k lf k2, |
■ • |
К • |
• - ) Х |
|
|
|
|
14 |
|
|
|
|
|
|
х |
| 4 ^ ( 4 |
к2, . |
. ., к„ . . |
. > , |
(28.23) |
|
где С(кь к2, |
...)— функция |
макроскопически |
большого |
числа |
|||
аргументов (по числу ячеек) и, разумеется, |
не может быть вы |
||||||
числена точно, даже если все |xFn,k> |
известны. Однако, когда |
||||||
неаддитивность мала, С можно |
аппроксимировать по |
теории |
возмущений с желаемой точностью. Таким образом, задача сво дится к определению l^n.kX из уравнения (28.22).
Полная система собственных векторов ячейки содержит со стояния с произвольным числом частиц (и, в частности, состоя
ния с jV „ » 1). Однако очевидно, |
что разложение |
(28.23) для |
любого реального (равновесного) |
состояния |Ч/> |
содержит в |
основном векторы |Ч/'п,к>, соответствующие равновесному чис лу частиц в ячейке N„ = nV0= (N/V)V0. Разумеется, в известных случаях флуктуации числа частиц в ячейке AN„ могут быть не
малы по сравнению |
с N„, и тогда |
их следует |
явно |
учесть. |
||
Во |
всяком |
случае, |
нас интересуют |
только решения |
уравне |
|
ния |
(28.22) |
для |
относительно небольшого |
числа |
частиц |
Nn -\-AN п-
Отметим, что неопределенность числа частиц в ячейке свя зана, по-видимому, не только с термодинамическими флуктуа циями. Как это видно из введения ячеистого базиса во вторич ном квантовании, номер ячейки приобретает смысл нового квантового числа. Поэтому вряд ли можно говорить о строгой локализации определенного числа частиц в ячейке.
Таким образом, вычисление энергии электронного газа сво дится к решению уравнения Шредингера (28.22) внутри одной ячейки, содержащей конечное число частиц. Разумеется, эта задача все еще очень сложна и, вообще говоря, не может быть решена точно, поскольку, во-первых, внутри ячейки в общем случае нет малого параметра и, во-вторых, ячейка может со держать все еще большое (хотя и конечное) число взаимодейст вующих частиц. Поэтому речь может идти лишь о получении оценок сверху и снизу для точного собственного значения мето дами, изложенными в § 27.
Уравнение (28.22) описывает систему, заключенную в ячей ку объемом Vo и содержащую конечное число частиц. Для такой системы уже сравнительно легко построить несложные подпространства УЛ, обладающие достаточной полнотой по от ношению к интересующим нас состояниям |Чгп,к>- Поэтому для оценки значений Еп_к сверху и снизу можно воспользо*
293
—; |
ч |
ваться проекциями |
и НШ соответствующих операторов на |
эти подпространства. Верхние и нижние граничные значения для Е n.k являются собственными значениями уравнений Шредингера:
1 Уп> к> = |
£ ? 'к | Фп,к> |
(28.24) |
и |
|
|
| Vn. к> |
=_Е?к I Wn, к. |
(28.25) |
Тогда, согласно результатам § 27, |
|
|
£ ? ;< ^ п .к < £ ? к . |
(28.26) |
|
Отсюда следует оценка для энергии большой системы объе |
||
ма У в аддитивном приближении |
|
|
2 £ Г к < £ п Дк < 2 Ё ^ к. |
(28.27) |
|
п —' |
а |
|
Важно, что для основного состояния пространственно однород ной системы (или системы в периодическом внешнем поле) E f 0 и E f 0не зависят от п:
ё а0д = 2 ё^о = ЛГЯЧ~Ef0 = NE^o/Nn. |
(28.28) |
П |
|
Видно, что при оценках энергии основного состояния полная энергия основного состояния большой системы пропорциональ на полному числу частиц в системе N. Отметим, что собствен
ные значения уравнений |
Шредингера (28.24) и (28.25) |
зави |
сят от числа частиц в ячейке и не зависят от полного |
числа |
|
частиц в системе. |
|
|
Аналогично для оценки снизу получаем |
|
|
ё а0д s |
2 «So = A^m /AU |
(28.29) |
|
П |
|
Тогда из уравнений (28.28) и (28.29) следует оценка для энер
гии, приходящейся на одну частицу в аддитивном |
прибли |
жении, |
|
/Nn < £ „ < l f 0/JV„, |
(28.30) |
где Е0— точное значение энергии основного состояния системы. В заключение параграфа сделаем одно существенное за
мечание. |
Очевидно, что |
уравнения |
(28.24) и (28.25) значи |
тельно проще исходного |
уравнения |
(28.22). Так, если подпро |
|
странство |
отвечает возбуждениям одной-двух пар над неко |
||
торым исходным состоянием | 0 „>. |
[это может быть детерми |
нант из плоских волн в ячейке (28.5)], то уравнения (28.24) и
294
(28.25) описывают в сущности именно движение этих пар. В этом случае приходится решать уравнения Шредингера, за висящие от координат двух — четырех квазичастиц, хотя исход
ное состояние |
|Оп> |
соответствует |
N n частицам |
в |
ячейке, |
|
число |
которых |
может |
быть много |
большим единицы. |
Более |
|
того, |
и эти уравнения |
нужно решать с заданной |
точностью, |
определяемой характером задачи, а отнюдь не со спектроскопи ческой точностью, как это нужно было делать Лёвдину в случае гелиеподобных ионов.
§ 29. ВЕРХНЯЯ И НИЖНЯЯ ОЦЕНКИ ДЛЯ ЭНЕРГИИ |
|
||
ОСНОВНОГО СОСТОЯНИЯ б о л ь ш о й |
с и с т ем ы |
|
|
Оценка сверху. Пусть |
— эрмитов проекционный опера |
||
тор, проектирующий пространство |
функций Ж», |
определенных |
|
в п-й ячейке, па подпространство |
функций $?„, |
определенных |
|
в этой ячейке. Легко видеть, что |
|
|
|
0 = |
П0„ |
(29.1) |
|
|
П |
|
|
также является эрмитовым проекционным оператором:
|
|
О2= |
О, |
0+ |
- |
О. |
|
|
|
(29.2) |
|
Эго определяется гем фактом, что операторы 0„ |
с различны |
||||||||||
ми п коммутируют. |
|
проекцию |
оператора |
/ч |
(28.16) |
ра |
|||||
Определим |
внешнюю |
Я |
|||||||||
венством |
|
|
|
|
|
|
|
|
|
|
|
|
|
■ 'Ч /-Ч / N У-Ч |
/ \ / Ч |
|
У«Ч |
« Ч |
|
|
|
|
|
|
|
Я = ОНО = 0Яад0 + W. |
|
|
|
|
|||||
Тогда, согласно неравенству (27.7), |
|
|
|
|
|
||||||
|
|
|
Е0< |
£ 0, |
|
|
|
|
(29.3) |
||
где Ео — паипизшее |
собственное |
значение гамильтониана |
/\ |
||||||||
Я. |
|||||||||||
Можно показать, что для |
ячеек, |
содержащих |
относительно |
||||||||
большое число частиц, вклад неаддитивного члена W мал по |
|||||||||||
сравнению с 0 Я ад0. |
Тогда |
Ей представляется |
в |
виде |
ряда |
||||||
\ |
|
|
|
|
|
|
|
|
|
|
|
по W, аналогичного обычному ряду теории возмущений в по |
|||||||||||
левой форме |
(разложение |
по |
связным |
диаграммам), |
по |
с |
другими «свободными» функциями Грина, которые теперь с са мого начала учитывают различные взаимодействия внутри ячейки. Задача, таким образом, сводится к диагонализации
гамильтониана 0 Я ад0. Собственные функции этого гамильто-
295
пиана, так же как и гамильтониана Н, принадлежат подпрост ранству ЯЛ=П ЯЛП. Но в подпространстве ЯЛ
ОЯадО = У о н по = у ОпЯпО„ = V н„. |
(29.4) |
Следовательно, гамильтониан ОЯадО можно заменить аддитив-
/ N
ным гамильтонианом 2 Я „, диагонализация которого не пред-
П
ставляет большого труда.
Оценка снизу. В представлении ячеистого базиса гамиль
тониан |
системы |
имеет |
вид |
(28.16) — (28.17). Нетрудно |
видеть, |
||
что для |
оценок |
снизу |
достаточно уметь |
строить внутренние |
|||
|
/ч |
|
/ч |
|
|
|
|
проекции Яп операторов Я„. Действительно, поскольку |
|
||||||
то |
|
|
Яп < |
Яп, |
|
|
|
|
Н ^ ^ H n + W = H. |
|
(29.5) |
||||
|
|
|
|||||
|
|
|
П |
|
' |
|
|
Следовательно, |
согласно результатам § 27, |
наинизшие собствен- |
|||||
|
|
|
/ч |
/Ч |
|
|
|
ные значения операторов Я и Я связаны неравенством |
|
||||||
|
|
|
|
Д0 > |
E f . |
|
(29.6) |
Как и в случае построения верхней оценки для Еа, Е |
в прин- |
||||||
|
|
|
|
|
|
|
/ч |
цице можно представить в виде ряда теории возмущений по W, и задача сводится к диагонализации аддитивного гамильтониа
на ЕЯП (см. § 27).
П
Однако при построении внутренней проекции, соглас-
/Ч
но (27.11), необходимо явное выражение для Д|/2, что в слу
чае задачи многих тел эквивалентно диагонализации |
гамильто- |
/Ч |
трудность |
ниана Н„. Построить такой оператор нелегко. Эту |
в задаче многих тел можно обойти, построив нижнюю проек цию другого вида. Так, для операторов вида
А = р+ р |
(29.7) |
можно несколько иначе, чем в выражении |
(27.11), определить |
/Ч |
/Ч |
внутреннюю проекцию к на подпространство ЯЛ = ОШ,'- |
|
A =p+Ofr |
(29.8) |
Покажем, что и в этом случае имеет место неравенство |
|
О < А < А. |
(29.9) |
296
Пусть |
О — эрмитов |
проекционный оператор. |
Оператор А, |
||||
очевидно, |
определен в том |
же |
пространстве^, |
что |
и опергГ- |
||
/Ч |
|
неотрицателен, т. е. для |
любого |ф > |
||||
тор А. Если оператор А |
|||||||
|
|
|
|
/ч |
(29.8) |
также |
|
из пространства Ш <cpi Л |ф > ^ 0 , то оператор^ |
|||||||
неотрицателен. Действительно, |
|
|
|
|
|||
<Ф I А I Ф> = <Ф I Р+Ор | Ф> == <Ф I Р+ОаР | ф > |
= |
||||||
|
= < ф | р+О+Ор | ф > = <Орф | Орф> > |
0. |
|
|
|||
Докажем теперь, что |
|
|
|
|
|
|
|
|
|
Л < |
Л. |
|
|
(29.10) |
|
|
|
/Ч |
/Ч |
поскольку |
|
|
|
Для этого заметим, что 0 ^ 1 , |
|
|
|
||||
<ф |Т _ |
ОI ф> = <ф I (Т— о2) | Ф> = <Ф | о — 6+) (Г— 5) I ф> = |
||||||
|
=»<(? — О)ф |(Т — 6 ) ф> > 0 . |
|
|
(29.11) |
|||
Отсюда следует, что для произвольного |ф > из |
простран |
||||||
ства Ш |
|
|
|
|
|
|
|
|
<Ф | А | ф> = <ф | р+0р | ф> = <рф | О| рф> < |
|
|
||||
|
/ \ |
|
/ ч |
/N А |
|
|
|
|
< <Р+ф| 1I РФ) = <Ф|Р+ Р|Ф>. |
|
|
|
что и доказывает неравенство (29.10).— Запишем теперь гамильтониан Н„ в представлении «ячеи
стых» плоских волн: |
|
|
, сч | |
v ^ '/!6s.0exp(ikr) |
внутри n-й ячейки; 0Q |
фп.к.о (Г, |
0 |
вне n-й ячейки. |
| |
При этом, если объем ячейки Vo достаточно велик, можно по ложить*
<кх, к21и | к3, к4> = и ( | ki — к31) б(к4 + к2 — к3 — к4),
где и — оператор взаимодействия электронов. Тогда
|
Нп = 2 |
<ki I Т + U | k3V & kfl„.ki - |
M V 9) 2 « (Ч) + |
|
|
|
|
|
4 |
|
|
+ (l/2V0)2 « (q )p n +qpn.q> |
(29.13) |
|
* |
В случае |
необходимости можно ввести |
поправки, |
учитывающие ма |
лость |
Vo- |
|
|
|
297
где оператор q-й компоненты флуктуации плотности в n-й ячей ке имеет вид
~~+
Pn,q — ^n.k—q.a^n.k.a; k,o
T и U — операторы кинетической энергии электрона и периоди ческого внешнего поля соответственно. При этом в операторы рождения а+ и поглощения а электрона в n-й ячейке явно вве дены спиновые индексы. Последний член в выражении (29.13) учитывает двухчастичное взаимодействие и представляется сум мой по q операторов (29.7), поскольку оператор взаимодейст вия в импульсном представлении d(q)—4ne2/q2 положительно определен. Выбирая далее набор эрмитовых проекционных опе
кая |
|
|
проекцию |
опера |
|
раторов 0 q |
, легко построить внутреннюю |
||||
тора Нп |
|
|
|
|
|
% = 2 |
<ki I ^ |
I кзЛ к.ап .к, |
(NJ2V0) ^ u ( q ) |
+ |
|
к ,к . |
|
|
q |
|
|
|
+ |
(1/2) ^ и (q) b .q O f рал. |
|
(29.14) |
|
|
|
q |
|
|
|
Перейдем теперь к верхним и нижним |
оценкам |
энергии |
|||
основного |
состояния |
электронного газа на |
компенсирующем |
фоне положительного заряда в двух предельных случаях. Что бы не усложнять выкладок, исключим из рассмотрения перио дическое внешнее поле, т. е. проведем вычисления в упрощен ной постановке задачи. Отметим, что предлагаемый здесь фор мализм диктуется соображениями удобства вычислений и не является единственно возможным. Тем не менее такое кон кретное вычисление является наглядной иллюстрацией эффек
тивное ги |
рассмотренного метода и указывает на малую роль |
||||
пеадднтивпон |
части |
энергии W даже в случае сильной связи |
|||
(rs^> 1). когда число частиц в ячейке мало. |
|
||||
Оценка |
сверху для энергии |
основного состояния! плотного |
|||
однородного |
электронного газа |
на положительном |
компенси |
||
рующем фоне (rs<cl, |
аддитивное приближение). В представле |
||||
нии ячеистых плоских волн гамильтониан n-й ячейки |
имеет вид |
||||
|
|
|
Яп — р.<* 6 (р) N п,р,0 "Ь |
|
|
+ (1/2У0) |
^ |
;+ |
|
|
|
и (q) fln/p.+q.CT, <2п/ра—ч ,о аЯп>р2,а аЯп,рг(29,сг,.15) |
d .P i .<*1 .P t.
где fln.p.a — оператор уничтожения частицы в состоянии (29.12);
р = Пк- е (р) = ра/2т, и (q) = 4ne2/q\
Яп,р,о #П+, р ,0 ^П,р,СГ.
298
Будем строить внешнюю проекцию оператора #„ на подпро странство ЭЛ, определяемое системой базисных векторов*:
( |
|0п> |
л |
л л |
л |
' |
' (29.16) |
( |
| Фп.Ч (Pi* Р 2) |
= ^ . P t + q O n . P j O n , —Pj—q a n lt —P, I ^ n ) > |
|
|||
где |
|
|
|
|
|
|
|
|0п> = |
П в ( 8р - (г)а+р I |
. . . > |
‘ |
(29.17) |
|
|
|
p |
|
|
|
|
— волновая функция n-й ячейки в приближении |
Хартри — Фо |
ка; р — химический потенциал n-й ячейки, в равновесии одина
ковый для всех ячеек; |
| ... > |
— волновая функция вакуума |
|||||
|
|
|
|
J |
1 при * > |
0; |
|
|
|
|
|
( 0 при х < |
0. |
||
Подпространство |
состояний |
ЯЛ 2 с |
базисом (29.16) соответст |
||||
вует, |
таким образом, |
конфигурациям, |
которые отличаются от |
||||
| Оп > |
не более чем возбуждением |
двух пар с суммарным им |
|||||
пульсом, равным нулю. |
|
|
|
|
|||
Оператор О®2 |
запишем в виде |
|
|
||||
|
0 f 2= |
22 |
I Фп.ч (Pi> Рг)> <Фп.ч (Pi, Р2) I + | On) <0„ |. (29.18) |
||||
|
|
Pi.P j.q |
|
|
|
|
|
Этот |
оператор |
эрмитов |
и проектирует |
пространство состояний |
n-й ячейки на подпространство ЯЛП. С помощью этого операто ра строим внешнюю проекцию
= О®11н З Т г. |
(29.19) |
Собственныевекторы оператора Н„ принадлежат |
подпрост |
ранству ЯЛ 2. Поэтому вектор основного состояния гамильтониа на (29.19) можно представить в виде
| ^ п> = |
|Оп> + |А¥„>, |
(29.20) |
где |
|
|
№ > - ( l / K 0) V |
Ф,(Р1, p2)|cpn.q(Pi,p2)>. |
(29.21) |
q.Pi.p2 |
|
|
На функцию Фч(рь рг) наложены условия: |
|
|
Фа(Pi. Р2) = Фч (Рг. Pi); |
Ф-q (рх, Рг) = Фч (— Рх, — Рг)! |
|
_________ Ф-q (— Рх, — Рг) = |
Фа(Рх, Рг), |
(29.22) |
* Здесь и далее под р понимаем (р, о ).
299