When Normality Tests Detect Equilibrium Distributions of Finite N-Body Systems
Jae Wan Shim
TL;DR
The paper addresses finite-$N$ deviations from Gaussianity in a one-dimensional, ideal-gas–like system by deriving a compact-support $q$-Gaussian via Havrda–Charvát entropy and evaluating five standard normality tests under a custom null across a wide $(N,n)$ grid. Using a large-scale Monte Carlo study, it shows that KS and AD maintain correct size under the custom null but have limited power to reject normality at larger $N$, while moment-based tests JB and SW exhibit much higher power, especially for smaller sample sizes, effectively detecting finite-size effects through skewness and kurtosis. As $N$ increases toward the Gaussian limit, all tests lose power, with ECDF-based tests requiring tens of thousands of observations to reach high power, whereas JB/SW reach high power with far fewer samples. The results provide practical guidance on test choice and required sample sizes for detecting finite-$N$ deviations in equilibrium distributions, with implications for physics and related fields where non-Gaussian finite-size effects are relevant.
Abstract
The particle number $N$ can be used as a quantitative gauge of non-Gaussianity. This idea extends to systems that are not literally finite by assigning them a notional $N $that captures the same deviation. For an ideal gas with $N$ insufficiently large for the thermodynamic limit, the velocity distribution that maximises Havrda-Charvát entropy departs markedly from the Maxwell-Boltzmann (Gaussian) form obtained in that limit. We explore how five standard normality tests-Kolmogorov-Smirnov, Anderson-Darling, Cramér-von Mises, Jarque-Bera and Shapiro-Wilk-respond to samples drawn from this finite-$N$ equilibrium distribution. A large-scale Monte-Carlo study maps the tests' statistical power across system size $N$ and sample size $n$, providing practical reference tables for deciding when finite-size effects remain detectable.
