## The Scaling of Chaos vs Periodicity: How Certain is it that an Attractor is Chaotic?

Madhura Joglekar, 1 Edward Ott, 1 and James A. Yorke 1

1

University of Maryland, College Park, Maryland 20742, USA

(Dated: 5th August 2014)

The character of the time-asymptotic evolution of physical systems can have complex, singular behavior with variation of a system parameter, particularly when chaos is involved. A perturbation of the parameter by a small amount glyph[epsilon1] can convert an attractor from chaotic to non-chaotic or viceversa. We call a parameter value where this can happen glyph[epsilon1] -uncertain. The probability that a random choice of the parameter is glyph[epsilon1] -uncertain commonly scales like a power law in glyph[epsilon1] . Surprisingly, two seemingly similar ways of defining this scaling, both of physical interest, yield different numerical values for the scaling exponent. We show why this happens and present a quantitative analysis of this phenomenon.

While low-dimensional chaotic attractors are common and fundamental in a vast range of physical phenomenon, it is predominantly the case that chaotic motions in such systems are 'structurally unstable'. This typically has the consequence that an arbitrarily small change of a system parameter can always be found that results in periodic behavior [1]. Physical models displaying chaotic attractors that are structurally unstable arise very often, e.g., in studies of plasma dynamics [2], Josephson junctions [3], chemical reactions [4] and many others. We also note that even the simple example of the one-dimensional quadratic map, observed uncertainty, applicable when all chaotic attractors of any size are considered [6], and this observation is at the heart of the resolution of the above-mentioned puzzle.

$$x _ { n + 1 } = C - x _ { n } ^ { \ 2 }, \quad \ \ \ ( 1 ) \quad \epsilon _ { - 1 }$$

displays this phenomenon, and, in this paper we will study this example as a convenient paradigm for such situations in general.

One reason for concern with this type of behavior is that physical systems have uncertainties in the values of their parameters, and one might therefore ask how confident one can be about the prediction of chaotic behavior from a model calculation (even when the model and its parameter dependence are precisely known and there is no noise). Despite the fundamental importance of this question, very little study has been done to quantitatively address it [5-7]. It is the purpose of this paper to readdress this general issue, and, in particular, to resolve a long-standing puzzle. This puzzle has to do with the scaling characterization of the fractal-like chaotic/periodic interweaving structure of parameter dependence associated with structural instability. In particular, studies on the quadratic map Eq.(1) have addressed scaling in two slightly different ways and obtained significantly different estimates of the scaling exponent [5-7]. The reason for this surprising discrepancy has remained unresolved. In one of the two ways of addressing scaling, attention was restricted to what might seem to be the most obvious source of uncertainty, namely, when a large (to be defined subsequently) chaotic attractor suddenly turns into a periodic attractor, as a parameter is varied. However, we find such transitions far too rare [5, 7] to account for the

Consider a dynamical system depending on a parameter C , and an attractor A C ( ) that, as C varies continuously in some range, can be uniquely associated with C . We say that a particular value of C is glyph[epsilon1] -uncertain with respect to chaos if A C ( ) is chaotic while either A C ( + ) glyph[epsilon1] or A C ( -glyph[epsilon1] ) or both are not chaotic. For example, for the case of the quadratic map, to which we henceforth restrict our considerations, it has been found that a random choice of C with uniform probability density yields glyph[epsilon1] -uncertainty with respect to chaos with a probability ¯ ( F glyph[epsilon1] ) that scales like a power of glyph[epsilon1] for small glyph[epsilon1] [6], ¯ ( F glyph[epsilon1] ) ∼ glyph[epsilon1] β . As a result of the fine-scaled interweaving of C values for which A C ( ) is chaotic and intervals of C values for which A C ( ) is periodic, the power law exponent β turns out to be less than 1, and the set of C values with A C ( ) chaotic has been called a 'fat-fractal' [5, 6] (more precisely, it is a Cantor set with positive Lebesgue measure [8]). We have repeated the numerical determination of ¯ ( F glyph[epsilon1] ) (Fig. 1) and find that for small glyph[epsilon1] ,

$$\bar { F } ( \epsilon ) \sim \epsilon ^ { \beta } \text{ with } \beta = 0. 3 9 2 \pm 0. 0 3 7. \quad \text{ (2)}$$

In performing this calculation we estimate ¯ ( F glyph[epsilon1] ) by first randomly choosing many C values with uniform probability in the range where the quadratic map has a unique bounded attractor, -1 / 4 ≤ C ≤ 2 . For each such C value we then compute the Lyapunov exponents for C , C + glyph[epsilon1] and C -glyph[epsilon1] , judging the corresponding attractors to be chaotic or not depending on whether the computed Lyapunov exponent is positive. We then estimate ¯ ( F glyph[epsilon1] ) as the fraction of those randomly chosen C -values that are computed to be glyph[epsilon1] -uncertain with respect to chaos. Reference [6] obtains a slightly larger β value of β ≈ 0 41 using . an glyph[epsilon1] -range with larger glyph[epsilon1] -values. We agree with their result in the range they tested, but, by pushing to small glyph[epsilon1] , obtain the result in Eq. (2) (See Fig. 1.)

We comment that Eq. (2) can be interpreted as implying a type of 'probabilistic stability' for chaos. That is,

FIG. 1: ¯ ( F glyph[epsilon1] ) vs glyph[epsilon1] .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0311v1-figure-9fa7f968-b81c-4399-b165-bbecb2f90a65.jpg)

while chaos occuring at some parameter value C may be structurally unstable in the sense that a chaos-destroying perturbation C → C + δC , can always be found with | δC | ≤ ˜ for any given ˜; chaos is still stable in the sense glyph[epsilon1] glyph[epsilon1] that, for ˜ small, such a chaos-destroying glyph[epsilon1] δC may have to be very carefully chosen, and the probability that a random choice of δC in | δC | ≤ ˜ destroys chaos approaches glyph[epsilon1] zero as ˜ is made smaller and smaller [9] (Eq. glyph[epsilon1] (2)). (Incidently, we note that quasiperiodicity appears to have this same type of structural-instability / probabilistic stability [6].)

We now recall the concept of a periodic window. A p -periodic window is an interval of the parameter, C ∗ ( p ) ≤ C ≤ C x ( p ) , such that at the beginning of the window, as C increases through C ∗ ( p ) , there is a bifurcation from a chaotic attractor to a periodic orbit attractor of period p , followed by a period-doubling cascade to chaos, followed by a sequence of band-mergings in each of which 2 m p separate pieces ( x -intervals) of the chaotic attractor pairwise merge into a 2 m -1 p -piece chaotic attractor, eventually forming a 'small' p -piece chaotic attractor, which subsequently terminates ('explodes') at the end of the window ( C = C x ( p ) ) through a crisis transition [10] to a larger chaotic attractor that is similar in size to the larger chaotic attractor just before the beginning of the window at C = C ∗ ( p ) . Thus, as is evident from viewing a bifurcation diagram for the quadratic map, 'small' chaotic attractors occur within windows only, and we call a chaotic attractor that is not contained in any window a 'large chaotic attractor'.

Now, instead of considering glyph[epsilon1] -uncertainty with respect to chaos, we consider glyph[epsilon1] -uncertainty with respect to the occurence of large chaotic attractors. That is, we consider C to be glyph[epsilon1] -uncertain if A C ( ) is a large chaotic attractor while either A C ( + ) or glyph[epsilon1] A C ( -glyph[epsilon1] ) or both are not large chaotic attractors. Reference [7] gives a detailed consideration and analysis of the scaling of glyph[epsilon1] uncertainty with respect to the occurence of large chaotic attractors for Eq.

(1) obtaining for the probability F 0 ( glyph[epsilon1] ) of glyph[epsilon1] -uncertainty

$$F _ { 0 } ( \epsilon ) \sim \epsilon ^ { \alpha } \text{ with } \alpha = 0. 5 1 \pm 0. 0 3. \quad \text{ (3)}$$

Thus, there are evidently two distinct scaling exponents, the β and α of Eqs. (2) and (3). We conjecture that the values of these exponents are universal for onedimensional maps with a quadratic maximum. For example, they would apply to physical situations like those in the plasma example in Ref. [2] and the chemical example in Ref. [4], where strong phase-space attraction leads to dynamics closely approximated by a one-dimensional map.

It is important to note that the chaos within windows has its own windows, which, in turn, have their own windows, and so on. Thus, it is useful to distinguish the 'order' of a window: We say a window is of order one (which we also call a 'primary' window), if it is not contained within another window, and we say a window is of order r &gt; 1, if it is contained within a window of order ( r -1), but it is not contained within a window of order ( r + 1). Thus, if a C value is not contained within a primary window, then it is not contained in any window (as is the case for a large chaotic attractor). For later reference, we define N r (∆) to be the number of windows of order r whose widths are greater than ∆ and we write this quantity for r = 1 (primary windows) as

$$N _ { 1 } ( \Delta ) = \sum _ { i = 1 } ^ { \infty } U ( \Delta _ { i } - \Delta ),$$

where ∆ 1 ≥ ∆ 2 ≥ ∆ 3 ≥ . . . denote the widths of the primary windows, and U denotes the unit step function ( U z ( ) = 1 for z &gt; 0, U z ( ) = 0 for z &lt; 0). If we consider all the windows of any order, then the number of these windows whose widths exceed ∆ is

$$\bar { N } ( \Delta ) = \sum _ { r = 1 } ^ { \infty } N _ { r } ( \Delta ).$$

Since glyph[epsilon1] 0 39 . &gt;&gt; glyph[epsilon1] 0 51 . for small glyph[epsilon1] , the results (2) and (3) imply that for small glyph[epsilon1] most of the parameter values that are glyph[epsilon1] -uncertain with respect to chaos lie in windows, and, in fact, as we will demonstrate elsewhere, for a randomly chosen parameter value ˜ C that is glyph[epsilon1] -uncertain with respect to chaos, the expectation value of the order of the lowest order window containing ˜ C approaches ∞ as glyph[epsilon1] → 0.

Preliminary to our analysis of the relationship between the exponents α and β , we need the following two results:

- (i) Self-similarity of windows: Reference [11] shows that the bifurcation structure and dynamics in windows of various orders and periods are self-similar. That is, considering x near 0, by use of uniform linear stretchings (magnification) in x and C , the p -times iterated map with the parameter ranging through the interval corresponding to a periodp

window, C ∗ ( p ) ≤ C ≤ C x ( p ) , very closely quantitatively replicates the behavior of the map in its full range, -1 / 4 ≤ C ≤ 2 . Furthermore, this selfsimilarity approximation is already extremely good even for the period-3 primary window, and becomes better and better for most r -order windows as r is increased. This is illustrated in Fig. 2, where we see that the bifurcation diagram in the full range -1 / 4 ≤ C ≤ 2 in Fig. 2a is virtually identical to a properly magnified, inverted ( x →-x ) version of the bifurcation diagram for the period-3 window, C ∗ (3) ≤ C ≤ C x (3) , blown up in the region near x = 0 (Fig. 2b).

- (ii) glyph[epsilon1] -uncertainty / window-width scaling equivalence: It is shown in the supplementary material that, for small window widths ∆, the scalings of ¯ (∆) and N N 1 (∆) are related to the small glyph[epsilon1] scalings of ¯ ( F glyph[epsilon1] ) and F 0 ( glyph[epsilon1] ) :

$$\bar { F } ( \epsilon ) \sim \epsilon \bar { N } ( \epsilon ) \text{ and } F _ { 0 } ( \epsilon ) \sim \epsilon N _ { 1 } ( \epsilon ). \quad \ \ ( 6 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Thus, from Eqs. (2) and (3), Eq. (6) yields

$$\bar { N } ( \Delta ) \sim \frac { 1 } { \Delta ^ { 1 - \beta } } \text{ and } N _ { 1 } ( \Delta ) \sim \frac { 1 } { \Delta ^ { 1 - \alpha } }. \quad \text{ (7)}$$

FIG. 2: (a) shows a bifurcation diagram for the quadratic map for -0 25 . ≤ C ≤ 2. (b) shows a blow-up of the bifurcation diagram in the period-3 window in the region near x = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0311v1-figure-60b7bbfa-77f7-4e00-8939-9971c82c0206.jpg)

By using Eq. (7), we reduce our analysis of the glyph[epsilon1] -uncertainty exponents α and β to an analysis of the scaling of the distribution of window-widths. Now, using the self-similarity of windows, we can express N r +1 (∆) in terms of N r (∆),

$$N _ { r + 1 } ( \Delta ) = \int _ { \Delta } ^ { 9 / 4 } ( \frac { - d N _ { r } ( \Delta ^ { \prime } ) } { d \Delta ^ { \prime } } ) \ N _ { 1 } ( \frac { 9 \Delta / 4 } { \Delta ^ { \prime } } ) \ d \Delta ^ { \prime }, \quad ( 8 )$$

where 9 4 / = 2 - -( 1 / 4) is the width of the C -range for the quadratic map. In Eq. (8), we have used the window self-similarity result to write the number of r +1 order windows of width ∆ that are contained in windows of width ∆ ′ as N 1 ( 9∆ 4 / ∆ ′ ). Letting u = ln(9 / (4∆)) and M u r ( ) = N r (∆) , Eq. (8) becomes,

$$M _ { r + 1 } ( u ) = \int _ { 0 } ^ { u } M _ { 1 } ( u - u ^ { \prime } ) \ \frac { d M _ { r } ( u ^ { \prime } ) } { d u ^ { \prime } } \ d u ^ { \prime }. \quad ( 9 )$$

Equation (9) is a convolution. It is therefore convenient to introduce the Laplace transform, ˆ M s r ( ) = ∫ ∞ 0 M u e r ( ) -su du , in terms of which Eq. (5) becomes ˆ M r +1 ( s ) = sM ˆ 1 ( s M ) ˆ r ( s ) (where we have made use of M r (0) = N r (9 / 4) = 0). Iterating this result, we obtain ˆ M s r ( ) in terms of ˆ M s 1 ( ),

$$\hat { M } _ { r } ( s ) = \frac { ( s \hat { M } _ { 1 } ( s ) ) ^ { r } } { s }. \quad \quad \quad ( 1 0 )$$

Introducing the Laplace transform ˆ ( ¯ M s ) of ¯ ( M u ) = ¯ (∆), we have from Eq. N (5) that ˆ ( ¯ M s ) = ∑ ∞ r =1 ˆ M s r ( ), which using Eq. (10) gives, upon summing the resultant geometric series,

$$\hat { \hat { M } } ( s ) = \frac { \hat { M } _ { 1 } ( s ) } { 1 - s \hat { M } _ { 1 } ( s ) }. \quad \quad ( 1 1 )$$

This expression is singular at values of s for which sM ˆ 1 ( s ) = 1 . From Eq. (4), M u 1 ( ) = ∑ ∞ i =1 U u ( -u i ) , where u i = ln(9 / (4∆ )) i . Thus, ˆ M s 1 ( ) = s -1 ∑ ∞ i =1 e -su i = ∑ ∞ i =1 (4∆ i / 9) s and we have that singularities of ˆ ( ¯ M s ) occur at s values satisfying 1 = ∑ ∞ i =1 δ i s , where δ i = (4∆ i / 9) are the normalized widths of the primary windows.

Now consider the inverse Laplace transform [12] of ˆ ( ¯ M s ) for large u (i.e., small ∆); we see that ¯ ( M u ) ∼ e γu where γ is the solution for s of 1 = ∑ i δ i s with the largest real part. It can be shown that this solution is real. Thus we obtain for ¯ (∆) N at small ∆, ¯ (∆) N ∼ ∆ -γ , which when compared with the first part of our result Eq. (7) shows that γ = 1 -β. We conclude that we can obtain the exponent β as the real positive root of

$$\sum _ { i = 1 } ^ { \infty } \delta _ { i } ^ { \gamma } = 1, \ \beta = 1 - \gamma.$$

We now use Eq. (12) to investigate the relationship between the exponents α and β . Writing Eq. (12) as ∑ I -1 i =1 δ i γ + ∑ ∞ i = I δ i γ = 1, for sufficiently large I , we can approximate the second summation by

$$\sum _ { i = I } ^ { \infty } \delta _ { i } ^ { \gamma } \approx \int _ { 0 } ^ { \delta _ { I } } ( - d N _ { 1 } ( \delta ) / d \delta ) \delta ^ { \gamma } d \delta. \quad \ ( 1 3 )$$

Since N 1 ( δ ) ∼ δ -(1 -α ) , the integrand is proportional to δ γ + α -2 , and the integral diverges unless γ + α &gt; 1 . Thus, the summation in (12) is infinity unless γ &gt; (1 -α ), and,

since each term in the sum decreases monotonically with increasing γ , we conclude that as γ increases past (1 -α ), the sum in (12) decreases monotonically.

For γ = 1, the sum is the normalized total length of all windows, which, by definition, is less than 1. We conclude that Eq. (12) has a single root for γ and that this root satisfies γ &lt; (1 -α . ) Thus, for γ = 1 -β , we must have that α &gt; β, in agreement with the numerical results α ≈ 0 51 . , β ≈ 0 39 . . Taking N 1 ( δ I ) = I ( δ I /δ ) 1 -α and performing the integration in (13), we obtain

$$\sum _ { i = 1 } ^ { I - 1 } \delta _ { i } ^ { 1 - \beta } + \frac { 1 - \alpha } { \alpha - \beta } I \delta _ { I } ^ { 1 - \beta } = 1. \quad \quad ( 1 4 )$$

For example, assuming availability of an estimate of α, one can use (14) to estimate β given numerical determinations of ( δ 1 , δ 2 , . . . , δ I ). (Note that the second term in (14) becomes smaller and smaller as I is increased and can be omitted for very large I .)

Figure 3 shows a plot of the estimated value of β obtained from (14) as a function of I for α = 0 51. . The result from Fig. 3 is β = 0 39 in good agreement with . the estimate from Fig. 1.

FIG. 3: Estimated value of β vs I . As I is increased, the estimated value of β converges to ≈ 0 39 . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0311v1-figure-32091578-2cb4-43a5-b415-052c4c49d7b3.jpg)

In conclusion, we have derived an analytic estimate Eq. (12) for β which yields good agreement with the numerical result in Eq. (2) and shows why α &gt; β. More generally, letting Q 0 and ¯ Q respectively denote the set of C values yielding large chaotic attractors and the set of C values yielding chaotic attractors of any size, one can view our work as using the self-similarity of windows to establish a quantitative link between the structure of these two sets. In particular, Eq. (12) relates the primary window widths { δ 1 , δ 2 , δ 3 , . . . } (a characterization of Q 0 ) to the exponent β (a characterization of ¯ ) Q. Although our considerations have focused on the quadratic map, we believe that the numerical results for the exponents α and β are universal for one-dimensional maps with a single quadratic maximum, and thus apply for situations as in [2], [4] where there is a strong phase-space contraction.

The work of E. O. was supported by ARO grant W911NF1210101.

- [1] That is, given a system parameter value C for which there is a chaotic attractor and a small value of ˜ glyph[epsilon1] &gt; 0, one can find a parameter perturbation C → C + δC with | δC | &lt; glyph[epsilon1] ˜ such that the attractor is no longer chaotic, and this can always be done no matter how small ˜ is. glyph[epsilon1]
- [2] J. M. Wersinger, J. M. Finn, and E. Ott, Phys. Fluids 23 , 1142 (1980).
- [3] B. A. Huberman, J. P. Crutchfield, and N. H. Packard, Appl. Phys. Lett. 37 , 750 (1980).
- [4] F. Argoul, A. Arneodo, P. Richetti, J. C. Roux, and H. L. Swinney, Acc. Chem. Res. 20 , 436 (1987).
- [5] J. D. Farmer, Phys. Rev. Lett. 55 , 351 (1985).
- [6] C. Grebogi, S. W. McDonald, E. Ott, and J. A. Yorke, Phys. Lett. A 110 , 1 (1985).
- [7] B. R. Hunt and E. Ott, J. Phys. A 30 , 7067 (1997).
- [8] M. V. Jacobson, Comm. Math. Phys. 81 , 39 (1981), has shown that C values for which A C ( ) is chaotic form a positive Lebesgue measure set, while J. Graczyk and G. Swiatek [Ann. of Math. 146, 1 (1997)] demonstrate that intervals of C values with periodic attractors are dense in the region, C F &lt; C &lt; 2, where chaos can occur (here C F denotes the 'Feigenbaum point' at which period doublings of periods 2 m accumulate).
- [9] We comment that this is a consequence of Lebesgue's density theorem which states that, for a Lebesgue measurable set Q (in this case the set of C values that yield chaotic attractors), for almost every point x in Q , the limit of glyph[epsilon1] → 0 of µ Q ( ∩ [ x -glyph[epsilon1], x + ]) (2 ) is one, where glyph[epsilon1] / glyph[epsilon1] µ denotes the Lebesgue measure.
- [10] C. Grebogi, E. Ott, and J. A. Yorke, Physica D 7(1) , 181 (1983).
- [11] J. A. Yorke, C. Grebogi, E. Ott, and L. TedeschiniLalli, Phys. Rev. Lett. 54 , 1095 (1985), see Sec. B of Supplementary Material.
- [12] In taking the Laplace transform ∫ ∞ 0 M u ( ) exp( -su du ) and the inverse Laplace transform (2 πi ) -1 ∫ p ˆ ( M s ) exp( su ds ) , where the contour P runs from - ∞ i + σ to + i ∞ + σ , we take Re s ( ) &gt; γ and σ &gt; γ , respectively.
- [13] There is some small error attached to our use of the approximate self-similarity of windows in our derivation of Eq. (12). As shown in Table 1 of Ref. [11] (reproduced in Sec. B of the supplementary material), however, the self-similarity assumption is quite good, and, as demonstrated in Sec. C of the supplementary material, we estimate that this small deviation from self-similarity is not significant compared to the error bounds of the numerical determination of β indicated in Eq. (2).
Madhura Joglekar Edward Ott James A. Yorke
The Scaling of Chaos vs Periodicity: How Certain is it that an Attractor is Chaotic?
The character of the time-asymptotic evolution of physical systems can have
complex, singular behavior with variation of a system parameter, particularly
when chaos is involved. A perturbation of the parameter by a small amount
$\epsilon$ can convert an attractor from chaotic to non-chaotic or vice-versa.
We call a parameter value where this can happen $\epsilon$-uncertain. The
probability that a random choice of the parameter is $\epsilon$-uncertain
commonly scales like a power law in $\epsilon$. Surprisingly, two seemingly
similar ways of defining this scaling, both of physical interest, yield
different numerical values for the scaling exponent. We show why this happens
and present a quantitative analysis of this phenomenon.