Quantitative range, short-range potentials, and off-equilibrium fluctuations: the open core of the long-time derivation of the Boltzmann equation from hard-sphere dynamics

In 2024-2025, Yu Deng, Zaher Hani, and Xiao Ma proved that Boltzmann's kinetic equation can be rigorously derived from the Newtonian dynamics of a rarefied hard-sphere gas for arbitrarily long times - any fixed time interval on which the limiting Boltzmann solution exists and stays bounded - extending Lanford's 1975 short-time validity theorem; together with a companion paper (arXiv:2503.01800) that adds the torus case and the passage to fluid equations, the authors claim completion of Hilbert's sixth problem for a rarefied hard-sphere gas, under strong regularity assumptions. Resolve the surviving open directions (A)-(C) below. Three linked directions remain open: (A) quantitative range - when the validity horizon is allowed to depend on the particle diameter $\varepsilon$, the current proof reaches only $t_{\rm fin}\sim\log|\log\varepsilon|$, the correlation estimates cover only cluster orders $s\leq|\log\varepsilon|$, and the convergence rate $\varepsilon^\theta$ is non-optimal (targets include $t_{\rm fin}\sim|\log\varepsilon|$, $s\leq\varepsilon^{-\kappa}$, and $\theta\to 1^-$; a known obstruction sits near $s\sim\varepsilon^{-1}$); (B) model generality - the extension to compactly supported short-range interaction potentials, where positive-duration and simultaneous multi-particle interactions must be controlled, is expected by the authors to be within reach but is unproved; (C) statistical depth - the fluctuating Boltzmann equation and the large-deviation principle of Bodineau-Gallagher-Saint-Raymond-Simonella hold only for short times, or for long times at exact equilibrium, and their long-time off-equilibrium derivation - linked by the authors to the optimal cumulant decay rates - is open.

Background

Consider a classical gas of $N$ identical hard spheres of diameter $\varepsilon$ moving in $\mathbb{R}^d$ ($d\geq 2$) under Newton's laws with elastic, energy- and momentum-conserving collisions. In the Boltzmann-Grad scaling one lets $N\to\infty$ and $\varepsilon\to 0$ with $N\varepsilon^{d-1}$ held fixed (a dilute gas), so that the mean free time - the average time between two successive collisions of a given particle - is of order one. Boltzmann proposed in 1872 that the one-particle distribution $f(t,x,v)$ then evolves by a closed equation, $(\partial_t + v\cdot\nabla_x)f = \iint_{\mathbb{R}^d\times S^{d-1}}\big((v-v_1)\cdot\omega\big)^+(f'f_1' - ff_1)\,d\omega\,dv_1,$ where $(v',v_1')$ are the post-collisional velocities of a pair with pre-collisional velocities $(v,v_1)$ and collision direction $\omega$. Closure rests on the Stosszahlansatz (molecular chaos): particles about to collide are statistically independent. Justifying this independence from Newtonian dynamics is the central problem of rigorous kinetic theory and constitutes the kinetic step of Hilbert's sixth problem (1900), which asks for the derivation of the equations of continuum mechanics from atomic mechanics. The exact dynamics is encoded in the BBGKY hierarchy - the chain of equations coupling the $s$-particle correlation functions (joint marginal densities) $f_s$ of the system; the failure of chaos is measured by the cumulants, i.e. the connected parts of the $f_s$ relative to products $\prod_{j=1}^s f(z_j)$. Lanford (1975) proved that, for random initial data with independent particles, the $f_s$ converge to products of the Boltzmann solution as $\varepsilon\to 0$, so the empirical measure of the gas converges to $f(t)$ - but only for a very short time, a small fraction of the mean free time, because his expansion of the hierarchy into collision trees converges only on such intervals. This short-time barrier stood for fifty years, with two kinds of exceptions: global-in-time validity for a near-vacuum gas cloud (Illner-Pulvirenti, where the collision rate is too small for correlations to build up), and extensions to compactly supported short-range interaction potentials (Gallagher-Saint-Raymond-Texier), still short-time. Pulvirenti-Simonella analyzed the correlation error at short times and proved decay rates of the cumulants proportional to the number of particles involved; in particular, for clusters of order $s\sim\varepsilon^{-1}$ the chaotic ansatz fails already at time zero. Beyond the law of large numbers, Bodineau-Gallagher-Saint-Raymond-Simonella developed the statistical dynamics of the hard-sphere gas in the same limit: central-limit-scale fluctuations of empirical observables around the Boltzmann approximation are governed by the fluctuating Boltzmann equation (a linearized kinetic equation with a forcing noise term), and probabilities of large deviations satisfy a functional Hamilton-Jacobi equation (an equation of Hamilton-Jacobi type for the rate functional, the functional measuring the cost of a fluctuation); these results hold for short times, or globally in time only at exact global Maxwellian equilibrium, where the linear structure can be exploited. In 2024-2025, Deng-Hani-Ma proved the long-time validity. Their Theorem 1: for the hard-sphere system started from a grand-canonical ensemble of independent, identically distributed particles (the particle number is Poisson-distributed at fixed density; positions and velocities are independent modulo the hard-core exclusion), if the Boltzmann solution exists on $[0,t_{\rm fin}]$ with bounds of Gaussian type in the velocities ($e^{\beta|v|^2}$-weighted sup-norms) and the initial datum has finite Boltzmann-type norms, then for $\varepsilon$ small enough, uniformly in $t\in[0,t_{\rm fin}]$ and in $s\leq|\log\varepsilon|$, $\Big\| f_s(t,z_s) - \prod_{j=1}^{s} f(t,z_j)\,\mathbbm{1}_{\mathcal{D}_s}(z_s)\Big\|_{L^1(\mathbb{R}^{2ds})} \leq \varepsilon^{\theta},$ where $\theta>0$ depends only on $d$ and $\mathbbm{1}_{\mathcal{D}_s}$ is the indicator of the hard-core domain $\mathcal{D}_s$ (configurations of $s$ particles with no overlapping spheres). The time $t_{\rm fin}$ is any fixed time within the life of the Boltzmann solution, so the validity time is as long as the limiting equation itself is well posed. The proof propagates a cumulant ansatz that keeps memory of the full collision history, via a partial time expansion; the smallness of the cumulants is reduced to combinatorial properties of 'molecules' (diagrams encoding the collision history of a cluster of particles), established by an elaborate cutting algorithm, with a uniform bound on the number of collisions in semi-dispersing billiards (billiard motion in a domain with convex scatterers, which is how the multi-particle hard-sphere dynamics is recast) due to Burago-Ferleger-Kononenko controlling the geometry. The strategy transplants the paradigm the first two authors developed for the wave kinetic equation of wave turbulence theory - the statistical kinetics of weakly nonlinear dispersive waves, whose limiting equation is the wave analog of the Boltzmann equation. A companion paper by the same authors extends the derivation to two- and three-dimensional tori and combines it with hydrodynamic-limit results to derive the compressible Euler and incompressible Navier-Stokes-Fourier equations from hard-sphere systems, claiming completion of Hilbert's sixth problem for a rarefied hard-sphere gas (under strong regularity assumptions, in an iterated limit: first Boltzmann-Grad, then a fast-relaxation limit in which collisions act much faster than transport). What remains open, and what this problem records, are the quantitative and model-generality limits of the long-time derivation itself, together with its extension from the law of large numbers to the fluctuation regime.

Problem Statement

Resolve the following three linked components, which are the surviving open directions that Deng, Hani, and Ma (arXiv:2408.07818, Remark 1.5 and Section 1.4) attach to their long-time validity theorem for the Boltzmann equation from hard-sphere dynamics in the Boltzmann-Grad limit. The setting of all components is that of their Theorem 1: the hard-sphere system of diameter $\varepsilon$ in $\mathbb{R}^d$ ($d\geq 2$), started from a grand-canonical ensemble of independent identically distributed particles under the Boltzmann-Grad scaling, with $f(t)$ the solution of the Boltzmann equation. (A) Quantitative range of validity. (A1) Time horizon: assuming the Boltzmann solution is global in time and uniformly bounded in the norms of Theorem 1, prove that the validity time can be taken to diverge as $\varepsilon\to 0$ faster than the currently available $t_{\rm fin}\sim\log|\log\varepsilon|$ (the authors suggest $t_{\rm fin}\sim|\log\varepsilon|$ as a target), or establish the optimal achievable rate by an upper bound. (A2) Cluster order: extend the correlation-defect estimate from $s\leq|\log\varepsilon|$ to $s\leq\varepsilon^{-\kappa}$ for some $\kappa>0$ (the authors state their method should reach this for sufficiently small $\kappa$ depending on $t_{\rm fin}$, but this is not carried out), and delimit the range of $s$ for which such an estimate can hold, in view of the expected failure near $s\sim\varepsilon^{-1}$ already at time zero. (A3) Rate: improve the convergence exponent $\theta$: for $s=1$, prove the expected optimal rate $\theta\to 1^-$; for $s\geq 2$, determine the optimal decay rate of the cumulants. (B) Short-range potentials. Extend the long-time validity theorem (validity for any fixed $t_{\rm fin}$ within the life of the Boltzmann solution, with the corresponding correlation-defect control) from hard spheres to particle systems interacting through compactly supported short-range potentials in the Boltzmann-Grad scaling, the class treated at short times by Gallagher-Saint-Raymond-Texier. This requires controlling interactions of positive duration and genuinely simultaneous interactions of more than two particles, for which the hard-sphere techniques (instantaneous collisions; Burago-Ferleger-Kononenko collision counting) have no direct analogue. (C) Long-time off-equilibrium fluctuations. Extend the fluctuating Boltzmann equation (central-limit-scale fluctuations of empirical observables around the Boltzmann approximation) and the associated large-deviation principle, proved by Bodineau-Gallagher-Saint-Raymond-Simonella only for short times, or for long times at exact global equilibrium, to the long-time, off-equilibrium regime; the authors state that the optimal cumulant decay rates of component (A3) for $s\geq 2$ are linked to the long-time derivation. The common objective of the three components is to determine the full domain of validity of the Boltzmann approximation of Newtonian dynamics in the Boltzmann-Grad limit: how long (A1), at which cluster order (A2) and precision (A3), for which interaction laws (B), and at which statistical order (C). Each component is a genuine open question as of September 2026, to be settled in its stated generality rather than in the strictly easier linearized, spatially homogeneous, or at-equilibrium special cases.

References

1. Yu Deng, Zaher Hani, Xiao Ma, 'Long time derivation of the Boltzmann equation from hard sphere dynamics', arXiv:2408.07818 (v1 August 2024, v3 July 2025), DOI 10.48550/arXiv.2408.07818, https://arxiv.org/abs/2408.07818
2. Yu Deng, Zaher Hani, Xiao Ma, 'Hilbert's sixth problem: derivation of fluid equations via Boltzmann's kinetic theory', arXiv:2503.01800 (March 2025), DOI 10.48550/arXiv.2503.01800, https://arxiv.org/abs/2503.01800
3. Thierry Bodineau, Isabelle Gallagher, Laure Saint-Raymond, Sergio Simonella, 'Derivation of the Boltzmann equation from hard-sphere dynamics (after Y. Deng, Z. Hani, and X. Ma)', Bourbaki Seminar, January 2026, arXiv:2602.04407, DOI 10.48550/arXiv.2602.04407, https://arxiv.org/abs/2602.04407
4. O. E. Lanford III, 'Time evolution of large classical systems', in Dynamical systems, theory and applications (Rencontres, Battelle Research Institute, Seattle, 1974), Lecture Notes in Physics 38, Springer, Berlin, 1975, pp. 1-111, DOI 10.1007/3-540-07171-7_1, https://doi.org/10.1007/3-540-07171-7_1
5. Isabelle Gallagher, Laure Saint-Raymond, Benjamin Texier, 'From Newton to Boltzmann: Hard Spheres and Short-range Potentials', Zurich Lectures in Advanced Mathematics 18, European Mathematical Society, 2014, DOI 10.4171/129, https://doi.org/10.4171/129
6. R. Illner, M. Pulvirenti, 'Global validity of the Boltzmann equation for two- and three-dimensional rare gas in vacuum: Erratum and improved result', Communications in Mathematical Physics 121 (1989), no. 1, 143-146, DOI 10.1007/BF01218628, https://doi.org/10.1007/BF01218628
7. M. Pulvirenti, S. Simonella, 'The Boltzmann-Grad limit of a hard sphere system: analysis of the correlation error', Inventiones mathematicae 207 (2017), no. 3, 1135-1237, DOI 10.1007/s00222-016-0682-4, https://doi.org/10.1007/s00222-016-0682-4
8. T. Bodineau, I. Gallagher, L. Saint-Raymond, S. Simonella, 'Statistical dynamics of a hard sphere gas: fluctuating Boltzmann equation and large deviations', Annals of Mathematics (2) 198 (2023), no. 3, 1047-1201, DOI 10.4007/annals.2023.198.3.3, https://doi.org/10.4007/annals.2023.198.3.3
9. T. Bodineau, I. Gallagher, L. Saint-Raymond, S. Simonella, 'Long-time derivation at equilibrium of the fluctuating Boltzmann equation', Annals of Probability 52 (2024), no. 1, 217-295, DOI 10.1214/23-AOP1656, https://doi.org/10.1214/23-aop1656
10. Nathalie Ayi, 'From Newton's Law to the Linear Boltzmann Equation Without Cut-Off', Communications in Mathematical Physics 350 (2017), no. 3, 1219-1274, DOI 10.1007/s00220-016-2821-6, https://doi.org/10.1007/s00220-016-2821-6
11. Corentin Le Bihan, 'Long Time Validity of the Linearized Boltzmann Equation for Hard Spheres: A Proof Without Billiard Theory', Archive for Rational Mechanics and Analysis 249 (2025), no. 4, Article 40, DOI 10.1007/s00205-025-02105-z, https://doi.org/10.1007/s00205-025-02105-z
12. Corentin Le Bihan, 'Long-Time Validity of the Linearized Boltzmann Uncut-Off and the Linearized Landau Equations From the Newton Law', Communications on Pure and Applied Mathematics 79 (2026), no. 7, 1673-1745, DOI 10.1002/cpa.70035, https://doi.org/10.1002/cpa.70035 (preprint arXiv:2408.03597)
13. Chenjiayue Qi, 'Global Solution of a Functional Hamilton-Jacobi Equation associated with a Hard Sphere Gas', arXiv:2409.02805 (2024), DOI 10.48550/arXiv.2409.02805, https://arxiv.org/abs/2409.02805
14. Shan Gao, 'Comment on "Hilbert's Sixth Problem: Derivation of Fluid Equations via Boltzmann's Kinetic Theory" by Deng, Hani, and Ma', arXiv:2504.06297 (2025), DOI 10.48550/arXiv.2504.06297, https://arxiv.org/abs/2504.06297
15. Ioakeim Ampatzoglou, Nataša Pavlović, William Warner, 'Derivation of the Higher Order Boltzmann Equation for Hard Spheres', arXiv:2409.05669 (2024), DOI 10.48550/arXiv.2409.05669, https://arxiv.org/abs/2409.05669
16. Yu Deng, Zaher Hani, 'Full derivation of the wave kinetic equation', Inventiones mathematicae 233 (2023), no. 2, 543-724, DOI 10.1007/s00222-023-01189-2, https://doi.org/10.1007/s00222-023-01189-2
17. David Hilbert, 'Mathematical problems' (lecture delivered to the International Congress of Mathematicians, Paris, 1900), Bulletin of the American Mathematical Society 37 (2000), 407-436 (reprint of the 1902 English translation), DOI 10.1090/s0273-0979-00-00881-8, https://doi.org/10.1090/s0273-0979-00-00881-8
18. D. Burago, S. Ferleger, A. Kononenko, 'Uniform Estimates on the Number of Collisions in Semi-Dispersing Billiards', Annals of Mathematics 147 (1998), no. 3, 695-708, DOI 10.2307/120962, https://doi.org/10.2307/120962