Source author record

Chaojie Yuan

Chaojie Yuan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.PR math.NA Molecular Networks Numerical Analysis

Catalog footprint

What is connected

2works

4topics

3close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Constrained Langevin approximation for the Togashi-Kaneko model of autocatalytic reactions

The Togashi Kaneko model (TK model), introduced by Togashi and Kaneko in 2001, is a simple stochastic reaction network that displays discreteness-induced transitions between meta-stable patterns. Here we study a constrained Langevin approximation (CLA) of this model. The CLA, obtained by Anderson et al. in 2019, is an obliquely reflected diffusion process on the positive orthant and hence it respects the constrain that chemical concentrations are never negative. We show that the CLA is a Feller process, is positive Harris recurrent, and converges exponentially fast to the unique stationary distribution. We also characterize the stationary distribution and show that it has finite moments. In addition, we simulate both the TK model and its CLA in various dimensions. For example, we describe how the TK model switches between meta-stable patterns in dimension 6. Our simulations suggest that, under the classical scaling, the CLA is a good approximation to the TK model in terms of both the stationary distribution and the transition times between patterns.

preprint2020arXiv

Variance of finite difference methods for reaction networks with non-Lipschitz rate functions

Parametric sensitivity analysis is a critical component in the study of mathematical models of physical systems. Due to its simplicity, finite difference methods are used extensively for this analysis in the study of stochastically modeled reaction networks. Different coupling methods have been proposed to build finite difference estimators, with the "split coupling," also termed the "stacked coupling," yielding the lowest variance in the vast majority of cases. Analytical results related to this coupling are sparse, and include an analysis of the variance of the coupled processes under the assumption of globally Lipschitz intensity functions [Anderson, SIAM Numerical Analysis, Vol. 50, 2012]. Because of the global Lipschitz assumption utilized in [Anderson, SIAM Numerical Analysis, Vol. 50, 2012], the main result there is only applicable to a small percentage of the models found in the literature, and it was conjectured that similar results should hold for a much wider class of models. In this paper we demonstrate this conjecture to be true by proving the variance of the coupled processes scales in the desired manner for a large class of non-Lipschitz models. We further extend the analysis to allow for time dependence in the parameters. In particular, binary systems with or without time-dependent rate parameters, a class of models that accounts for the vast majority of systems considered in the literature, satisfy the assumptions of our theory.