Detecting and characterizing phase synchronization in nonstationary dynamical systems

multiSyncPy: A Python package for assessing multivariate coordination dynamics

In order to support the burgeoning field of research into intra- and interpersonal synchrony, we present an open-source software package: multiSyncPy. Multivariate synchrony goes beyond the bivariate case and can be useful for quantifying how groups, teams, and families coordinate their behaviors, or estimating the degree to which multiple modalities from an individual become synchronized. Our package includes state-of-the-art multivariate methods including symbolic entropy, multidimensional recurrence quantification analysis, coherence (with an additional sum-normalized modification), the cluster-phase 'Rho' metric, and a statistical test based on the Kuramoto order parameter. We also include functions for two surrogation techniques to compare the observed coordination dynamics with chance levels and a windowing function to examine time-varying coordination for most of the measures. Taken together, our collation and presentation of these methods make the study of interpersonal synchronization and coordination dynamics applicable to larger, more complex and often more ecologically valid study designs. In this work, we summarize the relevant theoretical background and present illustrative practical examples, lessons learned, as well as guidance for the usage of our package - using synthetic as well as empirical data. Furthermore, we provide a discussion of our work and software and outline interesting further directions and perspectives. multiSyncPy is freely available under the LGPL license at: https://github.com/cslab-hub/multiSyncPy , and also available at the Python package index.

Windowed multiscale synchrony: modeling time-varying and scale-localized interpersonal coordination dynamics

Social interactions are pervasive in human life with varying forms of interpersonal coordination emerging and spanning different modalities (e.g. behaviors, speech/language, and neurophysiology). However, during social interactions, as in any dynamical system, patterns of coordination form and dissipate at different scales. Historically, researchers have used aggregate measures to capture coordination over time. While those measures (e.g. mean relative phase, cross-correlation, coherence) have provided a wealth of information about coordination in social settings, some evidence suggests that multiscale coordination may change over the time course of a typical empirical observation. To address this gap, we demonstrate an underutilized method, windowed multiscale synchrony, that moves beyond quantifying aggregate measures of coordination by focusing on how the relative strength of coordination changes over time and the scales that comprise social interaction. This method involves using a wavelet transform to decompose time series into component frequencies (i.e. scales), preserving temporal information and then quantifying phase synchronization at each of these scales. We apply this method to both simulated and empirical interpersonal physiological and neuromechanical data. We anticipate that demonstrating this method will stimulate new insights on the mechanisms and functions of synchrony in interpersonal contexts using neurophysiological and behavioral measures.

Autapses promote synchronization in neuronal networks

Neurological disorders such as epileptic seizures are believed to be caused by neuronal synchrony. However, to ascertain the causal role of neuronal synchronization in such diseases through the traditional approach of electrophysiological data analysis remains a controversial, challenging, and outstanding problem. We offer an alternative principle to assess the physiological role of neuronal synchrony based on identifying structural anomalies in the underlying network and studying their impacts on the collective dynamics. In particular, we focus on autapses - time delayed self-feedback links that exist on a small fraction of neurons in the network, and investigate their impacts on network synchronization through a detailed stability analysis. Our main finding is that the proper placement of a small number of autapses in the network can promote synchronization significantly, providing the computational and theoretical bases for hypothesizing a high degree of synchrony in real neuronal networks with autapses. Our result that autapses, the shortest possible links in any network, can effectively modulate the collective dynamics provides also a viable strategy for optimal control of complex network dynamics at minimal cost.

Multistable states in a system of coupled phase oscillators with inertia

We investigate the generalized Kuramoto model of globally coupled oscillators with inertia, in which oscillators with positive coupling strength are conformists and oscillators with negative coupling strength are contrarians. We consider the correlation between the coupling strengths of oscillators and the distributions of natural frequencies. Two different types of correlations are studied. It is shown that the model supports multistable synchronized states such as different types of travelling wave states, π state and another type of nonstationary state: an oscillating π state. The phase distribution oscillates in a confined region and the phase difference between conformists and contrarians oscillates around π periodically in the oscillating π state. The different types of travelling wave state may be characterized by the speed of travelling wave and the effective frequencies of oscillators. Finally, the bifurcation diagrams of the model in the parameter space are presented.

Detecting and characterizing high-frequency oscillations in epilepsy: a case study of big data analysis

We develop a framework to uncover and analyse dynamical anomalies from massive, nonlinear and non-stationary time series data. The framework consists of three steps: preprocessing of massive datasets to eliminate erroneous data segments, application of the empirical mode decomposition and Hilbert transform paradigm to obtain the fundamental components embedded in the time series at distinct time scales, and statistical/scaling analysis of the components. As a case study, we apply our framework to detecting and characterizing high-frequency oscillations (HFOs) from a big database of rat electroencephalogram recordings. We find a striking phenomenon: HFOs exhibit on-off intermittency that can be quantified by algebraic scaling laws. Our framework can be generalized to big data-related problems in other fields such as large-scale sensor data and seismic data analysis.

A phase-synchronization and random-matrix based approach to multichannel time-series analysis with application to epilepsy

We present a general method to analyze multichannel time series that are becoming increasingly common in many areas of science and engineering. Of particular interest is the degree of synchrony among various channels, motivated by the recognition that characterization of synchrony in a system consisting of many interacting components can provide insights into its fundamental dynamics. Often such a system is complex, high-dimensional, nonlinear, nonstationary, and noisy, rendering unlikely complete synchronization in which the dynamical variables from individual components approach each other asymptotically. Nonetheless, a weaker type of synchrony that lasts for a finite amount of time, namely, phase synchronization, can be expected. Our idea is to calculate the average phase-synchronization times from all available pairs of channels and then to construct a matrix. Due to nonlinearity and stochasticity, the matrix is effectively random. Moreover, since the diagonal elements of the matrix can be arbitrarily large, the matrix can be singular. To overcome this difficulty, we develop a random-matrix based criterion for proper choosing of the diagonal matrix elements. Monitoring of the eigenvalues and the determinant provides a powerful way to assess changes in synchrony. The method is tested using a prototype nonstationary noisy dynamical system, electroencephalogram (scalp) data from absence seizures for which enhanced cortico-thalamic synchrony is presumed, and electrocorticogram (intracranial) data from subjects having partial seizures with secondary generalization for which enhanced local synchrony is similarly presumed.

Effects of adaptive coupling on stochastic resonance of small-world networks

The phenomenon of stochastic resonance in networks with small-world connectivity is investigated when the coupling strength is adaptive. The effects of the fixed and adaptive couplings on stochastic resonance of the system are discussed. It is found that the resonance is a monotonically increasing function of the adaptive coupling strength, while there is a peak when the coupling strength is fixed. The resonance for the adaptive coupling can reach a much larger value than that for fixed coupling.

Amplitude-phase synchronization at the onset of permanent spatiotemporal chaos

Amplitude and phase synchronization due to multiscale interactions in chaotic saddles at the onset of permanent spatiotemporal chaos is analyzed using the Fourier-Lyapunov representation. By computing the power-phase spectral entropy and the time-averaged power-phase spectra, we show that the laminar (bursty) states in the on-off spatiotemporal intermittency correspond, respectively, to the nonattracting coherent structures with higher (lower) degrees of amplitude-phase synchronization across spatial scales.

Automated synchrogram analysis applied to heartbeat and reconstructed respiration

Phase synchronization between two weakly coupled oscillators has been studied in chaotic systems for a long time. However, it is difficult to unambiguously detect such synchronization in experimental data from complex physiological systems. In this paper we review our study of phase synchronization between heartbeat and respiration in 150 healthy subjects during sleep using an automated procedure for screening the synchrograms. We found that this synchronization is significantly enhanced during non-rapid-eye-movement (non-REM) sleep (deep sleep and light sleep) and is reduced during REM sleep. In addition, we show that the respiration signal can be reconstructed from the heartbeat recordings in many subjects. Our reconstruction procedure, which works particularly well during non-REM sleep, allows the detection of cardiorespiratory synchronization even if only heartbeat intervals were recorded.

Model-based detector and extraction of weak signal frequencies from chaotic data

Detecting a weak signal from chaotic time series is of general interest in science and engineering. In this work we introduce and investigate a signal detection algorithm for which chaos theory, nonlinear dynamical reconstruction techniques, neural networks, and time-frequency analysis are put together in a synergistic manner. By applying the scheme to numerical simulation and different experimental measurement data sets (Henon map, chaotic circuit, and NH(3) laser data sets), we demonstrate that weak signals hidden beneath the noise floor can be detected by using a model-based detector. Particularly, the signal frequencies can be extracted accurately in the time-frequency space. By comparing the model-based method with the standard denoising wavelet technique as well as supervised principal components analysis detector, we further show that the nonlinear dynamics and neural network-based approach performs better in extracting frequencies of weak signals hidden in chaotic time series.

Experimental observation of synchronous competition in the Chua system

We present experimental results for two different sinusoidal functions competing to phase synchronize a Chua oscillator. Our approach involves real-time observation of the synchronization process. It shows that depending on the amplitude and frequency values of the two sinusoidal functions, the Chua oscillator can stay phase synchronized to one or the other of the inputs all the time or can alternate synchronous states between them. Numerical simulations show good agreement with the experimental observations.

Characterization of synchrony with applications to epileptic brain signals

Measurement of synchrony in networks of complex or high-dimensional, nonstationary, and noisy systems such as the mammalian brain is technically difficult. We present a general method to analyze synchrony from multichannel time series. The idea is to calculate the phase-synchronization times and to construct a matrix. We develop a random-matrix-based criterion for proper choosing of the diagonal matrix elements. Monitoring of the eigenvalues and the determinant provides an effective way to assess changes in synchrony. The method is tested using a prototype nonstationary dynamical system, electroencephalogram (scalp) data from absence seizures for which enhanced synchrony is presumed, and electrocorticogram (intracranial) data from subjects having partial seizures with secondary generalization.

Experimental evidence for phase synchronization transitions in the human cardiorespiratory system

Transitions in the dynamics of complex systems can be characterized by changes in the synchronization behavior of their components. Taking the human cardiorespiratory system as an example and using an automated procedure for screening the synchrograms of 112 healthy subjects we study the frequency and the distribution of synchronization episodes under different physiological conditions that occur during sleep. We find that phase synchronization between heartbeat and breathing is significantly enhanced during non-rapid-eye-movement (non-REM) sleep (deep sleep and light sleep) and reduced during REM sleep. Our results suggest that the synchronization is mainly due to a weak influence of the breathing oscillator upon the heartbeat oscillator, which is disturbed in the presence of long-term correlated noise, superimposed by the activity of higher brain regions during REM sleep.

Quantum stochastic synchronization

We study, within the spin-boson dynamics, the synchronization of a quantum tunneling system with an external, time-periodic driving signal. As a main result, we find that at a sufficiently large system-bath coupling strength (i.e., for a friction strength alpha > 1) the thermal noise plays a constructive role in yielding forced synchronization. This noise-induced synchronization can occur when the driving frequency is larger than the zero-temperature tunneling rate. As an application evidencing the effect, we consider the charge transfer dynamics in molecular complexes.