# Statistical characterization of the standard map

## Abstract

The standard map, paradigmatic conservative system in the phase space, has been recently shown [Tirnakli and Borges (2016)] to exhibit interesting statistical behaviors directly related to the value of the standard map external parameter . A detailed numerical description is achieved in the present paper. More precisely, for large values of (e.g., ), the Lyapunov exponents are neatly positive over virtually the entire phase space, and, consistently with Boltzmann-Gibbs (BG) statistics, we verify , where is the -index for which the nonadditive entropy (with ) grows linearly with time before achieving its -dependent saturation value; characterizes the time increase of the sensitivity to the initial conditions, i.e., , where ; is the index associated with the -Gaussian distribution of the time average of successive iterations of the -coordinate; finally, characterizes the -exponential relaxation with time of the entropy towards its saturation value. In remarkable contrast, for small values of (e.g., ), the Lyapunov exponents are virtually zero over the entire phase space, and, consistently with -statistics, we verify , , and . The situation corresponding to intermediate values of , where both stable orbits and a chaotic sea are present, is discussed as well. The present results transparently illustrate when BG behavior or -statistical behavior are observed.

###### pacs:

05.10.-a , 05.45.Ac, 05.45.Pq## I Introduction

Statistical characterization of dynamical systems, particularly weakly chaotic low-dimensional dissipative maps, have been made according to its sensitivity to the initial conditions, the rate of increase of the entropy, the rate of relaxation to the attractor, and the distribution of a statistical variable Tsallis2009 (). For strongly chaotic systems, these behaviors present an exponential nature and are described by the Boltzmann-Gibbs (BG) statistical mechanical framework BeckSchlogl (). Weakly chaotic systems, that present zero largest Lyapunov exponent (LLE) , need a more general formalism to be properly handled. It has been observed that the nonextensive statistical mechanical framework is able to describe different features of various classes of complex systems Rapisarda-2004 (); Esquivel-2010 (); CMS-2010 (); Shao-2010 (); Sotolongo-2010 (); Mohanalin-2010 (); ALICE-2010 (); PHENIXbis (); PHENIX-2011 (); ALICE-2011 (); ATLAS-2011 (); Wong-2012 (); Marques-2013 (); Bogachev-2014 (); Combe-2015 (), including the logistic map Baldovin-Robledo-2002 (). If approaches zero, the sensitivity to the initial conditions and the rate of relaxation to the attractor follow -exponential behaviors, the distributions follow a -Gaussian form, and the rate of increase of entropy demands entropy, instead of BG entropy. All these behaviors are parameterized by the indices , , , , for sensitivity, entropy production, relaxation, and stationary distribution, respectively. It has been observed (and in some instances proved) that, for low-dimensional systems, , through the -generalized Pesin equality Tsallis2009 (). Consequently, three (possibly related) quantities remain to be addressed, that are currently referred to as the -triplet Tsallis2009 (). The -triplet has been experimentally observed in a variety of complex systems like in the ozone layer Ferri2010 () and in the Voyager I data for the solar wind Burlaga2013 (), among many others. Relations are expected to connect , , and , ultimately completely determined by the nonlinear dynamics of the system. The possible relations between the -triplet indices should in principle exhibit the fixed point (BG statistics). At least for low-dimensional systems, if one of them equals unity, the others also equal unity. This typically happens whenever . In fact, LLE plays a central role in the definition of which statistical framework is to be used in the description of nonlinear dynamical systems.

In the present work we address a low-dimensional conservative system, namely the standard map chirikov (); zaslavsky1 ():

(1) |

( and are taken as modulo ). Its phase space can exhibit strongly chaotic regions, characterized by positive local LLE , and weakly chaotic regions, with zero values of . The LLE has been numerically evaluated through the Benettin et al algorithm benettin (). It is calculated for each initial condition separately, as illustrated by Fig. 2 of tirnakli-borges-2016 () (see also BarangerLatoraRapisarda2002 ()). The extent of the phase space for each case (strongly or weakly chaotic) is determined by the control parameter . For sufficiently small values of , for instance , the entire phase space is practically stable (zero LLE), and if takes high values (), the points of the phase space have finite positive values of . In the present work, we evaluate the values of , , and for both strongly and weakly chaotic cases. We also revisit the stationary distributions, particularly the mixed case (strongly and weakly chaotic), that has been recently found for this system tirnakli-borges-2016 (). Finally, we have determined the relation of strongly and weakly chaotic indices, with respect to the stationary distributions of the map for arbitrary values of .

## Ii Sensitivity to initial conditions

The sensitivity to initial conditions is given by

(2) |

where is the temporal dependence of the discrepancy of two very close initial conditions at time . When ergodic behavior dominates over the whole phase space, Eq. (2) is expressed by the exponential dependence,

(3) |

and is the standard maximum Lyapunov exponent. Strongly chaotic systems are those which present . They are strongly sensitive to initial conditions. Insensitivity to initial conditions are described by . The marginal case is governed by a more subtle rule than Eq. (3) tsallis-plastino-zheng-1997 ():

(4) |

where , , is the so called –exponential function tsallis:quimica_nova (), and is a generalized Lyapunov exponent. If (with ) the system is weakly chaotic, in the sense that diverges slower than the exponential (strongly chaotic) case, — asymptotically, a power law. If the system is strongly chaotic, , and Eq. (4) recovers the standard exponential dependence, with .

The correctness of this description can be checked, for both chaotic regimes of standard map, in Figures 1 and 4. They show, for different values of , the averages of (where is the inverse function of the -exponential, and ) over realizations. Each realization starts with a randomly chosen pair of very close initial conditions, that are localized inside a particular cell of the equally partitioned phase space. In Fig. 1, that corresponds to the case of strong chaos , where LLE ), we consider decreasing initial discrepancies in Eq. (2) so as to obtain a well defined behavior for increasingly long times. In all cases (strongly and weakly chaotic chaos), after a transient, we verify a nontrivial property ananos-tsallis-2004 () namely that exists a special value of , noted (where av stands for average), which yields a linear dependence of with time. In other words, we verify .

For the strongly chaotic case, exponential sensitivity to the initial conditions is verified with , and the same generalized Lyapunov exponent characterizes the whole phase space no mater the value of or the particular cell to be chosen. Fig. 2 exemplifies this fact, as it shows that the slope of the temporal evolution of the averaged sensitivity to initial conditions is preserved over an intermediate regime where correlation coefficient equals one and (time steps), in accordance with latora-baranger1999 (). Consequently, generalized Lyapunov exponent does not depend on the spreading occupation index latora-baranger-tsallis-rapisarda-2000 () of averaging cells, and it makes sense choosing the random pairs of initial conditions to compute the realizations all over the phase space, i.e. .

For the weakly chaotic case, stability islands dominate the phase space. Consequently, the characterization of the spreading occupation index of cells is much more intricate. In this case, a “global characterization” of the map has tentatively been made Baldovin-2002 (); Baldovin-2004 (). But Fig. 3 (, where LLE = ) demonstrates that a common transient time can not be defined, and this fact tends to preclude a neat linear dependence of with time. We now do not choose the random pairs of initial conditions all over the phase space, but in particular cells, i.e., in local regions. Nevertheless, we identify regions that exhibit , whose value of depend on the initial cell, as shown in Figure 4. In fact, a huge set of randomly chosen cells exhibit analogous results of , but we have also found cells where linear behavior is satisfied for (Fig. 5).

## Iii The rate of entropy production

With respect to the entropy production per unit time, we may conveniently use the -entropy (, henceforth) Tsallis1988 (); Tsallis2009 ()

(5) |

The -entropy production is estimated dividing the phase space in equal cells, and randomly choosing initial conditions inside one of the cells (typically ). We follow the spreading of points within the phase space, and calculate Eq. (5) from the set of occupancy probabilities (). Fluctuations are of course present, but they can be reduced when we choose initial cells within which the initial conditions are chosen, and average over the realizations. The proper value of the entropic parameter is the special value of which makes the averaged -entropy production per unit time to be finite, which is ultimately related with the extensivity of the entropy. The -entropy production per unit time is consequently obtained as

(6) |

and must be calculated taking into account that the partitions of phase space, and are such as to obtain robust results.

When the LLE is largely positive and the phase space is dominated by a chaotic sea, the Boltzmann-Gibbs entropy () is expected to be the appropriate one (i.e. ). Our numerical results show, as expected, that the proper value of the entropic index is for (see Figure 6). Both indexes have been obtained comparing the nonlinearity measure of the polynomial fitting curves over the intermediate regime (after a transient and before saturation) were the correlation coefficient is constant and equals . Another interesting result that emerged is the coincidence of the slopes of the sensitivity and entropy functions of time (see Fig. 2 and Fig. 6), i.e., the -entropy production per unit time satisfies . These results are numerically compatible with a -generalization of a Pesin-like identity for ensemble averages in strongly chaotic conservative maps, and reinforce those in ananos-tsallis-2004 ().

Entropic characterization of weak chaos regime in conservative maps is much more subtle than in dissipative maps (see Casati-Tsallis-Baldovin-2005 (); RuizTsallis ()). First, let us notice that one of the practical difficulties to identify the weak chaos regime are the huge computer times required, because of the slow convergence of the finite time Larger Lyapunov Exponent to its limit value. Furthermore, we know that the phase space of weakly chaotic conservative maps is dominated by stability islands and regions with a huge diversity of spreading occupation indexes, and this makes a hard task to reduce fluctuations by averaging the entropy. We have already pointed out that spreading occupation index of cells, and consequently their transient time scales, are essential in the -sensitivity characterization of the map. Fig. 7 reveals the incidence of this phenomena on the entropy production, which is related to the sensibility to initial conditions, and we conclude that a smart averaging criterion must be adopted.

Initially we have tested the adequacy of the Smaller Alignment Index (SALI) Skokos () to identify weak chaos in the standard map, so as to reduce the computer times required for LLE computation. We observe that the SALI temporal evolution is a computationally cheap tool to detect chaotic behavior of trajectories, as it tends to zero in both ordered and chaotic orbits, but with completely different time rates: SALI decreases abruptly in strongly chaotic orbits and reaches the limit of accuracy of the computer after about 200 iterations and, in contrast, it decreases with time as a power law in ordered orbits. We have verified, for the first time to the best of our knowledge, that the SALI decreases with time as a power law even in the case of weakly chaotic orbits. Fig. 8 shows, for the standard map with , that the time required to characterize weak chaos of stability island is generically shorter if we study the SALI temporal evolution instead of the LLE limit computation.

We then obtain the temporal production of the averaged entropy , selecting the cells between those that present analogous SALI temporal evolution: the smoother the log-log SALI evolution is, the smother its entropy production. This criterion provides averaged -entropies, , that do not present the typical oscillations observed in conservative maps RuizTsallis (). Fig. 9 exemplifies this fact.

Summarizing, for weakly chaotic dissipative systems, particularly for the logistic map, it has been observed that there is a specific value of (for the logistic map, latora-baranger-tsallis-rapisarda-2000 ()) which yields finite entropy production per unit time, and hence () is the proper one to thermodynamically characterize the system. Verification of this behavior for a conservative system (the standard map) is presently done. Our numerical results show in both strongly and weakly chaotic regime, as expected, , and this points towards the conjecture of the generalized Pesin equality tsallis-plastino-zheng-1997 (). In particular, and for (strong chaotic regime, as ), and for (weakly chaotic regime, as , ). With respect to the weakly chaotic regime, a local characterization of generalized Pesin-like remains to be made, as and have both shown to be local properties.

## Iv The rate of relaxation of the entropy

The methods described in Sections II and III take an ensemble of initial conditions inside one of cells of the phase space, and thus , . In Ref. moura-tirnakli-lyra-2000 (), the authors have considered an ensemble of initial conditions spread over the entire phase space of a dissipative system (particularly, the logistic map), therefore, with the maximum possible value for the entropy . The time evolution brings the system to its attractor, and the number of occupied cells falls according to an exponential relation, in the case of a strongly chaotic system, or according to a -exponential, in the case of a weakly chaotic system. Later it was found a relation between and by analyzing the rate in which the entropy goes to its final value, beginning with , instead of ananos-borges-tsallis-oliveira-2002 ().

The method adopted in moura-tirnakli-lyra-2000 () cannot be identically applied to a conservative map, since there is no attractor, but that one introduced in ananos-borges-tsallis-oliveira-2002 () can be adapted to area-preserving models. We consider the variable and evaluate its average over an ensemble of initial conditions for . The average must be computed with a different criterion in strong and weak chaos regimes, according to the respective stability distributions on the phase space.

The strongly chaotic map is dominated all over the partitioned phase space by a chaotic sea, and we can analytically calculate . In that case, the average has been estimated over the most quickly spreading cells to minimize the time computation. On the contrary, the stability islands dominate the phase space in the weakly chaotic regime, and their respective value of is unknown, as it numerically depends on the cells where the initial conditions are randomly chosen. We must find a proper criterion for the numerical estimation of and , and we privilege those cells that cover a region were SALI behaves in a smoothly manner, such as to produce entropy in a smoothly manner too.

It is expected that , where is a relaxation time. Figures 10 and 11 show the time dependence of for the strongly chaotic case (), and the weakly chaotic case (), respectively. The former case displays an exponential relaxation, i.e., , consistently with BG framework. The weakly chaotic case presents a -exponential relaxation regime, with for .

## V Stationary distributions

Probability density distributions are among the most relevant items regarding the statistical description of dynamical systems. Gaussian distributions are typical for ergodic and mixing systems, for which LLE is positive tirnakli-beck-tsallis-2007 (); tirnakli-borges-2016 (), and -Gaussian distributions are typical for weakly chaotic systems for which LLE tirnakli-beck-tsallis-2007 (); tirnakli-beck-tsallis-2009 (); tirnakli-borges-2016 () and or chaotic layers are thin and LLE Ruiz-Bountis-tsallis-2012 ().

We define the -Gaussian distribution as

(7) |

where is the -mean value, is the -variance, is the normalization factor and is a parameter which characterizes the width of the distribution prato-tsallis-1999 ():

(8) |

(9) |

and where recovers the Gaussian distribution

Both Gaussian and -Gaussian distributions have recently been observed for the standard map tirnakli-borges-2016 (). We have analyzed the distribution of sums of iterates of the map, for representative values, as defined by

(10) |

where the average is both a time average over iterations and an ensemble average over initial conditions:

(11) |

The probability distributions of , namely , that emerge from the standard map with weakly chaotic phase space are -Gaussians, while Gaussian distributions are associated with the strongly chaotic regime. The scenario for intermediate values of is also quite rich. For these cases, the phase space displays regions with points that have positive LLE, — the chaotic sea —, and regions with ,— the stability islands. The probability distributions of for initial conditions taken inside the strongly chaotic regions are Gaussians, and those distributions for initial conditions that belong to the weakly chaotic regions are -Gaussians. One remarkable feature is that in all cases, regardless the value of .

If the initial conditions are spread over the entire phase space, thus the initial conditions are from both the chaotic sea and the stability islands, the distribution that emerges is a linear combination of a Gaussian and a -Gaussian. Consequently, the probability distribution of the standard map, for any arbitrary value of , can be modeled as

(12) |

where and are the probability densities of for the initial conditions taken inside the weak and strong chaos regions, respectively.

In Ref. tirnakli-borges-2016 () it was taken the linear combination of the probability distributions normalized to the maximum value , and the physical meaning of the interpolating parameter was not clear. Now we have taken the linear combination of the usual probability densities in Eq. (12). We find out the ratio of the areas for the strongly and weakly chaotic regions. In order to achieve this, we randomly generate initial conditions all over the phase space and count how many of those having non-zero Lyapunov and those having nearly zero Lyapunov. Dividing these numbers by the total number of initial conditions, we mimic the ratio of the areas. We classify an initial condition belonging to a weakly chaotic region if its Lyapunov exponent is smaller than a given threshold (time steps). Correspondingly, an initial condition is considered to belong to a strongly chaotic region whenever (time steps). These limiting values and are chosen according to what is explained in tirnakli-borges-2016 (): the phase-space ratio, i.e., [number of points with ] / [number of points with ] remains almost constant for . We call , and . There remains a small amount of initial conditions with intermediate Lyapunov exponent values () that are not taken into account in Eq. (12), and that are possibly negligible in the macroscopic limit. We have thus considered the approximated form of Eq. (12) as

(13) |

See the results of some representative values in the Table.

, | ||||

0.2 | 0.015 | 0 | 1 | |

0.6 | 0.015 | 0.0380 | 0.9416 | |

1 | 0.015 | 0.4826 | 0.5078 | |

2 | 0.010 | 0.7622 | 0.2367 | |

3 | 0.010 | 0.8825 | 0.1170 | |

10 | 1 | 0 |

Then, we compute the and values using the pdf of the system coming from the initial conditions with positive and with zero Lyapunov values, respectively. After determining these values now we have no fitting parameter and values are entirely determined by the intrinsic dynamics of the map. For some representative values of , the results are given in Fig. 12.

## Vi Concluding Remarks

The -triplet (, , ) is an important feature in the statistical description of dynamical systems that might be nonergodic and nonmixing, as a consequence of zero LLE. These indices define the rate of entropy production, the rate of relaxation, and the distribution of the stationary states. Positive LLE leads to ergodic, mixing systems, and strong chaos, and these indices collapse into , within the BG framework. Zero LLE may lead to breakdown of ergodicity, nonmixing dynamics, and weak chaos. The statistical description in many such circumstances appears to be associated with the -entropy, where the usual exponential and Gaussian functions turn into more general forms, namely the -exponential and the -Gaussian ones. This scenario has been previously observed and evaluated for dissipative low-dimensional systems.

We have now considered a paradigmatic low-dimensional conservative system, the standard map. We have evaluated the -triplet, and have verified that the sensitivity to initial conditions index satisfies , in both strongly and weakly chaotic regimes. Our results corroborate the expected situation, i.e., that if the phase space is completely strongly chaotic (LLE), and . These results are numerically compatible with a -generalization of a Pesin-like identity for ensemble averages, extending those in ananos-tsallis-2004 () to the case of conservative maps. If the entire phase space is weakly chaotic regime (), which happens for low values of , the standard map is characterized by , , and . We have shown that is a property that can be locally characterized by the finite temporal evolution of LLE or, alternatively, by the SALI, which has demonstrated to be a computationally cheaper tool to identify different chaotic regimes. The phase space for intermediate values of features regions of positive and regions of zero LLE. Ensemble averages of initial conditions taken within the strongly chaotic regions behave according to . In contrast, ensemble averages of initial conditions taken within the weakly chaotic regions behave according to the -triplet and, at least in what concerns the stationary distribution, with the same values of . This work extends what has been done in Ref. tirnakli-borges-2016 (). We have numerically found the proper superposition of distributions for the mixed case (coexistence of strongly and weakly chaotic regimes) without any additional fitting parameter.

The -triplet structure plays a central role on the statistical description of dynamical systems, and additional works addressing other conservative and dissipative maps shall bring further tests for the robustness of this framework, and they are very welcome.

## Acknowledgments

This work has been partially supported by CNPq and Faperj (Brazilian Agencies), and by TUBITAK (Turkish Agency) under the Research Project number 115F492. One of us (CT) also acknowledges partial financial support by the John Templeton Foundation (USA).

### References

- C. Tsallis, Introduction to Nonextensive Statistical Mechanics – Approaching a Complex World (Springer, New York, 2009).
- C. Beck and F. Schlögl, Thermodynamics of chaotic systems – An introduction (Cambridge Nonlinear Science, 1997).
- A. Rapisarda and V. Latora, in Nonextensive Entropy - Interdisciplinary Applications, Edited by Gell-Mann, C. Tsallis (Oxford University Press, New York, 2004).
- A. Esquivel and A. Lazarian, ApJ 710, 125 (2010).
- CMS Collaboration, Phys. Rev. Lett. 105, 022002 (2010).
- M. Shao, L. Yi, Z. Tang, H. Chen, C. Li and Z. Xu, J. Phys. G Nucl Partic 37, 085104 (2010).
- O. Sotolongo-Grau, D. RodrÃguez-PÃ©rez, J. C. Antoranz and O. Sotolongo-Costa, Phys. Rev. Lett. 105, 158105 (2010).
- Mohanalin, Beenamol, P. K. Kalra and N. Kumar, Comput. Math. Appl. 60, 2426 (2010).
- ALICE Collaboration, Phys. Lett. B 693, 53 (2010).
- PHENIX Collaboration, Phys. Rev. D 83, 052004 (2011).
- PHENIX Collaboration, Physical Review C 83, 024909 (2011).
- ALICE Collaboration, Eur. Phys. J. C C 71, 1655 (2011).
- ATLAS Collaboration, New J. Phys. 13, 053033 (2011).
- C. Y. Wong and G. Wilk, Phys. Rev. D 87, 114007 (2013)
- L. Marques, E. Andrade-II and A. Deppman, Phys. Rev. D 87, 114022 (2013)
- M. I. Bogachev, A. R. Kayumov and A. Bunde, PLoS ONE 9, e112534 (2014).
- G. Combe, V. Richefeu, M. Stasiak and A. P. F. Atman, Phys. Rev. Lett. 115, 238301 (2015).
- F. Baldovin and A. Robledo, Phys. Rev. E 66, 045104(R) (2002).
- G. L. Ferri, M. F. Reynoso Savio and A. Plastino, Physica A 389, 1829 (2010).
- L. F. Burlaga and N. F. Ness, Astrophys. J. 765, 35 (2013).
- B. V. Chirikov, Phys. Rep. 52, 264 (1979).
- G. M. Zaslavsky, Hamiltonian Chaos and Fractional Dynamics (Oxford University Press, New York, 2005).
- G. Benettin, L. Galgani, A. Giorgilli and J.-M. Strelcyn, Meccanica 15, 9 (1980).
- U. Tirnakli and E. P. Borges, Sci. Rep. 6, 23644 (2016).
- M. Baranger, V. Latora, and A. Rapisarda, Chaos, Solitons & Fractals 13, 471 (2002)
- C. Tsallis, A. R. Plastino and W.-M. Zheng, Chaos, Solitons & Fractals 8, 885 (1997)
- C. Tsallis, Quimica Nova 17, 468 (1994).
- G. F. J. Añaños and C. Tsallis, Phys. Rev. Lett. 93, 020601 (2004).
- V. Latora and M. Baranger, Phys. Rev. Lett. 82, 520 (1999).
- V. Latora, M. Baranger, A. Rapisarda and C. Tsallis, Phys. Lett. A 273, 97 (2000).
- F. Baldovin, Physica A 305, 124 (2002).
- F. Baldovin, in Nonextensive Entropy - Interdisciplinary Applications, Edited by Gell-Mann, C. Tsallis (Oxford University Press, New York, 2004).
- C. Tsallis, J. Stat. Phys. 52, 479 (1988).
- G. Casati, C. Tsallis and F. Baldovin, Europhys. Lett. 72, 355 (2005).
- G. Ruiz and C. Tsallis, Physica A 386, 720 (2007).
- Ch. Skokos, J. Phys. A: Math. Gen. 34, 10029 (2001).
- F. A. B. F. de Moura, U. Tirnakli and M. L. Lyra, Phys. Rev. E 62, 6361 (2000).
- E. P. Borges, C. Tsallis, G. F. J. Añaños, and P. M. C. de Oliveira, Phys. Rev. Lett. 89, 254103 (2002).
- U. Tirnakli, C. Beck and C. Tsallis, Phys. Rev. E 75, 040106(R) (2007).
- U. Tirnakli, C. Tsallis and C. Beck, Phys. Rev. E 79, 056209 (2009).
- G. Ruiz, T. Bountis and C. Tsallis, Int. J. Bif. and Chaos 22, 1250208 (2012).
- D. Prato, C. Tsallis, Phys. Rev. E 60, 2398 (1999).