Source author record

Constantino Tsallis

Constantino Tsallis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

61works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

61 published item(s)

preprint2026arXiv

Universal and non-universal facets of quantum critical phenomena unveiled along the Schmidt decomposition theorem

Critical phenomena have been extensively investigated both theoretically and experimentally in many fields, such as condensed matter physics, biology, e.g., brain criticality, and cosmology. In particular, the behaviour of response functions right at critical points (CPs) is highly topical. It turns out that in the frame of Boltzmann-Gibbs-von Neumann-Shannon approach, the extensive character of entropy breaks down at CPs. The latter implies diverging susceptibilities, which is at odds with experimental observations. Here, we investigate the influence of the spin magnitude $S$ on the quantum Grüneisen parameter $Γ^{0\text{K}}_{q}$ right at CPs for the 1D Ising model under a transverse magnetic field. Our findings are fourfold: $\textit{i}$) for higher $S$, $Γ^{0\text{K}}_{q}$ is increased, but remains finite, reflecting the enhancement of the Hilbert space dimensionality; $\textit{ii}$) the Schmidt decomposition theorem recovers the extensivity of the nonadditive $q$-entropy $S_q$ only for a $\textit{special}$ value of the entropic index $q$; $\textit{iii}$) the universality class in the frame of $S_q$ depends only on the symmetry of the system; $\textit{iv}$) we propose an experimental setup to explore finite-size effects in connection with the Hilbert space occupation at CPs. Our findings unveil both universal and non-universal aspects of quantum criticality in terms of $Γ^{0\text{K}}_{q}$ and $S_q$.

preprint2022arXiv

Complex network growth model: Possible isomorphism between nonextensive statistical mechanics and random geometry

In the realm of Boltzmann-Gibbs statistical mechanics there are three well known isomorphic connections with random geometry, namely (i) the Kasteleyn-Fortuin theorem which connects the $λ\to 1$ limit of the $λ$-state Potts ferromagnet with bond percolation, (ii) the isomorphism which connects the $λ\to 0$ limit of the $λ$-state Potts ferromagnet with random resistor networks, and (iii) the de Gennes isomorphism which connects the $n \to 0$ limit of the $n$-vector ferromagnet with self-avoiding random walk in linear polymers. We provide here strong numerical evidence that a similar isomorphism appears to emerge connecting the energy $q$-exponential distribution $\propto e_q^{-β_q \varepsilon}$ (with $q=4/3$ and $β_q ω_0 =10/3$) optimizing, under simple constraints, the nonadditive entropy $S_q$ with a specific geographic growth random model based on preferential attachment through exponentially-distributed weighted links, $ω_0$ being the characteristic weight.

preprint2022arXiv

Statistical mechanical approach of complex networks with weighted links

Systems which consist of many localized constituents interacting with each other can be represented by complex networks. Consistently, network science has become highly popular in vast fields focusing on natural, artificial and social systems. We numerically analyze the growth of $d$-dimensional geographic networks (characterized by the index $α_G\geq0$; $d = 1, 2, 3, 4$) whose links are weighted through a predefined random probability distribution, namely $P(w) \propto e^{-|w - w_c|/τ}$, $w$ being the weight $ (w_c \geq 0; \; τ> 0)$. In this model, each site has an evolving degree $k_i$ and a local energy $\varepsilon_i \equiv \sum_{j=1}^{k_i} w_{ij}/2$ ($i = 1, 2, ..., N$) that depend on the weights of the links connected to it. Each newly arriving site links to one of the pre-existing ones through preferential attachment given by the probability $Π_{ij}\propto \varepsilon_{i}/d^{\,α_A}_{ij} \;\;(α_A \ge 0)$, where $d_{ij}$ is the Euclidean distance between the sites. Short- and long-range interactions respectively correspond to $α_A/d>1$ and $0\leq α_A/d \leq 1$; $α_A/d \to \infty$ corresponds to interactions between close neighbors, and $α_A/d \to 0$ corresponds to infinitely-ranged interactions. The site energy distribution $p(\varepsilon)$ corresponds to the usual degree distribution $p(k)$ as the particular instance $(w_c,τ)=(2,0)$. We numerically verify that the corresponding connectivity distribution $p(\varepsilon)$ converges, when $α_A/d\to\infty$, to the weight distribution $P(w)$ for infinitely narrow distributions (i.e., $τ\to \infty, \,\forall w_c$) as well as for $w_c\to0, \, \forallτ$.

preprint2021arXiv

Along the lines of nonadditive entropies: $q$-prime numbers and $q$-zeta functions

The rich history of prime numbers includes great names such as Euclid, who first analytically studied the prime numbers and proved that there is an infinite number of them, Euler, who introduced the function $ζ(s)\equiv\sum_{n=1}^\infty n^{-s}=\prod_{p\,prime} \frac{1}{1- p^{-s}}$, Gauss, who estimated the rate at which prime numbers increase, and Riemann, who extended $ζ(s)$ to the complex plane $z$ and conjectured that all nontrivial zeros are in the $\mathbb{R}(z)=1/2$ axis. The nonadditive entropy $S_q=k\sum_ip_i\ln_q(1/p_i)\;(q\in\mathbb{R})$ involves the function $\ln_q z\equiv\frac{z^{1-q}-1}{1-q}\;(\ln_1 z=\ln z)$ that is interconnected to a $q$-generalized algebra, using $q$-numbers defined as $\langle x\rangle_q\equiv e^{\ln_q x}$ ($\langle x\rangle_1\equiv x$). The $q$-prime numbers are then defined as the $q$-natural numbers $\langle n\rangle_q\equiv e^{\ln_q n}\;(n=1,2,3,\dots)$, where $n$ is a prime number $p=2,3,5,7,\dots$ We show that, for any value of $q$, infinitely many $q$-prime numbers exist; for $q\le1$ they diverge for increasing prime number, whereas they converge for $q>1$; the standard prime numbers are recovered for $q=1$. For $q\le 1$, we generalize the $ζ(s)$ function as follows: $ζ_q(s)\equiv\langleζ(s)\rangle_q$ ($s\in\mathbb{R}$). We show that this function appears to diverge at $s=1+0$, $\forall q$. Also, we alternatively define, for $q\le 1$, $ζ_q^Σ(s)\equiv\sum_{n=1}^\infty\frac{1}{\langle n\rangle_q^s}=1+\frac{1}{\langle 2\rangle_q^s}+\dots $ and $ζ_q^Π(s)\equiv\prod_{p\,prime}\frac{1}{1-\langle p\rangle_q^{-s}}=\frac{1}{1-\langle 2\rangle_q^{-s}}\frac{1}{1-\langle 3\rangle_q^{-s}}\frac{1}{1-\langle 5\rangle_q^{-s}}\cdots$, which, for $q<1$, generically satisfy $ζ_q^Σ(s)<ζ_q^Π(s)$, in variance with the $q=1$ case, where of course $ζ_1^Σ(s)=ζ_1^Π(s)$.

preprint2021arXiv

Entropic extensivity and large deviations in the presence of strong correlations

The standard Large Deviation Theory (LDT) mirrors the Boltzmann-Gibbs (BG) factor which describes the thermal equilibrium of short-range Hamiltonian systems, the velocity distribution of which is Maxwellian. It is generically applicable to systems satisfying the Central Limit Theorem (CLT), among others. When we focus instead on stationary states of typical complex systems (e.g., classical long-range Hamiltonian systems), both the CLT and LDT need to be generalized. We focus here on a scale-invariant stochastic process involving strongly-correlated exchangeable random variables which, through the Laplace-de Finetti theorem, is known to yield a long-tailed $Q$-Gaussian $N\to\infty$ attractor in the space of distributions ($1 < Q<3)$. We present strong numerical indications that the corresponding LDT probability distribution is given by $P(N,z)=P_0\,e_q^{-r_q(z)N}=P_0[1-(1-q)r_q(z)N]^{1/(1-q)}$ with $q=2-1/Q \in (1,5/3)$. The rate function $r_q(z)$ seemingly equals the relative nonadditive $q_r$-entropy per particle, with $q_r \simeq \frac{7}{10} + \frac{6}{10}\frac{1}{Q-1}$, thus exhibiting a singularity at $Q=1$ and recovering the BG value $q_r=1$ in the $Q \to 3$ limit. Let us emphasize that the extensivity of $r_q(z)N$ appears to be verified, consistently with what is expected, from the Legendre structure of thermodynamics, for a total entropy. The present analysis of a relatively simple model somewhat mirroring spin-1/2 long-range-interacting ferromagnets (e.g., with strongly anisotropic XY coupling) might be helpful for a deeper understanding of nonequilibrium systems with global correlations and other complex systems.

preprint2021arXiv

Generalized Box-Muller method for generating q-Gaussian random deviates

Addendum: The generalized Box-Müller algorithm provides a methodology for generating q-Gaussian random variates. The parameter $-\infty<q\leq3$ is related to the shape of the tail decay; $q<1$ for compact-support including parabola $(q=0)$; $1<q\leq3$ for heavy-tail including Cauchy $(q=2)$. This addendum clarifies the transformation $q'=((3q-1)/(q+1))$ within the algorithm is due to a difference in the dimensions d of the generalized logarithm and the generalized distribution. The transformation is clarified by the decomposition of $q=1+2κ/(1+dκ)$, where the shape parameter $-1<κ\leq\infty$ quantifies the magnitude of the deformation from exponential. A simpler specification for the generalized Box- Müller algorithm is provided using the shape of the tail decay. Original: The q-Gaussian distribution is known to be an attractor of certain correlated systems, and is the distribution which, under appropriate constraints, maximizes the entropy Sq, basis of nonextensive statistical mechanics. This theory is postulated as a natural extension of the standard (Boltzmann-Gibbs) statistical mechanics, and may explain the ubiquitous appearance of heavy-tailed distributions in both natural and man-made systems. The q-Gaussian distribution is also used as a numerical tool, for example as a visiting distribution in Generalized Simulated Annealing. We develop and present a simple, easy to implement numerical method for generating random deviates from a q-Gaussian distribution based upon a generalization of the well known Box-Muller method. Our method is suitable for a larger range of q values, q<3, than has previously appeared in the literature, and can generate deviates from q-Gaussian distributions of arbitrary width and center. MATLAB code showing a straightforward implementation is also included.

preprint2020arXiv

Connecting complex networks to nonadditive entropies

Boltzmann-Gibbs statistical mechanics applies satisfactorily to a plethora of systems. It fails however for complex systems generically involving strong space-time entanglement. Its generalization based on nonadditive $q$-entropies adequately handles a wide class of such systems. We show here that scale-invariant networks belong to this class. We numerically study a $d$-dimensional geographically located network with weighted links and exhibit its 'energy' distribution per site at its quasi-stationary state. Our results strongly suggest a correspondence between the random geometric problem and a class of thermal problems within the generalised thermostatistics. The Boltzmann-Gibbs exponential factor is generically substituted by its $q$-generalisation, and is recovered in the $q=1$ limit when the nonlocal effects fade away. The present connection should cross-fertilise experiments in both research areas.

preprint2020arXiv

Dynamical robustness of discrete conservative systems: Harper and generalized standard maps

In recent years, statistical characterization of the discrete conservative dynamical systems (more precisely, paradigmatic examples of area-preserving maps such as the standard and the web maps) has been analyzed extensively and shown that, for larger parameter values for which the Lyapunov exponents are largely positive over the entire phase space, the probability distribution is a Gaussian, consistent with Boltzmann-Gibbs (BG) statistics. On the other hand, for smaller parameter values for which the Lyapunov exponents are virtually zero over the entire phase space, we verify this distribution appears to approach a $q$-Gaussian (with $q \simeq 1.935$), consistent with $q$-statistics. Interestingly, if the parameter values are in between these two extremes, then the probability distributions happen to exhibit a linear combination of these two behaviours. Here, we numerically show that the Harper map is also in the same universality class of the maps discussed so far. This constitutes one more evidence on the robustness of this behavior whenever the phase space consists of stable orbits. Then, we propose a generalization of the standard map for which the phase space includes many sticky regions, changing the previously observed simple linear combination behavior to a more complex combination.

preprint2019arXiv

A generalised model for asymptotically-scale-free geographical networks

We consider a generalised d-dimensional model for asymptotically-scale-free geographical networks. Central to many networks of this kind, when considering their growth in time, is the attachment rule, i.e. the probability that a new node is attached to one (or more) preexistent nodes. In order to be more realistic, a fitness parameter $η_i \in [0,1]$ for each node $i$ of the network is also taken into account to reflect the ability of the nodes to attract new ones. Our d-dimensional model takes into account the geographical distances between nodes, with different probability distribution for $η$ which sensibly modifies the growth dynamics. The preferential attachment rule is assumed to be $Π_i\propto k_i η_i r_{ij}^{-α_A} $ where $k_i$ is the connectivity of the $i$th pre-existing site and $α_A$ characterizes the importance of the euclidean distance r for the network growth. For special values of the parameters, this model recovers respectively the Bianconi-Barabási and the Barabási-Albert ones. The present generalised model is asymptotically scale-free in all cases, and its degree distribution is very well fitted with q-exponential distributions, which optimise the nonadditive entropy $S_q$, given by $p(k) \propto e_q^{-k/κ} \equiv 1/[1+(q-1)k/κ]^{1/(q-1)}$, with $(q,κ)$ depending uniquely only on the ratio $α_A/d$ and the fitness distribution. Hence this model constitutes a realization of asymptotically-scale-free geographical networks within nonextensive statistical mechanics, where $k$ plays the role of energy and $κ$ plays the role of temperature. General scaling laws are also found for q as a function of the parameters of the model.

preprint2019arXiv

Möbius transforms, cycles and q-triplets in statistical mechanics

In the realm of Boltzmann-Gibbs (BG) statistical mechanics and its q-generalisation for complex systems, we analyse observed sequences of q-triplets, or q-doublets if one of them is the unity, in terms of cycles of successive Möbius transforms of the line preserving unity ( q=1 corresponds to the BG theory). Such transforms have the form q --> (aq + 1-a)/[(1+a)q -a], where a is a real number; the particular cases a=-1 and a=0 yield respectively q --> (2-q) and q --> 1/q, currently known as additive and multiplicative dualities. This approach seemingly enables the organisation of various complex phenomena into different classes, named N-complete or incomplete. The classification that we propose here hopefully constitutes a useful guideline in the search, for non-BG systems whenever well described through q-indices, of new possibly observable physical properties.

preprint2016arXiv

Approach of complexity in nature: Entropic nonuniqueness

Boltzmann introduced in the 1870's a logarithmic measure for the connection between the thermodynamical entropy and the probabilities of the microscopic configurations of the system. His entropic functional for classical systems was extended by Gibbs to the entire phase space of a many-body system, and by von Neumann in order to cover quantum systems as well. Finally, it was used by Shannon within the theory of information. The simplest expression of this functional corresponds to a discrete set of $W$ microscopic possibilities, and is given by $S_{BG}= -k\sum_{i=1}^W p_i \ln p_i$ ($k$ is a positive universal constant; {\it BG} stands for {\it Boltzmann-Gibbs}). This relation enables the construction of BG statistical mechanics. The BG theory has provided uncountable important applications. Its application in physical systems is legitimate whenever the hypothesis of {\it ergodicity} is satisfied. However, {\it what can we do when ergodicity and similar simple hypotheses are violated?}, which indeed happens in very many natural, artificial and social complex systems. It was advanced in 1988 the possibility of generalizing BG statistical mechanics through a family of nonadditive entropies, namely $S_q=k\frac{1-\sum_{i=1}^W p_i^q}{q-1}$, which recovers the additive $S_{BG}$ entropy in the $q \to1$ limit. The index $q$ is to be determined from mechanical first principles. Along three decades, this idea intensively evolved world-wide (see Bibliography in \url{http://tsallis.cat.cbpf.br/biblio.htm}), and led to a plethora of predictions, verifications, and applications in physical systems and elsewhere. As expected whenever a {\it paradigm shift} is explored, some controversy naturally emerged as well in the community. The present status of the general picture is here described, starting from its dynamical and thermodynamical foundations, and ending with its most recent physical applications.

preprint2016arXiv

Dynamics and Statistics of the Fermi--Pasta--Ulam $β$--model with different ranges of particle interactions

In the present work we study the Fermi--Pasta--Ulam (FPU) $β$--model involving long--range interactions (LRI) in both the quadratic and quartic potentials, by introducing two independent exponents $α_1$ and $α_2$ respectively, which make the {forces decay} with distance $r$. Our results demonstrate that weak chaos, in the sense of decreasing Lyapunov exponents, and $q$--Gaussian probability density functions (pdfs) of sums of the momenta, occurs only when long--range interactions are included in the quartic part. More importantly, for $0\leq α_2<1$, we obtain extrapolated values for $q \equiv q_\infty >1$, as $N\rightarrow \infty$, suggesting that these pdfs persist in that limit. On the other hand, when long--range interactions are imposed only on the quadratic part, strong chaos and purely Gaussian pdfs are always obtained for the momenta. We have also focused on similar pdfs for the particle energies and have obtained $q_E$-exponentials (with $q_E>1$) when the quartic-term interactions are long--ranged, otherwise we get the standard Boltzmann-Gibbs weight, with $q=1$. The values of $q_E$ coincide, within small discrepancies, with the values of $q$ obtained by the momentum distributions.

preprint2016arXiv

Inter-occurrence times and universal laws in finance, earthquakes and genomes

A plethora of natural, artificial and social systems exist which do not belong to the Boltzmann-Gibbs (BG) statistical-mechanical world, based on the standard additive entropy $S_{BG}$ and its associated exponential BG factor. Frequent behaviors in such complex systems have been shown to be closely related to $q$-statistics instead, based on the nonadditive entropy $S_q$ (with $S_1=S_{BG}$), and its associated $q$-exponential factor which generalizes the usual BG one. In fact, a wide range of phenomena of quite different nature exist which can be described and, in the simplest cases, understood through analytic (and explicit) functions and probability distributions which exhibit some universal features. Universality classes are concomitantly observed which can be characterized through indices such as $q$. We will exhibit here some such cases, namely concerning the distribution of inter-occurrence (or inter-event) times in the areas of finance, earthquakes and genomes.

preprint2016arXiv

Nonlinear inhomogeneous Fokker-Planck equations: entropy and free-energy time evolution

We extend a recently introduced free-energy formalism for homogeneous Fokker-Planck equations to a wide, and physically appealing, class of inhomogeneous nonlinear Fokker-Planck equations. In our approach, the free-energy functional is expressed in terms of an entropic functional and an auxiliary potential, both derived from the coefficients of the equation. With reference to the introduced entropic functional, we discuss the entropy production in a relaxation process towards equilibrium. The properties of the stationary solutions of the considered Fokker-Planck equations are also discussed.

preprint2016arXiv

On the connection between linear combination of entropies and linear combination of extremizing distributions

We analyze the distribution that extremizes a linear combination of the Boltzmann--Gibbs entropy and the nonadditive $q$-entropy. We show that this distribution can be expressed in terms of a Lambert function. Both the entropic functional and the extremizing distribution can be associated with a nonlinear Fokker--Planck equation obtained from a master equation with nonlinear transition rates. Also, we evaluate the entropy extremized by a linear combination of a Gaussian distribution (which extremizes the Boltzmann--Gibbs entropy) and a $q$-Gaussian distribution (which extremizes the $q$-entropy). We give its explicit expression for $q=0$, and discuss the other cases numerically. The entropy that we obtain can be expressed, for $q=0$, in terms of Lambert functions, and exhibits a discontinuity in the second derivative for all values of $q<1$. The entire discussion is closely related to recent results for type-II superconductors and for the statistics of the standard map.

preprint2016arXiv

Sensitivity to initial conditions of a $d$-dimensional long-range-interacting quartic Fermi-Pasta-Ulam model: Universal scaling

We introduce a generalized $d$-dimensional Fermi-Pasta-Ulam (FPU) model in presence of long-range interactions, and perform a first-principle study of its chaos for $d=1,2,3$ through large-scale numerical simulations. The nonlinear interaction is assumed to decay algebraically as $d_{ij}^{-α}$ ($α\ge 0$), $\{d_{ij}\}$ being the distances between $N$ oscillator sites. Starting from random initial conditions we compute the maximal Lyapunov exponent $λ_{max}$ as a function of $N$. Our $N>>1$ results strongly indicate that $λ_{max}$ remains constant and positive for $α/d>1$ (implying strong chaos, mixing and ergodicity), and that it vanishes like $N^{-κ}$ for $0 \le α/d < 1$ (thus approaching weak chaos and opening the possibility of breakdown of ergodicity). The suitably rescaled exponent $κ$ exhibits universal scaling, namely that $(d+2) κ$ depends only on $α/d$ and, when $α/d$ increases from zero to unity, it monotonically decreases from unity to zero, remaining so for all $α/d >1$. The value $α/d=1$ can therefore be seen as a critical point separating the ergodic regime from the anomalous one, $κ$ playing a role analogous to that of an order parameter. This scaling law is consistent with Boltzmann-Gibbs statistics for $α/d > 1$, and possibly with $q$-statistics for $0 \le α/d < 1$.

preprint2016arXiv

Statistical mechanics for complex systems: On the structure of $q$-triplets

A plethora of natural, artificial and social complex systems exists which violate the basic hypothesis (e.g., ergodicity) of Boltzmann-Gibbs (BG) statistical mechanics. Many of such cases can be satisfactorily handled by introducing nonadditive entropic functionals, such as $S_q\equiv k\frac{1-\sum_{i=1}^W p_i^q}{q-1} \; \Bigl(q \in {\cal R}; \, \sum_{i=1}^W p_i=1 \Bigr)$, with $S_1=S_{BG}\equiv -k\sum_{i=1}^W p_i \ln p_i$. Each class of such systems can be characterized by a set of values $\{q\}$, directly corresponding to its various physical/dynamical/geometrical properties. A most important subset is usually referred to as the $q$-triplet, namely $(q_{sensitivity}, q_{relaxation}, q_{stationary\,state})$, defined in the body of this paper. In the BG limit we have $q_{sensitivity}=q_{relaxation}=q_{stationary\,state}=1$. For a given class of complex systems, the set $\{q\}$ contains only a few independent values of $q$, all the others being functions of those few. An illustration of this structure was given in 2005 [Tsallis, Gell-Mann and Sato, Proc. Natl. Acad. Sc. USA {\bf 102}, 15377; TGS]. This illustration enabled a satisfactory analysis of the Voyager 1 data on the solar wind. But the general form of these structures still is an open question. This is so, for instance, for the challenging $q$-triplet associated with the edge of chaos of the logistic map. We introduce here a transformation which sensibly generalizes the TGS one, and which might constitute an important step towards the general solution.

preprint2016arXiv

The limit distribution in the $q$-CLT for $q \ge 1$ is unique and can not have a compact support

In a paper by Umarov, Tsallis and Steinberg (2008), a generalization of the Fourier transform, called the $q$-Fourier transform, was introduced and applied for the proof of a $q$-generalized central limit theorem ($q$-CLT). Subsequently, Hilhorst illustrated (2009 and 2010) that the $q$-Fourier transform for $q>1$ is not invertible in the space of density functions. Indeed, using an invariance principle, he constructed a family of densities with the same $q$-Fourier transform and noted that "as a consequence, the $q$-central limit theorem falls short of achieving its stated goal". The distributions constructed there have compact support. We prove now that the limit distribution in the $q$-CLT is unique and can not have a compact support. This result excludes all the possible counterexamples which can be constructed using the invariance principle and fills the gap mentioned by Hilhorst.

preprint2015arXiv

A closer look at coupled logistic maps at the edge of chaos

We focus on a linear chain of $N$ first-neighbor-coupled logistic maps at their edge of chaos in the presence of a common noise. This model, characterised by the coupling strength $ε$ and the noise width $σ_{max}$, was recently introduced by Pluchino et al [Phys. Rev. E {\bf 87}, 022910 (2013)]. They detected, for the time averaged returns with characteristic return time $τ$, possible connections with $q$-Gaussians, the distributions which optimise, under appropriate constraints, the nonadditive entropy $S_q$, basis of nonextensive statistics mechanics. We have here a closer look on this model, and numerically obtain probability distributions which exhibit a slight asymmetry for some parameter values, in variance with simple $q$-Gaussians. Nevertheless, along many decades, the fitting with $q$-Gaussians turns out to be numerically very satisfactory for wide regions of the parameter values, and we illustrate how the index $q$ evolves with $(N, τ, ε, σ_{max})$. It is nevertheless instructive on how careful one must be in such numerical analysis. The overall work shows that physical and/or biological systems that are correctly mimicked by the Pluchino et al model are thermostatistically related to nonextensive statistical mechanics when time-averaged relevant quantities are studied.

preprint2015arXiv

A new entropy based on a group-theoretical structure

A multi-parametric version of the nonadditive entropy $S_{q}$ is introduced. This new entropic form, denoted by $S_{a,b,r}$, possesses many interesting statistical properties, and it reduces to the entropy $S_q$ for $b=0$, $a=r:=1-q$ (hence Boltzmann-Gibbs entropy $S_{BG}$ for $b=0$, $a=r \to 0$). The construction of the entropy $S_{a,b,r}$ is based on a general group-theoretical approach recently proposed by one of us \cite{Tempesta2}. Indeed, essentially all the properties of this new entropy are obtained as a consequence of the existence of a rational group law, which expresses the structure of $S_{a,b,r}$ with respect to the composition of statistically independent subsystems. Depending on the choice of the parameters, the entropy $S_{a,b,r}$ can be used to cover a wide range of physical situations, in which the measure of the accessible phase space increases say exponentially with the number of particles $N$ of the system, or even stabilizes, by increasing $N$, to a limiting value. This paves the way to the use of this entropy in contexts where a system "freezes" some or many of its degrees of freedom by increasing the number of its constituting particles or subsystems.

preprint2015arXiv

Conceptual Inadequacy of the Shore and Johnson Axioms for Wide Classes of Complex Systems

It is by now well known that the Boltzmann-Gibbs-von Neumann-Shannon logarithmic entropic functional ($S_{BG}$) is inadequate for wide classes of strongly correlated systems: see for instance the 2001 Brukner and Zeilinger's {\it Conceptual inadequacy of the Shannon information in quantum measurements}, among many other systems exhibiting various forms of complexity. On the other hand, the Shannon and Khinchin axioms uniquely mandate the BG form $S_{BG}=-k\sum_i p_i \ln p_i$; the Shore and Johnson axioms follow the same path. Many natural, artificial and social systems have been satisfactorily approached with nonadditive entropies such as the $S_q=k \frac{1-\sum_i p_i^q}{q-1}$ one ($q \in {\cal R}; \,S_1=S_{BG}$), basis of nonextensive statistical mechanics. Consistently, the Shannon 1948 and Khinchine 1953 uniqueness theorems have already been generalized in the literature, by Santos 1997 and Abe 2000 respectively, in order to uniquely mandate $S_q$. We argue here that the same remains to be done with the Shore and Johnson 1980 axioms. We arrive to this conclusion by analyzing specific classes of strongly correlated complex systems that await such generalization.

preprint2015arXiv

Convergence of the probability of large deviations in a model of correlated random variables having compact-support $Q$-Gaussians as limiting distributions

We consider correlated random variables $X_1,\dots,X_n$ taking values in $\{0,1\}$ such that, for any permutation $π$ of $\{1,\dots,n\}$, the random vectors $(X_1,\dots,X_n)$ and $(X_{π(1)},\dots,X_{π(n)})$ have the same distribution. This distribution, which was introduced by Rodríguez et al (2008) and then generalized by Hanel et al (2009), is scale-invariant and depends on a real parameter $ν>0$ ($ν\to\infty$ implies independence). Putting $S_n=X_1+\cdots+X_n$, the distribution of $S_n-n/2$ approaches a $Q$-Gaussian distribution with compact support ($Q=1-1/(ν-1)<1$) as $n$ increases, after appropriate scaling. In the present article, we show that the distribution of $S_n/n$ converges, as $n\to\infty$, to a beta distribution with both parameters equal to $ν$. In particular, the law of large numbers does not hold since, if $0\le x<1/2$, then $\mathbb{P}(S_n/n\le x)$, which is the probability of the event $\{S_n/n\le x\}$ (large deviation), does not converges to zero as $n\to\infty$. For $x=0$ and every real $ν>0$, we show that $\mathbb{P}(S_n=0)$ decays to zero like a power law of the form $1/n^ν$ with a subdominant term of the form $1/n^{ν+1}$. If $0<x\le 1$ and $ν>0$ is an integer, we show that we can analytically find upper and lower bounds for the difference between $\mathbb{P}(S_n/n\le x)$ and its ($n\to\infty$) limit. We also show that these bounds vanish like a power law of the form $1/n$ with a subdominant term of the form $1/n^2$.

preprint2015arXiv

From QCD-based hard-scattering to nonextensive statistical mechanical descriptions of transverse momentum spectra in high-energy $pp$ and $p\bar p$ collisions

Transverse spectra of both jets and hadrons obtained in high-energy $pp$ and $p\bar p $ collisions at central rapidity exhibit power-law behavior of $1/p_T^n$ at high $p_T$. The power index $n$ is 4-5 for jet production and is 6-10 for hadron production. Furthermore, the hadron spectra spanning over 14 orders of magnitude down to the lowest $p_T$ region in $pp$ collisions at LHC can be adequately described by a single nonextensive statistical mechanical distribution that is widely used in other branches of science. This suggests indirectly the possible dominance of the hard-scattering process over essentially the whole $p_T$ region at central rapidity in high-energy $pp$ and $p \bar p$ collisions. We show here direct evidences of such a dominance of the hard-scattering process by investigating the power indices of UA1 and ATLAS jet spectra over an extended $p_T$ region and the two-particle correlation data of the STAR and PHENIX Collaborations in high-energy $pp$ and $p \bar p$ collisions at central rapidity. We then study how the showering of the hard-scattering product partons alters the power index of the hadron spectra and leads to a hadron distribution that may be cast into a single-particle nonextensive statistical mechanical distribution. Because of such a connection, the nonextensive statistical mechanical distribution may be considered as a lowest-order approximation of the hard-scattering of partons followed by the subsequent process of parton showering that turns the jets into hadrons, in high energy $pp$ and $p\bar p$ collisions.

preprint2015arXiv

On the robustness of the $q$-Gaussian family

We introduce three deformations, called $α$-, $β$- and $γ$-deformation respectively, of a $N$-body probabilistic model, first proposed by Rodríguez et al. (2008), having $q$-Gaussians as $N\to\infty$ limiting probability distributions. The proposed $α$- and $β$-deformations are asymptotically scale-invariant, whereas the $γ$-deformation is not. We prove that, for both $α$- and $β$-deformations, the resulting deformed triangles still have $q$-Gaussians as limiting distributions, with a value of $q$ independent (dependent) on the deformation parameter in the $α$-case ($β$-case). In contrast, the $γ$-case, where we have used the celebrated $Q$-numbers and the Gauss binomial coefficients, yields other limiting probability distribution functions, outside the $q$-Gaussian family. These results suggest that scale-invariance might play an important role regarding the robustness of the $q$-Gaussian family.

preprint2015arXiv

Role of dimensionality in complex networks: Connection with nonextensive statistics

Deep connections are known to exist between scale-free networks and non-Gibbsian statistics. For example, typical degree distributions at the thermodynamical limit are of the form $P(k) \propto e_q^{-k/κ}$, where the $q$-exponential form $e_q^z \equiv [1+(1-q)z]^{\frac{1}{1-q}}$ optimizes the nonadditive entropy $S_q$ (which, for $q\to 1$, recovers the Boltzmann-Gibbs entropy). We introduce and study here $d$-dimensional geographically-located networks which grow with preferential attachment involving Euclidean distances through $r_{ij}^{-α_A} \; (α_A \ge 0)$. Revealing the connection with $q$-statistics, we numerically verify (for $d$ =1, 2, 3 and 4) that the $q$-exponential degree distributions exhibit, for both $q$ and $κ$, universal dependences on the ratio $α_A/d$. Moreover, the $q=1$ limit is rapidly achieved by increasing $α_A/d$ to infinity.

preprint2014arXiv

An introduction to nonadditive entropies and a thermostatistical approach of inanimate and living matter

The possible distinction between inanimate and living matter has been of interest to humanity since thousands of years. Clearly, such a rich question can not be answered in a single manner, and a plethora of approaches naturally do exist. However, during the last two decades, a new standpoint, of thermostatistical nature, has emerged. It is related to the proposal of nonadditive entropies in 1988, in order to generalise the celebrated Boltzmann-Gibbs additive functional, basis of standard statistical mechanics. Such entropies have found deep fundamental interest and uncountable applications in natural, artificial and social systems. In some sense, this perspective represents an epistemological paradigm shift. These entropies crucially concern complex systems, in particular those whose microscopic dynamics violate ergodicity. Among those, living matter and other living-like systems play a central role. We briefly review here this approach, and present some of its predictions, verifications and applications.

preprint2014arXiv

Boltzmann-Gibbs entropy is sufficient but not necessary for the likelihood factorization required by Einstein

In 1910 Einstein published a crucial aspect of his understanding of Boltzmann entropy. He essentially argued that the likelihood function of any system composed by two probabilistically independent subsystems {\it ought} to be factorizable into the likelihood functions of each of the subsystems. Consistently he was satisfied by the fact that Boltzmann (additive) entropy fulfills this epistemologically fundamental requirement. We show here that entropies (e.g., the $q$-entropy on which nonextensive statistical mechanics is based) which generalize the BG one through violation of its well known additivity can {\it also} fulfill the same requirement. This fact sheds light on the very foundations of the connection between the micro- and macro-scopic worlds.

preprint2014arXiv

Emergence of $q$-statistical functions in a generalized binomial distribution with strong correlations

We study a symmetric generalization $\mathfrak{p}^{(N)}_k(η, α)$ of the binomial distribution recently introduced by Bergeron et al, where $η\in [0,1]$ denotes the win probability, and $α$ is a positive parameter. This generalization is based on $q$-exponential generating functions ($e_{q^{gen}}^z \equiv [1+(1-q^{gen})z]^{1/(1-q^{gen})};\,e_{1}^z=e^z)$ where $q^{gen}=1+1/α$. The numerical calculation of the probability distribution function of the number of wins $k$, related to the number of realizations $N$, strongly approaches a discrete $q^{disc}$-Gaussian distribution, for win-loss equiprobability (i.e., $η=1/2$) and all values of $α$. Asymptotic $N\to \infty$ distribution is in fact a $q^{att}$-Gaussian $e_{q^{att}}^{-βz^2}$, where $q^{att}=1-2/(α-2)$ and $β=(2α-4)$. The behavior of the scaled quantity $k/N^γ$ is discussed as well. For $γ<1$, a large-deviation-like property showing a $q^{ldl}$-exponential decay is found, where $q^{ldl}=1+1/(ηα)$. For $η=1/2$, $q^{ldl}$ and $q^{att}$ are related through $1/(q^{ldl}-1)+1/(q^{att}-1)=1$, $\forall α$. For $γ=1$, the law of large numbers is violated, and we consistently study the large-deviations with respect to the probability of the $N\to\infty$ limit distribution, yielding a power law, although not exactly a $q^{LD}$-exponential decay. All $q$-statistical parameters which emerge are univocally defined by $(η, α)$. Finally we discuss the analytical connection with the Pólya urn problem.

preprint2014arXiv

Fermi-Pasta-Ulam model with long-range interactions: Dynamics and thermostatistics

We introduce and numerically study a long-range-interaction generalization of the one-dimensional Fermi-Pasta-Ulam (FPU) $β-$ model. The standard quartic interaction is generalized through a coupling constant that decays as $1/r^α$ ($α\ge 0$)(with strength characterized by $b>0$). In the $α\to\infty$ limit we recover the original FPU model. Through classical molecular dynamics computations we show that (i) For $α\geq 1$ the maximal Lyapunov exponent remains finite and positive for increasing number of oscillators $N$ (thus yielding ergodicity), whereas, for $0 \le α<1$, it asymptotically decreases as $N^{- κ(α)}$ (consistent with violation of ergodicity); (ii) The distribution of time-averaged velocities is Maxwellian for $α$ large enough, whereas it is well approached by a $q$-Gaussian, with the index $q(α)$ monotonically decreasing from about 1.5 to 1 (Gaussian) when $α$ increases from zero to close to one. For $α$ small enough, the whole picture is consistent with a crossover at time $t_c$ from $q$-statistics to Boltzmann-Gibbs (BG) thermostatistics. More precisely, we construct a "phase diagram" for the system in which this crossover occurs through a frontier of the form $1/N \propto b^δ/t_c^γ$ with $γ>0$ and $δ>0$, in such a way that the $q=1$ ($q>1$) behavior dominates in the $\lim_{N \to\infty} \lim_{t \to\infty}$ ordering ($\lim_{t \to\infty} \lim_{N \to\infty}$ ordering).

preprint2014arXiv

Nonlinear inhomogeneous Fokker-Planck equation within a generalized Stratonovich prescription

We deduce a nonlinear and inhomogeneous Fokker-Planck equation within a generalized Stratonovich, or stochastic $α$-, prescription ($α=0$, $1/2$ and $1$ respectively correspond to the Itô, Stratonovich and anti-Itô prescriptions). We obtain its stationary state $p_{st}(x)$ for a class of constitutive relations between drift and diffusion and show that it has a $q$-exponential form, $p_{st}(x) = N_q[1 - (1-q)βV(x)]^{1/(1-q)}$, with an index $q$ which does not depend on $α$ in the presence of any nonvanishing nonlinearity. This is in contrast with the linear case, for which the index $q$ is $α$-dependent.

preprint2014arXiv

Possible Implication of a Single Nonextensive $p_T$ Distribution for Hadron Production in High-Energy $pp$ Collisions

Multiparticle production processes in $pp$ collisions at the central rapidity region are usually considered to be divided into independent "soft" and "hard" components. The first is described by exponential (thermal-like) transverse momentum spectra in the low-$p_T$ region with a scale parameter $T$ associated with the temperature of the hadronizing system. The second is governed by a power-like distributions of transverse momenta with power index $n$ at high-$p_T$ associated with the hard scattering between partons. We show that the hard-scattering integral can be approximated as a nonextensive distribution of a quasi-power-law containing a scale parameter $T$ and a power index $n=1/(q -1)$, where $q$ is the nonextensivity parameter. We demonstrate that the whole region of transverse momenta presently measurable at LHC experiments at central rapidity (in which the observed cross sections varies by $14$ orders of magnitude down to the low $p_T$ region) can be adequately described by a single nonextensive distribution. These results suggest the dominance of the hard-scattering hadron-production process and the approximate validity of a "no-hair" statistical-mechanical description of the $p_T$ spectra for the whole $p_T$ region at central rapidity for $pp$ collisions at high-energies.

preprint2014arXiv

The transverse-momenta distributions in high-energy $pp$ collisions -- A statistical-mechanical approach

We analyze LHC available data measuring the distribution probability of transverse momenta~$p_T$ in proton-proton collisions at $\sqrt{s}=0.9\,\textrm{TeV}$ (CMS, ALICE, ATLAS) and $\sqrt{s}=7\,\textrm{TeV}$ (CMS, ATLAS). A remarkably good fitting can be obtained, along fourteen decades in magnitude, by phenomenologically using $q$-statistics for a {\it single} particle of a two-dimensional relativistic ideal gas. The parameters that have been obtained by assuming $\textrm{d}N/p_T\textrm{d}p_T\textrm{d}y \propto e_q^{-E_T/T}$ at mid-rapidity are, in all cases, $q \simeq 1.1$ and $T\simeq 0.13\,\textrm{GeV}$ (which satisfactorily compares with the pion mass). This fact suggests the approximate validity of a "no-hair" statistical-mechanical description of the hard-scattering hadron-production process in which the detailed mechanisms of parton scattering, parton cascades, parton fragmentation, running coupling and other information can be subsumed under the stochastic dynamics in the lowest-order description. In addition to that basic structure, a finer analysis of the data suggests a small oscillatory structure on top of the leading $q$-exponential. The physical origin of such intriguing oscillatory behavior remains elusive, though it could be related to some sort of fractality or scale-invariance within the system.

preprint2014arXiv

Thermodynamics is more powerful than the role to it reserved by Boltzmann-Gibbs statistical mechanics

We brief{}ly review the connection between statistical mechanics and thermodynamics. We show that, in order to satisfy thermodynamics and its Legendre transformation mathematical frame, the celebrated Boltzmann-Gibbs~(BG) statistical mechanics is suff{}icient but not necessary. Indeed, the $N\to\infty$ limit of statistical mechanics is expected to be consistent with thermodynamics. For systems whose elements are generically independent or quasi-independent in the sense of the theory of probabilities, it is well known that the BG theory (based on the additive BG entropy) does satisfy this expectation. However, in complete analogy, other thermostatistical theories (\emph{e.g.}, $q$-statistics), based on nonadditive entropic functionals, also satisfy the very same expectation. We illustrate this standpoint with systems whose elements are strongly correlated in a specific manner, such that they escape the BG realm.

preprint2013arXiv

Black hole thermodynamical entropy

As early as 1902, Gibbs pointed out that systems whose partition function diverges, e.g. gravitation, lie outside the validity of the Boltzmann-Gibbs (BG) theory. Consistently, since the pioneering Bekenstein-Hawking results, physically meaningful evidence (e.g., the holographic principle) has accumulated that the BG entropy $S_{BG}$ of a $(3+1)$ black hole is proportional to its area $L^2$ ($L$ being a characteristic linear length), and not to its volume $L^3$. Similarly it exists the \emph{area law}, so named because, for a wide class of strongly quantum-entangled $d$-dimensional systems, $S_{BG}$ is proportional to $\ln L$ if $d=1$, and to $L^{d-1}$ if $d>1$, instead of being proportional to $L^d$ ($d \ge 1$). These results violate the extensivity of the thermodynamical entropy of a $d$-dimensional system. This thermodynamical inconsistency disappears if we realize that the thermodynamical entropy of such nonstandard systems is \emph{not} to be identified with the BG {\it additive} entropy but with appropriately generalized {\it nonadditive} entropies. Indeed, the celebrated usefulness of the BG entropy is founded on hypothesis such as relatively weak probabilistic correlations (and their connections to ergodicity, which by no means can be assumed as a general rule of nature). Here we introduce a generalized entropy which, for the Schwarzschild black hole and the area law, can solve the thermodynamic puzzle.

preprint2013arXiv

Influence of the interaction range on the thermostatistics of a classical many-body system

We numerically study a one-dimensional system of $N$ classical localized planar rotators coupled through interactions which decay with distance as $1/r^α$ ($α\ge 0$). The approach is a first principle one (\textit{i.e.}, based on Newton's law), and yields the probability distribution of momenta. For $α$ large enough and $N\gg1$ we observe, for longstanding states, the Maxwellian distribution, landmark of Boltzmann-Gibbs thermostatistics. But, for $α$ small or comparable to unity, we observe instead robust fat-tailed distributions that are quite well fitted with $q$-Gaussians. These distributions extremize, under appropriate simple constraints, the nonadditive entropy $S_q$ upon which nonextensive statistical mechanics is based. The whole scenario appears to be consistent with nonergodicity and with the thesis of the $q$-generalized Central Limit Theorem.

preprint2013arXiv

Noise, Synchrony and Correlations at the Edge of Chaos

We study the effect of a weak random additive noise in a linear chain of N locally-coupled logistic maps at the edge of chaos. Maps tend to synchronize for a strong enough coupling, but if a weak noise is added, very intermittent fluctuations in the returns time series are observed. This intermittency tends to disappear when noise is increased. Considering the pdfs of the returns, we observe the emergence of fat tails which can be satisfactorily reproduced by $q$-Gaussians curves typical of nonextensive statistical mechanics. Interoccurrence times of these extreme events are also studied in detail. Similarities with recent analysis of financial data are also discussed.

preprint2012arXiv

Nonlinear Schroedinger Equation in the Presence of Uniform Acceleration

We consider a recently proposed nonlinear Schroedinger equation exhibiting soliton-like solutions of the power-law form $e_q^{i(kx-wt)}$, involving the $q$-exponential function which naturally emerges within nonextensive thermostatistics [$e_q^z \equiv [1+(1-q)z]^{1/(1-q)}$, with $e_1^z=e^z$]. Since these basic solutions behave like free particles, obeying $p=\hbar k$, $E=\hbar ω$ and $E=p^2/2m$ ($1 \le q<2$), it is relevant to investigate how they change under the effect of uniform acceleration, thus providing the first steps towards the application of the aforementioned nonlinear equation to the study of physical scenarios beyond free particle dynamics. We investigate first the behaviour of the power-law solutions under Galilean transformation and discuss the ensuing Doppler-like effects. We consider then constant acceleration, obtaining new solutions that can be equivalently regarded as describing a free particle viewed from an uniformly accelerated reference frame (with acceleration $a$) or a particle moving under a constant force $-ma$. The latter interpretation naturally leads to the evolution equation $i\hbar \frac{\partial}{\partial t}(\fracΦ{Φ_0}) = - \frac{1}{2-q}\frac{\hbar^2}{2m} \frac{\partial^2}{\partial x^2} [(\fracΦ{Φ_0})^{2-q}] + V(x)(\fracΦ{Φ_0})^{q}$ with $V(x)=max$. Remarkably enough, the potential $V$ couples to $Φ^q$, instead of coupling to $Φ$, as happens in the familiar linear case ($q=1$).

preprint2012arXiv

Reply to Comment on "Towards a large deviation theory for strongly correlated systems"

The paper that is commented by Touchette contains a computational study which opens the door to a desirable generalization of the standard large deviation theory (applicable to a set of $N$ nearly independent random variables) to systems belonging to a special, though ubiquitous, class of strong correlations. It focuses on three inter-related aspects, namely (i) we exhibit strong numerical indications which suggest that the standard exponential probability law is asymptotically replaced by a power-law as its dominant term for large $N$; (ii) the subdominant term appears to be consistent with the $q$-exponential behavior typical of systems following $q$-statistics, thus reinforcing the thermodynamically extensive entropic nature of the exponent of the $q$-exponential, basically $N$ times the $q$-generalized rate function; (iii) the class of strong correlations that we have focused on corresponds to attractors in the sense of the Central Limit Theorem which are $Q$-Gaussian distributions (in principle $1 < Q < 3$), which relevantly differ from (symmetric) Lévy distributions, with the unique exception of Cauchy-Lorentz distributions (which correspond to $Q = 2$), where they coincide, as well known. In his Comment, Touchette has agreeably discussed point (i), but, unfortunately, points (ii) and (iii) have, as we detail here, visibly escaped to his analysis. Consequently, his conclusion claiming the absence of special connection with $q$-exponentials is unjustified.

preprint2011arXiv

A scale-invariant probabilistic model based on Leibniz-like pyramids

We introduce a family of probabilistic {\it scale-invariant} Leibniz-like pyramids and $(d+1)$-dimensional hyperpyramids ($d=1,2,3,...$), characterized by a parameter $ν>0$, whose value determines the degree of correlation between $N$ $(d+1)$-valued random variables. There are $(d+1)^N$ different events, and the limit $ν\to\infty$ corresponds to independent random variables, in which case each event has a probability $1/(d+1)^N$ to occur. The sums of these $N$ $\,(d+1)$-valued random variables correspond to a $d-$dimensional probabilistic model, and generalizes a recently proposed one-dimensional ($d=1$) model having $q-$Gaussians (with $q=(ν-2)/(ν-1)$ for $ν\in [1,\infty)$) as $N\to\infty$ limit probability distributions for the sum of the $N$ binary variables [A. Rodr\'ıguez {\em et al}, J. Stat. Mech. (2008) P09006; R. Hanel {\em et al}, Eur. Phys. J. B {\bf 72}, 263 (2009)]. In the $ν\to\infty$ limit the $d-$dimensional multinomial distribution is recovered for the sums, which approach a $d-$dimensional Gaussian distribution for $N\to\infty$. For any $ν$, the conditional distributions of the $d-$dimensional model are shown to yield the corresponding joint distribution of the $(d-1)$-dimensional model with the same $ν$. For the $d=2$ case, we study the joint probability distribution, and identify two classes of marginal distributions, one of them being asymmetric and scale-invariant, while the other one is symmetric and only asymptotically scale-invariant. The present probabilistic model is proposed as a testing ground for a deeper understanding of the necessary and sufficient conditions for having $q$-Gaussian attractors in the $N\to\infty$ limit, the ultimate goal being a neat mathematical view of the causes clarifying the ubiquitous emergence of $q$-statistics verified in many natural, artificial and social systems.

preprint2011arXiv

Long memory constitutes a unified mesoscopic mechanism consistent with nonextensive statistical mechanics

We unify two paradigmatic mesoscopic mechanisms for the emergence of nonextensive statistics, namely the multiplicative noise mechanism leading to a {\it linear} Fokker-Planck (FP) equation with {\it inhomogenous} diffusion coefficient, and the non-Markovian process leading to the {\it nonlinear} FP equation with {\it homogeneous} diffusion coefficient. More precisely, we consider the equation $\frac{\partial p(x,t)}{\partial t}=-\frac{\partial}{\partial x}[F(x) p(x,t)] + 1/2D \frac{\partial^2}{\partial x^2} [ϕ(x,p)p(x,t)]$, where $D \in {\cal R}$ and $F(x)=-\partial V(x) /\partial x$, $V(x)$ being the potential under which diffusion occurs. Our aim is to find whether $ϕ(x,p)$ exists such that the inhomogeneous linear and the homogeneous nonlinear FP equations become unified in such a way that the (ubiquitously observed) $q$-exponentials remain as stationary solutions. It turns out that such solutions indeed exist for a wide class of systems, namely when $ϕ(x,p)=[A+BV(x)]^θ[p(x,t)]^η$, where $A$, $B$, $θ$ and $η$ are (real) constants. Our main result can be sumarized as follows: For $θ\neq 1$ and arbitrary confining potential $V(x)$, $p(x,\infty) \propto \lbrace 1-β(1-q)V(x)\rbrace ^{1/(1-q)} \equiv e_q^{-βV(x)}$, where $q= 1+ η/(θ-1)$. The present approach unifies into a single mechanism, essentially {\it long memory}, results currently discussed and applied in the literature.

preprint2011arXiv

Nonlinear Relativistic and Quantum Equations with a Common Type of Solution

Generalizations of the three main equations of quantum physics, namely, the Schrödinger, Klein-Gordon, and Dirac equations, are proposed. Nonlinear terms, characterized by exponents depending on an index $q$, are considered in such a way that the standard linear equations are recovered in the limit $q \rightarrow 1$. Interestingly, these equations present a common, soliton-like, travelling solution, which is written in terms of the $q$-exponential function that naturally emerges within nonextensive statistical mechanics. In all cases, the well-known Einstein energy-momentum relation is preserved for arbitrary values of $q$.

preprint2011arXiv

Restricted random walk model as a new testing ground for the applicability of q-statistics

We present exact results obtained from Master Equations for the probability function P(y,T) of sums $y=\sum_{t=1}^T x_t$ of the positions x_t of a discrete random walker restricted to the set of integers between -L and L. We study the asymptotic properties for large values of L and T. For a set of position dependent transition probabilities the functional form of P(y,T) is with very high precision represented by q-Gaussians when T assumes a certain value $T^*\propto L^2$. The domain of y values for which the q-Gaussian apply diverges with L. The fit to a q-Gaussian remains of very high quality even when the exponent $a$ of the transition probability g(x)=|x/L|^a+p with 0<p<<1 is different from 1, all though weak, but essential, deviation from the q-Gaussian does occur for $a\neq1$. To assess the role of correlations we compare the T dependence of P(y,T) for the restricted random walker case with the equivalent dependence for a sum y of uncorrelated variables x each distributed according to 1/g(x).

preprint2011arXiv

Sensitivity, entropy, and escape rates at the onset of chaos

We analytically link three properties of nonlinear dynamical systems, namely sensitivity to initial conditions, entropy production, and escape rate, in $z$-logistic maps for both positive and zero Lyapunov exponents. We unify these relations at chaos, where the Lyapunov exponent is positive, and at its onset, where it vanishes. Our result unifies, in particular, two already known cases, namely (i) the standard entropy rate in the presence of escape, valid for exponential functionality rates with strong chaos, and (ii) the Pesin-like identity with no escape, valid for the power-law behavior present at points such as the Feigenbaum one.

preprint2011arXiv

The nonadditive entropy $S_q$: A door open to the nonuniversality of the mathematical expression of the Clausius thermodynamic entropy in terms of the probabilities of the microscopic configurations

Clausius introduced, in the 1860s, a thermodynamical quantity which he named {\it entropy} $S$. This thermodynamically crucial quantity was proposed to be {\it extensive}, i.e., in contemporary terms, $S(N) \propto N$ in the thermodynamic limit $N \to\infty$. A decade later, Boltzmann proposed a functional form for this quantity which connects $S$ with the occurrence probabilities of the microscopic configurations (referred to as {\it complexions} at that time) of the system. This functional is, if written in modern words referring to a system with $W$ possible discrete states, $S_{BG}=-k_B \sum_{i=1}^W p_i \ln p_i$. The BG entropy is {\it additive}, meaning that, if A and B are two probabilistically independent systems, then $S_{BG}(A+B)=S_{BG}(A)+S_{BG}(B)$. The words, {\it extensive} and {\it additive}, were practically treated, for over more than one century, as almost synonyms, and $S_{BG}$ was considered to be the unique form that $S$ could take. In other words, the functional $S_{BG}$ was considered to be universal. It has become increasingly clear today that it is {\it not} so, and that those two words are {\it not} synonyms, but happen to coincide whenever we are dealing with paradigmatic Hamiltonians involving {\it short-range} interactions between their elements, presenting no strong frustration and other "pathologies". These facts constitute the basis of a generalization of the BG entropy and statistical mechanics, introduced in 1988, and frequently referred to as nonadditive entropy $S_q$ and nonextensive statistical mechanics, respectively. We briefly review herein these points, and exhibit recent as well as typical applications of these concepts in natural, artificial, and social systems, as shown through theoretical, experimental, observational and computational predictions and verifications.

preprint2011arXiv

Towards a large deviation theory for statistical-mechanical complex systems

The theory of large deviations constitutes a mathematical cornerstone in the foundations of Boltzmann-Gibbs statistical mechanics, based on the additive entropy $S_{BG}=- k_B\sum_{i=1}^W p_i \ln p_i$. Its optimization under appropriate constraints yields the celebrated BG weight $e^{-βE_i}$. An elementary large-deviation connection is provided by $N$ independent binary variables, which, in the $N\to\infty$ limit yields a Gaussian distribution. The probability of having $n \ne N/2$ out of $N$ throws is governed by the exponential decay $e^{-N r}$, where the rate function $r$ is directly related to the relative BG entropy. To deal with a wide class of complex systems, nonextensive statistical mechanics has been proposed, based on the nonadditive entropy $S_q=k_B\frac{1- \sum_{i=1}^W p_i^q}{q-1}$ ($q \in {\cal R}; \,S_1=S_{BG}$). Its optimization yields the generalized weight $e_q^{-β_q E_i}$ ($e_q^z \equiv [1+(1-q)z]^{1/(1-q)};\,e_1^z=e^z)$. We numerically study large deviations for a strongly correlated model which depends on the indices $Q \in [1,2)$ and $γ\in (0,1)$. This model provides, in the $N\to\infty$ limit ($\forall γ$), $Q$-Gaussian distributions, ubiquitously observed in nature ($Q\to 1$ recovers the independent binary model). We show that its corresponding large deviations are governed by $e_q^{-N r_q}$ ($\propto 1/N^{1/(q-1)}$ if $q>1$) where $q= \frac{Q-1}{γ(3-Q)}+1 \ge 1$. This $q$-generalized illustration opens wide the door towards a desirable large-deviation foundation of nonextensive statistical mechanics.

preprint2010arXiv

A generalization of the cumulant expansion. Application to a scale-invariant probabilistic model

As well known, cumulant expansion is an alternative way to moment expansion to fully characterize probability distributions provided all the moments exist. If this is not the case, the so called escort mean values (or q-moments) have been proposed to characterize probability densities with divergent moments [C. Tsallis et al, J. Math. Phys 50, 043303 (2009)]. We introduce here a new mathematical object, namely the q-cumulants, which, in analogy to the cumulants, provide an alternative characterization to that of the q-moments for the probability densities. We illustrate this new scheme on a recently proposed family of scale-invariant discrete probabilistic models [A. Rodriguez et al, J. Stat. Mech. (2008) P09006; R. Hanel et al, Eur. Phys. J. B 72, 263268 (2009)] having q-Gaussians as limiting probability distributions.

preprint2010arXiv

Generalization of symmetric $α$-stable Lévy distributions for $q>1$

The $α$-stable distributions introduced by Lévy play an important role in probabilistic theoretical studies and their various applications, e.g., in statistical physics, life sciences, and economics. In the present paper we study sequences of long-range dependent random variables whose distributions have asymptotic power law decay, and which are called $(q,α)$-stable distributions. These sequences are generalizations of i.i.d. $α$-stable distributions, and have not been previously studied. Long-range dependent $(q,α)$-stable distributions might arise in the description of anomalous processes in nonextensive statistical mechanics, cell biology, finance. The parameter $q$ controls dependence. If $q=1$ then they are classical i.i.d. with $α$-stable Lévy distributions. In the present paper we establish basic properties of $(q,α)$-stable distributions, and generalize the result of Umarov, Tsallis and Steinberg (2008), where the particular case $α=2, q\in [1,3),$ was considered, to the whole range of stability and nonextensivity parameters $α\in (0,2]$ and $q \in [1,3),$ respectively. We also discuss possible further extensions of the results that we obtain, and formulate some conjectures.

preprint2010arXiv

Limit distribution in the $q$-CLT for $q \ge 1$ can not have a compact support

In a recent paper Hilhorst \cite{Hilhorst2010} illustrated that the $q$-Fourier transform for $q>1$ is not invertible in the space of density functions. Using an invariance principle he constructed a family of densities with the same $q$-Fourier transform and claimed that $q$-Gaussians are not mathematically proved to be attractors. We show here that none of the distributions constructed in Hilhorst's counterexamples can be a limit distribution in the $q$-CLT, except the one whose support covers the whole real axis, which is precisely the $q$-Gaussian distribution.

preprint2009arXiv

Limit distributions of scale-invariant probabilistic models of correlated random variables with the q-Gaussian as an explicit example

Extremization of the Boltzmann-Gibbs (BG) entropy under appropriate norm and width constraints yields the Gaussian distribution. Also, the basic solutions of the standard Fokker-Planck (FP) equation (related to the Langevin equation with additive noise), as well as the Central Limit Theorem attractors, are Gaussians. The simplest stochastic model with such features is N to infinity independent binary random variables, as first proved by de Moivre and Laplace. What happens for strongly correlated random variables? Such correlations are often present in physical situations as e.g. systems with long range interactions or memory. Frequently q-Gaussians become observed. This is typically so if the Langevin equation includes multiplicative noise, or the FP equation to be nonlinear. Scale-invariance, i.e. exchangeable binary stochastic processes, allow a systematical analysis of the relation between correlations and non-Gaussian distributions. In particular, a generalized stochastic model yielding q-Gaussians for all q (including q>1) was missing. This is achieved here by using the Laplace-de Finetti representation theorem, which embodies strict scale-invariance of interchangeable random variables. We demonstrate that strict scale invariance together with q-Gaussianity mandates the associated extensive entropy to be BG.

preprint2009arXiv

Nonadditive entropy and nonextensive statistical mechanics - Some central concepts and recent applications

We briefly review central concepts concerning nonextensive statistical mechanics, based on the nonadditive entropy $S_q=k\frac{1-\sum_{i}p_i^q}{q-1} (q \in {\cal R}; S_1=-k\sum_{i}p_i \ln p_i)$. Among others, we focus on possible realizations of the $q$-generalized Central Limit Theorem, including at the edge of chaos of the logistic map, and for quasi-stationary states of many-body long-range-interacting Hamiltonian systems.

preprint2008arXiv

A closer look at time averages of the logistic map at the edge of chaos

The probability distribution of sums of iterates of the logistic map at the edge of chaos has been recently shown [see U. Tirnakli, C. Beck and C. Tsallis, Phys. Rev. E 75, 040106(R) (2007)] to be numerically consistent with a q-Gaussian, the distribution which, under appropriate constraints, maximizes the nonadditive entropy S_q, the basis of nonextensive statistical mechanics. This analysis was based on a study of the tails of the distribution. We now check the entire distribution, in particular its central part. This is important in view of a recent q-generalization of the Central Limit Theorem, which states that for certain classes of strongly correlated random variables the rescaled sum approaches a q-Gaussian limit distribution. We numerically investigate for the logistic map with a parameter in a small vicinity of the critical point under which conditions there is convergence to a q-Gaussian both in the central region and in the tail region, and find a scaling law involving the Feigenbaum constant delta. Our results are consistent with a large number of already available analytical and numerical evidences that the edge of chaos is well described in terms of the entropy S_q and its associated concepts.

preprint2008arXiv

Escort mean values and the characterization of power-law-decaying probability densities

Escort mean values (or $q$-moments) constitute useful theoretical tools for describing basic features of some probability densities such as those which asymptotically decay like {\it power laws}. They naturally appear in the study of many complex dynamical systems, particularly those obeying nonextensive statistical mechanics, a current generalization of the Boltzmann-Gibbs theory. They recover standard mean values (or moments) for $q=1$. Here we discuss the characterization of a (non-negative) probability density by a suitable set of all its escort mean values together with the set of all associated normalizing quantities, provided that all of them converge. This opens the door to a natural extension of the well known characterization, for the $q=1$ instance, of a distribution in terms of the standard moments, provided that {\it all} of them have {\it finite} values. This question would be specially relevant in connection with probability densities having {\it divergent} values for all nonvanishing standard moments higher than a given one (e.g., probability densities asymptotically decaying as power-laws), for which the standard approach is not applicable. The Cauchy-Lorentz distribution, whose second and higher even order moments diverge, constitutes a simple illustration of the interest of this investigation. In this context, we also address some mathematical subtleties with the aim of clarifying some aspects of an interesting non-linear generalization of the Fourier Transform, namely, the so-called $q$-Fourier Transform.

preprint2008arXiv

Nonadditive entropy: the concept and its use

The entropic form $S_q$ is, for any $q \neq 1$, {\it nonadditive}. Indeed, for two probabilistically independent subsystems, it satisfies $S_q(A+B)/k=[S_q(A)/k]+[S_q(B)/k]+(1-q)[S_q(A)/k][S_q(B)/k] \ne S_q(A)/k+S_q(B)/k$. This form will turn out to be {\it extensive} for an important class of nonlocal correlations, if $q$ is set equal to a special value different from unity, noted $q_{ent}$ (where $ent$ stands for $entropy$). In other words, for such systems, we verify that $S_{q_{ent}}(N) \propto N (N>>1)$, thus legitimating the use of the classical thermodynamical relations. Standard systems, for which $S_{BG}$ is extensive, obviously correspond to $q_{ent}=1$. Quite complex systems exist in the sense that, for them, no value of $q$ exists such that $S_q$ is extensive. Such systems are out of the present scope: they might need forms of entropy different from $S_q$, or perhaps -- more plainly -- they are just not susceptible at all for some sort of thermostatistical approach. Consistently with the results associated with $S_q$, the $q$-generalizations of the Central Limit Theorem and of its extended Lévy-Gnedenko form have been achieved. These recent theorems could of course be the cause of the ubiquity of $q$-exponentials, $q$-Gaussians and related mathematical forms in natural, artificial and social systems. All of the above, as well as presently available experimental, observational and computational confirmations -- in high energy physics and elsewhere --, are briefly reviewed. Finally, we address a confusion which is quite common in the literature, namely referring to distinct physical mechanisms {\it versus} distinct regimes of a single physical mechanism.

preprint2004arXiv

Preferential attachment growth model and nonextensive statistical mechanics

We introduce a two-dimensional growth model where every new site is located, at a distance $r$ from the barycenter of the pre-existing graph, according to the probability law $1/r^{2+α_G} (α_G \ge 0)$, and is attached to (only) one pre-existing site with a probability $\propto k_i/r^{α_A}_i (α_A \ge 0$; $k_i$ is the number of links of the $i^{th}$ site of the pre-existing graph, and $r_i$ its distance to the new site). Then we numerically determine that the probability distribution for a site to have $k$ links is asymptotically given, for all values of $α_G$, by $P(k) \propto e_q^{-k/κ}$, where $e_q^x \equiv [1+(1-q)x]^{1/(1-q)}$ is the function naturally emerging within nonextensive statistical mechanics. The entropic index is numerically given (at least for $α_A$ not too large) by $q = 1+(1/3) e^{-0.526 α_A}$, and the characteristic number of links by $κ\simeq 0.1+0.08 α_A$. The $α_A=0$ particular case belongs to the same universality class to which the Barabasi-Albert model belongs. In addition to this, we have numerically studied the rate at which the average number of links $<k_i>$ increases with the scaled time $t/i$; asymptotically, $<k_i > \propto (t/i)^β$, the exponent being close to $β={1/2}(1-α_A)$ for $0 \le α_A \le 1$, and zero otherwise. The present results reinforce the conjecture that the microscopic dynamics of nonextensive systems typically build (for instance, in Gibbs $Γ$-space for Hamiltonian systems) a scale-free network.

preprint2004arXiv

Quasi-stationary states in low-dimensional Hamiltonian systems

We address a simple connection between results of Hamiltonian nonlinear dynamical theory and thermostatistics. Using a properly defined dynamical temperature in low-dimensional symplectic maps, we display and characterize long-standing quasi-stationary states that eventually cross over to a Boltzmann-Gibbs-like regime. As time evolves, the geometrical properties (e.g., fractal dimension) of the phase space change sensibly, and the duration of the anomalous regime diverges with decreasing chaoticity. The scenario that emerges is consistent with the nonextensive statistical mechanics one.

preprint2003arXiv

Nonextensive statistical mechanics: Some links with astronomical phenomena

A variety of astronomical phenomena appear to not satisfy the ergodic hypothesis in the relevant stationary state, if any. As such, there is no reason for expecting the applicability of Boltzmann-Gibbs (BG) statistical mechanics. Some of these phenomena appear to follow, instead, nonextensive statistical mechanics. In the same manner that the BG formalism is based on the entropy $S_{BG}=-k \sum_i p_i \ln p_i$, the nonextensive one is based on the form $S_q=k(1-\sum_ip_i^q)/(q-1)$ (with $S_1=S_{BG}$). The stationary states of the former are characterized by an {\it exponential} dependence on the energy, whereas those of the latter are characterized by an (asymptotic) {\it power-law}. A brief review of this theory is given here, as well as of some of its applications, such as the solar neutrino problem, polytropic self-gravitating systems, galactic peculiar velocities, cosmic rays and some cosmological aspects. In addition to these, an analogy with the Keplerian elliptic orbits {\it versus} the Ptolemaic epicycles is developed, where we show that optimizing $S_q$ with a few constraints is equivalent to optimizing $S_{BG}$ with an infinite number of constraints.

preprint2002arXiv

Constructing a statistical mechanics for Beck-Cohen superstatistics

The basic aspects of both Boltzmann-Gibbs (BG) and nonextensive statistical mechanics can be seen through three different stages. First, the proposal of an entropic functional ($S_{BG} =-k\sum_i p_i \ln p_i$ for the BG formalism) with the appropriate constraints ($\sum_i p_i=1$ and $\sum_i p_i E_i = U$ for the BG canonical ensemble). Second, through optimization, the equilibrium or stationary-state distribution ($p_i = e^{-βE_i}/Z_{BG}$ with $Z_{BG}=\sum_j e^{-βE_j}$ for BG). Third, the connection to thermodynamics (e.g., $F_{BG}= -\frac{1}β\ln Z_{BG}$ and $U_{BG}=-\frac{\partial}{\partial β} \ln Z_{BG}$). Assuming temperature fluctuations, Beck and Cohen recently proposed a generalized Boltzmann factor $B(E) = \int_0^\infty dβf(β) e^{-βE}$. This corresponds to the second stage above described. In this letter we solve the corresponding first stage, i.e., we present an entropic functional and its associated constraints which lead precisely to $B(E)$. We illustrate with all six admissible examples given by Beck and Cohen.

preprint2001arXiv

Mixing and equilibration: Protagonists in the scene of nonextensive statistical mechanics

After a brief review of the present status of nonextensive statistical mechanics, we present a conjectural scenario where mixing (characterized by the entropic index $q_{mix} \le 1$) and equilibration (characterized by the entropic index $q_{eq} \ge 1$) play central and inter-related roles, and appear to determine {\it a priori} the values of the relevant indices of the formalism. Boltzmann-Gibbs statistical mechanics is recovered as the $q_{mix}=q_{eq}=1$ particular case.

preprint1999arXiv

Functions of linear operators: Parameter differentiation

We derive a useful expression for the matrix elements $[\frac{\partial f[A(t)]}{\partial t}]_{i j}$ of the derivative of a function $f[A(t)]$ of a diagonalizable linear operator $A(t)$ with respect to the parameter $t$. The function $f[A(t)]$ is supposed to be an operator acting on the same space as the operator $A(t)$. We use the basis which diagonalizes A(t), i.e., $A_{i j}=λ_i δ_{i j}$, and obtain $[\frac{\partial f[A(t)]}{\partial t}]_{i j}=[\frac{\partial A}{\partial t}]_ {i j}\frac{f(λ_j) - f(λ_i)} {λ_j - λ_i}$. In addition to this, we show that further elaboration on the (not necessarily simple) integral expressions given by Wilcox 1967 (who basically considered $f[A(t)]$ of the exponential type) and generalized by Rajagopal 1998 (who extended Wilcox results by considering $f[A(t)]$ of the $q$-exponential type where $\exp_q(x) \equiv [1+(1-q)x]^{1/(1-q)}$ with $q \in {\cal {R}}$; hence, $\exp_1 (x)=\exp(x))$ yields this same expression. Some of the lemmas first established by the above authors are easily recovered.

preprint1995arXiv

Generalized Simulated Annealing

We propose a new stochastic algorithm (generalized simulated annealing) for computationally finding the global minimum of a given (not necessarily convex) energy/cost function defined in a continuous D-dimensional space. This algorithm recovers, as particular cases, the so called classical ("Boltzmann machine") and fast ("Cauchy machine") simulated annealings, and can be quicker than both. Key-words: simulated annealing; nonconvex optimization; gradient descent; generalized statistical mechanics.