-
A survey on the Riemann-Lebesgue integrability in non-additive setting
Authors:
Anca Croitoru,
Alina Gavrilut,
Alina Iosif,
Anna Rita Sambucini
Abstract:
We present in this survey some results regarding Riemann_Lebesgue integrability with respect to arbitrary non-additive set functions.
We present in this survey some results regarding Riemann_Lebesgue integrability with respect to arbitrary non-additive set functions.
△ Less
Submitted 30 September, 2023;
originally announced October 2023.
-
A Computational Approach to Polynomial Conservation Laws
Authors:
Aurélien Desoeuvres,
Alexandru Iosif,
Christoph Lüders,
Ovidiu Radulescu,
Hamid Rahkooy,
Matthias Seiß,
Thomas Sturm
Abstract:
For polynomial ODE models, we introduce and discuss the concepts of exact and approximate conservation laws, which are the first integrals of the full and truncated sets of ODEs. For fast-slow systems, truncated ODEs describe the fast dynamics. We define compatibility classes as subsets of the state space, obtained by equating the conservation laws to constants. A set of conservation laws is compl…
▽ More
For polynomial ODE models, we introduce and discuss the concepts of exact and approximate conservation laws, which are the first integrals of the full and truncated sets of ODEs. For fast-slow systems, truncated ODEs describe the fast dynamics. We define compatibility classes as subsets of the state space, obtained by equating the conservation laws to constants. A set of conservation laws is complete when the corresponding compatibility classes contain a finite number of steady states. Complete sets of conservation laws can be used for model order reduction and for studying the multistationarity of the model. We provide algorithmic methods for computing linear, monomial, and polynomial conservation laws of polynomial ODE models and for testing their completeness. The resulting conservation laws and their completeness are either independent or dependent on the parameters. In the latter case, we provide parametric case distinctions. In particular, we propose a new method to compute polynomial conservation laws by comprehensive Gröbner systems and syzygies.
Keywords: First integrals, chemical reaction networks, polynomial conservation laws, syzygies, comprehensive Gröbner systems.
△ Less
Submitted 13 November, 2023; v1 submitted 30 December, 2022;
originally announced December 2022.
-
Reduction of Chemical Reaction Networks with Approximate Conservation Laws
Authors:
Aurélien Desoeuvres,
Alexandru Iosif,
Christoph Lüders,
Ovidiu Radulescu,
Hamid Rahkooy,
Matthias Seiß,
Thomas Sturm
Abstract:
Model reduction of fast-slow chemical reaction networks based on the quasi-steady state approximation fails when the fast subsystem has first integrals. We call these first integrals approximate conservation laws. In order to define fast subsystems and identify approximate conservation laws, we use ideas from tropical geometry. We prove that any approximate conservation law evolves slower than all…
▽ More
Model reduction of fast-slow chemical reaction networks based on the quasi-steady state approximation fails when the fast subsystem has first integrals. We call these first integrals approximate conservation laws. In order to define fast subsystems and identify approximate conservation laws, we use ideas from tropical geometry. We prove that any approximate conservation law evolves slower than all the species involved in it and therefore represents a supplementary slow variable in an extended system. By elimination of some variables of the extended system, we obtain networks without approximate conservation laws, which can be reduced by standard singular perturbation methods. The field of applications of approximate conservation laws covers the quasi-equilibrium approximation, well known in biochemistry. We discuss reductions of slow-fast as well as multiple timescale systems. Networks with multiple timescales have hierarchical relaxation. At a given timescale, our multiple timescale reduction method defines three subsystems composed of (i) slaved fast variables satisfying algebraic equations, (ii) slow driving variables satisfying reduced ordinary differential equations, and (iii) quenched much slower variables that are constant. The algebraic equations satisfied by fast variables define chains of nested normally hyberbolic invariant manifolds. In such chains, faster manifolds are of higher dimension and contain the slower manifolds. Our reduction methods are introduced algorithmically for networks with monomial reaction rates and linear, monomial or polynomial approximate conservation laws. We propose symbolic algorithms to reshape and rescale the networks such that geometric singular perturbation theory can be applied to them, test the applicability of the theory, and finally reduce the networks. As a proof of concept, we apply this method to a model of the TGF-b signaling pathway.
△ Less
Submitted 30 July, 2023; v1 submitted 27 December, 2022;
originally announced December 2022.
-
BioDynaMo: a general platform for scalable agent-based simulation
Authors:
Lukas Breitwieser,
Ahmad Hesam,
Jean de Montigny,
Vasileios Vavourakis,
Alexandros Iosif,
Jack Jennings,
Marcus Kaiser,
Marco Manca,
Alberto Di Meglio,
Zaid Al-Ars,
Fons Rademakers,
Onur Mutlu,
Roman Bauer
Abstract:
Motivation: Agent-based modeling is an indispensable tool for studying complex biological systems. However, existing simulators do not always take full advantage of modern hardware and often have a field-specific software design.
Results: We present a novel simulation platform called BioDynaMo that alleviates both of these problems. BioDynaMo features a general-purpose and high-performance simul…
▽ More
Motivation: Agent-based modeling is an indispensable tool for studying complex biological systems. However, existing simulators do not always take full advantage of modern hardware and often have a field-specific software design.
Results: We present a novel simulation platform called BioDynaMo that alleviates both of these problems. BioDynaMo features a general-purpose and high-performance simulation engine. We demonstrate that BioDynaMo can be used to simulate use cases in: neuroscience, oncology, and epidemiology. For each use case we validate our findings with experimental data or an analytical solution. Our performance results show that BioDynaMo performs up to three orders of magnitude faster than the state-of-the-art baseline. This improvement makes it feasible to simulate each use case with one billion agents on a single server, showcasing the potential BioDynaMo has for computational biology research.
Availability: BioDynaMo is an open-source project under the Apache 2.0 license and is available at www.biodynamo.org. Instructions to reproduce the results are available in supplementary information.
Contact: lukas.breitwieser@inf.ethz.ch, a.s.hesam@tudelft.nl, omutlu@ethz.ch, r.bauer@surrey.ac.uk
Supplementary information: Available at https://doi.org/10.5281/zenodo.4501515
△ Less
Submitted 5 February, 2021; v1 submitted 11 June, 2020;
originally announced June 2020.
-
Analysis of the Conradi-Kahle Algorithm for Detecting Binomiality on Biological Models
Authors:
Alexandru Iosif,
Hamid Rahkooy
Abstract:
We analyze the Conradi-Kahle Algorithm for detecting binomiality. We present experiments using two implementations of the algorithm in Macaulay2 and Maple on biological models and assess the performance of the algorithm on these models. We compare the two implementations with each other and with Gröbner bases computations up to their performance on these biological models.
We analyze the Conradi-Kahle Algorithm for detecting binomiality. We present experiments using two implementations of the algorithm in Macaulay2 and Maple on biological models and assess the performance of the algorithm on these models. We compare the two implementations with each other and with Gröbner bases computations up to their performance on these biological models.
△ Less
Submitted 14 December, 2019;
originally announced December 2019.
-
Efficiently and Effectively Recognizing Toricity of Steady State Varieties
Authors:
Dima Grigoriev,
Alexandru Iosif,
Hamid Rahkooy,
Thomas Sturm,
Andreas Weber
Abstract:
We consider the problem of testing whether the points in a complex or real variety with non-zero coordinates form a multiplicative group or, more generally, a coset of a multiplicative group. For the coset case, we study the notion of shifted toric varieties which generalizes the notion of toric varieties. This requires a geometric view on the varieties rather than an algebraic view on the ideals.…
▽ More
We consider the problem of testing whether the points in a complex or real variety with non-zero coordinates form a multiplicative group or, more generally, a coset of a multiplicative group. For the coset case, we study the notion of shifted toric varieties which generalizes the notion of toric varieties. This requires a geometric view on the varieties rather than an algebraic view on the ideals. We present algorithms and computations on 129 models from the BioModels repository testing for group and coset structures over both the complex numbers and the real numbers. Our methods over the complex numbers are based on Gröbner basis techniques and binomiality tests. Over the real numbers we use first-order characterizations and employ real quantifier elimination. In combination with suitable prime decompositions and restrictions to subspaces it turns out that almost all models show coset structure. Beyond our practical computations, we give upper bounds on the asymptotic worst-case complexity of the corresponding problems by proposing single exponential algorithms that test complex or real varieties for toricity or shifted toricity. In the positive case, these algorithms produce generating binomials. In addition, we propose an asymptotically fast algorithm for testing membership in a binomial variety over the algebraic closure of the rational numbers.
△ Less
Submitted 15 April, 2020; v1 submitted 9 October, 2019;
originally announced October 2019.
-
Properties of the Riemann-Lebesgue integrability in the non-additive case
Authors:
Domenico Candeloro,
Anca Croitoru,
Alina Gavrilut,
Alina Iosif,
Anna Rita Sambucini
Abstract:
We study Riemann-Lebesgue integrability of a vector function relative to an arbitrary non-negative set function. We obtain some classical integral properties. Results regarding the continuity properties of the integral and relationships among Riemann-Lebesgue, Birkhoff simple and Gould integrabilities are also established.
We study Riemann-Lebesgue integrability of a vector function relative to an arbitrary non-negative set function. We obtain some classical integral properties. Results regarding the continuity properties of the integral and relationships among Riemann-Lebesgue, Birkhoff simple and Gould integrabilities are also established.
△ Less
Submitted 10 May, 2019;
originally announced May 2019.
-
Multistationarity in the space of total concentrations for systems that admit a monomial parametrization
Authors:
Carsten Conradi,
Alexandru Iosif,
Thomas Kahle
Abstract:
We apply tools from real algebraic geometry to the problem of multistationarity of chemical reaction networks. A particular focus is on the case of reaction networks whose steady states admit a monomial parametrization. For such systems we show that in the space of total concentrations multistationarity is scale invariant: if there is multistationarity for some value of the total concentrations, t…
▽ More
We apply tools from real algebraic geometry to the problem of multistationarity of chemical reaction networks. A particular focus is on the case of reaction networks whose steady states admit a monomial parametrization. For such systems we show that in the space of total concentrations multistationarity is scale invariant: if there is multistationarity for some value of the total concentrations, then there is multistationarity on the entire ray containing this value (possibly for different rate constants) -- and vice versa. Moreover, for these networks it is possible to decide about multistationarity independent of the rate constants by formulating semi-algebraic conditions that involve only concentration variables. These conditions can easily be extended to include total concentrations. Hence quantifier elimination may give new insights into multistationarity regions in the space of total concentrations. To demonstrate this, we show that for the distributive phosphorylation of a protein at two binding sites multistationarity is only possible if the total concentration of the substrate is larger than either the total concentration of the kinase or the total concentration of the phosphatase. This result is enabled by the chamber decomposition of the space of total concentrations from polyhedral geometry. Together with the corresponding sufficiency result of Bihan et al. this yields a characterization of multistationarity up to lower dimensional regions.
△ Less
Submitted 13 August, 2019; v1 submitted 18 October, 2018;
originally announced October 2018.
-
A Geometric View of the Sieve of Eratosthenes
Authors:
Alexandru Iosif
Abstract:
We study the geometry of the Sieve of Eratosthenes. We introduce some concepts as Focals and Extremes. We find a symmetry in the distribution of the Focals (all the information about the primes is contained into a small set of numbers). We find that there is a geometric order in the Sieve and we give a formula for the greatest remainder that returns the same quotient.
We study the geometry of the Sieve of Eratosthenes. We introduce some concepts as Focals and Extremes. We find a symmetry in the distribution of the Focals (all the information about the primes is contained into a small set of numbers). We find that there is a geometric order in the Sieve and we give a formula for the greatest remainder that returns the same quotient.
△ Less
Submitted 30 December, 2022; v1 submitted 25 December, 2011;
originally announced December 2011.