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

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

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

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

Добавлен: 15.10.2024

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

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

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

Действительно, согласно вариационному принципу

fli — < ui I А | U l> ^ < u, | A | « ! > —

= < « i 1 О х Л О | u x > = < « ! | Л - 1 ых > = аг.

При этом использованы соотношения (27.6). Рассмотрим вспо­ могательную функцию ф= гг 1р1+ Й2Р2, удовлетворяющую усло­ виям

*C.u i I Ф^> “

I

?i _h <^1 I

P2 = 0;

< ф I ф > = I Pi |2 -I-1Р2 Г = 1;

Оф = ф .

Тогда, согласно вариационному принципу

 

я2 < < Ф М I Ф > = < О ф I А | О ф > = < Ф | Ъ+АО | ф > =

= <]ф I Л | Ф^> = •\WjPi

-f- W2P2 I А |

HlPl

Г u2p2^> =

=

«I | pi I +

Й2 I р.2 I2 <

а 2.

 

Если рассмотреть вспомогательную функцию, ортогональную

«1 и ы2. то аналогично получим

 

 

 

 

 

^ ^з-

 

 

Это рассуждение справедливо для любого k, что и доказывает

неравенство (27.7).

 

 

 

 

 

-S

 

 

 

 

 

 

 

 

Вариационный метод Ритца для оператора А эквивалентен

диагонализации оператора

А с

соответствующим

образом вы-

 

 

 

 

рассмотрим

экстре­

бранным оператором О. Действительно,

мальное значение интеграла

 

 

 

 

 

 

 

___

А

 

 

 

 

 

 

I = <м | А | м>/<ы | м>,

 

 

где и принадлежит

подпространству

ЭЛ,

выделенному

проек-

 

/ Ч

А .

 

 

 

 

 

циопным оператором О, т.А е.

ОЛ= й.

Введем

произвольную

функцию и и положим й = Ои. Тогда

 

 

 

 

 

, _

< и | 0+ ДО | и> _ < и | А | «>

 

 

 

< и | О | и >

 

< к | О | « >

 

 

Согласно вариационному

принципу,

6/=0

получаем

 

^

^

 

или

d

 

_

 

 

(Л —/О) м — 0,

Лц =

/и,

 

 

что и доказывает сделанное выше утверждение.

282


Определим в н у т р е н и юю

п р о е к ц и ю

эрмитова положи-

тельпо определенного оператора

А(А>0).

Внутренней

проек­

цией называют оператор вида

 

 

 

 

 

 

 

А =

А /гОА /г ,

 

 

 

который обладает

следующим

свойством.

Если имеет

место

^

 

собственных значений указанных

неравенство О< A < A , то для

операторов выполняется аналогичное соотношение

 

 

 

 

0 < а й< а й.

 

 

(27.8)

Операторное

неравенство

/ Ч

/ Ч

означает

неравенство

для

А<В

 

 

 

 

 

,--

ХЧ

 

для

диагональных матричных элементов <ср|ср> < <ср | В j ф>

всех возможных ср. Тогда верно и неравенство для собствен­

ных значений ак<Ьк. Очевидно,

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

справедливо. Оператор А 112

(«корень

из оператора») имеет

следующий смысл. Если

 

 

 

ч

|

ик~^>

|,

А =

k

 

 

 

то

А ',г = 2k а,‘к,г | ик> < и к |.

Неравенства (27.8) вытекают из положительной определенно-

-— ч

сти проекционного оператора О. Действительно,

< Ф | О | ф > = < 0 | О2 | ф > = < ф | 0+0 | ф > = < 0 ф | 0 ф > > 0.

Аналогичное рассмотрение матричных элементов от опера-

/ Ч -''Ч . / Ч

к

 

условию

л -~ч.

тора Р= 1—0 ^ 0 приводит

 

O ^ O ^ i, т. е.

 

 

получаем

0 < < ф |0 |ф Х < ф |ф > . Заменяя

ф на Л‘/2ф,

О< < Ф | А 4*О А''1

| ф >

< < ф | А | ф >,

что и доказывает утверждение (27.8).

Это свойство внутренней проекции оператора находит при­ менение в методе промежуточных гамильтонианов. Если

имеется система эрмитовых операторов

 

Я (1) < Я (2)< . . . < Я ,

(27.9)

(27.8)

сходящихся к данному гамильтониану Я, то неравенство

приводит к соответствующему расположению собственных зна­ чений:

4 ° < ^ 2)< . . . -< Ек,

(27.10)

что обеспечивает систему нижних оценок, сходящихся к точно­ му значению Eh. Соответственно при построении внешних про­

283


екций для гамильтониана Н получаем верхние оценки для энер­

гии Eh.

Если представить гамильтониан системы в виде

н = н 0 + и,

где U> 0 (положительно определенный оператор взаимодейст­ вия), то в качестве промежуточных гамильтонианов можно вы­ брать операторы вида

НМ

н 0ж и 4, Oin) U‘/а ,

(27.11)

где 0 <п>— последовательность проекционных

операторов, стре­

мящихся к единичному.

использовал Лёвдин

для вычисления

Подобный подход

энергетических уровней гелиеподобных ионов. При этом вычис­

лялись как

основные, так и

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

состояния

систем,

для которых гамильтониан

записывался

в виде

/X /Ч

H = H0+U,

 

£/ = — , Z = 1, 2, 3, . . .

 

 

 

Tia

 

 

 

 

Интересно

отметить, что уже

третье приближение

п —3

приво­

дит к значению энергии, отличающемуся от точного менее чем

на 1%.

Десятое же приближение

в случае

иона

лития (Li+,

Z = 3)

дает значения —7,28444

ат. ед.

для

верхней в

—7,279910 ат. ед. для нижней оценки полной энергии. Можно, по-видимому, утверждать, что в случае систем с конечным (не­ большим) числом частиц метод построения внешних и внут­

ренних проекций для оператора энергии

очень продуктивен.

Проекционный оператор можно построить с помощью орто-

нормированного набора

(конечного

или

бесконечного) векто­

ров состояний

 

 

 

 

(I <Pi>,

I Ф2> .

• •

.. I Ф«>1-

Тогда оператор О имеет вид

 

 

 

%

= 2i

I Ф / > < Ф ; I-

Этот оператор проектирует гильбертово пространство векторов состояний системы на подпространство Ш1, натянутое на базис­ ные векторы {|фг>}) причем <фг |фк> = 6г,„. Система функций {|ф,->}, вообще говоря, неполна. Добавляя к ней новые конфи­ гурации |фк> , будем получать все более полные системы. Та­ ким путем и получается последовательность проекционных

операторов 0 <">, которая стремится к единичному оператору, когда система приближается к полной.

284


Удобно говорить о полноте системы {I Фг>} по отношению к некоторому вектору состояний |Ч/ > (это может быть, напри-

мер, основное состояние, описываемое гамильтонианом Н). Тогда степень полноты системы по отношению к |ЧТ> можно характеризовать параметром

VsK= < ¥ |

| Ода^ > = 1 - < У 1Ош | ¥ > .

(27.12)

Параметр удд уменьшается при добавлении к

{|фг>}

новых

конфигураций и равен нулю для полной системы

{| фг>}■ Оче­

видно, что собственное значение промежуточного гамильтониа­

на

достаточно близко к точному

значению Eh лишь в том

случае, если

 

 

 

Y $ = 1 ~ < ^ I %

I ^ * > « 1 -

(27-13)

Мы не умеем точно вычислять параметр у эд.

Простая порядковая оценка показывает трудности, возни­ кающие при попытке использовать метод внешней и внутренней

проекций (Я и Н) в случае систем с большим числом частиц Nv

N -*■ оо, У^ - оо , п = N/V = const.

Чтобы уравнение Шредингера с промежуточным гамильтониа­

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

¥ > = |Оп> 4 -

(Pi, Pa) 'dt (Pi) d—q (Рг) Оп>. (27.14)

q , Р . р2

Волновая функция |ЧГ> записана здесь в импульсном пред­ ставлении; |Оп> — волновая функция системы в приближении

Хартри — Фока; V — объем системы; dp

и d_q— соответственно

операторы рождения

и уничтожения

пары.

Функции

вида

(27.14) принадлежат подпространству с базисом

 

( | 0 п >

~+

-

 

 

(27-15)

II Фч (Pi. Ра) > = d+ (Pi) d_q (р2) | On> .

 

Пара частиц рождается с суммарным

импульсом q и гибнет

с импульсом —q, так что закон сохранения

импульса

соблю­

дается. Второй член в выражении

(27.14) отражает виртуаль­

ные возбужденные состояния над системой одночастичных со­ стояний, описываемых волновой функцией |Оп> .

285


Чрезвычайно важно, что скалярное произведение < 0 11|Чго> в случае макроскопически, больших систем оказывается экспо­

ненциально малой величиной:

 

 

 

< 0 П| Ч'0> ~ е х р (-1 /),

(27.16)

где

|ЧЛ>>— волновая

функция,

описывающая

конечное число

парных возбуждений.

Легко видеть, что

 

q

2 I <Фч(р1, Р-з) I ^о> I2=

< lFo I О | VF0> ~ехр(—V).

, р. ■р2

 

 

 

Поэтому введенный выше параметр у®г близок к единице. Ины­

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

По этой причине среднее значение гамильтониана системы по состоянию (27.14) в пределе N—^oo равно энергии Хартри — Фока и не зависит от конкретного вида cDq (рь Рг). Эта функция

•определяет лишь вклад в < lF |# |xF > , не зависящий от числа

•частиц N. Нас же интересуют аддитивные свойства системы.

§ 28. ЭЛЕКТРОННЫЙ ГАЗ В ПЕРИОДИЧЕСКОМ ВНЕШНЕМ ПОЛЕ. ОПРЕДЕЛЕНИЕ И СВОЙСТВА ЯЧЕИСТОГО БАЗИСА

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

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

Н =

Г и (г,.)

(28.1)

 

 

i+ i

тде Mi — масса; pi — импульс i-й частицы; U(г,)— потенциаль­

ная энергия г-го электрона во внешнем поле;

N — число элек­

тронов в системе;

 

 

 

«(I г, — 17.|) = е2/\ г,-

гj |

 

—'Энергия взаимодействия электронов,

находящихся в точках

г,- и rj. Предполагается,

что система помещена в большой ку­

бический ящик объемом

V и в целом электронейтральиа.

Поскольку представляют интерес объемные свойства такой

системы, будет оценена,

в частности,

энергия,

приходящаяся

2 8 6