CC BY
43
ФИЗИКА
J
УДК 539.19+539.2
О. С. Е р к о в и ч, А. М. Руцкая
ТЕОРЕМА ВИРИАЛА И МАСШТАБНЫЕ СООТНОШЕНИЯ В МНОГОКОМПОНЕНТНЫХ СИСТЕМАХ ЗАРЯЖЕННЫХ ЧАСТИЦ
Настоящая работа посвящена исследованию единственности решения вариационной задачи в методе функционалов плотности. Исследованы масштабные соотношения, возникающие в системах, включающих частицы двух видов, а также в системах заряженных частиц, находящихся в электростатическом поле, созданном распределенным внешним зарядом. Рассмотрена возможность применения теоремы вириала. Показано, что для квантовых систем заряженных частиц, находящихся в поле распределенного заряда, не выполняется известное соотношение 2T = —V — W, связывающее между собой средние значения кинетической энергии системы T, потенциальной энергии взаимодействия частиц с внешним полем (V) и между собой (W).
E-mail: press@bmstu.ru
Ключевые слова: метод функционалов плотности, теорема вириала, масштабные соотношения.
При исследовании энергетических характеристик локализированных состояний систем взаимодействующих фермионов полезными оказываются соотношения, устанавливающие связь между различными вкладами в полную энергию системы, наиболее известным из которых является теорема вириала. Подобные соотношения, с одной стороны, облегчают выполнение численных расчетов, с другой — они могут выполнять функции дополнительного надежного критерия при определении границ применимости тех или иных приближенных теоретических методов. Наибольшую ценность эти критерии представляют при исследовании систем, для которых не имеется достаточно полных экспериментальных данных, к числу таких систем можно отнести, в первую очередь, нанопорошки и нанотрубки.
Настоящая работа посвящена исследованию масштабных соотношений для систем, включающих частицы двух видов, а именно систем заряженных частиц, находящихся в электростатическом поле, созданном распределенным внешним зарядом.
Решение поставленной задачи проводилось в рамках нерелятивистской квантовой механики с привлечением теории функционалов плотности.
Рассмотрим квантово-механическую систему, состоящую из двух подсистем, представляющих собой совокупности частиц двух видов
(например, электроны и ядра атомов), взаимодействующих между собой посредством кулоновских сил. Гамильтониан такой системы в отсутствие магнитных сил имеет вид
Н = Т + Т2 + VI + V + W. (1)
N-1 , 1 . N2 , 1 .
Операторы Т = ^ ( - ^АЛ и Т2 = ^ ( - — АЧ описывают
7=1 ^ ' 7=1 ^ '
N1 N2-1 1
кинетическую энергию частиц подсистем 1 и 2; У1 = ^^ ^^ —-—
г=2 ]<А 'Гг ^|
N2 N2-1 2
и У2 = У У -=тт — операторы потенциальной энергии взаимо-
7=2 и - (|
N1 N2
действия частиц подсистем 1 и 2; Ш = ^^ ^^ —-— — оператор
• 1 • 1 |Г * ¡71
7=1 7=1 -1 '
потенциальной энергии взаимодействия между подсистемами.
Основное состояние пространственно локализованной системы с гамильтонианом (1) описывается нормированной на единицу волновой функцией фдв.
Для исследования свойств этой системы проведем масштабное преобразование волновой функции основного состояния, введя в рассмотрение семейство нормированных на единицу функций
ф7 (гь ...,^N1; ¡ь..., ¡N2; а) = 1т/2Фд8 (тгъ..., т^м-; т?1,..., 7(к; а)
где а — совокупность спиновых координат частиц; 7 — вещественный параметр. Очевидно, что
ф7 (гь..., ; яи..., (¡N2; а)|7=1 = Фдз (ги ...^N1; ql,...,qN2;а).
Используя явный вид операторов Т, V и Ш, можно получить соотношения
<ф7 |Т1| ф7> = 72 <Ф iТ11 ф) ; (ф7 |Т21 Ф7> = 72 <Ф iТ21 ф> ; <ф7 | V! | ф7> = 7 <Ф | ф> ; <ф7 |V2| Ф7> = 7 <Ф ф> ; <ф7 |Ш| ф7> = 7 <Ф |Ш| ф>.
Воспользовавшись вариационным принципом Рэлея-Ритца
в
— <ф7 |Н |ф7 > =0,
"7 7=1
получим известную теорему вириала в виде
2 (7 + 72) = -V! - V - ".
Функции основного состояния фд8(г1,..., 1; <1,..., 2; системы с гамильтонианом Н соответствует т-частичная функция плотности пт(г1,..., <т), связанная с волновой функцией соотношением
(r!,...,rm) = CN ^ / ...,rjvi;gi,...,gw2;^)d3rm+i
... 1, ...
где суммирование проводится по всем спиновым переменным. Она удовлетворяет условию нормировки
nm (гЪ . . . , rm) d ri . . . d rm — Cm —
N!
т! (^ — т)!'
поэтому подвергнутой масштабному преобразованию функции ф^(г1,..., 1; ..., 2; будет соответствовать т-частичная функция плотности
пт7 (г1, . . . , гт) 7 пт (ТГ1, . . . , Тгт).
В соответствии с обобщенной теоремой Хоэнберга-Кона [1-7] волновая функция невырожденного основного состояния и, следовательно, кинетическая энергия, а также потенциальная энергия взаимодействия частиц системы с внешним полем и между собой являются однозначными универсальными функционалами т-частичной функции плотности основного состояния:
Ti[nm] — T2[nm] — Vi[nm] — V>[nm] — W [nm] —
Для функционалов Ti [nm], T2[«m], Vi [nm], V2 [nm] и W [nm] существуют масштабные соотношения
Ti[nm7 ] — Y 2Ti[nm];
T2[nm7 ] — Y 2T2[nm];
] |Ti| ^M); ] |T2| ^W); ] |Vi| ^H);
nm
nm ] |W| ^[«m]) .
Vi[nm7 ] — YVi[«m];
V [Пщ7 ] = 7 ^ [ищ];
Ш [ит7 ] = 7Ш [пт].
Тогда, используя принцип Рэлея-Ритца, получаем
2 (Т1 [Ищ] + Т2 [Пщ]) = - (V [Пщ] + V2 [Ищ] + Ш [Пщ]) .
Из этого следует, что соотношение пропорциональности связывает суммарные кинетические и потенциальные энергии составной системы, при этом условие
2Т1 [Пщ] = - (V [Пщ] + Ш [Пщ]), (2)
считающееся справедливым для любых систем заряженных частиц, в том числе и электронного газа, находящегося в поле кристаллической решетки, автоматически не выполняется [8-10].
В то же самое время неприменимость соотношения (2) для описания многокомпонентных систем приводит к необходимости получения аналогичных оценок для различных вкладов в полную энергию системы. При этом особый интерес представляют консервативные системы, состоящие из частиц одного типа, находящиеся в поле, созданном некоторым распределенным зарядом. Важнейшим частным случаем такой системы является электронный газ, находящийся в поле неподвижных атомных ядер. Отметим, что общепринято мнение о применимости к системе заряженных фермионов, движущихся в поле, созданном точечными неподвижными зарядами, теоремы вириала в виде (2) [4-10].
Возвратимся к анализу квантово-механической системы, состоящей из двух подсистем, представляющих собой совокупности частиц двух видов, взаимодействующих между собой посредством кулонов-ских сил. Гамильтониан такой системы имеет вид (1).
Учитывая, что частицы подсистем 1 и 2 различны, можно записать
фд5 (гь ..., ¡N1; дь. ..,(^2; а) = Фд5 (д; д; а) = ф (г ) р (д) х (а),
где г = г1,..., и д = ¡1,..., — совокупность координат подсистем 1 и 2. Условия нормировки волновой функции удобно задать в виде
<ф (г )| ф (г )> = 1; <р (д)| р (д)> = 1; <х (а)| х (а)> = 1,
откуда следует ||Ф (г, д )|| = 1.
Рассмотрим матричный элемент
/7 = <ф7 |Н |ф7> = 72Т + Т2 + 7^1 + V2 + Ш7.
Здесь
W7 = ( ф7 ) p (q ) X И
E
|g - g1
ф7 (r) p (g) x И ) =
= ( ф7 (r) p (g)
z
|r - gj1
i>j
ф7 (Г) P (g ) ) =
= 73N4 ф/R )p (g)
E
«>j
|g - g1
ф R p (g И =
= y 3Nl yu(r)p (g)
E
«>j
|g- g1
ф R p (g ) ) =
= y3n1 y i „/,*
[ ф* (R )p* (g) E
y |g - g1
ф R p (g) dr3 dg3 =
= Y
[ ф* (R)p* (g) E
R i - Ygj
ф R p (g) dr3dg3.
В соответствии с вариационным принципом Рэлея-Ритца
d
dY
— (ф7 |H |ф7 >
= 0.
7=1
Очевидно, что
d , , , , , dI7 ^ - -W7
— (ф7 |H |ф7 > = = 2ytTi + V + —7, dY -y -y
где
-w
-y
i ф* (R)p* (g) E
J i,j
Ri - Ygj
ф R p (g) dr 3dg 3+
[ ф* (R)p* (g) E
J i,j
-Ri - Ygj
3j
Yg2 - -Rigj) -^(-R)p (g) -r 3dg3-
-W7
-Y
=
7=1
+ / ф* R p* (g)
z
i,j
|gi-gj|
3j
g2-Rigj •^(R)p (g) dr 3dg 3=W+AW
z
z
z
z
z
z
z
откуда
2Ti + Vi + W + AWT = 0. (3)
Полученный результат позволяет получить оценочные соотношения, связывающие между собой различные вклады в полную энергию такой системы.
Отметим, что слагаемое
AW = i ф* (R) V* (д) £ 3 (j - %) Ф (R) V (?) dr 3dg3
i j 1 i *j 1
(4)
оказывается знакоопределенным: AWT ^ 0.
Таким образом, широко используемые результаты, представленные в классических монографиях [8-10], распространяющие теорему ви-риала в виде (2)
2Ti [nm] = - (V1 [nm ] + W [nm])
(доказанную для многоэлектронных атомов) на произвольные системы заряженных частиц (атомы, молекулы, твердые тела), требуют уточнения. Это обобщение оказывается справедливым только при выполнении условия
AW= / ф* (R) V* (д) ^ |Г -^ |3 (q? - %) Ф (R) V (?) dr 3dg3=0,
что возможно только для некоторых частных случаев.
Таким образом, можно сделать вывод о том, что в консервативных системах заряженных частиц, находящихся в поле распределенного внешнего заряда, средние значения кинетической энергии частиц и потенциальной энергии их взаимодействия между собой и с внешним полем не удовлетворяют известному соотношению 2T = — (V + W), т.е. теореме вириала.
В системах заряженных квантовых частиц, состоящих из двух подсистем, теорема вириала выполняется только для системы как целого:
2 (Ti [nm] + T? [nm]) = - (V [nm] + V? [nm] + W [nm]) ,
но не выполняется для отдельных подсистем.
Для подсистемы 1, находящейся в поле подсистемы 2, вклады в полную энергию подсистемы 1 определены соотношениями (3) и (4):
2T + Vi + W + AW = 0,
где слагаемое AW является знакоопределенным.
Полученные результаты могут быть применены в атомной и молекулярной физике, физике твердого тела и физике полимеров, а также при изучении нанообъектов.
СПИСОК ЛИТЕРАТУРЫ
1. Исследование корреляционных свойств многочастичных систем в рамках метода двухчастичных функционалов плотности / О.С. Еркович,
B.В.Комаров, А.М.Попова и др. // Вестник МГУ. Сер. 3: Физика, астрономия.
- 1991. - Т. 32., № 4. - С. 42-49.
2. Еркович О. С., К о м а р о в В. В., П о п о в а А. М., Борзилов В. А. Обобщенная теорема Хоэнберга-Кона в методе многочастичных функционалов плотности // Вестник МГУ. Сер. 3: Физика, астрономия. - 1994. - Т. 35, № 2. -
C. 33-39.
3. Е р к о в и ч О. С., Комаров В. В., Попова А. М. Анализ пространственной структуры электронного газа вблизи поверхности металла в методе двухчастичных функционалов плотности // Поверхность. - 1994. - № 10. -С. 23-28.
4. Е р к о в и ч О. С. Теорема вириала и масштабные соотношения в методе многочастичных функционалов плотности // Вестник МГУ. Сер. 3: Физика, астрономия. - 1996. - Т.37, № 5. - С. 70-72.
5. Е р к о в и ч О. С. Формулировка вариационного принципа в методе функционалов плотности // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. -2000. -№ 1(4). - С. 84-96.
6. Е р к о в и ч О. С. Метод многочастичных функционалов плотности: вид функционала кинетической энергии // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. - 2000. - № 2 (5). - С. 73-79.
7. Балашов В. В., Д о л и н о в В. К. Курс квантовой механики. - М.: Изд-во Моск. ун-та, 1982.-280 с.
8. Марч Н., Кон, Вашишта П., Лундквист С., Уильямс А., Барт У., Л э н г Н. Теория неоднородного электронного газа: Пер. с англ. / Под ред. С. Лундквиста и Н. Марча. - М.: Мир, 1987. - 400 с.
9. ParrR. G.,WeitaoYang. Density-functional theory of atoms and molecules.
- Oxford University Press, 1989.
10. D r e i z e r R. M., G r o s s E. K. U. Density functional theory. - Berlin: SpringerVerlag, 1990. - 303 p.
Статья поступила в редакцию 12.01.2010
Ольга Станиславовна Еркович родилась в 1962 г., окончила МГУ им. М.В. Ломоносова в 1984 г. Канд. физ.-мат. наук, доцент кафедры "Физика" МГТУ им. Н.Э. Баумана. Автор более 50 научных работ в области нерелятивистской квантовой механики.
O.S. Erkovich (b. 1962) graduated from the Lomonosov Moscow State University in 1984. Ph.D. (Phys.- Math), assoc. professor of "Physics" department of the Bauman State Technical University. Author of more than 50 publications in the field of nonrelativistic quantum mechanics.
Анна Михайловна Руцкая родилась в 1985 г., окончила МГТУ им. Н.Э. Баумана в 2008 г. Ассистент кафедры "Физика" МГТУ им. Н.Э. Баумана, магистр техники и технологий. Автор двух научных работ в области нерелятивистской квантовой механики.
A.M. Rutskaya (b. 1985) graduated from the Bauman Moscow State Technical University in 2008, holder of Master's degree in engineering and technology. Аssistant lecturer of "Physics" department of the Bauman Moscow State Technical University. Author of 2 publications in the field of nonrelativistic quantum mechanics.
