Source author record

Konstantin Mischaikow

Konstantin Mischaikow 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

18works
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

18 published item(s)

preprint2022arXiv

Extending combinatorial regulatory network modeling to include activity control and decay modulation

Understanding how the structure of within-system interactions affects the dynamics of the system is important in many areas of science. We extend a network dynamics modeling platform DSGRN, which combinatorializes both dynamics and parameter space to construct finite but accurate summaries of network dynamics, to new types of interactions. While the standard DSGRN assumes that each network edge controls the rate of abundance of the target node, the new edges may control either activity level or a decay rate of its target. While motivated by processes of post-transcriptional modification and ubiquitination in systems biology, our extension is applicable to the dynamics of any signed directed network.

preprint2022arXiv

Morse Graphs: Topological Tools for Analyzing the Global Dynamics of Robot Controllers

Understanding the global dynamics of a robot controller, such as identifying attractors and their regions of attraction (RoA), is important for safe deployment and synthesizing more effective hybrid controllers. This paper proposes a topological framework to analyze the global dynamics of robot controllers, even data-driven ones, in an effective and explainable way. It builds a combinatorial representation representing the underlying system's state space and non-linear dynamics, which is summarized in a directed acyclic graph, the Morse graph. The approach only probes the dynamics locally by forward propagating short trajectories over a state-space discretization, which needs to be a Lipschitz-continuous function. The framework is evaluated given either numerical or data-driven controllers for classical robotic benchmarks. It is compared against established analytical and recent machine learning alternatives for estimating the RoAs of such controllers. It is shown to outperform them in accuracy and efficiency. It also provides deeper insights as it describes the global dynamics up to the discretization's resolution. This allows to use the Morse graph to identify how to synthesize controllers to form improved hybrid solutions or how to identify the physical limitations of a robotic system.

preprint2021arXiv

Computing linear extensions for polynomial posets subject to algebraic constraints

In this paper we consider the classical problem of computing linear extensions of a given poset which is well known to be a difficult problem. However, in our setting the elements of the poset are multivariate polynomials, and only a small "admissible" subset of these linear extensions, determined implicitly by the evaluation map, are of interest. This seemingly novel problem arises in the study of global dynamics of gene regulatory networks in which case the poset is a Boolean lattice. We provide an algorithm for solving this problem using linear programming for arbitrary partial orders of linear polynomials. This algorithm exploits this additional algebraic structure inherited from the polynomials to efficiently compute the admissible linear extensions. The biologically relevant problem involves multilinear polynomials and we provide a construction for embedding it into an instance of the linear problem.

preprint2020arXiv

Contractibility of a persistence map preimage

This work is motivated by the following question in data-driven study of dynamical systems: given a dynamical system that is observed via time series of persistence diagrams that encode topological features of solutions snapshots, what conclusions can be drawn about solutions of the original dynamical system? In this paper we provide a definition of a persistence diagram for a point in $\mathbb{R}^N$ modeled on piecewise monotone functions. We then provide conditions under which time series of persistence diagrams can be used to guarantee the existence of a fixed point of the flow on $\mathbb{R}^N$ that generates the time series. To obtain this result requires an understanding of the preimage of the persistence map. The main theorem of this paper gives conditions under which these preimages are contractible simplicial complexes.

preprint2019arXiv

Conley index approach to sampled dynamics

The topological method for the reconstruction of dynamics from time series [K. Mischaikow, M. Mrozek, J. Reiss, A. Szymczak. Construction of Symbolic Dynamics from Experimental Time Series, Physical Review Letters, 82 (1999), 1144-1147] is reshaped to improve its range of applicability, particularly in the presence of sparse data and strong expansion. The improvement is based on a multivalued map representation of the data. However, unlike the previous approach, it is not required that the representation has a continuous selector. Instead of a selector, a recently developed new version of Conley index theory for multivalued maps [B. Batko and M. Mrozek. Weak index pairs and the Conley index for discrete multivalued dynamical systems, SIAM J. Applied Dynamical Systems 15 (2016), 1143-1162], [B.Batko. Weak index pairs and the Conley index for discrete multivalued dynamical systems. Part II: properties of the Index, SIAM J. Applied Dynamical Systems 16 (2017), 1587-1617] is used in computations. The existence of a continuous, single-valued generator of the relevant dynamics is guaranteed in the vicinity of the graph of the multivalued map constructed from data. Some numerical examples based on time series derived from the iteration of Hénon type maps are presented.

preprint2019arXiv

Interaction network analysis in shear thickening suspensions

Dense, stabilized, frictional particulate suspensions in a viscous liquid undergo increasingly strong continuous shear thickening (CST) as the solid packing fraction, $ϕ$, increases above a critical volume fraction, and discontinuous shear thickening (DST) is observed for even higher packing fractions. Recent studies have related shear thickening to a transition from mostly lubricated to predominantly frictional contacts with the increase in stress. The rheology and networks of frictional forces from two and three-dimensional simulations of shear-thickening suspensions are studied. These are analyzed using measures of the topology of the network, including tools of persistent homology. We observe that at low stress the frictional interaction networks are predominantly quasi-linear along the compression axis. With an increase in stress, the force networks become more isotropic, forming loops in addition to chain-like structures. The topological measures of Betti numbers and total persistence provide a compact means of describing the mean properties of the frictional force networks and provide a key link between macroscopic rheology and the microscopic interactions. A total persistence measure describing the significance of loops in the force network structure, as a function of stress and packing fraction, shows behavior similar to that of relative viscosity and displays a scaling law near the jamming fraction for both dimensionalities simulated.

preprint2015arXiv

Analysis of Kolmogorov Flow and Rayleigh-Bénard Convection using Persistent Homology

We use persistent homology to build a quantitative understanding of large complex systems that are driven far-from-equilibrium; in particular, we analyze image time series of flow field patterns from numerical simulations of two important problems in fluid dynamics: Kolmogorov flow and Rayleigh-Bénard convection. For each image we compute a persistence diagram to yield a reduced description of the flow field; by applying different metrics to the space of persistence diagrams, we relate characteristic features in persistence diagrams to the geometry of the corresponding flow patterns. We also examine the dynamics of the flow patterns by a second application of persistent homology to the time series of persistence diagrams. We demonstrate that persistent homology provides an effective method both for quotienting out symmetries in families of solutions and for identifying multiscale recurrent dynamics. Our approach is quite general and it is anticipated to be applicable to a broad range of open problems exhibiting complex spatio-temporal behavior.

preprint2015arXiv

Combinatorial Representation of Parameter Space for Switching Systems

We describe the theoretical and computational framework for the Dynamic Signatures for Genetic Regulatory Network (DSGRN) database. The motivation stems from urgent need to understand the global dynamics of biologically relevant signal transduction/gene regulatory networks that have at least 5 to 10 nodes, involve multiple interactions, and decades of parameters. The input to the database computations is a regulatory network, i.e.\ a directed graph with edges indicating up or down regulation, from which a computational model based on switching networks is generated. The phase space dimension equals the number of nodes. The associated parameter space consists of one parameter for each node (a decay rate), and three parameters for each edge (low and high levels of expression, and a threshold at which expression levels change). Since the nonlinearities of switching systems are piece-wise constant, there is a natural decomposition of phase space into cells from which the dynamics can be described combinatorially in terms of a state transition graph. This in turn leads to compact representation of the global dynamics called an annotated Morse graph that identifies recurrent and nonrecurrent. The focus of this paper is on the construction of a natural computable finite decomposition of parameter space into domains where the annotated Morse graph description of dynamics is constant. We use this decomposition to construct an SQL database that can be effectively searched for dynamic signatures such as bistability, stable or unstable oscillations, and stable equilibria. We include two simple 3-node networks to provide small explicit examples of the type information stored in the DSGRN database. To demonstrate the computational capabilities of this system we consider a simple network associated with p53 that involves 5-nodes and a 29-dimensional parameter space.

preprint2015arXiv

Discretization strategies for computing Conley indices and Morse decompositions of flows

Conley indices and Morse decompositions of flows can be found by using algorithms which rigorously analyze discrete dynamical systems. This usually involves integrating a time discretization of the flow using interval arithmetic. We compare the old idea of fixing a time step as a parameters to a time step continuously varying in phase space. We present an example where this second strategy necessarily yields better numerical outputs and prove that our outputs yield a valid Morse decomposition of the given flow.

preprint2015arXiv

Global Dynamics for Steep Sigmoidal Nonlinearities in Two Dimensions

We introduce a novel approach to obtaining mathematically rigorous results on the global dynamics of ordinary differential equations. Motivated by models of regulatory networks, we construct a state transition graph from a piecewise affine ordinary differential equation. We use efficient graph algorithms to compute an associated Morse graph that codifies the recurrent and gradient-like dynamics. We prove that for 2-dimensional systems, the Morse graph defines a Morse decomposition for the dynamics of any smooth differential equation that is sufficiently close to the original piecewise affine ordinary differential equation.

preprint2015arXiv

Topological data analysis of contagion maps for examining spreading processes on networks

Social and biological contagions are influenced by the spatial embeddedness of networks. Historically, many epidemics spread as a wave across part of the Earth's surface; however, in modern contagions long-range edges -- for example, due to airline transportation or communication media -- allow clusters of a contagion to appear in distant locations. Here we study the spread of contagions on networks through a methodology grounded in topological data analysis and nonlinear dimension reduction. We construct "contagion maps" that use multiple contagions on a network to map the nodes as a point cloud. By analyzing the topology, geometry, and dimensionality of manifold structure in such point clouds, we reveal insights to aid in the modeling, forecast, and control of spreading processes. Our approach highlights contagion maps also as a viable tool for inferring low-dimensional structure in networks.

preprint2014arXiv

Evolution of Force Networks in Dense Particulate Media

We introduce novel sets of measures with the goal of describing dynamical properties of force networks in dense particulate systems. The presented approach is based on persistent homology and allows for extracting precise, quantitative measures that describe the evolution of geometric features of the interparticle forces, without necessarily considering the details related to individual contacts between particles. The networks considered emerge from discrete element simulations of two dimensional particulate systems consisting of compressible frictional circular disks. We quantify the evolution of the networks for slowly compressed systems undergoing jamming transition. The main findings include uncovering significant but localized changes of force networks for unjammed systems, global (system-wide) changes as the systems evolve through jamming, to be followed by significantly less dramatic evolution for the jammed states. We consider both connected components, related in loose sense to force chains, and loops, and find that both measures provide a significant insight into the evolution of force networks. In addition to normal, we consider also tangential forces between the particles and find that they evolve in the consistent manner. Consideration of both frictional and frictionless systems leads us to the conclusion that friction plays a significant role in determining the dynamical properties of the considered networks. We find that the proposed approach describes the considered networks in a precise yet tractable manner, allowing to identify novel features which could be difficult or impossible to describe using other approaches.

preprint2014arXiv

Lattice Structures for Attractors II

The algebraic structure of the attractors in a dynamical system determine much of its global dynamics. The collection of all attractors has a natural lattice structure, and this structure can be detected through attracting neighborhoods, which can in principle be computed. Indeed, there has been much recent work on developing and implementing general computational algorithms for global dynamics, which are capable of computing attracting neighborhoods efficiently. Here we address the question of whether all of the algebraic structure of attractors can be captured by these methods.

preprint2013arXiv

Lattice Structures for Attractors I

We describe the basic lattice structures of attractors and repellers in dynamical systems. The structure of distributive lattices allows for an algebraic treatment of gradient-like dynamics in general dynamical systems, both invertible and noninvertible. We separate those properties which rely solely on algebraic structures from those that require some topological arguments, in order to lay a foundation for the development of algorithms to manipulate these structures computationally.

preprint2013arXiv

Reconstructing Functions from Random Samples

From a sufficiently large point sample lying on a compact Riemannian submanifold of Euclidean space, one can construct a simplicial complex which is homotopy-equivalent to that manifold with high confidence. We describe a corresponding result for a Lipschitz-continuous function between two such manifolds. That is, we outline the construction of a simplicial map which recovers the induced maps on homotopy and homology groups with high confidence using only finite sampled data from the domain and range, as well as knowledge of the image of every point sampled from the domain. We provide explicit bounds on the size of the point samples required for such reconstruction in terms of intrinsic properties of the domain, the co-domain and the function. This reconstruction is robust to certain types of bounded sampling and evaluation noise.

preprint2010arXiv

Topology-guided sampling of nonhomogeneous random processes

Topological measurements are increasingly being accepted as an important tool for quantifying complex structures. In many applications, these structures can be expressed as nodal domains of real-valued functions and are obtained only through experimental observation or numerical simulations. In both cases, the data on which the topological measurements are based are derived via some form of finite sampling or discretization. In this paper, we present a probabilistic approach to quantifying the number of components of generalized nodal domains of nonhomogeneous random processes on the real line via finite discretizations, that is, we consider excursion sets of a random process relative to a nonconstant deterministic threshold function. Our results furnish explicit probabilistic a priori bounds for the suitability of certain discretization sizes and also provide information for the choice of location of the sampling points in order to minimize the error probability. We illustrate our results for a variety of random processes, demonstrate how they can be used to sample the classical nodal domains of deterministic functions perturbed by additive noise and discuss their relation to the density of zeros.

preprint2007arXiv

Probabilistic validation of homology computations for nodal domains

Homology has long been accepted as an important computable tool for quantifying complex structures. In many applications, these structures arise as nodal domains of real-valued functions and are therefore amenable only to a numerical study based on suitable discretizations. Such an approach immediately raises the question of how accurate the resulting homology computations are. In this paper, we present a probabilistic approach to quantifying the validity of homology computations for nodal domains of random fields in one and two space dimensions, which furnishes explicit probabilistic a priori bounds for the suitability of certain discretization sizes. We illustrate our results for the special cases of random periodic fields and random trigonometric polynomials.