Scaling limits of general population processes - Wright-Fisher and branching processes in random environment Vincent Bansaye, Maria-Emilia Caballero, Sylvie Méléard arxiv:182.2362v1 [math.pr] 7 Feb 218 February 8, 218 Abstract Our motivation comes from the large population approximation of individual based models in population dynamics and population genetics. We propose a general method to investigate scaling limits of finite dimensional population size Markov chains to diffusion with jumps. The statements of tightness, identification and convergence in law are based on the convergence of suitable characteristics of the transition of the chain and strongly exploit the structure of the population processes defined recursively as sums of independent random variables. These results allow to reduce the convergence of characteristics of semimartingales to analytically tractable functional spaces. We develop two main applications. First, we extend the classical Wright-Fisher diffusion approximation to independent and identically distributed random environment. Second, we obtain the convergence in law of generalized Galton-Watson processes with interactions and random environment to the solution of stochastic differential equations with jumps. Key words: Tightness, diffusions with jumps, characteristics, semimartingales, Galton-Watson, Wright-Fisher, random environment MSC 21: 6J27, 6J75, 6F15, 6F5, 6F1, 92D25. Contents 1 Introduction 2 2 A criterion for tightness and convergence in law 5 2.1 The functional spaces and specific functions.................... 6 2.2 General statements.................................. 7 2.3 Proofs......................................... 8 CMAP, Ecole Polytechnique, CRS, route de Saclay, 91128 Palaiseau Cedex-France; E-mail: vincent.bansaye@polytechnique.edu UAM; E-mail: mariaemica@gmail.com CMAP, Ecole Polytechnique, CRS, route de Saclay, 91128 Palaiseau Cedex-France; E-mail: sylvie.meleard@polytechnique.edu 1
3 Wright-Fisher process with selection in Lévy environment 13 3.1 The discrete model.................................. 13 3.2 Tightness and identification............................. 16 3.3 Pathwise uniqueness and convergence in law.................... 19 3.4 Example........................................ 2 4 Continuous State Branching Process with interactions in Lévy environment 21 4.1 The discrete model.................................. 21 4.2 Tightness....................................... 24 4.3 Identification..................................... 25 4.4 Explosive CSBP with interactions and random environment........... 3 4.5 Conservative CSBP with interactions and random environment......... 33 5 Perspectives and multidimensional population models 35 6 Appendix 36 6.1 General construction of a discrete random variable satisfying A2........ 36 6.2 Taylor expansion for Galton-Watson process with cooperation.......... 4 6.3 Some technical results................................ 4 6.4 Stone-Weierstrass Theorem on locally compact space............... 42 1 Introduction Our motivation in this paper is the study of the scaling limits of discrete population models. The scaling parameter is and the successive population sizes Zn : n come from individual based models. We define a Markov chain recursively defined by an expression of the form Z n+1 = F Z n j=1 L i,nz n, E n, where F is a function estimating the number of individual events and for each z, e,, L i,n z, e : i, n 1 is a family of independent identically distributed random variables and En is the environment in generation n. In particular, this class of processes includes well known processes in population dynamics and population genetics, such as Galton-Watson and Wright-Fisher processes. It also takes into account the effect of random environment and can include a lot of additional ecological forces as competition, cooperation, sexe, as will be seen in this paper. For instance, in the Galton-Watson case, we have Zn, F z = z and L i,n z, e = L. The criterion for the convergence in law of the renormalized process Z[v t] / : t with a time scale v is known from the works of Lamperti [22, 23] and Grimwall [15]. A necessary and sufficient condition can indeed be given in terms of the reproduction random variable L for, which is equivalent to the convergence of the associated random walk to a Lévy process cf. [8]. The scaled process then converges as to a Continuous State Branching Process CSBP, which can be written as the unique strong solution of a Stochastic Differential Equation SDE. The parameters of this SDE are given by the limiting characteristic triplet of L. Our aim in this paper is to extend such criteria when the discrete process Z does not enjoy 2
the branching property. We obtain sufficient conditions on the asymptotic behavior of F and L, so that a suitable scaling of Z space and time changes converges in law. We develop here two main applications, which are our original motivations for this work. We generalize the convergence of Galton-Watson processes to CSBP both by adding independent and identically distributed random environments and by taking into account size dependent reproduction random variables. We also extend the Wright-Fisher diffusion approximation to the random environment framework, when the selection advantage is randomly fluctuating and can in particular provoke jumps in the frequencies of alleles. Many works have been devoted to the approximation of Markov processes. They essentially rely on tightness argument and identification of the martingale problem, see [11, 18] for two references on the topic. Several difficulties arise when considering the general framework of renormalized Galton-Watson processes and in particular when the limit is non-continuous and explosive. The branching property allowed different elegant general proofs, using the Laplace exponent [15] to get the finite dimensional convergence, but also the relation with the convergence of random walks to Lévy processes via Lamperti s transform [22, 8]. Lamperti has also introduced a powerful transform in the stable framework, see e.g. [23] and [26] and [5]. Other time changes have been successfully used to obtain scaling limits of discrete processes, in particular for diffusion approximations, see for instance [19] for branching processes in random environment and [3] for controlled branching processes. But such time change techniques are essentially restricted to branching processes or stable processes or diffusion approximations. We are here interested in processes where the branching property fails and major events may provoke jumps in the limit. The processes may even be explosive and are, in general, non-stable. To deal with such processes, we develop a general approach relying on the characteristic triplet of semimartingales Jacod-Shiryaev [18]. The law of a Markov process is determined by its transitions and thus the law of the scaled process Z[v.]/ is determined by the family of functions x Gx H = v E HZ1 / x for H continuous and bounded. Similarly, the asymptotic behavior of Gx H as for a large enough class of functions H captures the convergence of the scaled process. The convergence as of the characteristics of semimartingales defined by /H i [v t] G Z i allows to get tightness and identification of the scaled process [18] Chapter IX, under some regularity and domination assumptions. However controlling G Ṇ H for such a family of functions is difficult in general, even for Galton-Watson processes. We prove in this paper that one just needs to get uniform convergence of the characteristics G Ṇ H for H belonging to some rich enough functional space H, which is dense in the set of regular functions null at zero for a norm equivalent to H = u Hu/1 u 2. The choice of the space H will depend on the assumptions on the model. In our applications, we exploit the independence property of the variables L i,..,. : i 1. We plug at the level of the conditional increments the characterization of the law by the Laplace exponent, which is well adapted to the sum of non-negative independent random variables. 3
In the rest of the paper, we develop two applications to famous frameworks of population genetics and dynamics based on Wright-Fisher and Galton-Watson chains which were our original motivation and show the efficiency of our method. The first application studies a Wright-Fisher model with selection in a random environment which impacts the selective advantage. The environments are assumed to be independent and identically distributed and the associated random walk converges to a Lévy process. We obtain the convergence in law of the process evaluating the frequency of alleles A, jointly with the random walks of the environment, by using the functional space H = {u, w 1 e ku lw ; k, l }. We thus derive a diffusion with jumps in random environment, which generalizes the well known Wright-Fisher diffusion with selection and takes into account small random fluctuations and punctual dramatic advantages in the selective effects. The second application focuses on generalized Galton-Watson processes with reproduction law both density dependent and environment dependent. We obtain a result of convergence in law to the so called CSBP with interactions in Lévy environment introduced in [28, 16]. For such processes, characteristics are non-bounded and the result is deduced from the convergence of the compactified processes exp Z k /. To deal with the joint laws of the latter and the environment random walks, we use the functional space H = {v, w v k exp lw : k 1, l } {v, w 1 exp lw : l 1}. The results extend the criterion for the convergence of a sequence of Galton-Watson processes and the results we know in random environment [19, 4] or with interactions [1]. They are further applied to Galton-Watson processes with cooperation and to branching processes with logistic growth in random environment. The paper is organized as follows. In Section 2, we give general results for the tightness, the identification and the convergence in law of a scaled Markov process to a diffusion with jumps in R d. The functional space H is introduced in Section 2.1. Tightness and identification results are stated in Section 2.2 by assuming the uniform convergence and boundedness of characteristics G Ṇ H for any H H. Convergence requires an additional uniqueness assumption, which is here obtained from pathwise uniqueness in the applications, using standard techniques for non-negative SDE [17, 13]. Proofs of these general statements are given in Section 2.3. In Section 3, we apply our method to a Wright-Fisher model with selection in a random environment. We obtain in a suitable scaling limit a Wright-Fisher diffusion in random environment for which we prove uniqueness of solution. In Section 4 sections 4.1, 4.2, 4.3, we apply our method to Galton-Watson processes with reproduction law both density dependent and environment dependent. Section 4.4 is devoted to explosive CSBP with interactions and random environment and an application to Galton-Watson processes with cooperatives effects is developed. Section 4.5 is dedicated to conservative CSBP with interactions and random environment and an application to Galton-Watson processes with logistic competition and small environmental fluctuations is studied in details. Section 5 develops some perspectives of applications of tightness and identification results. Indeed, we expect that the method can be 4
applied in various contexts and in particular for structured populations models with mating functions F, competition or cooperation, etc. otation For x R d, we denote the euclidian norm by x. The functional norms are denoted by.. In particular the infinite norm of a bounded function f on a set U is denoted by f U,. The sets C b U, R and C c U, R denote the spaces of continuous real functions respectively bounded and with compact support defined on U. 2 A criterion for tightness and convergence in law Let X be a Borel subset of R d and U be a closed subset of R d containing a neighborhood of. Let us introduce a scaling parameter. For any, we consider a discrete time X -valued Markov chain Xk : k satisfying for any k, LX k+1 X k = x = LF x, where for any, F x, x X denotes a measurable family of X -valued random variables such that for any x X, the random variable F x x takes values in U. The natural filtration of the process X Xk+1 X k take values in U. is denoted by F k k. ote that the increments Our aim is the characterization of the convergence in law of the sequence of processes X [v.],, where v is a given sequence of positive real numbers going to infinity when tends to infinity. It is based on the criteria for tightness and identification of semimartingales by use of characteristics given in Jacod-Shiryaev [18], which consists in studying the asymptotic behavior of G x H = v E HF x x = v E HX k+1 X k X k = x, 1 for real valued bounded measurable functions H defined on U. Hypothesis H We first assume that the family of random variables Fx,x satisfies lim sup 1B,b c =. Gx b x X, This assumption avoids to get infinite jumps in the limit. We will see in the examples developed in the next sections that this condition affects both the population and the environment dynamics. Under Assumption H, we prove that the study of 1 can be reduced to a rich enough and tractable subclass H of functions H. The choice of H depends on the models and is illustrated in the examples. 5
2.1 The functional spaces and specific functions We consider a closed subset U of R d containing a neighborhood of. As usual, we write hu = ogu resp. hu gu when hu/gu tends to resp. to 1 as u tends to. We introduce the functional space d d Cb, 2 = C2 b, U, R = H C bu, R : Hu = α i u i + β i,j u i u j + o u 2, α i, β i,j R. The functions of Cb, 2 can be decomposed in a similar way with respect to any smooth function which behaves like identity at, as stated in the next lemma. The proof comes from the uniqueness of the second order Taylor expansion in a neighborhood of. Lemma 2.1. Let f = f 1,..., f d Cb, 2 d such that f i u u u i for i = 1,..., d. For any H Cb, 2, there exists a unique decomposition of the form i=1 i,j=1 H = d d α f i Hf i + β f i,j Hf i f j + H f, i=1 i,j=1 where H f = o f 2 is a continuous and bounded function and α f i H, βf i,j H, i, j = 1 d are real coefficients and β is a symmetric matrix. Let us now introduce a specific function h, which behaves as identity in a neighborhood of zero and plays a main role in our convergence results. In the whole paper, the function h satisfies h = h 1,, h d C 2 b, d ; h i u u u i ; h i u for u i = 1,..., d. 2 ote that in general this function is not a truncation function in the sense of Jacod-Shiryaev [18] since it may not coincide with identity function in a neighborhood of. Its choice will be driven by the processes we are considering and we are giving different choices of functions h in the next sections, for instance hx = 1 exp x on [ 1, when d = 1. We introduce a functional space H which contains the coordinates of h, their square products and "generates" the continuous functions with compact support in R d. More precisely, we make the following assumptions. Hypotheses H1 1. H is a subset of C 2 b, and hi, h i h j V ecth for i, j = 1,..., d. 2. For any g C c U, R with g =, there exists a sequence g n n Cb, 2 lim n g g n,u = and h 2 g n V ecth. such that Remark 2.2. In the examples of the next sections, Assumption H1.2 is proved with the use of the locally compact version of the Stone-Weierstrass Theorem. We refer to the Appendix for a precise statement. Moreover, contrary to the "convergence determining class" of [18], the functions of H will not be null or ou 2 in a neighborhood of. 6
2.2 General statements Hypotheses H 1 There exists G x H, x X, H H such that for any H H, 1. lim G x H G x H =. 2. sup x X sup x X G x H < +. Hence H 1 implies that the map x X G x H is measurable and bounded for any H H. Theorem 2.3. Assume that the sequence X is tight in X and that H, H1 and H 1 hold. Then the sequence of processes X[v.], is tight in D[,, X. We now introduce a truncation function in the sense of [18] : h = h 1,, h d, h C b U, R d, h x = x in a neighborhood of. 3 Obviously, h i hj C2 b,. The next assumption H2 is sufficient in addition to Hypotheses H1 and H 1 to get the identification of the limiting values and their representation as solutions of a stochastic differential equation. We use notations of Hypothesis H 1 and Lemma 2.1. Hypotheses H2 1. For any H H, the map x X G x H is continuous and extendable by continuity to X. 2. For any x X and any H H, G x H = d i=1 α h i Hb i x + d i,j=1 β h i,j Hc i,jx + V H h Kx, vµdv, 4 where α h i, β h i,j and Hh have been defined in Lemma 2.1, b i and σ i,j are measurable functions defined on X, V is a Polish space, µ is a σ-finite positive measure on V, K is a function from X V with values in U, V 1 K., v 2 µdv < + and c i,j x = d k=1 σ i,k xσ j,k x + h i h j Kx, vµdv. V The elements b, σ, V, µ, K will be specified in the applications. The proof of the next theorem will explicitly give their relation with the associated stochastic differential equation 5. 7
Theorem 2.4. Under the assumptions of Theorem 2.3 and the additional hypothesis H2, any limiting value of X[v.], is a semimartingale solution of the well-defined stochastic differential system t t t X t = X + bx s ds + σx s db s + h KX s, vñds, dv V t + Id h KX s, vds, dv, 5 V where X X, B is a d-dimensional Brownian motion, is a Poisson Point measure on R + V with intensity dsµdv, X, B and are independent and Ñ is the compensated martingale measure of, Id denotes the identity function. Finally the next additional assumption allows to obtain the convergence in law of the sequence of processes X[v.], in D[,, X. Hypothesis H3 The law of the initial condition X X being given, the uniqueness in law of the solution of 5 holds in D[,, X. Let us note that most often and at least in our applications, the uniqueness in law of 5 is deduced from the pathwise uniqueness. We are now in position to state the convergence result. Theorem 2.5. Assume that the sequence X converges in law in X to X and that H, H1, H 1, H2 and H3 hold. Then the sequence of processes X [v.] converges in law in D[,, X to the solution of 5. 2.3 Proofs From now on, we assume that hypotheses H, H1 and H 1 hold. We write U = U \{}. In the proofs, we use the space R b of continuous and bounded functions which are small enough close to : R b = {H C b U, R, Hu = o u 2 }. Using Lemma 2.1 and H1.1, we have We work with the norm defined for H R b such that sup u U of G x C 2 b, = V ecth + R b. 6 H h = sup u U Hu hu 2, Hu hu 2 for all x X and 1 imply that < +. In that case, the positivity and linearity G x H G x h 2 H h α h 2 H h, 7 where α h 2 = sup,x X G x h 2 < by H 1 since h 2 V ecth by H1.1. 8
2.3.1 Proof of Theorem 2.3 We first extend the assumptions H 1 to Cb, 2 in order to prove the tightness. We note that H 1.1 and H 1.2 extend immediately to H V ecth by linearity of H Gx H for any x X and 1. Lemma 2.6. For any x X, there exists a linear extension H C 2 b, G xh of G x such that H 1.1 and H 1.2 hold for any H C 2 b,. As a consequence, writing αh = sup,x X G x H, for any H C 2 b,, sup G x H sup Gx H = αh < +. 8 x X,x X Proof. Using 6 and linearity, we only have to prove the extension to R b. Let us first prove the result for the compact supported functions of R b. Let us consider H R b with compact support and show that the sequence Gx H converges when tends to infinity. The function H/ h 2 defined on U can be extended to a continuous function g on U with compact support and g =. Then by H1.2, there exists a sequence g n n of functions of Cb, 2 uniformly converging to g and such that H n = h 2 g n V ecth. Since H = h 2 g, H n H h when n. Moreover the sequence Gx H n converges to G xh n when tends to infinity for any fixed n and uniformly in x. Let us now consider two integers m and n. Using 7, we have sup G x H m Gx H n α h 2 H m H n h,x and letting go to infinity, we obtain that G x h 2 g n n is a Cauchy sequence. converges to a limit denoted by G x H, which satisfies sup x G x H <. Moreover Then it G x H G x H G x H G x H n + G x H n G x H n + G x H n G x H. Using that G x H G x H n α h 2 g g n and choosing n and then large enough ensure that the left hand side goes to : sup Gx H G x H. x X Let us now consider H R b. We introduce a non-decreasing sequence ϕ n n C 2 R d such that ϕ n 1, ϕ n x = 1 on B, n and ϕ n x = on B, n + 1 c. For any integer m n and x X and 1, G x Hϕ m G x Hϕ n H G x 1 B,n c H C n, where C n as n infinity by H. Letting tend to infinity, we obtain that for any x X, the sequence G x Hϕ n n is Cauchy and converges to some real number G x H. Moreover G x H G x Hϕ n C n H. It follows that for any H R b, G x H G x H G x H G x Hϕ n + G x Hϕ n G x Hϕ n + G x Hϕ n G x H 2C n H + G x Hϕ n G x Hϕ n Adding that Hϕ n R b and has compact support, G x Hϕ n G x Hϕ n goes to and then G x H G x H tends to as tends to infinity uniformly for x X. It proves H 1.1 and H 1.2 and 8. 9
We now prove that a σ-finite measure can be associated to G x with measurability in the x-dependence. It describes the jumps of the limiting process. Lemma 2.7. There exists a measurable family of σ-finite measures µ x : x X on U such that for any H R b, G x H = Huµ x du. 9 U For any x X, G x is then extended by 9 to any measurable and bounded function H on R d such that Hu = o u 2. Moreover lim sup G x 1 B,b c =. 1 b x X Proof. For any x R d and H C c U, R, the map H G x H is a positive linear operator. Adding that U is locally compact, Riesz Theorem leads to the existence of a σ-finite measure µ x on U such that for any H C c U, R, G x H = U Huµ x du. The extension of this identity to any H R b follows again from an approximation procedure, using ϕ n defined in the proof of Lemma 2.6. Indeed, on the one hand monotone convergence ensures that U Hϕ n µ x goes to U Hµ x. On the other hand, G x Hϕ n G x H C n H goes to. Finally 1 comes from H with a monotone approximation of 1 B,b c by elements of R b and the convergence of G to G. Let us now turn to the proof of Theorem 2.3. We first prove the convergence of conditional increments functionals, defined for any function H Cb,2 and t > by φ t H = [v t] k=1 E HX k [v X k 1 F k 1 t] 1 = G 11 v X k 1H, k=1 where the last identity comes from the Markov property. Proposition 2.8. For any function H Cb,2 and t >, lim sup φ t H t T t G X ds = [v s]h a.s. Proof. Using 8, we have 1 v Then [v t] k=1 G X k 1H = sup t T t φ t H t G ds G X[v s]h H ds = [v t] X[v s] v t t αh G ds+o. X[v s]h αh G X T sup G [v s]hds x H G x H + O x X and the conclusion follows from H1.1 proved in Lemma 2.6. v v 1
Let us now define on the canonical space D[,, X a triplet which characterizes the limiting values of the sequence X[v.],. Using the measurability and boundedness of x G x H for x X and H Cb,2 and the truncation function h introduced in 3, we define for any ω = ω s, s D[,, X the functionals B t ω = C ij t G ωs h 1,, G ωs h d ds, t t ω = G ωs h i h j ds, ν t ω, H = t G ωs H1 U ds = t U Huµ ωs du ds 12 for any H C b R d, R such that Hu = o u 2, where the last identity comes from 9. Following Chapters II. 2 & 3 in [18] adapted to the arrival space X instead of R d, the characteristic triplet associated with the semimartingale X is given for i, j {1,..., d} by Bt = Eh Uk F k 1 = φ t h 1,, φ t h D C,ij t = φ t H = k [v t] k [v t] k [v t] Eh i U k hj U k F k 1 Ehi U k F k 1 Ehj U k F k 1 EHUk F k 1, 13 where Uk = X k X k 1 and H is a continuous bounded function on Rd null in a neighborhood of. Proposition 2.8 implies the convergence of the characteristics, as stated below. Proposition 2.9. For any T > and any i, j = 1,, d and H a continuous bounded function null in a neighborhood of, we have the following almost-sure convergences : sup B,i t Bt i X[v.] ; 14 t T,ij ij sup C t C t X[v.] ; 15 t T φ t H ν t X[v.], H. 16 sup t T Proof. From Proposition 2.8, we obtain immediately the first and last convergences and φ t h i h j ij C t X[v.] a.s. sup t T So it remains to replace φ t h i hj,ij by C t. We have Eh k U k F k 1 1 v sup G x h k αh k. 1,x X v and αh k,ij < from 8. Hence the second term in C t vanishes, which yields the result. 11
To apply Theorem 3.9 IX p543 in [18] and get the tightness, we need to check the required strong majoration hypothesis and condition on big jumps. For any process X t, t with finite variation on any finite time interval, we denote by V arx t the total variation of X on [, t]. First, using that for any H C b,2, G.H is bounded ensures that there exists a positive constant A such that for any ω D[,, X, d d V arb i ω t + i=1 i,j=1 C ij t ω + ν tw, u u 2 1 A t. 17 Second, we can control the big jumps, using that ν t., 1 B,b c t G. 1 B,b c, which tends to as b tends to infinity from 1. We obtain lim sup b w D[,,X ν t w, 1 B,b c =. 18 The proof of Theorem 2.3 of X[v.], follows from 14-18 and the tightness of the initial condition, by an application of [18, Theorem 3.9, chapter IX, p543]. 2.3.2 Proofs of Theorems 2.4 and 2.5 Let us now assume the additional Hypothesis H2. We wish to identify the limiting values of X[v.], as solutions of the stochastic differential system 5. We first need to extend continuously the limiting characteristic triplet to the boundary. Lemma 2.1. i For any H R b, the map x X G x H is continuous and extendable by continuity to X. Moreover sup G x H αh < +. 19 x X ii For any H R b and x X, G x H = V HKx, vµdv 2 and in particular V 1 K., v 2 µdv is bounded on X. Proof. Let H R b. Using the sequence ϕ n of the proof of Lemma 2.6 and approximating ϕ n H for. h by H n V ecth R b as in the proof of Lemma 2.7, we obtain that G x H Gx H n H C n + ϕ n H H n h α h 2 sup x X, 1 tends to as n. Letting ensures that G. H n converges uniformly to G. H as n. Using H2.1 for H n, we obtain that G. H is continuous on X and extendable by continuity to X, which proves i. Moreover 8 yields 19 by continuity. For ii, we first consider H V ecth R b. Then α h H = β h H =, H h = H and H2.2 ensures that 2 holds for H. Let us now extend this identity to H R b with compact support. We note that H = h 2 g with g C c U, R and using Assumption H1.2, 12
we know that g is uniformly approximated by a sequence g n such that h 2 g n V ecth R b. Identity 4 implies that 2 holds for any h 2 g n and G x h 2 g n = h 2 g n Kx, vµdv. V We make n tend to infinity in both terms using 8 and the assumption V 1 Kx, v 2 µdv < +. The extension to R b follows again from a monotone approximation by the compact supported functions Hϕ n, which ends the proof. This lemma allows us to extend the definitions of the characteristics and the identity of 12 to any w D[,, X. Moreover i ensures that w D[,, X B t ω, C t ω, ν t ω, H is continuous and that the dominations 17 and 18 extend from X to X. We can then apply [18, Theorem 2.11, chapter IX, p53] on the closed set X for the identification. We obtain that any limiting value of the law of X [v.] is a solution of the martingale problem on the canonical space D[,, X with characteristic triplet B, C, ν, where C ij t = C ij t ν t., h i h j. Finally, using H2.2 for H {h i, hi hj } and 2, the characteristics in 12 write t B t w = bw s ds t D C ij t w = σ i,k w s σ j,k w s ds ν t w, H = t V k=1 HKw s, vµdvds, for any w D[,, X. By [18, Chapter III, Theorem 2.26 p157], the set of solutions of the martingale problem with characteristic triplet B, C, ν then coincides with the set of weak solutions of the stochastic differential equation 5. That proves Theorem 2.4. To conclude the proof of the convergence, we remark that uniqueness hypothesis H3 guarantees iii in [18, Theorem 3.21, chapter IX, p546]. The other points i vi of this theorem have been checked above and Theorem 2.5 follows. 3 Wright-Fisher process with selection in Lévy environment 3.1 The discrete model Let us consider the framework of the Wright-Fisher model: at each generation, the alleles of a fixed size population are sampled from the previous generation. We consider a population of individuals characterized by some allele which is either A or a. The number of individuals carrying allele A is a process Zk, k whose dynamics depends on the environment. When is fixed, we consider the coupled process describing both the discrete time dynamics 13
of the population process and the environment process. For each 1, it is recursively defined for k by Zk+1 = i=1 E k,i Z k /, E k, 21 Sk+1 = S k + E k, and S =, Z = [Z ], Z [, 1] is a finite random variable, E k k are independent and identically distributed with values in 1, + and the family of random variables E k,i z, w, z, w [, 1] 1, ; k 1, i 1 are independent and for each z, w [, 1] 1,. Moreover the random variables E k,i z, w; k 1, i 1 are identically distributed as a Bernoulli random variable E z, w defined by PE z, w = 1 = pz, w ; PE z, w = = 1 pz, w. We assume that Z, Ek,i z, w, z, w [, 1] 1, ; k 1, i 1 and Ek independent., k are We assume now that the random walk S [.] converges in law to a Lévy process with characteristics α E, β E, ν E. Following [18, chap.vii Corollary 3.6] p.415, it can be stated as follows, by considering a truncation function h E, i.e. h E defined on 1, + continuous and bounded and h E w = w in a neighborhood of. Assumption A. There exist α E R, σ E and a measure ν E on 1, + satisfying 1,+ w2 1ν E dw < + such that lim Eh EE = α E ; lim Eh2 EE = β E ; lim EfE = fwν E dw, 1, for any f null in a neighborhood of, continuous and bounded and β E = σ 2 E+ 1, h2 Ewν E dw. The small fluctuations of the environment are given by σ E, while the dramatic events are given by the jump measure ν E. egative jumps will correspond to dramatic disadvantages of allele A and an usual set of selective advantages is 1,, as illustrated in Section 3.4. The limiting environment process Y can thus be defined by Y t = α E t + t σ E db E s + t 1,+ h E wñ E ds, dw + t 1,+ w h E w E ds, dw,22 where B E is a Brownian motion and E is a Poisson point measure on R + 1, + independent of B E with intensity measure ν E. By construction, this Lévy process has jumps larger than 1. Let us first prove a consequence of Assumption A which will be needed in the proof of the next theorem. 14
Lemma 3.1. Let g C 3 [, 1] 1,, R bounded and satisfying gz, = for any z [, 1]. Then, under Assumption A, uniformly for z [, 1], with Egz, E B z g, g B z g = α E w z, + β E 2 g z, + ĝz, wν 2 w2 E dw 1, and ĝz, w = gz, w h E w g w z, h Ew 2 2 g 2 z,. w 2 Proof. Indeed, we can decompose as follows Egz, E = g w z, Eh EE + 1 2 g 2 w 2 z, Eh EE 2 + Eĝz, E. The two first terms converge uniformly as by a direct application of Assumption A. Moreover the last part of Assumption A can be extended to continuous function fw = ow 2 using a monotone approximation of f by functions null in a neighborhood of. Then the last term converges for fixed z and it remains to prove that the convergence is uniform on [, 1]. First, let us consider a compact set K = [, 1] [ 1 + ε, A] of [, 1] 1,. As g is C 3 [, 1] 1,, R, the function z, w ĝz, w h E w 2 = gz, w h Ew g w z, h E w 2 1 2 g z, 2 w2 and its first derivative with respect to z are well defined on [, 1] 1, [, 1] {} and extendable by continuity to [, 1] 1,. Thus the derivative of ĝz, w/h E w 2 with respect to z is bounded on K and using that Eh E E 2 is bounded by the second part of Assumption A, there exists C > such that for any 1, ĝz, E E 1 E [ 1+ε,A] E ĝz, E 1 E [ 1+ε,A] C z z. Moreover, using that all functions involved in the definition of ĝ are bounded, there exists C > such that ĝz, E E 1 E [ 1+ε,A] C PE [ 1 + ε, A] and by the last part of Assumption A, lim sup PE [ 1 + ε, A] = lim ν E 1, 1 + ε A, =. ε,a ε,a Combining the two last inequalities, we obtain that the family of functions Eĝ., E is uniformly equicontinuous on [, 1] and the convergence is uniform by Ascoli Theorem. 15
We now assume that p is a C 3 -function from [, 1] 1, to [, 1] and pz, = z for any z [, 1]. A main example is given by pz, w = z1 + w z1 + w + 1 z 23 and describes a case where the environment w acts as a selection factor on the probability to get the allele A. By construction, this selection factor is larger than 1. The particular case when the environment is non-random, i.e. Ek = s/ a.s. for some real number s 1, +, yields the classical Wright-Fisher process with weak selection. It is well known that in this case, the processes Z[.] converge in law to the Wright-Fisher diffusion with selection coefficient s whose equation is given by dz t = Z t 1 Z t db t + sz t 1 Z t dt. Here we generalize this result for a random environment. 3.2 Tightness and identification For the statement, we introduce the drift coefficient inherited from the fluctuations of the environment: p b 1 z = α E w z, + σ E 2 p z, + pz, w z h 2 w2 E w p z, ν E dw. 1, w Theorem 3.2. Under Assumption A, the sequence of processes Z[.] /, S [.] is tight in D[,, [, 1] R and any limiting value of this sequence is solution of the following well defined stochastic differential equation Z t = Z + + t 1, b 1 Z s ds + Y t = α E t + σ E B E t + t Zs 1 Z s db D s + σ E t pz t, w Z t Ñdt, dw; t 1, t h wñdt, dw + E p w Z s, db E s + 1, w h E wdt, dw, 24 where B D and B E are Brownian motions; is a point Poisson measure on R + 1, with intensity dtν E dw and Ñ is the compensated martingale measure of ; Z, B D, B E and are independent. Proof. We apply our results to the Markov chain X k = Z k, S k, k which takes values in X = [, 1] R. Let x = z, y X, we set Fx = Fz,y = 1 i=1 E iz, E, y + E and we have Fx 1 x = E i z, E z, E. 25 i=1 16
The state space of the random variables F x x is U = [ 1, 1] 1, +. We first prove that the assumptions H, H1, H 1 and H2 are satisfied with v =. i Let us first check H. We take b > and consider G x 1 B,b c = E1 B,b cf x x. Then E1 B,b cfx x 1 P E i z, E z > b/ 2 + P E > b/ 2 i=1 We observe that 1 i=1 E iz, E z [ 1, 1] a.s. Moreover the last part of Assumption A ensures that lim sup P E > b/ 2 ν[b/ 2 1,, which tends to as b +. Then sup,x [,1] 1, G x 1 B,b c and H is satisfied. ii We define the function h on U by hu, w = 1 e u, 1 e w and the space H is the subset of real functions on U defined as H = {u, w U H k,l u, w ; k, l }, with H k,l u, w = 1 e ku lw. We can apply the local Stone-Weierstrass Theorem on the locally compact Hausdorff space U {, } see Appendix 6.4 to the subalgebra V ecth, which is thus dense in the set of functions of C c U null in zero. Moreover V ecth is stable by multiplication by h 2. We deduce that H1.2 is satisfied, while H1.1 is obvious. Let us now prove that H 1 is satisfied. We need to study the limit of Gx H k,l as tends to infinity. Recall that Gx H k,l = EH k,l Fx x with x = z, y and Fx x given by 25. We have Gx H k,l = E 1 e k i=1 E iz,e z e le = 1 E E [e k Ez,E z E ] e le [ e = 1 E e k 1 z pz, E + e k z 1 pz, E ] le. Using the following Taylor expansion log e k 1 z p + e k z 1 p = k k2 z p + 2 2 p1 p + O1/ 3, with 3 O1/ 3 bounded uniformly in p, z [, 1], we obtain Gx H k,l = E 1 e kz pz,e le.e k2 2 pz,e 1 pz,e.e O1/ 2 = E 1 [1 A k,l 1 + B k, 1 + R k, ] z, E, 17
where 2 R k, z, w is uniformly bounded for z [, 1], w 1, and 1 and A k,l z, w = 1 exp kpz, w z lw ; By expansion, we obtain G x H k,l = E A k,l z, E 1 + O 1/ B k, z, w = k2 2 pz, w1 pz, w+o 1 2. k2 2 E pz, E 1 pz, E + O1/. 26 Using Lemma 3.1 both for z, w A k,l z, w and z, w pz, w1 pz, w z1 z, we obtain from 26 that Gx H k,l G x H k,l = B z A k,l k2 z1 z, 2 uniformly for x = z, y [, 1] R. Then H 1.1 is satisfied and for any x X, A k,l G x H k,l = α E w z, + β E 2 with 2 A k,l z, + w2 1,  k,l z, wν E dw k2 z1 z 27 2 and A k,l w p z, = k z, + l, w 2 A k,l w 2 z, = k 2 p w 2 z, k p w z, + l 2  k,l z, w = A k,l z, w h E w A k,l w z, h2 Ew 2 2 A k,l z,. w2 The assumptions on the function p allow us to conclude that H 1.2 is also satisfied. iii We now check H2. First, the continuity of G. on [, 1] R comes from the regularity of p and from the integrability assumption on ν E by Assumption A with use of Lebesgue s Theorem. ow, let us introduce the truncation function defined on U = [ 1, 1] 1, by u, w U, h u, w = u, h E w. With the notation of Section 2, recall that h 1 u, w = u and h2 u, w = h Ew. With the notation of Lemma 2.1, we have α h 1 H k,l = k, α h 2 H k,l = l, β h 11 H k,l = k2 2, βh 12 H k,l = β h 21 H k,l = kl 2, βh 22 H k,l = l2 2. Moreover, setting Kz, w = pz, w z, w, we note that A k,l = H k,l K and  k,l z, w = H k,l h Kx, w + kp 1 z, w k2 2 p 2z, w klh E wqz, w, 18
with p 1 z, w = pz, w z h E w p w z, h Ew 2 2 p 2 z,, w 2 2 p 2 z, w = pz, w z 2 h 2 p Ew w z, and qz, w = pz, w z he w p w z,. ow we set V = 1, and choose µ = ν E and for x = z, y [, 1] R, we define p b 1 x = b 1 z = α E w z, + β E 2 p z, + p 2 w2 1 z, wν E dw ; b 2 x = α E V σ 11 x = p z1 z ; σ 22 x = σ E ; σ 21 x = ; σ 12 x = σ E z,. w Then 27 rewrites G x H k,l = kb 1 x + lb 2 x k2 2 c 11x l2 2 c 22x klc 12 x + where, recalling that β E = σ 2 E + 1, h2 Ewν E dw, V H k,l h Kz, wν E dw, p 2 c 11 x = z1 z + β E z, + p 2 z, wν E dw w V = σ11x 2 + σ12x 2 + h 1 Kz, w 2 µdw, V c 22 x = β E = σ 22 x 2 + h 2 Kz, w 2 µdw, V p c 12 x = β E z, + h E wqz, wν E dw = σ 12 xσ 22 x + h 1 w h 2 Kz, wµdw. 1, V Thus H2 holds for any H = H k,l H. We can now apply Theorems 2.3 and 2.4 for tightness and identification and conclude. 3.3 Pathwise uniqueness and convergence in law To get uniqueness for 24 and convergence of the scaled process, we use here pathwise uniqueness results and need some additional regularity of the terms of 24. We follow Li-Pu [25] and assume that Assumption B. For any w 1, +, the function z pz, w is non-decreasing. Moreover there exist functions r, b r and b d such that for any z [, 1], b 1 z = b r z + b d z, with b d non-increasing and r non-negative, non-decreasing and concave,. 1/rzdz = +, and b r z b r z r z z for any z, z [, 1]. The monotonicity assumption on p is natural since pz, w denotes the probability to get an offspring carrying allele A in a population where z individuals carry allele A. The regularity of b 1 may impose some moment conditions on the measure ν E, as seen for an example in Section 3.4 below. Corollary 3.3. Under Assumptions A and B, the sequence of processes Z[.] /, S [.] converges in law in D[,, [, 1] R to the unique strong solution Z, Y of 24. 19
Proof. We show that pathwise uniqueness holds for 24, which ensures H3 and allows us to apply Theorem 2.5. The pathwise uniqueness of the process Y is well known. Let us focus on the first equation of 24 and prove the pathwise uniqueness of the process Z. We are in the conditions of applications of Theorem 3.2 in [25]. Indeed, we remark that the Brownian part of 24 writes p p 2 Zt 1 Z t dbt D + σ E w Z t, dbt E = Z t 1 Z t + σ 2 E w Z t, dw t = σz t dw t, with W Brownian motion since B D and B E are two independent Brownian motions. We prove easily that for any z 1, z 2 [, 1], σz 1 σz 2 2 L z 1 z 2 for some constant L >. Using ν E V V < with V = 1/2, 1 and z [, 1] pz, w z/w is uniformly Lipschitz for w V since its first derivative is bounded, all the required assumptions for [25] Theorem 3.2 are satisfied and we get the pathwise uniqueness of the solution of 24. 3.4 Example A main example is given by pz, w = z1 + w z1 + w + 1 z 28 and describes a case where the environment w acts as a selection factor on the probability to get the allele A. By construction, this selection factor is larger than 1. The particular case when the environment is non-random, i.e. Ek = s/ a.s. for some real number s 1, +, yields the classical Wright-Fisher process with weak selection. It is well known that in this case, the processes Z[.] converge in law to the Wright-Fisher diffusion with selection coefficient s whose equation is given by dz t = Z t 1 Z t db t + sz t 1 Z t dt. Here we generalize this result for a random environment. ote that in this case, we have Then b 1 z = α E z1 z σ E z 2 1 z + p w z, = z1 z ; 2 p w 2 z, = 2z2 1 z. 1, wz1 z zw + 1 h Ewz1 z ν E dw. 29 Under Assumption A and by Theorem 3.2, the sequence of processes Z[.] /, S [.] tight in D[,, [, 1] R and the first coordinate of any limiting value is solution of Z t = Z + t b 1 Z s ds + + t t t Zs 1 Z s dbs D + σ E Z s 1 Z s dbs E 1,+ wz s 1 Z s 1 + wz s Ñds, dw. 3 is 2
In particular if σ E = and ν E =, we recover the classical Wright-Fisher diffusion with deterministic selection α E. Our extension allows us to consider small random fluctuations asymptotically Brownian and punctual dramatic advantage of the selective effects. By Corollary 3.3, a sufficient condition for the convergence in law of Z[.] /, S [.] in D[,, [, 1] R to the unique strong solution of 3 is Indeed, writing kw, z = z +. w ν E dw < +. wz1 z zw + 1 h Ewz1 z = 1 + w 1 + wz 2 1 h Ew1 2z, the function z 1+w/1+wz 2 1 is decreasing for small w, kw, z = Ow 2 on [ 1/2, 1] and sup z [,1],w [1, kw, z Cw for some C >. Therefore, the function b 1 defined in 29 satisfies Assumption B. 4 Continuous State Branching Process with interactions in Lévy environment In this section, we are interested in large population approximations of population dynamics with random environment and interactions. We generalize in different directions the classical convergence of Galton-Watson processes to Continous State Branching processes CSBP, see for example [15, 22, 8]. We focus on models where the environment and the interaction mainly affect the mean of the reproduction law and thus modify the drift term of the CSBP by addition of stochastic and nonlinear terms. This allows us to construct a discrete population model approximating the so-called CSBP in Lévy environment with interactions BPLEI for large populations. The CSBPs in random environment or with interaction have recently been subject of a large attention. More precisely, we refer to [28] for existence of the solution of the associated SDE with general assumptions, [4, 3, 26] for approximations and study of some classes of CSBP in random environment without interactions, [2, 1, 29] for CSBP with interactions without random environment, while [3] considered diffusion approximations in the continuous context. The method we develop here, relying on Section 2, allows to obtain new statements both for convergence of discrete population models and for existence of solutions of SDE with jumps, as detailed in the following Theorems. 4.1 The discrete model Let us now describe our framework. We consider scaling the population size and a sequence v which tends to infinity with giving the time scale. As in Section 3, we introduce for any a sequence of independent identically distributed real-valued random variables E k k with same law as E. The asymptotical behavior of E k k is similar as the one in Section 3 Assumption A but the scaling parameter is v. As in the previous section, h E denotes a truncation function defined on 1, +. 21
Assumption A1. Let us consider α E R, σ E [, and ν E a measure on 1, such that 1 w 2 ν E dw <. 31 1, Writing β E = σ 2 E + 1, h2 Ewν E dw, we assume that lim v Eh E E = α E ; for any f nul in a neighborhood of zero. lim v Eh 2 EE = β E ; We also consider the associated random walk defined by lim v EfE = S =, S k+1 = S k + E k k. 1, fw ν E dw, We recall as in Section 3 that A1 is equivalent to the convergence of the random walk S [v.] to the Lévy process Y with characteristics α E, β E, ν E defined in 22. We reduce the set of jumps to 1, to avoid degenerated cases when a catastrophe below 1 could kill all the population in one generation. Let us fix. We assume that each individual reproduces independently at generation k with the same reproduction law which depends both on the underlying environment and on the population size. More precisely, given a population size n and an environment w, this law is the law of a random variable L n, w. We thus introduce random variables Z and L i,k n, w such that the family of random variables Z, L i,k n, w, n, w 1, +, E j ; i, k, j is independent and for each n, w 1, +, the random variables L i,k n, w are all distributed as L n, w for i, k. Moreover, we assume that the function L i,k defined on Ω 1, + endowed by the product σ-field is measurable. The population size Z k at generation k is recursively defined as follows, Z k Z = [Z ], Zk+1 = L i,k Z k, E k k, 32 i=1 with Z [, a fixed integer valued random variable. We investigate the convergence in law of the process 1 Z [v t], S [v t], t [,. We cannot apply directly our general result to Z[v.]. Indeed, the associated characteristics of the first component are not bounded. Moreover, scaling limits of Z can lead to explosive processes, as already happens in the Galton-Watson case. Therefore, we first study the convergence of the process Xk = exp Zk /, S k k 33 in DR +, [, 1] R where the state space of the first coordinate has been compactified. Following notation of Section 2, we introduce for x = exp z, y, 1] R the quantity F x = exp 1 [z] i=1 L i [z], E 1 z, y + E, 34 22
and observe that for any z /, conditionally on Xk Xk+1 is distributed as F x. = exp z, y, the random variable We apply now the theoretical framework developed in Section 2. Let us introduce χ =, 1] R and U = [ 1, 1] 1, and for u = v, w U, we define hu = hv, w = v, 1 exp w, h u = h v, w = v, h E w 35 respectively as the specific function and the truncation function. We choose the functional space H defined by H = {H k,l : k 1, l } {H l : l 1}, where for any u = v, w U, H k,l u = v k exp lw and H l u = 1 exp lw. The fact that H satisfies H1 is a consequence of the local Stone Weierstrass Theorem on [ 1, 1] [ 1, \ {, } cf. Appendix 6.4. For any k, l and x = exp z, y, 1] R, we have Gx H k,l = E H k,l exp 1 [z] L i [z], E 1 z exp z, E i=1 = e kz v E e 1 [z] k i=1 L i [z],e 1 1 e le. Using binomial expansion and setting P k z, w = E e k L [z],w 1 [z] 1 36 and A j,l z = v E Pj z, E e le. 37 we obtain by independence of the reproduction random variables conditionally on E that G x H k,l = e kz k j= for k 1, since k j= k j 1 j k =. Similarly, for l 1, k 1 k j A j,l z 38 j G x H l = v E1 exp le. Using Section 2, the asymptotic study of the sequence X [v.] can be determined via the convergence of the functions A j,l. More precisely, the uniform convergence and boundedness of exp kza j,l z ensure the tightness of X [v.] using Theorem 2.3. The continuity of the limiting functions is involved in the identification of the characteristic triplet of the limiting semimartingales. Finally, their representation as solutions of a stochastic differential equation 23
and the associated uniqueness yield the convergence Theorem 2.5. Let us note that this characterization of the convergence of the process in terms of characteristics of the joint distribution of the offspring size and the environment applies in particular to Galton-Watson processes. In that case, E =, L z, w = L and A j,l z becomes A j z = v Pj z, = v E e j [z] L 1 1. We observe that A j z v [z]ef j L 1/ as. The uniform convergence of e jz A j z is equivalent to the classical necessary and sufficient condition for convergence in law of Galton-Watson processes [15], see Section 4.4 for details. In the next sections, we generalize this convergence by considering a class of reproduction random variables taking into account interactions and random environment and converging to CSBP with interactions in Lévy environment. 4.2 Tightness Using the process X[v.] defined in 33, we first prove the tightness of the process by assuming uniform convergence of characteristics A j,l defined in 37. Assumption A1. For any 1 j k and l, there exists a bounded function A j,k,l such that e kz A j,l z A j,k,l z uniformly for z. We state a tightness criterion for the original scaled process in the state space [, ] R endowed with a distance d which makes it compact and then Polish, say dz 1, z 2 = exp z 1 exp z 2 for z 1, z 2 [, ] with the convention exp =. Theorem 4.1. Under Assumptions A1 and A1, the sequence of processes 1 Z [v t], S [v t], t [, is tight in DR +, [, ] R. Proof. We first recall that for l 1, v E 1 e le γl E 39 for 1 e lw = lh E w 1 2 l2 h 2 Ew + κw, where κw = ow 2 is continuous bounded and using Assumption A1. Then we can define G. on H l for l 1 as G x H l = γ E l, 4 24
for any x X =, 1] R. Let us now define G. for H k,l H and k 1, l. We set G x H k,l = k j= k 1 k j A j,k,l z. 41 j for x = e z, y X. Using Assumption A1 and 38, we obtain that lim G x H G x H = sup x X for any H H. Moreover G. H is bounded by A1 and Hypothesis H 1 is satisfied. The tightness of X [v.] is then a consequence of Theorem 2.3 and yields the result. 4.3 Identification We aim now at identifying the limiting values of X [v.] as diffusions with jumps. We are interested in models where the environment and the interactions affect the mean reproduction law and the limiting process will be CSBP with interactions in a Lévy environment BPLEI. We introduce a truncation function h D on the state space, +, the parameters α D R and σ D and a σ-finite measure ν D on, + such that 1 z 2 ν D dz < +. 42 We also consider a locally Lipschitz function g defined on R + such that e z z gz z. 43 The function g models the interaction between individuals. In the applications to population dynamics, the most relevant functions will be polynomial. We provide now the scaling assumption on the reproduction random variable L so that the limiting values of Z / can be identified to a BPLEI. This assumption will become more explicit and natural through the identification and examples of the next sections. Assumption A2. Setting for z, γz D = α D z 1 2 σ2 D z 2 + γz E = α E z 1 2 σ2 E z 2 + we assume that for any 1 j k and l,,+ 1,+ 1 e zv zh D v ν D dv, 1 e zw zh E w ν E dw, sup e kz A j,l z + jz gz + γd j z + γjz+l E γe l, 44 z where A j,l has been defined in 37. 25
Remark 4.2. i In Appendix 6.1, we provide an explicit construction of a family of random variables L z, e satisfying A2, in the case β D = and α D and ν D and α E, β E, ν E and g satisfying the moment and regularity conditions given above. ii We believe that the simple convergence in A2 is actually necessary for the convergence of the process Z / to a BPLEI. The uniformity can be obtained for many classes of reproduction laws via uniform continuity, using for example monotone or convex arguments, or boundedness of derivative on compact sets as in the previous section, see forthcoming examples. iii Finally, let us remark that we only need to prove the previous convergence for z / in A2, using the definition of A j,l z and the uniform continuity of the limit. It will be more convenient for examples. ote that under Assumption A2, Assumption A1 is satisfied with A j,k,l z = e kz jz gz γ D j z + γ E l γe jz+l, which is indeed bounded using 43 and the fact exp kzγ E jz+l is bounded, since γe jz+l C l,j z + z 2 + e jz/2 z 2 β E + e jz ν E 1, 1/2. Thus tightness holds. Moreover we can simplify the expression 41 of the limiting characteristic G x, which writes G x H k,l = e kz k j= for x = e z, y. For that purpose, we denote and use that k 1 k j jz gz γj D z + γl E j γe jz+l f z u = 1 e zu k k j= j 1 k j j = δ 1,k and k binomial expansion, k k j= j 1 k+1 e lu f z u k. By straightforward computation, we obtain that for k 3, 45 k j= j 1 k j j 2 = 2δ 2,k + δ 1,k and by a 1 k j f j u = 1 k+1 f 1 u k and k k j= j 1 k j f jz+l u = G x H k,l = 1 k e kz e lw f z w k ν E dw + z 1,+ and for k = 2,,+ f 1 v k ν D dv, 46 G x H 2,l = e {z 2z 2 β E + e lw f z w 2 z 2 h 2 Ew ν E dw + zβ D 1,+ + z f 2 1 v h 2 Dv } ν D dv, 47,+ and for k = 1, G x H 1,l = e z{ γ E l γe z+l zgz zγd 1 }. 48 26