Distribution-based bisimulation for labelled Markov processes

In this paper we propose a (sub)distribution-based bisimulation for labelled Markov processes and compare it with earlier definitions of state and event bisimulation, which both only compare states. In contrast to those state-based bisimulations, our distribution bisimulation is weaker, but corresponds more closely to linear properties. We construct a logic and a metric to describe our distribution bisimulation and discuss linearity, continuity and compositional properties.Comment: Accepted by FORMATS 201

Space-Time Evolution of the Oscillator, Rapidly moving in a random media

We study the quantum-mechanical evolution of the nonrelativistic oscillator, rapidly moving in the media with the random vector fields. We calculate the evolution of the level probability distribution as a function of time, and obtain rapid level diffusion over the energy levels. Our results imply a new mechanism of charmonium dissociation in QCD media.Comment: 32 pages, 13 figure

Process algebra modelling styles for biomolecular processes

We investigate how biomolecular processes are modelled in process algebras, focussing on chemical reactions. We consider various modelling styles and how design decisions made in the definition of the process algebra have an impact on how a modelling style can be applied. Our goal is to highlight the often implicit choices that modellers make in choosing a formalism, and illustrate, through the use of examples, how this can affect expressability as well as the type and complexity of the analysis that can be performed

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/

Annotations for Rule-Based Models

The chapter reviews the syntax to store machine-readable annotations and describes the mapping between rule-based modelling entities (e.g., agents and rules) and these annotations. In particular, we review an annotation framework and the associated guidelines for annotating rule-based models of molecular interactions, encoded in the commonly used Kappa and BioNetGen languages, and present prototypes that can be used to extract and query the annotations. An ontology is used to annotate models and facilitate their description

Biochemical Reaction Rules with Constraints

International audienceWe propose React(C), an expressive programming language for stochastic modeling and simulation in systems biology, that is based on biochemical reactions with constraints. We prove that React(C) can express the stochastic pi-calculus, in contrast to previous rule-based programming languages, and further illustrate the high expressiveness of React(C). We present a stochastic simulator for React(C) independently of the choice of the constraint language C. Our simulator must decide for a given reaction rule whether it can be applied to the current biochemical solution. We show that this decision problem is NP-complete for arbitrary constraint systems C, and that it can be solved in polynomial time for rules of bounded arity. In practice, we propose to solve this problem by constraint programming

GaSb quantum rings in GaAs/AlxGa1−xAs quantum wells

We report the results of continuous and time-resolved photoluminescence measurements on type-II GaSb quantum rings embedded within GaAs/AlxGa1-xAs quantum wells. A range of samples were grown with different well widths, compensation-doping concentrations within the wells, and number of quantum-ring layers. We find that each of these variants have no discernible effect on the radiative recombination, except for the very narrowest (5 nm) quantum well. In contrast, singleparticle numerical simulations of the sample predict changes in photoluminescence energy of up to 200 meV. This remarkable difference is explained by the strong Coulomb binding of electrons to rings that are multiply charged with holes. The resilience of the emission to compensation doping indicates that multiple hole occupancy of the quantum rings is required for efficient carrier recombination, regardless of whether these holes come from doping or excitation

Semiclassical theory of surface plasmons in spheroidal clusters

A microscopic theory of linear response based on the Vlasov equation is extended to systems having spheroidal equilibrium shape. The solution of the linearized Vlasov equation, which gives a semiclassical version of the random phase approximation, is studied for electrons moving in a deformed equilibrium mean field. The deformed field has been approximated by a cavity of spheroidal shape, both prolate and oblate. Contrary to spherical systems, there is now a coupling among excitations of different multipolarity induced by the interaction among constituents. Explicit calculations are performed for the dipole response of deformed clusters of different size. In all cases studied here the photoabsorption strength for prolate clusters always displays a typical double-peaked structure. For oblate clusters we find that the high--frequency component of the plasmon doublet can get fragmented in the medium size region ($N \sim 250$). This fragmentation is related to the presence of two kinds of three-dimensional electron orbits in oblate cavities. The possible scaling of our semiclassical equations with the valence electron number and density is investigated.Comment: 23 pages, 8 figures, revised version, includes discussion of scalin

The Attributed Pi Calculus with Priorities

International audienceWe present the attributed $\pi$-calculus for modeling concurrent systems with interaction constraints depending on the values of attributes of processes. The $\pi$-calculus serves as a constraint language underlying the $\pi$-calculus. Interaction constraints subsume priorities, by which to express global aspects of populations. We present a nondeterministic and a stochastic semantics for the attributed $\pi$-calculus. We show how to encode the $\pi$-calculus with priorities and polyadic synchronization $\pi$@ and thus dynamic compartments, as well as the stochastic $\pi$-calculus with concurrent objects spico. We illustrate the usefulness of the attributed $\pi$-calculus for modeling biological systems at two particular examples: Euglena’s spatial movement in phototaxis, and cooperative protein binding in gene regulation of bacteriophage lambda. Furthermore, population-based model is supported beside individual-based modeling. A stochastic simulation algorithm for the attributed $\pi$-calculus is derived from its stochastic semantics. We have implemented a simulator and present experimental results, that confirm the practical relevance of our approach

A High-Level Language for Rule-Based Modelling

Rule-based languages such as Kappa excel in their support for handling the combinatorial complexities prevalent in many biological systems, including signalling pathways. But Kappa provides little structure for organising rules, and large models can therefore be hard to read and maintain. This paper introduces a high-level, modular extension of Kappa called LBS-κ. We demonstrate the constructs of the language through examples and three case studies: a chemotaxis switch ring, a MAPK cascade, and an insulin signalling pathway. We then provide a formal definition of LBS-κ through an abstract syntax and a translation to plain Kappa. The translation is implemented in a compiler tool which is available as a web application. We finally demonstrate how to increase the expressivity of LBS-κ through embedded scripts in a general-purpose programming language, a technique which we view as generally applicable to other domain specific languages