Микроскопическая теория дисперсионных сил есть теория систем многих взаимодействующих молекул. Описание таких систем дает статистическая механика.
Для вывода необходимых соотношений статистической механики удобно придерживаться функционального подхода и использовать производящие свойства большой статистической суммы
ОО
S = X "ATT f • • • [ Dl . . . (1) . . . £ (ЛО ехр [~Р£/ (1, • • Ю] 0)
N—0 ‘ Рассматриваемой как функционал от функции
£ (0 = г%. ехр [—р« (01 (2)
Играющей роль обобщенной активности системы.
Здесь аргумент i обозначает совокупность переменных, относящихся к степеням свободы i—Й молекулы, а именно: ее сорт st — и пространственные координаты х;, включающие радиус-век-
Тор положения rt и, вообще говоря, углы ориентации 0г молекулы. Соответственно символ J di обозначает суммирование по сортам и интегрирование по пространственным координатам: j di = ^ J dxi— Далее, и (I — внешнее поле; zs. — актив-
Юность, связанная с химическим потенциалом ps соотношением
2S =exp({5jus.)/A3 (7) ‘(3)
; it
Где Л5. (Т) — известная функция температуры Т, возникающая в результате интегрирования по импульсам, канонически сопряженным пространственным координатам. Наконец, § = = /kBT — обратная температура (kB — постоянная Больц — мана); U (1, …, N) —энергия взаимодействия системы из N Молекул.
Функции распределения
00
Р<1……. п) = а-1 2 (Л/1П)! D («+»)… DNX
X G (1) — . — s ехр L—ре/ (1…………. Л/)] (4)
Описывающие вероятности распределения молекул, в методе производящего функционала могут быть компактно представлены в виде:
Р(1, .. ,/i) = C(l)…Ј(/i)B-,6n3/6Ј(l) … Б £(Л) (5)
Удобно ввести понятие о комплексе п, понимая под ним группу п молекул с заданными сортами и пространственными координатами. При этом р (1, …, п) = р {п.
Формулы (1)—(5) будут служить основой для дальнейшего.
Из (2) и (3) следует операторное соотношение
= J dxtt, (i) б/бс (0 (6)
Где производная по £ (i) берется при постоянстве температуры, а производная по — при постоянстве температуры и внешнего поля.
Применение этого соотношения к функциональному представлению (5) дает
J dx, Ip {п + 1} — р {/»> р (L)] = Р-1 — nSl {п} р {п} (7)
Где nSl{«} — число молекул сорта Sj в комплексе {п. Соотношение (7) выражает в наиболее общем виде теорему о сжимаемости.
Реалистическое описание межмолекулярных сил дает парный эффективный (учитывающий усредненно неаддитивные вклады) потенциал взаимодействия Ф (I, J). В терминах такого потенциала:
1 N
1/(1, …. A0 = — f £ Ф(/, /) (8)
^ »>/=1
Подстановка (8) в (1) позволяет представить экспоненту в (1) как сумму всевозможных произведений функций Майера F (I, /) = exp [—рФ (I, /)] — 1. Это приводит к известному для S выражению в виде суммы всех диаграмм, с вершинами которых сопоставляются обобщенные активности С. а с соединяющими вершины линиями — функции Майера /. Такие линии называются майеровскими связями.
Каждая из вершин в диаграммах для S является полевой. Под этим понимается, что в ней ведется суммирование по сортам и интегрирование по пространственным — координатам молекул.
Очевидно, что применение оператора £ (0 (0 к полевой вершине делает ее фиксированной. В такой вершине сорт и координаты молекулы, от которых зависит сопоставляемая с вершиной обобщенная активность принимают фиксированные значения s( и Xt.
Поскольку в функциональном методе функции распределения находят по формуле (5) путем последовательного применения указанного оператора к Е (предполагается, что все аргу
менты 1, 2, п различны), это дает простой способ вывода их диаграммных разложений. Таким путем для р (1, …, п) Получается выражение в виде суммы всех диаграмм, которые состоят из фиксированных вершин £ (1), …. £ (я), из полевых вершин из майеровских связей /ив которых каждая полевая вершина соединена (непосредственно или через другие полевые вершины) хотя бы с одной фиксированной вершиной. Указанное ограничение об отсутствии полевых вершин, не соединенных с фиксированными, очевидно, обязано множителю Е-1 в (5).
Построенные оперативно методом производящего функционала диаграммные разложения для функций распределения будут использованы ниже в качестве математического аппарата вычисления. Важно, что конечные результаты вычисления, в которых будут просуммированы существенные для рассматриваемых задач диаграммы во всех без исключения порядках, удастся представить в терминах физически наблюдаемых величин. Существование таких величин не связано со сходимостью разложений. Поэтому последующее применение диаграммных разложений как рабочего аппарата вычисления не будет ограничено областью их фактической сходимости. Оно справедливо и для конденсированных систем, в которых активности, играющие роль параметров диаграммных разложений, не являются малыми.
Описывая корреляционные и структурные свойства, функции распределения позволяют найти любые локальные характеристики системы. Соответствующие выражения следуют непосредственно из определения (4) для функций распределения.
Важнейшей локальной характеристикой капиллярных систем является тензор давлений. В общем случае нецентральных взаимодействий и молекул, имеющих ориентационные степени свободы, выражение для тензора давлений в некоторой про-
… дФ(1, 2) ^ / DD2———- ® <л2 Дг о |
Странственной точке г имеет вид:
PCr) = KBT^ РР Й
DrБ [г- Л (г2 — %) — rt] Р (1, 2)
Здесь 1 — единичный тензор;
Рр (г) — усредненная (проинтегрированная) по ориен-
Тациям одночастичная функция распределения
-+ -»
Компонента р в точке г производная по г2
Означает градиент по координате г2; (g) — символ, означающий прямое умножение векторов (в результате которого образуется тензор); б обозначает обычную б-функцию.
Тензор давлений удовлетворяет условию механического равновесия
Div Р = 0 (10)
Которое явствует из выражения (9) и определения (4). Физически условие равновесия выражает фундаментальный закон сохранения импульса в равновесных системах.
Условие равновесия в форме (10), в которой отсутствует член с градиентом внешнего поля, предполагает, что в выражении (9) для тензора давлений учитываются вклады от всех участков пространства, в которых присутствует вещество, способное оказывать силовое воздействие.
В случае центральных взаимодействий, когда Ф (I, J) зависит от пространственных координат только посредством расстояния можду молекулами, усреднению по угловым координатам в (9) подлежит лишь двухчастичная функция распределения. И тогда выражение (9) переходит (после замены инте-
Грирования по г2 интегрированием по R = г2 — гх) в обычное выражение Ирвинга—Кирквуда
Р Р, я
1
Х-^Э" Jdnpw(R—I|J?, 1-NR + R) (11)
О
Где индексы р и Q характеризуют компоненты, функции распределения усреднены по ориентациям молекул. Определяемый формулой Ирвинга—Кирквуда тензор давлений является симметричным.