Source author record

Diego Pazó

Diego Pazó appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Enlarged Kuramoto Model: Secondary Instability and Transition to Collective Chaos

The emergence of collective synchrony from an incoherent state is a phenomenon essentially described by the Kuramoto model. This canonical model was derived perturbatively, by applying phase reduction to an ensemble of heterogeneous, globally coupled Stuart-Landau oscillators. This derivation neglects nonlinearities in the coupling constant. We show here that a comprehensive analysis requires extending the Kuramoto model up to quadratic order. This "enlarged Kuramoto model" comprises three-body (nonpairwise) interactions, which induce strikingly complex phenomenology at certain parameter values. As the coupling is increased, a secondary instability renders the synchronized state unstable, and subsequent bifurcations lead to collective chaos. An efficient numerical study of the thermodynamic limit, valid for Gaussian heterogeneity, is carried out by means of a Fourier-Hermite decomposition of the oscillator density.

preprint2021arXiv

Comment on "The Winfree model with non-infinitesimal phase-response curve: Ott-Antonsen theory" [Chaos 30, 073139 (2020)]

In a recent paper [Chaos 30, 073139 (2020)] we analyzed an extension of the Winfree model with nonlinear interactions. The nonlinear coupling function Q was mistakenly identified with the non-infinitesimal phase-response curve (PRC). Here, we asses to what extent Q and the actual PRC differ in practice. By means of numerical simulations, we compute the PRCs corresponding to the Q functions previously considered. The results confirm a qualitative similarity between the PRC and the coupling function Q in all cases.

preprint2020arXiv

The Winfree model with non-infinitesimal phase-response curve: Ott-Antonsen theory

A novel generalization of the Winfree model of globally coupled phase oscillators, representing phase reduction under finite coupling, is studied analytically. We consider interactions through a non-infinitesimal (or finite) phase-response curve (PRC), in contrast to the infinitesimal PRC of the original model. For a family of non-infinitesimal PRCs, the global dynamics is captured by one complex-valued ordinary differential equation resorting to the Ott-Antonsen ansatz. The phase diagrams are thereupon obtained for four illustrative cases of non-infinitesimal PRC. Bistability between collective synchronization and full desynchronization is observed in all cases.

preprint2016arXiv

From quasiperiodic partial synchronization to collective chaos in populations of inhibitory neurons with delay

Collective chaos is shown to emerge, via a period-doubling cascade, from quasiperiodic partial synchronization in a population of identical inhibitory neurons with delayed global coupling. This system is thoroughly investigated by means of an exact model of the macroscopic dynamics, valid in the thermodynamic limit. The collective chaotic state is reproduced numerically with a finite population, and persists in the presence of weak heterogeneities. Finally, the relationship of the model's dynamics with fast neuronal oscillations is discussed.

preprint2015arXiv

Macroscopic description for networks of spiking neurons

A major goal of neuroscience, statistical physics and nonlinear dynamics is to understand how brain function arises from the collective dynamics of networks of spiking neurons. This challenge has been chiefly addressed through large-scale numerical simulations. Alternatively, researchers have formulated mean-field theories to gain insight into macroscopic states of large neuronal networks in terms of the collective firing activity of the neurons, or the firing rate. However, these theories have not succeeded in establishing an exact correspondence between the firing rate of the network and the underlying microscopic state of the spiking neurons. This has largely constrained the range of applicability of such macroscopic descriptions, particularly when trying to describe neuronal synchronization. Here we provide the derivation of a set of exact macroscopic equations for a network of spiking neurons. Our results reveal that the spike generation mechanism of individual neurons introduces an effective coupling between two biophysically relevant macroscopic quantities, the firing rate and the mean membrane potential, which together govern the evolution of the neuronal network. The resulting equations exactly describe all possible macroscopic dynamical states of the network, including states of synchronous spiking activity. Finally we show that the firing rate description is related, via a conformal map, with a low-dimensional description in terms of the Kuramoto order parameter, called Ott-Antonsen theory. We anticipate our results will be an important tool in investigating how large networks of spiking neurons self-organize in time to process and encode information in the brain.

preprint2014arXiv

Low-dimensional dynamics of populations of pulse-coupled oscillators

Large communities of biological oscillators show a prevalent tendency to self-organize in time. This cooperative phenomenon inspired Winfree to formulate a mathematical model that originated the theory of macroscopic synchronization. Despite its fundamental importance, a complete mathematical analysis of the model proposed by Winfree ---consisting of a large population of all-to-all pulse-coupled oscillators--- is still missing. Here we show that the dynamics of the Winfree model evolves into the so-called Ott-Antonsen manifold. This important property allows for an exact description of this high-dimensional system in terms of a few macroscopic variables, and the full investigation of its dynamics. We find that brief pulses are capable of synchronizing heterogeneous ensembles which fail to synchronize with broad pulses, specially for certain phase response curves. Finally, to further illustrate the potential of our results, we investigate the possibility of `chimera' states in populations of identical pulse-coupled oscillators. Chimeras are self-organized states in which the symmetry of a population is broken into a synchronous and an asynchronous part. Here we derive three ordinary differential equations describing two coupled populations, and uncover a variety of chimera states, including a new class with chaotic dynamics.

preprint2014arXiv

Synchronizing spatio-temporal chaos with imperfect models: a stochastic surface growth picture

We study the synchronization of two spatially extended dynamical systems where the models have imperfections. We show that the synchronization error across space can be visualized as a rough surface governed by the Kardar-Parisi-Zhang equation with both upper and lower bounding walls corresponding to nonlinearities and model discrepancies, respectively. Two types of model imperfections are considered: parameter mismatch and unresolved fast scales, finding in both cases the same qualitative results. The consistency between different setups and systems indicates the results are generic for a wide family of spatially extended systems.

preprint2013arXiv

On the angle between the first and the second Lyapunov vectors in spatio-temporal chaos

In a dynamical system the first Lyapunov vector (LV) is associated with the largest Lyapunov exponent and indicates ---at some point on the attractor--- the direction of maximal growth in tangent space. The LV corresponding to the second largest Lyapunov exponent generally points at a different direction, but tangencies between both vectors can in principle occur. Here we find that the probability density function (PDF) of the angle ψspanned by the first and the second LVs should be expected approximately symmetric around π/4 and peaked at 0 and π/2. Moreover, for small angles we uncover a scaling law for the PDF Q of ψ_l=\lnψwith the system size L: Q(ψ_l)=L^{-1/2} f(ψ_l L^{-1/2}). We give a theoretical argument that justifies this scaling form and also explains why it should be universal (irrespective of the system details) for spatio-temporal chaos in one spatial dimension.

preprint2013arXiv

Universal scaling of Lyapunov-exponent fluctuations in space-time chaos

Finite-time Lyapunov exponents of generic chaotic dynamical systems fluctuate in time. These fluctuations are due to the different degree of stability across the accessible phase-space. A recent numerical study of spatially-extended systems has revealed that the diffusion coefficient D of the Lyapunov exponents (LEs) exhibits a non-trivial scaling behavior, D(L) ~ L^{-γ}, with the system size L. Here, we show that the wandering exponent γcan be expressed in terms of the roughening exponents associated with the corresponding "Lyapunov-surface". Our theoretical predictions are supported by the numerical analysis of several spatially-extended systems. In particular, we find that the wandering exponent of the first LE is universal: in view of the known relationship with the Kardar-Parisi-Zhang equation, γcan be expressed in terms of known critical exponents. Furthermore, our simulations reveal that the bulk of the spectrum exhibits a clearly different behavior and suggest that it belongs to a possibly unique universality class, which has, however, yet to be identified.

preprint2012arXiv

Covariant hydrodynamic Lyapunov modes and strong stochasticity threshold in Hamiltonian lattices

We scrutinize the reliability of covariant and Gram-Schmidt Lyapunov vectors for capturing hydrodynamic Lyapunov modes (HLMs) in one-dimensional Hamiltonian lattices. We show that,in contrast with previous claims, HLMs do exist for any energy density, so that strong chaos is not essential for the appearance of genuine (covariant) HLMs. In contrast, Gram-Schmidt Lyapunov vectors lead to misleading results concerning the existence of HLMs in the case of weak chaos.

preprint2011arXiv

Characteristic Lyapunov vectors in chaotic time-delayed systems

We compute Lyapunov vectors (LVs) corresponding to the largest Lyapunov exponents in delay-differential equations with large time delay. We find that characteristic LVs, and backward (Gram-Schmidt) LVs, exhibit long-range correlations, identical to those already observed in dissipative extended systems. In addition we give numerical and theoretical support to the hypothesis that the main LV belongs, under a suitable transformation, to the universality class of the Kardar-Parisi-Zhang equation. These facts indicate that in the large delay limit (an important class of) delayed equations behave exactly as dissipative systems with spatiotemporal chaos.

preprint2011arXiv

Collective synchronization in the presence of reactive coupling and shear diversity

We analyze the synchronization dynamics of a model obtained from the phase reduction of the mean-field complex Ginzburg-Landau equation with heterogeneity. We present exact results that uncover the role of dissipative and reactive couplings on the synchronization transition when shears and natural frequencies are independently distributed. As it occurs in the purely dissipative case, an excess of shear diversity prevents the onset of synchronization, but this does not hold true if coupling is purely reactive. In this case the synchronization threshold turns out to depend on the mean of the shear distribution, but not on all the other distribution's moments.

preprint2011arXiv

Maximizing the statistical diversity of an ensemble of bred vectors by using the geometric norm

We show that the choice of the norm has a great impact on the construction of ensembles of bred vectors. The geometric norm maximizes (in comparison with other norms like the Euclidean one) the statistical diversity of the ensemble while, at the same time, enhances the growth rate of the bred vector and its projection on the linearly most unstable direction, i.e. the Lyapunov vector. The geometric norm is also optimal in providing the least fluctuating ensemble dimension among all the spectrum of q-norms studied. We exemplify our results with numerical integrations of a toy model of the atmosphere (the Lorenz-96 model), but our findings are expected to be generic for spatially extended chaotic systems.

preprint2011arXiv

Shear diversity prevents collective synchronization

Large ensembles of heterogeneous oscillators often exhibit collective synchronization as a result of mutual interactions. If the oscillators have distributed natural frequencies and common shear (or nonisochronicity), the transition from incoherence to collective synchronization is known to occur at large enough values of the coupling strength. However, here we demonstrate that shear diversity cannot be counterbalanced by diffusive coupling leading to synchronization. We present the first analytical results for the Kuramoto model with distributed shear, and show that the onset of collective synchronization is impossible if the width of the shear distribution exceeds a precise threshold.

preprint2011arXiv

The Kuramoto model with distributed shear

We uncover a solvable generalization of the Kuramoto model in which shears (or nonisochronicities) and natural frequencies are distributed and statistically dependent. We show that the strength and sign of this dependence greatly alter synchronization and yield qualitatively different phase diagrams. The Ott-Antonsen ansatz allows us to obtain analytical results for a specific family of joint distributions. We also derive, using linear stability analysis, general formulae for the stability border of incoherence.

preprint2010arXiv

Logarithmic bred vectors in spatiotemporal chaos: structure and growth

Bred vectors are a type of finite perturbation used in prediction studies of atmospheric models that exhibit spatially extended chaos. We study the structure, spatial correlations, and the growth- rates of logarithmic bred vectors (which are constructed by using a given norm). We find that, after a suitable transformation, logarithmic bred vectors are roughly piecewise copies of the leading Lyapunov vector. This fact allows us to deduce a scaling law for the bred vector growth rate as a function of their amplitude. In addition, we relate growth rates with the spectrum of Lyapunov exponents corresponding to the most expanding directions. We illustrate our results with simulations of the Lorenz '96 model.

preprint2010arXiv

Structure of characteristic Lyapunov vectors in anharmonic Hamiltonian lattices

In this work we perform a detailed study of the scaling properties of Lyapunov vectors (LVs) for two different one-dimensional Hamiltonian lattices: the Fermi-Pasta-Ulam and $Φ^4$ models. In this case, characteristic (also called covariant) LVs exhibit qualitative similarities with those of dissipative lattices but the scaling exponents are different and seemingly nonuniversal. In contrast, backward LVs (obtained via Gram-Schmidt orthonormalizations) present approximately the same scaling exponent in all cases, suggesting it is an artificial exponent produced by the imposed orthogonality of these vectors. We are able to compute characteristic LVs in large systems thanks to a `bit reversible' algorithm, which completely obviates computer memory limitations.

preprint2009arXiv

Existence of hysteresis in the Kuramoto model with bimodal frequency distributions

We investigate the transition to synchronization in the Kuramoto model with bimodal distributions of the natural frequencies. Previous studies have concluded that the model exhibits a hysteretic phase transition if the bimodal distribution is close to a unimodal one, due to the shallowness the central dip. Here we show that proximity to the unimodal-bimodal border does not necessarily imply hysteresis when the width, but not the depth, of the central dip tends to zero. We draw this conclusion from a detailed study of the Kuramoto model with a suitable family of bimodal distributions.

preprint2009arXiv

Stability of an [N/2]-dimensional invariant torus in the Kuramoto model at small coupling

When the natural frequencies are allocated symmetrically in the Kuramoto model there exists an invariant torus of dimension [N/2]+1 (N is the population size). A global phase shift invariance allows to reduce the model to $N-1$ dimensions using the phase differences, and doing so the invariant torus becomes [N/2]-dimensional. By means of perturbative calculations based on the renormalization group technique, we show that this torus is asymptotically stable at small coupling if N is odd. If N is even the torus can be stable or unstable depending on the natural frequencies, and both possibilities persist in the small coupling limit.