We derive the nonlinear Boltzmann equation from a grand-canonical Newtonian gas with initial pair exclusion. We treat stable, finite-range radial potentials that are C² except for an allowed repulsive singularity at the origin, and C¹ initial probability densities with spatially summable Gaussian bounds on the density and its spatial gradient. All fixed-order rescaled factorial marginals converge in L¹, uniformly throughout every finite interval on which the classical kinetic solution has a uniform Gaussian bound. The result allows attractive wells and dynamically formed clusters, without restrictions on the differential scattering cross-section.

We prove a finite-dimensional central limit theorem away from equilibrium for a deterministic grand-canonical hard-sphere gas in three dimensions. The initial probability density is smooth, with spatially summable Gaussian velocity bounds on the density and its spatial gradient. The limit holds on every finite interval on which the classical Boltzmann solution has uniform Gaussian velocity decay. The empirical measure is centered by its exact microscopic expectation. Its fluctuations converge to the Gaussian solution of the linear fluctuating Boltzmann equation.