Refine
Has Fulltext
- no (294) (remove)
Year of publication
Document Type
- Article (285)
- Monograph/Edited Volume (8)
- Other (1)
Is part of the Bibliography
- yes (294)
Keywords
- Complex networks (5)
- Event synchronization (4)
- Extreme rainfall (2)
- Synchronization (2)
- synchronization (2)
- 3D medical image analysis (1)
- African climate (1)
- Algebraic geometry (1)
- Amazon rainforest (1)
- Anisotropy (1)
Institute
- Institut für Physik und Astronomie (217)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (39)
- Institut für Geowissenschaften (20)
- Department Psychologie (13)
- Institut für Biochemie und Biologie (4)
- Department Linguistik (3)
- Institut für Informatik und Computational Science (2)
- Department Sport- und Gesundheitswissenschaften (1)
- Institut für Umweltwissenschaften und Geographie (1)
The main intention of this contribution is to discuss different nonlinear approaches to heart rate and blood pressure variability analysis for a better understanding of the cardiovascular regulation. We investigate measures of complexity which are based on symbolic dynamics, renormalised entropy and the finite time growth rates. The dual sequence method to estimate the baroreflex sensitivity and the maximal correlation method to estimate the nonlinear coupling between time series are employed for analysing bivariate data. The latter appears to be a suitable method to estimate the strength of the nonlinear coupling and the coupling direction. Heart rate and blood pressure data from clinical pilot studies and from very large clinical studies are analysed. We demonstrate that parameters from nonlinear dynamics are useful for risk stratification after myocardial infarction, for the prediction of life-threatening cardiac events even in short time series, and for modelling the relationship between heart rate and blood pressure regulation. These findings could be of importance for clinical diagnostics, in algorithms for risk stratification, and for therapeutic and preventive tools of next generation implantable cardioverter defibrillators.
Nonlinear detection of paleoclimate-variability transitions possibly related to human evolution
(2011)
Potential paleoclimatic driving mechanisms acting on human evolution present an open problem of cross-disciplinary scientific interest. The analysis of paleoclimate archives encoding the environmental variability in East Africa during the past 5 Ma has triggered an ongoing debate about possible candidate processes and evolutionary mechanisms. In this work, we apply a nonlinear statistical technique, recurrence network analysis, to three distinct marine records of terrigenous dust flux. Our method enables us to identify three epochs with transitions between qualitatively different types of environmental variability in North and East Africa during the (i) Middle Pliocene (3.35-3.15 Ma B. P.), (ii) Early Pleistocene (2.25-1.6 Ma B. P.), and (iii) Middle Pleistocene (1.1-0.7 Ma B. P.). A deeper examination of these transition periods reveals potential climatic drivers, including (i) large-scale changes in ocean currents due to a spatial shift of the Indonesian throughflow in combination with an intensification of Northern Hemisphere glaciation, (ii) a global reorganization of the atmospheric Walker circulation induced in the tropical Pacific and Indian Ocean, and (iii) shifts in the dominating temporal variability pattern of glacial activity during the Middle Pleistocene, respectively. A reexamination of the available fossil record demonstrates statistically significant coincidences between the detected transition periods and major steps in hominin evolution. This result suggests that the observed shifts between more regular and more erratic environmental variability may have acted as a trigger for rapid change in the development of humankind in Africa.
Linear methods of dimensionality reduction are useful tools for handling and interpreting high dimensional data. However, the cumulative variance explained by each of the subspaces in which the data space is decomposed may show a slow convergence that makes the selection of a proper minimum number of subspaces for successfully representing the variability of the process ambiguous. The use of nonlinear methods can improve the embedding of multivariate data into lower dimensional manifolds. In this article, a nonlinear method for dimensionality reduction, Isomap, is applied to the sea surface temperature and thermocline data in the tropical Pacific Ocean, where the El Nino-Southern Oscillation (ENSO) phenomenon and the annual cycle phenomena interact. Isomap gives a more accurate description of the manifold dimensionality of the physical system. The knowledge of the minimum number of dimensions is expected to improve the development of low dimensional models for understanding and predicting ENSO
We review the problem of estimating parameters and unobserved trajectory components from noisy time series measurements of continuous nonlinear dynamical systems. It is first shown that in parameter estimation techniques that do not take the measurement errors explicitly into account, like regression approaches, noisy measurements can produce inaccurate parameter estimates. Another problem is that for chaotic systems the cost functions that have to be minimized to estimate states and parameters are so complex that common optimization routines may fail. We show that the inclusion of information about the time-continuous nature of the underlying trajectories can improve parameter estimation considerably. Two approaches, which take into account both the errors-in-variables problem and the problem of complex cost functions, are described in detail: shooting approaches and recursive estimation techniques. Both are demonstrated on numerical examples
In this paper we apply symbolic transformations as a visualisation technique for analysing rhythm production. It is shown that qualitative information can be extracted from the experimental data. This approach may provide new insights into the organisation of temporal order by the brain on different levels of description. A simple phenomenological model for the explanation of the observed phenomena is proposed.
We analyse the X-ray light curves of compact objects using linear and nonlinear time series analysis methods. A Power Density Spectrum (PDS) describes the overall second order properties of the observed data well. To look beyond we propose the nonlinear Q-statistic to detect an asymmetry of the time series. This allows us to find relevant time scales. This method even grants a subclassification of the known states of X-ray sources.
In integrated medical considerations of the biological human system, both intellectual and motor performances in a similar manner are considered as a result of the function of the nervous system. Consequently, universal minimal dysfunctions of the central nervous system may lead to both intellectual and physical anomalies. Therefore, this study tests the hypothesis that there is a connection between the balance ability as a motor parameter and school success as an intellectual parameter. A postural measuring system based on the force-moment sensor technique was used to record the postural balance regulation of 773 children (circle divide 11 +/- 1 years). The school achievement of each child was determined by school grades. Data analysis was performed by linear as well as by nonlinear time series analyses. There are highly significant differences in balance regulation between good and poor pupils recognized by several linear and nonlinear parameters. Good pupils could be discriminated from pupils with bad results in learning to 80 %. The results support the hypothesis mentioned above. One possible explanation for the poor regulation of balance in bad learners could be a deficit in the neural maturity. In future, further developments will be targeted on higher discrimination levels, possibly in order to predict school success. On the other hand, the effects of special movement exercises on the neural development in childhood will be the focus in our further work
Hydrometric networks play a vital role in providing information for decision-making in water resource management. They should be set up optimally to provide as much information as possible that is as accurate as possible and, at the same time, be cost-effective. Although the design of hydrometric networks is a well-identified problem in hydrometeorology and has received considerable attention, there is still scope for further advancement. In this study, we use complex network analysis, defined as a collection of nodes interconnected by links, to propose a new measure that identifies critical nodes of station networks. The approach can support the design and redesign of hydrometric station networks. The science of complex networks is a relatively young field and has gained significant momentum over the last few years in different areas such as brain networks, social networks, technological networks, or climate networks. The identification of influential nodes in complex networks is an important field of research. We propose a new node-ranking measure – the weighted degree–betweenness (WDB) measure – to evaluate the importance of nodes in a network. It is compared to previously proposed measures used on synthetic sample networks and then applied to a real-world rain gauge network comprising 1229 stations across Germany to demonstrate its applicability. The proposed measure is evaluated using the decline rate of the network efficiency and the kriging error. The results suggest that WDB effectively quantifies the importance of rain gauges, although the benefits of the method need to be investigated in more detail.
This paper discusses translocation features of the 20S proteasome in order to explain typical proteasome length distributions. We assume that the protein transport depends significantly on the fragment length with some optimal length which is transported most efficiently. By means of a simple one-channel model, we show that this hypothesis can explain both the one- and the three-peak length distributions found in experiments. A possible mechanism of such translocation is provided by so-called fluctuation-driven transport.
Graphical models applying partial coherence to multivariate time series are a powerful tool to distinguish direct and indirect interdependencies in multivariate linear systems. We carry over the concept of graphical models and partialization analysis to phase signals of nonlinear synchronizing systems. This procedure leads to the partial phase synchronization index which generalizes a bivariate phase synchronization index to the multivariate case and reveals the coupling structure in multivariate synchronizing systems by differentiating direct and indirect interactions. This ensures that no false positive conclusions are drawn concerning the interaction structure in multivariate synchronizing systems. By application to the paradigmatic model of a coupled chaotic Roessler system, the power of the partial phase synchronization index is demonstrated
This work focuses on the dynamics of particles in a confined geometry with position-dependent diffusivity, where the confinement is modelled by a periodic channel consisting of unit cells connected by narrow passage ways. We consider three functional forms for the diffusivity, corresponding to the scenarios of a constant (D ₀), as well as a low (D ₘ) and a high (D d) mobility diffusion in cell centre of the longitudinally symmetric cells. Due to the interaction among the diffusivity, channel shape and external force, the system exhibits complex and interesting phenomena. By calculating the probability density function, mean velocity and mean first exit time with the Itô calculus form, we find that in the absence of external forces the diffusivity D d will redistribute particles near the channel wall, while the diffusivity D ₘ will trap them near the cell centre. The superposition of external forces will break their static distributions. Besides, our results demonstrate that for the diffusivity D d, a high dependence on the x coordinate (parallel with the central channel line) will improve the mean velocity of the particles. In contrast, for the diffusivity D ₘ, a weak dependence on the x coordinate will dramatically accelerate the moving speed. In addition, it shows that a large external force can weaken the influences of different diffusivities; inversely, for a small external force, the types of diffusivity affect significantly the particle dynamics. In practice, one can apply these results to achieve a prominent enhancement of the particle transport in two- or three-dimensional channels by modulating the local tracer diffusivity via an engineered gel of varying porosity or by adding a cold tube to cool down the diffusivity along the central line, which may be a relevant effect in engineering applications. Effects of different stochastic calculi in the evaluation of the underlying multiplicative stochastic equation for different physical scenarios are discussed.
When we fixate our gaze on a stable object, our eyes move continuously with extremely small involuntary and autonomic movements, that even we are unaware of during their occurrence. One of the roles of these fixational eye movements is to prevent the adaptation of the visual system to continuous illumination and inhibit fading of the image. These random, small movements are restricted at long time scales so as to keep the target at the centre of the field of view. In addition, the synchronisation properties between both eyes are related to binocular coordination in order to provide stereopsis. We investigated the roles of different time scale behaviours, especially how they are expressed in the different spatial directions (vertical versus horizontal). We also tested the synchronisation between both eyes. Results show different scaling behaviour between horizontal and vertical movements. When the small ballistic movements, i.e., microsaccades, are removed, the scaling behaviour in both axes becomes similar. Our findings suggest that microsaccades enhance the persistence at short time scales mostly in the horizontal component and much less in the vertical component. We also applied the phase synchronisation decay method to study the synchronisation between six combinations of binocular fixational eye movement components. We found that the vertical-vertical components of right and left eyes are significantly more synchronised than the horizontal-horizontal components. These differences may be due to the need for continuously moving the eyes in the horizontal plane in order to match the stereoscopic image for different viewing distances.
We study the effects of mutual and external chaotic phase synchronization in ensembles of bursting oscillators. These oscillators (used for modeling neuronal dynamics) are essentially multiple time scale systems. We show that a transition to mutual phase synchronization takes place on the bursting time scale of globally coupled oscillators, while on the spiking time scale, they behave asynchronously. We also demonstrate the effect of the onset of external chaotic phase synchronization of the bursting behavior in the studied ensemble by a periodic driving applied to one arbitrarily taken neuron. We also propose an explanation of the mechanism behind this effect. We infer that the demonstrated phenomenon can be used efficiently for controlling bursting activity in neural ensembles
We study phase synchronization effects of chaotic oscillators with a type-I intermittency behavior. The external and mutual locking of the average length of the laminar stage for coupled discrete and continuous in time systems is shown and the mechanism of this synchronization is explained. We demonstrate that this phenomenon can be described by using results of the parametric resonance theory and that this correspondence enables one to predict and derive all zones of synchronization
We demonstrate the existence of phase synchronization of two chaotic rotators. Contrary to phase synchronization of chaotic oscillators, here the Lyapunov exponents corresponding to both phases remain positive even in the synchronous regime. Such frequency locked dynamics with different ratios of frequencies are studied for driven continuous-time rotators and for discrete circle maps. We show that this transition to phase synchronization occurs via a crisis transition to a band-structured attractor.
An approach is presented for the reconstruction of phase synchronization phenomena in a chaotic CO2 laser from experimental data. We analyze this laser system in a regime able to phase synchronize with a weak sinusoidal forcing. Our technique recovers the synchronization diagram of the experimental system from only few measurement data sets, thus allowing the prediction of the regime of phase synchronization as well as nonsynchronization in a broad parameter space of forcing frequency and amplitude without further experiments
In the last decade, there has been an increasing interest in compensating thermally induced errors to improve the manufacturing accuracy of modular tool systems. These modular tool systems are interfaces between spindle and workpiece and consist of several complicatedly formed parts. Their thermal behavior is dominated by nonlinearities, delay and hysteresis effects even in tools with simpler geometry and it is difficult to describe it theoretically. Due to the dominant nonlinear nature of this behavior the so far used linear regression between the temperatures and the displacements is insufficient. Therefore, in this study we test the hypothesis whether we can reliably predict such thermal displacements via nonlinear temperature-displacement regression functions. These functions are estimated firstly from learning measurements using the alternating conditional expectation (ACE) algorithm and then tested on independent data sets. First, we analyze data that were generated by a finite element spindle model. We find that our approach is a powerful tool to describe the relation between temperatures and displacements for simulated data. Next, we analyze the temperature-displacement relationship in a silent real experimental setup, where the tool system is thermally forced. Again, the ACE-algorithm is powerful to estimate the deformation with high precision. The corresponding errors obtained by using the nonlinear regression approach are 10-fold lower in comparison to multiple linear regression analysis. Finally, we investigate the thermal behavior of a modular tool system in a working milling machine and get again promising results. The thermally inducedaccuracy using this nonlinear regression analysis. Therefore, this approach seems to be very useful for the development of new modular tool systems. errors can be estimated with 1-2 micrometer
Based on high-spatiotemporal-resolution data, the authors perform a climatological study of strong rainfall events propagating from southeastern South America to the eastern slopes of the central Andes during the monsoon season. These events account for up to 70% of total seasonal rainfall in these areas. They are of societal relevance because of associated natural hazards in the form of floods and landslides, and they form an intriguing climatic phenomenon, because they propagate against the direction of the low-level moisture flow from the tropics. The responsible synoptic mechanism is analyzed using suitable composites of the relevant atmospheric variables with high temporal resolution. The results suggest that the low-level inflow from the tropics, while important for maintaining sufficient moisture in the area of rainfall, does not initiate the formation of rainfall clusters. Instead, alternating low and high pressure anomalies in midlatitudes, which are associated with an eastward-moving Rossby wave train, in combination with the northwestern Argentinean low, create favorable pressure and wind conditions for frontogenesis and subsequent precipitation events propagating from southeastern South America toward the Bolivian Andes.
Changes in trabecular bone composition during development of osteoporosis are used as a model for bone loss in microgravity conditions during a space flight. Symbolic dynamics and measures of complexity are proposed and applied to assess quantitatively the structural composition of bone tissue from 3D data sets of human tibia bone biopsies acquired by a micro-CT scanner. In order to justify the newly proposed approach, the measures of complexity of the bone architecture were compared with the results of traditional 2D bone histomorphometry. The proposed technique is able to quantify the structural loss of the bone tissue and may help to diagnose and to monitor changes in bone structure of patients on Earth as well as of the space-flying personnel. © 2005 Elsevier Ltd. All rights reserved
Quantifying the roles of single stations within homogeneous regions using complex network analysis
(2018)
Regionalization and pooling stations to form homogeneous regions or communities are essential for reliable parameter transfer, prediction in ungauged basins, and estimation of missing information. Over the years, several clustering methods have been proposed for regional analysis. Most of these methods are able to quantify the study region in terms of homogeneity but fail to provide microscopic information about the interaction between communities, as well as about each station within the communities. We propose a complex network-based approach to extract this valuable information and demonstrate the potential of our approach using a rainfall network constructed from the Indian gridded daily precipitation data. The communities were identified using the network-theoretical community detection algorithm for maximizing the modularity. Further, the grid points (nodes) were classified into universal roles according to their pattern of within- and between-community connections. The method thus yields zoomed-in details of individual rainfall grids within each community.
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
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 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.
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.
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.
We present two different approaches to detect and quantify phase synchronization in the case of coupled non- phase coherent oscillators. The first one is based on the general idea of curvature of an arbitrary curve. The second one is based on recurrences of the trajectory in phase space. We illustrate both methods in the paradigmatic example of the Rossler system in the funnel regime. We show that the second method is applicable even in the case of noisy data. Furthermore, we extend the second approach to the application of chains of coupled systems, which allows us to detect easily clusters of synchronized oscillators. In order to illustrate the applicability of this approach, we show the results of the algorithm applied to experimental data from a population of 64 electrochemical oscillators
In this article we review the application of the synchronization theory to the analysis of multivariate biological signals. We address the problem of phase estimation from data and detection and quantification of weak interaction, as well as quantification of the direction of coupling. We discuss the potentials as well as limitations and misinterpretations of the approach
Starting from an initial wiring of connections, we show that the synchronizability of a network can be significantly improved by evolving the graph along a time dependent connectivity matrix. We consider the case of connectivity matrices that commute at all times, and compare several approaches to engineer the corresponding commutative graphs. In particular, we show that synchronization in a dynamical network can be achieved even in the case in which each individual commutative graphs does not give rise to synchronized behavior
In this paper, we study the complete synchronization of a class of time-varying delayed coupled chaotic systems using feedback control. In terms of Linear Matrix Inequalities, a sufficient condition is obtained through using a Lyapunov-Krasovskii functional and differential equation in equalities. The conditions can be easily verified and implemented. We present two simulation examples to illustrate the effectiveness of the proposed method.
In the present paper, two kinds of dynamical complex networks are considered. The first is that elements of every node have different time delays but all nodes in Such networks have the same time-delay vector. The second is that different nodes have different time-delay vectors, and the elements of each node also have different time delays. Corresponding synchronization theorems are established. Numerical examples show the efficiency of the derived theorems.
We present experimental and numerical evidence of synchronization of burst events in two different modulated CO2 lasers. Bursts appear randomly in each laser as trains of large amplitude spikes intercalated by a small amplitude chaotic regime. Experimental data and model show the frequency locking of bursts in a suitable interval of coupling strength. We explain the mechanism of this phenomenon and demonstrate the inhibitory properties of the implemented coupling.
We analyze synchronization between two interacting populations of different phase oscillators. For the important case of asymmetric coupling functions, we find a much richer dynamical behavior compared to that of symmetrically coupled populations of identical oscillators. It includes three types of bistabilities, higher order entrainment and the existence of states with unusual stability properties. All possible routes to synchronization of the populations are presented and some stability boundaries are obtained analytically. The impact of these findings for neuroscience is discussed.
Being one of the fundamental phenomena in nonlinear science, synchronization of oscillations has permanently remained an object of intensive research. Development of many asymptotic methods and numerical simulations has allowed an understanding and explanation of various phenomena of self-synchronization. But even in the classical case of coupled van der Pol oscillators a full description of all possible dynamical regimes, their mutual transitions and characteristics is still lacking. We present here a study of the phenomenon of mutual synchronization for two non-scalar- coupled non-identical limit-cycle oscillators and analyze phase, frequency and amplitude characteristics of synchronization regimes. A series of bifurcation diagrams that we obtain exhibit various regions of qualitatively different behavior. Among them we find mono-, bi- and multistability regions, beating and "oscillation death" ones; also a region, where one of the oscillators dominates the other one is observed. The frequency characteristics that we obtain reveal three qualitatively different types of synchronization: (i) on the mean frequency (the in-phase synchronization), (ii) with a shift from the mean frequency caused by a conservative coupling term (the anti-phase synchronization), and (iii) on the frequency of one of the oscillators (when one oscillator dominates the other). (C) 2003 Elsevier B.V. All rights reserved