Scale-free networks emerging from multifractal time series

Marcello A. Budroni,1 Andrea Baronchelli,2 and Romualdo Pastor-Satorras3 1Nonlinear Physical Chemistry Unit, Facult´edes Sciences, Universit´elibre de Bruxelles (ULB), CP 231 - Campus Plaine, 1050 Brussels, Belgium. 2Department of Mathematics - City, University of London - Northampton Square, London EC1V 0HB, UK 3Departament de F´ısica, Universitat Polit`ecnica de Catalunya, Campus Nord B4, 08034 Barcelona, Spain (Dated: November 5, 2018) Methods connecting dynamical systems and graph theory have attracted increasing interest in the past few years, with applications ranging from a detailed comparison of different kinds of dynamics to the characterisation of empirical data. Here we investigate the effects of the (multi)fractal properties of a time signal, common in sequences arising from chaotic or strange attractors, on the topology of a suitably projected network. Relying on the box counting formalism, we map boxes into the nodes of a network and establish analytic expressions connecting the natural measure of a box with its degree in the graph representation. We single out the conditions yielding to the emergence of a scale-free topology, and validate our findings with extensive numerical simulations.

PACS numbers: 89.75.Hc, 05.45.Tp, 02.50.Ey, 05.45.Ac

I. INTRODUCTION the nodes, which are connected if the signal transits from one state to another. This transformation procedure pre- Network science [1] has emerged in the last decades as serves the temporal information of the dynamics in the a transverse interpretative framework for understanding network, is stable to noise affecting real time series and the structure and function of a wide range of complex can account for the fundamental structure of the related systems, as well as the dynamic phenomena that takes attractors [21]. place on them, ranging from financial crises and traffic Using this kind of network projections, several aspects congestion to epidemic and social influence spreading [2– of dynamical systems have been cast in terms of topolog- 4]. ical network properties. Here we will focus in particular In the realm of dynamical systems [5], network tech- on the effects that the fractal and multifractal proper- niques have been applied to the analysis of nonlinear time ties of a time signal have on the topology of representa- series, with a particular focus on characterizing chaotic tive classes of projected network. Many time signal, par- dynamics [6]. The main idea of this methodology is to ticularly those arising from a chaotic or strange attrac- project the information of a time series into the topology tor [6], have a fractal structure, characterized by a sta- of a network. The key element of this approach resides tistically self-similar pattern in phase space. Moreover, in the identification of nodes and links in the network they can also show multifractal properties, described by from the time series information. Several alternatives a strongly heterogeneous probability of visiting different have been proposed in this context [7]. Thus, Zhang and neighborhoods in the phase space [22–24]. Using as a sim- Small [8] consider cycles of a pseudo-periodic time series ple example a transition network representation [19, 25], as the nodes of a network, which are connected by links framed within a box counting formalism [23] in which depending on the similarity between cycles. Lacasa at each box corresponds to a vertex, we find analytic ex- al. [9, 10], on the other hand, build on the concept of pressions tying the visitation probability of a box with visibility graphs, where nodes correspond to series data its associated degree. From these relations we obtain the points and two nodes are connected if a straight line can conditions under which a projected network grows with a be established between them without intersection with scale-free topology [26], characterized by a degree distri- any intermediate data height. bution of the form P (k) ∼ k−γ . Our results highlight the Another general approach is based on encoding topo- correspondence between an attractor’s structure and the logical information from the reconstructed phase space of topology of the associated projected network, and relate a time series into a proximity network. In these networks, the possible origin of a heterogeneous scale-free topology nodes represent segments of time series or vectors in the with the heterogeneous and hierarchical visitation prob- related reconstructed phase space, and links depend on a ability characterizing multifractal attractors. specific criterion determining adjacency in phase space. The present paper is organized as follows: In Sec. II we Cycle networks [11], correlation networks [12, 13], and re- briefly summarize the multifractal formalism for chaotic currence networks [14–17] are typical examples of prox- time series and attractors. In Sec. III we present a tran- arXiv:1612.07070v1 [cond-mat.stat-mech] 21 Dec 2016 imity networks. Finally, the related class of transition sition network mapping for general time series, based on networks [18–21] encompasses different models for map- the box-counting algorithm [19]. In Sec. IV we relate the ping time series into networks, in which the values on topological properties of the associated networks with the the time signal are mapped into a finite number of states multifractal properties of the original time series. This (regions in the phase space of the series), representing relation is directly mediated by the natural measure of 2 the series, defined as the probability that the time se- the same value of α, such that their number scales as −f(α) quence visits a given box in a partition of the substrate Nα(ε) ∼ ε . The function f(α), called the multifrac- of the series. Numerical checks of our theoretical pre- tal spectrum of the measure, defines the fractal dimen- dictions are detailed in Sec. V. Finally, we present our sion of the set of boxes with the given value α, and is in conclusions in Sec. VI. general a convex function with a single maximum. An equivalent, and numerically simpler, description can be obtained from the generalized dimensions Dq, defined as II. MULTIFRACTAL TIME SIGNALS [22]

For time signals arising from strange (chaotic) attrac- P q 1 log i pi(ε) Dq = lim , (3) tors, it is common that different regions of the phase q − 1 ε→0 log ε space are differently visited, and chaotic orbits spend most of their time in a small region of the support under- which, for q ≥ 0, fulfill D0 ≤ Dq ≤ D∞ ≡ αmin. For neath the chaotic attractor itself. This heterogeneity is a uniform measure D = const. = D . For multifractal at the basis of the so-called multifractal structure of the q 0 measures, Dq is a decreasing function of q which is re- strange attractor, which can be mathematically captured lated with the multifractal spectrum (α, f(α)) by means by the formalism presented below [22, 24]. of a Legendre transformation [22], defining a parametric Let us consider a real time signal, or a d-dimensional exponent α(q) that fulfills the equations chaotic attractor, given by the normalized sequence of points G = ~x : ~x ∈ [0, 1]d, t = 1, 2, . . . n , where the in- t t d dex t represents either time or the order in the sequence α(q) = (q − 1)Dq (4) of points that generate the attractor (e.g. the index of dq the iteration in an iterated map), while n  1 is the f(α(q)) = qα(q) + (1 − q)Dq (5) number of points in the signal. We consider a partition d −1/d of the set [0, 1] in M boxes of length ε = M . Box Numerically, the generalized dimensions can be esti- are labeled by the indexes i, with 1 ≤ i ≤ M. Let us mated from Eq. (3), by noticing that, for finite ε, (1) (2) (d) associate to each point ~xt = {xt , xt , . . . , xt } in the sequence an integer index X (q−1)Dq pi(ε) ∼ ε , (6) d i X j (r) k r−1 it = 1 + xt L L (1) r=1 allowing Dq to be determined from a linear regression of P pi(ε) for decreasing values of ε (increasing M) in a 1/d i in the range 1 ≤ it ≤ M, where L = M , and bzc is the log-log plot. floor function. In this sense, the signal or attractor can be interpreted as visiting the it-th box in the partition at time t. In heterogeneous fractals, not all the boxes will be equally visited. In general, during the n steps of III. NEWORK MAPPING the signal, the i-th box will be visited a number of times n , and the total number of boxes visited at least once i In order to construct a transition network representa- will be N(ε). This quantity coincides with the number tion, we follow the approach in Refs. [19, 25], and asso- of boxes of length ε needed to cover the fractal set, and ciate a (virtual) vertex to each box 1 ≤ i ≤ M in the thus we can define the box or capacity dimension of the partition of the chaotic attractor in phase space. Actual attractor [23], D , by the relation 0 vertices in the network are given by the set of boxes that N(ε) ∼ ε−D0 . (2) have been visited at least once by the signal, with a size N(ε). Edges in the network are established by associ- 0 Let us define the probability pi(ε) = limn→∞ ni/n, ating an undirected connection between vertices i and i 0 termed the natural measure, as the probability that the whenever the signal jumps between boxes i and i in two 0 chaotic map visits the i-th box of the N(ε) available dur- consecutive time steps, i.e. i ≡ it and i ≡ it+1. ing an infinitely long orbit. For an homogeneous struc- The resulting projected networks, which are charac- d ture in d dimensions, pi(ε) ∼ ε , while in the case of a terized by the coarse-grained scale ε, are connected by D0 uniform fractal of dimension D0 ≤ d, pi(ε) ∼ ε . In construction, and preserve temporal information of the more complex situations, however, the attractor exhibits generator of the signal. A completely random, stochas- a non-uniform fractal distribution, and we assume a gen- tic signal will lead to a fully connected network; on the αi eral form for the natural measure pi(ε) ∼ ε , where the other hand, for a limit cycle or periodic attractor, the exponent αi, taking values in the interval [αmin, αmax], projected network will be ring, with a number of nodes measures the strength of the local singularity of the mea- equal to the period of the cycle. The former case corre- sure at box i. In general, there will be many boxes with sponds to D0 = d, and the latter to D0 = 1. 3

IV. RELATING TOPOLOGY WITH For a homogeneous fractal set, Dq = D0 ∀q, and MULTIFRACTALITY hg(k)qi = const. On the other hand, for a multifractal strange attractor, since Dq < D0 for q > 0, we have that In the case of a multifractal time series, the topology the moments hg(k)qi, with q > 1, diverge as the number of the associated transition network described above can of boxes in the partition increases, i.e. as the network be related to the generalized dimensions Dq. We will size grows. This observation allows to extract conclu- consider in particular the degree distribution Pε(k) of a sions on the functional form of the degree distribution, projected network with a coarse-graining level ε, defined which will depend on the particular growth law g(k). We as the probability that a randomly chosen node has de- will consider analytically tractable forms in the following gree k, i.e., it is connected to k other nodes. To make Section. To avoid complications in the development, we explicit this relation, we observe that every node (box) will further assume in our analysis that the degree distri- i, will be characterized by a number of visits ni and a bution of the projected network is stable, meaning that degree ki. In the limit n → ∞, we assume that the rel- the effect of ε consists essentially in imposing an upper ative number of visits, i.e. the natural measure pi(ε), degree cut-off kc to a functional form P (k) independent and the degree of a node are stable quantities. The cor- of ε. responding degree distribution will thus depend only on the phase space discretization ε. On average, the two A. Exponential growth quantities pi(ε) and ki will be related, since obviously the more times a box is visited, the larger is expected to be the degree of the associated node. We can thus Let us consider first the case of an exponential (faster assume a functional relation between these two averaged than algebraic) growth of the number of visits in a box quantities, valid for sufficiently large ki and pi, of the with the associated node degree, i.e. form g(k) ∼ eβk, (12) pi(ε) ' Gε(ki), (7) where β > 0. The fact that heβki is constant and heβqki where Gε(z) is an increasing function of z. The function diverges for q > 1 is compatible with a degree distri- Gε(ki) depends in general on ε. Indeed, from the nor- bution that shows, at large values of k, the asymptotic malization of the natural measure, Eq. (7) implies that behavior P (k) ∼ e−αk, with α > β, that is, an ex- X X ponentially bounded degree distribution. Assuming a Gε(ki) ≡ N(ε) Pε(k)Gε(k) = N(ε)hGε(k)i = 1, stable degree distribution, for a non-zero ε the diver- i k gence of the exponential moments heβqki will be reflected P in a dependence on the network size N(ε), modulated where we have defined hF (k)i = k Pε(k)F (k). From here, we have by the largest degree in the network, or degree cut-off kc(ε) [27]. To estimate this value, we observe that, from −1 D0 1/β hGε(k)i = N(ε) ∼ ε . (8)  D0   −(D0−αi)/β Eq. (9), ki ∼ ln pi(ε)/ε ∼ ln ε . The We make the assumption, to be validated numerically largest value of ki will correspond to the minimum of αi, later on, that the ε dependence of the function Gε(k) αmin = D∞. Therefore, we have resides exclusively in a multiplicative prefactor, i.e. h i −(D0−D∞)/β Gε(k) = a(ε)g(k). From Eq. (8), we have that hGε(k)i = kc(ε) ∼ ln ε . (13) a(ε)hg(k)i ∼ εD0 . We make the additional assumption that the average of g(k), hg(k)i, is a constant, indepen- This expression allows to relate the network parameters dent of ε. From here, we obtain the relation α and β with the multifractal exponents D0 and D∞ by noticing that, in a network of finite size N(ε), the max- D0 ∞ pi(ε) ' ε g(ki). (9) imum degree is given by the condition P P (k) = k=kc(ε) This relation implies that we can express the multifractal 1/N(ε) [27]. With an exponential degree distribution P (k) = αe−αk we thus have, in the continuous degree properties of the attractor in terms of topological prop- ∞ approximation, R αe−αk dk = e−αkc(ε) = 1/N(ε), erties of the network. In fact, from Eq. (9), we have kc(ε) q from where we obtain X q X  D0  qD0 X q pi(ε) ' ε g(k) ' ε N(ε) Pε(k)g(k) h i h i 1/α −D0/α i i k kc(ε) ∼ ln N(ε) ∼ ln ε . (14) ' ε(q−1)D0 hg(k)qi, (10) Combining Eqs. (13) and (14), we obtain the relation where we have used N(ε) ∼ ε−D0 . From Eq. (6) we have also P p (ε)q ∼ ε(q−1)Dq . Combining this with Eq. (10), α D0 i i = . (15) we obtain the relation linking network and multifractal β D0 − D∞ properties, namely The properties of the network can also be used to ex- hg(k)qi ∼ ε−(q−1)(D0−Dq ). (11) tract information on the full set of generalized dimension 4 by building on relation Eq. (11). Indeed we can write the Finally, writing diverging moments in a finite network as Z kc(ε) δq δq−γ δq−γ+1 k (ε) hk i ∼ k dk ∼ k (ε) Z c e(qβ−α)kc(ε) c heβqki ∼ e(qβ−α)k dk ∼ qβ − α ∼ ε−D0[δq/(γ−1)−1], (23) (16) and comparing with Eq. (11), we obtain again the sim- which, taking into account Eq. (13), yields ple asymptotic expression for the generalized dimensions given by Eq. (18).

heβqki ∼ ε−(D0−D∞)(q−α/β) . (17) V. NUMERICAL EXPERIMENTS Combining Eq. (17) and Eq. (11) leads to the asymp- totic expression, valid for large q In order to check the validity of the predictions made q in Section IV, we have considered different multifractal D ∼ D , (18) q ∞ q − 1 time signals, generated by means of iterative maps. In particular we have studied three paradigmatic examples that recovers the result known for deterministic multi- of 1- and 2-dimensional chaotic attractors, namely the fractal measures [24]. logistic, the Duffing and the Henon map. The well-known logistic recurrence [6, 28] in dimension d = 1,

B. Algebraic growth xt+1 = µ xt(1 − xt) , (24) maps the interval x ∈ [0, 1] into itself when the control In the case of an algebraic growth of the number of parameter µ ranges between 0 and 4. This systems un- visits with the associated degree, we have dergoes a period-doubling bifurcation transition to chaos, g(k) ∼ kδ, (19) which sets-in at µ = 3.569456... . Multifractal chaotic regimes interspersed with periodic windows then occur in with δ > 0. From Eq. (11), we have that hkδi is finite, and the parameter interval µ ∈ [3.57, 4). Here we fix µ = 3.7. hkqδi, with q > 1, diverge in the limit of infinite network The two-dimensional Duffing map size. The fact that all higher degree moments diverge (1) (2) indicate that the degree distribution P (k) of the network xt+1 = xt , (25) 3 has long tails, which in the simplest case are compatible (2) (1) (2) h (2)i x = −b x + a x − x , (26) with a power law distribution of the form P (k) ∼ k−γ , t+1 t t t where we impose γ > δ + 1 to ensure a finite value of is a discrete representation of the Duffing oscillator, de- hkδi. scribing a forced oscillator coupled to a dissipative restor- Performing again a finite-size analysis, from Eq. (9) we ing force [6]. This map typically produces chaotic be- obtain a maximum degree haviours with the critical parameters a = 2.75 and (1) (2) −(D −D )/δ (1−D /D )/δ b = 0.2, and generates values of x and x in the k (ε) ∼ ε 0 ∞ ∼ N(ε) ∞ 0 . (20) c range [−1.71, 1.71]. Following our network projection al- On the other hand, the network relation P∞ P (k) = gorithm, each variable is thus shifted and normalized into k=kc(ε) the interval [0, 1]2. Finally, the Henon map [6, 29] in 1/N(ε) leads now, with a power-law degree distribution d = 2 is defined by the recurrence P (k) ∼ k−γ in the continuous degree approximation, to 2 [27] (1) h (1)i (2) xt+1 = 1 − a xt + xt , (27) 1/(γ−1) −D0/(γ−1) kc(ε) ∼ N(ε) ∼ ε . (21) (2) (2) xt+1 = b xt , (28) Combining Eqs. (20) and (21), we obtain the relation with the parameters a and b fixed to 1.4 and 0.3, respec- between network properties and multifractal exponents tively. For these values, starting from an initial point (1) (2) γ − 1 D0 (x0 , x0 ) the dynamics can either asymptotize to a frac- = . (22) tal attractor relying on the subset x(1) ∈ [−1.3, 1.3] and δ D0 − D∞ x(2) ∈ [−0.4, 0.4] or diverge to infinity. For other values In the case of an algebraic g(k) function, Eqs. (20) of a and b the map may be still chaotic, intermittent, or and (21) can be used to directly estimate δ and γ. This converge to a periodic orbit. Here we use the classical approach is more difficult in the case of an exponential parameter setting and, again, we transform the values of g(k), due to the much smaller range of variation of the each variable into the interval [0, 1]2. A one-dimensional logarithm of ε, see Eqs. (13) and (14). projection of the Henon map is also considered, obtained

D∞ = 0.96(1)

5 a) 0.51 b) a) 104 b) 102 0.79 δ = 1.04(2) β = 0.28(8) 0.505 101

) 2 ) q q 10 k

k 0

0.5 0.78 ( (

D 10 D

g 3 5 M = 10 g M = 2.5 x 10 104 106 0.495 10−1 0.77 100 105 2.25 x 106 D∞ = 0.49(1) D∞ = 0.75(1) 106 4 x 106 0.49 −2 1.01 1.02 1.03 1.04 1.01 1.02 1.03 1.04 10 100 101 102 103 104 0 10 20 30 q/(q-1) q/(q-1) k k 2 c) 0.78 d) 0.84 c) d) 10 102 δ = 1.10(5) β = 0.15(2) 0.775 101 0.83 ) ) q q

0 k 0 k 10 0.77 10 ( D ( D

3 g 4 0.82 g M = 10 M = 10 0.765 104 10−1 2.5 x 105 10−2 5 6 D∞ = 0.73(2) D∞ = 0.82(2) 10 10 106 −2 4 x 106 0.76 0.81 10 1.01 1.02 1.03 1.04 1.01 1.02 1.03 1.04 100 101 102 103 104 0 10 20 30 q/(q-1) q/(q-1) k k

FIG. 1. Generalized dimensions Dq as a function of q/(q − 1) FIG. 2. Rescaled average natural measure ε−D0 p¯(k) as a func- for the logistic map a), the Duffing map b), the d = 1 Henon tion of the degree k for the logistic map a), the Duffing map map c) and the d = 2 Henon map c). The corresponding b), the d = 1 Henon map c) and the d = 2 Henon map d). asymptotic values of D∞ can be extrapolated by means of a The functionp ¯(k) is computed for different values of M, con- linear regression of these plots following Eq. (18). The results sidering n = 1000 × M iterations. Dashed lines represent the are listed in Table I. estimated values of δ and β in each case, see Table I. by taking into account only one normalized variable of the map (both x(1) and x(2) give analogous results). g(k). From Fig. 2 we observe, for the logistic, Duffing and In Table I we present a summary of the multifractal d = 2 Henon maps, a perfect convergence of the rescaled properties of the chaotic time signals considered, com- average natural measure as a function of k, indicating puted by using the box counting formalism described in the validity of the assumption in Eq. (9). In the case of Sec. II. In particular, from Eq. (6), we compute the expo- the d = 1 projection of the Henon attractor, however, P the collapse ofp ¯ε(k) is not fulfilled. In this particular nent Dq by performing a linear regression of log i pi(ε) as a function of log ε, for ε = M 1/d, in a range of values of case, therefore, the predictions made in Sec. IV are not M between 103 and 4 × 106, depending of the particular expected to hold. attractor. The slope of this regression yields the factor From Fig. 2, we also observe that in the d = 1 cases (q − 1)Dq. According to Eq. (18), the asymptotic value (Figs. 2a and 2c), both the attractor of the logistic map D∞ is obtained by means of linear regressions of Dq as and that of the projected Henon map obey an algebraic a function of q/(q − 1), performed over suitable intervals behavior, g(k) ∼ kδ, while in the d = 2 cases (Figs. 2b of the variable q/(q − 1), see Fig. 1. and 2d), both the Duffing and the Henon systems show an exponential growth, g(k) ∼ eβk. A linear regression In Fig. 2 we plot the natural measurep ¯ε(k), averaged over all nodes of degree k, as a function of k, for different performed on the data in Fig. 2 provides an estimation of partitions of the multifractal attractors, i.e. different ε the exponents δ and β. In general, it can be noticed how (or M). In order to check the main assumption in Eq. (9), the linear profiles become more defined and stable while −D0 refining the statistics of g(k) by increasing M; we thus we plot the rescaled function ε p¯ε(k) as a function of fit data obtained with M = 106 in the algebraic case and the degree, using the fractal dimensions D0 quoted in 6 Table I. In this case, we expect all plots of each map for M = 4 × 10 for an exponential g(k). We thus find the different ε to collapse onto the single universal function exponents δ = 1.04(2) for the logistic map, β = 0.28(8) for the Duffing map, while δ = 1.10(5) and β = 0.15(2) are obtained for the d = 1 and the d = 2 Henon maps, respectively. Attractor d D0 D∞ β δ α γ Logistic 1 0.998(2) 0.49(1) — 1.04(2) — 3.14(1) In Fig. 3 we examine the topology of the projected net- works by plotting the cumulative degree distributions, Duffing 2 1.306(2) 0.75(1) 0.28(8) — 0.42(6) — P∞ Pcum(k) = P (q), for different values of M. As 1 1.000(1) 0.73(2) — 1.10(5) — 4.48(1) q=k Henon predicted in Sec. IV, the networks characterized by an 2 1.24(1) 0.82(2) 0.15(2) — 0.57(5) — algebraic growth g(k) (here d = 1 cases) exhibit power- law degree distributions. Panels a) and c) show how TABLE I. Properties of the different multifractal time signals for different network sizes (i.e. different values of M) and associated projected networks. all trends converge to a common power-law distribution 6 a) b) 3 4 a) b) 0 M = 10 M = 10 4 10 0 3 104 10 2.5 x 105 5 6 ) ) 10 10 ) ) ) ) −2 6 6 ε ε ( ( k 10 10 k 4 x 10 2.5 c c ( ( 3 −2 m m 10 k k u u ( (

c −4 c

10 g 2 g P P o o l s = 0.48(1) l 2 s = 0.54(1) 10−4 10−6 γ = 3.143(2) α = 0.42(6) 1.5

100 101 102 103 10 20 30 2 3 4 5 6 3 4 5 6 k k log(N(ε)) log(N(ε)) c) 3 d) 4 0 M = 10 M = 10 10 4 100 5 10 2.5 x 10 FIG. 4. Scaling of the degree cut-off, kc(), estimated as the 105 106 ) −2 6 ) 6 maximum degree in the network, as function of the network

k 10 10 k 4 x 10 ( ( 10−2

m m size N() for the logistic map, a), and the d = 1 Henon map, u u c 10−4 c b). P P 10−4 10−6 γ = 4.480(3) α = 0.57(5) 101 102 103 104 10 20 30 tic map and the d = 1 projection of the Henon attractor, k k we can proceed to check the behavior of the maximum degree k (ε) as a function of the network size N(ε). We FIG. 3. Cumulative degree distributions of the networks pro- c jected from the logistic map a), the Duffing map b), the d = 1 do not consider the relation between degree cut-off and Henon map c) and the d = 2 Henon map d), computed for dif- network size for maps with an exponential g(k), since the ferent values of M and with n = 1000×M iterations. Dashed very small span of network sizes obtained does not allow lines represent the estimated values of γ and/or α of the per- for a determination of the exponent in relation Eq. (13). tinent case, see Table I. Indeed, for the algebraic case, following Eq. (20), the be- s havior of the maximum degree is given by kc(ε) ∼ N(ε) , with an exponent s = (1−D∞/D0)/δ. The numerical ex- characterized by γ = 3.14(1) in the logistic networks and ponents obtained through a linear regression of log kc(ε) γ = 4.48(1) for the d = 1 Henon map. Interestingly, in as a function of log N(ε) are s = 0.48(1) and s = 0.54(1) this last case, the prediction of a scale-free degree distri- for the logistic and the d = 1 Henon map, respectively. bution holds, despite the fact that Eq. (9) is not fulfilled. From the values of δ, D0 and D∞ in Table I, our theoreti- This must attributed to the affect of an algebraic func- cal predictions are s = 0.48(3) and s = 0.25(3) for the lo- tion g(k), still present in the d = 1 Henon map. By gistic and the d = 1 Henon map, respectively. The agree- contrast, the networks resulting from the projection of ment between numerics and theory is very good for the the d = 2 maps follow short-tailed degree distributions, logistic map, but completely off in the d = 1 Henon case. compatible with an exponential behavior P (k) ∼ e−αk. The disagreement in this last case must be attributed to These plots present a poor statistics and we can only ex- the failure of Eq. (9). While the general trend towards a trapolate rough values for α, namely α = 0.42(6) and power-law degree distribution is ensured by the algebraic α = 0.57(5) for the Duffing and Henon maps, respec- form of the function g(k), the lack of the expected scal- tively. ing with ε (i.e. the prefactor) affects the scaling relations With these values characterizing the multifractal prop- deduced from Eq. (9). In the d = 1 Henon case yet an- erties of the maps and the topological properties obtained other of the approximations made breaks down, namely from the analysis of the projected transition networks, we the assumption of a stable degree distribution. This fact can validate our theoretical framework. We first check is checked in Fig. 5, which shows that the average degree the cross-relations given by Eqs. (15) and (22) for expo- hki of the projected networks is essentially independent of nential and algebraic cases, respectively. Regarding the ε for all the multifractal time signals considered, except algebraic examples we obtain (γ − 1)/δ = 2.06(7) and for the d = 1 Henon attractor, which exhibits instead a power-law increasing behavior. This fact indicates that D0/(D0 − D∞) = 1.96(5) for the logistic map (see val- ues from Table I). In this case, the identity Eq. (22) is d = 1 Henon networks belong to the class of accelerated well fulfilled within error bars. For the Henon attrac- networks [30]. This increasing average degree causes the tor in d = 1, we obtain instead (γ − 1)/δ = 3.2(2) and degree distribution to be non stable, introducing an ad- ditional scale, beyond the degree cut-off k . Surprisingly, D0/(D0 − D∞) = 3.7(3), again coinciding within error c bars. For exponential cases, the equality Eq. (15) leads however, the exponent relation Eq. (22) seems to still be fulfilled, at least within error bars. to the comparison α/β = 1.5(6) and D0/(D0 − D∞) = 2.34(5) for the Duffing map, and α/β = 3.8(8) and D0/(D0 − D∞) = 2.9(2) for the d = 2 Henon attrac- tor. In this case, the exponent relations are affected by VI. DISCUSSION stronger errors, but the trend is clearly towards a positive comparison. In this paper we have investigated the effect of the Finally, for an algebraic g(k), as observed in the logis- fractal and multifractal properties of a temporal signal on 7

the multifractal exponents of the generating signal. We 102 have verified the validity of our approach through ex- tensive numerical simulations, highlighting the excellent agreement observed in many cases, and discussing in de- tail the reasons why in some cases (e.g., the Henon map 1 10 in d = 1) the numerical experiments depart from theo- >

k retical predictions.

< In particular, we could conclude that a sufficient con- Henon d=1 100 dition to obtain a scale-free topology is that the natural Henon d=2 measure of a box must increase with the degree of the Duffing associated node in a algebraic fashion. In our numeri- Logistic cal experiments we have observed that this condition is −1 10 fulfilled in multifractal attractors in d = 1 with fractal 103 104 105 106 107 dimension D0 = 1. This fact leads us to conjecture that M scale-free networks can be observed in general multifrac- tal time series in which the fractal dimension is equal to FIG. 5. Average degree hki as a function of M for the different the euclidean dimension of the embedding phase space. maps considered. All of them, except the d = 1 projection Our work extends existing approaches bridging time of the Henon map are essentially independent of the network series analysis and network science by addressing the discretization. ubiquitous case of signals exhibiting multifractal prop- erties. By doing this, it enriches the set of interpretative the topology of the corresponding projected network. By tools available for a better characterization of empirical combining a transition network representation with the time-series. For this reason, we can envisage that it will box counting formalism, we have mapped temporal sig- be of interest also to the growing community of interdisci- nals into networks whose nodes are the boxes partitioning plinary researchers studying natural time series through the attractor of the temporal signal in the phase space, the lenses of network science. and links are established between successive pair of boxes between which the signal jumps. We have developed a mathematical framework connecting network topology to ACKNOWLEDGMENTS the multifractal properties of the generating signal. This formalism allows us to predict the functional form of the M.A.B. acknowledges financial support from FRS- network degree distribution on the basis of the relation, FNRS. R.P.-S. acknowledges financial support from the g(k), linking the natural measure of a box with the as- Spanish MINECO, under projects FIS2013-47282-C2- sociated node degree. We have focused on the prototyp- 2 and FIS2016-76830-C2-1-P, and EC FET-Proactive ical and general cases of an exponential and an algebraic Project MULTIPLEX (Grant No. 317532). R.P.-S. growth g(k), showing that the latter results in power law acknowledges additional financial support from ICREA degree distributions whose exponent γ is controlled by Academia, funded by the Generalitat de Catalunya.

[1] M. E. J. Newman, Networks: An introduction (Oxford Letters 29, 048903 (2012). University Press, Oxford, 2010). [9] L. Lacasa, B. Luque, F. Ballesteros, J. Luque, and J. C. [2] A. Barrat, M. Barth´elemy, and A. Vespignani, Dynami- Nu˜no,Proc. Natl. Acad. Sci. USA 105, 4972 (2008). cal Processes on Complex Networks (Cambridge Univer- [10] L. Lacasa and R. Toral, Phys. Rev. E 82, 036120 (2010). sity Press, Cambridge, 2008). [11] J. Zhang and M. Small, Phys. Rev. Lett. 96, 238701 [3] S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, (2006). Rev. Mod. Phys. 80, 1275 (2008). [12] Y. Yang and H. Yang, Physica A: Statistical Mechanics [4] R. Pastor-Satorras, C. Castellano, P. Van Mieghem, and and its Applications 387, 1381 (2008). A. Vespignani, Rev. Mod. Phys. , 925 (2015). [13] Z. Gao and N. Jin, Phys. Rev. E 79, 066303 (2009). [5] J. Guckenheimer and P. Holmes, Nonlinear oscillations, [14] N. Marwan, J. F. Donges, Y. Zou, R. V. Donner, and dynamical systems, and bifurcations of vector fields, Ap- J. Kurths, Physics Letters A 373, 4246 (2009). plied mathematical sciences (Springer, New York, 2002). [15] Z. Gao and N. Jin, Chaos 19, 033137 (2009). [6] E. Ott, Chaos in dynamical systems (Cambridge Univer- [16] R. V. Donner, Y. Zou, J. F. Donges, N. Marwan, and sity Press, Cambridge, U.K., 1993). J. Kurths, New Journal of Physics 12, 033025 (2010). [7] R. V. Donner, M. Small, J. F. Donges, N. Marwan, [17] R. Xiang, J. Zhang, X.-K. Xu, and M. Small, Chaos 22, Y. Zou, R. Xiang, and J. Kurths, International Jour- 013107 (2012). nal of Bifurcation and Chaos 21, 1019 (2011). [18] G. Nicolis, A. Garcia Cantu, and C. Nicolis, Interna- [8] P. Li, J. Zhang, X.-K. Xu, and M. Small, Chinese Physics tional Journal of Bifurcation and Chaos 15, 3467 (2005). 8

[19] A. H. Shirazi, G. R. Jafari, J. Davoudi, J. Peinke, [24] J. Feder, Fractals, Physics of Solids and Liquids M. R. R. Tabar, and M. Sahimi, Journal of Statisti- (Springer, New York, 1989). cal Mechanics: Theory and Experiment 2009, P07046 [25] M. A. Budroni, E. Tiezzi, and M. Rustici, Physica (2009). A: Statistical Mechanics and its Applications 389, 3883 [20] A. Campanharo, M. Sirer, R. Malmgren, F. Ramos, and (2010). L. Amaral, PLoS ONE 6, e23378 (2011). [26] A.-L. Barab´asiand R. Albert, Science 286, 509 (1999). [21] X. Sun, M. Small, Y. Zhao, and X. Xue, Chaos 24, [27] M. Bogu˜n´a, R. Pastor-Satorras, and A. Vespignani, 024402 (2014). Euro. Phys. J. B 38, 205 (2004). [22] T. C. Halsey, M. H. Jensen, L. P. Kadanoff, I. Procaccia, [28] R. M. May, Nature 261, 459 (1976). and B. I. Shraiman, Phys. Rev. A 33, 1141 (1986). [29] M. H´enon,Communications in Mathematical Physics 50, [23] K. Falconer, Fractal Geometry: Mathematical Founda- 69 (1976). tions and Applications, 2nd ed. (John Wiley & Sons, [30] S. N. Dorogovtsev and J. F. F. Mendes, Advances in Chichester, 2003). Physics 51, 1079 (2002).