Пусть

--- количество частиц (шариков),

--- размерность пространства. Ну, для наших целей, наверное, можно считать, что

.
Я использую то же самое пространство что и китайцы

,

Здесь

--- координата и скорость частицы с номером

,

На стр. 4 они дают формальное описание динамики шариков.
Уравнение (1.3) описывает свободное прямолинейное движение шариков между столкновениями.
Формулы (1.2) описывают скорости после удара.
Траекторию они называют "flow" и обозначают

в (1.4).
Точнее, они вводят "flow map"

.
Ну, это просто все траектории, как отображение начальных точек в конечные.
Далее они рассматривают все возможные начальные точки

и вводят на них вероятность по формуле (1.6), (1.7).
Точнее, у них там все немножко сложнее. Они рассматривают случайное количество частиц

и случайные начальные данные. Плюс, надо следить за тем, чтобы шарики не пересекались. Это накладывает некоторые ограничения на возможные значения

. Но для нас это все несущественно.
Главное, что они рассматривают ВСЕ ВОЗМОЖНЫЕ ТРАЕКТОРИИ из

частиц. Все они приходят в какие-то конечные точки

. На множестве этих конечных точек распределение вероятности на начальных данных индуцирует вероятность появления данной конечной точки в фазовом пространстве

по формуле (1.9).
Затем они рассматривают корреляционную функцию

--- совместное распределение

по формуле (1.11).
Например, при

величина

--- это в точности условная вероятность обнаружить первую частицу в состоянии

.
Теорема 1 утверждает, что эта величина удовлетворяет уравнению Больцмана (ну, понятное дело, приближенно, с оценкой погрешности итд).
Что здесь получается? Каждое состояние в

достижимо (с оговорками геометрического характера). Поэтому на всем этом пространстве имеется распределение вероятности. Вот по этой вероятности и усредняем. Получаем

. Дискретная динамика лишь инструмент для получения распределения вероятности.
А я хотел бы иметь дело с одной единственной траекторией. В пространстве

это одномерная линия. Вероятность обнаружить частицу в конкретном состоянии

имеет

-образный характер. Усредняя по маленькому шарику получаю гладкую среднюю плотность.
В сущности, вывод уравнения Больцмана как-то-так и выглядит. Рассматриваем маленький объемчик (шарик). За малое время

из него что-то вылетает, что-то в него залетает. Моделируем перераспределение по скоростям, составляем баланс итд. При этом нет усреднения по РАЗНЫМ начальным данным. Это все один и тот же процесс. При этом объемчик не такой уж и малый, чтобы пренебречь дискретностью, а все усреднение спрятано в столкновительный член. Сколько там чего столкнулось за это время и во что это все превратилось. Поэтому и время

не такое уж и малое, чтобы там набралась статистика столкновений. Разумеется, в этот момент и возникает вопрос об обосновании этой статистики. Почему детерминированная динамика допускает такое статистическое обоснование? Ну, да, мы видим, что процесс столкновений ИСКЛЮЧИТЕЛЬНО неустойчивый. В результате получается этакий датчик (псевдо)случайных чисел. Но на мой взгляд, мы имеем дело с двумя РАЗНЫМИ случайностями. В случае анасамбля начальных данных имеет место неустойчивость по начальным данным. А в случае одной траектории речь идет о неустойчивости и статистическом перемешивании в течение достаточно короткого промежутка времени

. Почему это так --- вопрос.