BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Michael Margaliot

Michael Margaliot contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.DSmath.OCmath.CASystems and ControlMolecular NetworksQuantitative MethodsComputational ComplexityData Structures and Algorithmsmath.COSubcellular Processes

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
10topics
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

19 published item(s)

preprint2025arXiv

Instability of equilibrium and convergence to periodic orbits in strongly 2-cooperative systems

We consider time-invariant nonlinear $n$-dimensional strongly $2$-cooperative systems, that is, systems that map the set of vectors with up to weak sign variation to its interior. Strongly $2$-cooperative systems enjoy a strong Poincare-Bendixson property: bounded solutions that maintain a positive distance from the set of equilibria converge to a periodic solution. For strongly $2$-cooperative systems whose trajectories evolve in a bounded and invariant set that contains a single unstable equilibrium, we provide a simple criterion for the existence of periodic trajectories. Moreover, we explicitly characterize a positive-measure set of initial conditions which yield solutions that asymptotically converge to a periodic trajectory. We demonstrate our theoretical results using two models from systems biology, the $n$-dimensional Goodwin oscillator and a $4$-dimensional biomolecular oscillator with RNA-mediated regulation, and provide numerical simulations that verify the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Negative feedback and oscillations in a model for mRNA translation

The ribosome flow model (RFM) is a phenomenological model for the unidirectional flow of particles along a 1D chain of $n$ sites. The RFM has been extensively used to study the dynamics of ribosome flow along a single mRNA molecule during translation. In this case, the particles model ribosomes and each site corresponds to a consecutive group of codons. Networks of interconnected RFMs have been used to model and analyze large-scale translation in the cell and, in particular, the effects of competition for shared resources. Here, we analyze the RFM with a negative feedback connection from the protein production rate to the initiation rate. This models, for example, the production of proteins that inhibit the translation of their own mRNA. Using tools from the theory of 2-cooperative dynamical systems, we provide a simple condition guaranteeing that the closed-loop system admits at least one non-trivial periodic solution. When this condition holds, we also explicitly characterize a large set of initial conditions such that any solution emanating from this set converges to a non-trivial periodic solution. Such a solution corresponds to a periodic pattern of ribosome densities along the mRNA, and to a periodic pattern of protein production.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Multiplicative and additive compounds via Kronecker products and Kronecker sums

Compound matrices play an important role in many fields of mathematics and have recently found new applications in systems and control theory. However, the explicit formulas for these compounds are non-trivial and not always easy to use. Here, we derive new formulas for the multiplicative and additive compounds of a matrix using Kronecker products and sums. This provides a new approach to matrix compounds based on the well-known and powerful theory of Kronecker products and sums. We demonstrate several applications of these new formulas, including deriving a new expression for the additive compound of the product of two matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Compound matrices in systems and control theory: a tutorial

The multiplicative and additive compounds of a matrix play an important role in several fields of mathematics including geometry, multi-linear algebra, combinatorics, and the analysis of nonlinear time-varying dynamical systems. There is a growing interest in applications of these compounds, and their generalizations, in systems and control theory. The goal of this tutorial paper is to provide a gentle and self-contained introduction to these topics with an emphasis on the geometric interpretation of the compounds, and to describe some of their recent applications including several non-trivial generalizations of positive systems, cooperative systems, contracting systems, and more.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Minimum effort decentralized control design for contracting network systems

We consider the problem of making a networked system contracting by designing minimal effort local controllers. Our method combines a hierarchical contraction characterization and a matrix-balancing approach to stabilizing a Metzler matrix via minimal diagonal perturbations. We demonstrate our approach by designing local controllers that render contractive a network of FitzHugh-Nagumo neurons with a general topology of interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On matrices whose exponential is a P-matrix

A matrix is called a P-matrix if all its principal minors are positive. P-matrices have found important applications in functional analysis, mathematical programming, and dynamical systems theory. We introduce a new class of real matrices denoted~$\EP$. A matrix is in~$\EP$ if and only if its matrix exponential is a P-matrix for all positive times. In other words, $A\in \EP$ if and only if the transition matrix of the linear system~$\dot x=Ax$ is a P-matrix for any positive time~$t$. We analyze the properties of this new class of matrices and describe an application of our theoretical results to opinion dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Verifying $k$-Contraction without Computing $k$-Compounds

Compound matrices have found applications in many fields of science including systems and control theory. In particular, a sufficient condition for $k$-contraction is that a logarithmic norm (also called matrix measure) of the $k$-additive compound of the Jacobian is uniformly negative. However, this may be difficult to check in practice because the $k$-additive compound of an $n\times n$ matrix has dimensions $\binom{n}{k}\times \binom{n}{k}$. For an $n\times n$ matrix $A$, we prove a duality relation between the $k$ and $(n-k)$ compounds of $A$. We use this duality relation to derive a sufficient condition for $k$-contraction that does not require the computation of any $k$-compounds. We demonstrate our results by deriving a sufficient condition for $k$-contraction of an $n$-dimensional Hopfield network that does not require to compute any compounds. In particular, for $k=2$ this sufficient condition implies that the network is $2$-contracting and this implies a strong asymptotic property: every bounded solution of the network converges to an equilibrium point, that may not be unique. This is relevant, for example, when using the Hopfield network as an associative memory that stores patterns as equilibrium points of the dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Behavior of Totally Positive Differential Systems Near a Periodic Solution

A time-varying nonlinear dynamical system is called a totally positive differential system (TPDS) if its Jacobian admits a special sign pattern: it is tri-diagonal with positive entries on the super- and sub-diagonals. If the vector field of a TPDS is T-periodic then every bounded trajectory converges to a T-periodic solution. In particular, when the vector field is time-invariant every bounded trajectory of a TPDS converges to an equlbrium. Here, we use the spectral theory of oscillatory matrices to analyze the behavior near a periodic solution of a TPDS. This yields information on the perturbation directions that lead to the fastest and slowest convergence to or divergence from the periodic solution. We demonstrate the theoretical results using a model from systems biology called the ribosome flow model.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Diagonal Stability of Discrete-time $k$-Positive linear Systems with Applications to Nonlinear Systems

A linear dynamical system is called $k$-positive if its dynamics maps the set of vectors with up to $k-1$ sign variations to itself. For $k=1$, this reduces to the important class of positive linear systems. Since stable positive linear time-invariant (LTI) systems always admit a diagonal quadratic Lyapunov function, i.e. they are diagonally stable, we may expect that this holds also for stable $k$-positive systems. We show that, in general, this is not the case both in the continuous-time (CT) and discrete-time (DT) case. We then focus on DT $k$-positive linear systems and introduce the new notion of DT $k$-diagonal stability. It is shown that this is a necessary condition for standard DT diagonal stability. We demonstrate an application of this new notion to the analysis of a class of DT nonlinear systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the Diagonal Stability of $k$-Positive Linear Systems

We consider $k$-positive linear systems, that is, systems that map the set of vectors with up to $k-1$ sign variations to itself. For $k=1$, this reduces to positive linear systems. It is well-known that stable positive linear time invariant (LTI) systems admit a diagonal Lyapunov function. This property has many important implications. A natural question is whether stable $k$-positive systemsalso admit a diagonal Lyapunov function. This paper shows that, in general, the answer is no. However, for both continuous-time and discrete-time $n$-dimensional systems that are $(n-1)$-positive we provide a sufficient condition for diagonal stability.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Variability in mRNA Translation: A Random Matrix Theory Approach

The rate of mRNA translation depends on the initiation, elongation, and termination rates of ribosomes along the mRNA. These rates depend on many "local" factors like the abundance of free ribosomes and tRNA molecules in the vicinity of the mRNA molecule. All these factors are stochastic and their experimental measurements are also noisy. An important question is how protein production in the cell is affected by this considerable variability. We develop a new theoretical framework for addressing this question by modeling the rates as identically and independently distributed random variables and using tools from random matrix theory to analyze the steady-state production rate. The analysis reveals a principle of universality: the average protein production rate depends only on the of the set of possible values that the random variable may attain. This explains how total protein production can be stabilized despite the overwhelming stochasticticity underlying cellular processes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Is my system of ODEs $k$-cooperative?

A linear dynamical system is called positive if its flow maps the non-negative orthant to itself. More precisely, it maps the set of vectors with zero sign variations to itself. A linear dynamical system is called $k$-positive if its flow maps the set of vectors with up to $k-1$ sign variations to itself. A nonlinear dynamical system is called $k$-cooperative if its variational system, which is a time-varying linear dynamical system, is $k$-positive. These systems have special asymptotic properties. For example, it was recently shown that strongly $2$-cooperative systems satisfy a strong Poincaré-Bendixson property. Positivity and~$k$-positivity are easy to verify in terms of the sign-pattern of the matrix in the dynamics. However, these sign conditions are not invariant under a coordinate transformation. A natural question is to determine if a given~$n$-dimensional system is $k$-positive up to a coordinate transformation. We study this problem for two special kinds of transformations: permutations and scaling by a signature matrix. For any $n\geq 4$ and~$k\in\{2,\dots, n-2\}$, we provide a graph-theoretical necessary and sufficient condition for $k$-positivity up to such coordinate transformations. We describe an application of our results to a specific class of Lotka-Volterra systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximizing average throughput in oscillatory biological synthesis systems: an optimal control approach

A dynamical system entrains to a periodic input if its state converges globally to an attractor with the same period. In particular, for a constant input the state converges to a unique equilibrium point for any initial condition. We consider the problem of maximizing a weighted average of the system's output along the periodic attractor. The gain of entrainment is the benefit achieved by using a non-constant periodic input relative to a constant input with the same time average. Such a problem amounts to optimal allocation of resources in a periodic manner. We formulate this problem as a periodic optimal control problem which can be analyzed by means of the Pontryagin maximum principle or solved numerically via powerful software packages. We then apply our framework to a class of occupancy models that appear frequently in biological synthesis systems and other applications. We show that, perhaps surprisingly, constant inputs are optimal for various architectures. This suggests that the presence of non-constant periodic signals, which frequently appear in biological occupancy systems, is a signature of an underlying time-varying objective functional being optimized.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Exponent of Several Classes of Oscillatory Matrices

Oscillatory matrices were introduced in the seminal work of Gantmacher and Krein. An $n\times n$ matrix $A$ is called oscillatory if all its minors are nonnegative and there exists a positive integer $k$ such that all minors of $A^k$ are positive. The smallest $k$ for which this holds is called the exponent of the oscillatory matrix $A$. Gantmacher and Krein showed that the exponent is always smaller than or equal to $n-1$. An important and nontrivial problem is to determine the exact value of the exponent. Here we use the successive elementary bidiagonal factorization of oscillatory matrices, and its graph-theoretic representation, to derive an explicit expression for the exponent of several classes of oscillatory matrices, and a nontrivial upper-bound on the exponent for several other classes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Random attraction in the TASEP model

The totally asymmetric simple exclusion process (TASEP) is a basic model of statistical mechanics that has found numerous applications. We consider the case of TASEP with a finite chain where particles may enter from the left and leave to the right at prescribed rates. This model can be formulated as a Markov process with a finite number of states. Due to the irreducibility of the process it is well-known that the probability distribution on the states is globally attracted to a unique equilibrium distribution. We extend this result to the more detailed level of individual trajectories. To do so we formulate TASEP as a random dynamical system. Our main result is that the trajectories from all possible initial conditions contract to each other yielding the existence of a random attractor that consists of a single trajectory almost surely. This implies that in the long run TASEP ``filters out'' any perturbation that changes the state of the particles along the chain. In order to prove our main result we first establish that any random dynamical system on a finite state space possesses both a global random pullback attractor and a global random forward attractor. This observation appears to be missing in the literature. We then provide sufficient and necessary conditions for these attractors to be singletons. Finally, we show that TASEP satisfies one of these conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A Generalization of Smillie's Theorem on Strongly Cooperative Tridiagonal Systems

Smillie (1984) proved an interesting result on the stability of nonlinear, time-invariant, strongly cooperative, and tridiagonal dynamical systems. This result has found many applications in models from various fields including biology, ecology, and chemistry. Smith (1991) has extended Smillie's result and proved entrainment in the case where the vector field is time-varying and periodic. We use the theory of linear totally nonnegative differential systems developed by Schwarz (1970) to give a generalization of these two results. This is based on weakening the requirement for strong cooperativity to cooperativity, and adding an additional observability-type condition.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Output Selection and Observer Design for Boolean Control Networks: A Sub-Optimal Polynomial-Complexity Algorithm

Using a graph-theoretic approach, we derive a new sufficient condition for observability of a Boolean control network (BCN). Based on this condition, we describe two algorithms: the first selects a set of nodes so that observing this set makes the BCN observable. The second algorithm builds an observer for the observable BCN. Both algorithms are sub-optimal, as they are based on a sufficient but not necessary condition for observability. Yet their time-complexity is linear in the length of the description of the BCN, rendering them feasible for large-scale networks. We discuss how these results can be used to provide a sub-optimal yet polynomial-complexity algorithm for the minimal observability problem in BCNs. Some of the theoretical results are demonstrated using a BCN model of the core network regulating the mammalian cell cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

A Polynomial-Time Algorithm for Solving the Minimal Observability Problem in Conjunctive Boolean Networks

Many complex systems in biology, physics, and engineering include a large number of state-variables, and measuring the full state of the system is often impossible. Typically, a set of sensors is used to measure part of the state-variables. A system is called observable if these measurements allow to reconstruct the entire state of the system. When the system is not observable, an important and practical problem is how to add a \emph{minimal} number of sensors so that the system becomes observable. This minimal observability problem is practically useful and theoretically interesting, as it pinpoints the most informative nodes in the system. We consider the minimal observability problem for an important special class of Boolean networks, called conjunctive Boolean networks (CBNs). Using a graph-theoretic approach, we provide a necessary and sufficient condition for observability of a CBN with $n$ state-variables, and an efficient~$O(n^2)$-time algorithm for solving the minimal observability problem. We demonstrate the usefulness of these results by studying the properties of a class of random CBNs.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Minimal Controllability of Conjunctive Boolean Networks is NP-Complete

Given a conjunctive Boolean network (CBN) with $n$ state-variables, we consider the problem of finding a minimal set of state-variables to directly affect with an input so that the resulting conjunctive Boolean control network (CBCN) is controllable. We give a necessary and sufficient condition for controllability of a CBCN; an $O(n^2)$-time algorithm for testing controllability; and prove that nonetheless the minimal controllability problem for CBNs is NP-hard.

0 citations0 reviewsNo code linkedOpen access