In this paper we consider stochastic processes taking values in a set of continuum graphs we call graphemes, defined as equivalence classes of sequences of vertices labelled by N embedded in an uncountable Polish space (with the cardinality of the continuum), together with an N x N connection matrix with entries 0 or 1 specifying the absence or presence of edges between pairs of vertices. In particular, we construct a Markov process on a Polish state space G of graphemes suitable to describe the time-space path of countable graphs. The class of dynamics we propose arises by specifying simple rules for the evolution of finite graphs and passing to the limit of infinite graphs. The evolution of graphemes is characterised by well-posed martingale problems, and leads to strong Markov processes with the Feller property.
We consider a system of interacting diffusions labeled by a geographic space that is given by the hierarchical group Ω N \Omega _N of order N ∈ N N\in \mathbb {N} . Individuals live in colonies and are subject to resampling and migration as long as they are active. Each colony has a seed-bank into which individuals can retreat to become dormant, suspending their resampling and migration until they become active again. The migration kernel has a hierarchical structure: individuals hop between colonies at a rate that depends on the hierarchical distance between the colonies. The seed-bank has a layered structure: when individuals become dormant they acquire a colour that determines the rate at which they become active again. The latter allows us to model seed-banks whose wake-up times have a fat tail. We analyse a system of coupled stochastic differential equations that describes the population in the large-colony-size limit. For fixed N ∈ N N\in \mathbb {N} , the system exhibits a dichotomy between coexistence (= locally multi-type equilibria) and clustering (= locally mono-type equilibria). To identify the clustering regime, i.e., the range of parameters controlling the migration and the seed-bank for which clustering prevails, we apply a necessary and sufficient criterion derived in (Electron. J. Probab. 27 (2022), Paper No. 18, 88, DOI 10.1214/21-ejp728. MR4377129) that is valid for any geographic space given by a countable Abelian group endowed with the discrete topology. The diffusion function controlling the resampling does not play a role in this criterion. We carry out a multi-scale renormalisation analysis in the hierarchical mean-field limit N → ∞ N\to \infty . In particular, we show that block averages on hierarchical space-time scale k ∈ N k \in \mathbb {N} perform a diffusion with a renormalised diffusion function that depends on k k . In the clustering regime, after an appropriate scaling with k k , this diffusion function converges to the Fisher-Wright diffusion function as k → ∞ k\to \infty , irrespective of the diffusion function controlling the resampling. Thus, the system exhibits full universality on large space-time scales in terms of the scaling limit. For several subclasses of parameters we identify the speed at which the scaled renormalised diffusion function converges to the Fisher-Wright diffusion as k → ∞ k\to \infty . This speed in turn determines the speed at which mono-type clusters grow in space and time. We show that the seed-bank reduces the speed compared to the model without seed-bank. This reduction is the result of a delicate interplay between migration and seed-bank, which in the limit as N → ∞ N\to \infty can be worked out explicitly.
For two Polish state spaces $E_X$ and $E_Y$, and an operator $G_X$, we obtain existence and uniqueness of a $G_X$-martingale problem provided there is a dual process $Y$ on $E_Y$ solving a $G_Y$-martingale problem. Duality here means the existence of a rich function $H$ and transition kernels $(\mu_t)_{t\geq 0}$ on $E_X$ such that $\mathbf E_y[H(x,Y_t)] = \int \mu_t(x,dx') H(x',y)$ for all $(x,y) \in E_X \times E_Y$ and $t\geq 0$. While duality is well-known to imply uniqueness of the $G_X$-martingale problem, we give here a set of conditions under which duality also implies existence without using approximation techniques. As examples, we treat branching and resampling models, such as Feller's branching diffusion and the Fleming-Viot superprocess.
In Athreya et al. (2021), models from population genetics were used to define stochastic dynamics in the space of graphons arising as continuum limits of dense graphs. In the present paper we exhibit an example of a simple neutral population genetics model for which this dynamics is a Markovian diffusion that can be characterized as the solution of a martingale problem. In particular, we consider a Markov chain in the space of finite graphs that resembles a Moran model with resampling and mutation. We encode the finite graphs as graphemes, which can be represented as a triple consisting of a vertex set (or more generally, a topological space), an adjacency matrix, and a sampling (Borel) measure. We equip the space of graphons with convergence of sample subgraph densities and show that the grapheme-valued Markov chain converges to a grapheme-valued diffusion as the number of vertices goes to infinity. We show that the grapheme-valued diffusion has a stationary distribution that is linked to the Griffiths-Engen-McCloskey (GEM) distribution. In a companion paper (Greven et al. 2023), we build up a general theory for obtaining grapheme-valued diffusions via genealogies of models in population genetics.
In the article Quenched large deviation principle for words in a letter sequence, Probab. Theory Relat. Fields 148, no. 3/4 (2010), 403-456 we derived a quenched large deviation principle for the empirical process of words obtained by cutting an i.i.d. sequence of letters according to an independent renewal process. We derived a representation of the associated rate function for stationary word processes in terms of certain specific relative entropies. Our proof of this representation is correct when the mean word length is finite, but is flawed when the mean word length is infinite. In this paper we fix the flaw in the proof. Along the way we derive new representations of the rate function that are interesting in their own right. A key ingredient in the proof is the observation that if the rate function in the annealed large deviation principle is finite at a stationary word process, then the letters in the tail of the long words in this process are typical.
We consider a system of interacting Fisher-Wright diffusions with seed-bank. Individuals carry type one of two types, live in colonies, and are subject to resampling and migration as long as they are active. Each colony has a structured seed-bank into which individuals can retreat to become dormant, suspending their resampling and migration until they become active again. As geographic space labelling the colonies we consider a countable Abelian group endowed with the discrete topology. In earlier work we showed that the system has a one-parameter family of equilibria controlled by the relative density of the two types. Moreover, these equilibria exhibit a dichotomy of coexistence (= locally multi-type equilibrium) versus clustering (= locally mono-type equilibrium). We identified the parameter regimes for which these two phases occur, and found that these regimes are different when the mean wake-up time of a dormant individual is finite or infinite. The goal of the present paper is to establish the finite-systems scheme, i.e., identify how a finite truncation of the system (both in the geographic space and in the seed-bank) behaves as both the time and the truncation level tend to infinity, properly tuned together. If the wake-up time has finite mean, then there is a single universality class for the scaling limit. On the other hand, if the wake-up time has infinite mean, then there are two universality classes depending on how fast the truncation level of the seed-bank grows compared to the truncation level of the geographic space.
When we cut an i.i.d. sequence of letters into words according to an independent renewal process, we obtain an i.i.d. sequence of words. In the annealed large deviation principle (LDP) for the empirical process of words, the rate function is the specific relative entropy of the observed law of words w.r.t. the reference law of words. In the present paper we consider the quenched LDP, i.e., we condition on a typical letter sequence. We focus on the case where the renewal process has an algebraic tail. The rate function turns out to be a sum of two terms, one being the annealed rate function, the other being proportional to the specific relative entropy of the observed law of letters w.r.t. the reference law of letters, with the former being obtained by concatenating the words and randomising the location of the origin. The proportionality constant equals the tail exponent of the renewal process. Earlier work by Birkner considered the case where the renewal process has an exponential tail, in which case the rate function turns out to be the first term on the set where the second term vanishes and to be infinite elsewhere. In a companion paper the annealed and the quenched LDP are applied to the collision local time of transient random walks, and the existence of an intermediate phase for a class of interacting stochastic systems is established.
We consider a system of interacting Fisher-Wright diffusions with seed-bank. Individuals live in colonies and are subject to resampling and migration as long as they are active. Each colony has a structured seed-bank into which individuals can retreat to become dormant, suspending their resampling and migration until they become active again. As geographic space labelling the colonies we consider a countable Abelian group $\mathbb{G}$ endowed with the discrete topology. The key example of interest is the Euclidean lattice $\mathbb{G}=\mathbb{Z}^d$. Our goal is to classify the long-time behaviour of the system in terms of the underlying model parameters. In particular, we want to understand in what way the seed-bank enhances genetic diversity. We introduce three models of increasing generality, namely, individuals become dormant: (1) in the seed-bank of their colony; (2) in the seed-bank of their colony while adopting a random colour that determines their wake-up time; (3) in the seed-bank of a random colony while adopting a random colour. The extension in (2) allows us to model wake-up times with fat tails while preserving the Markov property of the evolution. For each of the three models we show that the system converges to a unique equilibrium depending on a single density parameter that is determined by the initial state, and exhibits a dichotomy of coexistence (= locally multi-type equilibrium) versus clustering (= locally mono-type equilibrium) depending on the parameters controlling the migration and the seed-bank. The dichotomy between clustering and coexistence in model 1 is determined by migration only. In models (2) and (3), when the wake-up time has infinite mean, the dichotomy is determined by both the exchange with the seed-bank and migration. It turns out that the seed-bank affects the long-time behaviour both quantitatively and qualitatively.
In neutrally evolving populations subject to resampling and migration on geographic spaces, the longtime behaviour exhibits a dichotomy between clustering of types versus coexistence of types (that is between convergence to monotype equilibria versus multitype equilibria). In the classical setting, which of the two scenarios prevails depends on whether the difference between the positions of two lineages in the dual is a recurrent random walk or a transient random walk. In the present contribution, for two classes of models we present results on how this dichotomy is affected when the resampling has either high volatility or low volatility: (I) The spatial Cannings model with block resampling, where a positive fraction of the individuals in the next generation may inherit the type of a single individual in the previous generation. (II) The spatial Fleming–Viot model with seed-bank, where individuals may become dormant for a while, suspending resampling and migration, until they become active again. We will see that the dichotomy is shifted towards more clustering in class (I) and towards more coexistence in class (II). In particular, in class (II) we will see that for critically recurrent random walks (which are typical for two-dimensional geographic spaces) an infinite seed-bank turns clustering into coexistence. We will also explore the impact of allowing resampling that is controlled by a random environment. Along the way we discuss robustness, universality and critical dimension.
The paper has four goals. First, we want to generalize the classical concept of the branching property so that it becomes applicable for historical and genealogical processes (using the coding of genealogies by ($V$-marked) ultrametric measure spaces leading to state spaces $\mathbb{U}$ resp. $\mathbb{U}^V$). The processes are defined by well-posed martingale problems. In particular we want to complement the corresponding concept of infinite divisibility developed in \cite{infdiv} for this context. Second one of the two main points, we want to find a corresponding characterization of the generators of branching processes more precisely their martingale problems which is both easy to apply and general enough to cover a wide range of state spaces. As a third goal we want to obtain the branching property of the $\mathbb{U}$-valued Feller diffusion respectively $\mathbb{U}^V$-valued super random walk and the historical process on countable geographic spaces the latter as two examples of a whole zoo of spatial processes we could treat. The fourth goal is to show the robustness of the method and to get the branching property for genealogies marked with ancestral path, giving the line of descent moving through the ancestors and space, leading to path-marked ultra-metric measure spaces. This processes are constructed here giving our second major result. The starting point for all four points is the Feller diffusion model, the final goal the (historical) super random walk model.
We survey results on the description of stochastically evolving genealogies of populations and marked genealogies of multitype populations or spatial populations via tree-valued Markov processes on (marked) ultrametric measure spaces. In particular we explain the choice of state spaces and their topologies, describe the dynamics of genealogical Fleming-Viot and branching models by well-posed martingale problems, and formulate the typical results on the longtime behavior. Furthermore we discuss the basic techniques of proofs and sketch as two key tools of analysis the different forms of duality and the Girsanov transformation.
We consider the evolution of the genealogy of the population currently alive in a Feller branching diffusion model. In contrast to the approach via labeled trees in the continuum random tree world, the genealogies are modeled as equivalence classes of ultrametric measure spaces, the elements of the space $\mathbb{U}$. This space is Polish and has a rich semigroup structure for the genealogy. We focus on the evolution of the genealogy in time and the large time asymptotics conditioned both on survival up to present time and on survival forever. We prove existence, uniqueness and Feller property of solutions of the martingale problem for this genealogy valued, i.e., $\mathbb{U}$-valued Feller diffusion. We give the precise relation to the time-inhomogeneous $\mathbb{U}_1$-valued Fleming-Viot process. The uniqueness is shown via Feynman-Kac duality with the distance matrix augmented Kingman coalescent. Using a semigroup operation on $\mathbb{U}$, called concatenation, together with the branching property we obtain a L{\'e}vy-Khintchine formula for $\mathbb{U}$-valued Feller diffusion and we determine explicitly the L{\'e}vy measure on $\mathbb{U}\setminus\{0\}$. From this we obtain for $h>0$ the decomposition into depth-$h$ subfamilies, a representation of the process as concatenation of a Cox point process of genealogies of single ancestor subfamilies. Furthermore, we will identify the $\mathbb{U}$-valued process conditioned to survive until a finite time $T$. We study long time asymptotics, such as generalized quasi-equilibrium and Kolmogorov-Yaglom limit law on the level of ultrametric measure spaces. We also obtain various representations of the long time limits.
We consider the evolution of the genealogy of the population currently alive in a Feller branching diffusion model. In contrast to the approach via labeled trees in the continuum random tree world, the genealogies are modeled as equivalence classes of ultrametric measure spaces, the elements of the space $\mathbb{U}$. This space is Polish and has a rich semigroup structure for the genealogy. We focus on the evolution of the genealogy in time and the large time asymptotics conditioned both on survival up to present time and on survival forever. We prove existence, uniqueness and Feller property of solutions of the martingale problem for this genealogy valued, i.e., $\mathbb{U}$-valued Feller diffusion. We give the precise relation to the time-inhomogeneous $\mathbb{U}_1$-valued Fleming-Viot process. The uniqueness is shown via Feynman-Kac duality with the distance matrix augmented Kingman coalescent. Using a semigroup operation on $\mathbb{U}$, called concatenation, together with the branching property we obtain a L\'evy-Khintchine formula for $\mathbb{U}$-valued Feller diffusion and we determine explicitly the L\'evy measure on $\mathbb{U}\setminus\{0\}$. From this we obtain for $h>0$ the decomposition into depth-$h$ subfamilies, a representation of the process as concatenation of a Cox point process of genealogies of single ancestor subfamilies. Furthermore, we will identify the $\mathbb{U}$-valued process conditioned to survive until a finite time $T$. We study long time asymptotics, such as generalized quasi-equilibrium and Kolmogorov-Yaglom limit law on the level of ultrametric measure spaces. We also obtain various representations of the long time limits.
The paper reviews the results obtained for spatial population models and the evolution of the genealogies of these populations during the last decade by the author and his coworkers. The focus is on their large scale behaviour and on the analysis of universality classes of large scale behaviour via the methods of the hierarchical mean-field limit or via the spatial continuum limit and from another angel, the finite system scheme . We use genealogical information to analyze the type and location structure and vice versa. To apply this approach and to explain effects in biological situations we extend the classical model classes in new directions. Namely we look as population models here at: Fleming-Viot genealogies in continuum geographic space, Cannings models with block resampling (reducing diversity), Fisher-Wright diffusion with coloured seedbanks (enhancing diversity) and evolving genealogies of Fleming-Viot models with selection and rare mutation.
The goal of this work is to decompose random populations with a genealogy in subfamilies of a given degree of kinship and to obtain a notion of infinitely divisible genealogies. We model the genealogical structure of a population by (equivalence classes of) ultrametric measure spaces (um-spaces) as elements of the Polish space U which we recall. In order to then analyze the family structure in this coding we introduce an algebraic structure on um-spaces (a consistent collection of semigroups). This allows us to obtain a path of collections of subfamilies of fixed kinship h (described as ultrametric measure spaces), for every depth h as a measurable functional of the genealogy. Random elements in the semigroup are studied, in particular infinitely divisible random variables. Here we define infinite divisibility of random genealogies as the property that the h-tops can be represented as concatenation of independent identically distributed h-forests for every h and obtain a Levy-Khintchine representation of this object and a corresponding representation via a concatenation of points of a Poisson point process of h-forests. Finally the case of discrete and marked um-spaces is treated allowing to apply the results to both the individual based and most important spatial populations. The results have various applications. In particular the case of the genealogical (U-valued) Feller diffusion and genealogical (U V -valued) super random walk is treated based on the present work in [DG18b] and [GRG]. In the part II of this paper we go in a different direction and refine the study in the case of continuum branching populations, give a refined analysis of the Laplace functional and give a representation in terms of a Cox process on h-trees, rather than forests.
This paper studies the spatial coalescent on ^2. In our setting, the partition elements are located at the sites of ^2 and undergo local delayed coalescence and migration. That is, pairs of partition elements located at the same site coalesce into one partition element after exponential waiting times. In addition, the partition elements perform independent random walks. The system starts in either locally finite configurations or in configurations containing countably many partition elements per site. These two situations are relevant if the coalescent is used to study the scaling limits for genealogies in Moran models respectively interacting Fisher-Wright diffusions (or Fleming-Viot processes), which is the key application of the present work. Our goal is to determine the longtime behavior with an initial population of countably many individuals per site restricted to a box [-t^α/2, t^α/2]^2 ∩^2 and observed at time t^β with 1 ≥β≥α≥ 0. We study both asymptotics, as t→∞, for a fixed value of α as the parameter β∈[α,1] varies, and for a fixed β, as the parameter α∈ [0,β] varies. This exhibits the genealogical structure of the mono-type clusters arising in 2-dimensional Moran and Fisher-Wright systems. (... for more see the actual preprint)
In an earlier paper, we introduced and studied a system of hierarchically interacting measure-valued random processes which describes a large population of individuals carrying types and living in colonies labelled by the hierarchical group of order $N$. The individuals are subject to migration, resampling on all hierarchical scales simultaneously. Upon resampling, a random positive fraction of the population in a block of colonies inherits the type of a random single individual in that block, which is why we refer to our system as the hierarchical Cannings process. In the present paper, we study a version of the hierarchical Cannings process in random environment, namely, the resampling measures controlling the change of type of individuals in different blocks are chosen randomly with a given mean and are kept fixed in time (= the quenched setting). We give a necessary and sufficient condition under which a multi-type equilibrium is approached (= coexistence) as opposed to a mono-type equilibrium (= clustering). Moreover, in the hierarchical mean-field limit $N \to \infty$, with the help of a renormalization analysis, we obtain a full picture of the space-time scaling behaviour of block averages on all hierarchical scales simultaneously. We show that the $k$-block averages are distributed as the superposition of a Fleming-Viot diffusion with a deterministic volatility constant $d_k$ and a Cannings process with a random jump rate, both depending on $k$. In the random environment, $d_k$ turns out to be smaller than in the homogeneous environment of the same mean. We investigate how $d_k$ scales with $k$. This leads to five universality classes of cluster formation in the mono-type regime. We find that if clustering occurs, then the random environment slows down the growth of the clusters, i.e., enhances the diversity of types.
We study the evolution of genealogies of a population of individuals, whose type frequencies result in an interacting Fleming-Viot process on $\Z$. We construct and analyze the genealogical structure of the population in this genealogy-valued Fleming-Viot process as a marked metric measure space, with each individual carrying its spatial location as a mark. We then show that its time evolution converges to that of the genealogy of a continuum-sites stepping stone model on $\R$, if space and time are scaled diffusively. We construct the genealogies of the continuum-sites stepping stone model as functionals of the Brownian web, and furthermore, we show that its evolution solves a martingale problem. The generator for the continuum-sites stepping stone model has a singular feature: at each time, the resampling of genealogies only affects a set of individuals of measure $0$. Along the way, we prove some negative correlation inequalities for coalescing Brownian motions, as well as extend the theory of marked metric measure spaces (developed recently by Depperschmidt, Greven and Pfaffelhuber [DGP12]) from the case of probability measures to measures that are finite on bounded sets.
For a beneficial allele which enters a large unstructured population and eventually goes to fixation, it is known that the time to fixation is approximately $2\log(\alpha)/\alpha$ for a large selection coefficient $\alpha$. For a population that is distributed over finitely many colonies, with migration between these colonies, we detect various regimes of the migration rate $\mu$ for which the fixation times have different asymptotics as $\alpha \to \infty$. If $\mu$ is of order $\alpha$, the allele fixes (as in the spatially unstructured case) in time $\sim 2\log(\alpha)/\alpha$. If $\mu$ is of order $\alpha^\gamma, 0\leq \gamma \leq 1$, the fixation time is $\sim (2 + (1-\gamma)\Delta) \log(\alpha)/\alpha$, where $\Delta$ is the number of migration steps that are needed to reach all other colonies starting from the colony where the beneficial allele appeared. If $\mu = 1/\log(\alpha)$, the fixation time is $\sim (2+S)\log(\alpha)/\alpha$, where $S$ is a random time in a simple epidemic model. The main idea for our analysis is to combine a new moment dual for the process conditioned to fixation with the time reversal in equilibrium of a spatial version of Neuhauser and Krone's ancestral selection graph.
We consider a spatial multi-type branching model in which individuals migrate in geographic space according to random walks and reproduce according to a state-dependent branching mechanism which can be sub-, super- or critical depending on the local intensity of individuals of the different types. The model is a Lotka-Volterra type model with a spatial component and is related to two models studied in as well as to earlier work in and in . Our main focus is on the diffusion limit of small mass, locally many individuals and rapid reproduction. This system differs from spatial critical branching systems since it is not density preserving and the densities for large times do not depend on the initial distribution but mainly on the carrying capacities. We prove existence of the infinite particle model and the system of interacting diffusions as solutions of martingale problems or systems of stochastic equations. In the exchangeable case in which the parameters are not type dependent we show uniqueness of the solutions. For that purpose we establish a new exponential duality.