77,016 research outputs found
A finite state projection algorithm for the stationary solution of the chemical master equation
The chemical master equation (CME) is frequently used in systems biology to
quantify the effects of stochastic fluctuations that arise due to biomolecular
species with low copy numbers. The CME is a system of ordinary differential
equations that describes the evolution of probability density for each
population vector in the state-space of the stochastic reaction dynamics. For
many examples of interest, this state-space is infinite, making it difficult to
obtain exact solutions of the CME. To deal with this problem, the Finite State
Projection (FSP) algorithm was developed by Munsky and Khammash (Jour. Chem.
Phys. 2006), to provide approximate solutions to the CME by truncating the
state-space. The FSP works well for finite time-periods but it cannot be used
for estimating the stationary solutions of CMEs, which are often of interest in
systems biology. The aim of this paper is to develop a version of FSP which we
refer to as the stationary FSP (sFSP) that allows one to obtain accurate
approximations of the stationary solutions of a CME by solving a finite
linear-algebraic system that yields the stationary distribution of a
continuous-time Markov chain over the truncated state-space. We derive bounds
for the approximation error incurred by sFSP and we establish that under
certain stability conditions, these errors can be made arbitrarily small by
appropriately expanding the truncated state-space. We provide several examples
to illustrate our sFSP method and demonstrate its efficiency in estimating the
stationary distributions. In particular, we show that using a quantised tensor
train (QTT) implementation of our sFSP method, problems admitting more than 100
million states can be efficiently solved.Comment: 8 figure
A finite state projection algorithm for the stationary solution of the chemical master equation
The chemical master equation (CME) is frequently used in systems biology to
quantify the effects of stochastic fluctuations that arise due to biomolecular
species with low copy numbers. The CME is a system of ordinary differential
equations that describes the evolution of probability density for each
population vector in the state-space of the stochastic reaction dynamics. For
many examples of interest, this state-space is infinite, making it difficult to
obtain exact solutions of the CME. To deal with this problem, the Finite State
Projection (FSP) algorithm was developed by Munsky and Khammash (Jour. Chem.
Phys. 2006), to provide approximate solutions to the CME by truncating the
state-space. The FSP works well for finite time-periods but it cannot be used
for estimating the stationary solutions of CMEs, which are often of interest in
systems biology. The aim of this paper is to develop a version of FSP which we
refer to as the stationary FSP (sFSP) that allows one to obtain accurate
approximations of the stationary solutions of a CME by solving a finite
linear-algebraic system that yields the stationary distribution of a
continuous-time Markov chain over the truncated state-space. We derive bounds
for the approximation error incurred by sFSP and we establish that under
certain stability conditions, these errors can be made arbitrarily small by
appropriately expanding the truncated state-space. We provide several examples
to illustrate our sFSP method and demonstrate its efficiency in estimating the
stationary distributions. In particular, we show that using a quantised tensor
train (QTT) implementation of our sFSP method, problems admitting more than 100
million states can be efficiently solved.Comment: 8 figure
Massive Fields of Arbitrary Integer Spin in Symmetrical Einstein Space
We study the propagation of gauge fields with arbitrary integer spins in the
symmetrical Einstein space of any dimensionality. We reduce the problem of
obtaining a gauge-invariant Lagrangian of integer spin fields in such
background to an purely algebraic problem of finding a set of operators with
certain features using the representation of high-spin fields in the form of
some vectors of pseudo-Hilbert space. We consider such construction in the
linear order in the Riemann tensor and scalar curvature and also present an
explicit form of interaction Lagrangians and gauge transformations for massive
particles with spins 1 and 2 in terms of symmetrical tensor fields.Comment: 15 pages, latex, no figures,minor change
Algebra, coalgebra, and minimization in polynomial differential equations
We consider reasoning and minimization in systems of polynomial ordinary
differential equations (ode's). The ring of multivariate polynomials is
employed as a syntax for denoting system behaviours. We endow this set with a
transition system structure based on the concept of Lie-derivative, thus
inducing a notion of L-bisimulation. We prove that two states (variables) are
L-bisimilar if and only if they correspond to the same solution in the ode's
system. We then characterize L-bisimilarity algebraically, in terms of certain
ideals in the polynomial ring that are invariant under Lie-derivation. This
characterization allows us to develop a complete algorithm, based on building
an ascending chain of ideals, for computing the largest L-bisimulation
containing all valid identities that are instances of a user-specified
template. A specific largest L-bisimulation can be used to build a reduced
system of ode's, equivalent to the original one, but minimal among all those
obtainable by linear aggregation of the original equations. A computationally
less demanding approximate reduction and linearization technique is also
proposed.Comment: 27 pages, extended and revised version of FOSSACS 2017 pape
Skyrmions and the Nuclear Force
The derivation of the nucleon-nucleon force from the Skyrme model is
reexamined. Starting from previous results for the potential energy of
quasistatic solutions, we show that a calculation using the Born-Oppenheimer
approximation properly taking into account the mixing of nucleon resonances,
leads to substantial central attraction. We obtain a potential that is in
qualitative agreement with phenomenological potentials. We also study the
non-adiabatic corrections, such as the velocity dependent transition
potentials, and discuss their importance.Comment: 24 pages, UPR-0124M
- …