Arrow of time across five centuries of classical music

The concept of time series irreversibility—the degree by which the statistics of signals are not invariant under time reversal—naturally appears in nonequilibrium physics in stationary systems which operate away from equilibrium and produce entropy. This concept has not been explored to date in the realm of musical scores as these are typically short sequences whose time reversibility estimation could suffer from strong finite size effects which preclude interpretability. Here we show that the so-called horizontal visibility graph method—which recently was shown to quantify such statistical property even in nonstationary signals—is a method that can estimate time reversibility of short symbolic sequences, thus unlocking the possibility of exploring such properties in the context of musical compositions. Accordingly, we analyze over 8000 musical pieces ranging from the Renaissance to the early Modern period and show that, indeed, most of them display clear signatures of time irreversibility. Since by construction stochastic processes with a linear correlation structure (such as 1 / f noise) are time reversible, we conclude that musical compositions have a considerably richer structure, that goes beyond the traditional properties retrieved by the power spectrum or similar approaches. We also show that musical compositions display strong signs of nonlinear correlations, that nonlinearity is correlated to irreversibility, and that these are also related to asymmetries in the abundance of musical intervals, which we associate to the narrative underpinning a musical composition. These findings provide tools for the study of musical periods and composers, as well as criteria related to music appreciation and cognition

Exact scaling in the expansion-modification system

This work is devoted to the study of the scaling, and the consequent power-law behavior, of the correlation function in a mutation-replication model known as the expansion-modification system. The latter is a biology inspired random substitution model for the genome evolution, which is defined on a binary alphabet and depends on a parameter interpreted as a \emph{mutation probability}. We prove that the time-evolution of this system is such that any initial measure converges towards a unique stationary one exhibiting decay of correlations not slower than a power-law. We then prove, for a significant range of mutation probabilities, that the decay of correlations indeed follows a power-law with scaling exponent smoothly depending on the mutation probability. Finally we put forward an argument which allows us to give a closed expression for the corresponding scaling exponent for all the values of the mutation probability. Such a scaling exponent turns out to be a piecewise smooth function of the parameter.Comment: 22 pages, 2 figure

Transport Properties of the Diluted Lorentz Slab

We study the behavior of a point particle incident from the left on a slab of a randomly diluted triangular array of circular scatterers. Various scattering properties, such as the reflection and transmission probabilities and the scattering time are studied as a function of thickness and dilution. We show that a diffusion model satisfactorily describes the mentioned scattering properties. We also show how some of these quantities can be evaluated exactly and their agreement with numerical experiments. Our results exhibit the dependence of these scattering data on the mean free path. This dependence again shows excellent agreement with the predictions of a Brownian motion model.Comment: 14 pages of text in LaTeX, 7 figures in Postscrip

Evolution of scaling emergence in large-scale spatial epidemic spreading

Background: Zipf's law and Heaps' law are two representatives of the scaling concepts, which play a significant role in the study of complexity science. The coexistence of the Zipf's law and the Heaps' law motivates different understandings on the dependence between these two scalings, which is still hardly been clarified. Methodology/Principal Findings: In this article, we observe an evolution process of the scalings: the Zipf's law and the Heaps' law are naturally shaped to coexist at the initial time, while the crossover comes with the emergence of their inconsistency at the larger time before reaching a stable state, where the Heaps' law still exists with the disappearance of strict Zipf's law. Such findings are illustrated with a scenario of large-scale spatial epidemic spreading, and the empirical results of pandemic disease support a universal analysis of the relation between the two laws regardless of the biological details of disease. Employing the United States(U.S.) domestic air transportation and demographic data to construct a metapopulation model for simulating the pandemic spread at the U.S. country level, we uncover that the broad heterogeneity of the infrastructure plays a key role in the evolution of scaling emergence. Conclusions/Significance: The analyses of large-scale spatial epidemic spreading help understand the temporal evolution of scalings, indicating the coexistence of the Zipf's law and the Heaps' law depends on the collective dynamics of epidemic processes, and the heterogeneity of epidemic spread indicates the significance of performing targeted containment strategies at the early time of a pandemic disease.Comment: 24pages, 7figures, accepted by PLoS ON

Interaction of the IP3-Ca2+ and MAPK signaling systems in the Xenopus blastomere: a possible frequency encoding mechanism for the control of the Xbra gene expression

The intense periodic calcium activity experimentally observed in the Xenopus embryo at the Mid Blastula Transition stage is closely related to the competence of the embryonic cells of the marginal zone to respond to the posterior-mesodermal inducting signals from the Fibroblast Growth Factor (FGF). In this work we do a stability analysis and study numerically an extension of a mathematical model Previously introduced by us [Diaz, J., Baier, G., Martinez-Mekler, G., Pastor, N.. 2002. Interaction of the IP3-Ca2+ and the FGF-MAPK signaling pathways in the Xenopus laevis embryo: a qualitative approach to the mesodermal induction problem. Biophys. Chem. 97, 55-72] for the interaction of the Inositol 1,4,5-triphosphate-Calcium (IP3-Ca2+) and the Mitogen-Activated Protein Kinase (MAPK) signaling pathways at the Mid Blastula Transition stage or stage 8 of development. This allows us to consider the effect of the oscillatory calcium dynamics on the FGF input signal carried by the MAP kinase (ERK) into the nucleus. We find that this interaction of the pathways induces a limit cycle behavior for ERK with frequency-encoding characteristics. We believe that this periodic increase of the ERK levels in the nucleus is related to the ability of the cell to express posteriorizing mesodermal features induced by the FGF signal at stage 8. (c) 2004 Society for Mathematical Biology. Published by Elsevier Ltd. All rights reserved

Power spectrum crossover in sediments of a paleolake disturbed by volcanism

We study density fluctuations from sediments of a paleolake in central Mexico that was subjected to volcanic perturbations by means of computed tomography (CT) measurements on blocks chiselled out of mines at the lake's bed. The mine walls show laminations corresponding to the alternation of low density diatom sediments and high density volcanic ash depositions. We have previously shown that there is a range of scales where these fluctuations present a self-similar behavior [1]. Here we relate density correlation calculations to the power spectrum of the fluctuations. We show that a scaling region in the power spectrum coincides with the scaling region in the correlations produced by relaxation from intense volcanic perturbations to steady state fluctuations. There appears to be a kink-like crossover in the power spectrum from mid range scaling to a shorter range scale invariance. This, together with the density probability distribution of the fluctuations, draws attention to the dominant role of rare events. We believe that our analysis may be useful for the understanding of other phenomena with similar power spectrum properties, in which a scale invariance in the unperturbed system is altered by external perturbations that induce an additional scaling behavior

Role of a spatial distribution of IP3 receptors in the Ca2+ dynamics of the Xenopus embryo at the mid-blastula transition stage

Periodic calcium activity correlates temporally with the onset of gene expression in the embryo, suggesting a causal relation between these two events. Calcium transients are elicited by the action of fibroblast growth factor (FGF) through the activation of phospholipase C. In this work, we present a reaction-diffusion model that extends our previous results on the generation of calcium oscillations for a single and two coupled blastomere cells to a meridian of the Xenopus embryo at the mid-blastula transition. In the model, all cells are subject to the same amount of FGF and contain the same concentration of intracellular components, except for the amount of IP3 receptors (IP3R). A bell-shaped distribution of IP3R produces the correct shape of the calcium transients experimentally observed in the Xenopus blastula at stage 8 (mid-blastula transition stage). The model is also capable of predicting period and amplitude values close to the experimental values. In our model, calcium transients induce spatially localized ERK periodic transients that could activate specific nuclear genes, allowing for the regional differentiation of the cells in the zone under the influence of the calcium signal. (C) 2004 Wiley-Liss, Inc

Deterministic ratchets, circle maps, and current reversals

In this work we transform the deterministic dynamics of an overdamped tilting ratchet into a discrete dynamical map by looking stroboscopically at the continuous motion originally ruled by differential equations. We show that, for the simple and widely used case of periodic dichotomous driving forces, the resulting discrete map belongs to the class of circle homeomorphisms. This approach allows us to apply the well-known properties of such maps to derive the necessary and sufficient conditions that the ratchet potential must satisfy in order to have a vanishing current. Furthermore, as a consequence of the above, we show (i) that there is a class of periodic potentials which do not exhibit the rectification phenomenon in spite of their asymmetry and (ii) that current reversals occur in the deterministic case for a large class of ratchet potentials