Researcher profile

Per Sebastian Skardal

Per Sebastian Skardal contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 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.

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

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

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.