HYPER-FINITE METHODS FOR MULTI-DIMENSIONAL STOCHASTIC PROCESSES By MARK ALLAN REIMERS B.Sc, The University of Toronto, 1978 •Sc., The University of British Columbia, 1983 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in THE FACULTY OF GRADUATE STUDIES Department of Mathematics The University of British Columbia We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA October 1986 ©Mark Allan Reimers, 1986 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of The University of British Columbia 1956 Main Mall Vancouver, Canada V6T 1Y3 Date DE-6(3/81) ABSTRACT In this thesis we introduce Non-Standard Methods, in particular the use of hyperfinite difference equations, to the study of space-time random processes. We obtain a new existence theorem in the spirit of Keisler (1984) for the one dimensional heat equation forced non-linearly by white noise. We obtain several new results on the sample path properties of the Critical Branching Measure Diffusion, and show that in one dimension it has a density which satisfies a non-linearly forced heat equation. We also obtain results on the dimension of the support of the Fleming-Viot Process. Edwin Perkins ACKNOWLEDGEMENT I would like to thank Ed Perkins for many helpful conversations and of my friends for putting up with my complaining. iv TABLE OF CONTENTS page ABSTRACT ii ACKNOWLEDGEMENT iiTABLE OF CONTENTS v TABLE 1 1 CHAPTER ONE - Introduction 3 1.1 Why this Thesis1.2 SPDEs 4 1.3 The Dawson Critical Measure Valued Diffusion 7 1.4 The Fleming-Viot Process 9 CHAPTER TWO - Non-Standard Analysis and Probability 10 2.1 Some Definitions and Notations from Non-Standard Analysis 10 2.2 Non Standard White Noise 14 2.3 Adapted Stochastic Integrals 6 CHAPTER THREE - The Heat Equation with Non-Linear Stochastic Forcing 20 3.1 Scope 23.2 White Noise on the space U 23 3.3 Hyper-Finite Difference Equations 24 3.4 Some Useful Inequalities 26 3.5 Bounds on Moments of U 7 tx 3.6 Bounds on Moments of- Spatial Differences 31 3.7 Bounds on Moments of Temporal Differences 33 3.8 S-Continuity and the Standard Part 35 3.9 Solution of the SPDE 38 CHAPTER FOUR - The Dawson Critical Branching Diffusion 41 4.1 Introduction 4V 4.2 A Hyperfinite Difference Equation 43 4.3 The Coefficients Q 45 4.4 The Total Mass Process M 52 4.5 S-Continuity of the Process 5 4.6 Characterization by a Martingale Problem 62 4.7 New Results on the Dawson Measure Valued Diffusion 66 CHAPTER FIVE - The Critical Branching Diffusion in One Dimension 70 5.1 Introduction 75.2 The SPDE and the Measure Diffusion 71 CHAPTER SIX - The Support of the Fleming-Viot Process 73 6.1 Introduction and Construction 76.2 The Dimension of a Putative Support Set 75 6.3 A Useful Stochastic Differential Equation 80 6.4 Verification of Support 83 APPENDIX A — Some Inequalities Used in Chapter 3 ' 88 A.1 Purpose 8A.2 Some Identities 8A.3 Some Inequalities 91 APPENDIX B. - Internal Solutions to SPDEs in Higher Dimensions 97 REFERENCES 100 TABLE 1 of Notations Meaning Infinitesimal grid spacing in space Infinitesimal grid spacing in time At/Ax2 *-countable grid representing R1^ *-finite grid representing [0, t ] 2 Internal I.I.D. S-L random variables on T x X Borel subsets of X Coefficients for a discrete *-finite Green's formula; also density of an infinitesimal random walk Internal solution to stochastic *-finite difference equations Internal *-finite analogue of Internal *-finite analogue of A Internal solution to unforced *-finite analogue of the heat equation Lebesgue measure on R Internal measure on X , which assigns mass Axd at each grid point Space of positive finite measures on Rd !Space of probability measures on R^ Space of infinitely differentiable functions with compact support in R Space of C functions bounded on R Space of C, functions whose second b d derivatives are bounded on R Predictable square (increasing) process associated with a process x 3 CHAPTER ONE Introduction 1.1 Why this Thesis The aim of this work is to introduce nonstandard methods, and in particular the use of hyperfinite difference equations, to the theory of multi-dimensional stochastic processes. Non-standard analysis is a particularly appropriate tool when investigating random processes over a region in space, which evolve in time, especially when the support of these processes is confined to a region of infinitesimal volume, as is frequently the case. The evolution through time of all the processes discussed in this work g may be described by the heat operator A , but in the opinion of the o t author, nonstandard techniques may be equally fruitful in the analysis of processes whose development through time is described by different operators as well. 4 1.2 SPDEs The approach to the theory of SPDEs which we will follow has its home in the theory of multiparameter processes, and in particular in the theory of multi parameter stochastic integration that has been developed in recent years. Walsh (1986) contains a systematic treatment of this theory. This approach emphasizes sample path properties. An alternative approach considers SPDEs as stochastic evolutions on a space of functions, and emphasizes analytic properties. See Dawson (1975) and (1985) and the references there for further information on this approach. We will not follow it closely here. The type of SPDE we will be considering most often in this work, is Cl-1). |r- = Au + f (u) W^ , d t x where t e R+ and x e Rd . W is "white noise" on R+ x Rd ; that is, the tx derivative, in the sense of distributions (or generalized functions) of a random process W^ , indexed by sets A <= R+X.R^ •> such that 2 i) E(W^) = 0 ii) E(W^) = X(A) where X is Lebesgue measure iii). if A n B = 0 , then W^ is independent of Wg . For further information on white noise see Walsh (1986).- Chapter 1. f (u) is a real-valued function of the point values of u . Equation Cl-1). cannot possibly hold in the classical sense of an equation between the values of functions at every point in the domain. The W term is tx far too rough. Rather we usually interpret (1-1) in the weak sense. That is, OO £ if we multiply (.1-1) by a cc ( R ) (smooth, with compact support) function <j> (~x) and integrate over a rectangle CO,T] X A , where A contains supp cj> , then 5 r r (1-2) A u <j> (x) dx = Ox ( u A<J>(x)dx)ds Jo JA SX f(u )<j>(x)dW sx sx where the last integral on the right is the multiparameter stochastic integral in the sense of Ito discussed in Walsh (1986) Chapter 1. Equation (1-2) may be derived from (1-1) taken in the classical sense, by integration of the middle term by parts. An existence theory for (1-1) has been developed in the case d = 1 , when f : R -> R , is a Lipshitz - continuous function which grows at most linearly at infinity. Dawson (1972) established existence and uniqueness under these conditions, using a Hilbert space approach. Funaki (1983), established the same result, with joint continuity of sample paths in t and x . Walsh (1981) established a modulus of continuity in t and in x for solutions of an equation similar to (1-1) and investigated finer sample path properties, under the same conditions on f . In Chapter Three of this thesis an existence result is- established for (1-1) assuming only continuity and linear growth of f for d = 1 . may be regarded as a derivative of order d + 1 of a continuous function of I cl unbounded variation (the Brownian Sheet) on R x R . When d >_ 2 there is no hope of finding a continuous function u to satisfy the equation, even in the weak sense of (1-2) . The most we can hope for is to find a continuous process v , such that u may be regarded as a derivative, in the sense of The situation for d >_ 2 is entirely different. The term W tx in (1-1) 6 distributions, of order d - 1 , of v . In this case u will not, in general, have point values, and it is difficult to see what sense can be made of the term f(u ) occurring in the stochastic integral on the r.h.s. of sx (1-2). In order for our theory of stochastic integration to make sense of (1-2) we would need f(u ) to be an adapted continuous process. This is out of the question if f is supposed to be a real function of the (non-existent) point values of u . The functions f for which (1-2) can be reasonably expected to make sense are the constant functions. Walsh (1984) has shown existence and uniqueness of solutions to (1-2) in this case. We were not able to extend his results (see Appendix B). 7 1.3 The Dawson Critical Measure Valued Diffusion Measure-Valued Branching Processes (MB Processes) were first obtained by Jirina (1958) as a limit of a branching diffusion of a large number of particles. These processes were studied extensively by Watanabe (1968) and co-workers, and lately many details of the fine structure have been obtained by Dawson and Hochberg (1979). The name of Dawson is particularly associated with the case we shall study here, hence we refer to it often as the "Dawson Process". We will however most often make use of a martingale characterization of this process described in Roelly-Coppoletta (1986). A simple construction, for the case of a initial Lebesgue measure is as follows. Let particles be distributed initially on R^ or a portion thereof, according to a Poisson point process with intensity A . Suppose that thereafter each particle independently executes a Brownian motion on R . Also suppose that each particle independently undergoes critical branching with rate u , i.e. at fixed or exponentially distributed times, whose number in a unit time interval has expectation p , the particle dies or splits into two particles, each outcome being equally likely. If the particle splits, both daughter particles begin independent careers from the point of bifurcation. We now 1 d assign a mass, of' — to each particle, and obtain a random measure on R at A each time. Now suppose that both the initial density of particles A, (which is the reciprocal of the weight assigned to each particle) and the branching rate u, are allowed to go to infinity in such a way that —• is constant. Then there A d is a limiting process taking values in the space of positive measures on R Dawson (1972) indicated a connection in d = 1 between this process obtained as a limit of a particle system, and the solution to the SPDE 8 d-3) |E = + ^ „ , 9t . 2 tx 9x but this connection has not yet been made rigorous, since the function f : u —> /u is not Lipshitz, and hence (1-3) does not fall under the purview of the existence theory discussed in 1.2. With the existence theorem in chapter three, we may now close the gap in Chapter Five,' and establish the identity of the solution of (1-3) with the MB process in one dimension. As mentioned in 1.2, it is difficult to make sense of an equation like (1-3) in dimensions d >_ 2 . However a non-standard analogue of (1-3) may be constructed, and its solution may be shown to coincide with the Dawson process. This construction is used in Chapter Four to establish several previously unknown results on the sample paths of the Dawson Process. Specifically we show that i) any given Lebesgue null set is a.s. never charged, ii) the mass on any given Lebesgue set is a.s. a continuous function of time, iii). convergence of a sequence of bounded functions in measure, implies the a.s. convergence of the integrals of those functions with respect to the random measure, uniformly on finite time intervals. 9 1.4 The Fleming-Viot Process In Fleming and Viot (1979) a measure-valued process was introduced as a limit, under suitable scalings of time and space, of the Ohta-Kimura stepwise mutation model. Further results on the structure of the sample paths have been obtained by Dawson and Hochberg (1982). Among these results is the fact that for fixed times the support of the random measure has Hausdorff dimension not greater than 2, almost surely. In Chapter Six we establish this result for all times simultaneously using Non-standard methods to amplify some of the ideas of Dawson and Hochberg (1982). For further information on the Fleming-Viot process, the reader is referred to section 6.1. CHAPTER TWO Non-Standard Analysis and Probability 2.1 Some Definitions and Notations from Non-Standard Analysis For a real introduction the ideas of non-standard analysis, with a minimum of technical apparatus, we refer the reader to Cutland (1983). An internal object in the non-standard universe is one which may be referred to in the non-standard language. One of the consequences of the transfer principle is that internal objects described in the non-standard language inherit all the qualities of standard objects that are described in analoguous standard language. This is useful when dealing with hyper-finite collections, which may be treated as finite sets, though they are generally infinite. We will usually denote non-standard objects by capital letters or underlined letters. Unless otherwise noted, lower case roman letters will stand for standard objects. The embedding of a standard object into the non-standard universe will be denoted by an asterisk (*) to the left. We say that x_ e *R is infinitesimal if |x| < ^ for every n e N ; we denote this x ~ 0 ; x s y means x_ - y_ « 0 . We say x e *R is infinite, if |x| > n , for every n e N . If x e *R is finite (not infinite) then there is a unique x e R such that x_ « *x . We say that x is near standard and call x the standard part of JC , denoted °x or st(x) . These concepts may be extended, in the obvious fashion, to any space. In particular, if f e C[Rk;Rn) and F e *C(. Rk;, R") , then st(F) = f <=> ° F (x) = f(°x) for all nearstandard x e *R We shall require the.following axiom of saturation (see Cutland 1983, 1.9). 11 If {A } n neN is a decreasing family of non-empty internal sets, then n neN A is not empty. Two consequences of this are n a) Denumerable Comprehension. For every internal set A , and every function f : N -> A , there is an internal function F : *N -»- A which extends f . b) Infinitesimal Underflow. Let S be internal, S £ *R , and suppose for some a > 0 , x e S whenever 0 < °x <_ a . Then for some e & 0 , x e S whenever e <_ x_ <_ *a . A useful notion is S-continuity (S for standard). Definition 2.1.1. An internal function F : E c *Rn is called S-continuous iff x ~ y => F (x) PS F (y_) e . ns (* - Rn)' .. That this is the appropriate notion of continuity to link standard and non-standard,-.is shown by , Theorem 2.1.2. (Cutland (1983), Theorem 1.6). Let F : E <= *R *R be internal. Then °F exists and is continuous iff F is S-continuous on E . Nonstandard Probability Theory really came into its own after-the development of Loeb Measure. Theorem 2.1.3. (see Cutland (1983), Theorem 3.1). Every internal *finitely additive measure space (|~, F, y_) gives rise: to a classical o-additive measure space (|J, L(F), L (y_)) , such that F c L(F) , and if A e F , then L (A) = °y_(A) . If A e L(F) , then L(y_) (A) = inf{°y_(B) | A £ B , A e F} . Further, if L(y) (A) < 00 , then there is B e F , with L (y_) (B A A) =0 . The (complete) measure space (f, L(F), L(v_)) is called the Loeb extension of Cf, F, y_) . In the case of probability measures, we will denote the Loeb extension of (ft, F, P) by (Q, F, P) , contrary to our usual convention. 12 Theorem 2.1.4. (see Cutland (1983) Theorems 3.1, 3.5). If E is an internal field of subsets of [ , and F : [ •+ is an internal E-measurable function then the projection st ° F : [ + R u {°°} is L(E) measurable. Conversely if f : [ •-> is L(E) measurable, and [ is a-finite with respect to L(y_) (as is always the case here), then there is an internal E-measureable function F : [~ -> *R^ such that st ° F(x) = f(x) L(y) - a.e. Such an internal F is called a lifting of f . If f has finite support, i.e. L(y) ({xj f (x) ^ 0}) < 00 , then f admits a lifting F with the same property. If T £ *R^ if f : T R is continuous, then we may obtain a lifting F which is S-continuous and for which st ° F(x) = f(°x) , V 21 e r n st-1(Rd) (the near-standard points). Such an F is called a uniform lifting of f . Definition 2.1.5. An internal function F on (f, F, y) is called S-integraBle, iff F dp is finite, and F(x) dy + |F(x) I dp & 0 , for all {x: 1F(X) I < 6} {x: |F(X) | > H> infinitesimal 6 and infinite H . F is said to be S-IJ with respect to p_ , if |F|^ is S- integrable. Theorem 2.1.6. (see Cutland (1983) Theorem 3.9). If F is S -integrable on (f, F, p) , then for all A e F , F dy = F d L(y) A A Definition 2.1.7. An internal measure space (T, F, y_) is called a hyperfinite representation of a topological measure space (E, F, y) iff (i) T is a hyperfinite internal subset of *E (ii) F is the internal power set of [ (iii) A set B c E is y-measureable, iff st "*"(B) n [ is L (y) -measureable. In that case y(B) = L (y_) (st-1 (8) n D . The canonical example of a hyperfinite representation is the discrete representation of Lebesgue measure on R , which we will use frequently. Theorem 2.1.8. (see Cutland (1983) Theorem 4.1). Let Ax. be any infinitesimal, i = l,...,d. Let M. e *N\'N be infinite, l l i = 1, 2,...,d such that °M Ax. ^ 0 . Then let l I X = {k. Ax.,...,k.Ax.)Ik. e *Z, k.I < M.} . Define a measure A on X by 1 i j D ' l ' l — d setting the value of A_ on each point x_ e X to be II Ax- • Then i=l (X, P(X) , (A_)) is a hyperfinite representation of the rectangle {x| |x. | °M^Ax_^} c_ R^ , equipped with Lebesgue measure. In this thesis we will prove results about S-continuity of functions from hyperfinite grids X to *R . All the preceding definitions and * d theorems apply as if these functions were step functions on a rectangle in R Theorems about S-continuity on a hyperfinite grid X , can be translated into theorems about the weak convergence of a sequence of processes on a sequence of finite grids to a continuous limit. We will not explicitly make such a translation. 14 2.2 Non Standard White Noise Anderson (1976) introduced a hyperfinite representation of Brownian Motion, namely an infinitesimal random walk with spatial excursions of size /Kt in a time step At . Implicit in this construction was a non-standard representation of white noise on the line, as a sum of IID random variables £ , each of mean 0 and variance At , on a hyperfinite time-line of spacing At . Recently Andreas Stoll (1985) has generalized Anderson's construction to.arbitrary a-finite Radon spaces. For our purposes we only need representations of white noise on rectangles in R^ , or on all of R^ . Let X be a hyperfinite lattice as described above, and let ft be an internal space•on which are defined a family {£ } of I.I.D. S-L * X£X internal random variables, such that E(£ ) =0 and var(5 ) = 1 . For — 21 x. most of our applications we will need finiteness of all the higher moments of E, as well, x The existence of such a space ft may be shown by example. Let ft = {-1, 1} . Let F be the family of internal subsets of ft , and define A P_CA). = J—L for A e F . Then (ft, F, P) is an internal probability space, |ft| and we may define £ as the coordinate maps. We will most often be using a *countable lattice X = { (k^Ax, •. - ,k_.Ax) | kl'"*"'kd e * z^ • An exemplary space ft is then {-1, 1} .We may take F to be the *a-field generated by the *finite subsets of ft (closed under *a-unions and *a -intersections). The internal probability measure P_ may be defined as in the standard analogue: P_ is the unique *a-additive measure 15 on F such that, (<4 : x e X} are independent and for any x e X , P({co = l}) = P({w = -l}) = i . x x 2 The transfer principle guarantees that the Kolmogorov Extension Theorem carries over to the non-standard setting, and therefore that such a measure P_ exists. In this case again we will take the probability space (fi, F, P) to be the Loeb extension of (fi, F, P_) . Given such a space fi , and random variables £ , define for internal x r A ~ sets A c X , W(A) = } E, IT Ax ; the map W : A -»- W(A) is called — x . . i xeA — i=l d-dimensional S-white noise on X . Stoll (1986) shows Lemma 2.2.2. If °A_(A) < °° and _A(A A B) = 0 , then °W(A) = °W(B) P-a.s. Thus we may make Definition 2.2.3. For each Loeb measureable set A c X , with L (A_) (A) < °° , a standard random variable w(A) is well defined (up to a null set) by w(A) = °W(A) P -a.s. whenever A is internal and L(X_) (AAA) =0 . Theorem 2.2.4. (Stoll (1986) Theorem 2.5). The family (W]('A'I .|,.;L/'CA)••".<: ~}/ is a white noise on the Loeb extension of (X, X) with respect to the Loeb probability space (fi, F, P) . 16 2.3 Adapted Stochastic Integrals Let At and Ax be infinitesimals, and let T be a lattice of spacing At representing a line segment in R+ , and let X be a lattice of spacing Ax in each direction representing a rectangle in RD . Let 9, be an internal space supporting a collection E, , (t,x) e T x X , of I.I.D. tx - -2 S-L random variables, as described above. In this section we will use the notation dW to represent tx F v^At Axd . We define an internal filtration F , t e T on ft to be the tx t — algebra of internal sets generated by {W 0 < s < t, x e X} . On the Loeb sx — — space Cft, F, P) we define a filtration F , t e st(T) , by F = n a(F ) v W , where W is the collection of P-null sets. For t °t>t ± properties of this filtration see Hoover and Perkins (1983) §3. We say that an internal process utx^w^ a lifting of a process U2 I({|x| > H ; u^Cw) on R+ x Rd x fl , if for H e *NXN and n e N , t <_ n}) dX dX « 0 , and LC*. x X x P)C(Ct, x, u) | °U (u) ? u (w)}) = 0 . '-t -X — — — tX o. o - - _t, X We say that an internal process utx^ on ^ x ^ is F^-adapted, if TJ ,0<s<t,yeX,is F measureable. sy — t_ We shall call a process u^_ on R x R , F^_-adapted, if, for each t e R+ , u (jo). , (s e i-0, t], x e RD , a) e ft), is B([0,t]) x B( RD) x f sx t measureable. Since our filtration F is continuous and we do not wish to integrate discontinuous integrands, we shall not make distinctions between adapted, progressively measureable, optional, and predictable processes. We shall usually be working with F^adapted liftings of F^-adapted processes. Theorem 2.3.1. Suppose u^-x^w^ ;"'s an L F^_-adapted process on R x R x and that an internal F^ adapted process U (to) lifts u , and is — tx S-L2(T x X x n; L(A't *'A x.p)). Then for any t e R+ and t ss "t C2-1) a.s. u dw sx sx y y u dw 0<s<t x'x S5 sx where the integral on the left is in the sense of Ito. Proof: We first establish an isometry property: I I U dW 0<s<t xeX §- §-J I I I UsxdWsxV+ 2 X I I { I UsxdO 4 I Us' \o£s.<t (xeX §S s*j 0<3<£ s<s'<t\xeX 2- IxeX -dW x s (since dwsix is (conditionally) independent of Fg) = E/ y y [u2 dw2 + y u u ,aw dw 1 + 2x0 SA- v sx sx sx sx sx sx' \0<s<t xeX — — x'/x — — X " x:JteX I I E(U2 )AxdAt + 0 (using E(dW2 ) = AxdAt) 0<s<t xeX -- ^ 1 1 1 2 |u|| L (TxXxfi) Hence the mapping U > £ £ U dW acts isometrically from the ^ ^ v sx sx 0<s<t xeX --2 2 space of F -adapted' S-L (TxXxo,) processes to *L (fi) . Hence, if U and ' 2 U are any two S-L F^-adapted liftings of u , y I D dw - y yudw 0<s<t xeX — — 0<s<t xeX — — I I E|U. - U' |2 AxdAt 0<s<t xeX tx tx1 Hence °y £ U dW 0<s<t xeX a .s. o I I u' dW 0<s<t xeX §x §x A similar argument shows that if t1 RJ t , then a.s. y y u dw tv sx sx 0 , so that the r.h.s. of (2-1) is well defined t,<s<t xeX — up to a null set. We must show that the r.h.s. of (2-1) coincides with the Ito integral. Consider a process u of the form u. (o>) = Irj_" . . (t) • I. (x) • I_(u>.) , tx It ,t ) A K where R is an F measureable set in fi , and A is Lebesgue measureable 1 in R . Pick %2 ~ ^2 ' and P^c^- A internal in X such that LUx)(A A st 1(A)) .= .0; . By Theorem 3.2 in Hoover and Perkins (1983), we may ' find ' £ t - and an internal subset R of fi - , such that R e and P (R A R) =0 . _1 Let u\ (to) = I tx Then u dw = lv(w) sx sx K dw sx R [tl'tAt2)><A a .s. 1T)^ ' I I dw » (by. definition of the white t <s<tAt 2-noise dw. ) xeA tx I I 0<s<t xeX U dW sx sx Thus the left hand side and rhs of Theorem 2.3.1 coincide for simple functions u of the form indicated, and hence for linear combinations thereof. 2 Such linear combinations are clearly dense in the Hilbert space of L F^-adapted processes on, [0,..:t]-, as in the case-of the ordinary stochastic integral. Hence the two sides must coincide for all such u and U by the isometry property of the stochastic integral, and of the hyperfinite sum. CHAPTER THREE The Heat Equation with Non-Linear Stochastic Forcing 3.1 Scope In this chapter I will use hyperfinite methods to prove a weak existence theorem for solutions of equations of the form (3-1) |» f(u) * , dt 2 tx dx + * + on R x R . Here W is "white noise" on R x R . We suppose the function tx f to be continuous, and to satisfy a growth condition: there is a real K , such that (3-2) f2(u) £K(l+u2) for all u . We solve C3-1) subject to an initial condition C3-3J. u(0,x) = u (x) , where u^(x) is a bounded continuous function. With minor changes in notation f could be made to depend on t and x as well as u , and u^(x) could be taken as a random function. With more significant changes, the same construction can be made to work if (3-1) is modified by the introduction of an additional forcing term + g(t, x, u) , or by the introduction of a bounded, non-zero, continuous non-linear function 92u of u multiplying the term —— . With appropriate inequalities analogous 3x to those in Appendix A, a treatment very similar to the remainder of this chapter can be done for equation (3-1) on a strip R+ x [a, b] with Dirichlet or Neumann type boundary conditions specified at a and b . A random initial condition u independent of the white noise W may be 0 tx handled by enlarging the probability space 0, . As discussed in section 1.2, (3-1) is only solvable in a 'weak' sense, oo that is, if <|> e C ( R) , C3-4) utx<$> (x) dx - u d> (x) dx = Ox u Ad)(y)dyds sy ft 0 f(u )<j>(y)dw sy sy This is called the "weak" form of (3-1). The main theorem (3.9.2) of this chapter asserts that there exists a space Q , such that, for any f , a stochastic process u , jointly continuous in t and x , may be defined tx 00 on 9, , for which (3-4) holds for any (f> e C (R) . c Section 3.2 introduces the probability space ft and discusses the construction of a white noise on 0, . Section 3.3 exhibits a hyperfinite analogue of (3-1) and shows how it may be solved internally, for a solution U tx Section 3.4 contains several inequalities which are used subsequently. Proofs are delayed until Appendix A. Section 3.5 contains estimates on the moments of the internal solution Utx , which are necessary for 3.6 and 3.7. Section 3.6, obtains bounds on the moments of spatial differences U. - U . In section 3.7-we obtain bounds on the moments of temporal tx ty differences U\ - U .We use the results of 3.6 and 3.7 in 3.8 to show tx rx that U is, with probability 1, a lifting of a jointly continuous process tx In 3.9 we verify that the process u actually satisfies (3-4), hence is a weak solution of (3-1). 23 3.2 White Noise on the space ft Let a positive infinitesimal Ax be given. Let X be the set {kAx|k e *Z> . Now pick aa such that 0 < °a < -j , and let At = a Ax2 . Now suppose t^ > 0 is given and let t be any number of the form {kAt, k e *N} such that t_ ~ tc . Let T be {kAtlk e *N, k < t, / At} . -ft 1 — -f We have mentioned in section 2.2 that if ft is a *countable space r -,TxX (such as {-1, 1} ) on which are defined a family of I.I.D. random 2 variables {E, t 6 T, x e X} such that E(E ) = 0 & E (£ ) = 1 V(t,x) , tx- - — tx — tx --and such that 5^ possesses finite higher moments of all orders, then the random (internal) set function A > T E /AtAx induces a , • \ •» tx (t,x)eA — "white noise" {w^(oj) |A e B(['0, tf] !x R) } . Further, this white noise is adapted to the filtration derived as per section 2.3 from the internal filtration F generated by the values {E, s < t} . It is with respect to t sx - _ -this white noise that we shall solve (3-1). The method we shall use is similar in spirit to that employed in Keisler (1983) for stochastic O.D.E.'s,/ 24 3.3 Hyper-Finite Difference Equations Let F(u) be a uniform lifting of f(u), subject also to the growth condition (3-2). Let be a uniform lifting of u^ which is also uniformly bounded on X . Then consider the hyperfinite analogue of (3-1) , TTa +U U . - 2 U +U F(U ) £ t+At,x t,x t,x+Ax t,x t,x-Ax tx tx (3-5) At . 2 Ax /At AJ or equivalently, U,.,.., = U + At Ax 2[u . - 2 U +U . ] + F(U ) E — t+At,x t,x t,x+Ax t,x t,x-Ax tx *tx Ax We may solve (3-5) in principle inductively. The specification of U U / X and EQ ^ for x e X gives us enough information to find . Knowing U.. and £, , we may solve for U . , and so on. At,.- At,- 1 2At,-, Continuing in this manner, we define an internal process U (u) . We may find a closed form expression for this inductive definition as follows. From (3-5) it is clear that the value F(U ) £ enters into the definition tx tx of U , then U . . , U , U . , and then five t+At,x ' t+2At,x-Ax t+2At,x t+2At,x+Ax values of U , and so on. The coefficient QnAt with which F(U ) E t+jAt, X~Y tx tx enters into the definition of U • „ ,may be found as the solution of a .trHiAt,y difference equation: = 1 ; = 0 Vx ^ 0 ; , (n+l)At nAt L „ x^nAt ^ ^nAt (3-6 Q = a Q . + (l-2a)Q + a Q . *x • *x-Ax x x+Ax But we may recognize (3-6) as the difference equation governing the density of a Markov Process, in fact, a random walk. Thus we have t Lemma 3.3.1. The coefficients Q- are the "density" of an infinitesimal random walk on^ the lattice X , starting at x = 0 , at time At , and taking one step to the right (or left) in each'time period At with probability a . With this notation, we may write C3-7L u tx y l Q- - F(U ) L Lv ^x-y Sy 0<s<t yeX - i -I ? 1^ sy •* Ax + I Qt+At u y,X °l The second term on the r.h.s. of (3-7) is the deterministic solution to the hyperfinite heat equation (3-5), which we will designate hereafter U tx 3.4 Some Useful Inequalities We will find the following inequalities helpful in the next three sections. Proofs are deferred to Appendix A, since they do not illuminate the particular subject matter of this chapter. Lemma 3.4.1. There is a finite constant K , depending on a , such that for a 2 i t/A?t c *N , I (Q-) < K ^ v x — a J t xeX Lemma 3.4.2. There is a finite constant , depending on a , such that 2 for all t , I I CQ-j < K /t/At . " 0<s<t x€X S - a -Lemma 3.4.3. There is a constant K , depending on a , such that V z e X a I I CQ- - Q-A )2 < K |z| / Ax . 0<s<t -xeX x ~ a " Lemma 3.4.4. There is a constant K , depending on a , such that, for all a t e T , and r < t I I CQ-"") - Q- ")2 < K /(t-r)/At 0<s<r xeX We shall also require the following theorem of Burkholder (1973, Theorem 2.1.1., specialized slightly). Theorem 3.4.5. Let M , n e N be a hypermartingale and let <M> n e N n n be the associated predictable square function, and let p > 1 be finite. Then there is a finite constant K depending only on p , such that P E(MP) < K EC<M>P//2) + K E(max |M - M |P) . -n-p- n P- Q<k<n k+1 k 27 3.5 Bounds on Moments of U tx Pick q > 1 , and let R (t,x) = E|U |2q . Let H (t) = sup R (t,x) q - - - tx q - q - -Consider any fixed (t,x) / and let ji be the measure on [0, t] x X 2 defined by y({(s,y)}) = (Q-"§) — . Let lyl denote y([0, t] x X). - - X~Y Ax 1 1 Let c denote, in what follows, a finite constant, depending only on q and t , which may change its exact value from line to line. Now 1/2 R Ct,x). < c E q - - - - I I ECU ) Q±~S- E (f) 0<s<t y£X SY X"X sx Ax 2<5 r |2q + c U 1 tx1 In the summation above, the term corresponding to (s, y) may be regarded as a martingale increment with respect to the filtration F . Hence we may s apply Burkholder's inequality (3.4.5) to obtain R (t,x). < c E[ I I F2(U ) y(s,y)]q (3-8) + c E| sup| F(u ) |E /y(s,y) |2q + c lu |2q . — sy sy - - 1 tx1 0<s<t -1 — yeX The second term on the r.hs. of (3-8) is clearly equal to 2 2 ^ c E[ max F (U ) E y(s,y)] which is bounded by - 0<s<t ?y - -c E[ I I F2(u ) E2 y(s,y)]q . Thus " 0<s<t veX a S " R(t,x)<cE| I I F2(U )(1 + E2 ) Hil^-|q|y|g + c|u I2q 0<s<t YeX tl 51 |y 28 Applying Jensen's Inequality to the probability measure y(-)/|y| , we find R Ct,x) < C E[ I I |F(U )|2q(1 + \^\) Hitfl]|y|q + C|S \2* . q " " ~ ~0<s<t y£X ?y \V\ ' tx Now we use the facts that |y|q/|y| = | y |q 1 , that |y| is bounded uniformly for °t <_ t by Lemma 3.4.2, that the moments of £ are finite, 1 sy that 5 is independent of U , and that U is bounded uniformly sy sy tx in Ct,x) , since UQ is bounded, to obtain R(t,x)<cU+ I I E |Fltr;>|2q(l + E\K |2q) (Q-"-)2 7^) q - - 0<s<t yex ^ *-y A* < cCl 0<s<t I yeX K2q(l + E|U -' sy |2q) (Q^) *x-y Ax; <_ cCl + I I R (s,y) CQ- -)2 T^) 0<s<t yeX q - - x-y Ax < ca + l H (s) c I —fp- At) . 0<§<t q yeX AX Now using Ax = /At/a , and lemma 3.4.1 we find Lemma 3.5.1. r -1/2 (3-9) H (t) = sup R (t,y) <_c(l + I H(s)(t-s) At) , for all t in y q °fs<t [0, t^l , where c is a constant depending on a , q and tf . Now, iterating (3-9) and integrating by parts, H (t) <c I (1 + c I (1 + H (u))(s-u) 1/2At)(t-s)"1/2At 0<s<t 0<u<s < c (^tj+ I (1 + H(U))( I (s-u) 1/2(t-s)1/2At) At 0<u<t '* u<s<t ~ ' (3.10) <.c(l+ I (1 + H (u)) ( I s V2(tf-s) 1/2 At) At 1 c(l + I H (u) At) 0<u<t q where c is another constant depending only on a, q, and t^ . We now require a type of Gronwall's lemma. Lemma 3.5.2. There is a constant c depending only on q, maxlu |2q . a y 0y and tf , such that E_| | <_ c exp(ct) , for t <_ t , and for any x e Proof. We may take c to be the maximum of H (0)(= maxlu |)2q and the q y 1 Oy1 * + of equation (3-10). We proceed by induction, on n e Z . For n = 0 , the lemma holds for t = nAt . Suppose now that for n e N H CkAt) <_ c Cl + c At)k , for k = 0, 1, . . . , n-1 . Then HCnAt) <_ c (.1 + c £ (1 + c At)k At) (by (3-10)) 0<k<n = c ci + c ^—r-1— At cAt (summing a geometric series) = cCl + c At)n This is the induction step. By the transfer principle, we may conclude that this internal argument verifies H(nAt) ' <_ c (1+cAt) n , for all n e *N such that nAt <_ t . Then we notice that (1+cAt)n £ exp(cnAt) , (in fact they are infinitesimally close). • 31 3.6 Bounds on Moments of Spatial Differences Let U = U + V , where V represents the contribution from the 1—X L.X "CX random forcing. U is S-continuous (in fact S-smooth), and we are interested in the continuity of V . In this section we obtain estimates on the moments of the differences V - V tx ty Lemma 3.6.1. There is a constant c , depending only on q > 1 , max|u2q| , a , and t^ , such that, for all x, y in X , and for 0 <_ t <_ t , i 12q I I q -|vtx - v ± °lx - zl • Proof: We may write V - V as a sum of martingale increments with tx ty respect to the internal filtration F ; by (3-7), s V - V = y l (Q-~- - Q^)F(U ) £ (—)V2 £ a. 0<s<t zex V* V«' H We will estimate the 2q-th moment of this using Burkholder's Inequality CTheorem 3.4.5). We will designate by yCs,z) the measure on [0, t] x X which assigns t-s t-s 2 At to each point (s,z) the weight (Q- - - Q- -) — - - x-z y-z Ax Hlv^ " Vty|2q < c EC F2(u )dy(s,z))q sz - -[0,t)xX 2 2 + c E( max F (u ) £ y(s,z))q " 0<s<t ^ -S-Z " -32 < c E( [0,t]xX F2(u , dy_(s,z))qiijIq -- y + CE( [ F2(U ) K2 ^S'5),q|y|q sz sz < C E CF2(U )]q[i + |5|2q] T^-(§'5) -|y|q sz sz I I (Jensen's Inequality applied to -—r) < c K2q H (s) -^{s-'*] • |y|q q " y (using independence of U and E , and also finite moments of E ) sz sz sz <_ c I x - y | At (using lemmas 3.5.2 and 3.4.3, noting that — = ctAx ) • 33 3.7 Bounds on Moments of Temporal Differences In this section we obtain estimates on the moments of the differences V - V tx rx Lemma 3.7.1. There is a constant c , depending only on q > 1 , max |TJ | , a , and t , such that, for all t , r <_ t , and all x e X , uy r — — -r -I I2q | ,q/2 EV -V<ct-r —' tx rx' — '-Proof: We may suppose w.1.o.g. that r < t . We may write V - V as a sum of martingale differences with tx rx respect to the internal filtration F s S £x 0<s<r zeX IS ^ 25 ^ ^u ^. ~V x-z sz sz Ax r<s<t Z£A - -We will designate by y(s,z) the measure on [0, t) * X which assigns t-s r-s ^ At to each point (s,z) the weight (Q- - - Q- -) -— , if s < r , and the - - x-z x-z Ax - -2 weight ^2X_|^ Ax ' ^ - - - < *~ ' T^en usi-n9" Burkholder's Inequality CTheorem 3.4.5) I l2c* E V - V < c E( — tx rx — — J [0,t)xX F2(U )dy(s,z))q sz + c E( max F2(U ) £2 y(s,z))q - 0<s<t ~ ~ < C E (F2(U ))q(i + \K |2q) -^'5>.|y|q y (using Jensen's Inequality) < c E K2Q H (s) ^ • |ylq (using independence of U , and E, , and the sz sz f initeness of E I E, I 2q) — sz1 < c(t-r) q/2 (using lemma 3.5.2, and lemmas 3.4.2 and 3.4.4, with 34 35 3.8 S-Continuity and the Standard Part The main result in this section is that , the solution to the hyperfinite difference equations (3-5) is a.s. S-continuous. We shall obtain this by applying a non-standard version of Kolmogorov's Continuity Criterion: * Theorem 3.8.1. Let : fi x r -»- R be an internal process on a hyperfinite lattice F which represents a finite rectangle in R . If there exist positive real numbers 3 ,...,3 , y ,...,y , K such that for k = l,...,d h- d 1 d k I i k Uux _ U I <_ K I x - y | , whenever x, y e T are such that x - y lies along the k*^1 coordinate axis, then if 6, < f, /$, , k = l,...,d there k k k o is a set fi' c fi of Loeb Probability 1, a function 6 (w) , <5 (w) > 0 on fi' , ,6k and a constant c , such that for k=l,...,d U - U < c x - y ' x y' — ' - - 1 whenever x, y e T, |x-y| < 6(w) and x - y lies along the k^1 coordinate axis. In particular U is a.s. S-continuous on T . Proof: See Stoll (1984)t Lemma 3.2. The result he states is not as detailed as 3.8.1, but his proof is sufficient. Theorem 3.8.2. The hyperfinite process u-(-x^a)^ constructed by the solution of C3-5), is a.s. S-continuous on near standard points in T x X . Moreover, if < 4" and ^2 < \ and A c X is a rectangle whose sides have finite length, there is a set fi' <= fi of probability 1, a positive real function 6(a)). on fi' , and a constant c , depending on 3 , 3 max U , a , and 1 2 yeX °l V t , such that wefi' , x, yeA,t,r<t 36 1 y lUtx ' Ury' - c('t ~ ^1 + |x - y| ) , if |t-r| + |x-y| < 6 . Proof: Pick q e R+ such that 8 < -—— and 8 < q 2 . By lemmas 1 2q 2 2q 3.6.1 and 3.7.1, there are constants c such that I |2q I 12+(q/2 - 2) E V - V < c t - r ^ ; —1 tx rx1 — '- -' I I2q I i2+(q-2) —1 tx ty — 1 - -1 Hence by Theorem 3.8.1, the statement of the lemma is true with V in place of U on any set T x A , where A is an internal finite rectangle in X ,-but the near standard part of X is a a-union of such A . Now U = U + V and U is a lifting of a smooth function, which is-a-solution-to the heat equation. An examination of the explicit form (3-7) for U yields quickly that U. - U <c(t-r + x- y) tx ry — — — — — Hence the theorem is true for U = U + V . • In general the exponents j and ^ are best possible (see Walsh (1986), Corollary 3.4). We may allow slightly more general initial conditions, if we are prepared to relax the conclusion slightly. The arguments in lemmas 3.6.1 and 3.7.1 depend only on the boundedness of max |u I . Thus the boundedness of U y£X ' °T 0 is enough to ensure the S-continuity of U - U . However, in this case U itself will not be S-continuous in the monad of zero. If TJ is bounded but 0 discontinuous the conclusion of the theorem will have to be restricted to 37 °t, °r > 0 . Returning to the case when is continuous, we find that another way of phrasing the conclusion of Theorem 3.8.2 is that U is nearstandard almost surely in *C([0, tf] x R : R) . Hence we may define a process U(OJ) o as the standard part of U(OJ) , or equivalently u = U for all o o tx t, x t, x , a.s. It is clear that that set ft1 = {co :\/n e N, 3k e N,V x, y £ X , t, r<t,|x-y|<— , It - rl < — =>|u - U | < —} is in the -- — -f 1 - - ' k 1 - - — k tx £Zn a-algebra generated by the internal sets, hence is Loeb-measureable. Thus the process utx(o» = i °Utx^oO" for any (t,x) K (t,x) , if co e ft' 0 , if co i Q1 has sample paths in C([0, t ] x R : R) . We are of course really interested in solutions for all t . Thus take the hyperfinite time line T up to some infinite number L , and take ft larger enough to support a white noise on T x X . Construct the solution as before. All theorems proved previously hold true up until any finite time t . Thus we have Corollary 3.8.3. The solution on T x X , where T now is an infinitesimal grid representing R+ , constructed from the difference equations C3-5I is S-continuous in (°|x| < «>} n {°'t < °°} , a.s. o We note also that u is 0CF ) measureable for any s with s > t t-s - -hence by definition Cin section 2.3) u is F^_-adapted. 38 3.9 Solution of the SPDE We now show that the process u of section 3.8 is in fact a solution tx of (3-1) with respect to the white noise defined on 0, in sections 2.2 and 3.2. We must check condition (3-4). We need first a new definition. Definition 3.9.1. An internal function $ on an infinitesimal lattice X £ Rd , is called a lifting to order k of a function cj> : Rd -> R , if °6 (. . . (6 $)...) (x) = <J) (°x) for all x e X . Here 6 x. x. 3x, ...3X. - . *1 \ \ \ is the finite difference operator in the direction x.: (6 ($)(x) = i x. l [§(x,,...,x.+Ax,...,x ,) - $(x,,...,x )]/Ax . A consequence of this definition -± ~i -d l -d is that for nearstandard x, °[$(x+Ax) - 23>(x) +$ (x-Ax) ]/Ax2 = -^r- <j>(°x) 8x if $ is a lifting to order 2 of (j) . We note that every C function cj> has a canonical lifting to order k , * namely <f> restricted to the lattice. oo Now fix any e C ( R) , and let $ be a lifting to order 2 which is o exactly 0 on values of x whose standard parts x lie outside the closed support * (j) , (to avoid (unnecessary) concern over the convergence of -countable sums) . Then U • $ (x) is a (uniform) lifting of u d> (x) a.s. for any t ss t . Thus tx - tx | utx<Hx)dx - I u0x«Mx)dx a=S* °l (Utx - U0x) «(;, Ax 39 xeX 0<s<t u 4 - u s+At,x s,x At •At • $(x)-Ax ' I 1 0<s<t xeX U . 2 U + u SyX+Ax sx s,x-Ax Ax F(U ) sx $(x) + —TIZIZ- ? $(X) /At Ax --Ax At (3-11) I 0<s<t xeX 0(x+Ax) - 2$(x) +'$(x-Ax) I U _ - Ax At sx Ax + ° I I F(Ugx) $(x) /AtAx E 0<s<t xeX sx Now U C$ (x+Ax) - 2$(x) + $(x-Ax)]/Ax is a uniform lifting a.s. of u <f>" (x) sx sx Hence the first term on the r.h.s. of (3-11) is a.s. u c(>"(x)dx • R sx Now F(U ) $ (x) is a (uniform) F -adapted lifting of f(u )4>(x) sx — s sx Hence by Theorem 2.3.2, the second term on the r.h.s. of (3-11) is a.s. equal 't to 0 R f (u )<j)(x)dw . Thus we have, d) e C ( R) t e R sx sx c (3-12) u <i> (x) dx -tx R R u c|> (x) dx a=S ' Ox R u (j) (x) + sx 0 R f(u )<Mx)dw S ,X SX Theorem 3.9.2. There is a Loeb space , on which any equation of the form (3-1) has a solution jointly continuous in t and x with respect to the canonical white noise on fi . NOTE: We believe it is possible to extend Keisler's internal transformation principle to this Loeb space fi (see Keisler (1984)),. In this case, the equation (3-1) has a solution with respect to any white noise on R supported by 0, . Walsh (1986) has established uniqueness for the case when f i Lipshitz. Presumably this is false in general but we do not know a convenient counterexample. CHAPTER FOUR The Dawson Critical Branching Diffusion 4.1 Introduction The use of hyperfinite difference equations, which was successful in Chapter 3 for one spatial dimension may be attempted for higher dimensions as well. This approach does not succeed to the same extent as it does in one dimension, as is spelled out in Appendix B. However it does succeed with the Dawson Critical Branching Diffusion, as will be explained here. It is possible to represent the Dawson process non-standardly, by placing, on a hyperfinite grid X , representing a portion of R , d->"-l , an infinite number of particles, each executing an infinitesimal random walk, and undergoing branching. If we let U^ stand for the density _of„particles at any grid point x e X , at time t , then, if the initial, density ;. is- taken large enough, it is possible to show that U satisfies a hyperfinite difference equation of the form °*tu = A U + Z - tx where Z is an internal noise with E(Z. ) =0 , E (Z u tx' ' tx|Utx' _ d * — — AtAx Hence Z may be written as /u w where W is like an S-white noise in many respects. This noise W is however, a little awkward to work with, and therefore we adopt the simpler scheme of difference equations set out in 4.2. In 4.4 we examine the total mass of the process constructed in 4.2, and use this in 4.5 to establish some continuity results, which yield easily that the standard part is well-defined. In 4.6 we verify that this standard part does indeed coincide with the Dawson process. In 4.7 we obtain several new results about the pathwise regularity of the Dawson process, using nonstandard construction. 43 4.2 A Hyperfinite Difference Equation Let Ax be any infinitesimal, and let X be {xix = (k Ax,...,k Ax) , k. e Z , d > 1} so that X represents R . - - 1 d I W(: will treat here the construction of the Dawson Process only on the whole of R . A very similar treatment is possible if reflecting boundary conditions are imposed on several hyper-planes in Rd , or along the edges of a rectangle in R , but the inequalities are messier, and indeed depend on those for the unbounded domain. Let At be an infinitesimal, such that At/Axd £ 0 . This makes some parts of the treatment much easier. In case d = 1 we require 1 2 At < j& , as in Chapter 3. Let T be a hyperfinite time line of spacing At : T = {t : t = kAt , k e *Z+ , k £ M} . We will suppose 0 < ° (MAt) < °° . Let t = MAt . TxX It is easier, and it suffices for our purposes, to take fi = {-1, 1} and to let E (co) be the coordinate map, as-.outlined in section 2.2. Thus tx - •. . PCE = 11 = P{E = -1} = i . - tx tx 2 The analogue of equation (3-5) in higher dimensions is: C4-11 C6.U) = CAU) + A—- IE , t tx - tx \ ./ d 2At I tx or equivalently, U , = U + AtCAu) +1 / A 1 E , Ut+At,x utx tx % A 2 / tx (r: •) where 6 Is a finite difference analogue of — : t 9t (6tU)tx = ^ut+^t x - ut x^ / At , and A is the finite difference analogue of the Laplacian in R : C4-2) CAU) 2 /n (Ut(x ,...,x.+Ax,...,x ) x j_i=l - -1 -l -d + U4-/ A \ 2d U . t(x ,...,x.-Ax,...,x ) t(x,, .. .,x ) - -1 -i -d - -1 -d The term /u At u tx is substituted for simply / U At tx in (4-1), Ax - - * Ax-in order to ensure, that, a non-negative value at (t,x) , all of whose nearest neighbours are non-negative, will not become negative at the next time step. The values of U for which the linear term is taken are 4At 0 < U < Ax d ' For an initial condition for (4-1) we may use any non-negative internal def r . d function on X , requiring only that I u Ax , the total xeX Ox initial mass, be finite and that the mass on points of X which are not near -^standard be infinitesimal in sum. We may represent any finite positive d d Borel measure on R (we will denote the space of all such measures ( R )) by such, an internal function UQ (see Cutland (1983) Theorem 4.7 and preceeding remarks). 45 4.3 The Coefficients Q "u At U SY sy From (4-1) and (4-2) we observe that each value of / A —— £ . d 2 sy Ax --enters into the definition of subsequent U 's . We denote the coefficient tx t-s of the former term in the definition of the latter by 0~ ~ , observing that x-y these coefficients are homogeneous in space and time. We may then write the analogue of a Green's function formula: U At U sy £ tx K~y \ * A D 2 I sy y x-y Oy 0<s<t Y£/^ - - v Ax / -- yeX - -t+At Lemma 4.3.1. The coefficients Q- are the internal density for an internal infinitesimal random walk B on X . The standard part of B t t is a.s. d-dimensional Brownian motion of rate 2. Proof: As in Lemma 3.3.1 we may construct difference equations for the coefficients Q- from (4-1) and (4-2). We observe that these difference x equations correspond to a Markov process, B , on X , with parameter t e T where B^_ starts at 0 e X , and at each time step B takes a step to one of the 2d nearest neighbours of its current position, with infinitesimal 2 probability a = At/Ax for each of the 2d possibilities. The process stays put with probability 1 - 2da ~ 1 . Q-+^t = P(B = x) . X t ~ i th Let B denote the i coordinate of B . Since steps in opposite t t directions have equal probability, B1 is an internal martingale. Also E(:(Bt+A-t " B^).2 | Bs ; 0 <_ s <_ t) = 2aAx2 = 2At . Hence the (internal) predictable quadratic variation <Bl>t = 2t • Using Burkholder's Inequality i ±2 on the higher moments of , we may conclude that (B^_) is S-integrable t . Now a step in one direction excludes a step in any other direction during the same time step. Hence E[ (B-1 . - B3) (B1 „ - B1) ] = 0 if i ^ j . t+At t t+At t Successive steps are independent, so that the internal process <B1, B^> = 26..t . Now we invoke Hoover and Perkins (1983) Theorem 8.5 to assert that B has a.s. a standard part b , and we observe that b satisfies i) b^ is a martingale i = l,...,d ii) <b1, b3> = 26 t i, j = l,...,d t i: (again with reference to Hoover and Perkins (1983)). Now i) and ii) above characterize d-dimensional Brownian Motion. • We need a little more information about the Q's . o Q- 2 x" II || Lemma 4.3.2: If °t > 0 , —= P (°x) , where P (x) = 1 exp{ } :—^—' ~ i d t — t r— 4t Ax - » 4TT t Q-x is the density for Brownian motion of rate 2. Hence is S-continuous in Axd t for °t > 0 . Proof: We know from 4.3.1 that the distribution for B is the same as that for bQ^_ , which has density P0^_ , for any °t > 0 . The statement of the lemma will follow then, if Q-/Axd is s-continuous for °t > 0 . x We observe that the following equation holds: C4-4I (6tQx) = (A Q^)x Now it is clear from the definition of the random walk B that the coefficients Q- are symmetric with respect to change of sign on any of the d indices x Claim: For each t, x , Q- > Q- . _ for k = l,...,d , where e - - x — x+Axe, k - . k represents a unit vector in the x direction, which points away from 0 .. We prove this by *-finite induction on t . It is clearly true for t = At At x Q = P(B = x) =6- . Suppose the claim holds for some t e T . Then by x 0 - 0 - * (4-4) and (4-2), a C4-5) Qt+ t = (l-2da)Q- + a ( £ Q- . ) , and x x . . x.+ Axe j=l d t+At t , r t (4-6) Q _ . = (l-2da)Q- _ . + a ) n- _ Now if x 5^.0 then all terms appearing in (4-5) and (4-6) above lie -k in the same half-plane; each term in (4-6) is shifted by e Ax relative to, k and hence, by the induction assumption, is not greater than, the corresponding term in (4-5). Hence Q- > Q- , - . x - x+Axe, - k Now suppose (w.l.o.g. k = 1) x = 0 . We pick one of the two possibilities for $ and stick with it. Let y be a (k-1) tuple. By our assumptions on At , in any dimension d , (2d+l)a <_ 1 . Hence (4-7) Cl-(2d+l)a)Q§x > (l-(2d+l)a)Q^x > (1-(2d+2)a)Q|x^ + a p*^ Now adding a O- + a 07 + a Q~ to both sides of (4-7) and using Ax,y -Ax,y t t , Q- = Q- , we obtain' -Ax,y Ax,y (4-8) (l-2dcO + a(Qj + Q- ) >_ (l-2da)Cj£ + a (Q^ + Q| ) Oy Ax,y -Ax,y Ax,y Oy 2Ax,y Now by the induction assumption d t d t C4-9) a 1 Q5y+Axe. > « I 2Zx,y+Axe ' j=2 uy- j j=2 j Adding (4-8) and (4-9) we obtain t+At t+At 25y ^Zx.y * This establishes the claim for {Q^+^t|x e X} . Now we invoke the principle of induction under the transfer principle to establish the claim for all tel. Now by Lemma 4.3.1, for any internal rectangle A X d d xeA Ax (4-10) I —=- Axc P (x)dx , for °t > 0 In light of the monotonicity claim just established, this means °(Q-/Axd) = PQ (°x) if °t > 0 , and x e ns(X) . Therefore also is S-continuous in t for t > 0 , since P (x . d - - t Ax is continuous in t , in t > 0 . We prove the following inequality in order to obtain moment bounds in section 4.5. It has no independent interest. We introduce the notation A:Y- xeA VY Lemma 4.3.3. There is a constant K < °° and a positive infinitesimal At' 49 such that for all internal sets A , and all y e X , whenever t, r, s e T and r-s > t-r > At' , t-s r-s Q- - - Q- -A:y A:y (t-r) < K-- (r-s) Proof: Suppose first that °(t-r) > 0 Now clearly for fixed t, r and s , max IQ- - - Q- -1 = Q- - - Q- -A,y A:y- A^ Vt,r,s:y \'t,r,s:y where A = {x e X|Q- § < Q- -} y,t,r,s - x-y x-y The standard part of this set is easy to identify, using Lemma 4.3.2. It is the ball A ^ = {x e Rd| Mx-y||2 < 2d ^-s) (r-s) c ^r y,t,r,s 111 Jn2~ (t-r) r-s ' where y = °y , t = °t , r = °r and s = °s . The radius 6 of this - - - - t,r,s ball is computed by solving P^ ^(x) = P (x) -i.e. [4Tr(r-s)] d/2 exp[- 4(r-s) -] = [47T(t-s)] d/2 exp[- 4(t-s) -] d/2 r-s 4 (r-s) (t-s) .2 oJ (t-s) r-s . rt-s, <=> 6 = 2d logC ] t-s r-s Now (P (x-y) - P (x-y))dx r-s J t-s y,t,r,s 50 a-l' • t,r'S 1 , 2,«< „ d-l exp(-p /4(r-s))p dp , (4tr) d/2 J (r-s) d/2 d-11 (4TT) .d/2 TT) } t.r.s 1 , 2,., d-1 ' ' exp(-p /4(t-s))p dp (t-s) d/2 where S is the area of the surface, S, , , of the unit ball in R d-l d-1 Now let p1 = p^/tTs" in the first integral, and p' = p/2/r-s in the second, to obtain (4TT) d/2 j [ t,r,s/2/r-s ^ _p2 P e t,r,s/2/t-s d-1 -p Since the function f(p) = p e is bounded, we may find K to bound this integral by K 6 (~i t, r, s / r— vr-s vt-s K/log (1+ ) ( * r-s /t-s - /r-s, /t-Now if r-s >_ t-r we may bound this further, by KpEE ((!2}"(r's) _) r S /t-r(/t^s~ + /r-s) < K t-r r-s Thus, if °(t-r) > 0 (4-11) r-s>t-r=> sup \c£ -(t-r) Q£-s. I < K -7^T- , Vy e X AcX A:y A:y — (r-s) 51 Since the internal statement (4-11) is true for t - r > e for all real e > 0 , then by the principle of infinitesimal overflow it must hold for t - r > At' ss o . • We believe that this Lemma is true for t-r down to At but an internal proof of this is not easy. 52 4.4 The Total Mass Process M Now let M = ) U Ax be the total mass of the internal measure t Lv tx xeX — whose density U is obtained by solving (4-1) inductively. We have by definition, /U At U 0<s<t yeX xeX i~l y Ax r r r t-s / sx d In what follows c denotes a flexible constant whose values depend only on q and t , and which values may change from line to line. Lemma 4.4.1. There is a constant c depending on q > 1 , M , and t^ , such that \jt <_ t , E(M2q) < c eCt Proof: Using Theorem 3.4.5 (Burkholder's Inequality), the fact that v t-s ) Q- - = 1 and the fact that for any a, b, a A b < a v x-y — xeX - i E(M2q) < c E(| I I 1 • U AtAxd| ) - t — - 1 ^ L sy 1 0<s<t yeX -£ d ,2 • * (4-12) + c E( max U AtAx £ ) " 0<s<t ^ 5* yeX + c M2q 0 Now the second term in (4-12) above may be bounded by c E( max | M At|q) <_ c E| \ M At|q . Hence (4-12) becomes 0<s<t - 0<s<t -0<s<t r 0 E(M2q) < c EI T M At - t — - „ _ s < c E| I M Atl^ + c(M2q+D ~ ' 0<s<t 5 0 (4-13) < c y E(M2q)At + c(M2q+l) — « . s 0 0<s<t for all t <_ tf . We complete the proof with an appeal to the appropriate version of Gronwall' Lemma (Lemma 3.5.2). We may now estimate the differences M, - M as follows, in order to t r obtain continuity of the mass process. Let q > 2 , and w.l.o.g. take 0 <_ r < t <_ t • Then (4-14) E|M - M I25 < c E I J 1 U AtAxd|q t r' — - L i_ v sy 1 r<s<t yeX --(using Burkholder's Inequality and applying the same reasoning as in going from (4-12) to (4-13)) q < c E | I M " (t-r)q r<s<t ^ (t-r)1 - -< c E ( I Mq -^) (t-r) q s sx. s t-r r<s<t - - -(applying Jensen's Inequality to the probability measure —— ) on the interval r < s < t) t-r - - -54 < c I Ellfl) (t-r)q ~ r<s<t " S (t-r) - -(4-15) <_ c exp(c tf) (t-r)H , by Lemma 4.4.1. Thus we have Lemma 4.4.2. The total mass process M / defined above, is a.s. S-continuous on 0 < t < t . Proof: Apply the Kolmogorov Continuity Criterion (3.8.1) to (4-15). • In fact M is S-Holder continuous of order — - e for any E > 0 . 55 4.5 S-Coritinuity of the Process For internal sets A , and internal functions F let the notations X^ and X*" stand for T U Axd and J F(x) U Axd . We introduce the t u tx Lv - tx x£A — xeX class of (standard) functions C2 = {f e C2( Rd) | 3K < °° such that |f(x)| + I Af (x) | <_ K \/x e Rd} . Recall from section 3.9 that if an internal function F lifts to order 2 a function f e C, „ , then b, 2 °(AF)(x)= Af(°x) , Vx e ns(X) . Lemma 4.5.1. If F is an internal function such that |F(X) | + |AF(X) | <_ K , Vx e X , then the process XF is S-continuous on T a.s. Proof: For r < t e T , 7* - XF = Y F(x) (U - U )Axd t r y - tx rx xeX r s+At,x sx-i l IF(x) L ————J AtAx r<s<t x r r/Usx Usx i = y i F(X) (Au ) - A ; , ^ ^ L.-S-X VA A d 2At sx r<s<t x — — V A+Ax AtAxd (4-16) = I I (AF)(x) U Axd At , SX r<s<t x u sx . d + I [F(x)(^ A _p g /AtAxd ^ ^ L - sx 2 At sx r<s<t x 56 The first term in (4-16) above we recognize as I X At r<s<t -AF Now X- < K M . Hence if q > 1 s s -. -p j?. 2q _ 2q EX - X < c E( ) M At) -' t r1 - - L s r<s<t -+ c E( I I F2(x) u AxdAt)q r<s<t xeX (using Burkholder's Inequality for the second term on the right and incorporating the term involving "max" as we did in (4-12)) < c I E(M2q) (t-r) r<s<t 2q s t-r - -' + c I E( I F2(x) U Axd)qAt — \/ — SX r<s<t xeX r | F2,^ F2 2 (we now recognize the second terms as l E | X | At and X <_ K M ) r<s<t - - -(4-17) < c (t-r)2q + c (t-r)q , if t t , using Jensen's Inequality a second time. Now we apply the N.S. version of the Kolmogorov Continuity Criterion (Theorem 3.8.1). • 57 A Lemma 4.5.2. Um sup °X n = ° a-S* where n-*>° t e [0,tf] " A = {x e X| |x| >_ n} ; that is, X is nearstandard in M ( Rd) for all n t F t , a.s. Proof: Let H be infinite. Let F be an internal function such that 0<_F<_1,F=1 on A , F = 0 on A^ and AF is bounded. Then H H—1 — E(<> < ECX/"1) -1(1 Q-t+At) U Axd . Now by Lemma 4.3.1, for each y e ns(X) r -t-s r „t+At ) Q s 1 , so that e = ) Q- z 0 . By assumption UQ is nearstandardly concentrated, so that lim °( I U Axd) = 0 . n-X» y£A -n Hence, (4-18) Etxf) < I e U ~ 0 => xl s 0 a.s. t "yeX X °Z t F Now by Lemma 4.5.1 X is a.s. S-continuous. There is a countable S-dense subset of T for which (4-18) holds. Hence A °XF = OUteT,a.s. => sup °X H = 0 a.s. ~ " teC0,tf] £ A ° n P Therefore sup X —> 0 as n -> 00 , teCO,tf] £ A A n n+l , n but since X > X we must have convergence a.s. U 58 Theorem 4.5.3. X : T M ( R ) is a.s. S-continuous. Thus X is F nearstandard in *C( R ; ( Rd)) , and (°X)q = L(X ) 0 st-1 for all t t e [0, t ] a.s. Proof: Recall the weak topology on Mp( Rd) : + x <=> for all bounded continuous f , x (f) -»- x (f) . Let {(j) } be a countable collection of n Tk 2 cl CT _( R ) functions which constitute a convergence determining class for b, 2 M ( R ) . Let {$ } be liftings to order 2 of the {cj> } . Then by Theorem •t k k $ 4.5.1 each X^_ is a.s. S-continuous. Now by the Loeb construction (see section 2.2) (4-19) <h d L(X ) o st = °X k Therefore the l.h.s. of (4-19) is continuous a.s. for each <j) in the convergence determining class. Hence, almost surely L(X ) ° st-1 = L(X ) 0 st-1 t s for all t s s e T . * d From Lemmas 4.4.1 and 4.5.2, X is in fact nearstandard in ( R ) for all t , a.s., so that ° (X^) exists and equals L(Xfc) ° st~J-\^ t a.s. • Now we show a stronger form of continuity. * d Theorem 4.5.4. Let X be nearstandard in M ( R ) (as discussed at the ' ' 0 F end of section 4.2). Let A be internal. Then X is S-continuous for t °t > 0 on a coarser grid T' c T of infinitesimal spacing At' independent of A Proof: We take At' from Lemma 4.3.3. Let Y1= T Q-+At U Axd and let t Lv A:y Oy TEA - -59 U Y2 = t I I Qj.- ( Jv AtAx' 0<s<t yeX A:Y sy •^Axd) £ 2 sy Then XA = vJ; + Y2 t t t From Lemma 4.3.3, there is a K such that, if r < t < 2r in 7" , |Q- - Q- I < K(t-r)/r . Hence if t ~ r and °r > 0 - - - A:y A:y — - - -I |0j -Of | U Axd yeX A:x A:y °y < — (t-r) M — r - - 0 ~ 0 . Thus Y^ is S-continuous in T' n st-1(t>0) . Now suppose £<teT' ,t-r<l,0<y<lf and q >_ 2 . Then using Burkholder's Inequality and absorbing the terms involving a maximum into the summations, as we did in section 4.4, ((4-12) and (4-13)), i 2 2 i 2c2 E YI - Y < c E — t r1 — + C E + c E < C E 0<s<r - (t-r) Y ye Y (Q- - - Q- -) U AtA: u^ A:y A:y sy d q I v yeX X Y sy r-(t-r)'<s<r i I I (Q|.5)2 U AtA < ^. A:y sy r<s<t y£A u -± I 2-2Y Ct-r) M At 0<s< r-(t-r) (we invoke Lemma 4.3.3 and note (r-s) > (t-r) .) 60 + C E q 1 1 • M At r-(t-r)Y<s<r -(since |pjf - - -| < 1) 1 A:y A:y' — + c E T M At 1 r<s<t -q (4-20) < c((t-r)(2-2Y)q+ (t-r)Yq+ (t-r)q) (following the usual Jensen's inequality train of development, by now so familiar). Our faithful servant the constant c has now acquired dependence on Mo ' 01 ' -f ' and q ' kut not on t and - • 2 The best choice of y for (4-20) seems to be — . Using the Kolmogorov criterion (Theorem 3.8.1) we get S-continuity with a modulus -j - e for any e > 0 . • Remark: Of course if puts all its mass on a null set A , we get failure of S-continuity in the monad {t « o} . Corollary 4.5.5. Let F be any bounded internal function on X . Then F i — 1 the process X^_ is S-continuous on T n st (t>0) , a.s. Proof: Let K be a bound for |F| . Let XF = y y F(x) o~+At U Axd, and t v v ~ x-y Oy yeX xeX - *• -U let zf = I I • ( I F(x) QI~JT) (/u_AtAxd A Axd) ^ yeJ t w v - x-y sy 2 sy 0<s<t yeX xeX _ £ -£ z °£ F "F ~F . „. • -r< Then X* = X + Xfc . Now if t > r e T 61 |xF - xF| < K|X{f>0} - x{F>0}| + K|X{f<0} - x{F<0} 1 t r1 — 1 t r ' 1 t r Following the first part of Lemma 4.5.4, XF is S-continuous in T' n {°t > 0} . Now if q > 1 , E|XF - XF|2q < c E| I 1(1 F(x)(Q^ - Qx:|))2 U AtAxd|2q 0<s<r yeX xeX - £ - - -£ + c E| I I K2 u AtAxd|2q " r<s<t yeX Now C £ F(x) CQ- - - Q|~|) )2 < 2 K2 sup (Q| - - <g~-)2 < c — by Lemma 4.3.3. Hence following the second part of the argument in 4.5.4 we obtain the similar result I~F ~FI2q , ,2q/3 E X - X ^ < c(t-r) - t r — - -and thus ~X is S-continuous on T a.s. • Remark: We believe that X is Holder S-continuous of index 2 ~ z ^or any z > 0 . Remark: If we took tf to be hyperfinite, °t = °° , then 4.5.3, 4.5.4 and 4.5.5 will hold a.s. until all finite times, and hence a.s. on ns(T) . We will use this fact without further ado in section 4.7. 62 4.6 Characterization by a Martingale Problem We show now that the measure-valued process x , which we constructed in the last section as L(Xt) o st_1 is in fact the measure diffusion studied by Dawson and others. We will use one of the martingale characterizations given in Roelly Coppoletta (1986), which are refinements of the larger class of martingale problems investigated by Holley and Stroock (1978), Dawson and Kurtz (1982) and others. The characterization of the critical branching measure diffusion x^ given in Roelly-Coppoletta (1986; Theorem 1.3, condition iii)) is as follows: Theorem 4.6.1. x is the (unique) M ( Rd) valued process with weakly t F continuous paths, and the given initial condition such that for all $ e D(A) (the domain of the Laplacian A <j> def „ a) mt = Xt " xo x A$ ds is a continuous martingale. 0 s <mr> t t ,2 x • ds . 0 s We know already (Lemma 4.5.3) that our process is weakly continuous in M C R^) • We will verify a) and b) above for functions <|> e ( R^) only. F b / 2 The extension to D(A) (which is the closure under A of C, „( R ) in b, 2 C (R )) may be done by a limiting argument. It is in fact unnecessary to b 2 d Ms. Roelly-Coppoletta's argument in her Theorem 1.3 to go beyond C ( R ) . b, 2 For such a <{> let $ be a lifting to order 2, and suppose |$| <_ K . Then x<f> = r <M°x)d L(Xt) , t e T' , a.s. °t hs (X) (recall from 4.5.2 that the mass on X \ ns(X) is always infinitesimal a.s.) ° r d ) $(x) U Ax , by Loeb Lifting Theorem (Theorem 2.1.6) X — tx o $ = Xt Hence, <j> 4> <j) mt = xt " xo t A* x ds 0 5 is indistinguishable from the standard part of $ def $ $ r A* M = x - x - J xfAt - 0<s<t -= I $(x)(U - U )Axd - I I u (A$)(x) AxdAt xeX - Ox sx - -(U , A i ~ U ) s+Atx sx = I I *(£) rr AtAxd 0<s<t xeX At d Y l $(x) (A U ) AtAx 0<s<t xeX - -(using the infinitesimal analogue of integration-by-parts) (4-21) y y *(x)[7u AtAxd A HsX Axd]? « . LY - sx 2— sy 0<s<t xeX 1 $ Now is a martingale with respect to the internal filtration F^_ . Hence irr = M Vt e T , a.s. implies that irr is an F -martingale, by °t - " fc t Theorem 5.2 of Hoover and Perkins (1983). This verifies a) of Theorem 4.6.1. Further, since m^ is a.s. continuous in t , the predictable square function <m^> coincides a.s. with the square function [m^] . Now by ^ t another result of Hoover and Perkins (1983) (Theorem 6.7) [st(Mt)] = st[M ] , the standard part of the internal square function, and by Theorem 8.5 of the same paper, 64 CM ]T ~ <M > v t , a.s., $2 $ provided (M^) is S-integrable, where <M > is the internal predictable (with respect to F^_) square process. Now by Burkholder's Inequality, if p > 2 E( sup (M$)P) < c E " t<t £ ~ P " T l $2(x) U AxdAt u •• - sx 0<s<t xeX P/2 < c K E - p \ M At 0<s<t p/2 (where M is the total mass process) s which is finite by Lemma 4.4.1. $ 2 Now apply Chebychev's Inequality to conclude that (M ) is S-integrable as required. Now <M$> 0<s<t xeX (4-22) \ I $ (x) U AtA? 0<s<t xeX I I $2(x) I tr sx d (— — - U )AtAxd v - r 4At, 4 At sx 0<s<t xeX {u < —-i ~ sx A d Ax Now the second term above may be bounded by r, - j. ~v sx , 4Atl 0<s<t xeX — {u < TJ sx A d r I I i -AX { X < L} AtAx Ml2 I I u i /u AxdAt 0<s<t xeX — {u < —-r} { x > L> — Ax < K2 • t • (2L)d • ^ Axd (4-23) + K2 I X^'-L} At 0<s<t -By Lemma 4.5.2, the standard part of the second term in (4-23) goes to 0 a.s. as L -»• 00 . The first term is infinitesimal for all L finite. Hence, for each t 2 ,» . \ <k a.s. 0 ,$ a.s. ° v $ A a.s (4-24) <nr> = <M > = ) X At t t „ L s 11 - 0<s<t -2 f-S2 x ds , 0 s as required by condition b) of Theorem 4.6.1. Hence Theorem 4.6.2. The standard part xfc of the hyperfinite process X constructed by solving the difference equations (4-1) for an internal density U , is the critical branching measure diffusion considered by Dawson and others. 4.7 New Results on the Dawson Measure Valued Diffusion This section represents joint work with Ed Perkins, my thesis supervisor. Now that we know that the measure valued process x^_ which we constructed via solution of the hyperfinite difference equations (4-1) for an internal "density" TJ , is in fact a construction of the Dawson Process, we may use t-X this construction to ascertain some regularity properties of the process. First we need Lemma 4.7.1. Let A c X be a Loeb-measureable set of measure zero. Then on a set fi of probability 1 L(X )(A) =0 ,V t e T* n st_1(t>0) . Proof. Let {A } „ be a nested sequence of internal sets such that A c A n neN - n for each neN, and ° (y (A ) ) ->• 0 . Extend {A } to *N in such a way x n n that the A 's are still nested (that this may be done, follows from n w^-saturation). Let H be infinite in *N . Then u (A ) ss o . Suppose t e T' , x H — with °t > 0 . Now by Lemma 4.3.2, if K is a bound for p (x) , Q-t+At °-H i *€A„ Ax — rl ECX^H) = EC I Axd) = I U Axd < Kyx(AH)M0 = 0 . xeAH — yeX H — -t AH Thus for each t there is a set fi- of probability 1 on which X ~ 0 Let {t }, be a countable S-dense set in T' . Let fi' be the set of -k keN , H A H i probability one on which X is S-continuous for t e T . Let fiTT = n fi n n . Then on fi , (X ) = lim (x. ) = 0 for some k=l %c H t °t Vt ^ -k. - l H \J sequence {t } . Thus X x 0 Vt e T' , hence by infinitesimal underflow —. t ~~ 1 -for each to e ft , there is an infinitesimal e(co) such that rl A For n e *N , let Y = sup X n A 1 . Then Y is S-L for each n , and n teT' t the sequence {Y } is internal. Now E(Y ) x o , for any infinite n ne N H J H e *N . Hence ° (E (Y )) 4- 0 as n =° through N , and thus °Y 4-0 on - n n A a set ft of probability one. Now L(X )(A) <^L(X )(A )= °X n < • °Y V/t e T' , V/n e N • Hence L(X )(A) = 0 ,\/t e T' on ft . • t - A Hence, we may draw Corollary 4.7.2. If A c Rd is a Lebesgue null set (X(A) = 0) , then x (A) = 0 Vt > 0 , a.s. Proof: st_1CA) is a Loeb null set. • We introduce the notations x^ = t f dx. and x = x (A) d f Theorem 4.7.3. If f : R •+ R is a bounded measureable function, then x is a.s. continuous in t > 0 . In particular x^_ is a.s. continuous in t > 0 for any Lebesgue set A . Proof: Let F : X *R be a lifting with respect to L(y ) of f . Then x let A be {x : x e X °F(x) ^ f(°x)} which is a Loeb null set. By Lemma 4.7.1, L(X ) CA) = 0 \/t e T' a.s. Now F is S-L1 of (X,Xfc) V t a.s. (since M is finite Vt a.s.) . Hence xf = f °F(x)dL(X ) = °XF Vt e T' t v- ofc j - t t -F a.s. By Lemma 4.5.5 X is S-continuous a.s. t 68 Theorem 4.7.4. Let <j> be a bounded measureable function on R , and suppose a sequence {A } of uniformly bounded measureable functions k keM converges to 6 in Lebesgue measure. Then for any e > 0 sup Ix - x | —> 0 , He,tf] t t Proof: Let $ and $ k e N , be liftings of <)> and <j> , k e N JC JC respectively. Extend {$, }, to an internal sequence {$ } k e *N . k keN k Then if °t > 0 'E(X k) = ° J u" (*(X) - $ (x))Axd tv tx k — xeX u (<)>(x) - <j> (x))dx where u = st(U) . -i ° x Rd tx -> 0 as k -»• 00 , since u is bounded and L1 ( R) for tx t > 0 . * Hence if He N is infinite, (4-25) °E(X H) = 0 => X H ~ 0 a.s. Vt e T . $ $ H i -1 Now X^ and X^_ -are both S-continuous on T n st (t>0) a.s. by Lemma 5.5.2. Therefore since (4-25) holds for a countable S-dense subset $"$H of T , hence X ~ Q Vt e T' a.s., and def „ H Y„ = max X ^ 0 a.s. H teT' £ Therefore Y^ —> 0 as k ->- 00 through N . 69 Now we have seen (Lemma 4.7.3) that for k £ N (resp. x^_) and k $ st(Xt ) (resp. st(Xt)) are indistinguishable processes. Hence I k 4> I a.s. i ki ° P „ n sup | x^ - xY = o max X = Y —> 0 • t€[e,tfj t t teT' * k 70 CHAPTER 5 The Critical Branching Diffusion in One Dimension 5.1 Introduction In Chapter Three we obtained an existence theorem for SPDEs of the form 3u — = Au + f(u) wt 3t tx where f grows at most linearly at infinity, without the necessity of imposing a Lipshitz condition on f . As discussed in Chapter One, the one dimensional Dawson critical branching diffusion has been believed to satisfy the SPDE ,,. li 32u A /- • (.5-1) — = —- + Zu w dt „ 2 tx dX However the theory of such an equation has not been well-known since the function u -> is not Lipshitz. In this section we use the results of Chapter Three to show that the critical branching diffusion of Dawson does indeed satisfy (5-1) in one spatial dimension, and thus it has a.s. a jointly continuous density. 5.2 The SPDE and the Measure Diffusion TxX Let d = 1 , and let fi be the space {-1, 1} as described in section 4.2. Let U be the solution to the hyperfinite difference tx equations (4-1): U £ tx . tx (5-2) 6 U. = A U + (/rT- A —- -*-) " , t tx - tx tx 2 At v^x-and let X be the measure-valued process, whose internal density is U . In section 4.6 we verified that the process x^ = °(X^) was in fact the critical branching diffusion. In the case d = 1 however, we have the additional information that if U is S-continuous and S-L''" , then U 0 tx is S-continuous on TxX a.s. by Corollary 3.8.3. Furthermore u = °U satisfies (3-12), the weak form of an SPDE; to complete the identification of u as the solution to (5-1), we need only note that, since x 0 , then for u e ns T R) On the other hand, starting from the SPDE (5-1) we note that the 2 martingale problem is easily satisfied since for <f> e C ( R) c t Ao <Kx)dw 0 sx sx R is a martingale, whose increasing process is ft 2 u A (x)dsdx sx R If we also assume that uQ is L^ , then rt o u dsdx < 00 a.s. sx and a dominated convergence argument ensures that the martingale problem is 2 satisfied for <J> e C ( R) (which coincides with D(A) in this case) . b, 2. Thus we have Theorem 5.2.1: Let u_ (x) be continuous and L ( R) . Let w^ be the — 0 tx TxX white noise constructed in section 2.2 from the Loeb space of {-1, 1} onto R+ x R . There is a jointly continuous non-negative process u such that U0x = U0(X) ' and for all 6 e „( R) , and all t e R+ b, z C5-3) utx<t)(x)dx = U0x ^ ^ dx + R 32d) u —- (x) dxds sx „ 2 R dx 0 R /u * (x) dw sx sx Moreover any solution to (.5-3) is the jointly continuous density of the unique measure-valued solution to the martingale problem described in Theorem 4.6.1. Remark: The existence of a jointly continuous density holds for all realizations of the one-dimensional Dawson branching diffusion, not only the one constructed here. This follows from the fact that the existence of a jointly continuous density is a measureable property of the sample paths of this process Csee Cutler (1985)). CHAPTER SIX The Support of the Fleming-Viot Process 6.1 Introduction and Construction As mentioned in Chapter One, the Fleming-Viot process is the limiting case of a model used in theoretical genetics: the Ohta-Kimura model for d quantitative characters. Briefly, in this model, the total number of individuals is conserved, and the dynamics involve two processes, genetic "drift", and mutation. Mutation is modelled by a random walk on Z "Drift" is modelled by replacing individuals at random by new individuals whose genetic type matches another individual chosen at random from the rest of the population. If we denote the types by points k e Zd , then we may denote by p(t,k) , the number of individuals of type k alive at time t e R+ , divided by N , the total number of individuals (conserved), then pCt,-) forms a continuous-time countable state space Markov jump process with generator: C6-1) Lffp). = I l^je Z [y p(i)p(j) + D p(i) 9i.](f(p1D) - f (p)) where p1"1 (k) p(k) + i- , if k = j { p(k) - , if k = i P(k) otherwise, and f 1, if = 1 e -2d, if = 3 0, otherwise, and 74 y and D are positive constants describing the rates of "drift" and mutation respectively. Each of the N individuals takes one step ("mutates") in one of the 2d possible directions at Poisson times whose rate is D/N . Similarily each individual dies and is replaced by another whose type coincides with that of a given other of the N-1 individuals (i.e. an "offspring" of that 9 individual) according to a Poisson process with rate Y/N To construct the Fleming-Viot process, we re-scale time and space by 1/N2 and 1/N1/2 respectively: for A e B( Rd) let (6-2) XA = T p(N2t, k) — eA takes values in M^( Rd) , the space of probability measures on Rd . The result of Fleming and Viot (1979) is that in the limit as N •+ 00 , the measure-valued process defined by (6-2) converges in the space of (R ) -valued processes provided the initial measure converges in M ( Rd) . We refer to Fleming and Viot (1979) and Dawson and Hochberg (1982) for more details on this. However, a non-standard construction of the Fleming-Viot process is immediate from (6-2). We simply take N infinite and let the state space be *Z . Stated in non-standard language the result of Fleming and Viot is * + Theorem 6.1.1. For N infinite, and t e n.s.( R ) (6-2) defines a.s. a process with S-continuous paths in ns( ( R )) , if initially all but an * d r-infinitesimal fraction of the N ' particles are on nearstandard points in Z /vN . Note that this result asserts that although the parameter space in (6-2) is * + -1 + technically R , that X ~ X"_ in the weak topology whenever t ~ r e st (R ) , a.s. Hence X has a standard part, x , which is a measure-valued Strong Markov process called the Fleming-Viot Process. 75 6.2 The Dimension of a Putative Support Set We use here the terminology of Dawson and Hochberg (1982) and some of their results. What we aim to show in this chapter is that a non-standard construction makes much of their work more natural (and a good deal easier'.) as well as extending their results. The main result (Theorem 6.4.4) of this chapter asserts that the dimension of the support of the Fleming-Viot process is at most 2 for all times simultaneously, a.s. (Dawson and Hochberg were able to show this only for fixed times). In this section we derive the dimension of a set, that in section 6.4 will be shown to be the supporting set. Consider a particle (or individual) alive at some time t . As the process evolves and the particle wanders, at some time r > t this particle may disappear (to be replaced by another particle somewhere else), or else it may serve as the 'type-model' for the replacement of some other particle which disappears at time r . In this latter case we say that both particles at time r are 'descendants' of the original particle at time t . For ease of terminology we will say that the particle at any time s > t, up until the time of disappearance of that particle, is the descendant of that particle at time t . Note that ancestry is a transitive relation. Furthermore every particle at time t has a unique ancestor at any time s < t ; if we follow the paths of particles backward in time, they may converge, but they will never split. Two particles at a time r > t are said to have a common ancestor at time t , if they are both descendants of a given particle at time t . We will construct a supporting set for the mass of N particles at any time, by looking for a small (finite) set of ancestors, at an earlier time whose descendants comprise all of the N particles at the later time. 76 Consider any finite time interval [0, T] . Let e > 0 , and let {a^} and ^A^} be strictly decreasing sequences of positive numbers such a2/A that A e n n -> oo and n a2+E/A 0 as n -* 00 . (A convenient choice would be a = A (4+e)/ (4+2e) ^ ^ n n n We suppose w.l.o.g. that {T/A^} are all integers so that for each n , the (finite) sets {t |t = k A ,0<k< T/A } form a partition of n t JC n / Jc n — — n CO, T] . We let N be the number of ancestors at time t of the n,k n,k-l N particles alive at time t , . Let A . be the union of balls of n,k n,k radius a , centered at each of the N , ancestors at time t , , of n n,k n,k-l the system of N particles at time t . Let A be the union of n,k n,k a. smaller balls of radius —— , centered at the same points. (Technically all of the above are non-standard (internal) objects, but since they are all near-standard I will make the distinction between standard and internal only when necessary.) Let k (t) = Ct/A ] identify in which interval Ct , , t , ,), t n n n,k n,k+l lies. Then let (6-3) A. = u n st(A ) t i n,k (t) m=l n=m n which is a (standard) Borel set for each standard t in Co, T] . This will turn out to be our supporting set. The dimension of this set depends on the number of ancestors at time t , , of the N particles at the times t , (and n,k * n,k hence of all the particles alive between times t , and t , ) . n,k n,k Dawson and Hochberg treated this problem for fixed times by constructing an infinite particle system to describe the Fleming-Viot process at that particular time. Using our hyperfinite model, we can use the same system of N particles at all times. Dawson and Hochberg (1982; (6.23)) showed that N the distribution of the time T taken (in reverse time) to reduce the m number of ancestors of N particles to m only, had Laplace transform -STN , (6-4) E(e = H? . (1 + -f-T,) • k=m+l yk(k-l) (Their argument applies verbatim to our hyperfinite scheme.) We may use (6-4) to estimate P (N >_£.) = P(°TN/a > °A ) for integral -2- . From (6-3) n,k A c/A n A n n n (6-5) E(TN/A ) = I } =J^- 4~^,and - c/A L yk(k+l) Yn Yc n k= + 1 A A n n N N H A C/A_ ~ r..,.„. ,n2 k= A n n , c + 1 [Yk(k-l)] V 1 r 1 l-i 2 3 3' k- -£+ 1 3Y (k"D k A 1 n 2,c3~ 3J~ 23 Y ) N 3y c A 3/2 oA °A Now by Chebychev's Inequality, P( T^^ > + h ——j—) < — , n Y f3yc h hence, tak ing h = °/5c (yc-1)/SE~ , n (6-6) P(N , > T5-) < An (for c > Y"1) . n,k A -n 3C(YC-1) 78 2 + E Let 6 > 0 . Let y£(A) denote the Hausdorff x measure of a set Then PC max N , > 6/a2+£] k<T/A n'k — n < P(N > 6/a2+e) A n,l ' n ' n A < f- 2__ _ (FROM (6_6)) n n , n . > 3 ~2+7 (Y ~2+F " } a a n n (6-7) = A A 2 3 6 (Y6 ~2+F ~ 1} a a n n 0 as n -> 00 since A n 2+e a n 2+e Hence P( max a N , ,,>><5)">"0» ancj by taking a subsequence, if t£[0,T] n n'kn(t) necessary, we can ensure 0+c (6-8) max a N , #J_, -> 0 a.s. t£[0,T] n n'kn(t) oo Now for each m , A , ... is a covering of n A , . by N , fJ_. m,k (t) n,k (t) 1 m,k (t) m n=m n m balls of radius a . Hence (6-8) ensures that m y ( n A , ,,,) = 0 uniformly in t , a.s. z m n,k (t) 00 00 and hence y ( u n A , ,,..)= 0, for all t e [0, T] a.s. This is true e m=l m n-kn(t) 79 for any E > 0 . Lemma 6.2.1: The random set-valued frunction A defined by (6-3) has Hausdorff dimension at most 2 for all time, a.s. 6.3 A Useful Stochastic Differential Equation We study the numbers of descendants of a specified group of particles. At any time, we stop the Fleming-Viot process, designate n particles out of the N , and restart this (Markov) process. For this section w.l.o.g. suppose the process is re-started at time 0 . Let Y denote the mass associated with the designated n particles at time 0 and with their descendants at times t > 0 . Y = — (# of particles alive at time t t N which are descended from the original n) . Recall that we are including particles in the original n that have remained unchanged until time t , in this count. Now Y^_ changes only when i) one particle out of the NYfc disappears, and is replaced by a particle of a type-model not included in the NY^_ , or ii) one particle out of the N-NY^_ excluded particles disappears and is replaced according to the type of one of the NY^_ descendant particles. Disappearances and replacements happening entirely within the context of the NY^ designated particles or amongst the N-NY excluded particles, make no change to Y^_ . Now the replacement of any given particle according to the type model of any other given particle happens according to a Poisson process with rate y (see (6-1)). There are NY • (N-NY^_) possible ways for an event of type i) to occur each happening according to a Poisson process with rate y , and each causing a decrease of size — in Y. . There are (N-NY ) • NY possible ways for N t t t an event of type ii) to occur, each with the effect of increasing Y by and occurring at the same rate y as an event of type i). Thus Y is an (internal) martingale. The associated predictable increasing process is easy to compute, since the change to Y at any time is the sum of 2NYt(N-NY^) independent Poisson processes, each of rate y , whose predictable quadratic variations are each —j . Thus N <Y> = t 2y Y (1-Y )dt 0 * (6-9) d<Y>t = 2y Yt(l-Yfc)dt Hence by Hoover and Perkins (1983)"'" we may conclude that the paths of * o Y are a.s. nearstandard in the space C((0,°°); R) . Let y = Y . Then ,t from (6-9) d<Y>t = 2Y Yt(l-y )dt and hence b^_ = dy is o /2Y ytd-yt) t a standard Brownian motion (up till the time of extinction of y) . Thus by enlarging our probability space we may find a Brownian motion b such that (6-10) dyt = /2y yt(l-yt> dbt Furthermore we have the following lemma. Lemma 6.3.1: Suppose y = e where 0 < e < — . Then there are finite 0 2 1 * Note. Hoover and Perkins (1983) Theorem 8.5 refers to -discrete time processes X , such that if Ax is the change in X over an infinitesimal time step At , then sup °|Axj =0 a.s. By looking at our process Y only at intervals of At = 1/N2 , we may ensure sup °~|AY| = 0 a.s. and bring Y into the framework of the theorem quoted. 82 constants K^, K2 independent of e , such that P( 3s e [0, t] such that y < ^ e) < K. exp C- r^-] . s — 4 — 1 K^t Proof: Write yfc as the time change of a Brownian motion . As long as 3 5 Yfc stays in the range [- e, - E] the derivative of the time change must be at least — ^—^— • Now use the estimate -2yeU-4e) P( 3 s e [0, t] such that lb I > c) < 2 P(|b I > c) s — ' t1 — 2 c = 4 P(b >_ c) < K e 2t • 83 6.4 Verification of Support We now check that the random set Afc of section 6.2 does indeed support the measure xfc of Theorem 6.1.1 (the Fleming-Viot Process). Lemma 6.4.1: Consider at time t any particle p . Let s < t , and trace the path of the (unique) ancestor of p at each time r , for s <_ r < t . Then this path is an infinitesimal random walk whose standard part is a d-dimensional Brownian motion of rate 2D . Proof: Between the appearance and ultimate disappearance of any given particle, it takes a step of size — from its current position to any one /N of the 2d neighbouring positions according to a Poisson process with rate DN . Replacement (which is bifurcation of an ancestor particle) occurs according to a Poisson process which is independent of the motion of the particle. If we imagine such a motion continuing indefinitely and call this process B^_ , and its coordinates B^ , then clearly each B^ is an internal martingale, since steps to the right occur at the same rate as steps to the 1 2 left. Now <B1> = 2DNt( — ) = 2Dt , since B1 is the sum of 2 independent 1 /N Poisson processes, of rate 2DN , and of amplitude —— . The motions B1 and /N ~B? happen according to independent Poisson processes, hence <B1, B"5>t = 0 if i ^ j . Thus by Hoover and Perkins (1983) (once againl) B^ has a standard part b^_ a.s., and this b^ satisfies the characterization of d-dimensional Brownian motion: E(bb 0<r<s)=b ; <b> = 2DIt. • t' r — — s t 2 Again we must look at the process Bt at discrete intervals (of size At = 1/N) to put it in the framework of Theorem 8.5 of this paper. Referring to 6.2 for the definitions of the sets A , and A , . n,k n,k c Kl ~an/K2A Lemma 6.4.2: P(x ((°A ,))>£)< — °e n for finite constants ° n,k — e n,k V K2 • Proof: Consider the family of descendants of any one of the N particles n / JC at time t . . used in defining A If we trace back from time t , n,k-l 3 n,k n,k the movements of any one of the particles in this family, and its' progenitors, until we come to the position of the one ancestral particle at time t , , , n,k-l we find a motion of the kind described in Lemma 6.4.1, whose standard part is a Brownian Motion. Hence the displacement of any given particle at time t n,k from its' ancestor at time t ' , is distributed N(0, 2DA i) . For anv n,k-l n J particle p at time t , , n,k a„ -a2/K A E(IAC (P)) £P(BDA > -) < K e A . n n,k Hence N -a2/K A E(x U°A f)\< ECX^ CAC))= \ I EI (p.) < °K e n 2n °, n,k — t . n,k N . - ;o I — 1 t . n,k 1=1 A , n,k k and c K-, „ -a2/KnA PCx (C°A V)C) > e) < -i °e n 2 n . o n,k e fcn,k c K o -VVn Lemma 6.4.3: P(x ((A ,))>£)<— e o n,k e n,k+l Proof: As above, with larger K2 , since displacements are distributed N(0, 4D A I) . n Let E > 0 . From Lemmas 6.4.2 and 6.4.3 we may deduce P( max {x (°AC , ) v (°A° )} > T fcn,k n'k Si.k+l n'k 2 n K T -a2/K A < — enzn-vO as n <*> . - EA n Cas per usual, our trusty servants, the constants and K2 are changing values when necessary). Hence a fortiori PC max (x (°AC ) v x+. (°AC )} > f) -> 0 . n Now consider the possibility that the sets A , contain 1 - — of , n,k 2 the mass of the process x at times t , and t , , , but that more than t n,k n,k+l O E of mass lies outside the sets A , at some time in between t , and n,k n,k t , , - Let s , be the first such time, s is a stopping time for n,k+l n,k n,k the Markov process x^ and we may consider it restarted at time s , from t k its configuration at s , . One of the following must occur during the n,k interval Cs t . ,.) • n,k n,k+l a) at least one fourth of the mass e that lies outside A , initially n ,k a travels a distance — Cto re-enter A , ) 2 n,k b) the mass e that lies outside A , decreases by at least one fourth n,k e 3 Cthat is the — N particles have at most e N descendants at time t n,k+l) Consider case a). For each particle p at time t , , , the displacement from its n,k+l ancestor at time s , is distributed N(0, D(t , , - s )I) . Recall n,k n,k+l n,k from section 6.3 that the number of descendants at a future time of any given subset of particles, forms a martingale. Hence the expected mass at time t • , of descendants of the e N particles at time s , , which are n,k+l n n,k 3. outside balls of radius — centered on their ancestors, is equal to the initial mass e , times the probability that any one particle is outside such -a2/K2A a ball, which product is bounded by e K e n (as in Lemma 6.4.2). Hence using Chebychev's Inequality, "an/K2An P (for a fixed n and k, case a) occurs) < K^e T K1T "an/K2An Hence P(for some k < -— , s , < t , and case a) occurs) < —:— e A n,k n,k+l A n n which goes to zero as n -*• 00 . Now by Lemma 6.3.1 T P(for some k < — , s , < t , and case b) occurs) A n,k n,k+l n < — P ( 3 s < t , -s <A s.t. y = - E I v = e) — A - n,k+l n,k n ys 4 1 Y0 b' n (where y is the process mentioned in.6.3.1) -e/K.A T 2 n - A~K1 6 n which also goes to zero. Hence P( 3k < T^- such that x^ (°A °, ) < ^- and ( A C, ) < A t n,k 2 t n,k 2 n n,k n,k+l and 3 s e (t . , t ) such that x (°A C, ) > e) ->- 0 . n,k n,k+l s n,k 87 » C P Thus sup x ( A , ) —> 0 as n + te[0,T] t n'kn(t) By taking a subsequence, if necessary, we may ensure ° c a • s sup x ( A . .) —'—$ 0 . Again by taking a subsequence we may assert te[0,T] r n,Knm V sup x (°AC ) < nil tcCO.T] t n'kn(t) 00 a.s. Then sup x (A^) te[0,T] t o c <_ sup x ( u A . , . ) for any m e N , te[0,T] Z n=m n'Kn^> <_ sup I x (°A C ) te[0,T] n=m t n'kn(t) < y sup x (°A C ~TL t£[0,T] * n'kn(t) -> 0 as m -> 00 Thus we have Theorem 6.4.4: The random measures xfc , which are the realizations of the Fleming-Viot process, are supported for all times t on a set A of dimension at most 2 , a.s. APPENDIX A Some Inequalities Used in Chapter 3 A.l Purpose In Appendix A.3, we derive some inequalities involving the coefficients g Q- which were presented in Chapter 3.4. In section A.2 we derive some identities which provide a neat route to the inequalities in A.3. A.2 Some Identities s y y As we saw in Lemma 3.3.1, Q- = p(B- , = x where B- is an infinite-x-y - s-At simal random walk starting at y at time 0 , and thereafter taking steps of size Ax to the right or left with probability a , at each time interval 2 1 1 At , where a = At/Ax <_ — . The reason for not considering — < a <_ — will become clear later. As we saw in Lemma 4.3.1 the above remarks are true also in d-dimensions, where the random walk takes steps in any of 2d directions, each with 2 probability a = At/A^ . The identities in A.2 are true in d-dimensions although the d-dimensional versions will not be used The first identity is trivial. s Lemma A.2.1. 7 Q- = l,if —- £ N L x At x Ax e *Z 2 Lemma A.2.2. £ (Q|) = Q^-~At > if s/At £ *N x * ~ £ Z Ax 89 Proof: I (Q|)2 = ( I Q|)2 - I I Qx 9? x - x x y^x - -x - y^x (by A.l) = 1-^(Bs-At = ^ '^s-At^ x 1 - I ^_At = x) -^B^^-x) X (by symmetry; this step fails if we consider a reflecting random walk.) 1 - I £<B°_At = x) P(B°s.2At ^0|B°_At=x) x 1 - P(B° . f< 0) - 2s-2At 28-At y0 (by definition). • * * Lemma A.2.3. If z/Ax e Z , and s/At e N x/Axe Z - - -Proof: 7 (Q- - Q- ) ,.L * X x+z x/Axe Z - - -2 2 = I (Q-) - 2 y Q- Q- + I (Q§ ) ^ x ^ x x+z L x+z x- x - - - x - -2 - 2 l P(B°.4T - ,,P(B°.TT . ; • (by Lemma A.2.2) X (by symmetry) = 2 Q2S~At - 2 P(B° OA = z) *0 - 2s-2At = 2(2^-2^ u z r s Lemma A. 2.4. If — , — e N , J (Q£+S - Qf)2 = Q2£+2rAt + 2s-At _ f+2s-At x/Axe*Z x x 0 0 0 Proof: Y (Q-+- — Q-)2 v x x Ax x- x-- - x-2r+2,-At _ 2 l p(B0 = 0 = 0|B0 0 x - r+s-At - - r+2s-2At £+2-^ (by symmetry and A.2.2) 2r+2s-At ,„0 „. „2s-At 2 " " - 2 ^CBr+2s-2At = 0) + V A.3 Some Inequalities In this section I prove the four inequalities 3.4.1 through 3.4.4 of Chapter 3. Lemma A.3.1. (3.4.1). There is a constant K such that t I (QX) < K At/t , if ^ xeX * e N Proof: Clearly xeX x x/Axe Z x 0 By definition, ofk+1*At = P (IS, I < zr) , where S is 0 k — 2 k the sum of k I.I.D. random variables taking the values -1, 0, +1 with probabilities a , l-2a , a respectively. Var(S ) = 2ka . k Since °a > 0 , then by Corollary 2.2.3 of Bhattacharya and Rao (1976), (k+l)At -0 (2 ' /2kct 1 /2TT e dz "2 /2ka < — as k -> Now 2 ' /2ka 1 -z /2, e dz is asymptotically 1 x , 1, + o(—) as k •> 00 2/2irka /k •/2ka Thus there is a finite K such that (k+l)At < JC 0 "A ft By the transfer principle this must hold for all of N . Thus (checking t = At separately) 2t-At < _J±_ < At y0 - 2t-2At - t * Lemma A.3.2 (3.4.2). There is a constant K such that for t/At e N I I (Q-)2 1 K/t/At . 0<s<t xeX -n 1 Proof: Follows from A.3.1 and the fact that £ —f_ C Sn for all n e N , k=l /k * hence for all n e N . • 5 * Lemma A.3.3 (3.4.3). There is a finite constant K such that if -— e Z , I I Ax t * r t-ss^ I ? I and - e N , then ) ) (Q- - Q- ) < K. At 0<s<t X£A - -Proof: Suppose w.l.o.g. that z ^_ 0 . Let Bg be the infinitesimal random walk whose density is Q-+A^ • x Let J(t,x) = I I(Bs > x) (Bs+A " B ) . 0<s<t -Let L(t,x) = I I (B = x)Ax . 0<s<t -Let L(t,x) = I KB = x ; B = x + Ax)Ax - - „ . S - S+At -0<s<t Then L is the true occupation density ("local time"), for our random walk with 'pauses', while L is the discrete analogue of Brownian Local Time. Now E(L(t,x)) = 7 P(B = x ; B . - B = Ax)Ax --- ^ " - s - s+At s 0<s<t - -93 Y P(B = x)a • Ax u s -0<s<t = a E(L(t,x)) . In Perkins (1982) Lemma 3.1, Tanaka's formula is established by internal induction on t : (Bfc - x)+ = J(t,x) + L(t,x) Thus L(t,0) - L(t,x) = - (B - x)+ - J(t,0) + J(t,x) Now J(t,0) and J(t,x) are both internal martingales, hence E(L(t,0) - L(t,x)) = E[Bp - E[(Bt - x)+] (A.3.1) Now since 0 < °a < | , Q§ = V (Bs_At = 0) >_ P (Bs_At = x) = Q| , as shown in Lemma 4.3.2. Now by A.2.3 I 1 (Q- - Q~ )2 < 2 I Qn22"At - Q2s"At _ . v x x+z — „ u . 0 z 0<s<t xeX - - - 0<s<t < * i <er4t - c4t» At<s<2t U -= 2 I P(B = 0) - P(B = z) 0<s<2t-2At - -1 2 E(L(2t-At, 0) - L(2t-At, z)) |- E(L(2t-At,0) - L(2t-At,z) < - , by (A.3.1) ~~ a Ax Lemma A.3.4 (3.4.4). There is a finite constant K such that if r/At < t/At are in *N , 2 I I (Q!~- - Qr~~> < K/(t-r)/At 0<s<r xeX Proof: The l.h.s. above is equal to V r . s+(t-r) s.2 L L (Q~ ~ ~ ~ Q~' i which by A.2.4 is bounded by 0<s<r xeX x -(A-2) I (Q2s+2(t-r)-At + Q2s-At _ 2 Q(t-r)+2s-At) 0<s<r 0 0 First, suppose t-r < r and (t-r)/At is even. Then many of the terms in (A-2) cancel leaving t-r t-r 0<s<—— <s<t-r --2 2 -I Q2-"AT + I Q2S_AT t-r ° t-r 0 r<s<r+—— r+ — <s<t - - - 2 - 2 ° < 1 s Since a — - , the coefficients Q- are monotone decreasing as 3 0 s increases. Hence the sum of the second, third, and fourth terms in (A-3) above is bounded above by 0 . The first term in (A-3) is bounded by KV (t-r)/At by A.3.2, as required. t-r Now if t-r < r and ~~7~7~ is odd the sum in (A-2) is bounded by Y rn2s+2(t-r)-At n2s-At - (t-r)+2s1 , . . I LQ + Q - - 2 Q -J .. Applying the same 0<s<r v 2s~At i cancellation argument, this is bounded by \ Q - < K /t-r/it t-r 0<s<-— +At — 2 Now if t-r > r (i.e. r < t/2) , then I I <<^> - Q|)2 0<s<r xeX x x 0<s<r xeX y 2s+2(t-r)-At + Q2s-At ~ 0 0~ 0<s<r (by A.2.2) 2s-At < 2 I Q — L o 0<s<r — K AT (BY A,3-2) t-r < K —— . - At 1 Remark: If a = — then Lemmas A.3.3 and A.3.4 are false. The pattern of non-zero coefficients in the array {Q-} , s/At e *N , x/Ax e *Z , is a checkerboard pattern. Thus U is independent of U, , and of any U £x t-At,x J s; t-s x-y for which —-— + ~ is an odd hyper-integer. Thus the moment inequalities At Ax 96 on spatial and temporal differences will fail. This is the reason for not using the simplest finite difference scheme in Chapter 3. APPENDIX B Internal Solutions to SPDEs in Higher Dimensions The success of the hyperfinite difference equation approach in Chapters 3 and 4 to the existence of solutions to SPDEs in one dimension, and to the Dawson measure diffusion in higher dimensions, leads one to wonder if the use of hyperfinite difference equations might lead to a general theory of SPDEs in higher dimensions. So far at least, this hope has not borne fruit. The kind of equations that we would be led to consider after the analogy of those in Chapters 3 and 4, would be of the form (B-l) !^ = Au + f (u)dWfcx , t e R+ , x e Rd and the corresponding internal equation t'B-2) (6 U ) = (A U. ) + F(U ) 5 / AtAxd , for t -x t t- x tx tx t e T and x e X , hyperfinite grids representing R and R respectively. Now it is easy to see that an internal solution to (B-2) exists, by the usual inductive construction. What is not so clear is whether or not this internal solution U has a non-trivial standard part u , presumeably in some space of distributions. It is also not clear what it would mean for such a distribution-valued process u to be a solution of (B-l), in general, since non-linear operations on point values of distributions are undefined or discontinuous at best. It seems possible to make sense of (B-l) if f(u) is an operator valued function of u with values in a class of operators on a space of distributions, but to pursue this possibility would 98 take us too far afield from the ideas of this thesis. If we restrict ourselves to the case where f is a real valued function of a real variable, and that F is some natural lifting of f , such as f , then we may still ask whether or not the internal solution to (B-2) has an interesting standard part. Sadly, in several cases the answer seems to be 'no'. The first case we considered is when f is a continuous function of compact support, (and F is an S-continuous lifting of compact support). Then since At << Ax , changes to the internal solution U to (B-2) are infinitesimal at each step and an easy induction argument shows that the internal solution U is bounded by a finite constant. If we seek to tx estimate the variance of £ U Ax*3 for some finite internal rectangle xeA — (or other set) A , then we are naturally led to examine E(F (U )) . We may sy set up a difference equation for this quantity, and show that it is everywhere a.s. infinitesimal. Thus the variance of the integral over any finite region of U , is infinitesimal. Thus though the values taken by U lie almost always at either end of the (connected) support of F , these values balance very precisely on each monad. The balance point is infinitesimally close to the value of the deterministic solution to the heat equation with the given initial condition, at the standard point corresponding to the monad in question. The other cases which we examined were when F(U) = , when either 0<p<i,or ^<P^_1- We maY try to follow the development of section 4.4 to find the total mass. In the case 0 < p < - , we find that 2 E(M -Mn)2q is finite, at least in the case of reflecting boundary conditions - t u on a finite rectangle, simply by bounding U by 1 + U . However, a closer examination of the difference equations (B-2), will show that, . the values U are almost always infinite or infinitesimal. Thus the quantity £ U2p is actually infinitesimally __ smaller than £ , since almost all the mass comes from points x xeX where U is infinite. Thus we end up with, in fact E|Mt - M0|2q << E( 1 M^) , and thus the total mass is infinitesimally 0<s<t -close to MQ a.s. If we examine the predictable increasing process associated with £ (U - U )$(x) - £ xA$ At , as in section 4.5 we find X "cx ux — ^ . s 0<s<t = it to be £ J U2p $2(x)Axd , which is likewise infinitesimally close 0<s<t xeX --to 0 . Thus the process we get from solving (B-2) for F(U) = iP , 0 < p < turns out to have a deterministic standard part. If F(U) = UP , i < p <_ 1 , then, again U must almost always be infinite or infinitesimal. If we obtain an estimate of the variance of the total mass M , we find E(M - M )2 = Y Y U2p AtAxd . 0<s<t xeX — Whenever U is infinite U << U2p , and hence Y U2p Axd >> Y U Axd Lv sx Lv sx xeA — xeA — The moment bounds on cannot be obtained. An examination of the predictable ,* y XA$ ^ s 0<s<t -that this is a.s. infinite, and thus the quadratic variation is simply too larg increasing process for X* - X* - £ XM At as in section 4.5, shows t 0 S for U to be nearstandard as a measure-valued stochastic process. The only cases where the formulation (B-l) appears to yield anything are if a) F is essentially like F(x) = C , in which case we recover the linear theory of SPDEs in higher dimensions which was originally developed by Walsh (1986), or b) F(x) = c/x essentially, which is treated in Chapter 4. 100 REFERENCES R. M. Anderson, "A Non-Standard Representation for Brownian Motion and Ito Integration", Israel J. Math., 25, 15-46 (1976). K. B. Athreya and P. E. Ney, Branching Processes, Springer Verlag, Heidelberg, 1972. R. N. Bhattacharya and R. R. Rao, Normal Approximations and Asymptotic Expansions. John Wiley & Sons, New York, 1976. D. L. Burkholder, "Distribution Function Inequalities for Martingales", Ann. Prob., 1, 19-42 (1973). C. Cutler, "Some measure-theoretic and topological results for measure-valued and set-valued processes", Technical Report #49, Laboratory for Research in Statistics and Probability, Carleton U., (1985). N. J. Cutland, "Nonstandard Measure Theory and its Applications", Bull. London Math. Soc, 15, 529-589, (1983). D. A. Dawson, "Stochastic Evolution Equations". Math. Biosciences 15, 287-316, (1972). D. A. Dawson, "Stochastic Evolution Equations and Related Measure Processes". J. Multiv. Anal., 3, 1-52, (1975). D. A. Dawson and K. J. Hochberg, "The Carrying Dimension of a Stochastic Measure Diffusion", Ann. Prob., 7, 693-703, (1979). D. A. Dawson and K. J. Hochberg, "Wandering Random Measures in the Fleming-Viot model" , Ann. Prob., 10, 554-580, (1982). D. A. Dawson, and T. G. Kurtz. "Application of Duality to Measure-valued Processes", pp. 91-105 in Lecture Notes in Control and In. Sci. 42 Information Sciences, ed. W. Fleming and L. C. Gorostiza, Springer-Verlag (1982). D. A. Dawson, "Measure-Valued Processes: Construction, Qualitative Behaviour, and Stochastic Geometry", Technical Report #53, Laboratory for Research in Statistics and Probability, Carleton U., (1985). W. H. Fleming, and M. Viot, "Some Measure-Valued Markov Processes in Population Genetics Theory", Indiana Univ. Math. J., 28, 817-843 (1979). R. A. Holley, and D. W. Stroock, "Generalized Ornstein-Uhlenbeck Processes and Infinite Particle Branching Brownian Motions", Publ. R.I.M.S., Kyoto Univ., 14, 741-788 (1978). D. N. Hoover, and E. Perkins, "Non-Standard Construction of the Stochastic Integral and Applications to Stochastic Differential Equations", parts I and II, Trans. Amer. Math. Soc. 275, pp. 1-58, (1983). 101 H. J. Keisler, "An Infinitesimal Approach to Stochastic Analysis", Memoirs of the A.M.S., 297, (1984)* E. Perkins, "Weak Invariance Principles for Local Time", Z. Wahr. verw. Geb., 60, 437-451, (1982). S. Roelly-Coppoletta, "A Criterion of Convergence of Measure-Valued Processes: Application to Measure Branching Processes", Stochastics, 17, 43^65, (1986). A. Stoll, "A Non-Standard Construction of Levy Brownian Motion", Technical Report, Universitat Bochum (1984). A. Stoll, "A Nonstandard Construction of Levy Brownian Motion", Probab. Th. Rel. Fields, 71, 321-334 (1986). J. B. Walsh, "A Stochastic Model of Neural Response", Adv. Appl. Prob., 13, 231-281, (1981). J. B. Walsh, "Regularity Properties of a Stochastic Partial Differential Equation", pp. 257-290 in Seminar on Stochastic Process, Birkhauser, Boston, (1984) . J. B. Walsh, "Martingales with a multidimensional parameter and stochastic integrals in the plane", in Ecole d'Ete de St. Flour, Lecture Notes in Math. v. 1180, Springer-Verlag (1986). S. Watanabe, "A Limit Theory of Branching Processes and Continuous State Branching Processes", J. Math. Kyoto University, 8, 141-167, (1968). M. Jirina, "Stochastic Branching Processes with Continuous State Space", Czechoslovak Math. J., 8, (1958).
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Hyper-finite methods for multi-dimensional stochastic...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Hyper-finite methods for multi-dimensional stochastic processes Reimers, Mark Allan 1986-08-18
pdf
Page Metadata
Item Metadata
Title | Hyper-finite methods for multi-dimensional stochastic processes |
Creator |
Reimers, Mark Allan |
Publisher | University of British Columbia |
Date Issued | 1986 |
Description | In this thesis we introduce Non-Standard Methods, in particular the use of hyperfinite difference equations, to the study of space-time random processes. We obtain a new existence theorem in the spirit of Keisler (1984) for the one dimensional heat equation forced non-linearly by white noise. We obtain several new results on the sample path properties of the Critical Branching Measure Diffusion, and show that in one dimension it has a density which satisfies a non-linearly forced heat equation. We also obtain results on the dimension of the support of the Fleming-Viot Process. |
Subject |
Finite differences Stochastic processes |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2010-08-18 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0079656 |
URI | http://hdl.handle.net/2429/27515 |
Degree |
Doctor of Philosophy - PhD |
Program |
Mathematics |
Affiliation |
Science, Faculty of Mathematics, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-UBC_1987_A1 R44_6.pdf [ 4.15MB ]
- Metadata
- JSON: 831-1.0079656.json
- JSON-LD: 831-1.0079656-ld.json
- RDF/XML (Pretty): 831-1.0079656-rdf.xml
- RDF/JSON: 831-1.0079656-rdf.json
- Turtle: 831-1.0079656-turtle.txt
- N-Triples: 831-1.0079656-rdf-ntriples.txt
- Original Record: 831-1.0079656-source.json
- Full Text
- 831-1.0079656-fulltext.txt
- Citation
- 831-1.0079656.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0079656/manifest