4,927 research outputs found
Language-based Abstractions for Dynamical Systems
Ordinary differential equations (ODEs) are the primary means to modelling
dynamical systems in many natural and engineering sciences. The number of
equations required to describe a system with high heterogeneity limits our
capability of effectively performing analyses. This has motivated a large body
of research, across many disciplines, into abstraction techniques that provide
smaller ODE systems while preserving the original dynamics in some appropriate
sense. In this paper we give an overview of a recently proposed
computer-science perspective to this problem, where ODE reduction is recast to
finding an appropriate equivalence relation over ODE variables, akin to
classical models of computation based on labelled transition systems.Comment: In Proceedings QAPL 2017, arXiv:1707.0366
Challenges in Quantitative Abstractions for Collective Adaptive Systems
Like with most large-scale systems, the evaluation of quantitative properties
of collective adaptive systems is an important issue that crosscuts all its
development stages, from design (in the case of engineered systems) to runtime
monitoring and control. Unfortunately it is a difficult problem to tackle in
general, due to the typically high computational cost involved in the analysis.
This calls for the development of appropriate quantitative abstraction
techniques that preserve most of the system's dynamical behaviour using a more
compact representation. This paper focuses on models based on ordinary
differential equations and reviews recent results where abstraction is achieved
by aggregation of variables, reflecting on the shortcomings in the state of the
art and setting out challenges for future research.Comment: In Proceedings FORECAST 2016, arXiv:1607.0200
Reduction of dynamical biochemical reaction networks in computational biology
Biochemical networks are used in computational biology, to model the static
and dynamical details of systems involved in cell signaling, metabolism, and
regulation of gene expression. Parametric and structural uncertainty, as well
as combinatorial explosion are strong obstacles against analyzing the dynamics
of large models of this type. Multi-scaleness is another property of these
networks, that can be used to get past some of these obstacles. Networks with
many well separated time scales, can be reduced to simpler networks, in a way
that depends only on the orders of magnitude and not on the exact values of the
kinetic parameters. The main idea used for such robust simplifications of
networks is the concept of dominance among model elements, allowing
hierarchical organization of these elements according to their effects on the
network dynamics. This concept finds a natural formulation in tropical
geometry. We revisit, in the light of these new ideas, the main approaches to
model reduction of reaction networks, such as quasi-steady state and
quasi-equilibrium approximations, and provide practical recipes for model
reduction of linear and nonlinear networks. We also discuss the application of
model reduction to backward pruning machine learning techniques
Model reduction for stochastic CaMKII reaction kinetics in synapses by graph-constrained correlation dynamics.
A stochastic reaction network model of Ca(2+) dynamics in synapses (Pepke et al PLoS Comput. Biol. 6 e1000675) is expressed and simulated using rule-based reaction modeling notation in dynamical grammars and in MCell. The model tracks the response of calmodulin and CaMKII to calcium influx in synapses. Data from numerically intensive simulations is used to train a reduced model that, out of sample, correctly predicts the evolution of interaction parameters characterizing the instantaneous probability distribution over molecular states in the much larger fine-scale models. The novel model reduction method, 'graph-constrained correlation dynamics', requires a graph of plausible state variables and interactions as input. It parametrically optimizes a set of constant coefficients appearing in differential equations governing the time-varying interaction parameters that determine all correlations between variables in the reduced model at any time slice
Syntactic Markovian Bisimulation for Chemical Reaction Networks
In chemical reaction networks (CRNs) with stochastic semantics based on
continuous-time Markov chains (CTMCs), the typically large populations of
species cause combinatorially large state spaces. This makes the analysis very
difficult in practice and represents the major bottleneck for the applicability
of minimization techniques based, for instance, on lumpability. In this paper
we present syntactic Markovian bisimulation (SMB), a notion of bisimulation
developed in the Larsen-Skou style of probabilistic bisimulation, defined over
the structure of a CRN rather than over its underlying CTMC. SMB identifies a
lumpable partition of the CTMC state space a priori, in the sense that it is an
equivalence relation over species implying that two CTMC states are lumpable
when they are invariant with respect to the total population of species within
the same equivalence class. We develop an efficient partition-refinement
algorithm which computes the largest SMB of a CRN in polynomial time in the
number of species and reactions. We also provide an algorithm for obtaining a
quotient network from an SMB that induces the lumped CTMC directly, thus
avoiding the generation of the state space of the original CRN altogether. In
practice, we show that SMB allows significant reductions in a number of models
from the literature. Finally, we study SMB with respect to the deterministic
semantics of CRNs based on ordinary differential equations (ODEs), where each
equation gives the time-course evolution of the concentration of a species. SMB
implies forward CRN bisimulation, a recently developed behavioral notion of
equivalence for the ODE semantics, in an analogous sense: it yields a smaller
ODE system that keeps track of the sums of the solutions for equivalent
species.Comment: Extended version (with proofs), of the corresponding paper published
at KimFest 2017 (http://kimfest.cs.aau.dk/
Computationally-efficient stochastic cluster dynamics method for modeling damage accumulation in irradiated materials
An improved version of a recently developed stochastic cluster dynamics (SCD)
method {[}Marian, J. and Bulatov, V. V., {\it J. Nucl. Mater.} \textbf{415}
(2014) 84-95{]} is introduced as an alternative to rate theory (RT) methods for
solving coupled ordinary differential equation (ODE) systems for irradiation
damage simulations. SCD circumvents by design the curse of dimensionality of
the variable space that renders traditional ODE-based RT approaches inefficient
when handling complex defect population comprised of multiple (more than two)
defect species. Several improvements introduced here enable efficient and
accurate simulations of irradiated materials up to realistic (high) damage
doses characteristic of next-generation nuclear systems. The first improvement
is a procedure for efficiently updating the defect reaction-network and event
selection in the context of a dynamically expanding reaction-network. Next is a
novel implementation of the -leaping method that speeds up SCD
simulations by advancing the state of the reaction network in large time
increments when appropriate. Lastly, a volume rescaling procedure is introduced
to control the computational complexity of the expanding reaction-network
through occasional reductions of the defect population while maintaining
accurate statistics. The enhanced SCD method is then applied to model defect
cluster accumulation in iron thin films subjected to triple ion-beam
(, and \text{H\ensuremath{{}^{+}}})
irradiations, for which standard RT or spatially-resolved kinetic Monte Carlo
simulations are prohibitively expensive
Forward and Backward Bisimulations for Chemical Reaction Networks
We present two quantitative behavioral equivalences over species of a
chemical reaction network (CRN) with semantics based on ordinary differential
equations. Forward CRN bisimulation identifies a partition where each
equivalence class represents the exact sum of the concentrations of the species
belonging to that class. Backward CRN bisimulation relates species that have
the identical solutions at all time points when starting from the same initial
conditions. Both notions can be checked using only CRN syntactical information,
i.e., by inspection of the set of reactions. We provide a unified algorithm
that computes the coarsest refinement up to our bisimulations in polynomial
time. Further, we give algorithms to compute quotient CRNs induced by a
bisimulation. As an application, we find significant reductions in a number of
models of biological processes from the literature. In two cases we allow the
analysis of benchmark models which would be otherwise intractable due to their
memory requirements.Comment: Extended version of the CONCUR 2015 pape
- …