On the mathematical divergences emerging in the theory of critical phenomena within Boltzmann–Gibbs statistical mechanics

The standard Theory of Critical Phenomena within Boltzmann–Gibbs statistical mechanics is known to successfully predict, at any temperature differing from a continuous-phase-transition critical one, the (finite) values for thermodynamically relevant quantities such as isothermal magnetic susceptibility, specific heat, correlation length, Grüneisen parameter, and similar ones. However, at the precise critical point, empirically unreachable infinities emerge within that approach at the N→∞ thermodynamical limit. Such divergences, although mathematically correct, may be seen as poorly informative ones. Indeed, they make no distinction at all between say the Ising, XY, Heisenberg, and other models. Moreover, for say the XY model, they make no distinction whether the spins are, for instance, 1/2, 1 or 3/2, or whether the interactions are say first-neighbor or first-and-second-neighbor ones. The introduction of nonadditive entropic functionals within a consistently generalized statistical mechanics substantially improves the situation, i.e., it regularizes the theory providing informative finite values instead of mere infinities. This undoubtedly enriching fact is illustrated Soares et al., Phys. Rev. B 111, L060409 (2025) for the Grüneisen parameter at the zero-temperature quantum critical point of a spin-1/2 ferromagnetic Ising chain in the presence of a transverse magnetic field, by introducing the nonadditive entropic functional Sq with the unique value q=q⋆≡37−6=0.0828…, which precisely guarantees the thermodynamical extensivity of the entropy as mandated by the preservation of the Legendre-transform structure of classical thermodynamics.