Refine
Year of publication
Document Type
- Article (285)
- Postprint (11)
- Preprint (9)
- Monograph/Edited Volume (8)
- Other (1)
Keywords
- Complex networks (5)
- Event synchronization (4)
- precipitation (3)
- synchronization (3)
- Amazon rainforest (2)
- Extreme rainfall (2)
- Synchronization (2)
- channel (2)
- classification (2)
- climate networks (2)
Institute
- Institut für Physik und Astronomie (234)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (48)
- Institut für Geowissenschaften (21)
- Department Psychologie (17)
- Institut für Biochemie und Biologie (4)
- Department Linguistik (3)
- Extern (3)
- Institut für Informatik und Computational Science (2)
- Mathematisch-Naturwissenschaftliche Fakultät (2)
- Department Sport- und Gesundheitswissenschaften (1)
- Institut für Umweltwissenschaften und Geographie (1)
We use an index of chaotic synchronization based on the averaged coherence function for the quantitative analysis of the process of the complete synchronization loss in unidirectionally coupled oscillators and maps. We demonstrate that this value manifests different stages of the synchronization breaking. It is invariant to time delay and insensitive to small noise and distortions, which can influence the accessible signals at measurements. Peculiarities of the synchronization destruction in maps and oscillators are investigated
In the modern industrialized countries every year several hundred thousands of people die due to the sudden cardiac death. The individual risk for this sudden cardiac death cannot be defined precisely by common available, non-invasive diagnostic tools like Holter-monitoring, highly amplified ECG and traditional linear analysis of heart rate variability (HRV). Therefore, we apply some rather unconventional methods of nonlinear dynamics to analyse the HRV. Especially, some complexity measures that are basing on symbolic dynamics as well as a new measure, the renormalized entropy, detect some abnormalities in the HRV of several patients who have been classified in the low risk group by traditional methods. A combination of these complexity measures with the parameters in the frequency domain seems to be a promising way to get a more precise definition of the individual risk. These findings have to be validated by a representative number of patients.
During life bones constantly adapt their structure to their mechanical environment via a mechanically controlled process called bone remodeling. For trabecular bone, this process modifies the thickness of each trabecula leading occasionally to full resorption. We describe the irreversible dynamics of the trabecular thickness distribution (TTD) by means of a Markov chain discrete in space and time. By using thickness data from adult patients, we derive the transition probabilities in the chain. This allows a quantification, in terms of geometrical quantities, of the control of bone remodeling and thus to determine the evolution of the TTD with age.
This paper treats a problem of reconstructing ordinary differential equation from a single analytic time series with observational noise. We suppose that the noise is Gaussian (white). The investigation is presented in terms of classical theory of dynamical systems and modern time series analysis. We restrict our considerations on time series obtained as a numerical analytic solution of autonomous ordinary differential equation, solved with respect to the highest derivative and with polynomial right-hand side. In case of an approximate numerical solution with a rather small error, we propose a geometrical basis and a mathematical algorithm to reconstruct a low-order and low-power polynomial differential equation. To reduce the noise the given time series is smoothed at every point by moving polynomial averages using the least-squares method. Then a specific form of the least-squares method is applied to reconstruct the polynomial right-hand side of the unknown equation. We demonstrate for monotonous, periodic and chaotic solutions that this technique is very efficient
Using a special technique of data analysis, we have found out 34 grand minima of solar activity obtained from a 7,700 years long Δ14C record. The method used rests on a proper filtering of the Δ14C record and the extrapolation of verifiable results for the later history back in time. Additionally, we use a method of nonlinear dynamics, the recurrence rate, to back up the results. Our findings are not contradictory to the record of solar maxima resp. minima by Eddy [5], but constitute a considerable extension. Hence, it has become possible to look closer at the validity of models. This way, we have tested several models for solar activity, esp. the model of Barnes et al. [1]. There are hints for that the grand minima might solely be driven by the 209 year period found in the Δ14C record.
Using a special technique of data analysis, we have found out 34 grand minima of solar activity in a 7,700 years long C14 record. The method used rests on a proper filtering of the C14 record and the extrapolation of verifiable results for the later history back in time. Additionally, we have applied a method of nonlinear dynamics, the recurrence rate, to back up the results. Our findings are not contradictory to the record of grand minima by Eddy, but constitute a considerable extension. Hence, it has become possible to look closer at the validity of models. This way, we have tested esp. the model of Barnes et al. There are hints for that the grand minima might solely be driven by the 209--year period found in the C14 record.
Reconstruction of nonlinear time delay models from data by the use of optimal transformations
(1997)
Identifying causal relations from observational data sets has posed great challenges in data-driven causality inference studies. One of the successful approaches to detect direct coupling in the information theory framework is transfer entropy. However, the core of entropy-based tools lies on the probability estimation of the underlying variables. Herewe propose a data-driven approach for causality inference that incorporates recurrence plot features into the framework of information theory. We define it as the recurrence measure of conditional dependence (RMCD), and we present some applications. The RMCD quantifies the causal dependence between two processes based on joint recurrence patterns between the past of the possible driver and present of the potentially driven, excepting the contribution of the contemporaneous past of the driven variable. Finally, it can unveil the time scale of the influence of the sea-surface temperature of the Pacific Ocean on the precipitation in the Amazonia during recent major droughts.
Recurrence-plot-based time series analysis is widely used to study changes and transitions in the dynamics of a system or temporal deviations from its overall dynamical regime. However, most studies do not discuss the significance of the detected variations in the recurrence quantification measures. In this letter we propose a novel method to add a confidence measure to the recurrence quantification analysis. We show how this approach can be used to study significant changes in dynamical systems due to a change in control parameters, chaos-order as well as chaos-chaos transitions. Finally we study and discuss climate transitions by analysing a marine proxy record for past sea surface temperature. This paper is dedicated to the 25th anniversary of the introduction of recurrence plots.
Ventricular tachycardia or fibrillation (VT) as fatal cardiac arrhythmias are the main factors triggering sudden cardiac death. The objective of this recurrence quantification analysis approach is to find early signs of sustained VT in patients with an implanted cardioverter-defibrillator (ICD). These devices are able to safeguard patients by returning their hearts to a normal rhythm via strong defibrillatory shocks; additionally, they are able to store at least 1000 beat-to-beat intervals immediately before the onset of a life-threatening arrhythmia. We study the
Recurrence-plot-based measures of complexity and its application to heart-rate-variability data
(2002)
The knowledge of transitions between regular, laminar or chaotic behavior is essential to understand the underlying mechanisms behind complex systems. While several linear approaches are often insufficient to describe such processes, there are several nonlinear methods which however require rather long time observations. To overcome these difficulties, we propose measures of complexity based on vertical structures in recurrence plots and apply them to the logistic map as well as to heart rate variability data. For the logistic map these measures enable us not only to detect transitions between chaotic and periodic states, but also to identify laminar states, i.e. chaos-chaos transitions. The traditional recurrence quantification analysis fails to detect the latter transitions. Applying our new measures to the heart rate variability data, we are able to detect and quantify the laminar phases before a life-threatening cardiac arrhythmia occurs thereby facilitating a prediction of such an event. Our findings could be of importance for the therapy of malignant cardiac arrhythmias.
We study phase synchronization in a network motif with a starlike structure in which the central node's (the hub's) frequency is strongly detuned against the other peripheral nodes. We find numerically and experimentally a regime of remote synchronization (RS), where the peripheral nodes form a phase synchronized cluster, while the hub remains free with its own dynamics and serves just as a transmitter for the other nodes. We explain the mechanism for this RS by the existence of a free amplitude and also show that systems with a fixed or constant amplitude, such as the classic Kuramoto phase oscillator, are not able to generate this phenomenon. Further, we derive an analytic expression which supports our explanation of the mechanism.
We investigate noise-controlled resonant response of active media to weak periodic forcing, both in excitable and oscillatory regimes. In the excitable regime, we find that noise-induced irregular wave structures can be reorganized into frequency-locked resonant patterns by weak signals with suitable frequencies. The resonance occurs due to a matching condition between the signal frequency and the noise-induced inherent time scale of the media. m:1 resonant regions similar to the Arnold tongues in frequency locking of self-sustained oscillatory media are observed. In the self-sustained oscillatory regime, noise also controls the oscillation frequency and reshapes significantly the Arnold tongues. The combination of noise and weak signal thus could provide an efficient tool to manipulate active extended systems in experiments
The response of scale-free networks with community structure to external stimuli is studied. By disturbing some nodes with different strategies, it is shown that the robustness of this kind of network can be enhanced due to the existence of communities in the networks. Some of the response patterns are found to coincide with topological communities. We show that such phenomena also occur in the cat brain network which is an example of a scale-free like network with community structure. Our results provide insights into the relationship between network topology and the functional organization in complex networks from another viewpoint.
The results of the theoretical consideration of stochastic resonance in overdamped bistable oscillators are given. These results are founded not on the model of two states as in [McNamara B, Wiesenfeld K. Theory of stochastic resonance. Phys Rev A 1989;39:4854-69], but on splitting of motion into regular and random and the rigorous solution of the Fokker-Planck equation for the random component. We show that this resonance is caused by a change, under the influence of noise, of the system's effective stiffness and damping factor contained in the equation for the regular component. For a certain value of the noise intensity the effective stiffness is minimal, and this fact causes non-monotonic change of the output signal amplitude as the noise intensity changes. It is important that the location of the minimum and its value depend essentially on the signal frequency.
Two deterministic processes leading to roughening interfaces are considered. It is shown that the dynamics of linear perturbations of turbulent regimes in coupled map lattices is governed by a discrete version of the Kardar-Parisi-Zhang equation. The asymptotic scaling behavior of the perturbation field is investigated in the case of large lattices. Secondly, the dynamics of an order-disorder interface is modelled with a simple two-dimensional coupled map lattice, possesing a turbulent and a laminar state. It is demonstrated, that in some range of parameters the spreading of the turbulent state is accompanied by kinetic roughening of the interface.
Seismic quiescence as an indicator for large earthquakes in a system of self-organized criticality
(2000)
In this Letter, we show that coherence and phase synchronization analysis are sensitive but not specific in detecting the correct class of underlying dynamics. We propose procedures to increase specificity and demonstrate the power of the approach by application to paradigmatic dynamic model systems. (c) 2006 Elsevier B.V. All rights reserved
Ventricular tachycardia or fibrillation (VT-VF) as fatal cardiac arrhythmias are the main factors triggering sudden cardiac death. The objective of this study is to find early signs of sustained VT-VF in patients with an implanted cardioverter-defibrillator (ICD). These devices are able to safeguard patients by returning their hearts to a normal rhythm via strong defibrillatory shocks; additionally, they store the 1000 beat-to-beat intervals immediately before the onset of a life-threatening arrhythmia. We study these 1000 beat-to-beat intervals of 17 chronic heart failure ICD patients before the onset of a life-threatening arrhythmia and at a control time, i.e., without a VT-VF event. To characterize these rather short data sets, we calculate heart rate variability parameters from the time and frequency domain, from symbolic dynamics as well as the finite-time growth rates. We find that neither the time nor the frequency domain parameters show significant differences between the VT-VF and the control time series. However, two parameters from symbolic dynamics as well as the finite-time growth rates discriminate significantly both groups. These findings could be of importance in algorithms for next generation ICD's to improve the diagnostics and therapy of VT-VF.
We investigate the bifurcation structures in a two-dimensional parameter space (PS) of a parametrically excited system with two degrees of freedom both analytically and numerically. By means of the Renyi entropy of second order K-2, which is estimated from recurrence plots, we uncover that regions of chaotic behavior are intermingled with many complex periodic windows, such as shrimp structures in the PS. A detailed numerical analysis shows that, the stable solutions lose stability either via period doubling, or via intermittency when the parameters leave these shrimps in different directions, indicating different bifurcation properties of the boundaries. The shrimps of different sizes offer promising ways to control the dynamics of such a complex system.
Acoustic emission signals generated during high speed cutting of steel are investigated. The data are represen ted in time-folded form. Several methods from linear and nonlinear data analysis based on time- and frequency- domain are applied to the data and reveal signatures of the observed acoustic emission signal. These investiga tions are necessary for modeling the cutting process by means of differential equations.
Similar power laws for foreshock and aftershock sequences in a spring block model for earthquakes
(1999)
We have recently reported the phenomenon of doubly stochastic resonance [Phys. Rev. Lett. 85, 227 (2000)], a synthesis of noise-induced transition and stochastic resonance. The essential feature of this phenomenon is that multiplicative noise induces a bimodality and additive noise causes stochastic resonance behavior in the induced structure. In the present paper we outline possible applications of this effect and design a simple lattice of electronic circuits for the experimental realization of doubly stochastic resonance.
Understanding the functional dynamics of the mammalian brain is one of the central aims of modern neuroscience. Mathematical modeling and computational simulations of neural networks can help in this quest. In recent publications, a multilevel model has been presented to simulate the resting-state dynamics of the cortico-cortical connectivity of the mammalian brain. In the present work we investigate how much of the dynamical behavior of the multilevel model can be reproduced by a strongly simplified model. We find that replacing each cortical area by a single Rulkov map recreates the patterns of dynamical correlation of the multilevel model, while the outcome of other models and setups mainly depends on the local network properties, e. g. the input degree of each vertex. In general, we find that a simple simulation whose dynamics depends on the global topology of the whole network is far from trivial. A systematic analysis of different dynamical models and coupling setups is required.
We demonstrate the occurrence of regimes with singular continuous (fractal) Fourier spectra in autonomous dissipative dynamical systems. The particular example in an ODE system at the accumulation points of bifurcation sequences associated to the creation of complicated homoclinic orbits. Two different machanisms responsible for the appearance of such spectra are proposed. In the first case when the geometry of the attractor is symbolically represented by the Thue-Morse sequence, both the continuous-time process and its descrete Poincaré map have singular power spectra. The other mechanism owes to the logarithmic divergence of the first return times near the saddle point; here the Poincaré map possesses the discrete spectrum, while the continuous-time process displays the singular one. A method is presented for computing the multifractal characteristics of the singular continuous spectra with the help of the usual Fourier analysis technique.
Spatial recurrence plots
(2006)
We propose an extension of the recurrence plot concept to perform quantitative analyzes of roughness and disorder of spatial patterns at a fixed time. We introduce spatial recurrence plots (SRPs) as a graphical representation of the pointwise correlation matrix, in terms of a two-dimensional spatial return plot. This technique is applied to the study of complex patterns generated by coupled map lattices, which are characterized by measures of complexity based on SRPs. We show that the complexity measures we propose for SRPs provide a systematic way of investigating the distribution of spatially coherent structures, such as synchronization domains, in lattice profiles. This approach has potential for many more applications, e.g., in surface roughness analyzes
The South American Andes are frequently exposed to intense rainfall events with varying moisture sources and precipitation-forming processes. In this study, we assess the spatiotemporal characteristics and geographical origins of rainfall over the South American continent. Using high-spatiotemporal resolution satellite data (TRMM 3B42 V7), we define four different types of rainfall events based on their (1) high magnitude, (2) long temporal extent, (3) large spatial extent, and (4) high magnitude, long temporal and large spatial extent combined. In a first step, we analyze the spatiotemporal characteristics of these events over the entire South American continent and integrate their impact for the main Andean hydrologic catchments. Our results indicate that events of type 1 make the overall highest contributions to total seasonal rainfall (up to 50%). However, each consecutive episode of the infrequent events of type 4 still accounts for up to 20% of total seasonal rainfall in the subtropical Argentinean plains. In a second step, we employ complex network theory to unravel possibly non-linear and long-ranged climatic linkages for these four event types on the high-elevation Altiplano-Puna Plateau as well as in the main river catchments along the foothills of the Andes. Our results suggest that one to two particularly large squall lines per season, originating from northern Brazil, indirectly trigger large, long-lasting thunderstorms on the Altiplano Plateau. In general, we observe that extreme rainfall in the catchments north of approximately 20 degrees S typically originates from the Amazon Basin, while extreme rainfall at the eastern Andean foothills south of 20 degrees S and the Puna Plateau originates from southeastern South America.
We study the effects of parametric noise on a lattice network, which is locally modeled by a two-dimensional Rulkov map. We conclude that at some intermediate noise intensity, parametric noise can induce ordered circular patterns, which indicates the appearance of spatiotemporal coherence resonance in the studied lattice. With the observation of coherence-like manner in linear spatial cross-correlation, the coherence phenomena can be analyzed quantitatively.
Spatiotemporal dynamics of the Calvin cycle multistationarity and symmetry breaking instabilities
(2011)
The possibility of controlling the Calvin cycle has paramount implications for increasing the production of biomass. Multistationarity, as a dynamical feature of systems, is the first obvious candidate whose control could find biotechnological applications. Here we set out to resolve the debate on the multistationarity of the Calvin cycle. Unlike the existing simulation-based studies, our approach is based on a sound mathematical framework, chemical reaction network theory and algebraic geometry, which results in provable results for the investigated model of the Calvin cycle in which we embed a hierarchy of realistic kinetic laws. Our theoretical findings demonstrate that there is a possibility for multistationarity resulting from two sources, homogeneous and inhomogeneous instabilities, which partially settle the debate on multistability of the Calvin cycle. In addition, our tractable analytical treatment of the bifurcation parameters can be employed in the design of validation experiments.
We study numerically the behavior of the autocorrelation function (ACF) and the power spectrum of spiral attractors without and in the presence of noise. It is shown that the ACF decays exponentially and has two different time scales. The rate of the ACF decrease is defined by the amplitude fluctuations on small time intervals, i.e., when tau < tau(cor), and by the effective diffusion coefficient of the instantantaneous phase on large time intervals. it is also demonstrated that the ACF in the Poincare map also decreases according to the exponential law exp(-lambda(+)k), where lambda(+) is the positive Lyapunov exponent. The obtained results are compared with the theory of fluctuations for the Van der Pol oscillator
We employ a spectral decomposition method to analyze synchronization of a non-identical oscillator network. We study the case that a small parameter mismatch of oscillators is characterized by one parameter and phase synchronization is observed. We derive a linearized equation for each eigenmode of the coupling matrix. The parameter mismatch is reflected on inhomogeneous term in the linearized equation. We find that the oscillation of each mode is essentially characterized only by the eigenvalue of the coupling matrix with a suitable normalization. We refer to this property as spectral universality, because it is observed irrespective of network topology. Numerical results in various network topologies show good agreement with those based on linearized equation. This universality is also observed in a system driven by additive independent Gaussian noise.
In this paper we show that delay embedding produces spurious structures in a recurrence plot (RP) that are not present in the real attractor. We analyze typical sets of simulated data, such as white noise and data from the chaotic Rossler system to show the relevance of this effect. In the second part of the paper we show that the second order Renyi entropy and the correlation dimension are dynamical invariants that can be estimated from Recurrence Plots with arbitrary embedding dimension and delay
Stochastic bifurcations and coherencelike resonance in a self-sustained bistable noisy oscillator
(2010)
We investigate the influence of additive Gaussian white noise on two different bistable self-sustained oscillators: Duffing-Van der Pol oscillator with hard excitation and a model of a synthetic genetic oscillator. In the deterministic case, both oscillators are characterized with a coexistence of a stable limit cycle and a stable equilibrium state. We find that under the influence of noise, their dynamics can be well characterized through the concept of stochastic bifurcation, consisting in a qualitative change of the stationary amplitude distribution. For the Duffing-Van der Pol oscillator analytical results, obtained for a quasiharmonic approach, are compared with the result of direct computer simulations. In particular, we show that the dynamics is different for isochronous and anisochronous systems. Moreover, we find that the increase of noise intensity in the isochronous regime leads to a narrowing of the spectral line. This effect is similar to coherence resonance. However, in the case of anisochronous systems, this effect breaks down and a new phenomenon, anisochronous-based stochastic bifurcation occurs.
Charged dust grains in circumplanetary environments experience, beyond various deterministic forces, also stochastic perturbations caused, by fluctuations of the magnetic field, the charge of the grains, by chaotic rotation of aspherical grains, etc. Here we investigate the dynamics of a dust population in a circular orbit around a planet which is perturbed by a stochastic planetary magnetic field B', modeled by an isotropically Gaussian white noise. The resulting perturbation equations give rise to a modified diffusion of the inclinations i and eccentricities e. The diffusion coefficient is found to be D proportional to w^2 O /n^2 , where the gyrofrequency, the Kepler frequency, and the synodic frequency are denoted by w , O, and n, respectively. This behavior has been checked against numerical simulations. We have chosen dust grains (1 m in radius) ejected from Jupiter's satellite Europa in circular equatorial orbits around Jupiter and integrated numerically their trajectories over their typical lifetimes (100 years). The particles were exposed to a Gaussian fluctuating magnetic field B' with the same statistical properties as in the analytical treatment. These simulations have confirmed the analytical results. The theoretical studies showed the statistical properties of B' to be of decisive importance. To estimate them, we analyzed the magnetic field data obtained by the Galileo spacecraft magnetometer at Jupiter and found almost Gaussian fluctuations of about 5% of the mean field and exponentially decaying correlations. This results in a diffusion of orbital inclinations and eccentricities of the dust grains of about ten percent over the lifetime of the particles. For smaller dusty motes or for close-in particles (e.g., in Jovian gossamer rings) stochastics might well dominate the dynamics.
We study the noise-dependent dynamics in a chain of four very stiff excitable oscillators of the FitzHugh- Nagumo type locally coupled by inhibitor diffusion. We could demonstrate frequency- and noise-selective signal acceptance which is based on several noise-supported stochastic attractors that arise owing to slow variable diffusion between identical excitable elements. The attractors have different average periods distinct from that of an isolated oscillator and various phase relations between the elements. We explain the correspondence between the noise-supported stochastic attractors and the observed resonance peaks in the curves for the linear response versus signal frequency. (C) 2005 American Institute of Physics
Recent research using the complex network approach has revealed a rich and complicated network topology in the cortical connectivity of mammalian brains. It is of importance to understand the implications of such complex network structures in the functional organization of the brain activities. Here we study this problem from the viewpoint of dynamical complex networks. We investigate synchronization dynamics on the corticocortical network of the cat by modeling each node (cortical area) of the network with a sub-network of interacting excitable neurons. We find that the network displays clustered synchronization behavior, and the dynamical clusters coincide with the topological community structures observed in the anatomical network. Our results provide insights into the relationship between the global organization and the functional specialization of the brain cortex.
Fourier surrogate data are artificially generated time series, that - based on a resampling scheme - share the linear properties with an observed time series. In this paper we study a statistical surrogate hypothesis test to detect deviations from a linear Gaussian process with respect to asymmetry in time (Q-statistic). We apply this test to a Fourier representable function and obtain a representation of the asymmetry in time of the sample data, a characteristic for nonlinear processes, and the significance in terms of the Fourier coefficients. The main outcome is that we calculate the expected value of the mean and the standard deviation of the asymmetries of the surrogate data analytically and hence, no surrogates have to be generated. To illustrate the results we apply our method to the saw tooth function, the Lorenz system and to measured X-ray data of Cygnus X-1
Observational data of natural systems, as measured in medical measurements are typically quite different from those obtained in laboratories. Due to the peculiarities of these data, wellknown characteristics, such as power spectra or fractal dimension, often do not provide a suitable description. To study such data, we present here some measures of complexity, which are basing on symbolic dynamics. Firstly, a motivation for using symbolic dynamics and measures of complexity in data analysis based on the logistic map is given and next, two applications to medical data are shown. We demonstrate that symbolic dynamics is a useful tool for the risk assessment of patients after myocardial infarction as well as for the evaluation of th e architecture of human cancellous bone.
We analyse time series from a study on bimanual rhythmic movements in which the speed of performance (the external control parameter) was experimentally manipulated. Using symbolic transformations as a visualization technique we observe qualitative changes in the dynamics of the timing patterns. Such phase transitions are quantitatively described by measures of complexity. Using these results we develop an advanced symbolic coding which enables us to detect important dynamical structures. Furthermore, our analysis raises new questions concerning the modelling of the underlying human cognitive-motor system.