-
Epidemic models with varying infectivity on a refining spatial grid. I. The SI model
Authors:
Anicet Mougabe-Peurkor,
Étienne Pardoux,
Ténan Yeo
Abstract:
We consider a space-time SI epidemic model with infection age-dependent infectivity and non-local infections constructed on a grid of the torus $\mathbb{T}^1 =(0, 1]^d$, where the individuals may migrate from node to another. The migration processes in either of the two states are assumed to be Markovian. We establish a functional law of large numbers by letting jointly $N$ the initial approximate…
▽ More
We consider a space-time SI epidemic model with infection age-dependent infectivity and non-local infections constructed on a grid of the torus $\mathbb{T}^1 =(0, 1]^d$, where the individuals may migrate from node to another. The migration processes in either of the two states are assumed to be Markovian. We establish a functional law of large numbers by letting jointly $N$ the initial approximate number of individuals on each node go to infinity and $\varepsilon$ the mesh size of the grid go to zero. The limit is a system of parabolic PDE/integral equations. The constraint on the speed of convergence of the parameters $N$ and $\varepsilon$ is that $N\varepsilon^d \to \infty$ as $(N, \varepsilon)\to (+\infty, 0)$.
△ Less
Submitted 4 June, 2023;
originally announced June 2023.
-
Spatially dense stochastic epidemic models with infection-age dependent infectivity
Authors:
Guodong Pang,
Etienne Pardoux
Abstract:
We study an individual-based stochastic spatial epidemic model where the number of locations and the number of individuals at each location both grow to infinity. Each individual is associated with a random infection-age dependent infectivity function. Individuals are infected through interactions across the locations with heterogeneous effects. The epidemic dynamics can be described using a time-…
▽ More
We study an individual-based stochastic spatial epidemic model where the number of locations and the number of individuals at each location both grow to infinity. Each individual is associated with a random infection-age dependent infectivity function. Individuals are infected through interactions across the locations with heterogeneous effects. The epidemic dynamics can be described using a time-space representation for the the total force of infection, the number of susceptible individuals, the number of infected individuals that are infected at each time and have been infected for a certain amount of time, as well as the number of recovered individuals. We prove a functional law of large numbers for these time-space processes, and in the limit, we obtain a set of time-space integral equations. We then derive the PDE models from the limiting time-space integral equations, in particular, the density (with respect to the infection age) of the time-age-space integral equation for the number of infected individuals tracking the age of infection satisfies a linear PDE in time and age with an integral boundary condition. These integral equation and PDE limits can be regarded as dynamics on graphon under certain conditions.
△ Less
Submitted 12 September, 2023; v1 submitted 11 April, 2023;
originally announced April 2023.
-
PDE model for multi-patch epidemic models with migration and infection-age dependent infectivity
Authors:
Guodong Pang,
Etienne Pardoux
Abstract:
We study a stochastic epidemic model with multiple patches (locations), where individuals in each patch are categorized into three compartments, Susceptible, Infected and Recovered/Removed, and may migrate from one patch to another in any of the compartments. Each individual is associated with a random infectivity function which dictates the force of infectivity while the interactive infection pro…
▽ More
We study a stochastic epidemic model with multiple patches (locations), where individuals in each patch are categorized into three compartments, Susceptible, Infected and Recovered/Removed, and may migrate from one patch to another in any of the compartments. Each individual is associated with a random infectivity function which dictates the force of infectivity while the interactive infection process depends on the age of infection (elapsed time since infection). We prove a functional law of large number for the epidemic evolution dynamics including the aggregate infectivity process, the numbers of susceptible and recovered individuals as well as the number of infected individuals at each time that have been infected for a certain amount of time. From the limits, we derive a PDE model for the density of the number of infected individuals with respect to the infection age, which is a systems of linear PDE equations with a boundary condition that is determined by a set of integral equations.
△ Less
Submitted 18 August, 2023; v1 submitted 6 April, 2023;
originally announced April 2023.
-
Extinction time of an epidemic with infection age dependent infectivity
Authors:
Anicet Mougabe-Peurkor,
Ibrahima Drame,
Modeste N'zi,
Etienne Pardoux
Abstract:
This paper studies the distribution function of the time of extinction of a subcritical epidemic, when a large enough proportion of the population has been immunized and/or the infectivity of the infectious individuals has been reduced, so that the effective reproduction number is less than one. We do that for a SIR/SEIR model, where infectious individuals have an infection age dependent infectivi…
▽ More
This paper studies the distribution function of the time of extinction of a subcritical epidemic, when a large enough proportion of the population has been immunized and/or the infectivity of the infectious individuals has been reduced, so that the effective reproduction number is less than one. We do that for a SIR/SEIR model, where infectious individuals have an infection age dependent infectivity, as in the model introduced in the 1927 seminal paper of Kermack and McKendrick. Our main conclusion is that simplifying the model as an ODE SIR model, as it is largely done in the epidemics literature, introduces a biais toward shorter extinction time.
△ Less
Submitted 3 April, 2023;
originally announced April 2023.
-
Central limit theorem for a spatial stochastic epidemic model with mean field interaction
Authors:
Maxime Hauray,
Etienne Pardoux,
Yen V. Vuong
Abstract:
In this article, we study an interacting particle system in the context of epidemiology where the individuals (particles) are characterized by their position and infection state. We begin with a description at the microscopic level where the displacement of individuals is driven by mean field interactions and state-dependent diffusion, whereas the epidemiological dynamic is described by the Poisso…
▽ More
In this article, we study an interacting particle system in the context of epidemiology where the individuals (particles) are characterized by their position and infection state. We begin with a description at the microscopic level where the displacement of individuals is driven by mean field interactions and state-dependent diffusion, whereas the epidemiological dynamic is described by the Poisson processes with an infection rate based on the distribution of other nearby individuals, also of the mean-field type. Then under suitable assumptions, a form of law of large numbers has been established to show that the associated empirical measure to the above system converges to the law of the unique solution of a nonlinear McKean-Vlasov equation. As a natural follow-up question, we study the fluctuation of this stochastic system around its limit. We prove that this fluctuation process converges to a limit process, which can be characterized as the unique solution of a linear stochastic PDE. Unlike the existing literature using a coupling approach to prove the central limit theorem for interacting particle systems, the main idea in our proof is to use a semigroup formalism and some appropriate estimates to directly study the linearized evolution equation of the fluctuation process in a suitable weighted Sobolev space and follows a Hilbertian approach.
△ Less
Submitted 5 December, 2022;
originally announced December 2022.
-
Stochastic epidemic models with varying infectivity and susceptibility
Authors:
Raphaël Forien,
Guodong Pang,
Étienne Pardoux,
Arsene Brice Zotsa--Ngoufack
Abstract:
We study an individual-based stochastic epidemic model in which infected individuals become susceptible again following each infection. In contrast to classical compartment models, after each infection, the infectivity is a random function of the time elapsed since one's infection. Similarly, recovered individuals become gradually susceptible after some time according to a random susceptibility fu…
▽ More
We study an individual-based stochastic epidemic model in which infected individuals become susceptible again following each infection. In contrast to classical compartment models, after each infection, the infectivity is a random function of the time elapsed since one's infection. Similarly, recovered individuals become gradually susceptible after some time according to a random susceptibility function. We study the large population asymptotic behaviour of the model, by proving a functional law of large numbers (FLLN) and investigating the endemic equilibria properties of the limit. The limit depends on the law of the susceptibility random functions but only on the mean infectivity functions. The FLLN is proved by constructing a sequence of i.i.d. auxiliary processes and adapting the approach from the theory of propagation of chaos. The limit is a generalisation of a PDE model introduced by Kermack and McKendrick, and we show how this PDE model can be obtained as a special case of our FLLN limit.% for a particular set of infectivity and susceptibility random functions and initial conditions. For the endemic equilibria, if $ R_0 $ is lower than (or equal to) some threshold, the epidemic does not last forever and eventually disappears from the population, while if $ R_0 $ is larger than this threshold, the epidemic will not disappear and there exists an endemic equilibrium. The value of this threshold turns out to depend on the harmonic mean of the susceptibility a long time after an infection, a fact which was not previously known.
△ Less
Submitted 10 October, 2022;
originally announced October 2022.
-
Multi-patch multi-group epidemic model with varying infectivity
Authors:
Raphaël Forien,
Guodong Pang,
Étienne Pardoux
Abstract:
This paper presents a law of large numbers result, as the size of the population tends to infinity, of SIR stochastic epidemic models, for a population distributed over $L$ distinct patches (with migrations between them) and $K$ distinct groups (possibly age groups). The limit is a set of Volterra-type integral equations, and the result shows the effects of both spatial and population heterogeneit…
▽ More
This paper presents a law of large numbers result, as the size of the population tends to infinity, of SIR stochastic epidemic models, for a population distributed over $L$ distinct patches (with migrations between them) and $K$ distinct groups (possibly age groups). The limit is a set of Volterra-type integral equations, and the result shows the effects of both spatial and population heterogeneity. The novelty of the model is that the infectivity of an infected individual is infection age dependent. More precisely, to each infected individual is attached a random infection-age dependent infectivity function, such that the various random functions attached to distinct individuals are i.i.d.
The proof involves a novel construction of a sequence of i.i.d. processes to invoke the law of large numbers for processes in $D$, by using the solution of a MacKean-Vlasov type Poisson-driven stochastic equation (as in the propagation of chaos theory). We also establish an identity using the Feynman-Kac formula for an adjoint backward ODE. The advantage of this approach is that it assumes much weaker conditions on the random infectivity functions than our earlier work for the homogeneous model in [20], where standard tightness criteria for convergence of stochastic processes were employed. To illustrate this new approach, we first explain the new proof under the weak assumptions for the homogeneous model, and then describe the multipatch-multigroup model and prove the law of large numbers for that model.
△ Less
Submitted 23 June, 2022; v1 submitted 11 November, 2021;
originally announced November 2021.
-
Conditional propagation of chaos in a spatial stochastic epidemic model with common noise
Authors:
Yen V. Vuong,
Maxime Hauray,
Etienne Pardoux
Abstract:
We study a stochastic spatial epidemic model where the $N$ individuals carry two features: a position and an infection state, interact and move in $\R^d$. In this Markovian model, the evolution of the infection states are described with the help of the Poisson Point Processes , whereas the displacement of the individuals are driven by mean field advection, a (state dependence) diffusion and also a…
▽ More
We study a stochastic spatial epidemic model where the $N$ individuals carry two features: a position and an infection state, interact and move in $\R^d$. In this Markovian model, the evolution of the infection states are described with the help of the Poisson Point Processes , whereas the displacement of the individuals are driven by mean field advection, a (state dependence) diffusion and also a common noise, so that the spatial dynamic is a random process. We prove that when the number $N$ of individual goes to infinity, the conditional propagation of chaos holds : conditionnally to the common noise, the individuals are asymptotically independent and the stochastic dynamic converges to a "random" nonlinear McKean-Vlasov process. As a consequence, the associated empirical measure converges to a measure, which is solution of a stochastic mean-field PDE driven by the common noise.
△ Less
Submitted 4 November, 2021;
originally announced November 2021.
-
Recent Advances in Epidemic Modeling: Non-Markov Stochastic Models and their Scaling Limits
Authors:
Raphael Forien,
Guodong Pang,
Etienne Pardoux
Abstract:
In this survey paper, we review the recent advances in individual based non--Markovian epidemic models. They include epidemic models with a constant infectivity rate, varying infectivity rate or infection-age dependent infectivity, infection-age recovery rate (or equivalently, general law of infectious period), as well as varying susceptibility/immunity. We focus on the scaling limits with a large…
▽ More
In this survey paper, we review the recent advances in individual based non--Markovian epidemic models. They include epidemic models with a constant infectivity rate, varying infectivity rate or infection-age dependent infectivity, infection-age recovery rate (or equivalently, general law of infectious period), as well as varying susceptibility/immunity. We focus on the scaling limits with a large population, functional law of large numbers (FLLN) and functional central limit theorems (FCLT), while the large and moderate deviations for some Markovian epidemic models are also reviewed. In the FLLN, the limits are a set of Volterra integral equations, and in the FCLT, the limits are stochastic Volterra integral equations driven by Gaussian processes. We relate our deterministic limits to the results in the seminal papers by Kermack and McKendrick published in 1927, 1932 and 1933, where the varying infectivity and susceptibility/immunity were already considered. We also discuss some extensions, including models with heterogeneous population, spatial models and control problems, as well as open problems.
△ Less
Submitted 22 September, 2022; v1 submitted 15 June, 2021;
originally announced June 2021.
-
Functional law of large numbers and PDEs for epidemic models with infection-age dependent infectivity
Authors:
Guodong Pang,
Etienne Pardoux
Abstract:
We study epidemic models where the infectivity of each individual is a random function of the infection age (the elapsed time of infection). To describe the epidemic evolution dynamics, we use a stochastic process that tracks the number of individuals at each time that have been infected for less than or equal to a certain amount of time, together with the aggregate infectivity process. We establi…
▽ More
We study epidemic models where the infectivity of each individual is a random function of the infection age (the elapsed time of infection). To describe the epidemic evolution dynamics, we use a stochastic process that tracks the number of individuals at each time that have been infected for less than or equal to a certain amount of time, together with the aggregate infectivity process. We establish the functional law of large numbers (FLLN) for the stochastic processes that describe the epidemic dynamics. The limits are described by a set of deterministic integral equations, which has a further characterization using PDEs under some regularity conditions. The solutions are characterized with boundary conditions that are given by a system of Volterra equations. We also characterize the equilibrium points for the PDEs in the SIS model with infection-age dependent infectivity. To establish the FLLNs, we employ a useful criterion for weak convergence for the two-parameter processes together with useful representations for the relevant processes via Poisson random measures.
△ Less
Submitted 28 December, 2022; v1 submitted 7 June, 2021;
originally announced June 2021.
-
Functional central limit theorems for epidemic models with varying infectivity
Authors:
Guodong Pang,
Etienne Pardoux
Abstract:
In this paper, we prove functional central limit theorems (FCLTs) for a stochastic epidemic model with varying infectivity and general infectious periods recently introduced in Forien, Pang and Pardoux (2020).The infectivity process (total force of infection at each time) is composed of the independent infectivity random functions of each infectious individual, which starts at the time of infectio…
▽ More
In this paper, we prove functional central limit theorems (FCLTs) for a stochastic epidemic model with varying infectivity and general infectious periods recently introduced in Forien, Pang and Pardoux (2020).The infectivity process (total force of infection at each time) is composed of the independent infectivity random functions of each infectious individual, which starts at the time of infection. These infectivity random functions induce the infectious periods (as well as exposed, recovered or immune periods in full generality), whose probability distributions can be very general. The epidemic model includes the generalized non-Markovian SIR, SEIR, SIS, SIRS models with infection-age dependent infectivity. In the FCLTs for the generalized SIR and SEIR models, the limits of the diffusion-scaled fluctuations of the infectivity and susceptible processes are a unique solution to a two-dimensional Gaussian-driven stochastic Volterra integral equations, and then given these solutions, the limits for the infected (exposed/infectious) and recovered processes are Gaussian processes expressed in terms of the solutions to those stochastic Volterra integral equations. We also present the FCLTs for the generalized SIS and SIRS models.
△ Less
Submitted 11 September, 2022; v1 submitted 24 September, 2020;
originally announced September 2020.
-
Metastability between the clicks of Muller's ratchet
Authors:
Mauro Mariani,
Etienne Pardoux,
Aurélien Velleret
Abstract:
We prove the existence and uniqueness of a quasi-stationary distribution for three stochastic processes derived from the model of Muller's ratchet. This model was invented with the aim of evaluating the limitations of an asexual reproduction mode in preventing the accumulation of deleterious mutations through natural selection alone. The main considered model is non-classical, as it is a stochasti…
▽ More
We prove the existence and uniqueness of a quasi-stationary distribution for three stochastic processes derived from the model of Muller's ratchet. This model was invented with the aim of evaluating the limitations of an asexual reproduction mode in preventing the accumulation of deleterious mutations through natural selection alone. The main considered model is non-classical, as it is a stochastic diffusion evolving on an irregular set of infinite dimension with hard killing on an hyperplane. We are nonetheless able to prove exponential convergence in total variation to the quasi-stationary distribution even in this case. The parameters in this last convergence result are directly related to the core parameters of Muller's ratchet. The speed of convergence to the quasi-stationary distribution is deduced both for the infinite dimensional model and for approximations with a large yet finite number of potential mutations. Likewise, we give uniform moment estimates of the empirical distribution of mutations in the population under quasi-stationarity.
△ Less
Submitted 30 March, 2024; v1 submitted 29 July, 2020;
originally announced July 2020.
-
A SIR model on a refining spatial grid: Law of Large Numbers
Authors:
M. N'zi,
E. Pardoux,
T. Yeo
Abstract:
We study in this paper a compartmental SIR model for a population distributed in a bounded domain D of $\mathbb{R}^d$, d= 1, 2, or 3. We describe a spatial model for the spread of a disease on a grid of D. We prove two laws of large numbers. On the one hand, we prove that the stochastic model converges to the corresponding deterministic patch model as the size of the population tends to infinity.…
▽ More
We study in this paper a compartmental SIR model for a population distributed in a bounded domain D of $\mathbb{R}^d$, d= 1, 2, or 3. We describe a spatial model for the spread of a disease on a grid of D. We prove two laws of large numbers. On the one hand, we prove that the stochastic model converges to the corresponding deterministic patch model as the size of the population tends to infinity. On the other hand, by letting both the size of the population tend to infinity and the mesh of the grid go to zero, we obtain a law of large numbers in the supremum norm, where the limit is a diffusion SIR model in D.
△ Less
Submitted 22 July, 2020;
originally announced July 2020.
-
A Spatial Stochastic Epidemic Model: Law of Large Numbers and Central Limit Theorem
Authors:
Samuel Bowong Tsakou,
Alphonse Emakoua,
Etienne Pardoux
Abstract:
We consider an individual-based SIR stochastic epidemic model in continuous space. The evolution of the epidemic involves the rates of infection and cure of individuals. We assume that individuals move randomly on the two-dimensional torus according to independent Brownian motions. We define the empirical measures $μ^{S,N}$, $μ^{I,N}$ and $μ^{R,N}$ which describe the evolution of the position of t…
▽ More
We consider an individual-based SIR stochastic epidemic model in continuous space. The evolution of the epidemic involves the rates of infection and cure of individuals. We assume that individuals move randomly on the two-dimensional torus according to independent Brownian motions. We define the empirical measures $μ^{S,N}$, $μ^{I,N}$ and $μ^{R,N}$ which describe the evolution of the position of the susceptible, infected and removed individuals. We prove the convergence in propbability, as $N\rightarrow \infty$, of the sequence $(μ^{S,N},μ^{I,N})$ towards $(μ^{S},μ^{I})$ solution of a system of parabolic PDEs. We show that the sequence $(U^{N}=\sqrt{N}(μ^{S,N}-μ^{S}),V^{N}=\sqrt{N}(μ^{I,N}-μ^{I}))$ converges in law, as $N\rightarrow\infty$, towards a Gaussian distribution valued process, solution of a system of linear PDEs with highly singular Gaussian driving processes. In the case where the individuals do not move we also define and study the law of large numbers and central limit theorem of the sequence $(μ^{S,N},μ^{I,N})$.
△ Less
Submitted 13 July, 2020;
originally announced July 2020.
-
Epidemic models with varying infectivity
Authors:
Raphael Forien,
Guodong Pang,
Etienne Pardoux
Abstract:
We introduce an epidemic model with varying infectivity and general exposed and infectious periods, where the infectivity of each individual is a random function of the elapsed time since infection, those function being i.i.d. for the various individuals in the population. This approach models infection-age dependent infectivity, and extends the classical SIR and SEIR models. We focus on the infec…
▽ More
We introduce an epidemic model with varying infectivity and general exposed and infectious periods, where the infectivity of each individual is a random function of the elapsed time since infection, those function being i.i.d. for the various individuals in the population. This approach models infection-age dependent infectivity, and extends the classical SIR and SEIR models. We focus on the infectivity process (total force of infection at each time), and prove a functional law of large number (FLLN). In the deterministic limit of this LLN, the infectivity process and the susceptible process are determined by a two-dimensional deterministic integral equation. From its solutions, we then derive the exposed, infectious and recovered processes, again using integral equations. For the early phase, we study the stochastic model directly by using an approximate (non--Markovian) branching process, and show that the epidemic grows at an exponential rate on the event of non-extinction, which matches the rate of growth derived from the deterministic linearized equations. We also use these equations to derive the basic reproduction number $R_0$ during the early stage of an epidemic, in terms of the average individual infectivity function and the exponential rate of growth of the epidemic.
△ Less
Submitted 28 May, 2021; v1 submitted 27 June, 2020;
originally announced June 2020.
-
Multi-patch epidemic models with general exposed and infectious periods
Authors:
Guodong Pang,
Etienne Pardoux
Abstract:
We study multi-patch epidemic models where individuals may migrate from one patch to another in either of the susceptible, exposed/latent, infectious and recovered states. We assume that infections occur both locally with a rate that depends on the patch as well as "from distance" from all the other patches. The exposed and infectious periods have general distributions, and are not affected by the…
▽ More
We study multi-patch epidemic models where individuals may migrate from one patch to another in either of the susceptible, exposed/latent, infectious and recovered states. We assume that infections occur both locally with a rate that depends on the patch as well as "from distance" from all the other patches. The exposed and infectious periods have general distributions, and are not affected by the possible migrations of the individuals. The migration processes in either of the three states are assumed to be Markovian, and independent of the exposed and infectious periods. We establish a functional law of large number (FLLN) and a function central limit theorem (FCLT) for the susceptible, exposed/latent, infectious and recovered processes. In the FLLN, the limit is determined by a set of Volterra integral equations. In the special case of deterministic exposed and infectious periods, the limit becomes a system of ODEs with delays. In the FCLT, the limit is given by a set of stochastic Volterra integral equations driven by a sum of independent Brownian motions and continuous Gaussian processes with an explicit covariance structure.
△ Less
Submitted 29 August, 2022; v1 submitted 25 June, 2020;
originally announced June 2020.
-
Functional Limit Theorems for Non-Markovian Epidemic Models
Authors:
Guodong Pang,
Etienne Pardoux
Abstract:
We study non-Markovian stochastic epidemic models (SIS, SIR, SIRS, and SEIR), in which the infectious (and latent/exposing, immune) periods have a general distribution. We provide a representation of the evolution dynamics using the time epochs of infection (and latency/exposure, immunity). Taking the limit as the size of the population tends to infinity, we prove both a functional law of large nu…
▽ More
We study non-Markovian stochastic epidemic models (SIS, SIR, SIRS, and SEIR), in which the infectious (and latent/exposing, immune) periods have a general distribution. We provide a representation of the evolution dynamics using the time epochs of infection (and latency/exposure, immunity). Taking the limit as the size of the population tends to infinity, we prove both a functional law of large number (FLLN) and a functional central limit theorem (FCLT) for the processes of interest in these models. In the FLLN, the limits are a unique solution to a system of deterministic Volterra integral equations, while in the FCLT, the limit processes are multidimensional Gaussian solutions of linear Volterra stochastic integral equations. In the proof of the FCLT, we provide an important Poisson random measures representation of the diffusion-scaled processes converging to Gaussian components driving the limit process.
△ Less
Submitted 23 June, 2021; v1 submitted 6 March, 2020;
originally announced March 2020.
-
Fluctuations around a homogenised semilinear random PDE
Authors:
Martin Hairer,
Étienne Pardoux
Abstract:
We consider a semilinear parabolic partial differential equation in $\mathbf{R}_+\times [0,1]^d$, where $d=1, 2$ or $3$, with a highly oscillating random potential and either homogeneous Dirichlet or Neumann boundary condition. If the amplitude of the oscillations has the right size compared to its typical spatiotemporal scale, then the solution of our equation converges to the solution of a deter…
▽ More
We consider a semilinear parabolic partial differential equation in $\mathbf{R}_+\times [0,1]^d$, where $d=1, 2$ or $3$, with a highly oscillating random potential and either homogeneous Dirichlet or Neumann boundary condition. If the amplitude of the oscillations has the right size compared to its typical spatiotemporal scale, then the solution of our equation converges to the solution of a deterministic homogenised parabolic PDE, which is a form of law of large numbers. Our main interest is in the associated central limit theorem. Namely, we study the limit of a properly rescaled difference between the initial random solution and its LLN limit. In dimension $d=1$, that rescaled difference converges as one might expect to a centred Ornstein-Uhlenbeck process. However, in dimension $d=2$, the limit is a non-centred Gaussian process, while in dimension $d=3$, before taking the CLT limit, we need to subtract at an intermediate scale the solution of a deterministic parabolic PDE, subject (in the case of Neumann boundary condition) to a non-homogeneous Neumann boundary condition. Our proofs make use of the theory of regularity structures, in particular of the very recently developed methodology allowing to treat parabolic PDEs with boundary conditions within that theory.
△ Less
Submitted 5 May, 2021; v1 submitted 7 November, 2019;
originally announced November 2019.
-
Household epidemic models and McKean-Vlasov Poisson driven stochastic differential equations
Authors:
Raphael Forien,
Etienne Pardoux
Abstract:
This paper presents a new view of household epidemic models, where the interaction between the households is of mean field type. We thus obtain in the limit of infinitely many households a nonlinear Markov process solution of a McKean - Vlasov type Poisson driven SDE, and a propagation of chaos result. We also define a basic reproduction number R0, and show that if R0>1, then the nonlinear Markov…
▽ More
This paper presents a new view of household epidemic models, where the interaction between the households is of mean field type. We thus obtain in the limit of infinitely many households a nonlinear Markov process solution of a McKean - Vlasov type Poisson driven SDE, and a propagation of chaos result. We also define a basic reproduction number R0, and show that if R0>1, then the nonlinear Markov process has a unique non trivial ergodic invariant probability measure, whereas if R0<=1, it converges to 0 as t tends to infinity.
△ Less
Submitted 3 October, 2022; v1 submitted 5 July, 2019;
originally announced July 2019.
-
Moderate deviations and extinction of an epidemic
Authors:
Etienne Pardoux
Abstract:
Consider an epidemic model with a constant flux of susceptibles, in a situation where the corresponding deterministic epidemic model has a unique stable endemic equilibrium. For the associated stochastic model, whose law of large numbers limit is the deterministic model, the disease free equilibrium is an absorbing state, which is reached soon or later by the process. However, for a large populati…
▽ More
Consider an epidemic model with a constant flux of susceptibles, in a situation where the corresponding deterministic epidemic model has a unique stable endemic equilibrium. For the associated stochastic model, whose law of large numbers limit is the deterministic model, the disease free equilibrium is an absorbing state, which is reached soon or later by the process. However, for a large population size, i.e. when the stochastic model is close to its deterministic limit, the time needed for the stochastic perturbations to stop the epidemic may be enormous. In this paper, we discuss how the Central Limit Theorem, Moderate and Large Deviations allow us to give estimates of the extinction time of the epidemic, depending upon the size of the population.
△ Less
Submitted 22 May, 2019;
originally announced May 2019.
-
The height process of a continuous state branching process with interaction
Authors:
Zenghu Li,
Etienne Pardoux,
Anton Wakolbinger
Abstract:
For a generalized continuous state branching process with non-vanishing diffusion part, finite expectation and a directed ("left-to-right") interaction, we construct the height process of its forest of genealogical trees. The connection between this height process and the population size process is given by an extension of the second Ray--Knight theorem. This paper generalizes earlier work of the…
▽ More
For a generalized continuous state branching process with non-vanishing diffusion part, finite expectation and a directed ("left-to-right") interaction, we construct the height process of its forest of genealogical trees. The connection between this height process and the population size process is given by an extension of the second Ray--Knight theorem. This paper generalizes earlier work of the two last authors which was restricted to the case of continuous branching mechanisms. Our approach is different from that of Berestycki et al. There the diffusion part of the population process was allowed to vanish, but the class of interactions was more restricted.
△ Less
Submitted 12 November, 2020; v1 submitted 8 April, 2019;
originally announced April 2019.
-
Stochastic epidemics in a homogeneous community
Authors:
Tom Britton,
Etienne Pardoux
Abstract:
These notes describe stochastic epidemics in a homogenous community. Our main concern is stochastic compartmental models (i.e. models where each individual belongs to a compartment, which stands for its status regarding the epidemic under study : S for susceptible, E for exposed, I for infectious, R for recovered) for the spread of an infectious disease. In the present notes we restrict ourselves…
▽ More
These notes describe stochastic epidemics in a homogenous community. Our main concern is stochastic compartmental models (i.e. models where each individual belongs to a compartment, which stands for its status regarding the epidemic under study : S for susceptible, E for exposed, I for infectious, R for recovered) for the spread of an infectious disease. In the present notes we restrict ourselves to homogeneously mixed communities. We present our general model and study the early stage of the epidemic in chapter 1. Chapter 2 studies the particular case of Markov models, especially in the asymptotic of a large population, which leads to a law of large numbers and a central limit theorem. Chapter 3 considers the case of a closed population, and describes the final size of the epidemic (i.e. the total number of individuals who ever get infected). Chapter 4 considers models with a constant influx of susceptibles (either by birth, immigration of loss of immunity of recovered individuals), and exploits the CLT and Large Deviations to study how long it takes for the stochastic disturbances to stop an endemic situation which is stable for the deterministic epidemic model. The document ends with an Appendix which presents several mathematical notions which are used in these notes, as well as solutions to many of the exercises which are proposed in the various chapters.
△ Less
Submitted 5 February, 2020; v1 submitted 16 August, 2018;
originally announced August 2018.
-
Large Deviations of the Exit Measure through a Characteristic Boundary for a Poisson driven SDE
Authors:
Etienne Pardoux,
Brice Samegni-Kepgnou
Abstract:
Let O the basin of attraction of the unique stable equilibrium of a dynamical system, which is the law of large numbers limit of a Poissonian SDE. We consider the law of the exit point from O of that Poissonian SDE. We adapt the approach of M. Day (1990) for the same problem for a Brownian SDE. For that purpose, we will use the Large deviation for the Poissonian SDE reflected at the boundary of O…
▽ More
Let O the basin of attraction of the unique stable equilibrium of a dynamical system, which is the law of large numbers limit of a Poissonian SDE. We consider the law of the exit point from O of that Poissonian SDE. We adapt the approach of M. Day (1990) for the same problem for a Brownian SDE. For that purpose, we will use the Large deviation for the Poissonian SDE reflected at the boundary of O studied in our recent work Pardoux and Samegni (2018).
△ Less
Submitted 15 August, 2018;
originally announced August 2018.
-
Large Deviation Principle for Reflected Poisson driven SDEs in Epidemic Models
Authors:
Etienne Pardoux,
Brice Samegni-Kepgnou
Abstract:
We establish a large deviation principle for a reflected Poisson driven SDE. Our motivation is to study in a forthcoming paper the problem of exit of such a process from the basin of attraction of a locally stable equilibrium associated with its law of large numbers. Two examples are described in which we verify the assumptions that we make to establish the large deviation principle.
We establish a large deviation principle for a reflected Poisson driven SDE. Our motivation is to study in a forthcoming paper the problem of exit of such a process from the basin of attraction of a locally stable equilibrium associated with its law of large numbers. Two examples are described in which we verify the assumptions that we make to establish the large deviation principle.
△ Less
Submitted 14 August, 2018;
originally announced August 2018.
-
Approximation of the Height process of a CSBP with interaction
Authors:
Ibrahima Drame,
Etienne Pardoux
Abstract:
In this work, we first show that the properly rescaled height process of the genealogical tree of a continuous time branching process converges to the height process of the genealogy of a (possibly discontinuous) continuous state branching process. We then prove the same type of result for generalized branching processes with interaction.
In this work, we first show that the properly rescaled height process of the genealogical tree of a continuous time branching process converges to the height process of the genealogy of a (possibly discontinuous) continuous state branching process. We then prove the same type of result for generalized branching processes with interaction.
△ Less
Submitted 1 December, 2019; v1 submitted 19 June, 2017;
originally announced June 2017.
-
Convergence to a Continuous State Branching Process with jumps and Height Process
Authors:
Ibrahima Drame,
Etienne Pardoux
Abstract:
In this work, we study asymptotics of the genealogy of Galton-Watson processes. Thus we consider a offspring distribution such that the rescaled Galton-Watson processes converges to a continuous state branching process (CSBP) with jumps. After we show that the rescaled height (or exploration) process of the corresponding Galton-Watson family tree, converges in a functional sense, to the continuous…
▽ More
In this work, we study asymptotics of the genealogy of Galton-Watson processes. Thus we consider a offspring distribution such that the rescaled Galton-Watson processes converges to a continuous state branching process (CSBP) with jumps. After we show that the rescaled height (or exploration) process of the corresponding Galton-Watson family tree, converges in a functional sense, to the continuous height process that Le Gall and Le Jan introduced in 1998 on their paper "branching processes in Lévy processes : The exploration process".
△ Less
Submitted 18 June, 2017;
originally announced June 2017.
-
Large deviation principle for Poisson driven SDEs in epidemic models
Authors:
Etienne Pardoux,
Brice Samegni-Kepgnou
Abstract:
We consider a general class of epidemic models obtained by applying a random time change to a collection of Poisson processes and we show the large deviation principle for such models. We generalize to a more general situation the approach followed by Dolgoaschinnykh [3] in the case of the SIR epidemic model. Thanks to an assumption which is satisfied in many examples, we simplify the recent work…
▽ More
We consider a general class of epidemic models obtained by applying a random time change to a collection of Poisson processes and we show the large deviation principle for such models. We generalize to a more general situation the approach followed by Dolgoaschinnykh [3] in the case of the SIR epidemic model. Thanks to an assumption which is satisfied in many examples, we simplify the recent work by P. Kratz and E. Pardoux [8].
△ Less
Submitted 6 June, 2016;
originally announced June 2016.
-
Large deviations for infectious diseases models
Authors:
Peter Kratz,
Etienne Pardoux
Abstract:
We establish large deviation estimates for Poisson driven models of infectious disease, and apply those estimates to the time of exit from the basin of attraction of an endemic equilibrium. We apply our results to the classical SIRS model.
We establish large deviation estimates for Poisson driven models of infectious disease, and apply those estimates to the time of exit from the basin of attraction of an endemic equilibrium. We apply our results to the classical SIRS model.
△ Less
Submitted 8 February, 2016;
originally announced February 2016.
-
Non-binary branching process and Non-Markovian exploration process
Authors:
Ibrahima Dramé,
Etienne Pardoux,
Ahmadou Bamba Sow
Abstract:
We study the exploration (or height) process of a continuous time non-binary Galton-Watson random tree, in the subcritical, critical and supercritical cases. Thus we consider the branching process in continuous time (Z_{t})_{t\geq 0}, which describes the number of offspring alive at time t. We then renormalize our branching process and exploration process, and take the weak limit as the size of th…
▽ More
We study the exploration (or height) process of a continuous time non-binary Galton-Watson random tree, in the subcritical, critical and supercritical cases. Thus we consider the branching process in continuous time (Z_{t})_{t\geq 0}, which describes the number of offspring alive at time t. We then renormalize our branching process and exploration process, and take the weak limit as the size of the population tends to infinity. Finally we deduce a Ray-Knight representation.
△ Less
Submitted 5 February, 2016;
originally announced February 2016.
-
On the large time behaviour of the solution of an SDE driven by a Poisson Point Process
Authors:
Elma Nassar,
Etienne Pardoux
Abstract:
We study a stochastic differential equation driven by a Poisson point process, which models continuous changes in a population's environment, as well as the stochastic fixation of beneficial mutations that might compensate for this change. The fixation probability of a given mutation increases as the phenotypic lag X_t between the population and the optimum grows larger, and successful mutations a…
▽ More
We study a stochastic differential equation driven by a Poisson point process, which models continuous changes in a population's environment, as well as the stochastic fixation of beneficial mutations that might compensate for this change. The fixation probability of a given mutation increases as the phenotypic lag X_t between the population and the optimum grows larger, and successful mutations are assumed to fix instantaneously (leading to an adaptive jump). Our main result is that the process is transient (i.e., continued adaptation is impossible) if the rate of environmental change v exceeds a parameter m, which can be interpreted as the rate of adaptation in case every beneficial mutation gets fixed with probability 1. If v < m, the process is positive recurrent, while in the limiting case m=v, null recurrence or transience depends upon additional technical conditions. We show how our results can be extended to the case of a time varying rate of environmental change.
△ Less
Submitted 3 February, 2016;
originally announced February 2016.
-
Continuity of the Feynman-Kac formula for a generalized parabolic equation
Authors:
Etienne Pardoux,
Aurel Rascanu
Abstract:
It is well-known since the work of Pardoux and Peng [12] that Backward Stochastic Differential Equations provide probabilistic formulae for the solution of (systems of) second order elliptic and parabolic equations, thus providing an extension of the Feynman-Kac formula to semilinear PDEs, see also Pardoux and Rascanu [14]. This method was applied to the class of PDEs with a nonlinear Neumann boun…
▽ More
It is well-known since the work of Pardoux and Peng [12] that Backward Stochastic Differential Equations provide probabilistic formulae for the solution of (systems of) second order elliptic and parabolic equations, thus providing an extension of the Feynman-Kac formula to semilinear PDEs, see also Pardoux and Rascanu [14]. This method was applied to the class of PDEs with a nonlinear Neumann boundary condition first by Pardoux and Zhang [15]. However, the proof of continuity of the extended Feynman-Kac formula with respect to x (resp. to (t,x)) is not correct in that paper. Here we consider a more general situation, where both the equation and the boundary condition involve the (possibly multivalued) gradient of a convex function. We prove the required continuity. The result for the class of equations studied in [15] is a Corollary of our main results.
△ Less
Submitted 3 February, 2016;
originally announced February 2016.
-
Averaging for SDE-BSDE with null recurrent fast component Application to homogenization in a non periodic media
Authors:
K Bahlali,
A Elouaflin,
E Pardoux
Abstract:
We establish an averaging principle for a family of solutions$(X^{\varepsilon}, Y^{\varepsilon})$ $ :=$ $(X^{1,\,\varepsilon},\,X^{2,\,\varepsilon},\, Y^{\varepsilon})$ of a system of SDE-BSDEwith a null recurrent fast component $X^{1,\,\varepsilon}$. Incontrast to the classical periodic case, we can not rely on aninvariant probability and the slow forward component$X^{2,\,\varepsilon}$ cannot be…
▽ More
We establish an averaging principle for a family of solutions$(X^{\varepsilon}, Y^{\varepsilon})$ $ :=$ $(X^{1,\,\varepsilon},\,X^{2,\,\varepsilon},\, Y^{\varepsilon})$ of a system of SDE-BSDEwith a null recurrent fast component $X^{1,\,\varepsilon}$. Incontrast to the classical periodic case, we can not rely on aninvariant probability and the slow forward component$X^{2,\,\varepsilon}$ cannot be approximated by a diffusion process.On the other hand, we assume that the coefficients admit a limit in aCesaro sense. In such a case, the limit coefficients may havediscontinuity. We show that we can approximate the triplet$(X^{1,\,\varepsilon},\, X^{2,\,\varepsilon},\, Y^{\varepsilon})$ bya system of SDE-BSDE $(X^1, X^2, Y)$ where $X := (X^1, X^2)$ is aMarkov diffusion which is the unique (in law) weak solution of theaveraged forward component and $Y$ is the unique solution to the averaged backward component. This is done with a backward component whosegenerator depends on the variable $z$. Asapplication, we establish an homogenization result for semilinearPDEs when the coefficients can be neither periodic nor ergodic. Weshow that the averaged BDSE is related to the averaged PDE via aprobabilistic representation of the (unique) Sobolev $W\_{d+1,\text{loc}}^{1,2}(\R\_+\times\R^d)$--solution of the limitPDEs. Our approach combines PDE methods and probabilistic argumentswhich are based on stability property and weak convergence of BSDEsin the S-topology.
△ Less
Submitted 31 August, 2015;
originally announced August 2015.
-
A mixing tree-valued process arising under neutral evolution with recombination
Authors:
Andrej Depperschmidt,
Etienne Pardoux,
Peter Pfaffelhuber
Abstract:
The genealogy at a single locus of a constant size $N$ population in equilibrium is given by the well-known Kingman's coalescent. When considering multiple loci under recombination, the ancestral recombination graph encodes the genealogies at all loci in one graph. For a continuous genome $G$, we study the tree-valued process $(T^N_u)_{u\in G}$ of genealogies along the genome in the limit…
▽ More
The genealogy at a single locus of a constant size $N$ population in equilibrium is given by the well-known Kingman's coalescent. When considering multiple loci under recombination, the ancestral recombination graph encodes the genealogies at all loci in one graph. For a continuous genome $G$, we study the tree-valued process $(T^N_u)_{u\in G}$ of genealogies along the genome in the limit $N\to\infty$. Encoding trees as metric measure spaces, we show convergence to a tree-valued process with cadlag paths. In addition, we study mixing properties of the resulting process for loci which are far apart.
△ Less
Submitted 5 May, 2015;
originally announced May 2015.
-
Numerical methods in the context of compartmental models in epidemiology
Authors:
Peter Kratz,
Etienne Pardoux,
Brice Samegni Kepgnou
Abstract:
We consider compartmental models in epidemiology. For the study of the divergence of the stochastic model from its corresponding deterministic limit (i.e., the solution of an ODE) for long time horizon, a large deviations principle suggests a thorough numerical analysis of the two models. The aim of this paper is to present three such motivated numerical works. We first compute the solution of the…
▽ More
We consider compartmental models in epidemiology. For the study of the divergence of the stochastic model from its corresponding deterministic limit (i.e., the solution of an ODE) for long time horizon, a large deviations principle suggests a thorough numerical analysis of the two models. The aim of this paper is to present three such motivated numerical works. We first compute the solution of the ODE model by means of a non-standard finite difference scheme. Next we solve a constraint optimization problem via discrete-time dynamic programming: this enables us to compute the leading term in the large deviations principle of the time of extinction of a given disease. Finally, we apply the τ-leaping algorithm to the stochastic model in order to simulate its solution efficiently. We illustrate these numerical methods by applying them to two examples.
△ Less
Submitted 15 January, 2015;
originally announced January 2015.
-
Height and the total mass of the forest of genealogical trees of a large population with general competition
Authors:
Vi Le,
Etienne Pardoux
Abstract:
We consider branching processes with interaction in continuous time, both with values in the integers and in the reals (in the second case we restrict ourselves to continuous processes), which model the evolution of the size of a population. We assume that for large population size the interaction is of the type of a competition, which limits the size of the population. We discuss in which cases t…
▽ More
We consider branching processes with interaction in continuous time, both with values in the integers and in the reals (in the second case we restrict ourselves to continuous processes), which model the evolution of the size of a population. We assume that for large population size the interaction is of the type of a competition, which limits the size of the population. We discuss in which cases the interaction is strong enough so that the extinction time (or equivalently the height of the forest of genealogical trees) remains finite, as the number of ancestors tends to infinity, or even such that the length of the forest of genealogical trees (which in the case of continuous state is rather called its total mass) remains finite, as the ancestral population size tends to infinity.
△ Less
Submitted 21 October, 2014;
originally announced October 2014.
-
A Wong-Zakai theorem for stochastic PDEs
Authors:
Martin Hairer,
Étienne Pardoux
Abstract:
We prove a version of the Wong-Zakai theorem for one-dimensional parabolic nonlinear stochastic PDEs driven by space-time white noise. As a corollary, we obtain a detailed local description of solutions.
Dedicated to the memory of Kiyosi Itô on the occasion of the 100th anniversary of his birth.
We prove a version of the Wong-Zakai theorem for one-dimensional parabolic nonlinear stochastic PDEs driven by space-time white noise. As a corollary, we obtain a detailed local description of solutions.
Dedicated to the memory of Kiyosi Itô on the occasion of the 100th anniversary of his birth.
△ Less
Submitted 20 April, 2015; v1 submitted 10 September, 2014;
originally announced September 2014.
-
A path-valued Markov process indexed by the ancestral mass
Authors:
Etienne Pardoux,
Anton Wakolbinger
Abstract:
A family of Feller branching diffusions $Z^x$, $x \ge 0$, with nonlinear drift and initial value $x$ can, with a suitable coupling over the {\em ancestral masses} $x$, be viewed as a path-valued process indexed by $x$. For a coupling due to Dawson and Li, which in case of a linear drift describes the corresponding Feller branching diffusion, and in our case makes the path-valued process Markovian,…
▽ More
A family of Feller branching diffusions $Z^x$, $x \ge 0$, with nonlinear drift and initial value $x$ can, with a suitable coupling over the {\em ancestral masses} $x$, be viewed as a path-valued process indexed by $x$. For a coupling due to Dawson and Li, which in case of a linear drift describes the corresponding Feller branching diffusion, and in our case makes the path-valued process Markovian, we find an SDE solved by $Z$, which is driven by a random point measure on excursion space. In this way we are able to identify the infinitesimal generator of the path-valued process. We also establish path properties of $x\mapsto Z^x$ using various couplings of $Z$ with classical Feller branching diffusions.
△ Less
Submitted 20 February, 2015; v1 submitted 5 September, 2014;
originally announced September 2014.
-
Stochastic Variational Inequalities on Non-Convex Domains
Authors:
Rainer Buckdahn,
Lucian Maticiuc,
Etienne Pardoux,
Aurel Răşcanu
Abstract:
The objective of this work is to prove, in a first step, the existence and the uniqueness of a solution of the following multivalued deterministic differential equation: $dx(t)+\partial ^-\varphi (x(t))(dt)\ni dm(t),\ t>0$, $x(0)=x_0$, where $m:\mathbb{R}_+\rightarrow\mathbb{R}^d$ is a continuous function and $\partial^-\varphi$ is the Fréchet subdifferential of a semiconvex function $\varphi$; th…
▽ More
The objective of this work is to prove, in a first step, the existence and the uniqueness of a solution of the following multivalued deterministic differential equation: $dx(t)+\partial ^-\varphi (x(t))(dt)\ni dm(t),\ t>0$, $x(0)=x_0$, where $m:\mathbb{R}_+\rightarrow\mathbb{R}^d$ is a continuous function and $\partial^-\varphi$ is the Fréchet subdifferential of a semiconvex function $\varphi$; the domain of $\varphi$ can be non-convex, but some regularities of the boundary are required.
The continuity of the map $m\mapsto x:C([0,T];\mathbb{R}^{d})\rightarrow C([0,T] ;\mathbb{R}^{d})$, which associate the input function $m$ with the solution $x$ of the above equation, as well as tightness criteria allow to pass from the above deterministic case to the following stochastic variational inequality driven by a multi-dimensional Brownian motion: $X_t+K_t = ξ+\int_0^t F(s,X_{s})ds + \int_0^t G(s,X_s) dB_s,\; t\geq0$, $\;$ with $dK_{t}(ω)\in\partial^-\varphi( X_t (ω))(dt)$.
△ Less
Submitted 9 May, 2015; v1 submitted 7 July, 2014;
originally announced July 2014.
-
Invariant measure selection by noise: An Example
Authors:
Jonathan C. Mattingly,
Etienne Pardoux
Abstract:
We consider a deterministic system with two conserved quantities and infinity many invariant measures. However the systems possess a unique invariant measure when enough stochastic forcing and balancing dissipation are added. We then show that as the forcing and dissipation are removed a unique limit of the deterministic system is selected. The exact structure of the limiting measure depends on th…
▽ More
We consider a deterministic system with two conserved quantities and infinity many invariant measures. However the systems possess a unique invariant measure when enough stochastic forcing and balancing dissipation are added. We then show that as the forcing and dissipation are removed a unique limit of the deterministic system is selected. The exact structure of the limiting measure depends on the specifics of the stochastic forcing.
△ Less
Submitted 14 March, 2014;
originally announced March 2014.
-
Homogenization of a periodic semilinear elliptic degenerate PDE
Authors:
Etienne Pardoux,
Ahmadou Bamba Sow
Abstract:
In this paper a semilinear elliptic PDE with rapidly oscillating coefficients is homogenized. The novetly of our result lies in the fact that we allow the second order part of the differential operator to be degenerate in some portion of $\R^d$.\\ Our fully probabilistic method is based on the connection between PDEs and BSDEs with random terminal time and the weak convergence of a class of diffus…
▽ More
In this paper a semilinear elliptic PDE with rapidly oscillating coefficients is homogenized. The novetly of our result lies in the fact that we allow the second order part of the differential operator to be degenerate in some portion of $\R^d$.\\ Our fully probabilistic method is based on the connection between PDEs and BSDEs with random terminal time and the weak convergence of a class of diffusion processes.
△ Less
Submitted 6 May, 2013;
originally announced May 2013.
-
Binary trees, exploration processes, and an extented Ray--Knight Theorem
Authors:
Mamadou Ba,
Etienne Pardoux,
Ahmadou Bamba Sow
Abstract:
We study the bijection between binary Galton--Watson trees in continuous time and their exploration process, both in the sub- and in the supercritical cases. We then take the limit over renormalized quantities, as the size of the population tends to infinity. We thus deduce Delmas' generalization of the second Ray--Knight theorem.
We study the bijection between binary Galton--Watson trees in continuous time and their exploration process, both in the sub- and in the supercritical cases. We then take the limit over renormalized quantities, as the size of the population tends to infinity. We thus deduce Delmas' generalization of the second Ray--Knight theorem.
△ Less
Submitted 6 May, 2013;
originally announced May 2013.
-
From Brownian motion with a local time drift to Feller's branching diffusion with logistic growth
Authors:
Etienne Pardoux,
Anton Wakolbinger
Abstract:
We give a new proof for a Ray-Knight representation of Feller's branching diffusion with logistic growth in terms of the local times of a reflected Brownian motion $H$ with a drift that is affine linear in the local time accumulated by $H$ at its current level. In \cite{LPW}, such a representation was obtained by an approximation through Harris paths that code the genealogies of particle systems.…
▽ More
We give a new proof for a Ray-Knight representation of Feller's branching diffusion with logistic growth in terms of the local times of a reflected Brownian motion $H$ with a drift that is affine linear in the local time accumulated by $H$ at its current level. In \cite{LPW}, such a representation was obtained by an approximation through Harris paths that code the genealogies of particle systems. The present proof is purely in terms of stochastic analysis, and is inspired by previous work of Norris, Rogers and Williams \cite{NRW}.
△ Less
Submitted 6 May, 2013;
originally announced May 2013.
-
Trees under attack: a Ray-Knight representation of Feller's branching diffusion with logistic growth
Authors:
Vi Le,
Etienne Pardoux,
Anton Wakolbinger
Abstract:
We obtain a representation of Feller's branching diffusion with logistic growth in terms of the local times of a reflected Brownian motion $H$ with a drift that is affine linear in the local time accumulated by $H$ at its current level. As in the classical Ray-Knight representation, the excursions of $H$ are the exploration paths of the trees of descendants of the ancestors at time $t=0$, and the…
▽ More
We obtain a representation of Feller's branching diffusion with logistic growth in terms of the local times of a reflected Brownian motion $H$ with a drift that is affine linear in the local time accumulated by $H$ at its current level. As in the classical Ray-Knight representation, the excursions of $H$ are the exploration paths of the trees of descendants of the ancestors at time $t=0$, and the local time of $H$ at height $t$ measures the population size at time $t$ (see e.g. \cite{LG4}). We cope with the dependence in the reproduction by introducing a pecking order of individuals: an individual explored at time $s$ and living at time $t=H_s$ is prone to be killed by any of its contemporaneans that have been explored so far. The proof of our main result relies on approximating $H$ with a sequence of Harris paths $H^N$ which figure in a Ray-Knight representation of the total mass of a branching particle system. We obtain a suitable joint convergence of $H^N$ together with its local times {\em and} with the Girsanov densities that introduce the dependence in the reproduction.
△ Less
Submitted 6 May, 2013;
originally announced May 2013.
-
Branching processes with competition and generalized Ray Knight Theorem
Authors:
Mamadou Ba,
Etienne Pardoux
Abstract:
We consider a discrete model of population with interaction where the birth and death rates are non linear functions of the population size. After proceeding to renormalization of the model parameters, we obtain in the limit of large population that the population size evolves as a diffusion solution of the SDE Z^x_t =x+\int_0^t f(Z^x_s)ds+2\int_0^t\int_0^{Z^x_s}W(ds,du), where W(ds,du) is a time…
▽ More
We consider a discrete model of population with interaction where the birth and death rates are non linear functions of the population size. After proceeding to renormalization of the model parameters, we obtain in the limit of large population that the population size evolves as a diffusion solution of the SDE Z^x_t =x+\int_0^t f(Z^x_s)ds+2\int_0^t\int_0^{Z^x_s}W(ds,du), where W(ds,du) is a time space white noise on ([0,\infty))^2. We give a Ray-Knight representation of this diffusion in terms of the local times of a reflected Brownian motion H with a drift that depends upon the local time accumulated by H at its current level, through the function f'/2.
△ Less
Submitted 2 April, 2013;
originally announced April 2013.
-
Muller's ratchet clicks in finite time
Authors:
Julien Audiffren,
Etienne Pardoux
Abstract:
We consider the accumulation of deleterious mutations in an asexual population, a phenomenon known as Muller's ratchet, using the continuous time model proposed in Etheridge et all. \cite{epw}. We show that for any parameter $λ>0$ (the rate at which mutations occur), for any $α>0$ (the toxicity of the mutations) and for any size $N>0$ of the population, the ratchet clicks a.s. in finite time. That…
▽ More
We consider the accumulation of deleterious mutations in an asexual population, a phenomenon known as Muller's ratchet, using the continuous time model proposed in Etheridge et all. \cite{epw}. We show that for any parameter $λ>0$ (the rate at which mutations occur), for any $α>0$ (the toxicity of the mutations) and for any size $N>0$ of the population, the ratchet clicks a.s. in finite time. That is to say the minimum number of deleterious mutations in the population goes to infinity a.s.
△ Less
Submitted 8 March, 2013;
originally announced March 2013.
-
Random homogenisation of a highly oscillatory singular potential
Authors:
Martin Hairer,
Etienne Pardoux,
Andrey Piatnitski
Abstract:
In this article, we consider the problem of homogenising the linear heat equation perturbed by a rapidly oscillating random potential. We consider the situation where the space-time scaling of the potential's oscillations is \textit{not} given by the diffusion scaling that leaves the heat equation invariant. Instead, we treat the case where spatial oscillations are much faster than temporal oscill…
▽ More
In this article, we consider the problem of homogenising the linear heat equation perturbed by a rapidly oscillating random potential. We consider the situation where the space-time scaling of the potential's oscillations is \textit{not} given by the diffusion scaling that leaves the heat equation invariant. Instead, we treat the case where spatial oscillations are much faster than temporal oscillations. Under suitable scaling of the amplitude of the potential, we prove convergence to a deterministic heat equation with constant potential, thus completing the results previously obtained in \cite{MR2962093}.
△ Less
Submitted 8 March, 2013;
originally announced March 2013.
-
Lambda-lookdown model with selection
Authors:
Boubacar Bah,
Etienne Pardoux
Abstract:
The goal of this paper is to study the lookdown model with selection in the case of a population containing two types of individuals, with a reproduction model which is dual to the $Λ$-coalescent. In particular we formulate the infinite population "$Λ$-lookdown model with selection". When the measure $Λ$ gives no mass to 0, we show that the proportion of one of the two types converges, as the popu…
▽ More
The goal of this paper is to study the lookdown model with selection in the case of a population containing two types of individuals, with a reproduction model which is dual to the $Λ$-coalescent. In particular we formulate the infinite population "$Λ$-lookdown model with selection". When the measure $Λ$ gives no mass to 0, we show that the proportion of one of the two types converges, as the population size $N$ tends to infinity, towards the solution of a stochastic differential equation driven by a Poisson point process. We show that one of the two types fixates in finite time if and only if the $Λ$-coalescent comes down from infinity. We give precise asymptotic results in the case of the Bolthausen-Sznitman coalescent. We also consider the general case of a combination of the Kingman and the $Λ$-lookdown model.
△ Less
Submitted 21 October, 2014; v1 submitted 8 March, 2013;
originally announced March 2013.
-
Homogenization of a singular random one dimensional parabolic PDE with time varying coefficients
Authors:
E. Pardoux,
A. Piatnitski
Abstract:
The paper studies homogenization problem for a non-autonomous parabolic equation with a large random rapidly oscillating potential in the case of one dimensional spatial variable. We show that if the potential is a statistically homogeneous rapidly oscillating function of both temporal and spatial variables then, under proper mixing assumptions, the limit equation is deterministic and the conver…
▽ More
The paper studies homogenization problem for a non-autonomous parabolic equation with a large random rapidly oscillating potential in the case of one dimensional spatial variable. We show that if the potential is a statistically homogeneous rapidly oscillating function of both temporal and spatial variables then, under proper mixing assumptions, the limit equation is deterministic and the convergence in probability holds. To the contrary, for the potential having a microstructure only in one of these variables, the limit problem is stochastic and we only prove the convergence in law.
△ Less
Submitted 29 December, 2009;
originally announced December 2009.
-
Diffusion limit for many particles in a periodic stochastic acceleration field
Authors:
Yves Elskens,
Etienne Pardoux
Abstract:
The one-dimensional motion of any number $\cN$ of particles in the field of many independent waves (with strong spatial correlation) is formulated as a second-order system of stochastic differential equations, driven by two Wiener processes. In the limit of vanishing particle mass ${\mathfrak{m}} \to 0$, or equivalently of large noise intensity, we show that the momenta of all $N$ particles conv…
▽ More
The one-dimensional motion of any number $\cN$ of particles in the field of many independent waves (with strong spatial correlation) is formulated as a second-order system of stochastic differential equations, driven by two Wiener processes. In the limit of vanishing particle mass ${\mathfrak{m}} \to 0$, or equivalently of large noise intensity, we show that the momenta of all $N$ particles converge weakly to $N$ independent Brownian motions, and this convergence holds even if the noise is periodic. This justifies the usual application of the diffusion equation to a family of particles in a unique stochastic force field. The proof rests on the ergodic properties of the relative velocity of two particles in the scaling limit.
△ Less
Submitted 5 November, 2008;
originally announced November 2008.
-
Viscosity solutions for systems of parabolic variational inequalities
Authors:
Lucian Maticiuc,
Etienne Pardoux,
Aurel Răşcanu,
Adrian Zălinescu
Abstract:
In this paper, we first define the notion of viscosity solution for the following system of partial differential equations involving a subdifferential operator:\[\{[c]{l}\dfrac{\partial u}{\partial t}(t,x)+\mathcal{L}_tu(t,x)+f(t,x,u(t,x))\in\partialφ(u(t,x)),\quad t\in[0,T),x\in\mathbb{R}^d, u(T,x)=h(x),\quad x\in\mathbb{R}^d,\] where $\partialφ$ is the subdifferential operator of the proper co…
▽ More
In this paper, we first define the notion of viscosity solution for the following system of partial differential equations involving a subdifferential operator:\[\{[c]{l}\dfrac{\partial u}{\partial t}(t,x)+\mathcal{L}_tu(t,x)+f(t,x,u(t,x))\in\partialφ(u(t,x)),\quad t\in[0,T),x\in\mathbb{R}^d, u(T,x)=h(x),\quad x\in\mathbb{R}^d,\] where $\partialφ$ is the subdifferential operator of the proper convex lower semicontinuous function $φ:\mathbb{R}^k\to (-\infty,+\infty]$ and $\mathcal{L}_t$ is a second differential operator given by $\mathcal{L}_tv_i(x)={1/2}\operatorname {Tr}[σ(t,x)σ^*(t,x)\mathrm{D}^2v_i(x)]+< b(t,x),\nabla v_i(x)>$, $i\in\bar{1,k}$. We prove the uniqueness of the viscosity solution and then, via a stochastic approach, prove the existence of a viscosity solution $u:[0,T]\times\mathbb{R}^d\to\mathbb{R}^k$ of the above parabolic variational inequality.
△ Less
Submitted 23 February, 2010; v1 submitted 28 July, 2008;
originally announced July 2008.