делим систему на физически бесконечно малые объёмы и считаем энтропию как сумму по ним
Я тут ещё подумал, посмотрел источники. Очень странно всё получается в присутствии дальнодействующих сил. Настолько странно, что возможно этот подход не применим. В принципе такое деление нужно чтобы энтропия была экстенсивной величиной даже в неравновесном случае, больше ни для чего. И хотя формально экстенсивность получается в присутствии дальнодействующих сил, оказывается, что локальная энтропия зависит не только от локального состояния, что по сути губит весь профит от экстенсивности.
В реальных работах (например,
тут) берут просто одночастичную функцию распределения и считают энтропию Гиббса по ней, без дополнительного разбиения на малые объёмы.
Нашёл кстати названия игрушечных моделей подобных систем, может пригодится: SGR (Self Gravitating Ring) и HMF (Hamiltonian mean-field). (Они очень похожи, я даже не понимаю в чём между ними разница и есть ли она.)