Какое утверждение неверно? Что
![$\delta^{(3N)}(r - r_0(t))\delta^{(3N)}(p - p_0(t))$ $\delta^{(3N)}(r - r_0(t))\delta^{(3N)}(p - p_0(t))$](https://dxdy-04.korotkov.co.uk/f/3/c/4/3c4446387517a4104f1487843ba839fe82.png)
подчиняется уравнению Лиувилля? Перечитайте внимательнее.
Попробую объяснить.
(Оффтоп)
Жалко, что Cos(x-pi/2) спугнули. У него такие объяснения блестяще получаются.
1. О связи уравнения Лиувилля и динамики системы.
Пусть у нас есть
![$N$ $N$](https://dxdy-04.korotkov.co.uk/f/f/9/c/f9c4988898e7f532b9f826a75014ed3c82.png)
-частичная гамильтонова система. Пусть
![$X$ $X$](https://dxdy-01.korotkov.co.uk/f/c/b/f/cbfb1b2a33b28eab8a3e59464768e81082.png)
означает точку в фазовом пространстве размерности
![$2N,$ $2N,$](https://dxdy-02.korotkov.co.uk/f/5/e/4/5e4b27e7679e810261a2ad163a0a874782.png)
определяющую текущее состояние системы
![$X\to (p_1;q_1,p_2;q_2,\dots,p_N;q_N)$ $X\to (p_1;q_1,p_2;q_2,\dots,p_N;q_N)$](https://dxdy-01.korotkov.co.uk/f/4/8/f/48f148753032e0e44f1016da15921c2982.png)
.
В начальный момент значение функции
![$f(X_0,t=0)=W$ $f(X_0,t=0)=W$](https://dxdy-01.korotkov.co.uk/f/c/4/b/c4b8a3a285a34ffca8f6c996242045ea82.png)
в фиксированной точке
![$X_0$ $X_0$](https://dxdy-01.korotkov.co.uk/f/0/7/4/07478cd102054dc58a97f6fd8df8470582.png)
описывает вероятность найти систему в состоянии
![$X_0.$ $X_0.$](https://dxdy-03.korotkov.co.uk/f/6/3/5/635761736e5904d87d9838118db014ad82.png)
Тогда через время
![$t$ $t$](https://dxdy-01.korotkov.co.uk/f/4/f/4/4f4f4e395762a3af4575de74c019ebb582.png)
значение
![$W$ $W$](https://dxdy-01.korotkov.co.uk/f/8/4/c/84c95f91a742c9ceb460a83f9b5090bf82.png)
переползет в точку
![$X_0(t),$ $X_0(t),$](https://dxdy-02.korotkov.co.uk/f/1/7/a/17a8866a3ac0c3b3a92e580e4f2ca01482.png)
получающуюся из точки
![$X_0$ $X_0$](https://dxdy-01.korotkov.co.uk/f/0/7/4/07478cd102054dc58a97f6fd8df8470582.png)
сдвигом всех
![$p_i$ $p_i$](https://dxdy-01.korotkov.co.uk/f/0/d/1/0d19b0a4827a28ecffa01dfedf5f5f2c82.png)
и
![$q_i$ $q_i$](https://dxdy-02.korotkov.co.uk/f/9/2/9/9294da67e8fbc8ee3f1ac635fc79c89382.png)
вдоль траектории:
![$W=f(X_0(t)).$ $W=f(X_0(t)).$](https://dxdy-03.korotkov.co.uk/f/6/3/0/6304a8981edb4d87a29c87aa76ecf3b982.png)
Если эту процедуру проделать для всех
![$X,$ $X,$](https://dxdy-02.korotkov.co.uk/f/9/f/6/9f69f5d2d87593dc3e07849c256d40a982.png)
то получится, что
где
![$X(t)$ $X(t)$](https://dxdy-04.korotkov.co.uk/f/7/3/9/739fe69d785eb129eeb2746580e6437782.png)
- истинные механические траектории. Т.е. знание
![$f(X,t)$ $f(X,t)$](https://dxdy-02.korotkov.co.uk/f/9/b/7/9b7aeac03a5e7db9ce88e60091f8759382.png)
эквивалентно знанию всех траекторий (соединяем близкие точки с одинаковым
![$f$ $f$](https://dxdy-02.korotkov.co.uk/f/1/9/0/190083ef7a1625fbc75f243cffb9c96d82.png)
и получается траектория). Условием (необходимым и достаточным) того, что все так и будет является то, что
![$f$ $f$](https://dxdy-02.korotkov.co.uk/f/1/9/0/190083ef7a1625fbc75f243cffb9c96d82.png)
удовлетворяет уравнению Лиувилля для гамильтониана системы
![$H.$ $H.$](https://dxdy-03.korotkov.co.uk/f/6/2/0/620cf76d486c2bc6b701fc59aed46fc282.png)
В этом случае механические траектории будут характеристиками уравнения Лиувилля. Если уравнение, которому подчиняется
![$N$ $N$](https://dxdy-04.korotkov.co.uk/f/f/9/c/f9c4988898e7f532b9f826a75014ed3c82.png)
-частичная (полная) функция распределения отличается от уравнения Лиувилля, то бабушка надвое сказала, что вообще существует хоть какая механическая система, соответствующая такому уравнению.
2. О "малочастичных" функциях распределения.
Поскольку
![$N$ $N$](https://dxdy-04.korotkov.co.uk/f/f/9/c/f9c4988898e7f532b9f826a75014ed3c82.png)
-частичная функция содержит всё, толку от нее - ноль, мы ее никогда не узнаем. В большинстве случаев вместо нее используется "малочастичная" (одно- или двухчастичная) функция распределения. Она позволяет для некоторых величин, имеющих вид
получить средние,
совпадающие с получаемыми усреднением по полной функции распределения. Например, плотность в точке
![$r\;n(r,t)$ $r\;n(r,t)$](https://dxdy-03.korotkov.co.uk/f/2/d/f/2df3112e5b4930036ffa277e5dfecc5b82.png)
соответствует усреднению величины
![$\sum\limits_{1}^{N}\delta(q_i-r)$ $\sum\limits_{1}^{N}\delta(q_i-r)$](https://dxdy-04.korotkov.co.uk/f/f/f/d/ffde42b42eedc5ccac0dd099fb613ec882.png)
(векторы не ставлю, бо лень). Для одинаковых частиц функция распределения штука симметричная к перестановке аргументов. Тогда
![$f(p,q,t)=N\int f(X,t)dx_2\dots dx_N.$ $f(p,q,t)=N\int f(X,t)dx_2\dots dx_N.$](https://dxdy-02.korotkov.co.uk/f/d/5/a/d5afd04874017bb2e3c5811e0bd7996382.png)
Из этой формулы следует постоянство энтропии для гамильтоновых систем. Если одночастичную функцию определить как-нибудь по-другому, то средние по такой "новой" одночастичной и по "настоящей" функциям не совпадут.
3. О динамике "малочастичных" функций и уравнении Больцмана.
Поскольку одночастичная функция получается урезанием полной, то никакому замкнутому уравнению такая функция не удовлетворяет. Максимум, что можно получить - это уравнение, в левой части которого будет одночастичная функция, а в правой - двухчастичная. Этот процесс можно продолжить, и получить систему
![$N$ $N$](https://dxdy-04.korotkov.co.uk/f/f/9/c/f9c4988898e7f532b9f826a75014ed3c82.png)
(по числу частиц) зацепляющихся уравнений. Такая система будет эквивалентна уравнению Лиувилля для полной функции распределения. Если мы эту систему как-нибудь изменим, то изменится уравнение для полной функции распределения и соответствие статистической и обычной механики нарушится.
При выводе уравнения Больцмана двухчастичную функцию объявляют произведением одночастичных, обрывая цепочку то ли на первом, то ли на втором уравнении - смотря как считать. После этого восстановление механических уравнений по получившемуся аналогу уравнений Лиувилля дает результат, сильно отличный от исходной гамильтоновой смстемы. При этом уравнение Больцмана обладает забавным свойством - его правая часть обращается в ноль для равновесной функции распределения, а левая является уравнением Лиувилля для некоторой
фиктивной одночастичной механической системы. Поэтому благодаря правой части любое распределение эволюционирует к равновесному, после чего уравнение переходит в уравнение Лиувилля, сохраняющее во времени все, что надо.
Выводы.
Никакими легальными способами, не подменяющими исходную механическую гамильтонову систему чем-то другим невозможно получить термодинамические законы. Переход к квантовой механике положение не спасает, только вместо слов "функция распределения" надо произносить "матрица плотности", а все остальное будет тоже самое. Если у нас не получилось, что
![$\frac{dS}{dt}=0,$ $\frac{dS}{dt}=0,$](https://dxdy-03.korotkov.co.uk/f/a/f/c/afcc147d2404a0338f1b02aef3b09b3b82.png)
значит в каком-то месте мы подменили гамильтонову (Гайзенберговскую) систему чем-то ещё.