BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Nicolas Champagnat

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

ResearcherUnclaimed source record
math.PRPopulations and Evolutionmath.APmath.CAmath.DSmath.NAmath.OC

Catalog footprint

What is connected

28works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

28 published item(s)

preprint2025arXiv

Quasi-stationary distributions in reducible state spaces

We study quasi-stationary distributions and quasi-limiting behavior of Markov chains in general reducible state spaces with absorption. We propose a set of assumptions dealing with particular situations where the state space can be decomposed into three subsets between which communication is only possible in a single direction. These assumptions allow us to characterize the exponential order of magnitude and the exact polynomial correction, called polynomial convergence parameter, for the leading order term of the semigroup for large time. They also provide explicit convergence speeds to this leading order term. We apply these results to general Markov chains with finitely or denumerably many communication classes using a specific induction over the communication classes of the chain. We are able to explicitely characterize the polynomial convergence parameter, to determine the complete set of quasistationary distributions and to provide explicit estimates for the speed of convergence to quasi-limiting distributions in the case of finitely many communication classes. We conclude with an application of these results to the case of denumerable state spaces, where we are able to prove that, in general, there is existence of a quasi-stationary distribution without assuming irreducibility before absorption. This actually holds true assuming only aperiodicity, the existence of a Lyapunov function and the existence of a point in the state space from which the return time is finite with positive probability.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Filling the gap between individual-based evolutionary models and Hamilton-Jacobi equations

We consider a stochastic model for the evolution of a discrete population structured by a trait with values on a finite grid of the torus, and with mutation and selection. Traits are vertically inherited unless a mutation occurs, and influence the birth and death rates. We focus on a parameter scaling where population is large, individual mutations are small but not rare, and the grid mesh for the trait values is much smaller than the size of mutation steps. When considering the evolution of the population in a long time scale, the contribution of small sub-populations may strongly influence the dynamics. Our main result quantifies the asymptotic dynamics of sub-population sizes on a logarithmic scale. We establish that under the parameter scaling the logarithm of the stochastic population size process, conveniently normalized, converges to the unique viscosity solution of a Hamilton-Jacobi equation. Such Hamilton-Jacobi equations have already been derived from parabolic integro-differential equations and have been widely developed in the study of adaptation of quantitative traits. Our work provides a justification of this framework directly from a stochastic individual based model, leading to a better understanding of the results obtained within this approach. The proof makes use of almost sure maximum principles and careful controls of the martingale parts.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Stochastic approximation of quasi-stationary distributions for diffusion processes in a bounded domain

We study a random process with reinforcement, which evolves following the dynamics of a given diffusion process in a bounded domain and is resampled according to its occupation measure when it reaches the boundary. We show that its occupation measure converges to the unique quasi-stationary distribution of the diffusion process absorbed at the boundary of the domain. Our proofs use recent results in the theory of quasi-stationary distributions and stochastic approximation techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Moments of the frequency spectrum of a splitting tree with neutral Poissonian mutations

We consider a branching population where individuals live and reproduce independently. Their lifetimes are i.i.d. and they give birth at a constant rate b. The genealogical tree spanned by this process is called a splitting tree, and the population counting process is a homogeneous, binary Crump-Mode-Jagers process. We suppose that mutations affect individuals independently at a constant rate $\θ$ during their lifetimes, under the infinite-alleles assumption: each new mutation gives a new type, called allele, to his carrier. We study the allele frequency spectrum which is the numbers A(k, t) of types represented by k alive individuals in the population at time t. Thanks to a new construction of the coalescent point process describing the genealogy of individuals in the splitting tree, we are able to compute recursively all joint factorial moments of (A(k, t)) k$\ge$1. These moments allow us to give an elementary proof of the almost sure convergence of the frequency spectrum in a supercritical splitting tree.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the link between infinite horizon control and quasi-stationary distributions

We study infinite horizon control of continuous-time non-linear branching processes with almost sure extinction for general (positive or negative) discount. Our main goal is to study the link between infinite horizon control of these processes and an optimization problem involving their quasi-stationary distributions and the corresponding extinction rates. More precisely, we obtain an equivalent of the value function when the discount parameter is close to the threshold where the value function becomes infinite, and we characterize the optimal Markov control in this limit. To achieve this, we present a new proof of the dynamic programming principle based upon a pseudo-Markov property for controlled jump processes. We also prove the convergence to a unique quasi-stationary distribution of non-linear branching processes controlled by a Markov control conditioned on non-extinction.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Population processes with unbounded extinction rate conditioned to non-extinction

This article studies the quasi-stationary behaviour of population processes with unbounded absorption rate, including one-dimensional birth and death processes with catastrophes and multi-dimensional birth and death processes, modeling biological populations in interaction. To handle this situation, we develop original non-linear Lyapunov criteria. We obtain the exponential convergence in total variation of the conditional distributions to a unique quasi-stationary distribution, uniformly with respect to the initial distribution. Our results cover all one-dimensional birth and death processes which come down from infinity with catastrophe rate satisfying appropriate bounds, and multi-dimensional birth and death models with stronger intra-specific than inter-specific competition.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

From stochastic, individual-based models to the canonical equation of adaptive dynamics - In one step

We consider a model for Darwinian evolution in an asexual population with a large but non-constant populations size characterized by a natural birth rate, a logistic death rate modelling competition and a probability of mutation at each birth event. In the present paper, we study the long-term behavior of the system in the limit of large population $(K\to \infty)$ size, rare mutations $(u\to 0)$, and small mutational effects $(σ\to 0)$, proving convergence to the canonical equation of adaptive dynamics (CEAD). In contrast to earlier works, e.g. by Champagnat and Méléard, we take the three limits simultaneously, i.e. $u=u_K$ and $σ=σ_K$, tend to zero with $K$, subject to conditions that ensure that the time-scale of birth and death events remains separated from that of successful mutational events. This slows down the dynamics of the microscopic system and leads to serious technical difficulties that requires the use of completely different methods. In particular, we cannot use the law of large numbers on the diverging time needed for fixation to approximate the stochastic system with the corresponding deterministic one. To solve this problem we develop a "stochastic Euler scheme" based on coupling arguments that allows to control the time evolution of the stochastic system over time-scales that diverge with $K$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Links between deterministic and stochastic approaches for invasion in growth-fragmentation-death models

We present two approaches to study invasion in growth-fragmentation-death mod- els. The first one is based on a stochastic individual based model, which is a piecewise deterministic branching process with a continuum of types, and the second one is based on an integro-differential model. The invasion of the population is described by the survival probability for the former model and by an eigenproblem for the latter one. We study these two notions of invasion fitness, giving different characterizations of the growth of the population, and we make links between these two complementary points of view. In particular we prove that the two approaches lead to the same crite- rion of possible invasion. Based on Krein-Rutman theory, we also give a proof of the existence of a solution to the eigenproblem, which satisfies the conditions needed for our study of the stochastic model, hence providing a set of assumptions under which both approaches can be carried out. Finally, we motivate our work in the context of adaptive dynamics in a chemostat model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Monte Carlo methods for linear and non-linear Poisson-Boltzmann equation

The electrostatic potential in the neighborhood of a biomolecule can be computed thanks to the non-linear divergence-form elliptic Poisson-Boltzmann PDE. Dedicated Monte-Carlo methods have been developed to solve its linearized version (see e.g.Bossy et al 2009, Mascagni & Simonov 2004}). These algorithms combine walk on spheres techniques and appropriate replacements at the boundary of the molecule. In the first part of this article we compare recent replacement methods for this linearized equation on real size biomolecules, that also require efficient computational geometry algorithms. We compare our results with the deterministic solver APBS. In the second part, we prove a new probabilistic interpretation of the nonlinear Poisson-Boltzmann PDE. A Monte Carlo algorithm is also derived and tested on a simple test case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quasi-stationary distribution for multi-dimensional birth and death processes conditioned to survival of all coordinates

This article studies the quasi-stationary behaviour of multidimensional birth and death processes, modeling the interaction between several species, absorbed when one of the coordinates hits 0. We study models where the absorption rate is not uniformly bounded, contrary to most of the previous works. To handle this natural situation, we develop original Lyapunov function arguments that might apply in other situations with unbounded killing rates. We obtain the exponential convergence in total variation of the conditional distributions to a unique stationary distribution, uniformly with respect to the initial distribution. Our results cover general birth and death models with stronger intra-specific than inter-specific competition, and cases with neutral competition with explicit conditions on the dimension of the process.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exponential convergence to quasi-stationary distribution and Q-process

For general, almost surely absorbed Markov processes, we obtain necessary and sufficient conditions for exponential convergence to a unique quasi-stationary distribution in the total variation norm. These conditions also ensure the existence and exponential ergodicity of the $Q$-process (the process conditioned to never be absorbed). We apply these results to one-dimensional birth and death processes with catastrophes, multi-dimensional birth and death processes, infinite-dimensional population models with Brownian mutations and neutron transport dynamics absorbed at the boundary of a bounded domain.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Adaptation in a stochastic multi-resources chemostat model

We are interested in modeling the Darwinian evolution resulting from the interplay of phenotypic variation and natural selection through ecological interactions, in the specific scales of the biological framework of adaptive dynamics. Adaptive dynamics so far has been put on a rigorous footing only for direct competition models (Lotka-Volterra models) involving a competition kernel which describes the competition pressure from one individual to another one. We extend this to a multi-resources chemostat model, where the competition between individuals results from the sharing of several resources which have their own dynamics. Starting from a stochastic birth and death process model, we prove that, when advantageous mutations are rare, the population behaves on the mutational time scale as a jump process moving between equilibrium states (the polymorphic evolution sequence of the adaptive dynamics literature). An essential technical ingredient is the study of the long time behavior of a chemostat multi-resources dynamical system. In the small mutational steps limit this process in turn gives rise to a differential equation in phenotype space called canonical equation of adaptive dynamics. From this canonical equation and still assuming small mutation steps, we prove a rigorous characterization of the evolutionary branching points.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Strong solutions to stochastic differential equations with rough coefficients

We study strong existence and pathwise uniqueness for stochastic differential equations in $\RR^d$ with rough coefficients, and without assuming uniform ellipticity for the diffusion matrix. Our approach relies on direct quantitative estimates on solutions to the SDE, assuming Sobolev bounds on the drift and diffusion coefficients, and $L^p$ bounds for the solution of the corresponding Fokker-Planck PDE, which can be proved separately. This allows a great flexibility regarding the method employed to obtain these last bounds. Hence we are able to obtain general criteria in various cases, including the uniformly elliptic case in any dimension, the one-dimensional case and the Langevin (kinetic) case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Birth and death processes with neutral mutations

In this paper, we review recent results of ours concerning branching processes with general lifetimes and neutral mutations, under the infinitely many alleles model, where mutations can occur either at birth of individuals or at a constant rate during their lives. In both models, we study the allelic partition of the population at time t. We give closed formulae for the expected frequency spectrum at t and prove pathwise convergence to an explicit limit, as t goes to infinity, of the relative numbers of types younger than some given age and carried by a given number of individuals (small families). We also provide convergences in distribution of the sizes or ages of the largest families and of the oldest families. In the case of exponential lifetimes, population dynamics are given by linear birth and death processes, and we can most of the time provide general formulations of our results unifying both models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On Dirichlet eigenvectors for neutral two-dimensional Markov chains

We consider a general class of discrete, two-dimensional Markov chains modeling the dynamics of a population with two types, without mutation or immigration, and neutral in the sense that type has no influence on each individual's birth or death parameters. We prove that all the eigenvectors of the corresponding transition matrix or infinitesimal generator Π can be expressed as the product of "universal" polynomials of two variables, depending on each type's size but not on the specific transitions of the dynamics, and functions depending only on the total population size. These eigenvectors appear to be Dirichlet eigenvectors for Π on the complement of triangular subdomains, and as a consequence the corresponding eigenvalues are ordered in a specific way. As an application, we study the quasistationary behavior of finite, nearly neutral, two-dimensional Markov chains, absorbed in the sense that 0 is an absorbing state for each component of the process.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Splitting trees with neutral Poissonian mutations II: Largest and Oldest families

We consider a supercritical branching population, where individuals have i.i.d. lifetime durations (which are not necessarily exponentially distributed) and give birth (singly) at constant rate. We assume that individuals independently experience neutral mutations, at constant rate $θ$ during their lifetimes, under the infinite-alleles assumption: each mutation instantaneously confers a brand new type, called allele or haplotype, to its carrier. The type carried by a mother at the time when she gives birth is transmitted to the newborn. We are interested in the sizes and ages at time $t$ of the clonal families carrying the most abundant alleles or the oldest ones, as $t\to\infty$, on the survival event. Intuitively, the results must depend on how the mutation rate $θ$ and the Malthusian parameter $α>0$ compare. Hereafter, $N\equiv N_t$ is the population size at time $t$, constants $a,c$ are scaling constants, whereas $k,k'$ are explicit positive constants which depend on the parameters of the model. When $α>θ$, the most abundant families are also the oldest ones, they have size $cN^{1-θ/α}$ and age $t-a$. When $α<θ$, the oldest families have age $(α/θ)t+a$ and tight sizes; the most abundant families have sizes $k\log(N)-k'\log\log(N)+c$ and all have age $(θ-α)^{-1}\log(t)$. When $α=θ$, the oldest families have age $kt-k'\log(t)+a$ and tight sizes; the most abundant families have sizes $(k\log(N)-k'\log\log(N)+c)^2$ and all have age $t/2$. Those informal results can be stated rigorously in expectation. Relying heavily on the theory of coalescent point processes, we are also able, when $α\leqθ$, to show convergence in distribution of the joint, properly scaled ages and sizes of the most abundant/oldest families and to specify the limits as some explicit Cox processes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reconciling alternate methods for the determination of charge distributions: A probabilistic approach to high-dimensional least-squares approximations

We propose extensions and improvements of the statistical analysis of distributed multipoles (SADM) algorithm put forth by Chipot et al. in [6] for the derivation of distributed atomic multipoles from the quantum-mechanical electrostatic potential. The method is mathematically extended to general least-squares problems and provides an alternative approximation method in cases where the original least-squares problem is computationally not tractable, either because of its ill-posedness or its high-dimensionality. The solution is approximated employing a Monte Carlo method that takes the average of a random variable defined as the solutions of random small least-squares problems drawn as subsystems of the original problem. The conditions that ensure convergence and consistency of the method are discussed, along with an analysis of the computational cost in specific instances.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Splitting trees with neutral Poissonian mutations I: Small families

We consider a neutral dynamical model of biological diversity, where individuals live and reproduce independently. They have i.i.d. lifetime durations (which are not necessarily exponentially distributed) and give birth (singly) at constant rate b. Such a genealogical tree is usually called a splitting tree, and the population counting process (N_t;t\ge 0) is a homogeneous, binary Crump--Mode--Jagers process. We assume that individuals independently experience mutations at constant rate θduring their lifetimes, under the infinite-alleles assumption: each mutation instantaneously confers a brand new type, called allele, to its carrier. We are interested in the allele frequency spectrum at time t, i.e., the number A(t) of distinct alleles represented in the population at time t, and more specifically, the numbers A(k,t) of alleles represented by k individuals at time t, k=1,2,...,N_t. We mainly use two classes of tools: coalescent point processes and branching processes counted by random characteristics. We provide explicit formulae for the expectation of A(k,t) in a coalescent point process conditional on population size, which apply to the special case of splitting trees. We separately derive the a.s. limits of A(k,t)/N_t and of A(t)/N_t thanks to random characteristics. Last, we separately compute the expected homozygosity by applying a method characterizing the dynamics of the tree distribution as the origination time of the tree moves back in time, in the spirit of backward Kolmogorov equations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The evolutionary limit for models of populations interacting competitively with many resources

We consider a integro-differential nonlinear model that describes the evolution of a population structured by a quantitative trait. The interactions between traits occur from competition for resources whose concentrations depend on the current state of the population. Following the formalism of\cite{DJMP}, we study a concentration phenomenon arising in the limit of strong selection and small mutations. We prove that the population density converges to a sum of Dirac masses characterized by the solution $ϕ$ of a Hamilton-Jacobi equation which depends on resource concentrations that we fully characterize in terms of the function $ϕ$.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Limit theorems for conditioned multitype Dawson-Watanabe processes and Feller diffusions

A multitype Dawson-Watanabe process is conditioned, in subcritical and critical cases, on non-extinction in the remote future. On every finite time interval, its distribution is absolutely continuous with respect to the law of the unconditioned process. A martingale problem characterization is also given. Several results on the long time behavior of the conditioned mass process|the conditioned multitype Feller branching diffusion are then proved. The general case is first considered, where the mutation matrix which models the interaction between the types, is irreducible. Several two-type models with decomposable mutation matrices are analyzed too.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Polymorphic evolution sequence and evolutionary branching

We are interested in the study of models describing the evolution of a polymorphic population with mutation and selection in the specific scales of the biological framework of adaptive dynamics. The population size is assumed to be large and the mutation rate small. We prove that under a good combination of these two scales, the population process is approximated in the long time scale of mutations by a Markov pure jump process describing the successive trait equilibria of the population. This process, which generalizes the so-called trait substitution sequence, is called polymorphic evolution sequence. Then we introduce a scaling of the size of mutations and we study the polymorphic evolution sequence in the limit of small mutations. From this study in the neighborhood of evolutionary singularities, we obtain a full mathematical justification of a heuristic criterion for the phenomenon of evolutionary branching. To this end we finely analyze the asymptotic behavior of 3-dimensional competitive Lotka-Volterra systems.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Adaptive dynamics in logistic branching populations

We consider a trait-structured population subject to mutation, birth and competition of logistic type, where the number of coexisting types may fluctuate. Applying a limit of rare mutations to this population while keeping the population size finite leads to a jump process, the so-called `trait substitution sequence', where evolution proceeds by successive invasions and fixations of mutant types. The probability of fixation of a mutant is interpreted as a fitness landscape that depends on the current state of the population. It was in adaptive dynamics that this kind of model was first invented and studied, under the additional assumption of large population. Assuming also small mutation steps, adaptive dynamics' theory provides a deterministic ODE approximating the evolutionary dynamics of the dominant trait of the population, called `canonical equation of adaptive dynamics'. In this work, we want to include genetic drift in this models by keeping the population finite. Rescaling mutation steps (weak selection) yields in this case a diffusion on the trait space that we call `canonical diffusion of adaptive dynamics', in which genetic drift (diffusive term) is combined with directional selection (deterministic term) driven by the fitness gradient. Finally, in order to compute the coefficients of this diffusion, we seek explicit first-order formulae for the probability of fixation of a nearly neutral mutant appearing in a resident population. These formulae are expressed in terms of `invasibility coefficients' associated with fertility, defense, aggressiveness and isolation, which measure the robustness (stability w.r.t. selective strengths) of the resident type. Some numerical results on the canonical diffusion are also given.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Invasion and adaptive evolution for individual-based spatially structured populations

The interplay between space and evolution is an important issue in population dynamics, that is in particular crucial in the emergence of polymorphism and spatial patterns. Recently, biological studies suggest that invasion and evolution are closely related. Here we model the interplay between space and evolution starting with an individual-based approach and show the important role of parameter scalings on clustering and invasion. We consider a stochastic discrete model with birth, death, competition, mutation and spatial diffusion, where all the parameters may depend both on the position and on the trait of individuals. The spatial motion is driven by a reflected diffusion in a bounded domain. The interaction is modelled as a trait competition between individuals within a given spatial interaction range. First, we give an algorithmic construction of the process. Next, we obtain large population approximations, as weak solutions of nonlinear reaction-diffusion equations with Neumann's boundary conditions. As the spatial interaction range is fixed, the nonlinearity is nonlocal. Then, we make the interaction range decrease to zero and prove the convergence to spatially localized nonlinear reaction-diffusion equations, with Neumann's boundary conditions. Finally, simulations based on the microscopic individual-based model are given, illustrating the strong effects of the spatial interaction range on the emergence of spatial and phenotypic diversity (clustering and polymorphism) and on the interplay between invasion and evolution. The simulations focus on the qualitative differences between local and nonlocal interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Individual-based probabilistic models of adaptive evolution and various scaling approximations

We are interested in modelling Darwinian evolution, resulting from the interplay of phenotypic variation and natural selection through ecological interactions. Our models are rooted in the microscopic, stochastic description of a population of discrete individuals characterized by one or several adaptive traits. The population is modelled as a stochastic point process whose generator captures the probabilistic dynamics over continuous time of birth, mutation, and death, as influenced by each individual's trait values, and interactions between individuals. An offspring usually inherits the trait values of her progenitor, except when a mutation causes the offspring to take an instantaneous mutation step at birth to new trait values. We look for tractable large population approximations. By combining various scalings on population size, birth and death rates, mutation rate, mutation step, or time, a single microscopic model is shown to lead to contrasting macroscopic limits, of different nature: deterministic, in the form of ordinary, integro-, or partial differential equations, or probabilistic, like stochastic partial differential equations or superprocesses. In the limit of rare mutations, we show that a possible approximation is a jump process, justifying rigorously the so-called trait substitution sequence. We thus unify different points of view concerning mutation-selection evolutionary models.

0 citations0 reviewsNo code linkedOpen access