Source author record

Per Sebastian Skardal

Per Sebastian Skardal appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2026arXiv

Hodge Decomposition Guides the Optimization of Synchronization over Simplicial Complexes

Despite growing interest in synchronization dynamics over "higher-order" network models, optimization theory for such systems is limited. Here, we study a family of Kuramoto models inspired by algebraic topology in which oscillators are coupled over simplicial complexes (SCs) using their associated Hodge Laplacian matrices. We optimize such systems by extending the synchrony alignment function -- an optimization framework for synchronizing graph-coupled heterogeneous oscillators. Computational experiments are given to illustrate how this approach can effectively solve a variety of combinatorial problems including the joint optimization of projected synchronization dynamics onto lower- and upper-dimensional simplices within SCs. We also investigate the role of SC homology and develop bifurcation theory to characterize the extent to which optimal solutions are contained within (or spread across) the three Hodge subspaces. Our work extends optimization theory to the setting of higher-order networks, provides practical algorithms for Hodge-Laplacian-related dynamics including (but not limited to) Kuramoto oscillators, and paves the way for an emerging field that interfaces algebraic topology, combinatorial optimization, and dynamical systems.

preprint2022arXiv

Tiered synchronization in coupled oscillator populations with interaction delays and higher-order interactions

We study synchronization in large populations of coupled phase oscillators with time-delays, higher order interactions. With each of these effects individually giving rise to bistabiltiy between incoherence and synchronization via a subcriticality at the onset of synchronization and the development of a saddle node, we find that their combination yields another mechanism behind bistability, where supercriticality at onset may be maintained and instead the formation of two saddle nodes creates tiered synchronization, i.e., bistability between a weakly synchronized state and a strongly synchronized state. We demonstrate these findings by first deriving the low dimensional dynamics of the system and examining the system bifurcations using a stability and steady-state analysis.

preprint2021arXiv

Balanced Hodge Laplacians Optimize Consensus Dynamics over Simplicial Complexes

Despite the vast literature on network dynamics, we still lack basic insights into dynamics on higher-order structures (e.g., edges, triangles, and more generally, $k$-dimensional "simplices") and how they are influenced through higher-order interactions. A prime example lies in neuroscience where groups of neurons (not individual ones) may provide the building blocks for neurocomputation. Here, we study consensus dynamics on edges in simplicial complexes using a type of Laplacian matrix called a Hodge Laplacian, which we generalize to allow higher- and lower-order interactions to have different strengths. Using techniques from algebraic topology, we study how collective dynamics converge to a low-dimensional subspace that corresponds to the homology space of the simplicial complex. We use the Hodge decomposition to show that higher- and lower-order interactions can be optimally balanced to maximally accelerate convergence, and that this optimum coincides with a balancing of dynamics on the curl and gradient subspaces. We additionally explore the effects of network topology, finding that consensus over edges is accelerated when 2-simplices are well dispersed, as opposed to clustered together.

preprint2021arXiv

Chaos in Nonlinear Random Walks with Non-Monotonic Transition Probabilities

Random walks serve as important tools for studying complex network structures, yet their dynamics in cases where transition probabilities are not static remain under explored and poorly understood. Here we study nonlinear random walks that occur when transition probabilities depend on the state of the system. We show that when these transition probabilities are non-monotonic, i.e., are not uniformly biased towards the most densely or sparsely populated nodes, but rather direct random walkers with more nuance, chaotic dynamics emerge. Using multiple transition probability functions and a range of networks with different connectivity properties, we demonstrate that this phenomenon is generic. Thus, when such non-monotonic properties are key ingredients in nonlinear transport applications complicated and unpredictable behaviors may result.

preprint2021arXiv

Higher-order interactions improve optimal collective dynamics on networks

Collective behavior plays a key role in the function of a wide range of physical, biological, and neurological systems where empirical evidence has recently uncovered the prevalence of higher-order interactions, i.e., structures that represent interactions between more than just two individual units, in complex network structures. Here, we study the optimization of collective behavior in networks with higher-order interactions encoded in clique complexes. Our approach involves adapting the Synchrony Alignment Function framework to a new composite Laplacian matrix that encodes multi-order interactions including, e.g., both dyadic and triadic couplings. We show that as higher-order coupling interactions are equitably strengthened, so that overall coupling is conserved, the optimal collective behavior improves. We find that this phenomenon stems from the broadening of a composite Laplacian's eigenvalue spectrum, which improves the optimal collective behavior and widens the range of possible behaviors. Moreover, we find in constrained optimization scenarios that a nontrivial, ideal balance between the relative strengths of pair-wise and higher-order interactions leads to the strongest collective behavior supported by a network. This work provides insight into how systems balance interactions of different types to optimize or broaden their dynamical range of behavior, especially for self-regulating systems like the brain.

preprint2020arXiv

Bifurcation analysis and structural stability of simplicial oscillator populations

We present an analytical description for the collective dynamics of oscillator ensembles with higher-order coupling encoded by simplicial structure, which serves as an illustrative and insightful paradigm for brain function and information storage. The novel dynamics of the system, including abrupt desynchronization and multistability, are rigorously characterized and the critical points that correspond to a continuum of first-order phase transitions are found to satisfy universal scaling properties. More importantly, the underlying bifurcation mechanism giving rise to multiple clusters with arbitrary ensemble size is characterized using a rigorous spectral analysis of the stable cluster states. As a consequence of $SO_2$ group symmetry, we show that the continuum of abrupt desynchronization transitions result from the instability of a collective mode under the nontrivial antisymmetric manifold in the high dimensional phase space.

preprint2020arXiv

Spectrum of extensive multiclusters in the Kuramoto model with higher-order interactions

Globally coupled ensembles of phase oscillators serve as useful tools for modeling synchronization and collective behavior in a variety of applications. As interest in the effects of simplicial interactions (i.e., non-additive, higher-order interactions between three or more units) continues to grow we study an extension of the Kuramoto model where oscillators are coupled via three-way interactions that exhibits novel dynamical properties including clustering, multistability, and abrupt desynchronization transitions. Here we provide a rigorous description of the stability of various multicluster states by studying their spectral properties in the thermodynamic limit. Not unlike the classical Kuramoto model, a natural frequency distribution with infinite support yields a population of drifting oscillators, which in turn guarantees that a portion of the spectrum is located on the imaginary axes, resulting in neutrally stable or unstable solutions. On the other hand, a natural frequency distribution with finite support allows for a fully phase-locked state, whose spectrum is real and may be linearly stable or unstable.

preprint2019arXiv

Quasi-periodic dynamics and a Neimark-Sacker bifurcation in nonlinear random walks on complex networks

We study the dynamics of nonlinear random walks on complex networks. We investigate the role and effect of directed network topologies on long-term dynamics. While a period-doubling bifurcation to alternating patterns occurs at a critical bias parameter value, we find that some directed structures give rise to a different kind of bifurcation that gives rise to quasi-periodic dynamics. This does not occur for all directed network structure, but only when the network structure is sufficiently directed. We find that the onset of quasi-periodic dynamics is the result of a Neimark-Sacker bifurcation, where a pair of complex-conjugate eigenvalues of the system Jacobian passes through the unit circle, destabilizing the stationary distribution with high-dimensional rotations. We investigate the nature of these bifurcations, study the onset of quasi-periodic dynamics as network structure is tuned to be more directed, and present an analytically tractable case of a four-neighbor ring.

preprint2019arXiv

Synchronization of Network-Coupled Oscillators with Uncertain Dynamics

Synchronization of network-coupled dynamical units is important to a variety of natural and engineered processes including circadian rhythms, cardiac function, neural processing, and power grids. Despite this ubiquity, it remains poorly understood how complex network structures and heterogeneous local dynamics combine to either promote or inhibit synchronization. Moreover, for most real-world applications it is impossible to obtain the exact specifications of the system, and there is a lack of theory for how uncertainty affects synchronization. We address this open problem by studying the Synchrony Alignment Function (SAF), which is an objective measure for the synchronization properties of a network of heterogeneous oscillators with given natural frequencies. We extend the SAF framework to analyze network-coupled oscillators with heterogeneous natural frequencies that are drawn as a multivariate random vector. Using probability theory for quadratic forms, we obtain expressions for the expectation and variance of the SAF for given network structures. We conclude with numerical experiments that illustrate how the incorporation of uncertainty yields a more robust theoretical framework for enhancing synchronization, and we provide new perspectives for why synchronization is generically promoted by network properties including degree-frequency correlations, link directedness, and link weight delocalization.

preprint2016arXiv

Collective frequency variation in network synchronization and reverse PageRank

A wide range of natural and engineered phenomena rely on large networks of interacting units to reach a dynamical consensus state where the system collectively operates. Here we study the dynamics of self-organizing systems and show that for generic directed networks the collective frequency of the ensemble is {\it not} the same as the mean of the individuals' natural frequencies. Specifically, we show that the collective frequency equals a weighted average of the natural frequencies, where the weights are given by an out-flow centrality measure that is equivalent to a reverse PageRank centrality. Our findings uncover an intricate dependence of the collective frequency on both the structural directedness and dynamical heterogeneity of the network, and also reveal an unexplored connection between synchronization and PageRank, which opens the possibility of applying PageRank optimization to synchronization. Finally, we demonstrate the presence of collective frequency variation in real-world networks by considering the UK and Scandinavian power grids.

preprint2016arXiv

On controlling networks of limit-cycle oscillators

The control of network-coupled nonlinear dynamical systems is an active area of research in the nonlinear science community. Coupled oscillator networks represent a particularly important family of nonlinear systems, with applications ranging from the power grid to cardiac excitation. Here we study the control of network-coupled limit cycle oscillators, extending previous work that focused on phase oscillators. Based on stabilizing a target fixed point, our method aims to attain complete frequency synchronization, i.e., consensus, by applying control to as few oscillators as possible. We develop two types of control. The first type directs oscillators towards to larger amplitudes, while the second does not. We present numerical examples of both control types and comment on the potential failures of the method.

preprint2016arXiv

Optimal synchronization of directed complex networks

We study optimal synchronization of networks of coupled phase oscillators. We extend previous theory for optimizing the synchronization properties of undirected networks to the important case of directed networks. We derive a generalized synchrony alignment function that encodes the interplay between network structure and the oscillators' natural frequencies and serves as an objective measure for the network's degree of synchronization. Using the generalized synchrony alignment function, we show that a network's synchronization properties can be systematically optimized. This framework also allows us to study the properties of synchrony-optimized networks, and in particular, investigate the role of directed network properties such as nodal in- and out-degrees. For instance, we find that in optimally rewired networks the heterogeneity of the in-degree distribution roughly matches the heterogeneity of the natural frequency distribution, but no such relationship emerges for out-degrees. We also observe that a network's synchronization properties are promoted by a strong correlation between the nodal in-degrees and the natural frequencies of oscillators, whereas the relationship between the nodal out-degrees and the natural frequencies has comparatively little effect. This result is supported by our theory, which indicates that synchronization is promoted by a strong alignment of the natural frequencies with the left singular vectors corresponding to the largest singular values of the Laplacian matrix.

preprint2016arXiv

Spectral properties of the hierarchical product of graphs

The hierarchical product of two graphs represents a natural way to build a larger graph out of two smaller graphs with less regular and therefore more heterogeneous structure than the Cartesian product. Here we study the eigenvalue spectrum of the adjacency matrix of the hierarchical product of two graphs. Introducing a coupling parameter describing the relative contribution of each of the two smaller graphs, we perform an asymptotic analysis for the full spectrum of eigenvalues of the adjacency matrix of the hierarchical product. Specifically, we derive the exact limit points for each eigenvalue in the limits of small and large coupling, as well as the leading-order relaxation to these values in terms of the eigenvalues and eigenvectors of the two smaller graphs. Given its central roll in the structural and dynamical properties of networks, we study in detail the Perron-Frobenius, or largest, eigenvalue. Finally, as an example application we use our theory to predict the epidemic threshold of the Susceptible-Infected-Susceptible model on a hierarchical product of two graphs.

preprint2015arXiv

Control of coupled oscillator networks with application to microgrid technologies

The control of complex systems and network-coupled dynamical systems is a topic of vital theoretical importance in mathematics and physics with a wide range of applications in engineering and various other sciences. Motivated by recent research into smart grid technologies we study here control of synchronization and consider the important case of networks of coupled phase oscillators with nonlinear interactions--a paradigmatic example that has guided our understanding of self-organization for decades. We develop a method for control based on identifying and stabilizing problematic oscillators, resulting in a stable spectrum of eigenvalues, and in turn a linearly stable synchronized state. Interestingly, the amount of control, i.e., number of oscillators, required to stabilize the network is primarily dictated by the coupling strength, dynamical heterogeneity, and mean degree of the network, and depends little on the structural heterogeneity of the network itself.

preprint2015arXiv

Erosion of synchronization in networks of coupled oscillators

We report erosion of synchronization in networks of coupled phase oscillators, a phenomenon where perfect phase synchronization is unattainable in steady-state, even in the limit of infinite coupling. An analysis reveals that the total erosion is separable into the product of terms characterizing coupling frustration and structural heterogeneity, both of which amplify erosion. The latter, however, can differ significantly from degree heterogeneity. Finally, we show that erosion is marked by the reorganization of oscillators according to their node degrees rather than their natural frequencies.

preprint2015arXiv

Erosion of synchronization: Coupling heterogeneity and network structure

We study the dynamics of network-coupled phase oscillators in the presence of coupling frustration. It was recently demonstrated that in heterogeneous network topologies, the presence of coupling frustration causes perfect phase synchronization to become unattainable even in the limit of infinite coupling strength. Here, we consider the important case of heterogeneous coupling functions and extend previous results by deriving analytical predictions for the total erosion of synchronization. Our analytical results are given in terms of basic quantities related to the network structure and coupling frustration. In addition to fully heterogeneous coupling, where each individual interaction is allowed to be distinct, we also consider partially heterogeneous coupling and homogeneous coupling in which the coupling functions are either unique to each oscillator or identical for all network interactions, respectively. We demonstrate the validity of our theory with numerical simulations of multiple network models, and highlight the interesting effects that various coupling choices and network models have on the total erosion of synchronization. Finally, we consider some special network structures with well-known spectral properties, which allows us to derive further analytical results.

preprint2015arXiv

Frequency assortativity can induce chaos in oscillator networks

We investigate the effect of preferentially connecting oscillators with similar frequency to each other in networks of coupled phase oscillators (i.e., frequency assortativity). Using the network Kuramoto model as an example, we find that frequency assortativity can induce chaos in the macroscopic dynamics. By applying a mean-field approximation in combination with the dimension reduction method of Ott and Antonsen, we show that the dynamics can be described by a low dimensional system of equations. We use the reduced system to characterize the macroscopic chaos using Lyapunov exponents, bifurcation diagrams, and time-delay embeddings. Finally, we show that the emergence of chaos stems from the formation of multiple groups of synchronized oscillators, i.e., meta-oscillators.

preprint2015arXiv

Nonlinear Transport in the Stochastic Standard Map

We study a stochastically-driven standard map. The addition of a noise term destroys the invariant manifolds that organize the phase space which allows for more widespread transport than in the noiseless case. Using appropriately defined hitting times to quantify the dynamics, we identify two qualitatively different classes of transport: linear and nonlinear. Linear transport is primarily driven by the stochasticity in the system, while nonlinear transport results from a combination of the nonlinear dynamics and stochasticity and provides a significant speed-up in transport.

preprint2014arXiv

Coexisting chaotic and multi-periodic dynamics in a model of cardiac alternans

The spatiotemporal dynamics of cardiac tissue is an active area of research for biologists, physicists, and mathematicians. Of particular interest is the study of period-doubling bifurcations and chaos due to their link with cardiac arrhythmogenesis. In this paper we study the spatiotemporal dynamics of a recently developed model for calcium-driven alternans in a one dimensional cable of tissue. In particular, we observe in the cable coexistence of regions with chaotic and multi-periodic dynamics over wide ranges of parameters. We study these dynamics using global and local Lyapunov exponents and spatial trajectory correlations. Interestingly, near nodes -- or phase reversals -- low-periodic dynamics prevail, while away from the nodes the dynamics tend to be higher-periodic and eventually chaotic. Finally, we show that similar coexisting multi-periodic and chaotic dynamics can also be observed in a detailed ionic model.

preprint2014arXiv

Disorder induces explosive synchronization

We study explosive synchronization, a phenomenon characterized by first-order phase transitions between incoherent and synchronized states in networks of coupled oscillators. While explosive synchronization has been the subject of many recent studies, in each case strong conditions on either the heterogeneity of the network, its link weights, or its initial construction are imposed to engineer a first-order phase transition. This raises the question of how robust explosive synchronization is in view of more realistic structural and dynamical properties. Here we show that explosive synchronization can be induced in mildly heterogeneous networks by the addition of quenched disorder to the oscillators' frequencies, demonstrating that it is not only robust to, but moreover promoted by, this natural mechanism. We support these findings with numerical and analytical results, presenting simulations of a real neural network as well as a self-consistency theory used to study synthetic networks.

preprint2014arXiv

Optimal synchronization of complex networks

We study optimal synchronization in networks of heterogeneous phase oscillators. Our main result is the derivation of a synchrony alignment function that encodes the interplay between network structure and oscillators' frequencies and can be readily optimized. We highlight its utility in two general problems: constrained frequency allocation and network design. In general, we find that synchronization is promoted by strong alignments between frequencies and the dominant Laplacian eigenvectors, as well as a matching between the heterogeneity of frequencies and network structure.

preprint2014arXiv

Spatiotemporal Dynamics of Calcium-Driven Cardiac Alternans

We investigate the dynamics of spatially discordant alternans (SDA) driven by an instability of intracellular calcium cycling using both amplitude equations [P. S. Skardal, A. Karma, and J. G. Restrepo, Phys. Rev. Lett. 108, 108103 (2012)] and ionic model simulations. We find that, close to the alternans bifurcation, SDA is manifested as a smooth wavy modulation of the amplitudes of both repolarization and calcium transient (CaT) alternans, similarly to the well-studied case of voltage-driven alternans. In contrast, further away from the bifurcation, the amplitude of CaT alternans jumps discontinuously at the nodes separating out-of-phase regions, while the amplitude of repolarization alternans remains smooth. We show that node motion of discontinuous SDA patterns is strongly hysteretic even in homogeneous tissue due to the novel phenomenon of "unidirectional pinning": node movement can only be induced towards, but not away from, the pacing site in response to a change of pacing rate or physiological parameter. In addition, we show that the wavelength of discontinuous SDA patterns scales linearly with the conduction velocity restitution length scale, in contrast to the wavelength of smooth patterns that scales sub-linearly with this length scale. Those results are also shown to be robust against cell-to-cell fluctuations owing to the property that unidirectional node motion collapses multiple jumps accumulating in nodal regions into a single jump. Amplitude equation predictions are in good overall agreement with ionic model simulations. Finally, we briefly discuss physiological implications of our findings. In particular, we suggest that due to the tendency of conduction blocks to form near nodes, the presence of unidirectional pinning makes calcium-driven alternans potentially more arrhythmogenic than voltage-driven alternans.

preprint2012arXiv

Complex macroscopic behavior in systems of phase oscillators with adaptive coupling

Using recent dimensionality reduction techniques in large systems of coupled phase oscillators exhibiting bistability, we analyze complex macroscopic behavior arising when the coupling between oscillators is allowed to evolve slowly as a function of either macroscopic or local system properties. For example, we observe macroscopic excitability and intermittent synchrony in a system of time-delayed Kuramoto oscillators with Hebbian and anti-Hebbian learning. We demonstrate the robustness of our findings by considering systems with increasing complexity, including time-delayed oscillators with adaptive network structure and community interaction, as well as a system with bimodally distributed frequencies.

preprint2012arXiv

Effects of degree-frequency correlations on network synchronization: universality and full phase-locking

We introduce a model to study the effect of degree-frequency correlations on synchronization in networks of coupled oscillators. Analyzing this model, we find several remarkable characteristics. We find a stationary synchronized state that is (i) universal, i.e., the degree of synchrony, as measured by a global order parameter, is independent of network topology, and (ii) fully phase-locked, i.e., all oscillators become simultaneously phase-locked despite having different natural frequencies. This state separates qualitatively different behaviors for two other classes of correlations where, respectively, slow and fast oscillators can remain unsynchronized. We close by presenting analysis of the dynamics under arbitrary degree-frequency correlations.

preprint2012arXiv

Synchronization of Kuramoto oscillators in networks of networks

We study synchronization of Kuramoto oscillators in strongly modular networks in which the structure of the network inside each community is averaged. We find that the dynamics of the interacting communities can be described as an ensemble of coupled planar oscillators. In the limit of a large number of communities, we find a low dimensional description of the level of synchronization between the communities. In this limit, we describe bifurcations between incoherence, local synchrony, and global synchrony. We compare the predictions of this simplified model with simulations of heterogeneous networks in which the internal structure of each community is preserved and find excellent agreement. Finally, we investigate synchronization in networks where several layers of communities within communities may be present.

preprint2012arXiv

Unidirectional Pinning and Hysteresis of Spatially Discordant Alternans in Cardiac Tissue

Spatially discordant alternans is a widely observed pattern of voltage and calcium signals in cardiac tissue that can precipitate lethal cardiac arrhythmia. Using spatially coupled iterative maps of the beat-to-beat dynamics, we explore this pattern's dynamics in the regime of a calcium-dominated period-doubling instability at the single cell level. We find a novel nonlinear bifurcation associated with the formation of a discontinuous jump in the amplitude of calcium alternans at nodal lines separating discordant regions. We show that this jump unidirectionally pins nodal lines by preventing their motion away from the pacing site following a pacing rate decrease, but permitting motion towards this site following a rate increase. This unidirectional pinning leads to strongly history-dependent nodal line motion that is strongly arrhythmogenic.

preprint2011arXiv

Cluster Synchrony in Systems of Coupled Phase Oscillators with Higher-Order Coupling

We study the phenomenon of cluster synchrony that occurs in ensembles of coupled phase oscillators when higher-order modes dominate the coupling between oscillators. For the first time, we develop a complete analytic description of the dynamics in the limit of a large number of oscillators and use it to quantify the degree of cluster synchrony, cluster asymmetry, and switching. We use a variation of the recent dimensionality-reduction technique of Ott and Antonsen [Chaos {\bf 18}, 037113 (2008)] and find an analytic description of the degree of cluster synchrony valid on a globally attracting manifold. Shaped by this manifold, there is an infinite family of steady-state distributions of oscillators, resulting in a high degree of multi-stability in the cluster asymmetry. We also show how through external forcing the degree of asymmetry can be controlled, and suggest that systems displaying cluster synchrony can be used to encode and store data.

preprint2011arXiv

Hierarchical Synchrony of Phase Oscillators in Modular Networks

We study synchronization of sinusoidally coupled phase oscillators on networks with modular structure and a large number of oscillators in each community. Of particular interest is the hierarchy of local and global synchrony, i.e., synchrony within and between communities, respectively. Using the recent ansatz of Ott and Antonsen, we find that the degree of local synchrony can be determined from a set of coupled low-dimensional equations. If the number of communities in the network is large, a low-dimensional description of global synchrony can be also found. Using these results, we study bifurcations between different types of synchrony. We find that, depending on the relative strength of local and global coupling, the transition to synchrony in the network can be mediated by local or global effects.