Thierry E Huillet

A BRANCHING DIFFUSION MODEL OF SELECTION: FROM THE NEUTRAL WRIGHT-FISHER CASE TO THE ONE INCLUDING MUTATIONS

Keywords: Branching diffusion model of selection Diffusions, Doob transform, killing, branching, quasi-stationarity, Wright-Fisher model, neutral or with mutation, selection PACS classification: 87.23.Cc, 02.50.Ey, 87.23

We consider diffusion processes xt on the unit interval. Doobtransformation techniques consist of a selection of xt-paths procedure. The law of the transformed process is the one of a branching diffusion system of particles, each diffusing like a new process xt, superposing an additional drift to the one of xt. Killing and/or branching of xt-particles occur at some space-dependent rate λ. For this transformed process, so in the class of branching diffusions, the question arises as to whether the particle system is sub-critical, critical or super-critical. In the first two cases, extinction occurs with probability one.

We apply this circle of ideas to diffusion processes arising in population genetics. In this setup, the process xt is a Wright-Fisher (WF) diffusion, either neutral or with mutations.

We study a particular Doob transform which is based on the exponential function in the usual fitness parameter σ. We have in mind that this is an alternative way to introduce selection or fitness in both WF-like diffusions, leading to branching diffusion models ideas. For this Doob-transform model of fitness, the usual selection drift σx (1 -x) should be superposed to the one of xt to form xt which is the process that can branch, binarily.

In the first neutral case, there is a trade-off between branching events giving birth to new particles and absorption at the boundaries, killing the particles. Under our assumptions, the branching diffusion process gets eventually globally extinct in finite time with exponential tails.

In the second case with mutations, there is a trade-off between killing events removing some particles from the system and reflection at the boundaries where the particles survive. This branching diffusion process also gets eventually globally extinct but in very long finite time with power-law tails.

Our approach relies on the spectral expansion of the transition probability kernels of both xt and xt.

Introduction

We consider diffusion processes on the unit interval with in mind a series of elementary stochastic models arising chiefly in population dynamics. Special emphasis is put on Doob-transformation techniques of the diffusion processes under concern.

Most of the manuscript's content focuses on the specific Wright-Fisher (WF) diffusion model and some of its variations, describing the evolution of one two-locus colony undergoing random mating, possibly under the additional actions of mutation and selection. These models found their way over the last sixty years, chiefly in mathematical population genetics. We refer to the general monographs [START_REF] Crow | An introduction to population genetics theory[END_REF], [START_REF] Maruyama | Stochastic problems in population genetics[END_REF], [START_REF] Ewens | Mathematical population genetics. I. Theoretical introduction[END_REF], [START_REF] Durrett | Probability models for DNA sequence evolution[END_REF] and [START_REF] Gillespie | The Causes of Molecular Evolution[END_REF].

We now describe the content of this work in some more details.

Section 2 is devoted to generalities on one-dimensional diffusions on the unit interval [0, 1], say (x t ; t ≥ 0). Special emphasis is put on the Kolmogorov backward and forward equations, while stressing the crucial role played by the boundaries in such one-dimensional diffusion problems. Some questions such as the meaning of speed and scale functions, existence of an invariant measure, random time change... are addressed in the light of the Feller classification of boundaries. When the boundaries are absorbing, the important problem of evaluating additive functionals α along sample paths is then briefly discussed, emphasizing the prominent role played by the Green function of the model.

So far, we have dealt with a given diffusion process x t , and recalled the various ingredients for computing the expectations of various quantities of interest, summing up over the history of its paths. In this setup, there is no distinction among paths with different destinations, nor did we allow for annihilation or creation of paths inside the domain. The Doob transform of paths is an invitation to do so. This important class of transformations is a particular instance of a more general construction based on multiplicative functionals. We fix the background.

Roughly speaking, the Doob transformation of paths procedure allows to select sample paths x → y within any laps of time t, favoring large values of the ratio α (y) /α (x) for some specific functional α > 0 that fixes the selection problem under study. The process solving this selection of paths procedure belongs to a class of branching diffusion processes, where independent particles diffusing like a new process x t inside the interval are allowed to duplicate would the visited region of the state-space fulfill the selection of paths criterion or to die, if not. In the process, advantageous regions of the state-space are reinforced while unfavorable ones are left unexplored which is a reasonable physical way to look at selection of paths. The new process x t alluded to is obtained from x t just after superposing an additional suitable drift to the latter process. An important parameter is the statedependent rate λ at which killing and/or branching occur. Depending on λ and on the type of boundaries which {0, 1} are to x t , the full transformed process can have two stopping times: the time to absorption at the boundaries and the killing time inside the domain. Besides, there is or not an opportunity that the x t -particles duplicate, leading or not to a daughters particle system evolving independently starting from where their mother particle died. The killing/branching issues depend on the sign of λ.

It turns out that the same diffusion methods used in the previous discussion on simple diffusions apply to the transformed processes obtained after the induced change of measure. We develop this circle of ideas.

We next apply these general ideas to diffusion processes arising in population genetics.

In Section 3 we start recalling that WF diffusion models with various drifts are continuous space-time models which can be obtained as scaling limits of a biased discrete Galton-Watson model with a conservative number of offsprings over the generations. Sections 4 and 5 are devoted to a detailed study of both the neutral WF diffusion process (WFN), the WF diffusion with selection (WFS), the WF diffusion with mutations (WFM) and the WF diffusion with mutations and selection (WFMS) respectively.

In this context, our suggestion is the following one: we can view the action of selection (or fitness) on the evolution of the allele frequency distribution, either neutral or with mutations, as a functional deformation of the sample paths of the original process, say x t , favoring initial values x 0 = x with small α (x) and terminal values x t = y with large α (y) , for each t. In our construction, α (x) = e σx , σ > 0, is the chosen exponential fitness functional. Stated differently and more precisely, if p (x; t, y) is the transition probability density of x t (either WFN or WFM), our model of the action of fitness is p (x; t, y) fitness → p (x; t, y) = e σy e σx p (x; t, y) .

With this choice of α, the modification consists of selecting those paths x → y of x t for which e σ(y-x) is large. As a result of this transformation of paths, the usual positive selection drift σx (1 -x) has to be superposed to the one of x t to form the new process x t , but our functional definition of fitness also generates an additional branching multiplicative term, translating that a particle system pops in: The resulting transformed process is a (binary) branching diffusion of WF diffusions x t . We may call the obtained processes the branching neutral Wright-Fisher process and the branching Wright-Fisher process with mutations. This point of view seems to be new, to the best of the author's knowledge.

Because the spectral representation of both transition probability densities of WFN or WFM are known explicitly from the works of Crow and Kimura (see [START_REF] Kimura | Diffusion Models in Population Genetics[END_REF], [START_REF] Crow | An introduction to population genetics theory[END_REF] and [START_REF] Crow | Some genetic problems in natural populations[END_REF]), some easy consequences on the spectral structures of the branching transformed processes are available. For instance, it is possible to decide whether the BD process is sub-critical, critical or super-critical in the sense of ([1] and [START_REF] Asmussen | Some modified branching diffusion models[END_REF]).

In Section 6, we therefore give a detailed study of the binary branching diffusion process obtained while using the Doob transform α (x) = e σx when the starting point process is a WFN diffusion process. We end up with a branching particle system, each diffusing according to the WF model with a selection drift, but branching at a bounded rate b > 0. In this setup, the particles cannot get killed, rather they are allowed either to survive or to split: the transformed process is a pure binary branching diffusion. For this super-critical binary branching diffusion process, there is a trade-off between branching events giving birth to new particles and absorption at the boundaries, killing the particles. Thanks to the spectral representation of the WFN process, this problem is amenable to the results obtained in ([1] and [START_REF] Asmussen | Some modified branching diffusion models[END_REF]). Under our assumptions, the branching diffusion process turns out to be globally sub-critical: the branching diffusion process gets eventually globally extinct in finite exponential time. This requires the computation of the ground states associated with the smallest nonnegative eigenvalue of the infinitesimal generator of the transformed process which are here shown to be explicit. In particular, the expression of the quasi-stationary distribution of the particle system can be obtained in closed-form.

In Section 7, we study the binary branching diffusion process obtained while using the same Doob transform, when the starting point process is now a WF diffusion process with mutations, assuming reflecting boundaries. We end up in a branching particle system, each diffusing according to the WF model with a mutation and selection drift, but branching at quadratic rate λ, which is bounded from below and above. Although the particles are still allowed to split, they can now also get killed at the branching times: the transformed process is again a binary branching diffusion but with killing now allowed. In this setup, there is a competition between branching/ killing events and reflection at the boundaries where the particles survive. This problem is also amenable to the results obtained in ( [START_REF] Asmussen | Strong limit theorems for general supercritical branching processes with applications to branching diffusions[END_REF] and [START_REF] Asmussen | Some modified branching diffusion models[END_REF]) and we end up now in a globally critical branching particle system, each diffusing according to the WF model with a mutation and selection drift. This branching diffusion process turns out to be globally critical: it also gets eventually globally extinct but now in long finite time, with power-law tails.

Diffusion processes on the unit interval and Doob transforms

We start with generalities on one-dimensional diffusions with the WF model and its relatives in mind. For more technical details, we refer to [START_REF] Dynkin | Markov processes. Vols. I, II. Translated with the authorization and assistance of the author[END_REF], [START_REF] Ethier | Markov processes. Characterization and convergence[END_REF], [START_REF] Karlin | A second course in stochastic processes[END_REF] and [START_REF] Mandl | Analytical treatment of one-dimensional Markov processes[END_REF].

We also introduce Doob transforms as particular instances of the modification of the original diffusion process through a multiplicative functional.

2.1. One-dimensional diffusions on the interval [0, 1]. Let (w t ; t ≥ 0) be a standard one-dimensional Brownian (Wiener) motion. We consider a 1-dimensional Itô diffusion driven by (w t ; t ≥ 0) on the interval say [0, 1], see [START_REF] Itô | On stochastic differential equations[END_REF]. We assume it has locally Lipschitz continuous drift f (x) and local standard deviation (volatility) g (x), namely we consider the stochastic differential equation (SDE):

(1)

dx t = f (x t ) dt + g (x t ) dw t , x 0 = x ∈ I := (0, 1) .
The condition on f (x) and g (x) guarantees in particular that there is no point x * in I for which |f (x)| or |g (x)| would blow up and diverge as |x -x * | → 0.

The Kolmogorov backward infinitesimal generator of ( 1) is

G = f (x) ∂ x + 1 2 g 2 (x) ∂ 2
x . As a result, for all suitable ψ in the domain of the operator S t := e tG , u := u (x, t) = Eψ (x t∧τ x ) satisfies the Kolmogorov backward equation (KBE)

∂ t u = G (u) ; u (x, 0) = ψ (x) .
In the definition of the mathematical expectation u, we have t ∧ τ x := inf (t, τ x ) where τ x indicates a random time at which the process should possibly be stopped (absorbed), given the process was started in x. The description of this (adapted) absorption time is governed by the type of boundaries which ∂I := {0, 1} are to (x t ; t ≥ 0) . A classification of the boundaries exists, due to Feller (see [START_REF] Karlin | A second course in stochastic processes[END_REF] pp. 226): they can be either accessible (namely exit or regular), or inaccessible (namely entrance or natural).

2.2.

Natural coordinate, scale, speed measure, time change. For such Markovian diffusions, it is interesting to consider the G-harmonic coordinate ϕ ∈ C 2 belonging to the kernel of G, i.e. satisfying G (ϕ) = 0. For ϕ and its derivative ϕ ′ := dϕ/dy, with (x 0 , y 0 ) ∈ (0, 1), one finds

ϕ ′ (y) = ϕ ′ (y 0 ) e -2 y y 0 f (z) g 2 (z) dz ϕ (x) = ϕ (x 0 ) + ϕ ′ (y 0 ) x x0 e -2 y y 0 f (z) g 2 (z) dz dy.
One should choose a version of ϕ satisfying ϕ ′ (y) > 0, y ∈ I. The function ϕ kills the drift f of (x t ; t ≥ 0) in the sense that, considering the change of variable

y t = ϕ (x t ) , dy t = (ϕ ′ g) ϕ -1 (y t ) dw t , y 0 = ϕ (x)
. The drift-less diffusion (y t ; t ≥ 0) is often termed the diffusion in natural coordinates with state-space [ϕ (0) , ϕ (1)] =: ϕ (I). Its volatility is g (y) := (ϕ ′ g) ϕ -1 (y) . The function ϕ is often called the scale function.

Whenever ϕ (0) > -∞ and ϕ (1) < +∞, one can choose the integration constants defining ϕ (x) so that

ϕ (x) = x 0 e -2 y 0 f (z) g 2 (z) dz dy 1 0 e -2 y 0 f (z)
g 2 (z) dz dy , with ϕ (0) = 0 and ϕ (1) = 1. In this case, the state-space of (y t ; t ≥ 0) is again [0, 1] , the same as for (x t ; t ≥ 0) . Finally, considering the random time change t → θ t with inverse: θ → t θ defined by θ t θ = θ and θ = t θ 0 g 2 (y s ) ds, the novel diffusion (w θ := y t θ ; θ ≥ 0) is easily checked to be identical in law to a standard Brownian motion on ϕ (I). The random time t θ can be expressed as

t θ = θ 0 m ϕ -1 (w τ ) ϕ -1 ′ (w τ ) dτ
where m (x) := 1/ g 2 ϕ ′ (x) is the (positive) speed density at x = ϕ -1 (y). Both the scale function ϕ and the speed measure dµ = m (x) • dx are therefore essential ingredients to reduce the original stochastic process (x t ; t ≥ 0) to the standard Brownian motion (w θ ; θ ≥ 0). The Kolmogorov backward infinitesimal generator G may then be written in Feller form

G (•) = 1 2 d dµ d dϕ • .
Examples (from population genetics):

• Assume f (x) = 0 and g 2 (x) = x (1 -x). This is the neutral WF model discussed at length later. This diffusion is already in natural scale and ϕ

(x) = x, m (x) = [x (1 -x)] -1
. The speed measure is not integrable.

• With π 1 , π 2 > 0, assume f (x) = π 1 -(π 1 + π 2 )
x and g 2 (x) = x (1 -x). This is the WF model with mutation. The parameters π 1 , π 2 can be interpreted as mutation rates. The drift vanishes when x = π 1 /π (where π := π 1 + π 2 is the total mutation pressure) which is an attracting point for the dynamics. Here:

ϕ ′ (y) = ϕ ′ (y 0 ) y -2π1 (1 -y) -2π2 , ϕ (x) = ϕ (x 0 ) + ϕ ′ (y 0 ) x x0 y -2π1 (1 -y) -2π2 dy, with ϕ (0) = -∞ and ϕ (1) = +∞ if π 1 , π 2 > 1/2. The speed measure density is m (x) ∝ x 2π1-1 (1 -x)
2π2-1 and so is always integrable. After normalization to 1, m (x) is the beta(2π 1 , 2π 2 ) density.

• With σ ∈ R, assume a model with quadratic logistic drift f (x) = σx (1 -x) and local variance g 2 (x) = x (1 -x). This is the WF model with selection or selection. For this diffusion (see [START_REF] Kimura | On the probability of fixation of mutant genes in a population[END_REF]), ϕ (x) = 1-e -2σx

1-e -2σ and m (x) ∝ [x (1 -x)] -1 e 2σx is not integrable. Here, σ is a selection or fitness parameter.

Time change and subordination. We start from the diffusion (1) with infinitesimal generator

G = f ∂ x + 1 2 g 2 ∂ 2
x and consider the time change problem without passing first in natural coordinate. Let the random time change

t → θ t = t 0 g 2 (x s ) ds.
Its inverse: θ → t θ defined by θ t θ = θ is given by θ = t θ 0 g 2 (x s ) ds.

In this new stochastic time clock, the subordinated diffusion (y θ := x t θ ; θ ≥ 0) obeys the Langevin SDE with potential U (y) := -2

y 0 f (z) g 2 (z) dz dy θ = f g 2 (y θ ) dθ + dw θ , with backward infinitesimal generator G = g -2 G = f g 2 ∂ x + 1 2 ∂ 2
x (See [START_REF] Itô | Diffusion processes and their sample paths[END_REF], pp. 164-169).

We have • θ t = g 2 (x t ) meaning that at each point x t of the former motion, the motion of the path is accelerated or decelerated, depending on the rate g 2 (x t ) ≶ 1.

Note that conversely

• t θ = 1/g 2 (y θ ) . Under the time substitutions, the road maps of the paths of both (x t ; t ≥ 0) and (y θ ; θ ≥ 0) remain exactly the same. If a path of the former process is accelerated or decelerated by its squared volatility g 2 (its local variance) at each locality, then this process boils down to the latter one. Stated differently, if we measure time by the amount of squared volatility accumulated within each of its path, the process (x t ; t ≥ 0) becomes (y θ ; θ ≥ 0), both with statespace I.

2.3.

The transition probability density. Assume that f (x) and g (x) are now differentiable in I. Let then p (x; t, y) stand for the transition probability density function of x t at y given x 0 = x. Then p := p (x; t, y) is the smallest solution to the Kolmogorov forward (Fokker-Planck) equation (KFE):

(2)

∂ t p = G * (p) , p (x; 0, y) = δ y (x)
where

G * (•) = -∂ y (f (y) •) + 1 2 ∂ 2 y g 2 (y)
• is the adjoint of G (G * acts on the terminal value y whereas G acts on the initial value x). The way one can view this partial differential equation (PDE) depends on the type of boundaries that {0, 1} are.

Suppose for example that the boundaries • := 0 or 1 are both exit (or absorbing) boundaries. From the Feller classification of boundaries, this will be the case if ∀y 0 ∈ (0, 1):

(3) (i) m (y) / ∈ L 1 (y 0 , •) and (ii) ϕ ′ (y) y y0 m (z) dz ∈ L 1 (y 0 , •) ,
where a function f (y) ∈ L 1 (y 0 , •) if

• y0 |f (y)| dy < +∞. In this case, a sample path of (x t ; t ≥ 0) can reach • from the inside of I in finite time but cannot reenter. The sample paths are absorbed at •. There is an absorption at • at time τ x,• = inf (t > 0 : x t = • | x 0 = x) and P (τ x,• < ∞) = 1. Whenever both boundaries {0, 1} are absorbing, the diffusion x t should be stopped at τ x := τ x,0 ∧ τ x,1 . Would none of the boundaries {0, 1} be absorbing, then τ x = +∞. This occurs when the boundaries are inaccessible.

Examples of diffusion with exit boundaries are the WF model and the WF model with selection. In the WF model including mutations, the boundaries are entrance boundaries and so are not absorbing.

When the boundaries are absorbing, then p (x; t, y) is a sub-probability. Letting ρ t (x) := 1 0 p (x; t, y) dy, we clearly have ρ t (x) = P (τ x > t). Such models are non-conservative.

For one-dimensional diffusions, the transition density p (x; t, y) is reversible with respect to the speed density ( [START_REF] Karlin | A second course in stochastic processes[END_REF], Chapter 15, Section 13) and so detailed balance holds:

(4) m (x) p (x; t, y) = m (y) p (y; t, x) , 0 < x, y < 1.

The speed density m (y) satisfies G * (m) = 0. It may be written as a Gibbs measure with density: m (y) ∝ 1 g 2 (y) e -U(y) where the potential function U (y) reads:

(5)

U (y) = -2 y 0 f (z) g 2 (z) dz, 0 < y < 1
and with the measure dy g 2 (y) standing for the reference measure. Furthermore, if p (s, x; t, y) is the transition probability density from (s, x) to (t, y), s < t, then -∂ s p = G (p), with terminal condition p (t, x; t, y) = δ y (x) and so p (s, x; t, y) also satisfies the KBE when looking at it backward in time. The Feller evolution semigroup being time-homogeneous, one may as well observe that with p := p (x; t, y), operating the time substitution t -s → t, p itself solves the KBE (6)

∂ t p = G (p) , p (x; 0, y) = δ y (x) .
In particular, integrating over y,

∂ t ρ t (x) = G (ρ t (x)), with ρ 0 (x) = 1 (x ∈ (0, 1)).
p (x; t, y) being a sub-probability, we may define the normalized conditional probability density q (x; t, y) := p (x; t, y) /ρ t (x), now with total mass 1. We get

∂ t q = -∂ t ρ t (x) /ρ t (x) • q + G * (q) , q (x; 0, y) = δ y (x) .
The term b t (x) := -∂ t ρ t (x) /ρ t (x) > 0 is the time-dependent birth rate at which mass should be created to compensate the loss of mass of the original process due to absorption of (x t ; t ≥ 0) at the boundaries. In this creation of mass process, a diffusing particle started in x dies at rate b t (x) at point (t, y) where it is duplicated in two new independent particles both started at y (resulting in a global birth) evolving in the same diffusive way 1 . The birth rate function b t (x) depends here on x and t, not on y.

When the boundaries of x t are absorbing, the spectra of both -G and -G * are discrete (see [START_REF] Karlin | A second course in stochastic processes[END_REF] pp. 330): There exist positive eigenvalues (λ k ) k≥1 ordered in ascending sizes and eigenvectors

(v k , u k ) k≥1 of both -G * and -G satisfying -G * (v k ) = λ k v k and -G (y k ) = λ k u k such that, with u k , v k := 1 0 u k (x) v k (x) dx and b k := u k , v k -1 , the spectral representation (7) p (x; t, y) = k≥1 b k e -λ k t u k (x) v k (y)
holds.

Let λ 1 > λ 0 = 0 be the smallest non-null eigenvalue of the infinitesimal generator -G * (and of -G). Clearly, -1 t log ρ t (x) → t→∞ λ 1 and by L' Hospital rule therefore

b t (x) → t→∞ λ 1 .
Putting ∂ t q = 0 in the latter evolution equation, independently of the initial condition x

(8) q (x; t, y) → t→∞ q ∞ (y) = v 1 (y)
,

where v 1 is the eigenvector of -G * associated to λ 1 , satisfying -G * v 1 = λ 1 v 1 .
The limiting probability v 1 /norm (after a proper normalization) is called the quasistationary Yaglom limit law of (x t ; t ≥ 0) conditioned on being currently alive at all time t (see [START_REF] Yaglom | Certain limit theorems of the theory of branching random processes[END_REF]).

2.4.

Additive functionals along sample paths. Let (x t ; t ≥ 0) be the diffusion model defined by (1) on the interval I where both endpoints are assumed absorbing (exit). This process is thus transient and non-conservative. We wish to evaluate the nonnegative additive quantities

α (x) = E τ x 0 c (x s ) ds + d (x τ x ) ,
where the functions c and d are both assumed nonnegative on I and ∂I = {0, 1}. The functional α (x) ≥ 0 solves the Dirichlet problem:

-G (α) = c if x ∈ I α = d if x ∈ ∂I,
and α is a super-harmonic function for G, satisfying -G (α) ≥ 0.

Some examples:

1. Assume c = 1 and d = 0 : here, α = E (τ x ) is the mean time of absorption (average time spent in (0, 1) before absorption).

1 Consider a diffusion process with forward infinitesimal generator G * governing the evolution of p (x; t, y) . Suppose that a sample path of this process has some probability that it will be killed or create a new copy of itself, and that the killing and birth rates d and b depend on the current location y of the path. Then the process with the birth and death opportunities of a path has the infinitesimal generator λ (y) • +G * (•) , where λ (y) = b (y) -d (y). The rate can also depend on t and x.

2. Whenever both {0, 1} are exit boundaries, it is of interest to evaluate the probability that x t first hits [0, 1] (say) at 1, given x 0 = x. This can be obtained by choosing c = 0 and d (•) = 1 (• = 1) .

3. Let y ∈ I and put c = 1 2ε 1 (x ∈ (y -ε, y + ε)) and d = 0. As ε → 0, c converges weakly to δ y (x) and, α =:

g (x, y) = E lim 1 2ε τ x 0 1 (y-ε,y+ε) (x s ) ds = ∞ 0 p (x; s, y) ds is the Green function, solution to: -G (g) = δ y (x) if x ∈ I g = 0 if x ∈ ∂I.
g is therefore the mathematical expectation of the local time at y, starting from x (the sojourn time density at y). The solution is known to be (see [START_REF] Karlin | A second course in stochastic processes[END_REF], pp. 198 or [START_REF] Durrett | Probability models for DNA sequence evolution[END_REF], pp. 280)

g (x, y) = 2m (y) (ϕ (1) -ϕ (y)) ϕ (x) -ϕ (0) ϕ (1) -ϕ (0) if x ≤ y (9) g (x, y) = 2m (y) (ϕ (y) -ϕ (0)) ϕ (1) -ϕ (x) ϕ (1) -ϕ (0) if x ≥ y.
The Green function is of particular interest to solve the general problem of evaluating additive functionals α (x). Indeed, as is well-known, see [START_REF] Karlin | A second course in stochastic processes[END_REF] for example, the integral operator with respect to the Green kernel inverts the second order operator -G leading to

α (x) = I g (x, y) c (y) dy if x ∈ I α = d if x ∈ ∂I.
Under this form, α (x) appears as a potential function and any potential function is super-harmonic. Note that for all harmonic function h ≥ 0 satisfying -G (h) = 0,

α h (x) := I g (x, y) c (y) dy + h (x) is again super-harmonic because -G (α h ) = c ≥ 0.
2.5. Transformation of sample paths (Doob-transform) producing killing and/or branching. In the preceding Subsections, we have dealt with a given process and recalled the various ingredients for the expectations of various quantities of interest, summing over the history of paths. In this setup, there is no distinction among paths with different destinations nor did we allow for annihilation or creation of paths inside the domain before the process reached one of the boundaries. The Doob transform of paths allows to do so.

Consider a one-dimensional diffusion (x t ; t ≥ 0) as in [START_REF] Asmussen | Strong limit theorems for general supercritical branching processes with applications to branching diffusions[END_REF]. Let p (x; t, y) be its transition probability. Let α (x) ≥ 0 as x ∈ [0, 1] .

Define a new transformed stochastic process (x t ; t ≥ 0) by its transition probability

(10) p (x; t, y) = α (y) α (x) p (x; t, y) .
In this construction of (x t ; t ≥ 0) through a change of measure, sample paths x → y of (x t ; t ≥ 0) with a large value of the ratio α (y) /α (x) are favored. This is a selection of paths procedure due to Doob (see [START_REF] Dynkin | Markov processes. Vols. I, II. Translated with the authorization and assistance of the author[END_REF]).

The KFE for p clearly is ∂ t p = G * (p), with p (x; 0, y) = δ y (x) and G * (p) = α (y) G * (p/α (y)). The adjoint Kolmogorov backward operator of the transformed process is therefore by duality ( 11)

G (•) = 1 α (x) G (α (x) •) . Developing, with α ′ (x) := dα (x) /dx and G (•) := α ′ α g 2 ∂ x (•) + G (•), we get (12) G (•) = 1 α G (α) • + G (•) =: λ (x) • + G (•)
and the new KB operator can be obtained from the latter by adding a drift term

α ′ α g 2 ∂ x to
d x t = f ( x t ) dt + g ( x t ) dw t , x 0 = x ∈ (0, 1) ,
possibly killed or branching at rate λ (x) as soon as λ = 0. Whenever ( x t ; t ≥ 0) is killed, it enters conventionally into some coffin state {∂} added to the state-space.

Let us look at special cases:

(i) Suppose α ≥ 0 is such that -G (α) ≥ 0 (By α ≥ 0, we mean α > 0 in I, possibly with α (0) or α (1) equal 0). Then α is called a super-harmonic (or excessive) function for the process with infinitesimal generator G.

In this case, the rate λ (x) =: -d (x) satisfies λ (x) ≤ 0 and only killing occurs at rate d (x). Let τ x be the new absorption time at the boundaries of ( x t ; t ≥ 0) started at x (with τ x = ∞ would the boundaries be inaccessible to the new process x t ). Let τ x,∂ be the killing time of ( x t ; t ≥ 0) started at x (the hitting time of

∂), with τ x,∂ = ∞ if G (α) ≡ 0. Then τ x := τ x ∧ τ x,∂
is the novel stopping time for ( x t ; t ≥ 0) . The SDE for ( x t ; t ≥ 0), together with its global stopping time τ x characterize the new process (x t ; t ≥ 0) with full generator G to consider.

(ii) Suppose α ≥ 0 is such that -G (α) ≤ 0. Then α is called a sub-harmonic function for the process with generator G.

In this case, the rate λ (x) =: b (x) satisfies λ (x) ≥ 0 and only branching occurs at rate b (x). The transformed process (with infinitesimal backward generator G) accounts for a branching diffusion where a diffusing mother particle (with generator G and started at x) lives a random exponential time with constant rate 1. When the mother particle dies, it gives birth to a spatially dependent random number M (x) of particles, with mean µ (x) = 1 + λ (x) (where M (x)

d = 1 + ∆ (λ (x)
) and ∆ (λ (x)) is a geometrically distributed random variable on {0, 1, 2, ...} with mean λ (x)). Then M (x) independent daughter particles are started afresh where their mother particle died, with the event M (x) = 0 impossible; they move along a diffusion governed by G and reproduce, independently and so on for the subsequent particles.

(iii) If λ (x) =: b (x) is bounded above, λ may be put under the alternative form

λ (x) = λ * (µ (x) -1) ,
where λ * = sup x∈[0,1] λ (x) and 1 ≤ µ (x) ≤ 2. In this case, we can assume that M (x) can only take the values 1 or 2 with probability p 1 (x) and p 2 (x) respectively, with

p 1 (x) + p 2 (x) = 1. Then, µ (x) = E (M (x)) = p 1 (x) + 2p 2 (x) = 1 + p 2 (x) and λ (x) = λ * p 2 (x) . Note that µ (x)-1 = p 2 (x) > p 1 (x) = 2-µ (x)
. We get a binary branching process at rate λ * with the event to produce two particles being more likely than the one to produce a single one, whatever is x.

(iv) Whenever α is such that -G (α) has no constant sign, then both killing and branching can simultaneously occur at the death of the mother particle. λ (x) may be put under the form λ (x) = b (x) -d (x) where b (x) and d (x) are the birth (branching) and death (killing) components of λ (x).

(v) Suppose λ (x) is bounded below and let λ * = -inf x∈[0,1] λ (x) > 0. Then one may view λ as λ (x) = λ * (µ (x) -1)
, where µ (x) ≥ 0. In this case, branching occurs at rate λ * . When the mother particle dies, it gives birth to a spatially dependent random number M (x) of particles (where M (x) d = ∆ (µ (x)) and ∆ (µ (x)) is a geometrically distributed random variable on {0, 1, 2, ...} with mean µ (

x) = 1 + λ (x) /λ * ). With p m (x) = P (M (x) = m) = p 0 (x) q (x) m , m ≥ 0, p 0 (x) = 1 1+µ(x) , q 0 (x) = 1 -p 0 (x) λ (x) = λ *   m≥1 mp m (x) -1   = λ *   m≥2 (m -1) p m (x) -p 0 (x)   .
Thus, the decomposition λ (x) = b (x) -d (x) holds, where b and d can be read from

λ (x) = λ * µ (x) 2 1 + µ (x) - 1 1 + µ (x)
.

(vi) In some other examples, the killing/branching rate λ = G (α) /α is bounded above and below. Then λ may be put under the form

λ (x) = λ * (µ (x) -1)
,

where λ * = sup x∈[0,1] |λ (x)| and 0 ≤ µ (x) ≤ 2.
In this case, we can assume that M (x) can only take the values 0 or 2 with probability p 0 (x) and p 2 (x) respectively, with

p 0 (x) + p 2 (x) = 1. Then, µ (x) = E (M (x)) = 2p 2 (x) and λ (x) = λ * (2p 2 (x) -1) = λ * (p 2 (x) -p 0 (x)) ,
giving a simple decomposition of λ in the form λ (x) = b (x) -d (x) with the mother particle living a random exponential time now with constant rate λ * before giving birth to none or two descending particles (a binary branching process). Note that p 0 (x) ≥ p 2 (x) (respectively p 0 (x) ≥ p 2 (x)) when µ (x) ≤ 1 (µ (x) ≥ 1).

Examples of α. When (x t ; t ≥ 0) is non-conservative, consider

α (x) = E τ x 0 c (x s ) ds + d (x τ x ) ,
where the functions c and d are both assumed nonnegative on I and ∂I = {0, 1}.

Then α ≥ 0 solves the Dirichlet equation -Gα (x) = c (x) ≥ 0 on I (= d (x) on ∂I) and so α is super-harmonic or excessive. We refer to [START_REF] Huillet | On Wright-Fisher diffusion and its relatives[END_REF] for examples of Doob transforms based on such super-harmonic functions allowing to understand various conditionings of interest when the starting point process (x t ; t ≥ 0) is a neutral WF diffusion or a WF diffusion with selection.

Whenever α is super-harmonic for G, then

β = 1/α ≥ 0 is sub-harmonic for G = G + α ′ α g 2 ∂ x .
This results from the obvious identity

β -1 G (β) = -α -1 G (α) ,
showing that -Gα ≥ 0 entails -G (β) ≤ 0.

Whenever (x t ; t ≥ 0) is conservative and ergodic

1 t E x t 0 c (x s ) ds → t→∞ µ (c) := 1 0 c (y) dµ (y)
where dµ = m (y) dy is the invariant probability measure of (x t ; t ≥ 0) . Define We conclude that α is G-super-harmonic if ever c (x) ≥ µ (c), ∀x. ⋄ Background (multiplicative functional and path integral). The Doob transforms used here are particular instances of more general transformations based on multiplicative functionals. Let x t be the diffusion process [START_REF] Asmussen | Strong limit theorems for general supercritical branching processes with applications to branching diffusions[END_REF] governed by

G = f ∂ x + 1 2 g 2 ∂ 2
x with x 0 = x. Define the multiplicative functional M t as the solution of the differential equation

dM t = M t • (a (x t ) dt + b (x t ) dw t ) , M 0 = 1,
where a and b are arbitrary twice differentiable functions. Integrating, we get

M t = e t 0 (a-1 2 b 2 )(xs)ds+ t 0 b(xs)dws .
Let B be a Borel subset of I. Define a new process whose density p is obtained after a modification of the original one while using the multiplicative modulation factor M t as B p (x; t, y) dy

:= E x [M t 1 (x t ∈ B)] = B E x [M t | x t = y] p (x; t, y) dy.
Integrating M t over paths with fixed two endpoints x and y, E x [M t | x t = y] can be interpreted as the Radon-Nykodym derivative of p with respect to p, the density of x t . By duality, let

v (x, t) = E x [M t ψ (x t )] , v (x, 0) = ψ (x) .
Applying Itô calculus, we get

∂ t v = G (v) = (G + gb∂ x + a) (v) =: G + a (v) ,
where the modified backward infinitesimal generator G is obtained by adding a drift term gb∂ x to G to produce G and a multiplicative part a•. The adjoint KFE giving the evolution of p is thus

∂ t p = G * (p) = G * + a (p) , p (x; 0, y) = δ y (x) .
-(Cameron-Martin-Girsanov) For instance, when a = 0 and b = -f /g, the generator of the transformed diffusion is G = 1 2 g 2 ∂ 2

x killing the drift term of the original process governed by G. In this case,

M t = e -1 2 t 0 ( f g ) 2 (xs)ds+2 f g (xs)dws .
Clearly in this case M t is a martingale with E x (M t ) = 1, assuming b to be bounded. This construction kills the drift of the original process while using a change of measure.

-(Feynman-Kac) When b = 0, the generator of the transformed diffusion is G = G + a adding a multiplicative component a to the one G governing the original process. In this case

M t = e t 0 a(xs)ds is the exponential of the integrated rate. If v (x, t) = E x [M t ψ (x t )], v (x, 0) = ψ (x) , then v solves ∂ t v = G (v) = (G + a) (v) , v (x, 0) = ψ (x) . In particular, if v (x, t) = E x [M t ], v (x, 0) = 1 (x ∈ (0, 1)) , then v solves ∂ t v = G (v) = (G + a) (v) , v (x, 0) = 1 (x ∈ (0, 1)) .
-(Doob) Suppose now

dM t = M t α -1 (x t ) dα (x t ) , M 0 = 1.
This M t is a particular instance of the general M t introduced above. Indeed, applying Itô calculus,

α -1 (x t ) dα (x t ) = α -1 α ′ [f dt + gdw] + 1 2 α -1 α ′′ g 2 dt,
leading to

a = α -1 f α ′ + g 2 2 α ′′ = G (α) /α =: λ (x) b = α -1 α ′ g. Thus G = G + gb∂ x + a = G + α -1 α ′ g 2 ∂ x + λ (x)
as already observed earlier.

Now, from the differential generation of M t ,

M t = α (x t ) α (x) , M 0 = 1
only depends on the terminal and initial values of (x s ; 0 ≤ s ≤ t) and not on its intermediate values (such a particular Doob transformation is thus a gauge). Thus here

E x [M t | x t = y] = α(y) α(x)
consistently with the definition p (x; t, y) = α(y) α(x) p (x; t, y) .

A super-harmonic example.

Although this work chiefly focuses on Doob-transforms where branching is present in λ, let us give a significant example where the Doob transform just produces killing like in (i). Suppose (x t ; t ≥ 0) is a non-conservative diffusion. Let λ 1 be the smallest non-null eigenvalue of the infinitesimal generator G of (x t ; t ≥ 0) . Let α = u 1 be the corresponding eigenvector, that is satisfying -Gu 1 = λ 1 u 1 ≥ 0 with boundary conditions u 1 (0) = u 1 (1) = 0. Then c = λ 1 u 1 . The new KB operator associated to the transformed process (x t ; t ≥ 0) is ( 14)

G (•) = 1 α G (α) • + G (•) = -λ 1 • + G (•) ,
obtained while killing the sample paths of the process ( x t ; t ≥ 0) governed by G at constant death rate d = λ 1 . The transition probability of the transformed stochastic process (x t ; t ≥ 0) is p (x; t, y) = u 1 (y) u 1 (x) p (x; t, y) .

Define p (x; t, y) = e λ 1 t p (x; t, y) . It is the transition probability of the process ( x t ; t ≥ 0) governed by G; it corresponds to the original process (x t ; t ≥ 0) conditioned on never hitting the boundaries {0, 1} (the so-called Q-process of (x t ; t ≥ 0), see [START_REF] Lambert | Population dynamics and random genealogies[END_REF]). It is simply obtained from (x t ; t ≥ 0) by adding the additional drift term

u ′ 1 u1 g 2 to f
, where u 1 is the eigenvector of G associated to its smallest non-null eigenvalue. The determination of α = u 1 is a Sturm-Liouville problem. When t is large, to the dominant order

p (x; t, y) ∼ e -λ1t u 1 (x) v 1 (y) u 1 , v 1 ,
where v 1 is the Yaglom limit law of (x t ; t ≥ 0) . Therefore (15) p (x; t, y) ∼ e λ1t u 1 (y)

u 1 (x) e -λ1t u 1 (x) v 1 (y) u 1 , v 1 = u 1 (y) v 1 (y) u 1 , v 1 .
Thus the limit law of the Q-process ( x t ; t ≥ 0) is the normalized Hadamard product of the eigenvectors u 1 and v 1 associated respectively to G and G * . On the other hand, the limit law of ( x t ; t ≥ 0) is directly given by ( 16) p (x; t, y)

→ t→∞ p (y) = 1 Zg 2 (y) e 2 y 0 f (z)+ u ′ 1 u 1 g 2 (z) g 2 (z) dz = u 2 1 (y) Zg 2 (y) e 2 y 0 f (z) g 2 (z) dz ,
where Z is the appropriate normalizing constant. Comparing ( 15) and ( 16)

v 1 (y) = u 1 (y) g 2 (y) e 2 y 0 f (z) g 2 (z) dz = u 1 (y) m (y) .
The eigenvector v 1 associated to G * is therefore equal to the eigenvector u 1 associated to G times the speed density of (x t ; t ≥ 0) .

When dealing for example with the neutral WF diffusion (see Section 4 for additional details), it is known that λ 1 = 1 with u 1 = x (1 -x) and v 1 ≡ 1. The Q-process ( x t ; t ≥ 0) in this case obeys ( 17)

d x t = (1 -2 x t ) dt + x t (1 -x t )dw t ,
with an additional stabilizing drift toward 1/2:

f (x) = u ′ 1 u1 g 2 (x) = 1 -2x.
The limit law of the Q-process ( x t ; t ≥ 0) in this case is 6y (1 -y). The latter conditioning is more stringent than the Yaglom conditioning and so the limiting law has more mass away from the boundaries (compare with the uniform quasistationary Yaglom limit (8) with v 1 ≡ 1).

The Wright-Fisher example

In this Section, we briefly and informally recall that the celebrated WF diffusion process with or without a drift may be viewed as a scaling limit of a simple two alleles discrete space-time branching process preserving the total number N of individuals in the subsequent generations (see [START_REF] Karlin | A second course in stochastic processes[END_REF], [START_REF] Ethier | Markov processes. Characterization and convergence[END_REF], for example).

P (ν N = k N ) = N ! • N -N N n=1 k n ! .
This distribution can be obtained by conditioning N independent Poisson distributed random variables on summing to N . Assume subsequent iterations of this reproduction law are independent so that the population is with constant size for all generations.

Let N r (n) be the offspring number of the n first individuals at the discrete generation r ∈ N 0 corresponding to (say) allele A 1 (the remaining number N -N r (n) counts the number of alleles A 2 at generation r). This sibship process is a discretetime Markov chain with binomial transition probability given by:

P (N r+1 (n) = k ′ | N r (n) = k) = N k ′ k N k ′ 1 - k N N -k ′ .
Assume next that n = [N x] where x ∈ (0, 1) . Then, as well-known, the dynamics of the continuous space-time re-scaled process x t := N [N t] (n) /N , t ∈ R + can be approximated for large N , to the leading term in N -1 , by a Wright-Fisher-Itô diffusion on [0, 1] (the purely random genetic drift case):

(18)

dx t = x t (1 -x t )dw t , x 0 = x.
Here (w t ; t ≥ 0) is a standard Wiener process. For this scaling limit process, a unit laps of time t = 1 corresponds to a laps of time N for the original discrete-time process; thus time is measured in units of N . If the initial condition is x = N -1 , x t is the diffusion approximation of the offspring frequency of a singleton at generation

[N t].
Equation ( 18) is a 1-dimensional diffusion as in ( 1) on [0, 1] , with zero drift f (x) = 0 and volatility g (x) = x (1 -x). This diffusion is already in natural coordinate and so ϕ (x) = x. The scale function is x and the speed measure [x (1 -x)] -1 dx. One can check that both boundaries are exit in this case: The stopping time is τ x = τ x,0 ∧ τ x,1 where τ x,0 is the extinction time and τ x,1 the fixation time. The corresponding infinitesimal generators are

G (•) = 1 2 x (1 -x) ∂ 2 x (•) and G * (•) = 1 2 ∂ 2 y (y (1 -y) •) .

3.2.

Non-neutral cases. Two alleles WF models (with non-null drifts) are classically obtained by considering the binomial transition probabilities bin(N, p N ) :

P (N r+1 (n) = k ′ | N r (n) = k) = N k ′ p N k N k ′ 1 -p N k N N -k ′ where p N (x) : x ∈ (0, 1) → (0, 1)
is now some state-dependent probability (which is different from the identity x) reflecting some deterministic evolutionary drift from the allele A 1 to the allele A 2 .

For each r, we have

E (N r+1 (n) | N r (n) = k) = N p N k N σ 2 (N r+1 (n) | N r (n) = k) = N p N k N 1 -p N k N
which is amenable to a diffusion approximation in terms of

x t := N [N t] (n) /N , t ∈ R + under suitable conditions.
For instance, taking

p N (x) = (1 -π 2,N ) x + π 1,N (1 -x)
where (π 1,N , π 2,N ) are small (N -dependent) mutation probabilities from A 2 to A 1 (respectively

A 1 to A 2 ). Assuming (N • π 1,N , N • π 2,N ) → N →∞
(π 1 , π 2 ), leads after scaling to the drift of WF model with positive mutations rates (π 1 , π 2 ).

Taking

p N (x) = (1 + s 1,N ) x 1 + s 1,N x + s 2,N (1 -x)
where s i,N > 0 are small N -dependent selection parameter satisfying N •s i,N → N →∞ σ i > 0, i = 1, 2, leads, after scaling, to the WF model with selective drift f (x) = σx (1 -x), where σ := σ 1 -σ 2 . Typically, the drift f (x) is a large N approximation of the bias: N (p N (x) -x) . The WF diffusion with selection is thus:

(19) dx t = σx t (1 -x t ) dt + x t (1 -x t )dw t
where time is measured in units of N. Letting θ t = N t define a new time-scale with inverse t θ = θ/N , the time-changed process y θ = x θ/N now obeys the SDE

dy θ = sy θ (1 -y θ ) dθ + 1 N y θ (1 -y θ )dw θ ,
with a small diffusion term. Here s = s 1 -s 2 and time θ is the usual time-clock.

The WF diffusion with selection [START_REF] Mandl | Analytical treatment of one-dimensional Markov processes[END_REF] tends to drift to • = 1 (respectively

• = 0) if allele A 1 is selectively advantageous over A 2 : σ 1 > σ 2 (respectively σ 1 < σ 2 )
in the following sense: if σ > 0 (respectively < 0), the fixation probability at • = 1, which is [START_REF] Kimura | On the probability of fixation of mutant genes in a population[END_REF]

P (τ x,1 < τ x,0 ) = 1 -e -2σx
1 -e -2σ , increases (decreases) with σ taking larger (smaller) values.

The usual way to look at the WF diffusion with mutation and selection is to compose the two above mechanisms p N (x) corresponding to mutation and selection respectively. In the scaling limit, one obtains the standard WF diffusion model including mutations and selection as: [START_REF] Mano | Duality, ancestral and diffusion processes in models with selection[END_REF]

dx t = [(π 1 -(π 1 + π 2 ) x t ) + σx t (1 -x t )] dt + x t (1 -x t )dw t .

The neutral WF model

In this Section, we particularize the general ideas developed in the introductory Section 2 to the neutral WF diffusion [START_REF] Lambert | Population dynamics and random genealogies[END_REF] and draw some straightforward conclusions most of which are known which illustrate the use of Doob transforms.

4.1.

Explicit solutions of the neutral KFE. As shown by Kimura in ( [START_REF] Kimura | Stochastic processes and distribution of gene frequencies under natural selection[END_REF]), the Kolmogorov forward (and backward) equation is exactly solvable in this case, using spectral theory. The solutions involve a series expansion in terms of eigen-functions of the KB infinitesimal generator with discrete eigenvalues spectrum.

Let

λ k = k (k + 1) /2, k ≥ 0. There exist u k = u k (x) and v k = v k (y) solving the eigenvalue problem: -G (u k ) = λ k u k and -G * (v k ) = λ k v k . With v k , u k = 1 0 u k (x) v k (x)
dx, the transition probability density p (x; t, y) of the neutral WF models admits the spectral expansion

p (x; t, y) = k≥1 b k e -λ k t u k (x) v k (y) where b k = 1 v k , u k
where u k (x) are the Gegenbauer polynomials rescaled on [0, 1] and normalized to have value 1 at x = 0. In particular, u 0

(x) = x, u 1 (x) = x -x 2 , u 2 (x) = x-3x 2 +2x 3 , u 3 (x) = x-6x 2 +10x 3 -5x 4 , u 4 (x) = x-10x 2 +30x 3 -35x 4 +14x 5 , ... Next, v k (y) = m (y) u k (y) where m (y) = 1/ (y (1 -y)) is the speed density of the neutral WF diffusion. For instance, v 0 (y) = 1 1-y , v 1 (y) = 1, v 2 (y) = 1 -2y, v 3 (y) = 1 -5y + 5y 2 , v 4 (y) = 1 -9y + 21y 2 -14y 3 ,...
Although λ 0 = 0 really constitutes an eigenvalue, only v 0 (y) is not a polynomial and the spectral expansion of p should start at k = 1, expressing that p is a subprobability. When k ≥ 1, from their definition, the u k (x) polynomials satisfy

u k (0) = u k (1) = 0 in such a way that v k (y) = m (y) • u k (y) , k ≥ 1 is a polynomial with degree k -1.
The series expansion for p (x; t, y) solves the KFE of the WF model.

We have P (τ x > t) = 1 0 P (x t ∈ dy) and so

ρ t (x) := P (τ x > t) = k≥1 1 0 v k (y) dy v k , u k e -λ k t u k (x)
is the exact tail distribution of the absorption time.

Since v 1 (y) = 1, to the leading order in t, for large time

P (x t ∈ dy) = 6e -t • x (1 -x) dy + O e -3t
which is independent of y. Integrating over y, ρ t (x) := P (τ x > t) ∼ 6e -t • x (1 -x) so that the conditional probability

4.2.

Additive functionals for the neutral WF and Doob transforms. Let (x t ; t ≥ 0) be the neutral WF diffusion model defined by [START_REF] Lambert | Population dynamics and random genealogies[END_REF] on the interval I = [0, 1] where both endpoints are absorbing (exit). Consider the additive quantities

α (x) = E τ x 0 c (x s ) ds + d (x τ x ) ,
where functions c and d are both nonnegative.

With G = 1 2 x (1 -x) ∂ 2 x , α (x) solves: -G (α) = c if x ∈ I α = d if x ∈ ∂I.
Therefore α is a super-harmonic function for G.

Take c = lim ε↓0 1 2ε 1 (y-ε,y+ε) (x) =: δ y (x) and d = 0, when y ∈ I : in this case, α := g (x, y) is the Green function (the mean local time at y given the process started at x). The solution takes the simple form

g (x, y) = 2 x y if x < y g (x, y) = 2 1 -x 1 -y if x > y.
The Green function solves the above general problem of evaluating additive functionals α (x):

α (x) = I g (x, y) c (y) dy if x ∈ I α = d if x ∈ ∂I.
There are many interesting choices of c therefore leading to α, allowing to compute for example the mean time till absorption for the neutral WF diffusion, the probability to hit state 1 before 0... For each choice of α, it is interesting to study the transformed process (x t ; t ≥ 0) whose transition probability is given by

p (x; t, y) = α (y) α (x) p (x; t, y) ,
in terms of the original process transition probability p (x; t, y) . This allows for example to understand the neutral WF process conditioned on exit at some boundary and to evaluate for this new process interesting average additive functionals such as the mean time needed to hit the exit boundary...For detailed similar examples arising in the context of WF diffusions and related ones, see [START_REF] Huillet | On Wright-Fisher diffusion and its relatives[END_REF].

The WF model with selection

Now we briefly focus on the diffusion process [START_REF] Mandl | Analytical treatment of one-dimensional Markov processes[END_REF]. Let (v k (y)) k≥1 be the Gegenbauer eigen-polynomials of the KF operator corresponding to the neutral WF diffusion [START_REF] Lambert | Population dynamics and random genealogies[END_REF], so with eigenvalues λ k = k (k + 1) /2, k ≥ 1. Define the oblate spheroidal wave functions on [0, 1] as

(22) w σ k (y) = ′ l≥1 f l k v l (y) ,
where f l k obey the three-term recurrence defined in [START_REF] Mano | Duality, ancestral and diffusion processes in models with selection[END_REF]. In the latter equality, the l summation is over odd (even) values if k is even (odd).

Define v σ k (y) = e σy w σ k (y) and u σ k (x) = 1 m(x) v σ k (x) where m (x) = e 2σx / (x (1 -x))
is the speed measure density of the WF model with selection [START_REF] Mandl | Analytical treatment of one-dimensional Markov processes[END_REF].

The system (u σ k (x) , v σ k (x)) k≥1 constitute a system of eigen-functions for the WF with selection generators -G and -G * with eigenvalues λ σ k implicitly defined in [START_REF] Mano | Duality, ancestral and diffusion processes in models with selection[END_REF], thus with -G (u

σ k ) = λ σ k u σ k and -G * (v σ k ) = λ σ k v σ k .
The eigen-function expansion of the transition probability density of the WF model with selection is thus, [START_REF] Kimura | Diffusion Models in Population Genetics[END_REF]:

(23) p (x; t, y) = k≥1 b σ k e -λ σ k t u σ k (x) v σ k (y) where b σ k = v σ k , u σ k -1 .
The WF model with selection can be viewed as a perturbation problem of the neutral WF model (see [START_REF] Maruyama | Stochastic problems in population genetics[END_REF]). There exist perturbation developments of λ σ k around λ k with respect to σ 2 , [START_REF] Kimura | Diffusion Models in Population Genetics[END_REF]. They are valid and useful for small σ.

The WF diffusion process x t with selection ( 19) is non-conservative, with finite hitting time τ x of one of the boundaries. Following the general arguments developed in Section 2, the Yaglom limit of x t conditioned on τ x > t is the normalized version of (24) v σ 1 (y) = e σy w σ 1 (y) .

The limit law of x t conditioned on never hitting the boundaries in the remote future is the normalized version of (25)

u σ 1 (y) v σ 1 (y) = 1 m (y) v σ 1 (x) 2 = y (1 -y) w σ 1 (y) 2 .
Because the latter conditioning is more stringent than the former, the probability mass of ( 25) is more concentrated inside the interval than (24). Compare with the statements at the end of Section 2 concerning the neutral WF diffusion.

Doob transform of the neutral WF model: sub-critical BD

In this Section, we define the branching WF diffusion model with selection while applying a Doob transform to the neutral WF model, based on the sub-harmonic additive functional α (x) = e σx , say with σ > 0. We then study in detail the obtained branching process.

The starting point is thus the neutral WF diffusion:

dx t = x t (1 -x t )dw t , x 0 = x ∈ (0, 1) . For this model, G = 1 2 x (1 -x) ∂ 2
x and both boundaries are exit. With λ k = k (k + 1) /2, k ≥ 0, its transition density p (x; t, y) admits the spectral representation

(26) p (x; t, y) = k≥1 b k e -λ k t u k (x) v k (y) ,
in terms of the Gegenbauer eigen-polynomials (see Subsection 4.1). We shall consider the following transformation of paths on the neutral WF model: Let α (x) = e σx , σ > 0 and consider

G (•) = α -1 G (α•) = G (•) + b (x) • . We now have G (α) = 1 2 σ 2 x (1 -x) e σx and so b (x) = G (α) /α = σ 2 2 x (1 -x) ≥ 0. Note that -G (α) ≤ 0 indicating that α is sub-harmonic for G.
In this case study, one selects sample paths of (x t ; t ≥ 0) with large α (y) and we claim that this is an alternative interesting way to introduce selection in the neutral WF diffusion process.

The dynamics of ( x t ; t ≥ 0) governed by G is easily seen to be the standard WF with selection dynamics [START_REF] Mandl | Analytical treatment of one-dimensional Markov processes[END_REF]

d x t = σ x t (1 -x t ) dt + x t (1 -x t )dw t , subject to additional quadratic branching at rate b (x) = 1 2 σ 2 x (1 -x) inside I. We indeed have G (•) = e -σx G (e σx •) = b (x) • + G (•) , where G =: f ∂ x + 1 2 g 2 ∂ 2 x = σx (1 -x) ∂ x + 1 2 x (1 -x) ∂ 2 x
is the KBE operator of the dynamics ( x t ; t ≥ 0) . Recall that x t is transient and so hits one of the boundaries {0, 1} in finite time τ x .

To summarize, in our branching diffusion way to look at selection, we move from the neutral WF diffusion (x t ; t ≥ 0) to the standard WF diffusion with selection ( x t ; t ≥ 0) but subject to additional branching at rate b (x) .

Remark. With β (x) := α (x) -1 = e -σx , we clearly have

G (β (x)) = 0
and β is an harmonic function for G and as a result, Doob-transforming G by β, we get

β -1 G (β•) = (αβ) -1 G (αβ•) = G (•)
which is the infinitesimal generator of the original neutral WF martingale. ⋄

The birth (creating) rate b ≥ 0 in G is bounded from above on (0, 1). It may be put into the canonical form b

(x) = b * (µ (x) -1) where b * = max x∈[0,1] (b (x)) = σ 2 8 > 0 and (27) µ (x) = 1 + 4x (1 -x) ,
whose range is the interval [START_REF] Asmussen | Strong limit theorems for general supercritical branching processes with applications to branching diffusions[END_REF][START_REF] Asmussen | Some modified branching diffusion models[END_REF]

as x ∈ [0, 1] .
The density of the transformed process is p (x; t, y) = α(y) α(x) p (x; t, y) . It is exactly known because so is p is from (26).

The transformed process (with infinitesimal backward generator G) accounts for a branching diffusion (BD) where a diffusing mother particle (with generator G and started at x) lives a random exponential time with constant rate b * . When the mother particle dies, it gives birth to a spatially dependent random number M (x) of particles (with mean µ (x)). M (x) independent daughter particles are started where their mother particle died; they move along a WF diffusion with selection and reproduce, independently and so on.

Because µ (x) is bounded above by 2 and larger than 1 (indicating a super-critical branching process), we actually get a BD with binary scission whose random offspring number satisfies ('w.p.' meaning 'with probability')

M (x) = 0 w.p. p 0 (x) = 0 M (x) = 1 w.p. p 1 (x) = 2 -µ (x)
M (x) = 2 w.p. p 2 (x) = µ (x) -1, with p 2 (x) ≥ p 1 (x) (the event that 2 particles are generated in a splitting event is more probable than a single one).

For such a transformed process, the trade-off is as follows: there is a competition between the boundaries {0, 1} which are absorbing for the particle system and the number of particles N t (x) in the system at each time t, which may grow due to binary branching events (or remain steady when M (x) = 1).

The density p of the transformed process has the following interpretation

(28) p (x; t, y) = E   Nt(x) n=1 p (n) (x; t, y)   ,
where p (n) (x; t, y) is the density at (t, y) of the nth alive particle descending from the ancestral one (Eve), started at x. In the latter formula, the sum vanishes if N t (x) = 0. A particle is alive at time t if it came to birth before t and has not been yet absorbed by the boundaries.

Let ρ t (x) = (0,1) p (x; t, y) dy. Then ρ t (x) is the expected number of particle alive at time t. We have

∂ t ρ t (x) = G (ρ t (x)) , ρ 0 (x) = 1 (x ∈ (0, 1)) .
Remark. From the Feynman-Kac formula, p in (28) is also But then q (x; t, y) := p (x; t, y) /ρ t (x) obeys the forward PDE ∂ t q (x; t, y) = -∂ t ρ t (x) ρ t (x) + b (y) q (x; t, y) + G * (q (x; t, y))

as a result of ∂ t p (x; t, y) = G * (p (x; t, y)). We have (29) q (x; t, y) =

E Nt(x) n=1 p (n) (x; t, y) E [N t (x)]
showing that q (x; t, y) is the average presence density at (t, y) of the system of particles all descending from Eve started at x.

Clearly -log ρ t (x) t → t→∞ λ 1 = 1 (and therefore also -∂tρ t (x) ρ t (x) by L' Hospital rule), because

ρ t (x) = 1 α (x) k≥1 b k e -λ k t u k (x) 1 0 α (y) v k (y) dy.
The expected number of particles in the system decays globally and exponentially at rate λ 1 .

The BD transformed process therefore admits an integrable Yaglom limit q ∞ , solution to

-G * (q ∞ ) = (λ 1 + b (y)) q ∞ or -G * (q ∞ ) = λ 1 q ∞ . With v 1 (y) = 1, the
first eigenvector of -G * associated to the smallest positive eigenvalue λ 1 = 1, q ∞ is of the product form (30) q ∞ (y) = C * e σy v 1 (y) = σe σy e σ -1 .

The arbitrary multiplicative constant C * was chosen in such a way that q ∞ (y) is a probability.

By analogy with the Yaglom construction, this limiting probability q ∞ can be called the quasi-stationary Yaglom average density at (t, y) for the BD particle system (it is also the ground state for G * ).

There is also a natural eigenvector φ ∞ of the backward operator -G, satisfying -G φ ∞ = λ 1 φ ∞ (the ground state for G). It is explicitly here

(31) φ ∞ (x) = C α (x) u 1 (x) = 6 (e σ -1) σ e -σx x (1 -x) .
The arbitrary multiplicative constant C = 6/C * was chosen in such a way that 1 0 q ∞ (y) φ ∞ (y) dy = 1. Note that the spectral structures of both G * and G are easily obtainable from the ones of G * and G thanks to the Doob transform structure.

In the terminology of [START_REF] Pinsky | Positive harmonic functions and diffusion[END_REF], both operators G (•) + λ 1 • and its adjoint are critical 2 . In this context, the constant λ 1 is called the generalized principal eigenvalue. The eigen-functions φ ∞ , q ∞ are their associated ground states. We note that we have the L 1 -product property (See [START_REF] Pinsky | Positive harmonic functions and diffusion[END_REF], Subsection 4.9).

1 0 φ ∞ (x) q ∞ (x) dx = 6 1 0 u 1 (x) v 1 (x) dx = 1 < ∞.
Remark. Using the Feynman-Kac representation of ρ t (x) , we get We have the x log x condition:

- 1 t log E x e
(32)

1 0 l (x) φ ∞ (x) q ∞ (x) dx = 48 log 2 1 0 x (1 -x) u 1 (x) v 1 (x) dx < ∞.
We conclude (following [START_REF] Asmussen | Strong limit theorems for general supercritical branching processes with applications to branching diffusions[END_REF] and [START_REF] Asmussen | Some modified branching diffusion models[END_REF]) that, as a result of the condition (32) being trivially satisfied, global extinction holds in the following sense:

(i) P (N t (x) = 0) → t→∞ 1, uniformly in x.
(ii) there exists a constant γ > 0 :

e λ 1 t [1 -P (N t (x) = 0)] → t→∞ γφ ∞ (x) , uniformly in x. 2 G (•)+λ 1 • (G * (•)+λ 1 •
) is said to be critical if there exists some function φ ∞ ∈ C 2 (respectively q ∞ ∈ C 2 ), strictly positive in (0, 1) , such that: G φ ∞ + λ 1 φ ∞ = 0 (respectively G * (q ∞ ) + λ 1 q ∞ = 0) and the operators do not possess a minimal positive Green function.

(iii) For all bounded measurable function ψ on I :

E   Nt(x) n=1 ψ x (n) t | N t (x) > 0   → t→∞ γ -1 (0,1)
ψ (y) q ∞ (y) dy.

From (i), it is clear that the process gets ultimately extinct with probability 1. In the trade-off between pure branching and absorption at the boundaries, all particles get eventually absorbed and the global BD process turns out be sub-critical (even though µ (x) = EM (x) > 1 for all x ∈ (0, 1)): Probability mass escapes out of I although the BD survives with positive probability.

In the statement (ii) , the quantity 1 -P (N t (x) = 0) = P (N t (x) > 0) is also P (T (x) > t) where T (x) is the global extinction time of the particle system descending from an Eve particle started at x. The number -λ 1 is the usual Malthus exponential decay rate parameter. From (ii) , φ ∞ (x) has a natural interpretation in terms of the propensity of the particle system to survive to its extinction fate: the so-called reproductive value in demography.

(iii

) with ψ = 1 reads E [N t (x) | N t (x) > 0] → t→∞ γ -1
giving an interpretation of the constant γ (which may be hard to evaluate in practise).

The ground states of G + λ 1 and its adjoint are thus φ ∞ , q ∞ and explicit here.

It is useful to consider the process whose infinitesimal generator is given by the Doob-transform

φ -1 ∞ G + λ 1 φ ∞ • = φ -1 ∞ G + b + λ 1 φ ∞ • , because
product-criticality is preserved under this transformation. The ground states associated to this new operator and its dual are 1, φ ∞ q ∞ . Developing, we obtain a process whose infinitesimal generator is

G + φ ′ ∞ φ ∞ g 2 ∂ x = G + u ′ 1 u 1 g 2 ∂ x ,
with no multiplicative part. In our case study, we get 1 2 x (1 -x) ∂ 2

x + (1 -2x) ∂ x adding a stabilizing drift towards 1/2 to the original neutral WF model. The associated diffusion process is positive recurrent and so its invariant measure φ ∞ q ∞ = 6u 1 v 1 = 6y (1 -y) is integrable with mass 1. It is the beta(2, 2) limit law of the Q-process (see [START_REF] Kolmogorov | Étude de l'équation de la diffusion avec croissance de la quantité de matière et son application à un problème biologique[END_REF] and the comments at the end of Section 2 relative to the neutral WF diffusion).

Remarks.

(i) At time t, let x (n) t Nt(x) n=1
denote the positions of the BD particle system. Let

u (x, t; z) = E Nt(x) n=1 z ψ x (n) t
stand for the functional generating function (|z| ≤ 1) of the measure-valued branching particle system. u (x, t; z) obeys the nonlinear (quadratic) Kolmogorov-Petrovsky-Piscounoff PDE, [START_REF] Kolmogorov | Étude de l'équation de la diffusion avec croissance de la quantité de matière et son application à un problème biologique[END_REF]:

∂ t u (x, t; z) = b * θ (x, u (x, t; z)) + G (u (x, t; z)) ; u (x, 0; z) = z ψ(x) , where θ (x, z) = E z M(x) -z = p 2 (x) z 2 + p 1 (x) z -z or θ (x, z) = 4x (1 -x) z (z -1)
is the shifted probability generating function of the branching law of M (x) . Thus, the nonlinear part reads b * θ (x, u (x, t; z)) = b (x) u (x, t; z) (u (x, t; z) -1) which is quadratic in u.

In particular, if u (x, t)

:= ∂ z u (x, t; z) z=1 = E Nt(x) n=1 ψ x (n) t , u (x, t) obeys the linear backward PDE ∂ t u (x, t) = b (x) u (x, t) + G (u (x, t)) ; u (x, 0) = ψ (x) involving G (•) = G (•) + b (x) •. We have the Feynman-Kac interpretation u (x, t) = E x e t∧ τ x 0 b( xs)ds ψ ( x t ) .
The latter evolution equation is the backward version of the forward PDE giving the evolution of p (x; t, y) as ∂ t p (x; t, y) = G * (p (x; t, y)) , p (x; 0, y) = δ x (y) .

(ii) Let us look at the branching diffusion process governed by G would time be measured using the time substitution

θ t = t 0 g 2 ( x s ) ds = t 0 x s (1 -x s ) ds
for each of the particles that came to birth before t.

Then G → G := 1 x(1-x) G = σ∂ x + 1 2 ∂ 2 x + 1 2 σ 2 • .
In particular, each motion y θ = x t θ is a Brownian motion with constant drift (a Gaussian process). This new G is the one of absorbing Brownian motion with drift σ on [0, 1] , including branching at constant rate 1 2 σ 2 . The Sturm-Liouville problem for G admits the eigenvalues

λ k = k 2 +σ 2 2
, k ≥ 1 with eigen-states u k (x) ∝ e -σx sin (kπx) and v k (y) ∝ e σy sin (kπy) . The spectral gap is λ 1 = 1+σ 2 2 > 0 and the time-changed branching diffusion also becomes eventually extinct, sub-critically: The time substitution changes the spectral structure of the model but not its qualitative features. ⋄

Doob transform of the WF model with mutations: critical BD

In this Section, we start from the WF model with mutations. Using the same Doob transform based on the additive functional α (x) = e σx to introduce selection, we end up with a WF diffusion process with killing and branching describing the effect of selection on the WF model in the presence of mutations. We show that in this setup, the resulting branching diffusion process is no longer sub-critical; rather, it turns out to be critical.

Suppose the starting point model is now the WF diffusion with mutations:

dx t = (π 1 -πx t ) dt + x t (1 -x t )dw t , x 0 = x ∈ (0, 1) , with π := π 1 + π 2 . For this model, G = (π 1 -πx) ∂ x + 1 2 x (1 -x) ∂ 2
x and both boundaries are chosen as being entrance (reflecting) 3 . The WF diffusion process with mutations is now ergodic. With

λ k = k (k -1 + π) 2 , k ≥ 0,
its transition density p (x; t, y) now admits the discrete spectral representation

(33) p (x; t, y) = k≥0 b k e -λ k t u k (x) v k (y) .
Here, u k (x) are the Jacobi polynomials rescaled on [0, 1] and normalized to have value 1 at x = 0. In particular, u 0

(x) = 1, u 1 (x) = 1 -π π2 x, u 2 (x) = 1 -2(1+π) π2 x + (1+π)(2+π) π2(1+π2) x 2 ,... Next, v k (y) = m (y) u k (y) where m (y) = Γ (2π) Γ (2π 1 ) Γ (2π 2 ) y 2π1-1 (1 -y) 2π2-1
is the speed density of the ergodic WF diffusion mutations (its normalized invariant measure). Note that the k = 0 term in (33) is precisely m (y) as required.

Because the transition probability density of the WF diffusion with mutations has also a discrete spectral representation, this model is amenable to a similar analysis than the neutral WF diffusion.

Proceeding as for the neutral case, we shall consider the following transformation of paths for the WF model with mutations: Let α (x) = e σx and consider a transformed process with infinitesimal generator

G (•) = α -1 G (α•) . The multiplicative part of G is now λ (x) = G (α) /α = σ (π 1 -πx) + σ 2 2 x (1 -x) .
Note that now α is neither sub-harmonic nor super-harmonic for the infinitesimal generator G including mutations because the sign of λ (x) varies as x varies.

In this case study, one selects sample paths of the WF diffusion model with mutations (x t ; t ≥ 0) with large terminal values of α (y) . The dynamics of ( x t ; t ≥ 0) is easily seen to be the WF with mutation and selection dynamics of the type (20)

d x t = [(π 1 -π x t ) + σ x t (1 -x t )] dt + x t (1 -x t )dw t ,
subject to additional quadratic killing and branching at rate λ (x) inside I. We indeed have

G (•) = e -σx G (e σx •) = λ (x) • + G (•) , where G = [(π 1 -πx) + σx (1 -x)] ∂ x + 1 2 x (1 -x) ∂ 2
x is the KBE operator of the dynamics ( x t ; t ≥ 0) . 3 When both the mutation rates u 1 and u 2 are greater than 1/2, the boundaries are entrance.

When either u 1 or u 2 is smaller than 1/2 the corresponding boundary is regular and one needs to specify whether it is reflecting or absorbing or a mixture of the two. We force here the regular boundaries to be entrance.

To summarize, in our branching diffusion way to look at the action of selection, we move from the WF diffusion with mutations (x t ; t ≥ 0) to the standard WF diffusion with mutation and selection ( x t ; t ≥ 0) , but subject to additional killing/branching at rate λ (x) . Note that when π > σ/2, λ * = σ (π 1 ∨ π 2 ) whereas when π < σ/2, λ * = λ (x * )∨σπ 2 where x * = 1/2 -π/σ > 0.

The density of the transformed process is p (x; t, y) = α(y) α(x) p (x; t, y) . It is exactly known because p is known from (33).

The transformed process (with infinitesimal backward generator G) accounts for a branching diffusion (BD) where a diffusing mother particle (with generator G and started at x) lives a random exponential time with constant rate λ * . When the mother particle dies, it gives birth to a spatially dependent random number M (x) of particles (with mean µ (x)). If M (x) = 0, M (x) independent daughter particles are started where their mother particle died; they move along a WF diffusion with mutation and selection (with generator G) and reproduce independently, and so on. For such a transformed process, the trade-off is of a different nature: there is a competition between the boundaries {0, 1} which are now reflecting for the system of particles and the number of particles N t (x) in the system at each time t, which may grow or diminish due either to branching or killing events. In the presence of mutations, the particles are no longer killed once they hit the boundaries, suggesting that there should be a greater amount of them alive in the system. However, in this new model, there is an opportunity to kill the particles inside the definition domain, when they branch. The question now being: does the new trade-off result in global extinction or global explosion of the particle system? We will now show that critical global extinction occurs.

The density p of the transformed process again has the interpretation (28), where p (n) (x; t, y) is the density at (t, y) of the nth alive particle descending from the ancestral one (Eve), started at x. In the latter formula, the sum vanishes if N t (x) = 0. A particle is alive at time t if it came to birth before t and has not yet been killed by a killing event.

Let ρ t (x) = (0,1) p (x; t, y) dy. Then ρ t (x) is the expected number of particle alive at time t. We have ∂ t ρ t (x) = G (ρ t (x)) , ρ 0 (x) = 1 (x ∈ (0, 1)) .

But then q (x; t, y) := p (x; t, y) /ρ t (x) obeys the forward PDE ∂ t q (x; t, y) = -∂ t ρ t (x) ρ t (x) + b (y) q (x; t, y) + G * (q (x; t, y))

as a result of ∂ t p (x; t, y) = G * (p (x; t, y)). We again have (29), with q (x; t, y) the average presence density at (t, y) of the system of particles all descending from Eve started at x.

Clearly -log ρ t (x) t → t→∞ λ 0 = 0 (and therefore also -∂tρ t (x) ρ t (x) ), because

ρ t (x) = 1 α (x) k≥0
b k e -λ k t u k (x) The BD transformed process therefore admits an integrable Yaglom limit q ∞ , solution to -G * (q ∞ ) = λ (y) q ∞ or -G * (q ∞ ) = 0. With v 0 (y) = m (y) , the first eigenvector of -G * associated to the smallest positive eigenvalue λ 0 = 0 (the equilibrium density of the WF diffusion with mutations), q ∞ is of the product form (35) q ∞ (y) = e σy m (y)

1 0 e σy m (y) dy .

This explicit limiting probability q ∞ is the Yaglom limiting average presence density at (t, y) for the BD system of particles (it is also the ground state for G * ).

There is also a natural eigenvector φ ∞ of the backward operator -G, satisfying -G φ ∞ = 0 (the ground state for G). It is explicitly here Both operators G (•) and its adjoint are again critical. The constant λ 0 = 0 is the new generalized principal eigenvalue; The eigen-functions φ ∞ , q ∞ are the new associated ground states.

c 0 (

 0 (x s ) ds -tµ (c) = α (x) . Thus α (x) := ∞ E x (c (x s )) -µ (c)) ds solves the Poisson equation -Gα (x) = c (x) := c (x) -µ (c) .

3. 1 .

 1 The neutral Wright-Fisher model. Consider a discrete-time Galton Watson branching process preserving the total number of individuals in each generation. We start with N individuals. The initial reproduction law is defined as follows: Let |k N | := N m=1 k m = N and k N := (k 1 , ..., k N ) be integers. Assume the first-generation random offspring numbers ν N := (ν N (1) , ..., ν N (N )) admit the following joint exchangeable polynomial distribution on the discrete simplex |k N | = N :

( 21 )

 21 P (x t ∈ dy | τ x > t) ∼ t→∞ dy is asymptotically uniform in the Yaglom limit. As time passes by, given absorption did not occur in the past, x t d → x ∞ (as t → ∞) which is a uniformly distributed random variable on [0, 1].

p

  (x; t, y) = E x e t∧ τ x 0 b( xs)ds | x t = y p (x; t, y) and ρ t (x) = E x e t∧ τ x 0 b( xs)ds . ⋄

  x) . ⋄ With p m (x) = P (M (x) = m), let l (x) = m≥1 p m (x) m log m = 2 log 2p 2 (x) .

Remark.

  With β = α -1 = e -σx , again G (β) = 0 and β -1 G (β•) = G (•) is the infinitesimal generator of the original WF model, now with mutations. ⋄ The birth (creating) and death (annihilating) rate λ in G is bounded from above and below on (0, 1). It may now be put into the canonical form λ (x) = λ * (µ (x) -1) where λ * = max x∈[0,1] (|λ (x)|)and (34) µ (x) = 1 + λ (x) λ * whose range belongs to the interval [0, 2] as x ∈ [0, 1] .

2 M

 2 Because µ (x) is bounded above by 2 and larger than 0, we actually get a BD with binary scission whose random offspring number satisfiesM (x) = 0 w.p. p 0 (x) = 1 -µ (x) /(x) = 1 w.p. p 1 (x) = 0 M (x) = 2 w.p. p 2 (x) = µ (x) /2. Note that λ (x) = λ * (p 2 (x) -p 0 (x)) =: b (x) -d (x)identifying the birth and death components of the full multiplicative rate λ (x).

1 0α 1 0α 1 0e

 111 (y) v k (y) dy.The expected number of particles in the system decays globally at rate λ 1 towards the non-zero limiting value ρ ∞ (x) := α (x)-1 b 0 u 0 (x) (y) v 0 (y) dy = e -σx σy m (y) dy.

e 1 0e

 1 σy m (y) dy = e -σx σy m (y) dy.

We note that we have the L 1 -product property

Clearly the ground states of -G * and -G are defined up to arbitrary multiplicative constants. Note that we chose these constants in such a way that 1 0 q ∞ (y) dy = 1 and

Because p 2 (x) is a degree two polynomial in x, we have the condition:

(37)

We conclude (following [START_REF] Asmussen | Strong limit theorems for general supercritical branching processes with applications to branching diffusions[END_REF] and [START_REF] Asmussen | Some modified branching diffusion models[END_REF]) that, as a result of the condition (37) being trivially satisfied, global extinction holds critically, in the following sense:

There exists a finite positive constant :

(iii) For all bounded measurable function ψ on I :

ψ (y) q ∞ (y) dy.

From (i), it is clear that the process gets ultimately extinct with probability 1. In the trade-off between killing-branching and reflection at the boundaries, all particles get eventually absorbed but the global BD process turns out be critical. Thus, the killing part of λ (x) is strong enough to avoid the explosion of the number of particles inside the unit interval, resulting in an overall critical process where global extinction still holds.

In the statement (ii) , 1 -P (N t (x) = 0) = P (T (x) > t) where T (x) is the global extinction time of the particle system. The Pareto tails of T (x) decay like t -1 , thus algebraically slowly: the time till extinction in this critical model is much longer than in the previous neutral sub-critical case (with exponential tails). From (ii) , φ ∞ (x) has again a natural interpretation in terms of the propensity of the particle system to survive to its extinction fate.

(iii) with ψ = 1 reads

µ giving an interpretation of the constant µ. The constant µ is also ( [START_REF] Asmussen | Some modified branching diffusion models[END_REF], page 287)

and so is explicitly available in our case.

The ground states of G + λ 0 and its adjoint are thus φ ∞ , q ∞ and explicit here. It is also useful to consider the process whose infinitesimal generator is given by the Doob-transform φ

criticality is preserved under this transformation. The ground states associated to this new operator and its dual are 1, φ ∞ q ∞ . Developing, we obtain a process whose infinitesimal generator is

with no multiplicative part. The associated diffusion process is the starting point WF diffusion with mutations, which is positive recurrent and so its invariant measure φ ∞ q ∞ = u 0 v 0 = m (y) is integrable.

Remark. The functional generating function u (x, t; z) = E

of the measure-valued branching particle system obeys now the nonlinear (quadratic) PDE:

is the shifted probability generating function of the branching law of M (x) .


Romain Millot

Philippe Négrel

Chemical weathering of granitic rock: experiments and Pb-Li isotopes tracing

In a recent study (Négrel et al. 2010, Chem. Geol. Vol. 274) focusing on the lead geochemistry and Pb-isotope ratios of groundwaters along a small (53 km 2 ) endoreic granitic catchment in India (Masheshwaram, Andhra Pradesh), we have shown that most of the lead in the groundwaters is of geogenic origin. Combining a weathering model and field observations, we were able to define a two-step weathering process that includes a control on the Pb-isotope ratios by accessory phases and by the main minerals from the granite in a second step of weathering.

In order to go further and to better characterize water/rock interactions, we performed laboratory experiments with the two main granite bedrocks from this site. The aim of the present work is to better constrain the processes of water/rock interactions both in terms of source (dissolution of different primary minerals) and extent of weathering, by measuring Pb isotope signatures in addition with Li isotope signatures.

Laboratory experiments consisted in measuring the evolution through time of major and trace elements, as well as Pb and Li isotopic compositions of a rainwater solution in equilibrium with a granite powder. Experiments were carried out at 25°C with a solution/powder mass ratio of 10 considering 15 mL of reference solution TMRAIN-95 and 1.5 g of powdered granite placed in screw-top Teflon® PFA beakers. The beakers were kept in a temperaturecontrolled oven, which temperature was maintained within 5% of target temperature over the total duration of the experiments. Aliquots of the solution (after filtration at 0.2 µm) in contact with the granite powder were periodically sampled (from weeks up to 2 years) and analyzed for lead and lithium isotopic compositions.

The results show that a radiogenic contribution of lead is observed during the experiments, in agreement with the field observations, and that the light lithium isotope (6Li) is preferentially retained during uptake of Li into secondary minerals.

The results of these experiments will be discussed in the frame of the relative proportion of granite weathering (dissolution of primary minerals) to secondary mineral formation.