Цепочка Боголюбова — Борна — Грина — Кирквуда — Ивона

Цепочка Боголюбова — Борна — Грина — Кирквуда — Ивона (сокращённо ББГКИ, также иерархия ББГКИ) — это система зацепляющихся интегро-дифференциальных уравнений, описывающая эволюцию неравновесных статистических ансамблей в классической и квантовой статистической механике. Она связывает одночастичные, двухчастичные и более высокие функции распределения, позволяя последовательно учитывать корреляции между частицами. Цепочка названа по именам физиков, внёсших вклад в её разработку: Николая Николаевича Боголюбова, Макса Борна, Герберта Сиднея Грина, Джона Гэмблина Кирквуда и Жака Ивона (иногда также упоминается Йозеф Майер). Она является фундаментальным инструментом для изучения кинетических процессов, гидродинамики, теории плазмы и фазовых переходов.

История

Истоки цепочки ББГКИ восходят к работам конца XIX — начала XX века, когда Джозайя Уиллард Гиббс заложил основы статистической механики. Однако систематическое описание неравновесных процессов потребовало учёта многочастичных взаимодействий. В 1935 году Жак Ивон, французский физик, вывел уравнения для функций распределения в классической жидкости. В 1946 году американский физик Джон Гэмблин Кирквуд независимо получил аналогичные уравнения, используя метод проекционных операторов. В том же году Макс Борн и Герберт Сидней Грин опубликовали работу, в которой представили иерархию уравнений для газов и жидкостей. Решающий вклад внёс советский математик и физик Николай Николаевич Боголюбов, который в 1946 году в монографии «Проблемы динамической теории в статистической физике» предложил строгий вывод цепочки на основе функциональных производных и метода адиабатического включения взаимодействия. Боголюбов также разработал метод расцепления цепочки, позволяющий получать замкнутые кинетические уравнения, такие как уравнение Больцмана.

Формулировка

Классический случай

Для системы из \(N\) одинаковых частиц в фазовом пространстве (координаты \(\mathbf{q}_i\) и импульсы \(\mathbf{p}_i\)) вводится функция распределения \(f_N(\mathbf{q}_1, \mathbf{p}_1, \dots, \mathbf{q}_N, \mathbf{p}_N, t)\), нормированная на единицу. Эволюция системы определяется гамильтонианом: \[ H = \sum_{i=1}^N \frac{\mathbf{p}_i^2}{2m} + \sum_{i<j} \Phi(|\mathbf{q}_i - \mathbf{q}_j|), \] где \(\Phi\) — потенциал парного взаимодействия. Уравнение Лиувилля для \(f_N\) имеет вид: \[ \frac{\partial f_N}{\partial t} + \{f_N, H\} = 0, \] где \(\{ \cdot , \cdot \}\) — скобки Пуассона.

Интегрируя уравнение Лиувилля по части переменных, получают уравнения для \(s\)-частичных функций распределения \(f_s\) (\(s = 1, 2, \dots, N-1\)). В общем виде уравнение для \(f_s\) содержит интеграл от \(f_{s+1}\), что и образует цепочку: \[ \frac{\partial f_s}{\partial t} + \sum_{i=1}^s \frac{\mathbf{p}_i}{m} \cdot \frac{\partial f_s}{\partial \mathbf{q}_i} - \sum_{i<j}^s \frac{\partial \Phi_{ij}}{\partial \mathbf{q}_i} \cdot \frac{\partial f_s}{\partial \mathbf{p}_i} = \sum_{i=1}^s \int \frac{\partial \Phi_{i,s+1}}{\partial \mathbf{q}_i} \cdot \frac{\partial f_{s+1}}{\partial \mathbf{p}_i} \, d\mathbf{q}_{s+1} d\mathbf{p}_{s+1}. \]

Здесь левая часть описывает свободное движение и взаимодействие между частицами внутри группы из \(s\) частиц, а правая — влияние внешних частиц через парные силы. Для \(s = 1\) уравнение принимает вид: \[ \frac{\partial f_1}{\partial t} + \frac{\mathbf{p}_1}{m} \cdot \frac{\partial f_1}{\partial \mathbf{q}_1} = \int \frac{\partial \Phi_{12}}{\partial \mathbf{q}_1} \cdot \frac{\partial f_2}{\partial \mathbf{p}_1} \, d\mathbf{q}_2 d\mathbf{p}_2, \] что соответствует кинетическому уравнению с интегралом столкновений, выраженным через двухчастичную функцию.

Квантовый случай

В квантовой статистической механике цепочка ББГКИ формулируется для матриц плотности. Для системы из \(N\) частиц вводится матрица плотности \(\rho_N(t)\), удовлетворяющая уравнению фон Неймана: \[ i\hbar \frac{\partial \rho_N}{\partial t} = [H, \rho_N]. \] Путём взятия частичного следа по переменным частиц получают уравнения для \(s\)-частичных матриц плотности \(\rho_s\). Структура цепочки аналогична классической, но содержит коммутаторы вместо скобок Пуассона.

Методы расцепления

Цепочка ББГКИ является незамкнутой — для нахождения \(f_s\) требуется знание \(f_{s+1}\). Для практического применения её необходимо «расцепить», то есть выразить \(f_{s+1}\) через функции меньшего порядка. Основные подходы:

• Приближение молекулярного хаоса (Stosszahlansatz): предполагает, что перед столкновением частицы статистически независимы: \(f_2(\mathbf{q}_1, \mathbf{p}_1, \mathbf{q}_2, \mathbf{p}_2, t) \approx f_1(\mathbf{q}_1, \mathbf{p}_1, t) f_1(\mathbf{q}_2, \mathbf{p}_2, t)\). Это приводит к уравнению Больцмана для разреженных газов.
• Метод Боголюбова: основан на разложении по малому параметру — плотности или силе взаимодействия. В первом порядке получается уравнение Больцмана, во втором — уравнения с учётом тройных корреляций (например, уравнение Чепмена — Энскога).
• Приближение среднего поля: для систем с дальнодействием (например, плазма) корреляции учитываются через самосогласованное поле, что приводит к уравнению Власова.
• Проекционные операторы: метод Мори (1965) и Цванцига (1960) позволяет формально замкнуть цепочку, вводя память и случайные силы.

Применение

Цепочка ББГКИ лежит в основе многих областей физики:

• Кинетическая теория газов: вывод уравнения Больцмана, уравнений Навье — Стокса и гидродинамических уравнений из микроскопической динамики.
• Физика плазмы: обоснование уравнения Власова для бесстолкновительной плазмы и уравнений Ландау для столкновительной.
• Теория жидкостей: описание корреляционных функций в равновесных и неравновесных жидкостях, в частности, через уравнения Орнштейна — Цернике.
• Фазовые переходы: изучение критических явлений и образования кластеров в конденсированных средах.
• Химическая физика: моделирование химических реакций в растворах и биологических системах.

Критика и ограничения

Основная трудность применения цепочки ББГКИ — необходимость расцепления, которое вносит приближения, не всегда обоснованные для плотных сред или систем с сильными корреляциями. В квантовом случае цепочка становится ещё более сложной из-за некоммутативности операторов и принципа тождественности частиц. Кроме того, для систем с дальнодействующими потенциалами (например, гравитация) разложение по парным взаимодействиям может быть некорректным. В современных исследованиях для преодоления этих ограничений используются численные методы (молекулярная динамика) и теория функционала плотности.

Источники

• Боголюбов Н. Н. Проблемы динамической теории в статистической физике. — М.: Гостехиздат, 1946.
• Born M., Green H. S. A General Kinetic Theory of Liquids. — Cambridge University Press, 1949.
• Kirkwood J. G. The Statistical Mechanical Theory of Transport Processes // Journal of Chemical Physics. — 1946. — Vol. 14, № 3. — P. 180–201.
• Yvon J. La théorie statistique des fluides et l'équation d'état. — Actualités scientifiques et industrielles, 1935.
• Ландау Л. Д., Лифшиц Е. М. Статистическая физика. Часть 1. — М.: Наука, 1976.
• Зубарев Д. Н. Неравновесная статистическая термодинамика. — М.: Наука, 1971.