-
Chaotic renormalization group flow and entropy gradients over Haros graphs
Authors:
Jorge Calero-Sanz,
Bartolo Luque,
Lucas Lacasa
Abstract:
Haros graphs have been recently introduced as a set of graphs bijectively related to real numbers in the unit interval. Here we consider the iterated dynamics of a graph operator $\cal R$ over the set of Haros graphs. This operator was previously defined in the realm of graph-theoretical characterisation of low-dimensional nonlinear dynamics, and has a renormalization group (RG) structure. We find…
▽ More
Haros graphs have been recently introduced as a set of graphs bijectively related to real numbers in the unit interval. Here we consider the iterated dynamics of a graph operator $\cal R$ over the set of Haros graphs. This operator was previously defined in the realm of graph-theoretical characterisation of low-dimensional nonlinear dynamics, and has a renormalization group (RG) structure. We find that the dynamics of $\cal R$ over Haros graphs is complex and includes unstable periodic orbits or arbitrary period and non-mixing aperiodic orbits, overall portraiting a chaotic RG flow. We identify a single RG stable fixed point whose basin of attraction is the set of rational numbers, associate periodic RG orbits with (pure) quadratic irrationals and aperiodic RG orbits with (non-mixing) families of non-quadratic algebraic irrationals and trascendental numbers. Finally, we show that the entropy gradients inside periodic RG orbits are constant. We discuss the possible physical interpretation of such chaotic RG flow and speculate on the entropy-constant periodic orbits as a possible confirmation of a (quantum field-theoretic) $c$-theorem applied inside the invariant set of a RG flow.
△ Less
Submitted 13 December, 2022;
originally announced December 2022.
-
Haros graphs: an exotic representation of real numbers
Authors:
Jorge Calero-Sanz,
Bartolo Luque,
Lucas Lacasa
Abstract:
This paper introduces Haros graphs, a construction which provides a graph-theoretical representation of real numbers in the unit interval reached via paths in the Farey binary tree. We show how the topological structure of Haros graphs yields a natural classification of the reals numbers into a hierarchy of families. To unveil such classification, we introduce an entropic functional on these graph…
▽ More
This paper introduces Haros graphs, a construction which provides a graph-theoretical representation of real numbers in the unit interval reached via paths in the Farey binary tree. We show how the topological structure of Haros graphs yields a natural classification of the reals numbers into a hierarchy of families. To unveil such classification, we introduce an entropic functional on these graphs and show that it can be expressed, thanks to its fractal nature, in terms of a generalised de Rham curve. We show that this entropy reaches a global maximum at the reciprocal of the Golden number and otherwise displays a rich hierarchy of local maxima and minima that relate to specific families of irrationals (noble numbers) and rationals, overall providing an exotic classification and representation of the reals numbers according to entropic principles. We close the paper with a number of conjectures and outline a research programme on Haros graphs.
△ Less
Submitted 7 July, 2022;
originally announced July 2022.
-
Enhancing self-supervised monocular depth estimation with traditional visual odometry
Authors:
Lorenzo Andraghetti,
Panteleimon Myriokefalitakis,
Pier Luigi Dovesi,
Belen Luque,
Matteo Poggi,
Alessandro Pieropan,
Stefano Mattoccia
Abstract:
Estimating depth from a single image represents an attractive alternative to more traditional approaches leveraging multiple cameras. In this field, deep learning yielded outstanding results at the cost of needing large amounts of data labeled with precise depth measurements for training. An issue softened by self-supervised approaches leveraging monocular sequences or stereo pairs in place of exp…
▽ More
Estimating depth from a single image represents an attractive alternative to more traditional approaches leveraging multiple cameras. In this field, deep learning yielded outstanding results at the cost of needing large amounts of data labeled with precise depth measurements for training. An issue softened by self-supervised approaches leveraging monocular sequences or stereo pairs in place of expensive ground truth depth annotations. This paper enables to further improve monocular depth estimation by integrating into existing self-supervised networks a geometrical prior. Specifically, we propose a sparsity-invariant autoencoder able to process the output of conventional visual odometry algorithms working in synergy with depth-from-mono networks. Experimental results on the KITTI dataset show that by exploiting the geometrical prior, our proposal: i) outperforms existing approaches in the literature and ii) couples well with both compact and complex depth-from-mono architectures, allowing for its deployment on high-end GPUs as well as on embedded devices (e.g., NVIDIA Jetson TX2).
△ Less
Submitted 12 August, 2019; v1 submitted 8 August, 2019;
originally announced August 2019.
-
On a dynamical approach to some prime number sequences
Authors:
Lucas Lacasa,
Bartolo Luque,
Ignacio Gómez,
Octavio Miramontes
Abstract:
In this paper we show how the cross-disciplinary transfer of techniques from Dynamical Systems Theory to Number Theory can be a fruitful avenue for research. We illustrate this idea by exploring from a nonlinear and symbolic dynamics viewpoint certain patterns emerging in some residue sequences generated from the prime number sequence. We show that the sequence formed by the residues of the primes…
▽ More
In this paper we show how the cross-disciplinary transfer of techniques from Dynamical Systems Theory to Number Theory can be a fruitful avenue for research. We illustrate this idea by exploring from a nonlinear and symbolic dynamics viewpoint certain patterns emerging in some residue sequences generated from the prime number sequence. We show that the sequence formed by the residues of the primes modulo $k$ are maximally chaotic and, while lacking forbidden patterns, display a non-trivial spectrum of Renyi entropies which suggest that every block of size $m>1$, while admissible, occurs with different probability. This non-uniform distribution of blocks for $m>1$ contrasts Dirichlet's theorem that guarantees equiprobability for $m=1$. We then explore in a similar fashion the sequence of prime gap residues. This sequence is again chaotic (positivity of Kolmogorov-Sinai entropy), however chaos is weaker as we find forbidden patterns for every block of size $m>1$. We relate the onset of these forbidden patterns with the divisibility properties of integers, and estimate the densities of gap block residues via Hardy-Littlewood $k$-tuple conjecture. We use this estimation to argue that the amount of admissible blocks is non-uniformly distributed, what supports the fact that the spectrum of Renyi entropies is again non-trivial in this case. We complete our analysis by applying the Chaos Game to these symbolic sequences, and comparing the IFS attractors found for the experimental sequences with appropriate null models.
△ Less
Submitted 22 February, 2018;
originally announced February 2018.
-
Emergence of linguistic laws in human voice
Authors:
Ivan Gonzalez Torre,
Bartolo Luque,
Lucas Lacasa,
Jordi Luque,
Antoni Hernandez-Fernandez
Abstract:
Linguistic laws constitute one of the quantitative cornerstones of modern cognitive sciences and have been routinely investigated in written corpora, or in the equivalent transcription of oral corpora. This means that inferences of statistical patterns of language in acoustics are biased by the arbitrary, language-dependent segmentation of the signal, and virtually precludes the possibility of mak…
▽ More
Linguistic laws constitute one of the quantitative cornerstones of modern cognitive sciences and have been routinely investigated in written corpora, or in the equivalent transcription of oral corpora. This means that inferences of statistical patterns of language in acoustics are biased by the arbitrary, language-dependent segmentation of the signal, and virtually precludes the possibility of making comparative studies between human voice and other animal communication systems. Here we bridge this gap by proposing a method that allows to measure such patterns in acoustic signals of arbitrary origin, without needs to have access to the language corpus underneath. The method has been applied to six different human languages, recovering successfully some well-known laws of human communication at timescales even below the phoneme and finding yet another link between complexity and criticality in a biological system. These methods further pave the way for new comparative studies in animal communication or the analysis of signals of unknown code.
△ Less
Submitted 9 October, 2016;
originally announced October 2016.
-
Canonical Horizontal Visibility Graphs are uniquely determined by their degree sequence
Authors:
Bartolo Luque,
Lucas Lacasa
Abstract:
Horizontal visibility graphs (HVGs) are graphs constructed in correspondence with number sequences that have been introduced and explored recently in the context of graph-theoretical time series analysis. In most of the cases simple measures based on the degree sequence (or functionals of these such as entropies over degree and joint degree distributions) appear to be highly informative features f…
▽ More
Horizontal visibility graphs (HVGs) are graphs constructed in correspondence with number sequences that have been introduced and explored recently in the context of graph-theoretical time series analysis. In most of the cases simple measures based on the degree sequence (or functionals of these such as entropies over degree and joint degree distributions) appear to be highly informative features for automatic classification and provide nontrivial information on the associated dynam- ical process, working even better than more sophisticated topological metrics. It is thus an open question why these seemingly simple measures capture so much information. Here we prove that, under suitable conditions, there exist a bijection between the adjacency matrix of an HVG and its degree sequence, and we give an explicit construction of such bijection. As a consequence, under these conditions HVGs are unigraphs and the degree sequence fully encapsulates all the information of these graphs, thereby giving a plausible reason for its apparently unreasonable effectiveness.
△ Less
Submitted 17 May, 2016;
originally announced May 2016.
-
Walking on exoplanets: Is Star Wars right?
Authors:
Fernando J. Ballesteros,
Bartolo Luque
Abstract:
As the number of detected extrasolar planets increases, exoplanet databases become a valuable resource, confirming some details about planetary formation, but also challenging our theories with new unexpected properties.
As the number of detected extrasolar planets increases, exoplanet databases become a valuable resource, confirming some details about planetary formation, but also challenging our theories with new unexpected properties.
△ Less
Submitted 26 April, 2016;
originally announced April 2016.
-
Speech earthquakes: scaling and universality in human voice
Authors:
Jordi Luque,
Bartolo Luque,
Lucas Lacasa
Abstract:
Speech is a distinctive complex feature of human capabilities. In order to understand the physics underlying speech production, in this work we empirically analyse the statistics of large human speech datasets ranging several languages. We first show that during speech the energy is unevenly released and power-law distributed, reporting a universal robust Gutenberg-Richter-like law in speech. We f…
▽ More
Speech is a distinctive complex feature of human capabilities. In order to understand the physics underlying speech production, in this work we empirically analyse the statistics of large human speech datasets ranging several languages. We first show that during speech the energy is unevenly released and power-law distributed, reporting a universal robust Gutenberg-Richter-like law in speech. We further show that such earthquakes in speech show temporal correlations, as the interevent statistics are again power-law distributed. Since this feature takes place in the intra-phoneme range, we conjecture that the responsible for this complex phenomenon is not cognitive, but it resides on the physiological speech production mechanism. Moreover, we show that these waiting time distributions are scale invariant under a renormalisation group transformation, suggesting that the process of speech generation is indeed operating close to a critical point. These results are put in contrast with current paradigms in speech processing, which point towards low dimensional deterministic chaos as the origin of nonlinear traits in speech fluctuations. As these latter fluctuations are indeed the aspects that humanize synthetic speech, these findings may have an impact in future speech synthesis technologies. Results are robust and independent of the communication language or the number of speakers, pointing towards an universal pattern and yet another hint of complexity in human speech.
△ Less
Submitted 5 August, 2014;
originally announced August 2014.
-
On the thermodynamic origin of metabolic scaling
Authors:
Fernando J. Ballesteros,
Vicent J. Martínez,
Bartolo Luque,
Lucas Lacasa,
Enric Valor,
Andrés Moya
Abstract:
The origin and shape of metabolic scaling has been controversial since Kleiber found that basal metabolic rate of animals seemed to vary as a power law of their body mass with exponent 3/4, instead of 2/3, as a surface-to-volume argument predicts. The universality of exponent 3/4 -claimed in terms of the fractal properties of the nutrient network- has recently been challenged according to empirica…
▽ More
The origin and shape of metabolic scaling has been controversial since Kleiber found that basal metabolic rate of animals seemed to vary as a power law of their body mass with exponent 3/4, instead of 2/3, as a surface-to-volume argument predicts. The universality of exponent 3/4 -claimed in terms of the fractal properties of the nutrient network- has recently been challenged according to empirical evidence that observed a wealth of robust exponents deviating from 3/4. Here we present a conceptually simple thermodynamic framework, where the dependence of metabolic rate with body mass emerges from a trade-off between the energy dissipated as heat and the energy efficiently used by the organism to maintain its metabolism. This balance tunes the shape of an additive model from which different effective scalings can be recovered as particular cases, thereby reconciling previously inconsistent empirical evidence in mammals, birds, insects and even plants under a unified framework. This model is biologically motivated, fits remarkably well the data, and also explains additional features such as the relation between energy lost as heat and mass, the role and influence of different climatic environments or the difference found between endotherms and ectotherms.
△ Less
Submitted 16 January, 2018; v1 submitted 14 July, 2014;
originally announced July 2014.
-
Quasiperiodic graphs at the onset of chaos
Authors:
Bartolo Luque,
Marta Cordero-Gracia,
Mariola Gómez,
Alberto Robledo
Abstract:
We examine the connectivity fluctuations across networks obtained when the horizontal visibility (HV) algorithm is used on trajectories generated by nonlinear circle maps at the quasiperiodic transition to chaos. The resultant HV graph is highly anomalous as the degrees fluctuate at all scales with amplitude that increases with the size of the network. We determine families of Pesin-like identitie…
▽ More
We examine the connectivity fluctuations across networks obtained when the horizontal visibility (HV) algorithm is used on trajectories generated by nonlinear circle maps at the quasiperiodic transition to chaos. The resultant HV graph is highly anomalous as the degrees fluctuate at all scales with amplitude that increases with the size of the network. We determine families of Pesin-like identities between entropy growth rates and generalized graph-theoretical Lyapunov exponents. An irrational winding number with pure periodic continued fraction characterizes each family. We illustrate our results for the so-called golden, silver and bronze numbers.
△ Less
Submitted 11 December, 2013;
originally announced December 2013.
-
Phase transitions in Number Theory: from the Birthday Problem to Sidon Sets
Authors:
Bartolo Luque,
Ivan G. Torre,
Lucas Lacasa
Abstract:
In this work, we show how number theoretical problems can be fruitfully approached with the tools of statistical physics. We focus on g-Sidon sets, which describe sequences of integers whose pairwise sums are different, and propose a random decision problem which addresses the probability of a random set of k integers to be g-Sidon. First, we provide numerical evidence showing that there is a cros…
▽ More
In this work, we show how number theoretical problems can be fruitfully approached with the tools of statistical physics. We focus on g-Sidon sets, which describe sequences of integers whose pairwise sums are different, and propose a random decision problem which addresses the probability of a random set of k integers to be g-Sidon. First, we provide numerical evidence showing that there is a crossover between satisfiable and unsatisfiable phases which converts to an abrupt phase transition in a properly defined thermodynamic limit. Initially assuming independence, we then develop a mean field theory for the g-Sidon decision problem. We further improve the mean field theory, which is only qualitatively correct, by incorporating deviations from independence, yielding results in good quantitative agreement with the numerics for both finite systems and in the thermodynamic limit. Connections between the generalized birthday problem in probability theory, the number theory of Sidon sets and the properties of q-Potts models in condensed matter physics are briefly discussed.
△ Less
Submitted 12 October, 2013;
originally announced October 2013.
-
Horizontal Visibility graphs generated by type-I intermittency
Authors:
Ángel M. Núñez,
Bartolo Luque,
Lucas Lacasa,
José Patricio Gómez,
Alberto Robledo
Abstract:
The type-I intermittency route to (or out of) chaos is investigated within the Horizontal Visibility graph theory. For that purpose, we address the trajectories generated by unimodal maps close to an inverse tangent bifurcation and construct, according to the Horizontal Visibility algorithm, their associated graphs. We show how the alternation of laminar episodes and chaotic bursts has a fingerpri…
▽ More
The type-I intermittency route to (or out of) chaos is investigated within the Horizontal Visibility graph theory. For that purpose, we address the trajectories generated by unimodal maps close to an inverse tangent bifurcation and construct, according to the Horizontal Visibility algorithm, their associated graphs. We show how the alternation of laminar episodes and chaotic bursts has a fingerprint in the resulting graph structure. Accordingly, we derive a phenomenological theory that predicts quantitative values of several network parameters. In particular, we predict that the characteristic power law scaling of the mean length of laminar trend sizes is fully inherited in the variance of the graph degree distribution, in good agreement with the numerics. We also report numerical evidence on how the characteristic power-law scaling of the Lyapunov exponent as a function of the distance to the tangent bifurcation is inherited in the graph by an analogous scaling of the block entropy over the degree distribution. Furthermore, we are able to recast the full set of HV graphs generated by intermittent dynamics into a renormalization group framework, where the fixed points of its graph-theoretical RG flow account for the different types of dynamics. We also establish that the nontrivial fixed point of this flow coincides with the tangency condition and that the corresponding invariant graph exhibit extremal entropic properties.
△ Less
Submitted 21 January, 2013;
originally announced January 2013.
-
Nestedness in mutualistic networks
Authors:
Rudolf P. Rohr,
Miguel A. Fortuna,
Bartolo Luque,
Jordi Bascompte
Abstract:
James et al. (2012) presented simulations that apparently falsify the analytical result by Bastolla et al. (2009), who showed that nested mutualistic interactions decrease interspecific competition and increase biodiversity in model ecosystems. This contradiction, however, mainly stems from the incorrect application of formulas derived for fully connected networks to empirical, sparse networks.
James et al. (2012) presented simulations that apparently falsify the analytical result by Bastolla et al. (2009), who showed that nested mutualistic interactions decrease interspecific competition and increase biodiversity in model ecosystems. This contradiction, however, mainly stems from the incorrect application of formulas derived for fully connected networks to empirical, sparse networks.
△ Less
Submitted 16 January, 2013;
originally announced January 2013.
-
Phase transition in the Countdown problem
Authors:
Lucas Lacasa,
Bartolo Luque
Abstract:
Here we present a combinatorial decision problem, inspired by the celebrated quiz show called the countdown, that involves the computation of a given target number T from a set of k randomly chosen integers along with a set of arithmetic operations. We find that the probability of winning the game evidences a threshold phenomenon that can be understood in the terms of an algorithmic phase transiti…
▽ More
Here we present a combinatorial decision problem, inspired by the celebrated quiz show called the countdown, that involves the computation of a given target number T from a set of k randomly chosen integers along with a set of arithmetic operations. We find that the probability of winning the game evidences a threshold phenomenon that can be understood in the terms of an algorithmic phase transition as a function of the set size k. Numerical simulations show that such probability sharply transitions from zero to one at some critical value of the control parameter, hence separating the algorithm's parameter space in different phases. We also find that the system is maximally efficient close to the critical point. We then derive analytical expressions that match the numerical results for finite size and permit us to extrapolate the behavior in the thermodynamic limit.
△ Less
Submitted 13 June, 2012;
originally announced June 2012.
-
Feigenbaum graphs at the onset of chaos
Authors:
Bartolo Luque,
Lucas Lacasa,
Alberto Robledo
Abstract:
We analyze the properties of the self-similar network obtained from the trajectories of unimodal maps at the transition to chaos via the horizontal visibility (HV) algorithm. We first show that this network is uniquely determined by the encoded sequence of positions in the dynamics within the Feigenbaum attractor and it is universal in that it is independent of the shape and nonlinearity of the ma…
▽ More
We analyze the properties of the self-similar network obtained from the trajectories of unimodal maps at the transition to chaos via the horizontal visibility (HV) algorithm. We first show that this network is uniquely determined by the encoded sequence of positions in the dynamics within the Feigenbaum attractor and it is universal in that it is independent of the shape and nonlinearity of the maps in this class. We then find that the network degrees fluctuate at all scales with an amplitude that increases as the size of the network grows. This suggests the definition of a graph-theoretical Lyapunov exponent that measures the expansion rate of trajectories in network space. On good agreement with the map's counterpart, while at the onset of chaos this exponent vanishes, the subexponential expansion and contraction of network degrees can be fully described via a Tsallis-type scalar deformation of the expansion rate, that yields a discrete spectrum of non-null generalized exponents. We further explore the possibility of defining an entropy growth rate that describes the amount of information created along the trajectories in network space. Making use of the trajectory distributions in the map's accumulation point and the scaling properties of the associated network, we show that such entropic growth rate coincides with the spectrum of graph-theoretical exponents, what appears as a set of Pesin-like identities in the network.
△ Less
Submitted 31 October, 2012; v1 submitted 9 May, 2012;
originally announced May 2012.
-
Quasiperiodic graphs: structural design, scaling and entropic properties
Authors:
Bartolo Luque,
Fernando J. Ballesteros,
Ángel M. Núñez,
Alberto Robledo
Abstract:
A novel class of graphs, here named quasiperiodic, are constructed via application of the Horizontal Visibility algorithm to the time series generated along the quasiperiodic route to chaos. We show how the hierarchy of mode-locked regions represented by the Farey tree is inherited by their associated graphs. We are able to establish, via Renormalization Group (RG) theory, the architecture of the…
▽ More
A novel class of graphs, here named quasiperiodic, are constructed via application of the Horizontal Visibility algorithm to the time series generated along the quasiperiodic route to chaos. We show how the hierarchy of mode-locked regions represented by the Farey tree is inherited by their associated graphs. We are able to establish, via Renormalization Group (RG) theory, the architecture of the quasiperiodic graphs produced by irrational winding numbers with pure periodic continued fraction. And finally, we demonstrate that the RG fixed-point degree distributions are recovered via optimization of a suitably defined graph entropy.
△ Less
Submitted 16 March, 2012;
originally announced March 2012.
-
Analytical properties of horizontal visibility graphs in the Feigenbaum scenario
Authors:
Bartolo Luque,
Lucas Lacasa,
Fernando J. Ballesteros,
Alberto Robledo
Abstract:
Time series are proficiently converted into graphs via the horizontal visibility (HV) algorithm, which prompts interest in its capability for capturing the nature of different classes of series in a network context. We have recently shown [1] that dynamical systems can be studied from a novel perspective via the use of this method. Specifically, the period-doubling and band-splitting attractor cas…
▽ More
Time series are proficiently converted into graphs via the horizontal visibility (HV) algorithm, which prompts interest in its capability for capturing the nature of different classes of series in a network context. We have recently shown [1] that dynamical systems can be studied from a novel perspective via the use of this method. Specifically, the period-doubling and band-splitting attractor cascades that characterize unimodal maps transform into families of graphs that turn out to be independent of map nonlinearity or other particulars. Here we provide an in depth description of the HV treatment of the Feigenbaum scenario, together with analytical derivations that relate to the degree distributions, mean distances, clustering coefficients, etc., associated to the bifurcation cascades and their accumulation points. We describe how the resultant families of graphs can be framed into a renormalization group scheme in which fixed-point graphs reveal their scaling properties. These fixed points are then re-derived from an entropy optimization process defined for the graph sets, confirming a suggested connection between renormalization group and entropy optimization. Finally, we provide analytical and numerical results for the graph entropy and show that it emulates the Lyapunov exponent of the map independently of its sign.
△ Less
Submitted 12 January, 2012;
originally announced January 2012.
-
Feigenbaum graphs: a complex network perspective of chaos
Authors:
Bartolo Luque,
Lucas Lacasa,
Fernando J. Ballesteros,
Alberto Robledo
Abstract:
The recently formulated theory of horizontal visibility graphs transforms time series into graphs and allows the possibility of studying dynamical systems through the characterization of their associated networks. This method leads to a natural graph-theoretical description of nonlinear systems with qualities in the spirit of symbolic dynamics. We support our claim via the case study of the period…
▽ More
The recently formulated theory of horizontal visibility graphs transforms time series into graphs and allows the possibility of studying dynamical systems through the characterization of their associated networks. This method leads to a natural graph-theoretical description of nonlinear systems with qualities in the spirit of symbolic dynamics. We support our claim via the case study of the period-doubling and band-splitting attractor cascades that characterize unimodal maps. We provide a universal analytical description of this classic scenario in terms of the horizontal visibility graphs associated with the dynamics within the attractors, that we call Feigenbaum graphs, independent of map nonlinearity or other particulars. We derive exact results for their degree distribution and related quantities, recast them in the context of the renormalization group and find that its fixed points coincide with those of network entropy optimization. Furthermore, we show that the network entropy mimics the Lyapunov exponent of the map independently of its sign, hinting at a Pesin-like relation equally valid out of chaos.
△ Less
Submitted 6 September, 2011;
originally announced September 2011.
-
Detecting series periodicity with horizontal visibility graphs
Authors:
Angel M. Núñez,
Lucas Lacasa,
Eusebio Valero,
Jose Patricio Gómez,
Bartolo Luque
Abstract:
The horizontal visibility algorithm has been recently introduced as a mapping between time series and networks. The challenge lies in characterizing the structure of time series (and the processes that generated those series) using the powerful tools of graph theory. Recent works have shown that the visibility graphs inherit several degrees of correlations from their associated series, and therefo…
▽ More
The horizontal visibility algorithm has been recently introduced as a mapping between time series and networks. The challenge lies in characterizing the structure of time series (and the processes that generated those series) using the powerful tools of graph theory. Recent works have shown that the visibility graphs inherit several degrees of correlations from their associated series, and therefore such graph theoretical characterization is in principle possible. However, both the mathematical grounding of this promising theory and its applications are on its infancy. Following this line, here we address the question of detecting hidden periodicity in series polluted with a certain amount of noise. We first put forward some generic properties of horizontal visibility graphs which allow us to define a (graph theoretical) noise reduction filter. Accordingly, we evaluate its performance for the task of calculating the period of noisy periodic signals, and compare our results with standard time domain (autocorrelation) methods. Finally, potentials, limitations and applications are discussed.
△ Less
Submitted 8 August, 2011;
originally announced August 2011.
-
Time series irreversibility: a visibility graph approach
Authors:
Lucas Lacasa,
Ángel M. Núñez,
Édgar Roldán,
Juan M. R. Parrondo,
Bartolo Luque
Abstract:
We propose a method to measure real-valued time series irreversibility which combines two differ- ent tools: the horizontal visibility algorithm and the Kullback-Leibler divergence. This method maps a time series to a directed network according to a geometric criterion. The degree of irreversibility of the series is then estimated by the Kullback-Leibler divergence (i.e. the distinguishability) be…
▽ More
We propose a method to measure real-valued time series irreversibility which combines two differ- ent tools: the horizontal visibility algorithm and the Kullback-Leibler divergence. This method maps a time series to a directed network according to a geometric criterion. The degree of irreversibility of the series is then estimated by the Kullback-Leibler divergence (i.e. the distinguishability) between the in and out degree distributions of the associated graph. The method is computationally effi- cient, does not require any ad hoc symbolization process, and naturally takes into account multiple scales. We find that the method correctly distinguishes between reversible and irreversible station- ary time series, including analytical and numerical studies of its performance for: (i) reversible stochastic processes (uncorrelated and Gaussian linearly correlated), (ii) irreversible stochastic pro- cesses (a discrete flashing ratchet in an asymmetric potential), (iii) reversible (conservative) and irreversible (dissipative) chaotic maps, and (iv) dissipative chaotic maps in the presence of noise. Two alternative graph functionals, the degree and the degree-degree distributions, can be used as the Kullback-Leibler divergence argument. The former is simpler and more intuitive and can be used as a benchmark, but in the case of an irreversible process with null net current, the degree-degree distribution has to be considered to identifiy the irreversible nature of the series.
△ Less
Submitted 8 August, 2011;
originally announced August 2011.
-
Horizontal visibility graphs: exact results for random time series
Authors:
Bartolo Luque,
Lucas Lacasa,
Fernando Ballesteros,
Jordi Luque
Abstract:
The visibility algorithm has been recently introduced as a mapping between time series and complex networks. This procedure allows to apply methods of complex network theory for characterizing time series. In this work we present the horizontal visibility algorithm, a geometrically simpler and analytically solvable version of our former algorithm, focusing on the mapping of random series (series…
▽ More
The visibility algorithm has been recently introduced as a mapping between time series and complex networks. This procedure allows to apply methods of complex network theory for characterizing time series. In this work we present the horizontal visibility algorithm, a geometrically simpler and analytically solvable version of our former algorithm, focusing on the mapping of random series (series of independent identically distributed random variables). After presenting some properties of the algorithm, we present exact results on the topological properties of graphs associated to random series, namely the degree distribution, clustering coefficient, and mean path length. We show that the horizontal visibility algorithm stands as a simple method to discriminate randomness in time series, since any random series maps to a graph with an exponential degree distribution of the shape P(k) = (1/3)(2/3)**(k-2), independently of the probability distribution from which the series was generated. Accordingly, visibility graphs with other P(k) are related to non-random series. Numerical simulations confirm the accuracy of the theorems for finite series. In a second part, we show that the method is able to distinguish chaotic series from i.i.d. theory, studying the following situations: (i) noise-free low-dimensional chaotic series, (ii) low-dimensional noisy chaotic series, even in the presence of large amounts of noise, and (iii) high-dimensional chaotic series (coupled map lattice), without needs for additional techniques such as surrogate data or noise reduction methods. Finally, heuristic arguments are given to explain the topological properties of chaotic series and several sequences which are conjectured to be random are analyzed.
△ Less
Submitted 24 February, 2010;
originally announced February 2010.
-
The Visibility Graph: a new method for estimating the Hurst exponent of fractional Brownian motion
Authors:
Lucas Lacasa,
Bartolo Luque,
Jordi Luque,
Juan Carlos Nuno
Abstract:
Fractional Brownian motion (fBm) has been used as a theoretical framework to study real time series appearing in diverse scientific fields. Because its intrinsic non-stationarity and long range dependence, its characterization via the Hurst parameter H requires sophisticated techniques that often yield ambiguous results. In this work we show that fBm series map into a scale free visibility graph…
▽ More
Fractional Brownian motion (fBm) has been used as a theoretical framework to study real time series appearing in diverse scientific fields. Because its intrinsic non-stationarity and long range dependence, its characterization via the Hurst parameter H requires sophisticated techniques that often yield ambiguous results. In this work we show that fBm series map into a scale free visibility graph whose degree distribution is a function of H. Concretely, it is shown that the exponent of the power law degree distribution depends linearly on H. This also applies to fractional Gaussian noises (fGn) and generic f^(-b) noises. Taking advantage of these facts, we propose a brand new methodology to quantify long range dependence in these series. Its reliability is confirmed with extensive numerical simulations and analytical developments. Finally, we illustrate this method quantifying the persistent behavior of human gait dynamics.
△ Less
Submitted 7 January, 2009;
originally announced January 2009.
-
The first digit frequencies of primes and Riemann zeta zeros tend to uniformity following a size-dependent generalized Benford's law
Authors:
Bartolo Luque,
Lucas Lacasa
Abstract:
Prime numbers seem to distribute among the natural numbers with no other law than that of chance, however its global distribution presents a quite remarkable smoothness. Such interplay between randomness and regularity has motivated sci- entists of all ages to search for local and global patterns in this distribution that eventually could shed light into the ultimate nature of primes. In this wo…
▽ More
Prime numbers seem to distribute among the natural numbers with no other law than that of chance, however its global distribution presents a quite remarkable smoothness. Such interplay between randomness and regularity has motivated sci- entists of all ages to search for local and global patterns in this distribution that eventually could shed light into the ultimate nature of primes. In this work we show that a generalization of the well known first-digit Benford's law, which addresses the rate of appearance of a given leading digit d in data sets, describes with astonishing precision the statistical distribution of leading digits in the prime numbers sequence. Moreover, a reciprocal version of this pattern also takes place in the sequence of the nontrivial Riemann zeta zeros. We prove that the prime number theorem is, in the last analysis, the responsible of these patterns. Some new relations concerning the prime numbers distribution are also deduced, including a new approximation to the counting function pi(n). Furthermore, some relations concerning the statistical conformance to this generalized Benford's law are derived. Some applications are finally discussed.
△ Less
Submitted 20 November, 2008;
originally announced November 2008.
-
From time series to complex networks: the visibility graph
Authors:
Lucas Lacasa,
Bartolo Luque,
Fernando Ballesteros,
Jordi Luque,
Juan Carlos Nuno
Abstract:
In this work we present a simple and fast computational method, the visibility algorithm, that converts a time series into a graph. The constructed graph inherits several properties of the series in its structure. Thereby, periodic series convert into regular graphs, and random series do so into random graphs. Moreover, fractal series convert into scale-free networks, enhancing the fact that pow…
▽ More
In this work we present a simple and fast computational method, the visibility algorithm, that converts a time series into a graph. The constructed graph inherits several properties of the series in its structure. Thereby, periodic series convert into regular graphs, and random series do so into random graphs. Moreover, fractal series convert into scale-free networks, enhancing the fact that power law degree distributions are related to fractality, something highly discussed recently. Some remarkable examples and analytical tools are outlined in order to test the method's reliability. Many different measures, recently developed in the complex network theory, could by means of this new approach characterize time series from a new point of view.
△ Less
Submitted 6 October, 2008;
originally announced October 2008.
-
Number theoretic example of scale-free topology inducing self-organized criticality
Authors:
Bartolo Luque,
Octavio Miramontes,
Lucas Lacasa
Abstract:
In this work we present a general mechanism by which simple dynamics running on networks become self-organized critical for scale free topologies. We illustrate this mechanism with a simple arithmetic model of division between integers, the division model. This is the simplest self-organized critical model advanced so far, and in this sense it may help to elucidate the mechanism of self-organiza…
▽ More
In this work we present a general mechanism by which simple dynamics running on networks become self-organized critical for scale free topologies. We illustrate this mechanism with a simple arithmetic model of division between integers, the division model. This is the simplest self-organized critical model advanced so far, and in this sense it may help to elucidate the mechanism of self-organization to criticality. Its simplicity allows analytical tractability, characterizing several scaling relations. Furthermore, its mathematical nature brings about interesting connections between statistical physics and number theoretical concepts. We show how this model can be understood as a self-organized stochastic process embedded on a network, where the onset of criticality is induced by the topology.
△ Less
Submitted 17 September, 2008; v1 submitted 12 September, 2008;
originally announced September 2008.
-
Phase transition and computational complexity in a stochastic prime number generator
Authors:
Lucas Lacasa,
Bartolo Luque,
Octavio Miramontes
Abstract:
We introduce a prime number generator in the form of a stochastic algorithm. The character of such algorithm gives rise to a continuous phase transition which distinguishes a phase where the algorithm is able to reduce the whole system of numbers into primes and a phase where the system reaches a frozen state with low prime density. In this paper we firstly pretend to give a broad characterizati…
▽ More
We introduce a prime number generator in the form of a stochastic algorithm. The character of such algorithm gives rise to a continuous phase transition which distinguishes a phase where the algorithm is able to reduce the whole system of numbers into primes and a phase where the system reaches a frozen state with low prime density. In this paper we firstly pretend to give a broad characterization of this phase transition, both in terms of analytical and numerical analysis. Critical exponents are calculated, and data collapse is provided. Further on we redefine the model as a search problem, fitting it in the hallmark of computational complexity theory. We suggest that the system belongs to the class NP. The computational cost is maximal around the threshold, as common in many algorithmic phase transitions, revealing the presence of an easy-hard-easy pattern. We finally relate the nature of the phase transition to an average-case classification of the problem.
△ Less
Submitted 19 December, 2007;
originally announced December 2007.
-
Phase transition in a stochastic prime number generator
Authors:
Bartolo Luque,
Lucas Lacasa,
Octavio Miramontes
Abstract:
We introduce a stochastic algorithm that acts as a prime number generator. The dynamics of such algorithm gives rise to a continuous phase transition which separates a phase where the algorithm is able to reduce a whole set of integers into primes and a phase where the system reaches a frozen state with low prime density. We present both numerical simulations and an analytical approach in terms…
▽ More
We introduce a stochastic algorithm that acts as a prime number generator. The dynamics of such algorithm gives rise to a continuous phase transition which separates a phase where the algorithm is able to reduce a whole set of integers into primes and a phase where the system reaches a frozen state with low prime density. We present both numerical simulations and an analytical approach in terms of an annealed approximation, by means of which the data are collapsed. A critical slowing down phenomenon is also outlined.
△ Less
Submitted 13 July, 2007;
originally announced July 2007.
-
Self-overlap as a method of analysis in Ising models
Authors:
A. Ferrera,
B. Luque,
L. Lacasa,
E. Valero
Abstract:
The damage spreading method (DS) provided a useful tool to obtain analytical results of the thermodynamics and stability of the 2D Ising model --amongst many others--, but it suffered both from ambiguities in its results and from large computational costs. In this paper we propose an alternative method, the so called self-overlap method, based on the study of correlation functions measured at su…
▽ More
The damage spreading method (DS) provided a useful tool to obtain analytical results of the thermodynamics and stability of the 2D Ising model --amongst many others--, but it suffered both from ambiguities in its results and from large computational costs. In this paper we propose an alternative method, the so called self-overlap method, based on the study of correlation functions measured at subsequent time steps as the system evolves towards its equilibrium. Applying markovian and mean field approximations to a 2D Ising system we obtain both analytical and numerical results on the thermodynamics that agree with the expected behavior. We also provide some analytical results on the stability of the system. Since only a single replica of the system needs to be studied, this method would seem to be free from the ambiguities that afflicted DS. It also seems to be numerically more efficient and analytically simpler.
△ Less
Submitted 14 May, 2007;
originally announced May 2007.
-
Bonabeau hierarchy models revisited
Authors:
Lucas Lacasa,
Bartolo Luque
Abstract:
What basic processes generate hierarchy in a collective? The Bonabeau model provides us a simple mechanism based on randomness which develops self-organization through both winner/looser effects and relaxation process. A phase transition between egalitarian and hierarchic states has been found both analytically and numerically in previous works. In this paper we present a different approach: by…
▽ More
What basic processes generate hierarchy in a collective? The Bonabeau model provides us a simple mechanism based on randomness which develops self-organization through both winner/looser effects and relaxation process. A phase transition between egalitarian and hierarchic states has been found both analytically and numerically in previous works. In this paper we present a different approach: by means of a discrete scheme we develop a mean field approximation that not only reproduces the phase transition but also allows us to characterize the complexity of hierarchic phase. In the same philosophy, we study a new version of the Bonabeau model, developed by Stauffer et al. Several previous works described numerically the presence of a similar phase transition in this later version. We find surprising results in this model that can be interpreted properly as the non-existence of phase transition in this version of Bonabeau model, but a changing in fixed point structure.
△ Less
Submitted 11 November, 2005;
originally announced November 2005.
-
Small-Worlds, Mazes and Random Walks
Authors:
Bartolo Luque,
Miramontes Octavio
Abstract:
We establish a relationship between the Small-World behavior found in complex networks and a family of Random Walks trajectories using, as a linking bridge, a maze iconography. Simple methods to generate mazes using Random Walks are discussed along with related issues and it is explained how to interpret mazes as graphs and loops as shortcuts. Small-World behavior was found to be non-logarithmic…
▽ More
We establish a relationship between the Small-World behavior found in complex networks and a family of Random Walks trajectories using, as a linking bridge, a maze iconography. Simple methods to generate mazes using Random Walks are discussed along with related issues and it is explained how to interpret mazes as graphs and loops as shortcuts. Small-World behavior was found to be non-logarithmic but power-law in this model, we discuss the reason for this peculiar scaling
△ Less
Submitted 18 November, 2002;
originally announced November 2002.
-
Self-Organized Critical Random Boolean Networks
Authors:
Bartolo Luque,
Fernando J. Ballesteros,
Enrique M. Muro
Abstract:
Standard Random Boolean Networks display an order-disorder phase transition. We add to the standard Random Boolean Networks a disconnection rule which couples the control and order parameters. By this way, the system is driven to the critical line transition. Under the influence of perturbations the system point out self-organized critical behavior. Several numerical simulations have been done a…
▽ More
Standard Random Boolean Networks display an order-disorder phase transition. We add to the standard Random Boolean Networks a disconnection rule which couples the control and order parameters. By this way, the system is driven to the critical line transition. Under the influence of perturbations the system point out self-organized critical behavior. Several numerical simulations have been done and compared with a proposed analytical treatment.
△ Less
Submitted 22 February, 2001;
originally announced February 2001.
-
Small-world behavior in a system of mobile elements
Authors:
Susanna C. Manrubia,
Jordi Delgado,
Bartolo Luque
Abstract:
We analyze the propagation of activity in a system of mobile automata. A number r L^d of elements move as random walkers on a lattice of dimension d, while with a small probability p they can jump to any empty site in the system. We show that this system behaves as a Dynamic Small-World (DSW) and present analytic and numerical results for several quantities. Our analysis shows that the persisten…
▽ More
We analyze the propagation of activity in a system of mobile automata. A number r L^d of elements move as random walkers on a lattice of dimension d, while with a small probability p they can jump to any empty site in the system. We show that this system behaves as a Dynamic Small-World (DSW) and present analytic and numerical results for several quantities. Our analysis shows that the persistence time T* (equivalent to the persistence size L* of small-world networks) scales as T* ~ (r p)^(-t), with t = 1/(d+1).
△ Less
Submitted 5 February, 2001;
originally announced February 2001.
-
Measuring Mutual Information in Random Boolean Networks
Authors:
Bartolo Luque,
Antonio Ferrera
Abstract:
During the last few years an area of active research in the field of complex systems is that of their information storing and processing abilities. Common opinion has it that the most interesting beaviour of these systems is found ``at the edge of chaos'', which would seem to suggest that complex systems may have inherently non-trivial information proccesing abilities in the vicinity of sharp ph…
▽ More
During the last few years an area of active research in the field of complex systems is that of their information storing and processing abilities. Common opinion has it that the most interesting beaviour of these systems is found ``at the edge of chaos'', which would seem to suggest that complex systems may have inherently non-trivial information proccesing abilities in the vicinity of sharp phase transitions. A comprenhensive, quantitative understanding of why this is the case is however still lacking. Indeed, even ``experimental'' (i.e., often numerical) evidence that this is so has been questioned for a number of systems. In this paper we will investigate, both numerically and analitically, the behavior of Random Boolean Networks (RBN's) as they undergo their order-disorder phase transition. We will use a simple mean field approximation to treat the problem, and without lack of generality we will concentrate on a particular value for the connectivity of the system. In spite of the simplicity of our arguments, we will be able to reproduce analitically the amount of mutual information contained in the system as measured from numerical simulations.
△ Less
Submitted 3 September, 1999;
originally announced September 1999.
-
Statistical measures of complexity for strongly interacting systems
Authors:
Ricard V. Sole,
Bartolo Luque
Abstract:
In recent studies, new measures of complexity for nonlinear systems have been proposed based on probabilistic grounds, as the LMC measure (Phys. Lett. A {\bf 209} (1995) 321) or the SDL measure (Phys. Rev. E {\bf 59}
(1999) 2). All these measures share an intuitive consideration: complexity seems to emerge in nature close to instability points, as for example the phase transition points charact…
▽ More
In recent studies, new measures of complexity for nonlinear systems have been proposed based on probabilistic grounds, as the LMC measure (Phys. Lett. A {\bf 209} (1995) 321) or the SDL measure (Phys. Rev. E {\bf 59}
(1999) 2). All these measures share an intuitive consideration: complexity seems to emerge in nature close to instability points, as for example the phase transition points characteristic of critical phenomena. Here we discuss these measures and their reliability for detecting complexity close to critical points in complex systems composed of many interacting units. Both a two-dimensional spatially extended problem (the 2D Ising model) and a $\infty$-dimensional (random graph) model (random Boolean networks) are analysed. It is shown that the LMC and the SDL measures can be easily generalized to extended systems but fails to detect real complexity.
△ Less
Submitted 27 August, 1999;
originally announced September 1999.
-
Phase Transition in Random Networks with Multiple States
Authors:
Ricard V. Sole,
Bartolo Luque,
Stuart Kauffman
Abstract:
The critical boundaries separating ordered from chaotic behavior in randomly wired S-state networks are calculated. These networks are a natural generalization of random Boolean nets and are proposed as on extended approach to genetic regulatory systems, sets of cells in different states or collectives of agents engaged into a set of S possible tasks. A order parameter for the transition is comp…
▽ More
The critical boundaries separating ordered from chaotic behavior in randomly wired S-state networks are calculated. These networks are a natural generalization of random Boolean nets and are proposed as on extended approach to genetic regulatory systems, sets of cells in different states or collectives of agents engaged into a set of S possible tasks. A order parameter for the transition is computed and analysed. The relevance of these networks to biology, their relationships with standard cellular automata and possible extensions are outlined.
△ Less
Submitted 23 July, 1999;
originally announced July 1999.
-
Lyapunov Exponents in Random Boolean Networks
Authors:
Bartolo Luque,
Ricard V. Sole
Abstract:
A new order parameter approximation to Random Boolean Networks (RBN) is introduced, based on the concept of Boolean derivative. A statistical argument involving an annealed approximation is used, allowing to measure the order parameter in terms of the statistical properties of a random matrix. Using the same formalism, a Lyapunov exponent is calculated, allowing to provide the onset of damage sp…
▽ More
A new order parameter approximation to Random Boolean Networks (RBN) is introduced, based on the concept of Boolean derivative. A statistical argument involving an annealed approximation is used, allowing to measure the order parameter in terms of the statistical properties of a random matrix. Using the same formalism, a Lyapunov exponent is calculated, allowing to provide the onset of damage spreading through the network and how sensitive it is to minimal perturbations. Finally, the Lyapunov exponents are obtained by means of different approximations: through distance method and a discrete variant of the Wolf's method for continuous systems.
△ Less
Submitted 30 June, 1999;
originally announced July 1999.