Search

The Arecchi-Bonifacio (or Maxwell-Bloch) model is the benchmark for the description of active optical media. However, in the presence of a fast relaxation of the atomic polarization, its implementation is a challenging task even in the simple ring-laser configuration, due to the presence of multiple time scales. In this Article we show that the dynamics is nearly Hamiltonian over time scales much longer than those of the cavity losses. More precisely, we prove that it can be represented as a pseudo spatio-temporal pattern generated by a nonlinear wave equation equipped with a Toda potential. The existence of two constants of motion (identified as pseudo energies), thereby, elucidates the reason why it is so hard to simplify the original model: the adiabatic elimination of the polarization must be accurate enough to describe the dynamics correctly over unexpectedly long time scales. Finally, since the nonlinear wave equation with Toda potential can be simulated on much longer times than the previous models, this opens up the route to the numerical (and theoretical) investigation of realistic setups.

Rhythmic activities that alternate between coherent and incoherent phases are ubiquitous in chemical, ecological, climate, or neural systems. Despite their importance, general mechanisms for their emergence are little understood. In order to fill this gap, we present a framework for describing the emergence of recurrent synchronization in complex networks with adaptive interactions. This phenomenon is manifested at the macroscopic level by temporal episodes of coherent and incoherent dynamics that alternate recurrently. At the same time, the dynamics of the individual nodes do not change qualitatively. We identify asymmetric adaptation rules and temporal separation between the adaptation and the dynamics of individual nodes as key features for the emergence of recurrent synchronization. Our results suggest that asymmetric adaptation might be a fundamental ingredient for recurrent synchronization phenomena as seen in pattern generators, e.g., in neuronal systems.

We show that many delay-based reservoir computers considered in the literature can be characterized by a universal master memory function (MMF). Once computed for two independent parameters, this function provides linear memory capacity for any delay-based single-variable reservoir with small inputs. Moreover, we propose an analytical description of the MMF that enables its efficient and fast computation. Our approach can be applied not only to reservoirs governed by known dynamical rules such as Mackey-Glass or Ikeda-like systems but also to reservoirs whose dynamical model is not available. We also present results comparing the performance of the reservoir computer and the memory capacity given by the MMF., Comment: To be published

Real-world systems can be strongly influenced by time delays occurring in self-coupling interactions, due to unavoidable finite signal propagation velocities. When the delays become significantly long, complicated high-dimensional phenomena appear and a simple extension of the methods employed in low-dimensional dynamical systems is not feasible. We review the general theory developed in this case, describing the main destabilization mechanisms, the use of visualization tools, and commenting on the most important and effective dynamical indicators as well as their properties in different regimes. We show how a suitable approach, based on a comparison with spatio-temporal systems, represents a powerful instrument for disclosing the very basic mechanism of long-delay systems. Various examples from different models and a series of recent experiments are reported.

Phase transitions in equilibrium and nonequilibrium systems play a major role in the natural sciences. In dynamical networks, phase transitions organize qualitative changes in the collective behavior of coupled dynamical units. Adaptive dynamical networks feature a connectivity structure that changes over time, co-evolving with the nodes' dynamical state. In this Letter, we show the emergence of two distinct first-order nonequilibrium phase transitions in a finite size adaptive network of heterogeneous phase oscillators. Depending on the nature of defects in the internal frequency distribution, we observe either an abrupt single-step transition to full synchronization or a more gradual multi-step transition. This observation has a striking resemblance to heterogeneous nucleation. We develop a mean-field approach to study the interplay between adaptivity and nodal heterogeneity and describe the dynamics of multicluster states and their role in determining the character of the phase transition. Our work provides a theoretical framework for studying the interplay between adaptivity and nodal heterogeneity.

The acclaimed Maxwell-Bloch (or Arecchi-Bonifacio) equations are a valid dynamical model, effectively describing wave propagation in nonlinear optical media: from the amplification in input-output devices to multimode instabilities arising in laser systems. However, the inherent spatial variability of the physical observables represents an obstacle to fast simulations and analysis, especially whenever networks of active elements have to be considered. In this paper, we propose an approach which, stripping the spatial dependence of its role as a generator of dynamical richness, allows for a compelling simple portrait. It leads to (a few) ordinary differential equations in input-output configurations, complemented by a time-delayed feedback in closed-loop setups. Such scheme reproduces accurately the dynamics, paving the way to a plain treatment of the wealth of phenomena described by the Maxwell-Bloch equations., Comment: 3 figures

In this paper we give a profound insight into the computation capability of delay-based reservoir computing via an eigenvalue analysis. We concentrate on the task-independent memory capacity to quantify the reservoir performance and compare these with the eigenvalue spectrum of the dynamical system. We show that these two quantities are deeply connected, and thus the reservoir computing performance is predictable by analyzing the small signal response of the reservoir. Our results suggest that any dynamical system used as a reservoir can be analyzed in this way. We apply our method exemplarily to a photonic laser system with feedback and compare the numerically computed recall capabilities with the eigenvalue spectrum. Optimal performance is found for a system with the eigenvalues having real parts close to zero and off-resonant imaginary parts., New Journal Submission

An equilibrium of a delay differential equation (DDE) is absolutely stable, if it is locally asymptotically stable for all delays. We present criteria for absolute stability of DDEs with discrete time-delays. In the case of a single delay, the absolute stability is shown to be equivalent to asymptotic stability for sufficiently large delays. Similarly, for multiple delays, the absolute stability is equivalent to asymptotic stability for hierarchically large delays. Additionally, we give necessary and sufficient conditions for a linear DDE to be hyperbolic for all delays. The latter conditions are crucial for determining whether a system can have stabilizing or destabilizing bifurcations by varying time delays., 21 pages, 2 figures

First Published in SIAM Journal on Applied Dynamical Systems in Volume 20, Issue 4 (2021), published by the Society for Industrial and Applied Mathematics (SIAM).Copyright © by SIAM. Unauthorized reproduction of this article is prohibited. Multiplex networks are networks composed of multiple layers such that the number of nodes in all layers is the same and the adjacency matrices between the layers are diagonal. We consider the special class of multiplex networks where the adjacency matrices for each layer are simultaneously triagonalizable. For such networks, we derive the relation between the spectrum of the multiplex network and the eigenvalues of the individual layers. As an application, we propose a generalized master stability approach that allows for a simplified, low-dimensional description of the stability of synchronized solutions in multiplex networks. We illustrate our result with a duplex network of FitzHugh--Nagumo oscillators. In particular, we show how interlayer interaction can lead to stabilization or destabilization of the synchronous state. Finally, we give explicit conditions for the stability of synchronous solutions in duplex networks of linear diffusive systems.

Adaptive networks change their connectivity with time, depending on their dynamical state. While synchronization in structurally static networks has been studied extensively, this problem is much more challenging for adaptive networks. In this Letter, we develop the master stability approach for a large class of adaptive networks. This approach allows for reducing the synchronization problem for adaptive networks to a low-dimensional system, by decoupling topological and dynamical properties. We show how the interplay between adaptivity and network structure gives rise to the formation of stability islands. Moreover, we report a desynchronization transition and the emergence of complex partial synchronization patterns induced by an increasing overall coupling strength. We illustrate our findings using adaptive networks of coupled phase oscillators and FitzHugh-Nagumo neurons with synaptic plasticity.

Power grid networks, as well as neuronal networks with synaptic plasticity, describe real-world systems of tremendous importance for our daily life. The investigation of these seemingly unrelated types of dynamical networks has attracted increasing attention over the past decade. In this paper, we provide insight into the fundamental relation between these two types of networks. For this, we consider well-established models based on phase oscillators and show their intimate relation. In particular, we prove that phase oscillator models with inertia can be viewed as a particular class of adaptive networks. This relation holds even for more general classes of power grid models that include voltage dynamics. As an immediate consequence of this relation, we discover a plethora of multicluster states for phase oscillators with inertia. Moreover, the phenomenon of cascading line failure in power grids is translated into an adaptive neuronal network.

We present an index for the local sensitivity of spatiotemporal structures in coupled oscillatory systems based on the asymptotic scaling of local-in-space, finite-time Lyapunov Exponents. For a system of nonlocally-coupled R\"{o}ssler oscillators, we show that deviations of this index reflect the sensitivity to noise and the onset of spatial chaos for the patterns where coherence and incoherence regions coexist., Comment: 15 pages, 6 figures

We investigate the stability of synchronization in networks of dynamical systems with strongly delayed connections. We obtain strict conditions for synchronization of periodic and equilibrium solutions. In particular, we show the existence of a critical coupling strength $\kappa_{c}$, depending only on the network structure, isolated dynamics and coupling function, such that for large delay and coupling strength $\kappa, Comment: Submitted to DCDS-B

We study an excitable active rotator with slowly adapting nonlinear feedback and noise. Depending on the adaptation and the noise level, this system may display noise-induced spiking, noise-perturbed oscillations, or stochastic busting. We show how the system exhibits transitions between these dynamical regimes, as well as how one can enhance or suppress the coherence resonance, or effectively control the features of the stochastic bursting. The setup can be considered as a paradigmatic model for a neuron with a slow recovery variable or, more generally, as an excitable system under the influence of a nonlinear control mechanism. We employ a multiple timescale approach that combines the classical adiabatic elimination with averaging of rapid oscillations and stochastic averaging of noise-induced fluctuations by a corresponding stationary Fokker-Planck equation. This allows us to perform a numerical bifurcation analysis of a reduced slow system and to determine the parameter regions associated with different types of dynamics. In particular, we demonstrate the existence of a region of bistability, where the noise-induced switching between a stationary and an oscillatory regime gives rise to stochastic bursting.

Trabajo presentado en el IFISC Poster Party (online).-- The IFISC Poster Party is an annual activity where PhD students and postdoctoral researchers of IFISC present their research in a poster format.-- Nonlinear Photonics.

The time-delay-based reservoir computing setup has seen tremendous success in both experiment and simulation. It allows for the construction of large neuromorphic computing systems with only few components. However, until now the interplay of the different timescales has not been investigated thoroughly. In this manuscript, we investigate the effects of a mismatch between the time-delay and the clock cycle for a general model. Typically, these two time scales are considered to be equal. Here we show that the case of equal or resonant time-delay and clock cycle could be actively detrimental and leads to an increase of the approximation error of the reservoir. In particular, we can show that non-resonant ratios of these time scales have maximal memory capacities. We achieve this by translating the periodically driven delay-dynamical system into an equivalent network. Networks that originate from a system with resonant delay-times and clock cycles fail to utilize all of their degrees of freedom, which causes the degradation of their performance., S.Y. acknowledges the financial support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Project 411803875. A.R. and K.L. acknowledge support from the Deutsche Forschungsgemeinschaft (Germany) in the framework of the CRC910. F.S. acknowledges financial support provided by the Deutsche Forschungsgemeinschaft through the IRTG 1740. A.R. acknowledges the Spanish State Research Agency, through the Severo Ochoa and Maria de Maeztu Program for Centers and Units of Excellence in R&D (MDM-2017-0711). This paper was developed within the scope of the IRTG 1740 / TRP 2015/50122-0, funded by the DFG / FAPESP .

Adaptive dynamical networks appear in various real-word systems. One of the simplest phenomenological models for investigating basic properties of adaptive networks is the system of coupled phase oscillators with adaptive couplings. In this paper, we investigate the dynamics of this system. We extend recent results on the appearance of hierarchical frequency-multi-clusters by investigating the effect of the time-scale separation. We show that the slow adaptation in comparison with the fast phase dynamics is necessary for the emergence of the multi-clusters and their stability. Additionally, we study the role of double antipodal clusters, which appear to be unstable for all considered parameter values. We show that such states can be observed for a relatively long time, i.e., they are metastable. A geometrical explanation for such an effect is based on the emergence of a heteroclinic orbit., 16 pages, 8 figures

This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in: Berner, R., Yanchuk, S., Maistrenko, Y., & Schöll, E. (2021). Generalized splay states in phase oscillator networks. In Chaos: An Interdisciplinary Journal of Nonlinear Science (Vol. 31, Issue 7, p. 073128). AIP Publishing. https://doi.org/10.1063/5.0056664 and may be found at https://doi.org/10.1063/5.0056664., Networks of coupled phase oscillators play an important role in the analysis of emergent collective phenomena. In this article, we introduce generalized m-splay states constituting a special subclass of phase-locked states with vanishing mth order parameter. Such states typically manifest incoherent dynamics, and they often create high-dimensional families of solutions (splay manifolds). For a general class of phase oscillator networks, we provide explicit linear stability conditions for splay states and exemplify our results with the well-known Kuramoto-Sakaguchi model. Importantly, our stability conditions are expressed in terms of just a few observables such as the order parameter or the trace of the Jacobian. As a result, these conditions are simple and applicable to networks of arbitrary size. We generalize our findings to phase oscillators with inertia and adaptively coupled phase oscillator models. Published under an exclusive license by AIP Publishing.

The recently introduced complex active optical network (LANER) generalizes the concept of laser system to a collection of links, building a bridge with random-laser physics and quantum-graphs theory. So far, LANERs have been studied with a linear approach. Here, we develop a nonlinear formalism in the perspective of describing realistic experimental devices. The propagation along active links is treated via suitable rate equations, which require the inclusion of an auxiliary variable: the population inversion. Altogether, the resulting mathematical model can be viewed as an abstract network, its nodes corresponding to the (directed) fields in the physical links. The dynamical equations differ from standard network models in that, they are a mixture of differential delay (for the active links) and algebraic equations (for the passive links). The stationary states of a generic setup with a single active medium are thoroughly discussed, showing that the role of the passive components can be combined into a single transfer function that takes into account the corresponding resonances., 11 pages, 7 figures

In this article, we analyze a nonlocal ring network of adaptively coupled phase oscillators. We observe a variety of frequency-synchronized states such as phase-locked, multicluster and solitary states. For an important subclass of the phase-locked solutions, the rotating waves, we provide a rigorous stability analysis. This analysis shows a strong dependence of their stability on the coupling structure and the wavenumber which is a remarkable difference to an all-to-all coupled network. Despite the fact that solitary states have been observed in a plethora of dynamical systems, the mechanisms behind their emergence were largely unaddressed in the literature. Here, we show how solitary states emerge due to the adaptive feature of the network and classify several bifurcation scenarios in which these states are created and stabilized.

Localized states are a universal phenomenon observed in spatially distributed dissipative nonlinear systems. Known as dissipative solitons, auto-solitons, spot or pulse solutions, these states play an important role in data transmission using optical pulses, neural signal propagation, and other processes. While this phenomenon was thoroughly studied in spatially extended systems, temporally localized states are gaining attention only recently, driven primarily by applications from fiber or semiconductor lasers. Here we present a theory for temporal dissipative solitons (TDS) in systems with time-delayed feedback. In particular, we derive a system with an advanced argument, which determines the profile of the TDS. We also provide a complete classification of the spectrum of TDS into interface and pseudo-continuous spectrum. We illustrate our theory with two examples: a generic delayed phase oscillator, which is a reduced model for an injected laser with feedback, and the FitzHugh-Nagumo neuron with delayed feedback. Finally, we discuss possible destabilization mechanisms of TDS and show an example where the TDS delocalizes and its pseudo-continuous spectrum develops a modulational instability., 5 pages, 4 figures, submitted (revision) supplementary material available at https://doi.org/10.6084/m9.figshare.8269757 source code for reproducing figures 1a,b, 2a, 3a,b,c available at https://doi.org/10.6084/m9.figshare.8241674

We investigate the relation between the dynamics of a single oscillator with delayed self-feedback and a feed-forward ring of such oscillators, where each unit is coupled to its next neighbor in the same way as in the self-feedback case. We show that periodic solutions of the delayed oscillator give rise to families of rotating waves with different wave numbers in the corresponding ring. In particular, if for the single oscillator the periodic solution is resonant to the delay, it can be embedded into a ring with instantaneous couplings. We discover several cases where the stability of a periodic solution for the single unit can be related to the stability of the corresponding rotating wave in the ring. As a specific example we demonstrate how the complex bifurcation scenario of simultaneously emerging multi-jittering solutions can be transferred from a single oscillator with delayed pulse feedback to multi-jittering rotating waves in a sufficiently large ring of oscillators with instantaneous pulse coupling. Finally, we present an experimental realization of this dynamical phenomenon in a system of coupled electronic circuits of FitzHugh-Nagumo type.

We study the dynamics of an array of nearest-neighbor coupled spatially distributed systems each generating a periodic sequence of short pulses. We demonstrate that unlike a solitary system generating a train of equidistant pulses, an array of such systems can produce a sequence of clusters of closely packed pulses, with the distance between individual pulses depending on the coupling phase. This regime associated with the formation of locally coupled pulse trains bounded due to a balance of attraction and repulsion between them is different from the pulse bound states reported earlier in different laser, plasma, chemical, and biological systems. We propose a simplified analytical description of the observed phenomenon, which is in a good agreement with the results of direct numerical simulations of a model system describing an array of coupled mode-locked lasers., 6 pages, 4 figures. The supplemental material contains: 1.) expressions for the coefficients used in the main paper 2.) video which shows the fixed points and basins of attraction for the bound states in the system of two coupled mode-locked lasers for various coupling phases

Phase response curve (PRC) is an extremely useful tool for studying the response of oscillatory systems, e.g. neurons, to sparse or weak stimulation. Here we develop a framework for studying the response to a series of pulses which are frequent or/and strong so that the standard PRC fails. We show that in this case, the phase shift caused by each pulse depends on the history of several previous pulses. We call the corresponding function which measures this shift the phase response function (PRF). As a result of the introduction of the PRF, a variety of oscillatory systems with pulse interaction, such as neural systems, can be reduced to phase systems. The main assumption of the classical PRC model, i.e. that the effect of the stimulus vanishes before the next one arrives, is no longer a restriction in our approach. However, as a result of the phase reduction, the system acquires memory, which is not just a technical nuisance but an intrinsic property relevant to strong stimulation. We illustrate the PRF approach by its application to various systems, such as Morris-Lecar, Hodgkin-Huxley neuron models, and others. We show that the PRF allows predicting the dynamics of forced and coupled oscillators even when the PRC fails.

We perform bifurcation analysis of plane wave solutions in a one-dimensional complex cubic-quintic Ginzburg--Landau equation with delayed feedback. Our study reveals how multistability and snaking behavior of plane waves emerge as time delay is introduced. For intermediate values of the delay, bifurcation diagrams are obtained by a combination of analytical and numerical methods. For large delays, using an asymptotic approach we classify plane wave solutions into strongly unstable, weakly unstable, and stable. The results of analytical bifurcation analysis are in agreement with those obtained by direct numerical integration of the model equation.

We investigate the dynamics of the pendulum suspended on the forced Duffing oscillator. The detailed bifurcation analysis in two parameter space (amplitude and frequency of excitation) which presents both oscillating and rotating periodic solutions of the pendulum has been performed. We identify the areas with low number of coexisting attractors in the parameter space as the coexistence of different attractors has a significant impact on the practical usage of the proposed system as a tuned mass absorber., Comment: Accepted

Interaction via pulses is common in many natural systems, especially neuronal. In this article we study one of the simplest possible systems with pulse interaction: a phase oscillator with delayed pulsatile feedback. When the oscillator reaches a specific state, it emits a pulse, which returns after propagating through a delay line. The impact of an incoming pulse is described by the oscillator's phase reset curve (PRC). In such a system we discover an unexpected phenomenon: for a sufficiently steep slope of the PRC, a periodic regular spiking solution bifurcates with several multipliers crossing the unit circle at the same parameter value. The number of such critical multipliers increases linearly with the delay and thus may be arbitrary large. This bifurcation is accompanied by the emergence of numerous "jittering" regimes with non-equal interspike intervals (ISIs). Each of these regimes corresponds to a periodic solution of the system with a period roughly proportional to the delay. The number of different "jittering" solutions emerging at the bifurcation point increases exponentially with the delay. We describe the combinatorial mechanism that underlies the emergence of such a variety of solutions. In particular, we show how a periodic solution exhibiting several distinct ISIs can imply the existence of multiple other solutions obtained by rearranging of these ISIs. We show that the theoretical results for phase oscillators accurately predict the behavior of an experimentally implemented electronic oscillator with pulsatile feedback.

We investigate the dynamic mechanisms for low frequency fluctuations in semiconductor lasers subject to delayed optical feedback, using the Lang-Kobayashi model. This system of delay differential equations displays pronounced envelope dynamics, ranging from erratic, so called low frequency fluctuations to regular pulse packages, if the time scales of fast oscillations and envelope dynamics are well separated. We investigate the parameter regions where low frequency fluctuations occur and compute their Lyapunov spectrum. Using geometric singular perturbation theory, we study this intermittent chaotic behavior and characterize these solutions as bursting slow-fast oscillations., Comment: 22 pages, 7 figures

Oscillatory systems with time-delayed pulsatile feedback appear in various applied and theoretical research areas, and received a growing interest in the last years. For such systems, we report a remarkable scenario of destabilization of a periodic regular spiking regime. In the bifurcation point numerous regimes with non-equal interspike intervals emerge simultaneously. We show that this bifurcation is triggered by the steepness of the oscillator's phase resetting curve and that the number of the emerging, so-called "jittering" regimes grows exponentially with the delay value. Although this appears as highly degenerate from a dynamical systems viewpoint, the "multi-jitter" bifurcation occurs robustly in a large class of systems. We observe it not only in a paradigmatic phase-reduced model, but also in a simulated Hodgkin-Huxley neuron model and in an experiment with an electronic circuit.

We show how a variety of stable spatio-temporal periodic patterns can be created in 2D-lattices of coupled oscillators with non-homogeneous coupling delays. The results are illustrated using the FitzHugh-Nagumo coupled neurons as well as coupled limit cycle (Stuart-Landau) oscillators. A “hybrid dispersion relation” is introduced, which describes the stability of the patterns in spatially extended systems with large time-delay.

Dynamical systems with complex delayed interactions arise commonly when propagation times are significant, yielding complicated oscillatory instabilities. In this Letter, we introduce a class of systems with multiple, hierarchically long time delays, and using a suitable space-time representation we uncover features otherwise hidden in their temporal dynamics. The behaviour in the case of two delays is shown to ''encode'' two-dimensional spiral defects and defects turbulence. A multiple scale analysis sets the equivalence to a complex Ginzburg-Landau equation, and a novel criterium for the attainment of the long-delay regime is introduced. We also demonstrate this phenomenon for a semiconductor laser with two delayed optical feedbacks.

© 2014, Springer Science+Business Media New York. We show that for large coupling delays the synchronizability of delay-coupled networks of identical units relates in a simple way to the spectral properties of the network topology. The master stability function used to determine stability of synchronous solutions has a universal structure in the limit of large delay: it is rotationally symmetric around the origin and increases monotonically with the radius in the complex plane. We give details of the proof of this structure and discuss the resulting universal classification of networks with respect to their synchronization properties. We illustrate this classification by means of several prototype network topologies.

For a system of globally pulse-coupled phase-oscillators, we derive conditions for stability of the completely synchronous state and all possible two-cluster states and explain how the different states are naturally connected via bifurcations. The coupling is modeled using the phase-response-curve (PRC), which measures the sensitivity of each oscillator's phase to perturbations. For large systems with a PRC, which turns to zero at the spiking threshold, we are able to find the parameter regions where multiple stable two-cluster states coexist and illustrate this by an example. In addition, we explain how a locally unstable one-cluster state may form an attractor together will its homoclinic connections. This leads to the phenomenon of intermittent, asymptotic synchronization with abating beats away from the perfect synchrony., 12 pages. 6 figures

We study locking of the modulation frequency of a relative periodic orbit in a general $S^1$-equivariant system of ordinary differential equations under an external forcing of modulated wave type. Our main result describes the shape of the locking region in the three-dimensional space of the forcing parameters: intensity, wave frequency, and modulation frequency. The difference of the wave frequencies of the relative periodic orbit and the forcing is assumed to be large and differences of modulation frequencies to be small. The intensity of the forcing is small in the generic case and can be large in the degenerate case, when the first order averaging vanishes. Applications are external electrical and/or optical forcing of selfpulsating states of lasers., Comment: 5 figures

We study chaotic synchronization in networks with time-delayed coupling. We introduce the notion of strong and weak chaos, distinguished by the scaling properties of the maximum Lyapunov exponent within the synchronization manifold for large delay times, and relate this to the condition for stable or unstable chaotic synchronization, respectively. In simulations of laser models and experiments with electronic circuits, we identify transitions from weak to strong and back to weak chaos upon monotonically increasing the coupling strength.

Stability of synchronization in delay-coupled networks of identical units generally depends in a complicated way on the coupling topology. We show that for large coupling delays synchronizability relates in a simple way to the spectral properties of the network topology. The master stability function used to determine stability of synchronous solutions has a universal structure in the limit of large delay: it is rotationally symmetric around the origin and increases monotonically with the radius in the complex plane. This allows a universal classification of networks with respect to their synchronization properties and solves the problem of complete synchronization in networks with strongly delayed coupling., Comment: 4 pages, 1 figure, 1 table. Phys. Rev. Lett. in print, v2 has minor changes in response to PRL referee comments

The goal of this paper is twofold. In the first part we discuss a general approach to determine Lyapunov exponents from ensemble- rather than time-averages. The approach passes through the identification of locally stable and unstable manifolds (the Lyapunov vectors), thereby revealing an analogy with generalized synchronization. The method is then applied to a periodically forced chaotic oscillator to show that the modulus of the Lyapunov exponent associated to the phase dynamics increases quadratically with the coupling strength and it is therefore different from zero already below the onset of phase-synchronization. The analytical calculations are carried out for a model, the generalized special flow, that we construct as a simplified version of the periodically forced Rossler oscillator., Comment: Submitted to Physica D