Source author record

Constantino Tsallis

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

61works
19topics
4close collaborators

Published work

61 publications

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 SS on the quantum Grüneisen parameter Γq0KΓ^{0\text{K}}_{q} right at CPs for the 1D Ising model under a transverse magnetic field. Our findings are fourfold: i\textit{i}) for higher SS, Γq0KΓ^{0\text{K}}_{q} is increased, but remains finite, reflecting the enhancement of the Hilbert space dimensionality; ii\textit{ii}) the Schmidt decomposition theorem recovers the extensivity of the nonadditive qq-entropy SqS_q only for a special\textit{special} value of the entropic index qq; iii\textit{iii}) the universality class in the frame of SqS_q depends only on the symmetry of the system; iv\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 Γq0KΓ^{0\text{K}}_{q} and SqS_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 λ→1λ\to 1 limit of the λλ-state Potts ferromagnet with bond percolation, (ii) the isomorphism which connects the λ→0λ\to 0 limit of the λλ-state Potts ferromagnet with random resistor networks, and (iii) the de Gennes isomorphism which connects the n→0n \to 0 limit of the nn-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 qq-exponential distribution ∝eq−βqε\propto e_q^{-β_q \varepsilon} (with q=4/3q=4/3 and βqω0=10/3β_q ω_0 =10/3) optimizing, under simple constraints, the nonadditive entropy SqS_q with a specific geographic growth random model based on preferential attachment through exponentially-distributed weighted links, ω0ω_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 dd-dimensional geographic networks (characterized by the index αG≥0α_G\geq0; d=1,2,3,4d = 1, 2, 3, 4) whose links are weighted through a predefined random probability distribution, namely P(w)∝e−∣w−wc∣/τP(w) \propto e^{-|w - w_c|/τ}, ww being the weight (wc≥0;  τ>0) (w_c \geq 0; \; τ> 0). In this model, each site has an evolving degree kik_i and a local energy εi≡∑j=1kiwij/2\varepsilon_i \equiv \sum_{j=1}^{k_i} w_{ij}/2 (i=1,2,...,Ni = 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∝εi/dij αA    (αA≥0)Π_{ij}\propto \varepsilon_{i}/d^{\,α_A}_{ij} \;\;(α_A \ge 0), where dijd_{ij} is the Euclidean distance between the sites. Short- and long-range interactions respectively correspond to αA/d>1α_A/d>1 and 0≤αA/d≤10\leq α_A/d \leq 1; αA/d→∞α_A/d \to \infty corresponds to interactions between close neighbors, and αA/d→0α_A/d \to 0 corresponds to infinitely-ranged interactions. The site energy distribution p(ε)p(\varepsilon) corresponds to the usual degree distribution p(k)p(k) as the particular instance (wc,τ)=(2,0)(w_c,τ)=(2,0). We numerically verify that the corresponding connectivity distribution p(ε)p(\varepsilon) converges, when αA/d→∞α_A/d\to\infty, to the weight distribution P(w)P(w) for infinitely narrow distributions (i.e., τ→∞, ∀wcτ\to \infty, \,\forall w_c) as well as for wc→0, ∀τw_c\to0, \, \forallτ.

preprint2021arXiv

Along the lines of nonadditive entropies: qq-prime numbers and qq-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)≡∑n=1∞n−s=∏p prime11−p−sζ(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)ζ(s) to the complex plane zz and conjectured that all nontrivial zeros are in the R(z)=1/2\mathbb{R}(z)=1/2 axis. The nonadditive entropy Sq=k∑ipiln⁡q(1/pi)  (q∈R)S_q=k\sum_ip_i\ln_q(1/p_i)\;(q\in\mathbb{R}) involves the function ln⁡qz≡z1−q−11−q  (ln⁡1z=ln⁡z)\ln_q z\equiv\frac{z^{1-q}-1}{1-q}\;(\ln_1 z=\ln z) that is interconnected to a qq-generalized algebra, using qq-numbers defined as ⟨x⟩q≡eln⁡qx\langle x\rangle_q\equiv e^{\ln_q x} (⟨x⟩1≡x\langle x\rangle_1\equiv x). The qq-prime numbers are then defined as the qq-natural numbers ⟨n⟩q≡eln⁡qn  (n=1,2,3,… )\langle n\rangle_q\equiv e^{\ln_q n}\;(n=1,2,3,\dots), where nn is a prime number p=2,3,5,7,…p=2,3,5,7,\dots We show that, for any value of qq, infinitely many qq-prime numbers exist; for q≤1q\le1 they diverge for increasing prime number, whereas they converge for q>1q>1; the standard prime numbers are recovered for q=1q=1. For q≤1q\le 1, we generalize the ζ(s)ζ(s) function as follows: ζq(s)≡⟨ζ(s)⟩qζ_q(s)\equiv\langleζ(s)\rangle_q (s∈Rs\in\mathbb{R}). We show that this function appears to diverge at s=1+0s=1+0, ∀q\forall q. Also, we alternatively define, for q≤1q\le 1, ζqΣ(s)≡∑n=1∞1⟨n⟩qs=1+1⟨2⟩qs+…ζ_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)≡∏p prime11−⟨p⟩q−s=11−⟨2⟩q−s11−⟨3⟩q−s11−⟨5⟩q−s⋯ζ_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<1q<1, generically satisfy ζqΣ(s)<ζqΠ(s)ζ_q^Σ(s)<ζ_q^Π(s), in variance with the q=1q=1 case, where of course ζ1Σ(s)=ζ1Π(s)ζ_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 QQ-Gaussian N→∞N\to\infty attractor in the space of distributions (1<Q<3)1 < Q<3). We present strong numerical indications that the corresponding LDT probability distribution is given by P(N,z)=P0 eq−rq(z)N=P0[1−(1−q)rq(z)N]1/(1−q)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∈(1,5/3)q=2-1/Q \in (1,5/3). The rate function rq(z)r_q(z) seemingly equals the relative nonadditive qrq_r-entropy per particle, with qr≃710+6101Q−1q_r \simeq \frac{7}{10} + \frac{6}{10}\frac{1}{Q-1}, thus exhibiting a singularity at Q=1Q=1 and recovering the BG value qr=1q_r=1 in the Q→3Q \to 3 limit. Let us emphasize that the extensivity of rq(z)Nr_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 −∞<q≤3-\infty<q\leq3 is related to the shape of the tail decay; q<1q<1 for compact-support including parabola (q=0)(q=0); 1<q≤31<q\leq3 for heavy-tail including Cauchy (q=2)(q=2). This addendum clarifies the transformation q′=((3q−1)/(q+1))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κ)q=1+2κ/(1+dκ), where the shape parameter −1<κ≤∞-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 qq-entropies adequately handles a wide class of such systems. We show here that scale-invariant networks belong to this class. We numerically study a dd-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 qq-generalisation, and is recovered in the q=1q=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 qq-Gaussian (with q≃1.935q \simeq 1.935), consistent with qq-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∈[0,1]η_i \in [0,1] for each node ii 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∝kiηirij−αAΠ_i\propto k_i η_i r_{ij}^{-α_A} where kik_i is the connectivity of the iith pre-existing site and αAα_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 SqS_q, given by p(k)∝eq−k/κ≡1/[1+(q−1)k/κ]1/(q−1)p(k) \propto e_q^{-k/κ} \equiv 1/[1+(q-1)k/κ]^{1/(q-1)}, with (q,κ)(q,κ) depending uniquely only on the ratio αA/dα_A/d and the fitness distribution. Hence this model constitutes a realization of asymptotically-scale-free geographical networks within nonextensive statistical mechanics, where kk 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 WW microscopic possibilities, and is given by SBG=−k∑i=1Wpiln⁡piS_{BG}= -k\sum_{i=1}^W p_i \ln p_i (kk 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 Sq=k1−∑i=1Wpiqq−1S_q=k\frac{1-\sum_{i=1}^W p_i^q}{q-1}, which recovers the additive SBGS_{BG} entropy in the q→1q \to1 limit. The index qq 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α_1 and α2α_2 respectively, which make the {forces decay} with distance rr. Our results demonstrate that weak chaos, in the sense of decreasing Lyapunov exponents, and qq--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≤α2<10\leq α_2<1, we obtain extrapolated values for q≡q∞>1q \equiv q_\infty >1, as N→∞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 qEq_E-exponentials (with qE>1q_E>1) when the quartic-term interactions are long--ranged, otherwise we get the standard Boltzmann-Gibbs weight, with q=1q=1. The values of qEq_E coincide, within small discrepancies, with the values of qq 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 SBGS_{BG} and its associated exponential BG factor. Frequent behaviors in such complex systems have been shown to be closely related to qq-statistics instead, based on the nonadditive entropy SqS_q (with S1=SBGS_1=S_{BG}), and its associated qq-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 qq. 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 qq-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 qq-Gaussian distribution (which extremizes the qq-entropy). We give its explicit expression for q=0q=0, and discuss the other cases numerically. The entropy that we obtain can be expressed, for q=0q=0, in terms of Lambert functions, and exhibits a discontinuity in the second derivative for all values of q<1q<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 dd-dimensional long-range-interacting quartic Fermi-Pasta-Ulam model: Universal scaling

We introduce a generalized dd-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,3d=1,2,3 through large-scale numerical simulations. The nonlinear interaction is assumed to decay algebraically as dij−αd_{ij}^{-α} (α≥0α\ge 0), {dij}\{d_{ij}\} being the distances between NN oscillator sites. Starting from random initial conditions we compute the maximal Lyapunov exponent λmaxλ_{max} as a function of NN. Our N>>1N>>1 results strongly indicate that λmaxλ_{max} remains constant and positive for α/d>1α/d>1 (implying strong chaos, mixing and ergodicity), and that it vanishes like N−κN^{-κ} for 0≤α/d<10 \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)κ(d+2) κ depends only on α/dα/d and, when α/dα/d increases from zero to unity, it monotonically decreases from unity to zero, remaining so for all α/d>1α/d >1. The value α/d=1α/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α/d > 1, and possibly with qq-statistics for 0≤α/d<10 \le α/d < 1.

preprint2016arXiv

Statistical mechanics for complex systems: On the structure of qq-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 Sq≡k1−∑i=1Wpiqq−1  (q∈R; ∑i=1Wpi=1)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 S1=SBG≡−k∑i=1Wpiln⁡piS_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}\{q\}, directly corresponding to its various physical/dynamical/geometrical properties. A most important subset is usually referred to as the qq-triplet, namely (qsensitivity,qrelaxation,qstationary state)(q_{sensitivity}, q_{relaxation}, q_{stationary\,state}), defined in the body of this paper. In the BG limit we have qsensitivity=qrelaxation=qstationary state=1q_{sensitivity}=q_{relaxation}=q_{stationary\,state}=1. For a given class of complex systems, the set {q}\{q\} contains only a few independent values of qq, 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 qq-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 qq-CLT for q≥1q \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 qq-Fourier transform, was introduced and applied for the proof of a qq-generalized central limit theorem (qq-CLT). Subsequently, Hilhorst illustrated (2009 and 2010) that the qq-Fourier transform for q>1q>1 is not invertible in the space of density functions. Indeed, using an invariance principle, he constructed a family of densities with the same qq-Fourier transform and noted that "as a consequence, the qq-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 qq-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 NN 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σ_{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 qq-Gaussians, the distributions which optimise, under appropriate constraints, the nonadditive entropy SqS_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 qq-Gaussians. Nevertheless, along many decades, the fitting with qq-Gaussians turns out to be numerically very satisfactory for wide regions of the parameter values, and we illustrate how the index qq evolves with (N,τ,ε,σmax)(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 SqS_{q} is introduced. This new entropic form, denoted by Sa,b,rS_{a,b,r}, possesses many interesting statistical properties, and it reduces to the entropy SqS_q for b=0b=0, a=r:=1−qa=r:=1-q (hence Boltzmann-Gibbs entropy SBGS_{BG} for b=0b=0, a=r→0a=r \to 0). The construction of the entropy Sa,b,rS_{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 Sa,b,rS_{a,b,r} with respect to the composition of statistically independent subsystems. Depending on the choice of the parameters, the entropy Sa,b,rS_{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 NN of the system, or even stabilizes, by increasing NN, 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 (SBGS_{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 SBG=−k∑ipiln⁡piS_{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 Sq=k1−∑ipiqq−1S_q=k \frac{1-\sum_i p_i^q}{q-1} one (q∈R; S1=SBGq \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 SqS_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 QQ-Gaussians as limiting distributions

We consider correlated random variables X1,…,XnX_1,\dots,X_n taking values in {0,1}\{0,1\} such that, for any permutation ππ of {1,…,n}\{1,\dots,n\}, the random vectors (X1,…,Xn)(X_1,\dots,X_n) and (Xπ(1),…,Xπ(n))(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ν>0 (ν→∞ν\to\infty implies independence). Putting Sn=X1+⋯+XnS_n=X_1+\cdots+X_n, the distribution of Sn−n/2S_n-n/2 approaches a QQ-Gaussian distribution with compact support (Q=1−1/(ν−1)<1Q=1-1/(ν-1)<1) as nn increases, after appropriate scaling. In the present article, we show that the distribution of Sn/nS_n/n converges, as n→∞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≤x<1/20\le x<1/2, then P(Sn/n≤x)\mathbb{P}(S_n/n\le x), which is the probability of the event {Sn/n≤x}\{S_n/n\le x\} (large deviation), does not converges to zero as n→∞n\to\infty. For x=0x=0 and every real ν>0ν>0, we show that P(Sn=0)\mathbb{P}(S_n=0) decays to zero like a power law of the form 1/nν1/n^ν with a subdominant term of the form 1/nν+11/n^{ν+1}. If 0<x≤10<x\le 1 and ν>0ν>0 is an integer, we show that we can analytically find upper and lower bounds for the difference between P(Sn/n≤x)\mathbb{P}(S_n/n\le x) and its (n→∞n\to\infty) limit. We also show that these bounds vanish like a power law of the form 1/n1/n with a subdominant term of the form 1/n21/n^2.

preprint2015arXiv

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

Transverse spectra of both jets and hadrons obtained in high-energy pppp and ppˉp\bar p collisions at central rapidity exhibit power-law behavior of 1/pTn1/p_T^n at high pTp_T. The power index nn 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 pTp_T region in pppp 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 pTp_T region at central rapidity in high-energy pppp and ppˉ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 pTp_T region and the two-particle correlation data of the STAR and PHENIX Collaborations in high-energy pppp and ppˉ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 pppp and ppˉp\bar p collisions.

preprint2015arXiv

On the robustness of the qq-Gaussian family

We introduce three deformations, called αα-, ββ- and γγ-deformation respectively, of a NN-body probabilistic model, first proposed by Rodríguez et al. (2008), having qq-Gaussians as N→∞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 qq-Gaussians as limiting distributions, with a value of qq independent (dependent) on the deformation parameter in the αα-case (ββ-case). In contrast, the γγ-case, where we have used the celebrated QQ-numbers and the Gauss binomial coefficients, yields other limiting probability distribution functions, outside the qq-Gaussian family. These results suggest that scale-invariance might play an important role regarding the robustness of the qq-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)∝eq−k/κP(k) \propto e_q^{-k/κ}, where the qq-exponential form eqz≡[1+(1−q)z]11−qe_q^z \equiv [1+(1-q)z]^{\frac{1}{1-q}} optimizes the nonadditive entropy SqS_q (which, for q→1q\to 1, recovers the Boltzmann-Gibbs entropy). We introduce and study here dd-dimensional geographically-located networks which grow with preferential attachment involving Euclidean distances through rij−αA  (αA≥0)r_{ij}^{-α_A} \; (α_A \ge 0). Revealing the connection with qq-statistics, we numerically verify (for dd =1, 2, 3 and 4) that the qq-exponential degree distributions exhibit, for both qq and κκ, universal dependences on the ratio αA/dα_A/d. Moreover, the q=1q=1 limit is rapidly achieved by increasing αA/dα_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 qq-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 qq-statistical functions in a generalized binomial distribution with strong correlations

We study a symmetric generalization pk(N)(η,α)\mathfrak{p}^{(N)}_k(η, α) of the binomial distribution recently introduced by Bergeron et al, where η∈[0,1]η\in [0,1] denotes the win probability, and αα is a positive parameter. This generalization is based on qq-exponential generating functions (eqgenz≡[1+(1−qgen)z]1/(1−qgen); e1z=ez)e_{q^{gen}}^z \equiv [1+(1-q^{gen})z]^{1/(1-q^{gen})};\,e_{1}^z=e^z) where qgen=1+1/αq^{gen}=1+1/α. The numerical calculation of the probability distribution function of the number of wins kk, related to the number of realizations NN, strongly approaches a discrete qdiscq^{disc}-Gaussian distribution, for win-loss equiprobability (i.e., η=1/2η=1/2) and all values of αα. Asymptotic N→∞N\to \infty distribution is in fact a qattq^{att}-Gaussian eqatt−βz2e_{q^{att}}^{-βz^2}, where qatt=1−2/(α−2)q^{att}=1-2/(α-2) and β=(2α−4)β=(2α-4). The behavior of the scaled quantity k/Nγk/N^γ is discussed as well. For γ<1γ<1, a large-deviation-like property showing a qldlq^{ldl}-exponential decay is found, where qldl=1+1/(ηα)q^{ldl}=1+1/(ηα). For η=1/2η=1/2, qldlq^{ldl} and qattq^{att} are related through 1/(qldl−1)+1/(qatt−1)=11/(q^{ldl}-1)+1/(q^{att}-1)=1, ∀α\forall α. For γ=1γ=1, the law of large numbers is violated, and we consistently study the large-deviations with respect to the probability of the N→∞N\to\infty limit distribution, yielding a power law, although not exactly a qLDq^{LD}-exponential decay. All qq-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α1/r^α (α≥0α\ge 0)(with strength characterized by b>0b>0). In the α→∞α\to\infty limit we recover the original FPU model. Through classical molecular dynamics computations we show that (i) For α≥1α\geq 1 the maximal Lyapunov exponent remains finite and positive for increasing number of oscillators NN (thus yielding ergodicity), whereas, for 0≤α<10 \le α<1, it asymptotically decreases as N−κ(α)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 qq-Gaussian, with the index q(α)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 tct_c from qq-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∝bδ/tcγ1/N \propto b^δ/t_c^γ with γ>0γ>0 and δ>0δ>0, in such a way that the q=1q=1 (q>1q>1) behavior dominates in the lim⁡N→∞lim⁡t→∞\lim_{N \to\infty} \lim_{t \to\infty} ordering (lim⁡t→∞lim⁡N→∞\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α=0, 1/21/2 and 11 respectively correspond to the Itô, Stratonovich and anti-Itô prescriptions). We obtain its stationary state pst(x)p_{st}(x) for a class of constitutive relations between drift and diffusion and show that it has a qq-exponential form, pst(x)=Nq[1−(1−q)βV(x)]1/(1−q)p_{st}(x) = N_q[1 - (1-q)βV(x)]^{1/(1-q)}, with an index qq which does not depend on αα in the presence of any nonvanishing nonlinearity. This is in contrast with the linear case, for which the index qq is αα-dependent.

preprint2014arXiv

Possible Implication of a Single Nonextensive pTp_T Distribution for Hadron Production in High-Energy pppp Collisions

Multiparticle production processes in pppp 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-pTp_T region with a scale parameter TT associated with the temperature of the hadronizing system. The second is governed by a power-like distributions of transverse momenta with power index nn at high-pTp_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 TT and a power index n=1/(q−1)n=1/(q -1), where qq 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 1414 orders of magnitude down to the low pTp_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 pTp_T spectra for the whole pTp_T region at central rapidity for pppp collisions at high-energies.

preprint2014arXiv

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

We analyze LHC available data measuring the distribution probability of transverse momenta~pTp_T in proton-proton collisions at s=0.9 TeV\sqrt{s}=0.9\,\textrm{TeV} (CMS, ALICE, ATLAS) and s=7 TeV\sqrt{s}=7\,\textrm{TeV} (CMS, ATLAS). A remarkably good fitting can be obtained, along fourteen decades in magnitude, by phenomenologically using qq-statistics for a {\it single} particle of a two-dimensional relativistic ideal gas. The parameters that have been obtained by assuming dN/pTdpTdy∝eq−ET/T\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≃1.1q \simeq 1.1 and T≃0.13 GeVT\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 qq-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→∞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.}, qq-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 SBGS_{BG} of a (3+1)(3+1) black hole is proportional to its area L2L^2 (LL being a characteristic linear length), and not to its volume L3L^3. Similarly it exists the \emph{area law}, so named because, for a wide class of strongly quantum-entangled dd-dimensional systems, SBGS_{BG} is proportional to ln⁡L\ln L if d=1d=1, and to Ld−1L^{d-1} if d>1d>1, instead of being proportional to LdL^d (d≥1d \ge 1). These results violate the extensivity of the thermodynamical entropy of a dd-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 NN classical localized planar rotators coupled through interactions which decay with distance as 1/rα1/r^α (α≥0α\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≫1N\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 qq-Gaussians. These distributions extremize, under appropriate simple constraints, the nonadditive entropy SqS_q upon which nonextensive statistical mechanics is based. The whole scenario appears to be consistent with nonergodicity and with the thesis of the qq-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 qq-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 eqi(kx−wt)e_q^{i(kx-wt)}, involving the qq-exponential function which naturally emerges within nonextensive thermostatistics [eqz≡[1+(1−q)z]1/(1−q)e_q^z \equiv [1+(1-q)z]^{1/(1-q)}, with e1z=eze_1^z=e^z]. Since these basic solutions behave like free particles, obeying p=ℏkp=\hbar k, E=ℏωE=\hbar ω and E=p2/2mE=p^2/2m (1≤q<21 \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 aa) or a particle moving under a constant force −ma-ma. The latter interpretation naturally leads to the evolution equation iℏ∂∂t(ΦΦ0)=−12−qℏ22m∂2∂x2[(ΦΦ0)2−q]+V(x)(ΦΦ0)qi\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)=maxV(x)=max. Remarkably enough, the potential VV couples to ΦqΦ^q, instead of coupling to ΦΦ, as happens in the familiar linear case (q=1q=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 NN 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 NN; (ii) the subdominant term appears to be consistent with the qq-exponential behavior typical of systems following qq-statistics, thus reinforcing the thermodynamically extensive entropic nature of the exponent of the qq-exponential, basically NN times the qq-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 QQ-Gaussian distributions (in principle 1<Q<31 < Q < 3), which relevantly differ from (symmetric) Lévy distributions, with the unique exception of Cauchy-Lorentz distributions (which correspond to Q=2Q = 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 qq-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)(d+1)-dimensional hyperpyramids (d=1,2,3,...d=1,2,3,...), characterized by a parameter ν>0ν>0, whose value determines the degree of correlation between NN (d+1)(d+1)-valued random variables. There are (d+1)N(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)N1/(d+1)^N to occur. The sums of these NN  (d+1)\,(d+1)-valued random variables correspond to a d−d-dimensional probabilistic model, and generalizes a recently proposed one-dimensional (d=1d=1) model having q−q-Gaussians (with q=(ν−2)/(ν−1)q=(ν-2)/(ν-1) for ν∈[1,∞)ν\in [1,\infty)) as N→∞N\to\infty limit probability distributions for the sum of the NN 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−d-dimensional multinomial distribution is recovered for the sums, which approach a d−d-dimensional Gaussian distribution for N→∞N\to\infty. For any νν, the conditional distributions of the d−d-dimensional model are shown to yield the corresponding joint distribution of the (d−1)(d-1)-dimensional model with the same νν. For the d=2d=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 qq-Gaussian attractors in the N→∞N\to\infty limit, the ultimate goal being a neat mathematical view of the causes clarifying the ubiquitous emergence of qq-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 ∂p(x,t)∂t=−∂∂x[F(x)p(x,t)]+1/2D∂2∂x2[ϕ(x,p)p(x,t)]\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∈RD \in {\cal R} and F(x)=−∂V(x)/∂xF(x)=-\partial V(x) /\partial x, V(x)V(x) being the potential under which diffusion occurs. Our aim is to find whether ϕ(x,p)ϕ(x,p) exists such that the inhomogeneous linear and the homogeneous nonlinear FP equations become unified in such a way that the (ubiquitously observed) qq-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)]ηϕ(x,p)=[A+BV(x)]^θ[p(x,t)]^η, where AA, BB, θθ and ηη are (real) constants. Our main result can be sumarized as follows: For θ≠1θ\neq 1 and arbitrary confining potential V(x)V(x), p(x,∞)∝{1−β(1−q)V(x)}1/(1−q)≡eq−β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)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 qq, are considered in such a way that the standard linear equations are recovered in the limit q→1q \rightarrow 1. Interestingly, these equations present a common, soliton-like, travelling solution, which is written in terms of the qq-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 qq.

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=∑t=1Txty=\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∗∝L2T^*\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 aa 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≠1a\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 zz-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 SqS_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} SS. This thermodynamically crucial quantity was proposed to be {\it extensive}, i.e., in contemporary terms, S(N)∝NS(N) \propto N in the thermodynamic limit N→∞N \to\infty. A decade later, Boltzmann proposed a functional form for this quantity which connects SS 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 WW possible discrete states, SBG=−kB∑i=1Wpiln⁡piS_{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 SBG(A+B)=SBG(A)+SBG(B)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 SBGS_{BG} was considered to be the unique form that SS could take. In other words, the functional SBGS_{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 SqS_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 SBG=−kB∑i=1Wpiln⁡piS_{BG}=- k_B\sum_{i=1}^W p_i \ln p_i. Its optimization under appropriate constraints yields the celebrated BG weight e−βEie^{-βE_i}. An elementary large-deviation connection is provided by NN independent binary variables, which, in the N→∞N\to\infty limit yields a Gaussian distribution. The probability of having n≠N/2n \ne N/2 out of NN throws is governed by the exponential decay e−Nre^{-N r}, where the rate function rr 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 Sq=kB1−∑i=1Wpiqq−1S_q=k_B\frac{1- \sum_{i=1}^W p_i^q}{q-1} (q∈R; S1=SBGq \in {\cal R}; \,S_1=S_{BG}). Its optimization yields the generalized weight eq−βqEie_q^{-β_q E_i} (eqz≡[1+(1−q)z]1/(1−q); e1z=ez)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∈[1,2)Q \in [1,2) and γ∈(0,1)γ\in (0,1). This model provides, in the N→∞N\to\infty limit (∀γ\forall γ), QQ-Gaussian distributions, ubiquitously observed in nature (Q→1Q\to 1 recovers the independent binary model). We show that its corresponding large deviations are governed by eq−Nrqe_q^{-N r_q} (∝1/N1/(q−1)\propto 1/N^{1/(q-1)} if q>1q>1) where q=Q−1γ(3−Q)+1≥1q= \frac{Q-1}{γ(3-Q)}+1 \ge 1. This qq-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>1q>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,α)(q,α)-stable distributions. These sequences are generalizations of i.i.d. αα-stable distributions, and have not been previously studied. Long-range dependent (q,α)(q,α)-stable distributions might arise in the description of anomalous processes in nonextensive statistical mechanics, cell biology, finance. The parameter qq controls dependence. If q=1q=1 then they are classical i.i.d. with αα-stable Lévy distributions. In the present paper we establish basic properties of (q,α)(q,α)-stable distributions, and generalize the result of Umarov, Tsallis and Steinberg (2008), where the particular case α=2,q∈[1,3),α=2, q\in [1,3), was considered, to the whole range of stability and nonextensivity parameters α∈(0,2]α\in (0,2] and q∈[1,3),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 qq-CLT for q≥1q \ge 1 can not have a compact support

In a recent paper Hilhorst \cite{Hilhorst2010} illustrated that the qq-Fourier transform for q>1q>1 is not invertible in the space of density functions. Using an invariance principle he constructed a family of densities with the same qq-Fourier transform and claimed that qq-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 qq-CLT, except the one whose support covers the whole real axis, which is precisely the qq-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 Sq=k1−∑ipiqq−1(q∈R;S1=−k∑ipiln⁡pi)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 qq-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 qq-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=1q=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=1q=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 qq-Fourier Transform.

preprint2008arXiv

Nonadditive entropy: the concept and its use

The entropic form SqS_q is, for any q≠1q \neq 1, {\it nonadditive}. Indeed, for two probabilistically independent subsystems, it satisfies Sq(A+B)/k=[Sq(A)/k]+[Sq(B)/k]+(1−q)[Sq(A)/k][Sq(B)/k]≠Sq(A)/k+Sq(B)/kS_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 qq is set equal to a special value different from unity, noted qentq_{ent} (where entent stands for entropyentropy). In other words, for such systems, we verify that Sqent(N)∝N(N>>1)S_{q_{ent}}(N) \propto N (N>>1), thus legitimating the use of the classical thermodynamical relations. Standard systems, for which SBGS_{BG} is extensive, obviously correspond to qent=1q_{ent}=1. Quite complex systems exist in the sense that, for them, no value of qq exists such that SqS_q is extensive. Such systems are out of the present scope: they might need forms of entropy different from SqS_q, or perhaps -- more plainly -- they are just not susceptible at all for some sort of thermostatistical approach. Consistently with the results associated with SqS_q, the qq-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 qq-exponentials, qq-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 rr from the barycenter of the pre-existing graph, according to the probability law 1/r2+αG(αG≥0)1/r^{2+α_G} (α_G \ge 0), and is attached to (only) one pre-existing site with a probability ∝ki/riαA(αA≥0\propto k_i/r^{α_A}_i (α_A \ge 0; kik_i is the number of links of the ithi^{th} site of the pre-existing graph, and rir_i its distance to the new site). Then we numerically determine that the probability distribution for a site to have kk links is asymptotically given, for all values of αGα_G, by P(k)∝eq−k/κP(k) \propto e_q^{-k/κ}, where eqx≡[1+(1−q)x]1/(1−q)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α_A not too large) by q=1+(1/3)e−0.526αAq = 1+(1/3) e^{-0.526 α_A}, and the characteristic number of links by κ≃0.1+0.08αAκ\simeq 0.1+0.08 α_A. The αA=0α_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 <ki><k_i> increases with the scaled time t/it/i; asymptotically, <ki>∝(t/i)β<k_i > \propto (t/i)^β, the exponent being close to β=1/2(1−αA)β={1/2}(1-α_A) for 0≤αA≤10 \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 SBG=−k∑ipiln⁡piS_{BG}=-k \sum_i p_i \ln p_i, the nonextensive one is based on the form Sq=k(1−∑ipiq)/(q−1)S_q=k(1-\sum_ip_i^q)/(q-1) (with S1=SBGS_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 SqS_q with a few constraints is equivalent to optimizing SBGS_{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 (SBG=−k∑ipiln⁡piS_{BG} =-k\sum_i p_i \ln p_i for the BG formalism) with the appropriate constraints (∑ipi=1\sum_i p_i=1 and ∑ipiEi=U\sum_i p_i E_i = U for the BG canonical ensemble). Second, through optimization, the equilibrium or stationary-state distribution (pi=e−βEi/ZBGp_i = e^{-βE_i}/Z_{BG} with ZBG=∑je−βEjZ_{BG}=\sum_j e^{-βE_j} for BG). Third, the connection to thermodynamics (e.g., FBG=−1βln⁡ZBGF_{BG}= -\frac{1}β\ln Z_{BG} and UBG=−∂∂βln⁡ZBGU_{BG}=-\frac{\partial}{\partial β} \ln Z_{BG}). Assuming temperature fluctuations, Beck and Cohen recently proposed a generalized Boltzmann factor B(E)=∫0∞dβf(β)e−βEB(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)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 qmix≤1q_{mix} \le 1) and equilibration (characterized by the entropic index qeq≥1q_{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 qmix=qeq=1q_{mix}=q_{eq}=1 particular case.

preprint1999arXiv

Functions of linear operators: Parameter differentiation

We derive a useful expression for the matrix elements [∂f[A(t)]∂t]ij[\frac{\partial f[A(t)]}{\partial t}]_{i j} of the derivative of a function f[A(t)]f[A(t)] of a diagonalizable linear operator A(t)A(t) with respect to the parameter tt. The function f[A(t)]f[A(t)] is supposed to be an operator acting on the same space as the operator A(t)A(t). We use the basis which diagonalizes A(t), i.e., Aij=λiδijA_{i j}=λ_i δ_{i j}, and obtain [∂f[A(t)]∂t]ij=[∂A∂t]ijf(λj)−f(λi)λj−λi[\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)]f[A(t)] of the exponential type) and generalized by Rajagopal 1998 (who extended Wilcox results by considering f[A(t)]f[A(t)] of the qq-exponential type where exp⁡q(x)≡[1+(1−q)x]1/(1−q)\exp_q(x) \equiv [1+(1-q)x]^{1/(1-q)} with q∈Rq \in {\cal {R}}; hence, exp⁡1(x)=exp⁡(x))\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.

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.