High expectations on phase locking: Better quantifying the concentration of circular data.

The degree to which unimodal circular data are concentrated around the mean direction can be quantified using the mean resultant length, a measure known under many alternative names, such as the phase locking value or the Kuramoto order parameter. For maximal concentration, achieved when all of the data take the same value, the mean resultant length attains its upper bound of one. However, for a random sample drawn from the circular uniform distribution, the expected value of the mean resultant length achieves its lower bound of zero only as the sample size tends to infinity. Moreover, as the expected value of the mean resultant length depends on the sample size, bias is induced when comparing the mean resultant lengths of samples of different sizes. In order to ameliorate this problem, here, we introduce a re-normalized version of the mean resultant length. Regardless of the sample size, the re-normalized measure has an expected value that is essentially zero for a random sample from the circular uniform distribution, takes intermediate values for partially concentrated unimodal data, and attains its upper bound of one for maximal concentration. The re-normalized measure retains the simplicity of the original mean resultant length and is, therefore, easy to implement and compute. We illustrate the relevance and effectiveness of the proposed re-normalized measure for mathematical models and electroencephalographic recordings of an epileptic seizure.

Chimera states in multiplex networks: Chameleon-like across-layer synchronization.

Different across-layer synchronization types of chimera states in multilayer networks have been discovered recently. We investigate possible relations between them, for example, if the onset of some synchronization type implies the onset of some other type. For this purpose, we use a two-layer network with multiplex inter-layer coupling. Each layer consists of a ring of non-locally coupled phase oscillators. While oscillators in each layer are identical, the layers are made non-identical by introducing mismatches in the oscillators' mean frequencies and phase lag parameters of the intra-layer coupling. We use different metrics to quantify the degree of various across-layer synchronization types. These include phase-locking between individual interacting oscillators, amplitude and phase synchronization between the order parameters of each layer, generalized synchronization between the driver and response layer, and the alignment of the incoherent oscillator groups' position on the two rings. For positive phase lag parameter mismatches, we get a cascaded onset of synchronization upon a gradual increase of the inter-layer coupling strength. For example, the two order parameters show phase synchronization before any of the interacting oscillator pairs does. For negative mismatches, most synchronization types have their onset in a narrow range of the coupling strength. Weaker couplings can destabilize chimera states in the response layer toward an almost fully coherent or fully incoherent motion. Finally, in the absence of a phase lag mismatch, sufficient coupling turns the response dynamics into a replica of the driver dynamics with the phases of all oscillators shifted by a constant lag.

Seizure forecasting: Where do we stand?

A lot of mileage has been made recently on the long and winding road toward seizure forecasting. Here we briefly review some selected milestones passed along the way, which were discussed at the International Conference for Technology and Analysis of Seizures-ICTALS 2022-convened at the University of Bern, Switzerland. Major impetus was gained recently from wearable and implantable devices that record not only electroencephalography, but also data on motor behavior, acoustic signals, and various signals of the autonomic nervous system. This multimodal monitoring can be performed for ultralong timescales covering months or years. Accordingly, features and metrics extracted from these data now assess seizure dynamics with a greater degree of completeness. Most prominently, this has allowed the confirmation of the long-suspected cyclical nature of interictal epileptiform activity, seizure risk, and seizures. The timescales cover daily, multi-day, and yearly cycles. Progress has also been fueled by approaches originating from the interdisciplinary field of network science. Considering epilepsy as a large-scale network disorder yielded novel perspectives on the pre-ictal dynamics of the evolving epileptic brain. In addition to discrete predictions that a seizure will take place in a specified prediction horizon, the community broadened the scope to probabilistic forecasts of a seizure risk evolving continuously in time. This shift of gears triggered the incorporation of additional metrics to quantify the performance of forecasting algorithms, which should be compared to the chance performance of constrained stochastic null models. An imminent task of utmost importance is to find optimal ways to communicate the output of seizure-forecasting algorithms to patients, caretakers, and clinicians, so that they can have socioeconomic impact and improve patients' well-being.

Phase irregularity: A conceptually simple and efficient approach to characterize electroencephalographic recordings from epilepsy patients.

The severe neurological disorder epilepsy affects almost 1% of the world population. For patients who suffer from pharmacoresistant focal-onset epilepsy, electroencephalographic (EEG) recordings are essential for the localization of the brain area where seizures start. Apart from the visual inspection of the recordings, quantitative EEG signal analysis techniques proved to be useful for this purpose. Among other features, regularity versus irregularity and phase coherence versus phase independence allowed characterizing brain dynamics from the measured EEG signals. Can phase irregularities also characterize brain dynamics? To address this question, we use the univariate coefficient of phase velocity variation, defined as the ratio of phase velocity standard deviation and the mean phase velocity. Beyond that, as a bivariate measure we use the classical mean phase coherence to quantify the degree of phase locking. All phase-based measures are combined with surrogates to test null hypotheses about the dynamics underlying the signals. In the first part of our analysis, we use the Rössler model system to study our approach under controlled conditions. In the second part, we use the Bern-Barcelona EEG database which consists of focal and nonfocal signals extracted from seizure-free recordings. Focal signals are recorded from brain areas where the first seizure EEG signal changes can be detected, and nonfocal signals are recorded from areas that are not involved in the seizure at its onset. Our results show that focal signals have less phase variability and more phase coherence than nonfocal signals. Once combined with surrogates, the mean phase velocity proved to have the highest discriminative power between focal and nonfocal signals. In conclusion, conceptually simple and easy to compute phase-based measures can help to detect features induced by epilepsy from EEG signals. This holds not only for the classical mean phase coherence but even more so for univariate measures of phase irregularity.

What Models and Tools can Contribute to a Better Understanding of Brain Activity?

Despite impressive scientific advances in understanding the structure and function of the human brain, big challenges remain. A deep understanding of healthy and aberrant brain activity at a wide range of temporal and spatial scales is needed. Here we discuss, from an interdisciplinary network perspective, the advancements in physical and mathematical modeling as well as in data analysis techniques that, in our opinion, have potential to further advance our understanding of brain structure and function.

Chimeras confined by fractal boundaries in the complex plane.

Complex-valued quadratic maps either converge to fixed points, enter into periodic cycles, show aperiodic behavior, or diverge to infinity. Which of these scenarios takes place depends on the map's complex-valued parameter c and the initial conditions. The Mandelbrot set is defined by the set of c values for which the map remains bounded when initiated at the origin of the complex plane. In this study, we analyze the dynamics of a coupled network of two pairs of two quadratic maps in dependence on the parameter c. Across the four maps, c is kept the same whereby the maps are identical. In analogy to the behavior of individual maps, the network iterates either diverge to infinity or remain bounded. The bounded solutions settle into different stable states, including full synchronization and desynchronization of all maps. Furthermore, symmetric partially synchronized states of within-pair synchronization and across-pair synchronization as well as a symmetry broken chimera state are found. The boundaries between bounded and divergent solutions in the domain of c are fractals showing a rich variety of intriguingly esthetic patterns. Moreover, the set of bounded solutions is divided into countless subsets throughout all length scales in the complex plane. Each individual subset contains only one state of synchronization and is enclosed within fractal boundaries by c values leading to divergence.

Seizure Cycles in Focal Epilepsy.

Importance: Focal epilepsy is characterized by the cyclical recurrence of seizures, but, to our knowledge, the prevalence and patterns of seizure cycles are unknown. Objective: To establish the prevalence, strength, and temporal patterns of seizure cycles over timescales of hours to years. Design, Setting, and Participants: This retrospective cohort study analyzed data from continuous intracranial electroencephalography (cEEG) and seizure diaries collected between January 19, 2004, and May 18, 2018, with durations up to 10 years. A total of 222 adults with medically refractory focal epilepsy were selected from 256 total participants in a clinical trial of an implanted responsive neurostimulation device. Selection was based on availability of cEEG and/or self-reports of disabling seizures. Exposures: Antiseizure medications and responsive neurostimulation, based on clinical indications. Main Outcomes and Measures: Measures involved (1) self-reported daily seizure counts, (2) cEEG-based hourly counts of electrographic seizures, and (3) detections of interictal epileptiform activity (IEA), which fluctuates in daily (circadian) and multiday (multidien) cycles. Outcomes involved descriptive characteristics of cycles of IEA and seizures: (1) prevalence, defined as the percentage of patients with a given type of seizure cycle; (2) strength, defined as the degree of consistency with which seizures occur at certain phases of an underlying cycle, measured as the phase-locking value (PLV); and (3) seizure chronotypes, defined as patterns in seizure timing evident at the group level. Results: Of the 222 participants, 112 (50%) were male, and the median age was 35 years (range, 18-66 years). The prevalence of circannual (approximately 1 year) seizure cycles was 12% (24 of 194), the prevalence of multidien (approximately weekly to approximately monthly) seizure cycles was 60% (112 of 186), and the prevalence of circadian (approximately 24 hours) seizure cycles was 89% (76 of 85). Strengths of circadian (mean [SD] PLV, 0.34 [0.18]) and multidien (mean [SD] PLV, 0.34 [0.17]) seizure cycles were comparable, whereas circannual seizure cycles were weaker (mean [SD] PLV, 0.17 [0.10]). Across individuals, circadian seizure cycles showed 5 peaks: morning, mid-afternoon, evening, early night, and late night. Multidien cycles of IEA showed peak periodicities centered around 7, 15, 20, and 30 days. Independent of multidien period length, self-reported and electrographic seizures consistently occurred during the days-long rising phase of multidien cycles of IEA. Conclusions and Relevance: Findings in this large cohort establish the high prevalence of plural seizure cycles and help explain the natural variability in seizure timing. The results have the potential to inform the scheduling of diagnostic studies, the delivery of time-varying therapies, and the design of clinical trials in epilepsy.

Remote pacemaker control of chimera states in multilayer networks of neurons.

Networks of coupled nonlinear oscillators allow for the formation of nontrivial partially synchronized spatiotemporal patterns, such as chimera states, in which there are coexisting coherent (synchronized) and incoherent (desynchronized) domains. These complementary domains form spontaneously, and it is impossible to predict where the synchronized group will be positioned within the network. Therefore, possible ways to control the spatial position of the coherent and incoherent groups forming the chimera states are of high current interest. In this work we investigate how to control chimera patterns in multiplex networks of FitzHugh-Nagumo neurons, and in particular we want to prove that it is possible to remotely control chimera states exploiting the multiplex structure. We introduce a pacemaker oscillator within the network: this is an oscillator that does not receive input from the rest of the network but is sending out information to its neighbors. The pacemakers can be positioned in one or both layers. Their presence breaks the spatial symmetry of the layer in which they are introduced and allows us to control the position of the incoherent domain. We demonstrate how the remote control is possible for both uni- and bidirectional coupling between the layers. Furthermore we show which are the limitations of our control mechanisms when it is generalized from single-layer to multilayer networks.

Seizure Onset Zone Lateralization Using a Non-linear Analysis of Micro vs. Macro Electroencephalographic Recordings During Seizure-Free Stages of the Sleep-Wake Cycle From Epilepsy Patients.

The application of non-linear signal analysis techniques to biomedical data is key to improve our knowledge about complex physiological and pathological processes. In particular, the use of non-linear techniques to study electroencephalographic (EEG) recordings can provide an advanced characterization of brain dynamics. In epilepsy these dynamics are altered at different spatial scales of neuronal organization. We therefore apply non-linear signal analysis to EEG recordings from epilepsy patients derived with intracranial hybrid electrodes, which are composed of classical macro contacts and micro wires. Thereby, these electrodes record EEG at two different spatial scales. Our aim is to test the degree to which the analysis of the EEG recorded at these different scales allows us to characterize the neuronal dynamics affected by epilepsy. For this purpose, we retrospectively analyzed long-term recordings performed during five nights in three patients during which no seizures took place. As a benchmark we used the accuracy with which this analysis allows determining the hemisphere that contains the seizure onset zone, which is the brain area where clinical seizures originate. We applied the surrogate-corrected non-linear predictability score (ψ), a non-linear signal analysis technique which was shown previously to be useful for the lateralization of the seizure onset zone from classical intracranial EEG macro contact recordings. Higher values of ψ were found predominantly for signals recorded from the hemisphere containing the seizure onset zone as compared to signals recorded from the opposite hemisphere. These differences were found not only for the EEG signals recorded with macro contacts, but also for those recorded with micro wires. In conclusion, the information obtained from the analysis of classical macro EEG contacts can be complemented by the one of micro wire EEG recordings. This combined approach may therefore help to further improve the degree to which quantitative EEG analysis can contribute to the diagnostics in epilepsy patients.

Two populations of coupled quadratic maps exhibit a plentitude of symmetric and symmetry broken dynamics.

We numerically study a network of two identical populations of identical real-valued quadratic maps. Upon variation of the coupling strengths within and across populations, the network exhibits a rich variety of distinct dynamics. The maps in individual populations can be synchronized or desynchronized. Their temporal evolution can be periodic or aperiodic. Furthermore, one can find blends of synchronized with desynchronized states and periodic with aperiodic motions. We show symmetric patterns for which both populations have the same type of dynamics as well as chimera states of a broken symmetry. The network can furthermore show multistability by settling to distinct dynamics for different realizations of random initial conditions or by switching intermittently between distinct dynamics for the same realization. We conclude that our system of two populations of a particularly simple map is the most simple system that can show this highly diverse and complex behavior, which includes but is not limited to chimera states. As an outlook to future studies, we explore the stability of two populations of quadratic maps with a complex-valued control parameter. We show that bounded and diverging dynamics are separated by fractal boundaries in the complex plane of this control parameter.

Controlling chimera states via minimal coupling modification.

We propose a method to control chimera states in a ring-shaped network of nonlocally coupled phase oscillators. This method acts exclusively on the network's connectivity. Using the idea of a pacemaker oscillator, we investigate which is the minimal action needed to control chimeras. We implement the pacemaker choosing one oscillator and making its links unidirectional. Our results show that a pacemaker induces chimeras for parameters and initial conditions for which they do not form spontaneously. Furthermore, the pacemaker attracts the incoherent part of the chimera state, thus controlling its position. Beyond that, we find that these control effects can be achieved with modifications of the network's connectivity that are less invasive than a pacemaker, namely, the minimal action of just modifying the strength of one connection allows one to control chimeras.

Inferring directed networks using a rank-based connectivity measure.

Inferring the topology of a network using the knowledge of the signals of each of the interacting units is key to understanding real-world systems. One way to address this problem is using data-driven methods like cross-correlation or mutual information. However, these measures lack the ability to distinguish the direction of coupling. Here, we use a rank-based nonlinear interdependence measure originally developed for pairs of signals. This measure not only allows one to measure the strength but also the direction of the coupling. Our results for a system of coupled Lorenz dynamics show that we are able to consistently infer the underlying network for a subrange of the coupling strength and link density. Furthermore, we report that the addition of dynamical noise can benefit the reconstruction. Finally, we show an application to multichannel electroencephalographic recordings from an epilepsy patient.

Mean field phase synchronization between chimera states.

We study two-layer networks of identical phase oscillators. Each individual layer is a ring network for which a non-local intra-layer coupling leads to the formation of a chimera state. The number of oscillators and their natural frequencies is in general different across the layers. We couple the phases of individual oscillators in one layer to the phase of the mean field of the other layer. This coupling from the mean field to individual oscillators is done in both directions. For a sufficient strength of this inter-layer coupling, the phases of the mean fields lock across the two layers. In contrast, both layers continue to exhibit chimera states with no locking between the phases of individual oscillators across layers, and the two mean field amplitudes remain uncorrelated. Hence, the networks' mean fields show phase synchronization which is analogous to the one between low-dimensional chaotic oscillators. The required coupling strength to achieve this mean field phase synchronization increases with the mismatches in the network sizes and the oscillators' natural frequencies.

Robustness and versatility of a nonlinear interdependence method for directional coupling detection from spike trains.

The detection of directional couplings between dynamics based on measured spike trains is a crucial problem in the understanding of many different systems. In particular, in neuroscience it is important to assess the connectivity between neurons. One of the approaches that can estimate directional coupling from the analysis of point processes is the nonlinear interdependence measure L. Although its efficacy has already been demonstrated, it still needs to be tested under more challenging and realistic conditions prior to an application to real data. Thus, in this paper we use the Hindmarsh-Rose model system to test the method in the presence of noise and for different spiking regimes. We also examine the influence of different parameters and spike train distances. Our results show that the measure L is versatile and robust to various types of noise, and thus suitable for application to experimental data.

Generalized synchronization between chimera states.

Networks of coupled oscillators in chimera states are characterized by an intriguing interplay of synchronous and asynchronous motion. While chimera states were initially discovered in mathematical model systems, there is growing experimental and conceptual evidence that they manifest themselves also in natural and man-made networks. In real-world systems, however, synchronization and desynchronization are not only important within individual networks but also across different interacting networks. It is therefore essential to investigate if chimera states can be synchronized across networks. To address this open problem, we use the classical setting of ring networks of non-locally coupled identical phase oscillators. We apply diffusive drive-response couplings between pairs of such networks that individually show chimera states when there is no coupling between them. The drive and response networks are either identical or they differ by a variable mismatch in their phase lag parameters. In both cases, already for weak couplings, the coherent domain of the response network aligns its position to the one of the driver networks. For identical networks, a sufficiently strong coupling leads to identical synchronization between the drive and response. For non-identical networks, we use the auxiliary system approach to demonstrate that generalized synchronization is established instead. In this case, the response network continues to show a chimera dynamics which however remains distinct from the one of the driver. Hence, segregated synchronized and desynchronized domains in individual networks congregate in generalized synchronization across networks.

Coupling strength versus coupling impact in nonidentical bidirectionally coupled dynamics.

The understanding of interacting dynamics is important for the characterization of real-world networks. In general, real-world networks are heterogeneous in the sense that each node of the network is a dynamics with different properties. For coupled nonidentical dynamics symmetric interactions are not straightforwardly defined from the coupling strength values. Thus, a challenging issue is whether we can define a symmetric interaction in this asymmetric setting. To address this problem we introduce the notion of the coupling impact. The coupling impact considers not only the coupling strength but also the energy of the individual dynamics, which is conveyed via the coupling. To illustrate this concept, we follow a data-driven approach by analyzing signals from pairs of coupled model dynamics using two different connectivity measures. We find that the coupling impact, but not the coupling strength, correctly detects a symmetric interaction between pairs of coupled dynamics regardless of their degree of asymmetry. Therefore, this approach allows us to reveal the real impact that one dynamics has on the other and hence to define symmetric interactions in pairs of nonidentical dynamics.

Ictal time-irreversible intracranial EEG signals as markers of the epileptogenic zone.

OBJECTIVE: To show that time-irreversible EEG signals recorded with intracranial electrodes during seizures can serve as markers of the epileptogenic zone. METHODS: We use the recently developed method of mapping time series into directed horizontal graphs (dHVG). Each node of the dHVG represents a time point in the original intracranial EEG (iEEG) signal. Statistically significant differences between the distributions of the nodes' number of input and output connections are used to detect time-irreversible iEEG signals. RESULTS: In 31 of 32 seizure recordings we found time-irreversible iEEG signals. The maximally time-irreversible signals always occurred during seizures, with highest probability in the middle of the first seizure half. These signals spanned a large range of frequencies and amplitudes but were all characterized by saw-tooth like shaped components. Brain regions removed from patients who became post-surgically seizure-free generated significantly larger time-irreversibilities than regions removed from patients who still had seizures after surgery. CONCLUSIONS: Our results corroborate that ictal time-irreversible iEEG signals can indeed serve as markers of the epileptogenic zone and can be efficiently detected and quantified in a time-resolved manner by dHVG based methods. SIGNIFICANCE: Ictal time-irreversible EEG signals can help to improve pre-surgical evaluation in patients suffering from pharmaco-resistant epilepsies.

All together now: Analogies between chimera state collapses and epileptic seizures.

Conceptually and structurally simple mathematical models of coupled oscillator networks can show a rich variety of complex dynamics, providing fundamental insights into many real-world phenomena. A recent and not yet fully understood example is the collapse of coexisting synchronous and asynchronous oscillations into a globally synchronous motion found in networks of identical oscillators. Here we show that this sudden collapse is promoted by a further decrease of synchronization, rather than by critically high synchronization. This strikingly counterintuitive mechanism can be found also in nature, as we demonstrate on epileptic seizures in humans. Analyzing spatiotemporal correlation profiles derived from intracranial electroencephalographic recordings (EEG) of seizures in epilepsy patients, we found a pronounced decrease of correlation at the seizure onsets. Applying our findings in a closed-loop control scheme to models of coupled oscillators in chimera states, we succeed in both provoking and preventing outbreaks of global synchronization. Our findings not only advance the understanding of networks of coupled dynamics but can open new ways to control them, thus offering a vast range of potential new applications.

Resected Brain Tissue, Seizure Onset Zone and Quantitative EEG Measures: Towards Prediction of Post-Surgical Seizure Control.

BACKGROUND: Epilepsy surgery is a potentially curative treatment option for pharmacoresistent patients. If non-invasive methods alone do not allow to delineate the epileptogenic brain areas the surgical candidates undergo long-term monitoring with intracranial EEG. Visual EEG analysis is then used to identify the seizure onset zone for targeted resection as a standard procedure. METHODS: Despite of its great potential to assess the epileptogenicty of brain tissue, quantitative EEG analysis has not yet found its way into routine clinical practice. To demonstrate that quantitative EEG may yield clinically highly relevant information we retrospectively investigated how post-operative seizure control is associated with four selected EEG measures evaluated in the resected brain tissue and the seizure onset zone. Importantly, the exact spatial location of the intracranial electrodes was determined by coregistration of pre-operative MRI and post-implantation CT and coregistration with post-resection MRI was used to delineate the extent of tissue resection. Using data-driven thresholding, quantitative EEG results were separated into normally contributing and salient channels. RESULTS: In patients with favorable post-surgical seizure control a significantly larger fraction of salient channels in three of the four quantitative EEG measures was resected than in patients with unfavorable outcome in terms of seizure control (median over the whole peri-ictal recordings). The same statistics revealed no association with post-operative seizure control when EEG channels contributing to the seizure onset zone were studied. CONCLUSIONS: We conclude that quantitative EEG measures provide clinically relevant and objective markers of target tissue, which may be used to optimize epilepsy surgery. The finding that differentiation between favorable and unfavorable outcome was better for the fraction of salient values in the resected brain tissue than in the seizure onset zone is consistent with growing evidence that spatially extended networks might be more relevant for seizure generation, evolution and termination than a single highly localized brain region (i.e. a "focus") where seizures start.