-
Asymptotic analysis for stationary distributions of scaled reaction networks
Authors:
Linard Hoessly,
Carsten Wiuf,
Panqiu Xia
Abstract:
We study stationary distributions in the context of stochastic reaction networks. In particular, we are interested in complex balanced reaction networks and reduction of such networks by assuming a set of species (called non-interacting species) are degraded fast (and therefore essentially absent in the network), implying some reaction rates are large compared to others. Technically, we assume the…
▽ More
We study stationary distributions in the context of stochastic reaction networks. In particular, we are interested in complex balanced reaction networks and reduction of such networks by assuming a set of species (called non-interacting species) are degraded fast (and therefore essentially absent in the network), implying some reaction rates are large compared to others. Technically, we assume these reaction rates are scaled by a common parameter $N$ and let $N\to\infty$. The limiting stationary distribution as $N\to\infty$ is compared to the stationary distribution of the reduced reaction network obtained by algebraic elimination of the non-interacting species. In general, the limiting stationary distribution might differ from the stationary distribution of the reduced reaction network. We identify various sufficient conditions for when these two distributions are the same, including when the reaction network is detailed balanced and when the set of non-interacting species consists of intermediate species. In the latter case, the limiting stationary distribution essentially retains the form of the complex balanced distribution. This finding is particularly surprising given that the reduced reaction network might be non-weakly reversible and exhibit unconventional kinetics.
△ Less
Submitted 3 February, 2024;
originally announced February 2024.
-
Complex balanced distributions for chemical reaction networks
Authors:
Linard Hoessly,
Carsten Wiuf,
Panqiu Xia
Abstract:
Stationary distributions of continuous time Markov chains (CTMCs) are often a main interest, but hard to find. We consider CTMCs modelling reaction networks, and characterise complex balanced distributions (provided they exist) of reaction networks with arbitrary transition functions through conditions on the cycles of their corresponding digraph. The proof works by constructing a dynamically equi…
▽ More
Stationary distributions of continuous time Markov chains (CTMCs) are often a main interest, but hard to find. We consider CTMCs modelling reaction networks, and characterise complex balanced distributions (provided they exist) of reaction networks with arbitrary transition functions through conditions on the cycles of their corresponding digraph. The proof works by constructing a dynamically equivalent reaction network of disjoint cycles. We further derive a sufficient condition for the existence of a complex balanced distribution, and give precise conditions on when it is necessary. The sufficient condition holds for mass-action kinetics or if the digraph consists of only cyclic connected components. Hence, we fully characterise the existence and form of complex balanced distributions. Moreover, these results can be used to design complex balanced reaction networks with a given stationary distribution.
△ Less
Submitted 10 January, 2023;
originally announced January 2023.
-
On the sum of chemical reactions
Authors:
Linard Hoessly,
Carsten Wiuf,
Panqiu Xia
Abstract:
It is standard in chemistry to represent a sequence of reactions by a single overall reaction, often called a complex reaction in contrast to an elementary reaction. Photosynthesis $6 \text{CO}_2+6 \text{H}_2\text{O} \to \ \text{C}_6\text{H}_{12}\text{O}_6$ $+\ 6 \text{O}_2$ is an example of such complex reaction. We introduce a mathematical operation that corresponds to summing two chemical react…
▽ More
It is standard in chemistry to represent a sequence of reactions by a single overall reaction, often called a complex reaction in contrast to an elementary reaction. Photosynthesis $6 \text{CO}_2+6 \text{H}_2\text{O} \to \ \text{C}_6\text{H}_{12}\text{O}_6$ $+\ 6 \text{O}_2$ is an example of such complex reaction. We introduce a mathematical operation that corresponds to summing two chemical reactions. Specifically, we define an associative and non-communicative operation on the product space $\mathbb{N}_0^n\times \mathbb{N}_0^n$ (representing the reactant and the product of a chemical reaction, respectively). The operation models the overall effect of two reactions happening in succession, one after the other. We study the algebraic properties of the operation and apply the results to stochastic reaction networks, in particular to reachability of states, and to reduction of reaction networks.
△ Less
Submitted 25 March, 2022; v1 submitted 10 May, 2021;
originally announced May 2021.
-
Interplay between COVID-19 vaccines and social measures for ending the SARS-CoV-2 pandemic
Authors:
Glenn Young,
Pengcheng Xiao,
Kenneth Newcomb,
Edwin Michael
Abstract:
The development and authorization of COVID-19 vaccines has provided the clearest path forward to eliminate community spread hence end the ongoing SARS-CoV-2 pandemic. However, the limited pace at which the vaccine can be administered motivates the question, to what extent must we continue to adhere to social intervention measures such as mask wearing and social distancing? To address this question…
▽ More
The development and authorization of COVID-19 vaccines has provided the clearest path forward to eliminate community spread hence end the ongoing SARS-CoV-2 pandemic. However, the limited pace at which the vaccine can be administered motivates the question, to what extent must we continue to adhere to social intervention measures such as mask wearing and social distancing? To address this question, we develop a mathematical model of COVID-19 spread incorporating both vaccine dynamics and socio-epidemiological parameters. We use this model to study two important measures of disease control and eradication, the effective reproductive number $R_t$ and the peak intensive care unit (ICU) caseload, over three key parameters: social measure adherence, vaccination rate, and vaccination coverage. Our results suggest that, due to the slow pace of vaccine administration, social measures must be maintained by a large proportion of the population until a sufficient proportion of the population becomes vaccinated for the pandemic to be eradicated. By contrast, with reduced adherence to social measures, hospital ICU cases will greatly exceed capacity, resulting in increased avoidable loss of life. These findings highlight the complex interplays involved between vaccination and social protective measures, and indicate the practical importance of continuing with extent social measures while vaccines are scaled up to allow the development of the herd immunity needed to end or control SARS-CoV-2 sustainably.
△ Less
Submitted 6 March, 2021;
originally announced March 2021.
-
Proton Conducting Graphene Oxide Coupled Neuron Transistors for Brain-Inspired Cognitive Systems
Authors:
Changjin Wan,
Liqiang Zhu,
Yanghui Liu,
Ping Feng,
Zhaoping Liu,
Hailiang Cao,
Peng Xiao,
Yi Shi,
Qing Wan
Abstract:
Neuron is the most important building block in our brain, and information processing in individual neuron involves the transformation of input synaptic spike trains into an appropriate output spike train. Hardware implementation of neuron by individual ionic/electronic hybrid device is of great significance for enhancing our understanding of the brain and solving sensory processing and complex rec…
▽ More
Neuron is the most important building block in our brain, and information processing in individual neuron involves the transformation of input synaptic spike trains into an appropriate output spike train. Hardware implementation of neuron by individual ionic/electronic hybrid device is of great significance for enhancing our understanding of the brain and solving sensory processing and complex recognition tasks. Here, we provide a proof-of-principle artificial neuron based on a proton conducting graphene oxide (GO) coupled oxide-based electric-double-layer (EDL) transistor with multiple driving inputs and one modulatory input terminal. Paired-pulse facilitation, dendritic integration and orientation tuning were successfully emulated. Additionally, neuronal gain control (arithmetic) in the scheme of rate coding is also experimentally demonstrated. Our results provide a new-concept approach for building brain-inspired cognitive systems.
△ Less
Submitted 20 October, 2015;
originally announced October 2015.