A recent paper by T. Dauxois entitled "Non-Gaussian distributions under scrutiny" is submitted to scrutiny. Several comments on its content are made, which constitute, at the same time, a brief state-of-the-art review of nonextensive statistical mechanics, a current generalization of the Boltzmann-Gibbs theory. Some inadvertences and misleading sentences are pointed out as well.

We give a closer look at the Central Limit Theorem (CLT) behavior in quasi-stationary states of the Hamiltonian Mean Field model, a paradigmatic one for long-range-interacting classical many-body systems. We present new calculations which show that, following their time evolution, we can observe and classify three kinds of long-standing quasi-stationary states (QSS) with different correlations. The frequency of occurrence of each class depends on the size of the system. The different microsocopic nature of the QSS leads to different dynamical correlations and therefore to different results for the observed CLT behavior.

We present a molecular dynamics test of the Central Limit Theorem (CLT) in a paradigmatic long-range-interacting many-body classical Hamiltonian system, the HMF model. We calculate sums of velocities at equidistant times along deterministic trajectories for different sizes and energy densities. We show that, when the system is in a chaotic regime (specifically, at thermal equilibrium), ergodicity is essentially verified, and the Pdfs of the sums appear to be Gaussians, consistently with the standard CLT. When the system is, instead, only weakly chaotic (specifically, along longstanding metastable Quasi-Stationary States), nonergodicity (i.e., discrepant ensemble and time averages) is observed, and robust $q$-Gaussian attractors emerge, consistently with recently proved generalizations of the CLT.

Recent measurements of the temperature fluctuations of the cosmic microwave background (CMB) radiation from the WMAP satellite provide indication of a non-Gaussian behavior. Although the observed feature is small, it is detectable and analyzable. Indeed, the temperature distribution P^{CMB}(Delta T) of these data can be quite well fitted by the anomalous probability distribution emerging within nonextensive statistical mechanics, based on the entropy S_q = k (1 - \int dx [P(x)]^q)/(q - 1) (where in the limit case q -> 1 we obtain the Boltzmann-Gibbs entropy S_1 = - k \int dx P(x) ln[P(x)]). For the CMB frequencies analysed, \nu= 40.7, 60.8, and 93.5 GHz, P^{CMB}(Delta T) is well described by P_q(Delta T) \propto 1/[1 + (q-1) B(\nu) (Delta T)^2]^{1/(q-1)}, with q = 1.04 \pm 0.01, the strongest non-Gaussian contribution coming from the South-East sector of the celestial sphere. Moreover, Monte Carlo simulations exclude, at the 99% confidence level, P_1(Delta T) \propto e^{- B(\nu) (Delta T)^2} to fit the three-year WMAP data.

In a recent letter (PRL 95, 190601 (2005)), T.M. Rocha Filho and coworkers address the very interesting issue of the entropic form to be used for Hamiltonians with long-range interactions. In our opinion the letter misses several points quite debated in the recent literature which are of fundamental importance for a complete discussion of this problem. Moreover it contains several statements which are not true or not corroborated by any evidence. In this comment we discuss these arguments and severely question the generality of the conclusions of this letter.

Temperature maps of the Cosmic Microwave Background (CMB) radiation, as those obtained by the Wilkinson Microwave Anisotropy Probe (WMAP), provide one of the most precise data sets to test fundamental hypotheses of modern cosmology. One of these issues is related to the statistical properties of the CMB temperature fluctuations, which would have been produced by Gaussian random density fluctuations when matter and radiation were in thermal equilibrium in the early Universe. We analysed here the WMAP data and found that the distribution of the CMB temperature fluctuations P^{CMB}(Delta T) can be quite well fitted by the anomalous temperature distribution emerging within nonextensive statistical mechanics. This theory is based on the nonextensive entropy S_q = k (1 - \int dx [P_q(x)]^q) /(q-1), with the Boltzmann-Gibbs expression as the limit case q -> 1. For the frequencies investigated (\nu= 40.7, 60.8, and 93.5 GHz), we found that P^{CMB}(Delta T) is well described by P_q(Delta T) \propto 1/[1 + (q-1) B(\nu) (Delta T)^2]^{1/(q-1)}, with q = 1.055 \pm 0.002, which exclude, at the 99% confidence level, exact Gaussian temperature distributions P^{Gauss}(Delta T) \propto e^{- B(\nu) Delta T^2}, corresponding to the q -> 1 limit, to properly represent the CMB temperature fluctuations measured by WMAP.

The analysis of cosmic rays fluxes as a function of energy reveals a {\it knee} slightly below $10^{16}$ eV and an {\it ankle} close to $10^{19}$ eV. Their physical origins remain up to now quite enigmatic; in particular, no elementary process is known which occurs at energies close to $10^{16}$ eV. We propose a phenomenological approach along the lines of nonextensive statistical mechanics, a formalism which contains Boltzmann-Gibbs statistical mechanics as a particular case. The knee then appears as a crossover between two fractal-like thermal regimes, the crossover being caused by process occurring at energies ten million times lower than that of the knee, in the region of the quark hadron transition ($\simeq 10^{9}$ eV). This opens the door to an unexpected standpoint for further clarifying the phenomenon.

We study the dynamics of a Hamiltonian system of N classical spins with infinite-range interaction. We present numerical results which confirm the existence of metaequilibrium Quasi Stationary States (QSS), characterized by non-Gaussian velocity distributions, anomalous diffusion, L\'evy walks and dynamical correlation in phase-space. We show that the Thermodynamic Limit (TL) and the Infinite-Time Limit (ITL) do not commute. Moreover, if the TL is taken before the ITL the system does not relax to the Boltzmann-Gibbs equilibrium, but remains in this new equilibrium state where nonextensive thermodynamics seems to apply.

We study the dynamics of a system of N classical spins with infinite-range interaction. We show that, when the thermodynamic limit is taken before the infinite-time limit, the system does not relax to the Boltzmann-Gibbs equilibrium, but exhibits different equilibrium properties, characterized by stable non-Gaussian velocity distributions and dynamical power-law correlation in phase-space.

Under certain conditions, the rate of increase of the statistical entropy of a simple, fully chaotic, conservative system is known to be given by a single number, characteristic of this system, the Kolmogorov-Sinai entropy rate. This connection is here generalized to a simple dissipative system, the logistic map, and especially to the chaos threshold of the latter, the edge of chaos. It is found that, in the edge-of-chaos case, the usual Boltzmann-Gibbs-Shannon entropy is not appropriate. Instead, the non-extensive entropy $S_q\equiv \frac{1-\sum_{i=1}^W p_i^q}{q-1}$, must be used. The latter contains a parameter q, the entropic index which must be given a special value $q^*\ne 1$ (for q=1 one recovers the usual entropy) characteristic of the edge-of-chaos under consideration. The same q^* enters also in the description of the sensitivity to initial conditions, as well as in that of the multifractal spectrum of the attractor.

We show that the observational data recently provided by Giovanelli et al. (1996 a, b) and discussed by Bahcall and Oh (1996) concerning the velocity distribution of clusters of galaxies can be naturally fitted by a statistical distribution which generalizes the Maxwell-Boltzmann one (herein recovered for the entropic index $q=1$). Indeed a recent generalization of the Boltzmann-Gibbs thermostatistics suggests for this problem that the probability function is given, within a simple phenomenological model, by \begin{eqnarray} P(>v) & = & \frac{ \int_{v}^{v_{max}} dv \left[ 1-(1-q)(v/v_0)^2 \right]^\frac{q}{1-q} } {\int_{0}^{v_{max}} dv \left[ 1-(1-q)(v/v_0)^2 \right]^\frac{q}{1-q}}, \nonumber \\ v_{max} & \equiv & \left\{ \begin{array}{ll} v_0 (1-q)^{-1/2} & \mbox{if $q<1$} \\ \infty & \mbox{if $q \geq 1$.} \end{array}\right. \nonumber \end{eqnarray} A remarkably good fitting with the data is obtained for $q=0.23^{+0.07}_{-0.05}$ and $v_0=490\pm 5$ km s$^{-1}$.