BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Thomas G. Kurtz

Thomas G. Kurtz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.PRmath.DSmath.NAPopulations and Evolution

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2023arXiv

Looking forwards and backwards: dynamics and genealogies of locally regulated populations

We introduce a broad class of spatial models to describe how spatially heterogeneous populations live, die, and reproduce. Individuals are represented by points of a point measure, whose birth and death rates can depend both on spatial position and local population density, defined via the convolution of the point measure with a nonnegative kernel. We pass to three different scaling limits: an interacting superprocess, a nonlocal partial differential equation (PDE), and a classical PDE. The classical PDE is obtained both by first scaling time and population size to pass to the nonlocal PDE, and then scaling the kernel that determines local population density; and also (when the limit is a reaction-diffusion equation) by simultaneously scaling the kernel width, timescale and population size in our individual based model. A novelty of our model is that we explicitly model a juvenile phase: offspring are thrown off in a Gaussian distribution around the location of the parent, and reach (instant) maturity with a probability that can depend on the population density at the location at which they land. Although we only record mature individuals, a trace of this two-step description remains in our population models, resulting in novel limits governed by a nonlinear diffusion. Using a lookdown representation, we retain information about genealogies and, in the case of deterministic limiting models, use this to deduce the backwards in time motion of the ancestral lineage of a sampled individual. We observe that knowing the history of the population density is not enough to determine the motion of ancestral lineages in our model. We also investigate the behaviour of lineages for three different deterministic models of a population expanding its range as a travelling wave: the Fisher-KPP equation, the Allen-Cahn equation, and a porous medium equation with logistic growth.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Finite Markov chains coupled to general Markov processes and an application to metastability I

We consider a diffusion given by a small noise perturbation of a dynamical system driven by a potential function with a finite number of local minima. The classical results of Freidlin and Wentzell show that the time this diffusion spends in the domain of attraction of one of these local minima is approximately exponentially distributed and hence the diffusion should behave approximately like a Markov chain on the local minima. By the work of Bovier and collaborators, the local minima can be associated with the small eigenvalues of the diffusion generator. Applying a Markov mapping theorem, we use the eigenfunctions of the generator to couple this diffusion to a Markov chain whose generator has eigenvalues equal to the eigenvalues of the diffusion generator that are associated with the local minima and establish explicit formulas for conditional probabilities associated with this coupling. The fundamental question then becomes to relate the coupled Markov chain to the approximate Markov chain suggested by the results of Freidlin and Wentzel. In Part II of this work, we provide a complete analysis of this relationship in the special case of a double-well potential in one dimension. More generally, the coupling can be constructed for a general class of Markov processes and any finite set of eigenvalues of the generator.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Finite Markov chains coupled to general Markov processes and an application to metastability II

We consider a diffusion given by a small noise perturbation of a dynamical system driven by a potential function with a finite number of local minima. The classical results of Freidlin and Wentzell show that the time this diffusion spends in the domain of attraction of one of these local minima is approximately exponentially distributed and hence the diffusion should behave approximately like a Markov chain on the local minima. By the work of Bovier and collaborators, the local minima can be associated with the small eigenvalues of the diffusion generator. In Part I of this work, by applying a Markov mapping theorem, we used the eigenfunctions of the generator to couple this diffusion to a Markov chain whose generator has eigenvalues equal to the eigenvalues of the diffusion generator that are associated with the local minima and established explicit formulas for conditional probabilities associated with this coupling. The fundamental question now becomes to relate the coupled Markov chain to the approximate Markov chain suggested by the results of Freidlin and Wentzel. In this paper, we take up this question and provide a complete analysis of this relationship in the special case of a double-well potential in one dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Central limit theorems and diffusion approximations for multiscale Markov chain models

Ordinary differential equations obtained as limits of Markov processes appear in many settings. They may arise by scaling large systems, or by averaging rapidly fluctuating systems, or in systems involving multiple time-scales, by a combination of the two. Motivated by models with multiple time-scales arising in systems biology, we present a general approach to proving a central limit theorem capturing the fluctuations of the original model around the deterministic limit. The central limit theorem provides a method for deriving an appropriate diffusion (Langevin) approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Weak and strong solutions of general stochastic models

Typically, a stochastic model relates stochastic "inputs" and, perhaps, controls to stochastic "outputs". A general version of the Yamada-Watanabe and Engelbert theorems relating existence and uniqueness of weak and strong solutions of stochastic equations is given in this context. A notion of {\em compatibility} between inputs and outputs is critical in relating the general result to its classical forebears. The relationship between the compatibility condition and the usual formulation of stochastic differential equations driven by semimartingales is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Error analysis of tau-leap simulation methods

We perform an error analysis for numerical approximation methods of continuous time Markov chain models commonly found in the chemistry and biochemistry literature. The motivation for the analysis is to be able to compare the accuracy of different approximation methods and, specifically, Euler tau-leaping and midpoint tau-leaping. We perform our analysis under a scaling in which the size of the time discretization is inversely proportional to some (bounded) power of the norm of the state of the system. We argue that this is a more appropriate scaling than that found in previous error analyses in which the size of the time discretization goes to zero independent of the rest of the model. Under the present scaling, we show that midpoint tau-leaping achieves a higher order of accuracy, in both a weak and a strong sense, than Euler tau-leaping; a result that is in contrast to previous analyses. We present examples that demonstrate our findings.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Poisson representations of branching Markov and measure-valued branching processes

Representations of branching Markov processes and their measure-valued limits in terms of countable systems of particles are constructed for models with spatially varying birth and death rates. Each particle has a location and a "level," but unlike earlier constructions, the levels change with time. In fact, death of a particle occurs only when the level of the particle crosses a specified level $r$, or for the limiting models, hits infinity. For branching Markov processes, at each time $t$, conditioned on the state of the process, the levels are independent and uniformly distributed on $[0,r]$. For the limiting measure-valued process, at each time $t$, the joint distribution of locations and levels is conditionally Poisson distributed with mean measure $K(t)\times\varLambda$, where $\varLambda$ denotes Lebesgue measure, and $K$ is the desired measure-valued process. The representation simplifies or gives alternative proofs for a variety of calculations and results including conditioning on extinction or nonextinction, Harris's convergence theorem for supercritical branching processes, and diffusion approximations for processes in random environments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Product-form stationary distributions for deficiency zero chemical reaction networks

We consider stochastically modeled chemical reaction systems with mass-action kinetics and prove that a product-form stationary distribution exists for each closed, irreducible subset of the state space if an analogous deterministically modeled system with mass-action kinetics admits a complex balanced equilibrium. Feinberg's deficiency zero theorem then implies that such a distribution exists so long as the corresponding chemical network is weakly reversible and has a deficiency of zero. The main parameter of the stationary distribution for the stochastically modeled system is a complex balanced equilibrium value for the corresponding deterministically modeled system. We also generalize our main result to some non-mass-action kinetics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Separation of time-scales and model reduction for stochastic reaction networks

A stochastic model for a chemical reaction network is embedded in a one-parameter family of models with species numbers and rate constants scaled by powers of the parameter. A systematic approach is developed for determining appropriate choices of the exponents that can be applied to large complex networks. When the scaling implies subnetworks have different time-scales, the subnetworks can be approximated separately providing insight into the behavior of the full network through the analysis of these lower dimensional approximations.

0 citations0 reviewsNo code linkedOpen access