Уравнение Больцмана и H-теорема в функциональной формулировке классической механики Текст научной статьи по специальности «Математика»

УДК 517.958

УРАВНЕНИЕ БОЛЬЦМАНА И H-ТЕОРЕМА В ФУНКЦИОНАЛЬНОЙ ФОРМУЛИРОВКЕ КЛАССИЧЕСКОЙ МЕХАНИКИ

А. С. Трушечкин1,2

1 Математический институт им. В. А. Стеклова РАН,

119991, Москва, ул. Губкина, 8.

2 Национальный исследовательский ядерный университет «МИФИ»,

115409, Москва, Каширское ш., 31.

E-mails: [email protected]

Предлагается процедура получения уравнения Больцмана из уравнения Лиувил-ля в пределе, отличном от термодинамического. Она основывается на цепочках Боголюбова, функциональной формулировке классической механики и различении двух масштабов пространства-времени — макро- и микроскопического.

В соответствии с функциональным подходом к механике начальное состояние системы частиц формируется на основе измерений, которые имеют погрешности. Следовательно, можно говорить о точности, с которой задана начальная функция плотности вероятности в уравнении Лиувилля. Допустим, измерительные приборы прослеживают изменения физических величин лишь на макромасштабе, много большем, чем характерный радиус взаимодействия частиц (микромасштаб). Тогда соответствующую начальную функцию плотности нельзя использовать в качестве начального данного для уравнения Лиувилля, поскольку последнее представляет собой описание динамики на микромасштабе и в него явно входит потенциал взаимодействия между частицами (с характерным радиусом взаимодействия). Тем не менее для макроскопической начальной функции плотности можно получить уравнение Больцмана, воспользовавшись уравнением Лиувилля и идеологией цепочек Боголюбова, если предположить, что начальные условия для микроскопических функций плотности задаются макроскопической функцией. Показано, что для полученного уравнения верна H-теорема о возрастании энтропии.

Ключевые слова: статистическая механика, физическая кинетика, уравнение Больцмана, уравнение Лиувилля, цепочка уравнений Боголюбова.

Введение. Работа посвящена проблеме вывода кинетического уравнения Больцмана из уравнений микроскопической динамики (уравнения Лиувилля). Основной интерес заключается в том, что уравнение Лиувилля обратимо по времени (симметрично относительно замены t на — t), в то время как уравнение Больцмана необратимо. Напротив, для него справедлива так называемая H-теорема о возрастании энтропии. Проблема согласования необратимых по времени уравнений макроскопической динамики (Больцмана, Навье— Стокса и др.) с обратимыми уравнениями классической и квантовой микроскопической динамики (Ньютона, Гамильтона, Лиувилля, Шрёдингера) носит название проблемы необратимости [1].

Изящный вариант вывода уравнения Больцмана из уравнения Лиувилля предложен Н. Н. Боголюбовым [2], который ввёл иерархию времён. Но данный вывод приводит к расходимостям в членах разложения по степеням плотности более высокого порядка [3], что требует дальнейших исследований.

Антон Сергеевич Трушечкин (к.ф.-м.н.), научный сотрудник, отд. математической физики1; доцент, каф. системного анализа2.

В данной работе представлен несколько иной вариант вывода уравнения Больцмана из уравнения Лиувилля. Идеи метода цепочек Боголюбова здесь дополняются идеями недавно предложенной в работах И. В. Воловича функциональной формулировки механики [4], а также идеей различения двух масштабов пространства-времени — микро- и макроскопического (кинетического).

Другая попытка вывести уравнение типа Боголюбова—Больцмана в рамках функционального подхода к механике предпринята в работе [5]. Но полученное в этой работе кинетическое уравнение для системы из двух частиц по-прежнему (как и уравнение Лиувилля, из которого оно выведено) является обратимым. В настоящей работе получено необратимое уравнение Больцмана, для которого верна Н-теорема.

1. Уравнение Лиувилля. Пусть даны N частиц, находящихся в области

0 С К3 конечного объёма V. Их состояние в произвольный момент времени І описывается функцией / (х\, Х2, . . . , ,І), где Хі = (ді,Рі), Яі Є О (координата), рг Є К3 (импульс), і = 1, 2, ...,N, причём уУ / имеет смысл плотности вероятности нахождения системы частиц около соответствующей фазовой точки. Динамика функции / задаётся уравнением Лиувилля

I = {Н/).

где

N 2 N , | N

Н = £ 2т + £ *(^) + £ и(Яі)

і= 1 і,3 = 1 і=1

і>3

— гамильтониан системы, т > 0 — масса одной частицы, {■, ■} — скобка Пуассона. Первое слагаемое в Н соответствует кинетической энергии системы, Ф(г) —потенциал взаимодействия частиц, и(ді) — потенциал взаимодействия частицы с границей области О (стенкой сосуда).

Полагаем, что функция Ф(г) непрерывно дифференцируема, ограничена снизу, Ф(г) ^ +ж при г ^ 0, Ф(г) ^ 0 при г ^ ж. Также для простоты будем предполагать, что Ф(г) монотонно убывает с ростом г, а функция и непрерывно дифференцируема, и(д) ^ +ж при д ^ дО, где дО — граница области О.

Далее, у > 0 — малый параметр, который имеет смысл отношения характерного радиуса взаимодействия частиц г0 к характерному расстоянию

1 изменения одночастичной функции распределения /т, т. е. у = г0/1. Предел у ^ 0 соответствует стремлению к нулю характерного радиуса взаимодействия частиц в сравнении с характерным расстоянием изменения одночастичной функции распределения (или, эквивалентно, стремлению характерного расстояния изменения /1 к бесконечности в сравнении с характерным радиусом взаимодействия частиц).

Следуя Боголюбову, определим частичные функции распределения /3, в = 1, 2,..., N - 1:

/э(Х1, . . . ,Хз,г) = VN 3 / (Х1 ,...,ХN ,І)(ІХа+1 ...dХN,

■УпУ—

где Оу = С х М3 — фазовое пространство одной частицы, /м = /.

2. Микро- и макроскопический масштабы. Вывод Боголюбова уравнения Больцмана начинается с задачи Коши для уравнения Лиувилля:

Здесь необходимо вспомнить, что начальная функция распределения /0 не задана «объективно», а формируется на основе измерений [6,7]. Оказывается, существенным здесь является то, с какой степенью точности проведены эти измерения. Если они позволяют проследить изменение физических величин на масштабах порядка радиуса взаимодействия, то решение задачи Коши для уравнения Лиувилля теоретически возможно.

Нас интересует случай, когда измерения позволяют прослеживать изменение физических величин (и начальной функции распределения) лишь на масштабах, превышающих радиус взаимодействия частиц (и этот случай почти всегда и имеет место на практике, когда мы находимся на кинетическом уровне рассмотрения) [2,8]. Условие малости параметра ц = т0/1 означает то, что чувствительность наших приборов намного грубее, чем радиус взаимодействия. Характерное расстояние изменения одночастичной функции распределения I как раз непосредственно связано с чувствительностью прибора, можно сказать, I — характерное расстояние, на котором прибор улавливает изменения функции распределения. Масштаб порядка радиуса взаимодействия го будем называть «микромасштабом», масштаб порядка характерного изменения одночастичной функции распределения I — «макромасштабом» (точнее было бы сказать «кинетическом масштабом»). Следовательно, мы не знаем N-частичную функцию распределения /0(х1,... ,Хм), поскольку она должна содержать информацию о корреляциях частиц на масштабе порядка радиуса их взаимодействия. Причём эта информация является существенной. Например, когда расстояние между всеми частицами большое (|д — ^ | ^ г0 при всех г = ]), корреляции между ними могут отсутствовать, т. е. выполнено свойство

г=1

Именно это свойство и установит прибор, т. к. он прослеживает изменение функций распределения лишь на масштабах, много превышающих Г0. Однако

(вследствие того, что Ф(г) ^ при г ^ 0).

Итак, на масштабе порядка радиуса взаимодействия корреляции между частицами становятся существенными, но прибор не может их установить. Следовательно, мы не знаем начальной многочастичной функции распределения /0 и задача Коши (1) не имеет непосредственного физического смысла.

/ (Х1, ...,хм, 0) = / (Х1, .. .,хм),

(1)

N

/0 (Х1,Х2,...,Х3, 0) = Л /1^, 0).

если мы распространим это свойство на всё фазовое пространство , то система будет обладать бесконечной средней энергией

3. Постановка задачи для уравнения Лиувилля. Прибор позволяет нам установить начальную одночастичную функцию распределения /0(х1). Можно сказать, что /10 — «макроскопическое» распределение вероятностей, поскольку меняется на расстояниях, много превышающих радиус взаимодействия частиц. Будем предполагать, что /10 является непрерывной вместе со своими производными по каждому аргументу.

Определим другую задачу для уравнения Лиувилля для семейства функций /м = /м(х1, ..., ХМ) в пределе ц — 0, N — то, Nц2 — 0:

' ^ = {Н/м},

/1(Х1,0) — Л0(х1),

< ^-^4[/2(х1,х2,г — Дг) — /1(х1 ,г — Дг)Д(х2,г — Дг)} — 0, (2)

I Н (х1,...,хм )/м (Х1,...,ХМ ,г) (1х1 ...(1хм < то,

Vм

, /м (х1,...,хм ,г) = /м (х^1) ,...,хцм ),г), Уг е Р (1,..., N),

где Дг — 0, Дг/ц — то в рассматриваемом пределе, 8 > 0. Здесь исполь-(2)

зовано обозначение ;(Х1,Х2) —двухчастичный гамильтонов поток, опреде-

лённый на всех фазовых точках (Х1,Х2) таких, что Ц1 = д2. Таким образом, (2)

)^(х1,х2) = ^(х1*,х2*) (р — произвольная функция), где фазовая точка (хи,х21) получена сдвигом фазовой точки (х1,х2) вдоль гамильтонова потока

Н = ± 2т+ф( ^) + £ и («)

.= 1 ^ .= 1

на г. Через Р(1,..., N) обозначено множество перестановок натуральных чисел от 1 до N.

Отметим, что здесь рассматривается не термодинамический предел N—то, V ——то), N/V = еоп81, а предел бесконечного числа частиц в конечном объёме при бесконечно малом радиусе взаимодействия.

Третье условие в (2) говорит об ослаблении корреляций на бесконечно далёких от каждого столкновения временах. Оно является аналогом начального условия для двухчастичной функции распределения. Начальные условия для функций распределения большего числа частиц для вывода уравнения Больцмана не нужны, т. к. последнее учитывает только парные взаимодействия. Они нужны при выводе поправок к уравнению Больцмана, которые мы в данной работе не рассматриваем. Четвёртое условие в (2) —условие конечности энергии, пятое условие — условие симметричности многочастичных функций распределения относительно перестановок частиц.

Поясним условия Дг — 0 и Дг/ц — то при ц — 0 в третьем соотношении (2). Здесь Дг — промежуток времени в макроскопической шкале (порядка 1/й, где и — средняя скорость движения молекул), Дг/ц — в микроскопической (порядка г0/и, параметр ц = г0/1, напомним, и является масштабирующим множителем, связывающим эти шкалы). Другими словами, эти условия говорят о том, что рассматривается промежуток времени, который в микроскопической шкале в пределе бесконечно велик, а в макроскопической —

бесконечно мал. Образно говоря, то, что на макроскопическом временном масштабе времени длится «мгновение», на микроскопическом — «бесконечно долго». Аналогично с соответствующими областями пространства (шарами радиуса Дги): что на макромасштабе — «точка», то на микромасштабе — «бесконечное пространство». Таким образом, в нашем выводе явно присутствуют два различных масштаба пространства-времени.

Уравнение Лиувилля есть не что иное, как уравнение динамики на микромасштабе, поскольку в него явно входит потенциал взаимодействия между частицами с характерным радиусом г>. А в уравнении Больцмана одновременно участвуют микро- и макроскопический масштабы. С одной стороны, столкновение в уравнении Больцмана рассматривается как мгновенное и точечное, с другой — в нём участвует сечение рассеяния, которое определяется как асимптотика при времени, стремящемся к бесконечности. В первом случае (при рассмотрении динамики одночастичной функции распределения /1) мы имеем дело с макроскопическим масштабом, во втором (при рассмотрении динамики столкновения отдельных частиц) — с микроскопическим.

4. Основная теорема.

Теорема. Пусть семейство функций /м = /м (х1,..., хм ,г) удовлетворяет задаче (2). Пусть дополнительно Ф(г) в уравнении Лиувилля при г — то имеет асимптотику

Нш г1 Ф(г) = С = 0,

Г—

где ^ > 2. Тогда функция /1(х1,г) удовлетворяет уравнению д/1 Р1 д/1 , дид/1 , , (1'

дг т дд1 дд1 др1

Бг/1 = 1 / ——— [/1(я1,р[,г)/1(д1,р,2,г) — /1(д1,Р1 ,г)/1(дьр2,г)]йафь V,] в т

(4)

при ц — 0, N — то, Nц2 — 0. Здесь V = V/(Nц2); Б = М2 х М3; р\ и Р2 — импульсы после столкновения двух частиц, которые до столкновения имели импульсы р1 и р2 соответственно; йа = рйрйр —дифференциальное сечение рассеяния; р > 0 — прицельный параметр; р е [0, 2п) — полярный угол; р1 и р'2 являются функциями от р1, р2 и р; зависимость полностью определяется двухчастичным гамильтонианом.

Соотношение (3) представляет собой не что иное, как уравнение Больцмана, для которого имеет место Н-теорема:

йБ (г) _ й йг йг

/ / (х, г) 1п / (х,г) йх

1Пу

= А +о(1)

(для тех функций, для которых интеграл сходится), где А ^ 0, А = 0(1), при ц — 0, N — то, Nц2 — 0.

Заключение. Итак, мы получили уравнение Больцмана из уравнения Лиу-вилля в конечном объёме. Мы воспользовались цепочками Боголюбова, функциональной формулировкой механики и различением двух масштабов пространства-времени: микро- и макроскопического. Главное отличие нашего вывода от вывода Боголюбова состоит в том, что вместо задачи Коши (1) мы

исходим из задачи (2) для уравнения Лиувилля. В постановке этой задачи нет начальной многочастичной функции распределения, поскольку она, как правило, не может быть установлена измерительными приборами. Вместо этого начальное условие для многочастичной («микроскопической») функции распределения определяется из условия ослабления корреляций, т. е. на основе одночастичной («макроскопической») функции.

В теоретической и математической физике принята картина, согласно которой динамика системы полностью определена, если заданы начальные условия. Однако сами начальные условия понимаются как нечто внешнее по отношению к уравнениям математической физики («Законы природы хранят молчание относительно всего, что касается состояния мира в данный момент» [9]). Здесь мы предлагаем другую картину: начальные данные для рассматриваемого уровня природы задаются из более высокого уровня.

Автор благодарен И. В. Воловичу за постановку задачи и полезные замечания.

Работа выполнена при частичной финансовой поддержке Российского фонда фундаментальных исследований (проекты №№ 11-01-00828-а, 09-01-12161-офи-м), гранта Президента РФ (проект № НШ-7675.2010.1) и программы ОМН РАН.

Автор выражает благодарность за поддержку на конференции «Математическая физика и ее приложения - 2010» лабораторией математической физики СамГУ, грантами АВЦП 3341 и 10854 и контрактом ФЦП 2173.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Козлов В. В. Тепловое равновесие по Гиббсу и Пуанкаре. М., Ижевск: Институт компьютерных исследований, 2002. 320 с. [Kozlov V. V. Thermal equilibrium in the sense of Gibbs and Poincare. Moscow, Izhevsk: Institut Komp’yuternykh Issledovanij, 2002. 320 pp.]

2. Боголюбов Н. Н. Проблемы динамической теории в статистической физике. М., Л.: Гостехиздат, 1946. 119 с. [Bogolyubov N. N. Problems of Dynamical Theory in Statistical Physics. Moscow, Leningrad: Gostekhizdat, 1946. 119 pp.]

3. Боголюбов Н. Н. Кинетические уравнения и функции Грина в статистической механике / В сб.: Собрание научных трудов в двенадцати томах. Т. V. М.: Наука, 2005. С. 616638. [Bogolyubov N. N. Kinetic equations and Green functions in statistical mechanics / In: Collection of scientific works in twelve volumes. Vol. V. Moscow: Nauka, 2005. Pp. 616-638].

4. Волович И. В. Проблема необратимости и функциональная формулировка классической механики// Вестн. Сам. гос. ун-та. Естественнонаучн. сер., 2008. №8/1(67). С. 35-55, arXiv: 0907.2445 [cond-mat.stat-mech]. [Volovich I. V. Time Irreversibility Problem and Functional Formulation of Classical Mechanics // Vestn. Samar. Gos. Univ. Estest-vennonauchn. Ser., 2008. no. 8/1(67). Pp. 35-55].

5. Волович И. В. Уравнения Боголюбова и функциональная механика // ТМФ, 2010. Т. 164, №3. С. 354-362; англ. пер.: Volovich I. V. Bogoliubov equations and functional mechanics // Theoret. and Math. Phys., 2010. Vol. 164, no. 3. Pp. 1128-1135.

6. Trushechkin A. S., Volovich I. V. Functional classical mechanics and rational numbers // p-Adic Numbers, Ultrametric Analysis, and Applications, 2009. Т. 1, №4. С. 361-367, arXiv: 0910.1502 [math-ph].

7. Трушечкин А. С. Необратимость и роль измерительного прибора в функциональной формулировке классической механики// ТМФ, 2010. Т. 164, №3. С. 435-440; англ. пер.: Trushechkin A. S. Irreversibility and the role of an instrument in the functional formulation of classical mechanics// Theoret. and Math. Phys., 2010. Vol. 164, no. 3. Pp. 1198-1201.

8. Лифшиц Е.М., Питаевский Л. П. Физическая кинетика. М.: Физматлит, 2002. 536 с. [Lifshitz E. M., Pitaevskii L. M. Physical Kinetics. Moscow: Fizmatlit, 2002. 536 pp.]

9. Wigner E. The Unreasonable Effectiveness of Mathematics in the Natural Sciences // Commun. Pure Appl. Math., 1960. Vol. 13, no. 1. Pp. 1-14; русск. пер.: Вигнер Е. Непо-

стижимая эффективность математики в естественных науках / В сб.: Этюды о симметрии. М.: Мир, 1971. С. 182-198.

Поступила в редакцию 21/Х11/2010; в окончательном варианте — 21/11/2011.

MSC: 82C05, 82C40

BOLTZMANN EQUATION AND H-THEOREM IN THE FUNCTIONAL FORMULATION OF CLASSICAL MECHANICS

A. S. Trushechkin1,2

1 Steklov Mathematical Institute, Russian Academy of Sciences,

8, Gubkina st., Moscow, 119991, Russia.

2 National Research Nuclear University “MEPhI”,

31, Kashirskoe sh., Moscow, 115409, Russia.

E-mails: [email protected]

We propose a procedure for obtaining the Boltzmann equation from the Liouville equation in a non-thermodynamic limit. It is based on the BBGKY hierarchy, the functional formulation of classical mechanics, and the distinguishing between two scales of space-time, i.e., macro- and microscale. According to the functional approach to mechanics, a state of a system of particles is formed from the measurements, which have errors. Hence, one can speak about accuracy of the initial probability density function in the Liouville equation. Let’s assume that our measuring instruments can observe the variations of physical values only on the macroscale, which is much greater than the characteristic interaction radius (microscale). Then the corresponfing initial density function cannot be used as initial data for the Liouville equation, because the last one is a description of the microscopic dynamics, and the particle interaction potential (with the characteristic interaction radius) is contained in it explicitly. Nevertheless, for a macroscopic initial density function we can obtain the Boltzmann equation using the BBGKY hierarchy, if we assume that the initial data for the microscopic density functions are assigned by the macroscopic one. The H-theorem (entropy growth) is valid for the obtained equation.

Key words: statistical mechanics, physical kinetics, Boltzmann equation, Liouville equation, BBGKY hierarchy.

Original article submitted 21/XII/2010; revision submitted 21/II/2011.

Anton S. Trushechkin (Ph.D. (Phys. & Math.)), Researcher, Dept. of Mathematical Physics1; Associate Professor, Dept. of System Analysis2.