Действительно, согласно вариационному принципу
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, |
Лц = |
/и, |
|
|
что и доказывает сделанное выше утверждение.
Определим в н у т р е н и юю |
п р о е к ц и ю |
эрмитова положи- |
тельпо определенного оператора |
А(А>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. Соответственно при построении внешних про
екций для гамильтониана Н получаем верхние оценки для энер
гии 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 <">, которая стремится к единичному оператору, когда система приближается к полной.
Удобно говорить о полноте системы {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) отражает виртуаль |
ные возбужденные состояния над системой одночастичных со стояний, описываемых волновой функцией |Оп> .
Чрезвычайно важно, что скалярное произведение < 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 и в целом электронейтральиа. |
Поскольку представляют интерес объемные свойства такой |
системы, будет оценена, |
в частности, |
энергия, |
приходящаяся |