<ph f="cmr"> </ph><ph f="cmbx">A superprocess involving both branching and coalescing</ph>

Xiaowen Zhou1 1 E-mail address: zhou@alcor.concordia.ca

Department of Mathematics and Statistics, Concordia University, 7141 Sherbrooke St. W., Montreal, H4B 1R6, Canada

Keywords: coalescing simple random walk; coalescing Brownian motion, duality; super process with coalescing Brownian spatial motion; Arratia flow; Laplace functional; Feller's branching process; Bessel process AMS Subject Classification: 60G57; 60J65; 60G80; 60K35

1 Introduction

In this paper we mainly consider the following branching-coalescing particle system which can be described intuitively as follows. A collection of particles with masses execute coalescing Brownian motions. In the mean while the masses for these particles evolve according to independent Feller's branching processes. Upon coalescing those particles involved merge into one particle with their respective masses added up.
The above-mentioned particle system can be described using a measure-valued process Z   .
More precisely, the support of Z t   represents the locations of those particles at time t   , and the measure Z t   assigns to each supporting point stands for the mass for the corresponding particle.
This processes Z   , referred as the superprocess with coalescing Brownian spatial motion (SCSM), was first introduced in [5] . It arises as a scaling limit of another measure-valued process, which was referred in [3] as the superprocess with dependent spatial motion (SDSM). As to SDSM, it arises as a high density limit of a critical branching particle system in which the motion of each particle is subjected to both an independent Brownian motion and a common white noise applied to all the particles. More precisely, the movement of the i th   particle is governed by equation d x i ( t ) = σ ( x i ( t ) ) d B i ( t ) + ∫ R h ( y − x i ( t ) ) W ( d y , d t ) ,   where ( B i )   is a collection of independent Brownian motions which is independent of the white noise W   ; see [3] . A similar model was also studied in [12] .
It was shown in Theorem 4.2 of [5] that, after appropriate time-space scaling, SDSM converges weakly to SCSM. A functional dual for SCSM was given in Theorem 3.4 of [5] . In addition, using coalescing Brownian motions and excursions for Feller's branching process, a construction of SDSM was found in [5] , an idea that initially came from [4] . In this paper we always denote such a SCSM as Z   .
One of the most interesting problems in the study of a measure-valued process is to recover a certain dual relationship concerning the measure-valued process. Such a dual relationship often leads to the uniqueness of the measure-valued process; see [10] for some classical examples on super Brownian motion and related processes. It is not hard to show the existence of Z   as a high density limit of the branching-coalescing particle system. The main goal of this paper is to propose a new way of characterizing the measure-valued process Z   via duality, in which the self duality for coalescing Brownian motions plays a key role. To this end, we first prove a rather general duality on two coalescing Brownian motions running in the opposite directions. We derive this duality from an analogous, essentially combinatorial, fact about coalescing simple random walk. With this duality we can express certain Laplace functionals for Z   in terms of systems of coalescing Brownian motions.
We could carry out some explicit computation thanks to the above-mentioned duality. In particular, we first show that, starting with a possibly diffuse initial finite measure Z 0   , Z t   collapses into a discrete measure with a finite support as soon as t > 0   . Then we can identify Z t   interchangeably with a finite collection of spatially distributed particles with masses. When there is such a particle at a fixed location, we obtain the Laplace transform of its mass. The total number of particles in Z t   decreases in t   due to both branching and coalescing. When there is only one particle left at time t   , we also recover the joint distribution of its location and its mass. Eventually, all the particles will die out. We further find the distribution of the location where the last particle disappears. Coincidentally, super Brownian motion shares the same near extinction behavior.
Connections between superprocesses and stochastic flows have been noticed before. In [9] a superprocess was obtained from the empirical measure of a coalescing flow. Arratia flow serves as a fundamental example of coalescing flow. In this paper we point out several connections between Z   and the Arratia flow. More precisely, the support of Z t   at a fixed time t > 0   can be identified with a Cox process whose intensity measure is determined by the Arratia flow. A version of Z t   can be constructed using the Arratia flow. The general Laplace functional for Z   can also be expressed in terms of the Arratia flow.
Replacing the Feller's branching process by the square of Bessel process, we discuss a more general model at the end of this paper. The mass-dimension evolution of such a model can also be characterized by coalescing Brownian motions.
The rest of this paper is arranged as follows. As a preliminary, we first state and prove a dual relationship on coalescing Brownian motions in Section  2 . In Section  3 , we define the process Z   as a weak limit of the empirical measure for the branching-coalescing particle system. Then we proceed to prove the duality between Z   and coalescing Brownian motions. The uniqueness of Z   follows from such a duality immediately. We continue to study several properties of this process in Section  4 . We further discuss the connections between the Arratia flow and Z   in Section  5 .
At the end, we propose a more general model and establish its duality in Section 6.

2 Coalescing Brownian motions and their duality

An m   -dimensional coalescing Brownian motion can be described as follows. Consider a system of m   indexed particles with locations in R   that evolve as follows. Each particle moves according to an independent standard Brownian motion on R   until two particles are at the same location.
At this moment a coalescence occurs and the particle of higher index starts to move together with the particle of lower index. We say the particle with higher index is attached to the particle with lower index, which is still free. The particle system then continues its evolution in the same fashion. Note that indices are not essential here, the collection of locations of the particles is Markovian in its own right, but it will be convenient to think of the process as taking values in R m   rather than subsets of R   with at most m   elements. For definiteness, throughout this section we will further assume that the particles are indexed in increasing order of their initial positions: it it clear that the dynamics preserve this ordering. Call the resulting Markov process X = ( X 1 , … , X m )   .
Write 1 { B } ( . )   for the indicator function of a set B   . The distribution of X ( t )   is uniquely specified by knowing for each choice of y 1 < y 2 < … < y n   the joint probabilities of which “balls” X 1 ( t ) , X 2 ( t ) , … , X m ( t )   lie in which of the “boxes” ] y 1 , y 2 ] , ] y 2 , y 3 ] , … , ] y n − 1 , y n ]   . That is, the distribution of X ( t )   is determined by the joint distribution of the indicators I i j → ( t , y ) : = 1 { X i ( t ) ∈ ] y j , y j + 1 ] }   for 1 ≤ i ≤ m   , 1 ≤ j ≤ n − 1   and y = ( y 1 , … , y n )   .
Suppose now that Y : = ( Y 1 , … , Y n )   is another coalescing Brownian motion. The distribution of Y ( t )   is uniquely specified by knowing for each choice of x 1 < x 2 < … < x n   the distribution of the indicators I i j ← ( t , x ) : = 1 { x i ∈ ] Y j ( t ) , Y j + 1 ( t ) ] }   for 1 ≤ i ≤ m   , 1 ≤ j ≤ n − 1   and x = ( x 1 , … , x m )   .
The next “balls-in-boxes” duality is crucial in characterizing the distribution of the measure-valued process concerned in this paper.
Theorem 2.1. Suppose in the notation above that X = ( X 1 , … , X m )   is an m   -dimensional coalescing Brownian motion and Y = ( Y 1 , … , Y n )   is an n   -dimensional coalescing Brownian motion. Then for each t ≥ 0   the joint distribution of the m × ( n − 1 )   -dimensional random array ( I i j → ( t , Y ( 0 ) ) )   coincides with that of the m × ( n − 1 )   -dimensional random array ( I i j ← ( t , X ( 0 ) ) )   .
Theorem  2.1 generalizes Theorem 1.1 in [16] . Some other more elaborate dualities on coalescing-reflecting Brownian systems can be found in [14] and [13] .
We first prove the counterpart of Theorem  2.1 for continuous time simple coalescing random walks, which is interesting in its own right. Notice that X   is a coalescing Brownian motion if and only if X i   is a ( ℱ t X )   -Brownian motion for each 1 ≤ i ≤ m   , and ( X j − X i ) / 2   is a ( ℱ t X )   -Brownian motion stopped at 0   , where ( ℱ t X )   denotes the filtration generated by X   . Then Theorem  2.1 follows from a straight forward martingale argument proof of the convergence of scaled random walk to Brownian motion.
A p   -simple random walk on Z   is a continuous time simple random walk that makes jumps at unit rate, and when it makes a jump from some site it jumps to the right neighbor with probability p   and to the left neighbor with probability 1 − p   . An m   -dimensional p   -simple coalescing random walk is defined in the same way as the coalescing Brownian motion at the beginning of this section. When p = 1 2   we just call this particle system a simple coalescing random walk.
Some notation is useful to keep track of the interactions among the particles in the coalescing system. Let P m   denote the set of interval partitions of the totality of indices N m : = { 1 , … , m }   .
That is, an element π   of P m   is a collection π = { A 1 ( π ) , … , A h ( π ) }   of disjoint subsets of N m   such that ∪ i A i ( π ) = N m   and a < b   for all a ∈ A i   , b ∈ A j   , i < j   . The sets A 1 ( π ) , … A h ( π )   consisting of consecutive indices are the intervals of the partition π   . The integer h   is the length of π   and is denoted by l ( π )   . Equivalently, we can think of P m   as a set of equivalence relations on N m   and write i ∼ π j   if i   and j   belong to the same interval of π ∈ P m   . Of course, if i ∼ π j   , then i ∼ π k ∼ π j   for all i ≤ k ≤ j   .
Given π ∈ P m   , define α i ( π ) : = min A i ( π )   to be the left-hand end-point of the i t h   interval A i ( π )   . Put Z π m : = { ( x 1 , … , x m ) ∈ Z m : x 1 ≤ … ≤ x m and x i = x j if i ∼ π j }   and Z ^ π m : = { ( x 1 , … , x m ) ∈ Z m : x 1 ≤ … ≤ x m and x i = x j if and only if i ∼ π j } .   Note that Z m   is the disjoint union of the sets Z ^ π m   , π ∈ P m   .
Write X = ( X 1 , … , X m )   for the coalescing random walk. If X ( t ) ∈ Z ^ π m   , then the free particles at time t   have indices α 1 ( π ) , … , α l ( π ) ( π )   and the i t h   particle at time t   is attached to the free particle with index min { j : 1 ≤ j ≤ m , j ∼ π i } = max { α k ( π ) : α k ( π ) ≤ i } .   In order to write down the generator of X   , we require a final piece of notation. Let { e i k : 1 ≤ i ≤ k }   be the set of coordinate vectors in Z k   ; that is, e i k   is the vector that has the i t h   coordinate 1   and all the other coordinates 0   . For π ∈ P m   , define a map K π : Z π m → Z l ( π )   by K π ( x ) = K π ( x 1 , … , x m ) : = ( x α 1 ( π ) , … , x α l ( π ) ( π ) )   Notice that K π   is a bijection between Z π m   and { x ∈ Z l ( π ) : x 1 ≤ x 2 ≤ … ≤ x l ( π ) }   , and we write K π − 1   for the inverse of K π   . For brevity, we will sometimes write x π   for K π ( x )   .
Write B ( Z m )   for the collection of all bounded functions on Z m   . The generator G   of X   is the operator G : B ( Z m ) → B ( Z m )   given by
G f ( x ) : = p ∑ i = 1 l ( π ) f ∘ K π − 1 ( x π + e i l ( π ) ) + ( 1 − p ) ∑ i = 1 l ( π ) f ∘ K π − 1 ( x π − e i l ( π ) ) − l ( π ) f ∘ K π − 1 ( x π ) , f ∈ B ( Z m ) , x ∈ Z ^ π m , π ∈ P m .
This expression is well-defined, because if x ∈ Z ^ π m   , then x π   , x π + e i l ( π )   and x π − e i l ( π )   are all in { x ∈ Z l ( π ) : x 1 ≤ x 2 ≤ … ≤ x l ( π ) }   .
Note: From now on we will suppress the dependence on dimension and write e i l ( π )   as e i   .
Write Z ′ : = Z + 1 2 = { i + 1 2 : i ∈ Z }   . An n   -dimensional q   -simple coalescing random walk on Z ′ n   and its generator H   can be defined in the obvious way. Such a process, with q = 1 − p   , will serve as the process dual to the p   -simple coalescing random walk on Z m   in the following way.
Fix x ∈ Z m   with x 1 ≤ … ≤ x m   and y ∈ Z ′ n   with y 1 ≤ … ≤ y n   . Put I i j → ( t , y ) : = 1 { X i ( t ) ∈ ] y j , y j + 1 ] }   and I i j ← ( t , x ) : = 1 { x i ∈ ] Y j ( t ) , Y j + 1 ( t ) ] }   for 1 ≤ i ≤ m   and 1 ≤ j ≤ n − 1   .
Lemma 2.2. Suppose in the notation above that X = ( X 1 , … , X m )   is an m   -dimensional Z m   -valued p   -simple coalescing random walk and Y = ( Y 1 , … , Y n )   is an n   -dimensional Z ′ n   -valued ( 1 − p )   -simple coalescing random walk. Then for each t ≥ 0   the joint distribution of the m × ( n − 1 )   -dimensional random array ( I i j → ( t , Y ( 0 ) ) )   coincides with that of the m × ( n − 1 )   -dimensional random array ( I i j ← ( t , X ( 0 ) ) )   .
  • Proof. For a function g : { 0 , 1 } m ( n − 1 ) → R   , a vector x ∈ Z m   with x 1 ≤ … ≤ x m   , and a vector y ∈ Z ′ n   with y 1 ≤ … ≤ y n   , set g ¯ ( x ; y ) : = g ( 1 { ] y 1 , y 2 ] } ( x 1 ) , … , 1 { ] y n − 1 , y n ] } ( x 1 ) , … , 1 { ] y 1 , y 2 ] } ( x m ) , … , 1 { ] y n − 1 , y n ] } ( x m ) ) .   We may assume that X   and Y   are defined on the same probability space ( Ω , ℱ , P )   . We need to show that
    P [ g ¯ ( X ( t ) ; Y ( 0 ) ) ] = P [ g ¯ ( X ( 0 ) ; Y ( t ) ) ] . (2.1)
    For x ∈ Z m   , put g ¯ x ( ⋅ ) : = g ¯ ( x ; ⋅ )   , and for y ∈ Z ′ n   , put g ¯ y ( ⋅ ) : = g ¯ ( ⋅ ; y )   . In order to establish ( 2.1 ), it suffices by a standard argument (cf. Section 4.4 in [6] ) to show that
    G ( g ¯ y ) ( x ) = H ( g ¯ x ) ( y ) (2.2)
    (recall that G   and H   are the generators of X   and Y   , respectively). Fix x ∈ Z ^ π m   and y ∈ Z ^ ϖ ′ n   for some π ∈ P m   and ϖ ∈ P n   . Put I + : = { i : 1 ≤ i ≤ l ( π ) , x α i ( π ) + 1 2 = y α j ( ϖ ) for some 1 ≤ j ≤ l ( ϖ ) }   and I − : = { i : 1 ≤ i ≤ l ( π ) , x α i ( π ) − 1 2 = y α j ( ϖ ) for some 1 ≤ j ≤ l ( ϖ ) } .   Similarly, put J − : = { j : 1 ≤ j ≤ l ( ϖ ) , y α j ( ϖ ) − 1 2 = x α i ( π ) for some 1 ≤ i ≤ l ( π ) }   and J + : = { j : 1 ≤ j ≤ l ( ϖ ) , y α j ( ϖ ) + 1 2 = x α i ( π ) for some 1 ≤ i ≤ l ( π ) } .   Recall that x α 1 ( π ) < … < x α l ( π ) ( π )   and y α 1 ( ϖ ) < … < y α l ( ϖ ) ( ϖ )   . Therefore, for each i ∈ I +   there is a unique j ∈ J −   such that x α i ( π ) + 1 2 = y α j ( ϖ )   and vice versa. Fix such a pair ( i , j )   .
    Observe that x ′ : = x + ∑ k ∈ A i ( π ) e k m = K π − 1 ( x π + e i )   and y ′ : = y − ∑ k ∈ A j ( ϖ ) e k n = K π − 1 ( y ϖ − e j ) .   We are going to verify that
    ( 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) ) = ( 1 { ] y j ′ ′ , y j ′ + 1 ′ ] } ( x i ′ ) ) (2.3)
    by considering all the possible scenarios.
    Given any i ′ ∈ A i ( π )   we have:
    • ∙   for j ′ = α j ( ϖ ) − 1   ,
      1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ + 1 ) = 0 = 1 { ] y j ′ , y j ′ + 1 − 1 ] } ( x i ′ ) = 1 { ] y j ′ ′ , y j ′ + 1 ′ ] } ( x i ′ ) ,
    • ∙   for α j ( ϖ ) ≤ j ′ < max A j ( ϖ )   ,
      1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ + 1 ) = 0 = 1 { ] y j ′ − 1 , y j ′ + 1 − 1 ] } ( x i ′ ) = 1 { ] y j ′ ′ , y j ′ + 1 ′ ] } ( x i ′ ) ,
    • ∙   for j ′ = max A j ( ϖ )   ,
      1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ + 1 ) = 1 = 1 { ] y j ′ − 1 , y j ′ + 1 ] } ( x i ′ ) = 1 { ] y j ′ ′ , y j ′ + 1 ′ ] } ( x i ′ ) ,
    • ∙   and for j ′ < α j ( ϖ ) − 1   or j ′ > max A j ( ϖ )   ,
      1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ + 1 ) = 0 = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ) = 1 { ] y j ′ ′ , y j ′ + 1 ′ ] } ( x i ′ ) .
    Moreover, given any i ′ ∉ A i ( π )   , we have x i ′ ≠ x α i ( π )   . Hence
    • ∙   for j ′ = α j ( ϖ ) − 1   , 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ) = 1 { ] y j ′ , y j ′ + 1 − 1 ] } ( x i ′ ) = 1 { ] y j ′ ′ , y j ′ + 1 ′ ] } ( x i ′ ) ,  
    • ∙   for j ′ = max A j ( ϖ )   , 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ) = 1 { ] y j ′ − 1 , y j ′ + 1 ] } ( x i ′ ) = 1 { ] y j ′ ′ , y j ′ + 1 ′ ] } ( x i ′ ) ,  
    • ∙   for α j ( ϖ ) ≤ j ′ < max A j ( ϖ )   , 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ) = 1 { ] y j ′ − 1 , y j ′ + 1 − 1 ] } ( x i ′ ) = 1 { ] y j ′ ′ , y j ′ + 1 ′ ] } ( x i ′ ) ,  
    • ∙   and for j ′ < α j ( ϖ ) − 1   or j ′ > max A j ( ϖ )   , 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ′ ) = 1 { ] y j ′ , y j ′ + 1 ] } ( x i ′ ) = 1 { ] y j ′ ′ , y j ′ + 1 ′ ] } ( x i ′ ) .  
    Combining the above observations yields ( 2.3 ).
    Therefore, g ¯ y ∘ K π − 1 ( x π + e i ) = g ¯ x ∘ K π − 1 ( y ϖ − e j ) .   Furthermore, it is easy to see for i ′ ∉ I +   that g ¯ y ∘ K π − 1 ( x π + e i ′ ) = g ¯ y ∘ K π − 1 ( x π )   and for j ′ ∉ J −   that g ¯ x ∘ K ϖ − 1 ( y ϖ − e j ′ ) = g ¯ x ∘ K ϖ − 1 ( y ϖ ) .   Similarly, for any i ∈ I −   there exists a unique j ∈ J +   such that x α i ( π ) − 1 2 = y α j ( ϖ )   and vice versa. For such a pair ( i , j )   we have g ¯ y ∘ K π − 1 ( x π − e i ) = g ¯ x ∘ K π − 1 ( y ϖ + e j ) .   Furthermore, we see for i ′ ∉ I −   that g ¯ y ∘ K π − 1 ( x π − e i ′ ) = g ¯ y ∘ K π − 1 ( x π )   and for j ′ ∉ J +   that g ¯ x ∘ K ϖ − 1 ( y ϖ + e j ′ ) = g ¯ x ∘ K ϖ − 1 ( y ϖ ) .   Lastly, note that g ¯ y ∘ K π − 1 ( x π ) = g ¯ ( x ; y ) = g ¯ x ∘ K ϖ − 1 ( y ϖ )   and so
    G ( g ¯ y ) ( x ) − H ( g ¯ x ) ( y ) = p ∑ i = 1 l ( π ) ( g ¯ y ∘ K π − 1 ( x π + e i ) − g ¯ y ∘ K π − 1 ( x π ) ) + ( 1 − p ) ∑ i = 1 l ( π ) ( g ¯ y ∘ K π − 1 ( x π − e i ) − g ¯ y ∘ K π − 1 ( x π ) ) − p ∑ j = 1 l ( ϖ ) ( g ¯ x ∘ K ϖ − 1 ( y ϖ − e i ) − g ¯ x ∘ K ϖ − 1 ( y ϖ ) ) − ( 1 − p ) ∑ j = 1 l ( ϖ ) ( g ¯ x ∘ K ϖ − 1 ( y ϖ + e i ) − g ¯ x ∘ K ϖ − 1 ( y ϖ ) ) = p ∑ i ∈ I + g ¯ y ∘ K π − 1 ( x π + e i ) − p ∑ j ∈ J − g ¯ x ∘ K ϖ − 1 ( y ϖ − e j ) + ( 1 − p ) ∑ i ∈ I − g ¯ y ∘ K π − 1 ( x π − e i ) − ( 1 − p ) ∑ j ∈ J + g ¯ x ∘ K ϖ − 1 ( y ϖ + e j ) = 0 ,
    as required. □
Remark 2.3. For discrete time coalescing simple random walks such a duality is evident from Fig. 7 in [13] . But the duality seems to be less apparent for continuous time coalescing simple random walk.

3 Existence and uniqueness

A construction of Z   was given in [5] using Feller's branching excursions. In this paper we adopt a weak convergence approach, which is commonly used in the study of measure-valued processes.
Recall that a nonnegative valued process ξ   is a Feller's branching process with initial value x ≥ 0   if it is the unique strong solution to the following stochastic differential equation ξ ( t ) = x + ∫ 0 t γ ξ ( s ) d B ( s ) ,   where γ   is a positive constant and B   is a one-dimensional Brownian motion. ξ ( t )   is a martingale.
It has a Laplace transform
P [ exp { − λ ξ ( t ) } ] = exp { − 2 λ x 2 + λ γ t } ; (3.1)
its extinction probability is given by P { ξ ( t ) = 0 } = exp { − 2 x γ t } ;   see Section II.1 and II.5 in [10] .
Observe that independent Feller's branching processes are additive; i.e. if ξ   and η   are two independent Feller's branching processes (with the same parameter γ   ), then ξ + η   is also a Feller's branching process. This fact will be used repeatedly in our discussions.
Write M F ( R )   for the space of finite measures on R   equipped with the topology of weak convergence. Given any finite measure Z 0   on R   , put z ¯ : = Z 0 ( R )   . For any positive integer m   , let ( ξ 1 ( m ) , … , ξ m ( m ) )   be a collection of m   independent Feller's branching processes each with initial value z ¯ / m   . Choose ( x 1 , … , x n )   to be i.i.d. samples from distribution Z 0 / z ¯   . Let ( X 1 ( m ) , … , X m ( m ) )   be an m   -dimensional coalescing Brownian motion starting at ( x 1 , … , x m )   .
Moreover, we always assume that ( ξ 1 ( m ) , … , ξ m ( m ) )   and ( X 1 ( m ) , … , X m ( m ) )   are independent. Let δ x   denote the point mass at x ∈ R   . Then Z t ( m ) : = ∑ i = 1 m ξ i ( m ) ( t ) δ X i ( m ) ( t )   defines a M F ( R )   -valued process. From now on we will suppress the dependence of m   in ξ i ( m )   and X i ( m )   .
Recall that a collection of processes { Z α , α ∈ I }   with sample paths in D ( M F ( R ) )   is C-relatively compact if it is relatively compact and all its weak limits are a.s. continuous. The proof of the next lemma is standard; see, e.g. the proofs for Lemma 3.2 in [16] and Proposition II.4.2 in [10] .
Lemma 3.1. { Z ( m ) }   is C-relatively compact.
  • Proof. We first check the compact containment condition. For any ε > 0   and T > 0   , choose a compact set K 0 ⊂ D ( R )   such that P { X 1 ∈ K 0 c } < ε 2   . Let K : = { x t : x ∈ K 0 , t ≤ T }   . Then K   is compact in R   , and P { X 1 ( t ) ∈ K c , ∃ t ≤ T } ≤ P { X 1 ∈ K 0 c } < ε 2 .   Write N m : = # I K : = # { 1 ≤ i ≤ m : X i ( t ) ∈ K c , ∃ t ≤ T } ,   where # I K   denotes the cardinality of the indices set I K   . Conditioning on N m   , by the additivity for Feller's branching processes, we see that ∑ i ∈ I K ξ i   is a Feller's branching Process with initial value N m z ¯ / m   . Then by Doob's maximal inequality, P { sup 0 ≤ t ≤ T ∑ i ∈ I K ξ i ( t ) > ε | N m } ≤ N m z ¯ m ε .   Therefore,
    P { sup 0 ≤ t ≤ T Z t ( m ) ( K c ) > ε } ≤ P { sup 0 ≤ t ≤ T ∑ i ∈ I K ξ i ( t ) > ε } ≤ P [ N m ] z ¯ m ε ≤ z ¯ ε .
    For any f ∈ C b 2 ( R )   put Z t ( m ) ( f ) : = ∫ − ∞ ∞ f ( x ) Z t ( m ) ( d x )   . Now we are going to show that { Z . ( m ) ( f ) }   is C-relatively compact in D ( R )   . By Itô's formula, we have
    Z t ( m ) ( f ) = ∑ i = 1 m [ z ¯ m f ( x i ) + ∫ 0 t f ( X i ( s ) ) d ξ i ( s ) + ∫ 0 t ξ i ( s ) f ′ ( X i ( s ) ) d X i ( s ) + 1 2 ∫ 0 t ξ i ( s ) f ′ ′ ( X i ( s ) ) d s ] .
    The additivity for ( ξ i )   gives P [ sup 0 ≤ s ≤ t ∑ i = 1 m ξ i ( s ) ] < ∞ and P [ sup 0 ≤ s ≤ t ∑ i , j = 1 m ξ i ( s ) ξ j ( s ) ] < ∞ , t > 0 ,   then ∑ i = 1 m ∫ 0 t ξ i ( s ) f ′ ′ ( X i ( s ) ) d s   is C-relatively compact following from the Arzela-Ascoli theorem and Proposition VI.3.26 of [7] .
    Note that 〈 ∑ i = 1 m ∫ 0 . ξ i ( s ) f ′ ( X i ( s ) ) d X i ( s ) 〉 t = ∑ i , j = 1 m ∫ 0 t ξ i ( s ) ξ j ( s ) f ′ ( X i ( s ) ) f ′ ( X j ( s ) ) d 〈 X i , X j 〉 s ,   where 〈 X i , X j 〉 s = s − T i j ∧ s   and T i j : = inf { s ≥ 0 : X i ( s ) = X j ( s ) }   . By Arzela-Ascoli theorem again, { 〈 ∑ i = 1 m ∫ 0 . ξ i ( s ) f ′ ( X i ( s ) ) d X i ( s ) 〉 . }   is C-relatively compact. Theorem VI.4.13 and Proposition VI.3.26 in [7] then imply that the collection of martingales { ∑ i = 1 m ∫ 0 . ξ i ( s ) f ′ ( X i ( s ) ) d X i ( s ) }   is C-relatively compact.
    Similarly, { ∑ i = 1 m ∫ 0 . f ( X i ( s ) ) d ξ i ( s ) }   is also C-relatively compact. Moreover, 1 m ∑ i = 1 m f ( x i ) → Z 0 ( f ) a.s. .   { Z ( m ) ( f ) }   is thus C-relatively compact. Consequently, by Theorem II.4.1 in [10] we can conclude that { Z ( m ) }   is C-relatively compact.
    □
Write Z   for the weak limit of { Z ( m ) }   . The Laplace functional of Z   can be obtained from the duality in Theorem  2.1 . As a result, its uniqueness is settled.
In the sequel we always write ( Y 1 , … , Y 2 n )   for a coalescing Brownian motion starting at ( y 1 , … , y 2 n )   with y 1 ≤ … ≤ y 2 n   .
Theorem 3.2. { Z ( m ) }   has a unique weak limit Z   in C ( M F ( R ) )   . Given a j > 0 , j = 1 , … , n   , for any y 1 ≤ y 2 ≤ … ≤ y 2 n   and any t > 0   , we have
P [ exp { − ∑ j = 1 n a j Z t ( ] y 2 j − 1 , y 2 j ] ) } ] = P [ exp { − ∫ − ∞ ∞ Z 0 ( d x ) 2 ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) 2 + γ t ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) } ] . (3.2)
  • Proof. We first condition on ( ξ 1 ( t ) , … , ξ m ( t ) )   . By Theorem  2.1 we can show that
    P [ exp { − ∑ j = 1 n ∑ i = 1 m a j ξ i ( t ) 1 { ] y 2 j − 1 , y 2 j ] } ( X i ( t ) ) } | ( ξ i ( t ) ) ] = P [ exp { − ∑ j = 1 n ∑ i = 1 m a j ξ i ( t ) 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x i ) } | ( ξ i ( t ) ) ] , (3.3)
    where ( Y 1 , … , Y 2 n )   is independent of ( X i )   and ( ξ i )   .
    Now take expectations on both sides of ( 3.3 ) and then condition on ( x i )   and ( Y i ( t ) )   . Since ξ 1 ( t ) , … , ξ m ( t )   are independent of each other, and they are independent of ( x i )   and ( Y i ( t ) )   , it follows from ( 3.1 ) that
    P [ exp { − ∑ j = 1 n a j Z t ( m ) ( ] y 2 j − 1 , y 2 j ] ) } ] = P [ P [ exp { − ∑ j = 1 n ∑ i = 1 m a j ξ i ( t ) 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x i ) } | ( x i ) , ( Y i ( t ) ) ] ] = P [ m ∏ i = 1 exp { − 2 z ¯ ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x i ) m { 2 + γ t ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x i ) } } ] = P [ ( 1 z ¯ ∫ − ∞ ∞ Z 0 ( d x ) exp { − 2 z ¯ ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) m ( 2 + γ t ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) ) } ) m ] . (3.4)
    Let Z   be any weak limit of { Z ( m ) }   . Let m → ∞   in ( 3.4 ). Then
    lim m → ∞ P [ exp { − ∑ j = 1 n a j Z t ( m ) ( ] y 2 j − 1 , y 2 j ] ) } ] = lim m → ∞ P [ ( 1 − ∫ − ∞ ∞ Z 0 ( d x ) 2 ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) m ( 2 + γ t ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) ) ) m ] = P [ exp { − ∫ − ∞ ∞ Z 0 ( d x ) 2 ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) 2 + γ t ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) } ] . (3.5)
    We still need to make sure that
    P [ exp { − ∑ j = 1 n a j Z t ( ] y 2 j − 1 , y 2 j ] ) } ] = lim m → ∞ P [ exp { − ∑ j = 1 n a j Z t ( m ) ( ] y 2 j − 1 , y 2 j ] ) } ] . (3.6)
    To this end we can suppose that y 1 < y 2 < … < y 2 n   . Then for small enough ε > 0   , similar to ( 3.4 ) we have
    P [ exp { − ∑ j = 1 n a j Z t ( m ) ( ] y 2 j − 1 + ε , y 2 j − ε ] ) } ] − P [ exp { − ∑ j = 1 n a j Z t ( m ) ( ] y 2 j − 1 − ε , y 2 j + ε ] ) } ] ≤ 1 − P [ exp { − ∑ j = 1 n a j Z t ( m ) ( ] y 2 j − 1 − ε , y 2 j − 1 + ε ] ∪ ] y 2 j − ε , y 2 j + ε ] ) } ] = 1 − P [ ( 1 z ¯ ∫ − ∞ ∞ Z 0 ( d x ) exp { − 2 z ¯ ∑ j = 1 n a j 1 { ] Y 2 j − 1 ′ ( t ) , Y 2 j − 1 ′ ′ ( t ) ] ∪ ] Y 2 j ′ ( t ) , Y 2 j ′ ′ ( t ) ] } ( x ) m ( 2 + γ t ∑ j = 1 n a j 1 { ] Y 2 j − 1 ′ ( t ) , Y 2 j − 1 ′ ′ ( t ) ] ∪ ] Y 2 j ′ ( t ) , Y 2 j ′ ′ ( t ) ] } ( x ) ) } ) m ] ≤ 1 − P { ∩ j = 1 2 n { Y j ′ ( t ) = Y j ′ ′ ( t ) } } , (3.7)
    where ( Y 1 ′ , Y 1 ′ ′ , … , Y 2 n ′ , Y 2 n ′ ′ )   is a coalescing Brownian motion starting at ( y 1 − ε , y 1 + ε , … , y 2 n − ε , y 2 n + ε )   . Clearly, ( 3.7 ) converges (uniformly in m   ) to 0   as ε → 0 +   . So, ( 3.6 ) holds.
    It is clear that the distribution of Z   is uniquely determined by ( 3.2 ). So, Z   is the unique weak limit of { Z ( m ) }   .
    □
The moments of Z   can be obtained immediately from ( 3.2 ).
Proposition 3.3. For any y 1 ≤ y 2 ≤ … ≤ y 2 n   and t > 0   , we have P [ Z t ( ∑ j = 1 n a j 1 { ] y 2 j − 1 , y 2 j ] } ) ] = P [ Z 0 ( ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ) ]   and
P [ Z t 2 ( ∑ j = 1 n a j 1 { ] y 2 j − 1 , y 2 j ] } ) ] = P [ Z 0 2 ( ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ) ] + P [ γ t Z 0 ( ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ) ] .
Remark 3.4. Another consequence of duality ( 3.2 ) is that Z   is a Markov process; see Theorem 3.3 in [16] for a proof on a similar model.
Martingale problem is often used to characterize a superprocess. Z   is the solution to the martingale problem (see [5] ): for any φ ∈ C 2 ( R )   , M t ( φ ) = Z t ( φ ) − Z 0 ( φ ) − 1 2 ∫ 0 t Z s ( φ ′ ′ ) d s , t ≥ 0 ,   is a continuous martingale relative to ( ℱ t ) t ≥ 0   with quadratic variation process 〈 M ( φ ) 〉 t = γ ∫ 0 t Z s ( φ 2 ) d s + ∫ 0 t d s ∫ Δ φ ′ ( x ) φ ′ ( y ) Z s ( d x ) Z s ( d y ) ,   where Δ = { ( x , x ) : x ∈ R }   .
But a remarkable feature of such a martingale problem is that its solution is not unique.
For example, let ξ 1   and ξ 2   be two independent branching processes each with initial value 1   .
Let B 1   and B 2   be two independent Brownian motions. Assume that ( ξ 1 , ξ 2 )   and ( B 1 , B 2 )   are independent. Then Z t ′ : = ξ 1 ( t ) δ B 1 ( t ) + ξ 2 ( t ) δ B 1 ( t )   is another solution to this martingale problem; also see [16] for a similar counter example.

4 Some properties

Our first result in this section is a straight forward consequence of Theorem  3.2 .
Proposition 4.1. For any y 1 ≤ y 2 ≤ … ≤ y 2 n   and t > 0   , we have
P [ exp { − λ ∑ j = 1 n Z t ( ] y 2 j − 1 , y 2 j ] ) } ] = P [ exp { − 2 λ 2 + λ γ t ∑ j = 1 n Z 0 ( ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] ) } ] , λ > 0 . (4.1)
  • Proof. Observe that ∑ j = 1 n 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x )   is either 0   or 1   , then ( 4.1 ) follows readily from ( 3.2 ).
    □
Proposition  4.1 allows us to carry out some explicit computations on Z   . First, by letting λ → ∞   we can easily see that P [ Z t ( ( − ∞ , ∞ ) ) = 0 ] = exp { − 2 z ¯ γ t } .   We are going to further study the probability that Z t   does not charge on an arbitrary finite interval. For any x , y , a   and b   , put x ~ : = x − y 2 , y ~ : = x + y 2 , a ~ : = a + b 2 and b ~ : = b − a 2 .  
Proposition 4.2. Given a < b   and t > 0   , we have
P { Z t ( ] a , b ] ) = 0 } = 1 2 π t ∫ − ∞ ∞ d x ∫ 0 ∞ d y exp { − 2 Z 0 ( ] x ~ , y ~ ] ) t γ − ( x − a ~ ) 2 2 t } ( exp { − ( y − b ~ ) 2 2 t } − exp { − ( y + b ~ ) 2 2 t } ) + 2 2 π t ∫ b ~ ∞ d x exp { − x 2 2 t } . (4.2)
  • Proof. Let λ → ∞   in ( 4.1 ) we have
    P { Z t ( ] a , b ] ) = 0 } = P [ exp { − 2 γ t Z 0 ( ] Y 1 ( t ) , Y 2 ( t ) ] ) } ] = P [ exp { − 2 γ t Z 0 ( ] Y 1 ( t ) , Y 2 ( t ) ] ) } ; Y 1 ( t ) ≠ Y 2 ( t ) ] + P { Y 1 ( t ) = Y 2 ( t ) } ,
    where ( Y 1 , Y 2 )   is a coalescing Brownian motion starting from ( a , b )   .
    To find the distribution of ( Y 1 , Y 2 )   , one could rotate the coordinate system anti-clockwise by π / 4   . Under the new coordinate system ( Y 1 , Y 2 )   becomes a process ( Y 1 ′ , Y 2 ′ )   such that Y 1 ′   is aBrownian motion starting at a ~   , Y 2 ′   is a Brownian motion starting at b ~   and stopped at 0   , and Y 1 ′   and Y 2 ′   are independent. So ( 4.2 ) just follows from the reflection principle for Brownian motion.
    □
Write S t   for the support of Z t   . Intuitively, starting with particles of a total initial mass Z 0 ( R )   , as soon as t > 0   the particles near − ∞   and ∞   will die out due to branching. Z t   is then expected to be supported by a finite set because of coalescence. The next two results concern the cardinality of S t   .
Proposition 4.3. Given a < b   and t > 0   , we have
P [ # S t ∩ ] a , b ] ] = b − a π t − 1 2 π t 2 ∫ a b d z ∫ − ∞ ∞ d x ∫ 0 ∞ d y y exp { − 2 Z 0 ( ] x ~ , y ~ ] ) t γ − ( x − 2 z ) 2 + y 2 2 t } . (4.3)
  • Proof. It is easy to see from ( 4.2 ) that for any z ∈ R   ,
    P { Z t ( d z ) ≠ 0 } = d z π t − d z 2 π t 2 ∫ − ∞ ∞ d x ∫ 0 ∞ d y y exp { − 2 Z 0 ( ] x ~ , y ~ ] ) t γ − ( x − 2 z ) 2 + y 2 2 t } . (4.4)
    Then ( 4.3 ) is obtained by taking integrations on both sides of ( 4.4 ) from a   to b   .
    □
Proposition 4.4. With probability 1   , # S t < ∞ , ∀ t > 0   .
  • Proof. Given s > 0   , we first claim that P [ # S s ] < ∞   if Z 0   has a bounded support. Suppose that Z 0 ( ] − b , b ] ) = 1   for some b > 0   . Then by ( 4.4 ),
    P [ # S s ] = 1 2 π s 2 ∫ − ∞ ∞ d z ∫ − ∞ ∞ d x ∫ 0 ∞ d y y ( 1 − exp { − 2 Z 0 ( ] x ~ , y ~ ] ) s γ } ) exp { − ( x − 2 z ) 2 + y 2 2 s } = 1 2 π s s ∫ − ∞ ∞ d x ∫ 0 ∞ d y y ( 1 − exp { − 2 Z 0 ( ] x ~ , y ~ ] ) s γ } ) exp { − y 2 2 s } ≤ 1 2 π s s ∫ 0 ∞ d y ∫ − y − 2 b y + 2 b d x y exp { − y 2 2 s } < ∞ .
    Our claim is proved.
    Now given any integer j   , let η j ( s )   be the Feller's branching process with initial value η j ( 0 ) : = Z 0 ( ] j , j + 1 ] )   . Since ∑ j = − ∞ ∞ P { η j ( s ) ≠ 0 } = ∑ j = − ∞ ∞ ( 1 − exp { − 2 η j ( 0 ) γ s } ) ≤ ∑ j = − ∞ ∞ 2 η j ( 0 ) γ s = 2 z ¯ γ s ,   by Borel-Cantelli lemma we have that, with probability 1   , η j ( s ) ≠ 0   for only finitely many values of j   .
    Therefore, for any t > 0   , with probability 1   , Z t / 2   must have a bounded support. The Markov property for Z   , together with the claim from the first part of the proof, implies that # S t < ∞   a.s..
    Finally, by the Markov property for Z   we conclude that P { # S t < ∞ , ∀ t > 0 } = 1   . □
By Proposition  4.4 , as soon as t > 0   , S t   becomes a finite set. For any z ∈ S t   , we associate it with a particle located at z   with mass Z t ( { z } )   . We can thus identify Z t   interchangeably with a collection of spatially distributed particles with masses. As time goes on, the total number of particles decreases either because two “alive” particles coalesce into one particle, or because each particle disappears due to its branching.
Since # S t < ∞   , a small neighborhood of z   contains at most one particle in Z t   . When there is such a particle, we want to find the distribution of its mass. Formally, we are looking for P [ exp { − λ Z t ( { z } ) } ; Z t ( { z } ) > 0 ] .  
Proposition 4.5. For any z ∈ R   and t > 0   , we have
P [ exp { − λ Z t ( d z ) } ; Z t ( d z ) > 0 ] = d z 2 π t 2 ∫ − ∞ ∞ d x ∫ 0 ∞ d y y ( exp { − 2 λ Z 0 ( ] x ~ , y ~ ] ) 2 + λ γ t } − exp { − 2 Z 0 ( ] x ~ , y ~ ] ) t γ } ) exp { − ( x − 2 z ) 2 + y 2 2 t } . (4.5)
  • Proof. We fix ( ξ i ( t ) )   first. Apply Theorem  2.1 to exp { − λ ∑ i = 1 m ξ i ( t ) 1 { ] a , b ] } ( X i ( t ) ) } 1 { ∑ i = 1 m ξ i ( t ) 1 { ] a , b ] } ( X i ( t ) ) > 0 } .   Then condition on ( Y 1 ( t ) , Y 2 ( t ) )   and take an expectation with respect to ( ξ i ( t ) )   . Similar to the proof for Theorem  3.2 we have that
    P [ exp { − λ ∑ i = 1 m ξ i ( t ) 1 { ] a , b ] } ( X i ( t ) ) } ; ∑ i = 1 m ξ i ( t ) 1 { ] a , b ] } ( X i ( t ) ) > 0 ] = P [ exp { − λ ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) } ; ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) > 0 ] = P [ exp { − λ ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) } ] − P [ ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) = 0 ] = P [ exp { − 2 λ Z 0 ( m ) ( ] Y 1 ( t ) , Y 2 ( t ) ] ) 2 + λ γ t } ] − P [ exp { − 2 Z 0 ( m ) ( ] Y 1 ( t ) , Y 2 ( t ) ] ) γ t } ] . (4.6)
    Therefore,
    P [ exp { − λ Z t ( ] a , b ] ) } ; Z t ( ] a , b ] ) > 0 ] = P [ exp { − 2 λ Z 0 ( ] Y 1 ( t ) , Y 2 ( t ) ] ) 2 + λ γ t } ] − P [ exp { − 2 Z 0 ( ] Y 1 ( t ) , Y 2 ( t ) ] ) γ t } ] = 1 2 π t ∫ − ∞ ∞ d x ∫ 0 ∞ d y ( exp { − 2 λ Z 0 ( ] x ~ , y ~ ] ) 2 + λ γ t } − exp { − 2 Z 0 ( ] x ~ , y ~ ] ) t γ } ) exp { − ( x − a ~ ) 2 2 t } ( exp { − ( y − b ~ ) 2 2 t } − exp { − ( y + b ~ ) 2 2 t } ) . (4.7)
    So, ( 4.5 ) is obtained by letting b → a +   .
    □
At a fixed time t > 0   , with a positive probability there can be only one particle (with a positive mass) left. When this happens, we are interested in the joint distribution of the mass and the location of that particle. More precisely, we want to find P [ exp { − λ Z t ( R ) } ; Z t ( R ) ≠ 0 , S t ⊂ d z ] .  
Proposition 4.6. For any z ∈ R   and t > 0   , we have
P [ exp { − λ Z t ( R ) } ; Z t ( R ) ≠ 0 , S t ⊂ d z ] = d z 2 π t 2 ∫ − ∞ ∞ d x ∫ 0 ∞ d y y exp { − 2 λ Z 0 ( ] x ~ , y ~ ] ) 2 + λ γ t − 2 Z 0 ( ] x ~ , y ~ ] c ) γ t − ( x − 2 z ) 2 + y 2 2 t } − d z π t exp { − 2 z ¯ γ t } . (4.8)
  • Proof. Put B : = { ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] c } ( x i ) = 0 } for x i : = X i ( 0 ) .   It follows from Theorem  2.1 that
    P [ exp { − λ ∑ i = 1 m ξ i ( t ) } ; ∑ i = 1 m ξ i ( t ) 1 { ] a , b ] } ( X i ( t ) ) > 0 , ∑ i = 1 m ξ i ( t ) 1 { ] a , b ] c } ( X i ( t ) ) = 0 ] = P [ exp { − λ ∑ i = 1 m ξ i ( t ) } ; ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) > 0 , B ] = P [ exp { − λ ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) } ; ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) > 0 , B ] = P [ exp { − λ ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) } ; B ] − P [ ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) = 0 ; B ] = P [ exp { − 2 λ z ¯ ∑ i = 1 m 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) m ( 2 + λ γ t ) } exp { − 2 z ¯ ∑ i = 1 m 1 { ] Y 1 ( t ) , Y 2 ( t ) ] c } ( x i ) m γ t } ] − P { ∑ i = 1 m ξ i ( t ) = 0 } , (4.9)
    where in obtaining the last equation we have used the fact that, given ( Y 1 ( t ) , Y 2 ( t ) )   and ( x i )   , exp { − λ ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] } ( x i ) } and event { ∑ i = 1 m ξ i ( t ) 1 { ] Y 1 ( t ) , Y 2 ( t ) ] c } ( x i ) = 0 }   are independent.
    Now letting m → ∞   in ( 4.9 ) we have
    P [ exp { − λ Z t ( R ) } ; Z t ( R ) ≠ 0 , S t ⊂ ( a , b ) ] = P [ exp { − 2 λ Z 0 ( ] Y 1 ( t ) , Y 2 ( t ) ] ) 2 + λ γ t − 2 Z 0 ( ] Y 1 ( t ) , Y 2 ( t ) ] c ) γ t } ] − exp { − 2 z ¯ γ t } = 1 2 π t ∫ − ∞ ∞ d x ∫ 0 ∞ d y exp { − 2 λ Z 0 ( ] x ~ , y ~ ] ) 2 + λ γ t − 2 Z 0 ( ] x ~ , y ~ ] c ) γ t } exp { − ( x − a ~ ) 2 2 t } ( exp { − ( y − b ~ ) 2 2 t } − exp { − ( y + b ~ ) 2 2 t } ) + 2 2 π t exp { − 2 γ t } ∫ b ~ ∞ d x exp { − x 2 2 t } − exp { − 2 z ¯ γ t } . (4.10)
    Finally, ( 4.8 ) is obtained by letting b → a +   in ( 4.10 ).
    □
Remark 4.7. Let λ = 0   in ( 4.10 ). We then obtain a result on the range of S t   .
The total number of particles in Z   will decrease one by one. Put τ : = inf { s ≥ 0 : # S s = 1 } .   Then τ < ∞   is the first time when there is exactly one particle left. The distribution of τ   is given in the following Proposition.
Proposition 4.8.
P { τ < t } = ∫ − ∞ ∞ d z ( 1 2 π t 2 ∫ − ∞ ∞ d x ∫ 0 ∞ d y y exp { − 2 Z 0 ( ] x ~ , y ~ ] c ) γ t − ( x − 2 z ) 2 + y 2 2 t } − 1 π t exp { − 2 z ¯ γ t } ) + 1 − exp { 2 z ¯ γ t } . (4.11)
  • Proof. Observe that
    P { τ < t } = ∫ − ∞ ∞ P { Z t ( R ) ≠ 0 , S t ⊂ d z } + P { Z t ( R ) = 0 } ,
    then ( 4.11 ) follows from Proposition  4.6 .
    □
Let T : = inf { t ≥ 0 : Z t ( R ) = 0 } .   T   is the time when all the particles disappear. Its distribution can be found easily. P { T ≤ t } = P { Z t ( R ) = 0 } = exp { − 2 z ¯ γ t } .   Let F   denote the location of the last particle immediately before extinction, i.e. { F } = S T −   .
We could recover the explicit distribution for F   .
Proposition 4.9. F   has the same distribution as X T   , where X   is a Brownian motion with initial distribution Z 0 / z ^   , and X   and T   are independent.
  • Proof. First assume that Z 0 = ∑ i = 1 m a i δ x i   with a i > 0   and ∑ i = 1 m a i = a   . Then Z t : = ∑ i = 1 m ξ i ( t ) δ X i ( t ) ,   where ξ i ( 0 ) = a i > 0   and ( X 1 , … , X m )   is a coalescing Brownian motion starting at ( x 1 , … , x m )   .
    Write T i : = inf { t ≥ 0 : ξ i ( t ) = 0 } , i = 1 , … , m   . Then T : = max 1 ≤ i ≤ m T i   . Therefore, F = ∑ i = 1 m X i ( T i − ) 1 { T = T i } = ∑ i = 1 m X i ( T ) 1 { T = T i } .   Our first observation is that P { T i ≤ t } = P { ξ i ( t ) = 0 } = exp { − 2 a i γ t } .   Then P { T = T i } = a i / a   , and F = X i ( T )   with probability a i / a   . Our second observation is that conditional on { T = T i }   , the distribution for T   is the same as its unconditional distribution. So, F   has the same distribution as the random variable obtained by running a Brownian motion X   with initial distribution P { X ( 0 ) = x i } = a i / a , i = 1 , … , m   , and stopping it independently at time T   . As a result, F   has the desired distribution.
    By conditioning on Z ε   and letting ε → 0 +   , the conclusion in the proposition also follows for any general initial measure Z 0   .
    □
Remark 4.10. This near extinction behavior is the same as the super Brownian motion (see Theorem 1 in [15] ).

5 Connections with the Arratia flow

Arratia flow is a stochastic flow which describes the evolution of a continuous family of coalescing Brownian motions on R   . We refer to [1] for a detailed account and [2] for a survey on stochastic flows. By definition, the Arratia flow is a collection { φ ( s , t , x ) : 0 ≤ s ≤ t , x ∈ R }   of random variables such that
  • ∙   the random map ( s , t , x ) ↦ φ ( s , t , x )   is jointly measurable,
  • ∙   for each s   and x   , the map t ↦ φ ( s , t , x )   , t ≥ s   , is continuous,
  • ∙   for each s   and t   with s ≤ t   , the map x ↦ φ ( s , t , x )   is non-decreasing and right-continuous,
  • ∙   for s ≤ t ≤ u   , φ ( t , u , ⋅ ) ∘ φ ( s , t , ⋅ ) = φ ( s , u , ⋅ )   ,
  • ∙   for u > 0   , ( s , t , x ) ↦ φ ( s + u , t + u , x )   has the same distribution as φ   ,
  • ∙   for x 1 < … < x m   the process ( φ ( 0 , t , x 1 ) , … , φ ( 0 , t , x m ) ) t ≥ 0   has the same distribution as a coalescing Brownian motion starting at ( x 1 , … , x m )   .
Fix t > 0   , it is known that { φ ( 0 , t , x ) : x ∈ R }   , the image of R   under map φ ( 0 , t , . )   , is a discrete set (see [1] ). Let … < x − 1 * < x 0 * < x 1 * < …   be a sequence of random variables such that
{ φ ( 0 , t , x ) : x ∈ R } = { … , x − 1 * , x 0 * , x 1 * , … } . (5.1)
Since Brownian motion has continuous sample paths, the Arratia flow is order-preserving; i.e.
φ ( 0 , t , x 1 ) ≤ … ≤ φ ( 0 , t , x m )   whenever x 1 ≤ , … , ≤ x m   . Set Π i : = sup { x : φ ( 0 , t , x ) = x i * } .   Write φ − 1 ( 0 , t , x )   for the pre-image of x   under φ ( 0 , t , . )   Then ( Π i )   determines a partition on R   such that φ − 1 ( 0 , t , x i * ) = [ Π i − 1 , Π i [   .
Not surprisingly, the Arratia flow is closely connected to the process Z   studied in the previous sections. We first consider its support S t   . Since S t   is a discrete set, we can identify it with a simple point process by placing a unit mass on each point of S t   . For any y 1 ≤ y 2 ≤ … ≤ y 2 n   , by Proposition  4.1 ,
P { Z t ( ∪ j = 1 n ] y 2 j − 1 , y 2 j ] ) = 0 } = P [ exp { − 2 γ t Z 0 ( ∪ j = 1 n ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] ) } ] . (5.2)
We thus get the following characterization of the avoidance function for S t   .
P { S t ∩ ∪ j = 1 n ] y 2 j − 1 , y 2 j ] = ∅ } = P [ exp { − 2 γ t Z 0 ( ∪ j = 1 n ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] ) } ] . (5.3)
Consequently the distribution of S t   is uniquely determined by ( 5.3 ); see Theorem 3.3 in [8] .
( 5.3 ) suggests a connection between S t   and the Arratia flow. Let M t ( d y )   be a random measure on R   such that M t ( ∪ j = 1 n ] y 2 j − 1 , y 2 j ] ) = 2 γ t ∑ j = 1 n Z 0 ( ] φ ( 0 , t , y 2 j − 1 ) , φ ( 0 , t , y 2 j ) ] ) , y 1 ≤ y 2 ≤ … ≤ y 2 n .   Then S t   can be identified with a Cox process with a finite random intensity measure M t   . ( 5.2 ) also leads to a result on the occupation time for Z   . For any Borel set B   in R   , ∫ 0 t d s P { Z s ( B ) = 0 } = ∫ 0 t d s P [ exp { − M t ( B ) } ] .   A particle representation for Z t   is available by using the image of the Arratia flow as a skeleton.
Given ( x i * )   as in ( 5.1 ), let ( … , ξ − 1 , ξ 0 , ξ 1 , … )   be independent non-negative random variables such that P [ exp { − λ ξ i } | ( x i * ) ] = exp { − 2 λ Z 0 ( φ − 1 ( 0 , t , x i * ) ) 2 + λ γ t } .   Then
Z t = D ∑ i = − ∞ ∞ ξ i δ x i * . (5.4)
To see this, define Z t ( m ) : = ∑ i = − m 2 m m 2 m ξ i ( m ) ( t ) δ φ ( 0 , t , i / 2 m ) ,   where ( ξ i ( m ) ) i = − m 2 m m 2 m   is a sequence of independent Feller's branching processes with initial values ( Z 0 ( [ ( i − 1 ) / 2 m , i / 2 m [ ) ) i = − m 2 m m 2 m   , and in addition, ( ξ i ( m ) )   is independent of { φ ( s , t , x ) }   .
For any a j ≥ 0 , j = 1 , … , n   , and y 1 ≤ … ≤ y 2 n   , by the same argument as in the proof for Theorem  3.2 , we have
lim m → ∞ P [ exp { − ∑ j = 1 n a j Z t ( m ) ( ] y 2 j − 1 , y 2 j ] ) } ] = lim m → ∞ P [ m 2 m ∏ i = − m 2 m exp { − 2 Z 0 ( [ ( i − 1 ) / 2 m , i / 2 m [ ) ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( i / 2 m ) 2 + γ t ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( i / 2 m ) } ] = lim m → ∞ P [ exp { − ∑ i = − m 2 m m 2 m 2 Z 0 ( [ ( i − 1 ) / 2 m , i / 2 m [ ) ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( i / 2 m ) 2 + γ t ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( i / 2 m ) } ] = P [ exp { − ∫ − ∞ ∞ Z 0 ( d x ) 2 ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) 2 + γ t ∑ j = 1 n a j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) } ] = P [ exp { − ∑ j = 1 n a j Z t ( ] y 2 j − 1 , y 2 j ] ) } ] .
Therefore, Z t ( m ) → D Z t .   Further, by the definition of ( x i * )   and the additive property for Feller's branching processes we obtain that Z t ( m ) = ∑ i = − ∞ ∞ ∑ Π i − 1 ≤ j / 2 m < Π i ξ j ( m ) ( t ) δ x i * → D ∑ i = − ∞ ∞ ξ i δ x i * .   Putting these together gives ( 5.4 ). This interplay is remarkable. On one hand, Z   can be constructed using the Arratia flow; on the other hand, Z   tells us how an initial measure Z 0   is transported over time under both the Arratia flow and the branching.
The Laplace functional for Z t   can also be expressed in terms of ( x i * )   and ( Π i )   . Given any nonnegative bounded continuous function f   , for y j = j / 2 n   , Theorem  3.2 yields
P [ exp { − ∑ j = − n 2 n n 2 n f ( y j ) Z t ( ] y j − 1 , y j ] ) } ] = P [ exp { − ∫ − ∞ ∞ Z 0 ( d x ) 2 ∑ j = − n 2 n n 2 n f ( y j ) 1 { ] φ ( 0 , t , y j − 1 ) , φ ( 0 , t , y j ) ] } ( x ) 2 + γ t ∑ j = − n 2 n n 2 n f ( y j ) 1 { ] φ ( 0 , t , y j − 1 ) , φ ( 0 , t , y j ) ] } ( x ) } ] . (5.5)
Again, let m → ∞   . It follows that
P [ exp { − λ ∫ − ∞ ∞ f ( x ) Z t ( d x ) } ] = P [ exp { − ∫ − ∞ ∞ Z 0 ( d x ) 2 λ ∑ i = − ∞ ∞ f ( Π i ) 1 { ] x i * , x i + 1 * ] } ( x ) 2 + λ γ t ∑ i = − ∞ ∞ f ( Π i ) 1 { ] x i * , x i + 1 * ] } ( x ) } ] .

6 A more general model

Evans observes that what is really at work in the proof for Theorem  3.2 is the additivity for the Feller's branching processes. He then suggested that we could use the square of Bessel processes (BESQ) to describe the evolution of masses. We are going to carry it out in this section.
For x ≥ 0   and δ ≥ 0   the square of δ   -dimensional Bessel process starting at x   , denoted by BESQ δ ( x )   , is a non-negative valued process ξ   which solves the following stochastic differential equation ξ t = x + 2 ∫ 0 t ξ s d B s + δ t ,   where B   is a one-dimensional Brownian motion. The Laplace transform for ξ   is given by
P [ exp { − λ ξ t } ] = 1 ( 1 + 2 λ t ) δ 2 exp { − λ x 1 + 2 λ t } . (6.1)
Notice that the Feller's branching process is just a BESQ 0   . We refer to Chapter XI in [11] for a more detailed introduction on the Bessel processes.
It is easy to see from ( 6.1 ) that BESQ δ ( x )   is additive in both δ   and x   ; i.e. if { ξ i , i = 1 , … , m }   is a sequence of independent processes such that each ξ i   is a BESQ δ i ( x i )   . Then ∑ i = 1 m ξ i   is a BESQ ∑ i = 1 m δ i ( ∑ i = 1 m x i )   .
Now we are going to modify the process Z   defined in Section 3 by letting the masses of the particles be governed by the BESQ processes. Since the dimension is an additional parameter for BESQ, we need to introduce another measure-valued process Δ   to describe the evolution of the dimension.
As in Section 3, we first consider two systems of interacting particles. Given a finite measure Z 0   on R   , for any m   , choose x 1 , … , x m   to be i.i.d. random variables with a common distribution Z ¯ 0 : = Z 0 / Z 0 ( R )   . Let ( X 1 , … , X m )   be an m   -dimensional coalescing Brownian motion starting at ( x 1 , … , x m )   .
Given another finite measure Δ 0   on R +   with a finite “mean” μ : = ∫ 0 ∞ δ Δ 0 ( d δ )   , let δ 1 , … , δ m   be i.i.d. random variables with a common distribution Δ ¯ : = Δ 0 / Δ 0 ( R )   . We further suppose that and ( X i )   and ( δ i )   are independent.
Put δ ¯ : = Δ 0 ( R )   and z ¯ : = Z 0 ( R )   . Let ( ξ 1 , … , ξ m )   be a collection of m   independent BESQ δ i δ ¯ / m ( z ¯ / m )   processes. Then Z t ( m ) : = ∑ i = 1 m ξ i ( t ) δ X i ( t )   and Δ t ( m ) : = 1 m ∑ i = 1 m δ i δ ¯ δ X i ( t )   define two M F ( R )   -valued processes.
Similar to Lemma  3.1 we can show that both { Z ( m ) }   and { Δ ( m ) }   are C-relatively compact in D ( M F ( R ) )   . They have unique weak limits by Theorem  6.1 , which we will prove shortly. Let Z   and Δ   be the weak limits for { Z ( m ) }   and { Δ ( m ) }   . Intuitively, { ( Z 0 ( B ) , Δ 0 ( B ) ) : B ∈ ℬ ( R ) }   describes the initial mass-dimension distribution on R   , and { ( Z t ( B ) , Δ t ( B ) ) : B ∈ ℬ ( R ) , 0 ≤ t < ∞ }   describes the simultaneous mass-dimension evolution for such a model, which we call a super square of Bessel process with spatial coalescing Brownian motion.
For any nonnegative constants α j , β j , j = 1 , … , n   and t > 0   , put I t ( x ) : = ∑ j = 1 n α j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x )   and J t ( x ) : = ∑ j = 1 n β j 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x ) ,   where, as usual, ( Y 1 , … , Y 2 n )   is an 2 n   -dimensional coalescing Brownian motion starting at ( y 1 , … , y 2 n )   . The next result determines the joint distribution for ( Z t ( B ) , Δ t ( B ) ) , B ∈ ℬ ( R )   .
Theorem 6.1. For any α j ≥ 0 , β j ≥ 0 , j = 1 , … , n   and any y 1 ≤ … ≤ y 2 n   , we have
P [ exp { − ∑ j = 1 n α j Z t ( ] y 2 j − 1 , y 2 j ] ) − ∑ j = 1 n β j Δ t ( ] y 2 j − 1 , y 2 j ] ) } ] = P [ exp { − ∫ − ∞ ∞ Z 0 ( d x ) ( μ 2 z ¯ ln ( 1 + 2 t I t ( x ) ) + I t ( x ) 1 + 2 t I t ( x ) + μ z ¯ J t ( x ) ) } ] . (6.2)
  • Proof. To prove ( 6.2 ), again, we first fix ( δ 1 , … , δ m )   and ( ξ 1 , … , ξ m )   . It follows that
    P [ exp { − ∑ j = 1 n α j Z t ( m ) ( ] y 2 j − 1 , y 2 j ] ) − ∑ j = 1 n β j Δ t ( m ) ( ] y 2 j − 1 , y 2 j ] ) } ] = P [ exp { − ∑ j = 1 n ( ∑ i = 1 m α j ξ i ( t ) 1 { ] y 2 j − 1 , y 2 j ] } ( X i ( t ) ) + ∑ i = 1 m β j δ i δ ¯ m 1 { ] y 2 j − 1 , y 2 j ] } ( X i ( t ) ) ) } ] = P [ exp { − ∑ j = 1 n ( ∑ i = 1 m α j ξ i ( t ) 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x i ) + ∑ i = 1 m β j δ i δ ¯ m 1 { ] Y 2 j − 1 ( t ) , Y 2 j ( t ) ] } ( x i ) ) } ] = P [ exp { − ∑ i = 1 m ξ i ( t ) I t ( x i ) − ∑ i = 1 m δ i δ ¯ m J t ( x i ) } ] . (6.3)
    We then fix ( x i )   and ( δ i )   , and take expectations with respect to ( ξ 1 , … , ξ m )   . By ( 6.1 ) the right hand side of ( 6.3 ) is equal to
    P [ m ∏ i = 1 ( 1 + 2 t I t ( x i ) ) − δ i δ ¯ 2 m exp { − z ¯ I t ( x i ) m ( 1 + 2 t I t ( x i ) ) − δ i δ ¯ J t ( x i ) m } ] = P [ ( ∫ 0 ∞ Z ¯ 0 ( d x ) ∫ 0 ∞ Δ ¯ 0 ( d δ ) ( 1 + 2 t I t ( x ) ) − δ δ ¯ 2 m exp { − z ¯ I t ( x ) m ( 1 + 2 t I t ( x ) ) − δ δ ¯ J t ( x ) m } ) m ] . (6.4)
    Let m → ∞   in ( 6.4 ). We finally have
    P [ exp { − ∑ j = 1 n α j Z t ( ] y 2 j − 1 , y 2 j ] ) − ∑ j = 1 n β j Δ t ( ] y 2 j − 1 , y 2 j ] ) } ] = lim m → ∞ P [ exp { − ∑ j = 1 n α j Z t ( m ) ( ] y 2 j − 1 , y 2 j ] ) − ∑ j = 1 n β j Δ t ( m ) ( ] y 2 j − 1 , y 2 j ] ) } ] = lim m → ∞ P [ { ∫ R × R + Z ¯ 0 ( d x ) Δ ¯ 0 ( d δ ) ( 1 − δ δ ¯ 2 m ln ( 1 + 2 t I t ( x ) ) ) ( 1 − z ¯ I t ( x ) m ( 1 + 2 t I t ( x ) ) − δ δ ¯ J t ( x ) m ) } m ] = lim m → ∞ P [ { ∫ − ∞ ∞ Z ¯ 0 ( d x ) ∫ 0 ∞ Δ ¯ 0 ( d δ ) ( 1 − δ δ ¯ 2 m ln ( 1 + 2 t I t ( x ) ) − z ¯ I t ( x ) m ( 1 + 2 t I t ( x ) ) − δ δ ¯ J t ( x ) m ) } m ] = P [ exp { − ∫ − ∞ ∞ Z ¯ 0 ( d x ) ∫ 0 ∞ Δ ¯ 0 ( d δ ) ( δ δ ¯ 2 ln ( 1 + 2 t I t ( x ) ) + z ¯ I t ( x ) 1 + 2 t I t ( x ) + δ δ ¯ J t ( x ) ) } ] = P [ exp { − ∫ − ∞ ∞ Z 0 ( d x ) ( μ 2 z ¯ ln ( 1 + 2 t I t ( x ) ) + I t ( x ) 1 + 2 t I t ( x ) + μ z ¯ J t ( x ) ) } ] .
    □
Remark 6.2. Notice that Δ   is just the process Z   in Theorem  3.2 with γ = 0   .
The generalized model considered in this section will not die out for μ > 0   . Many of the properties in Section 3 and 4 can be discussed in a similar fashion. But we leave the details to the interested readers.
Acknowledgement: The author is grateful to Steven Evans for a suggestion that results in Section 5 of this paper. The author also thanks Carl M u ¨   ller for a helpful comment.
References

  1. R. Arratia, Coalescing Brownian motions on the line, Ph.D. thesis, University of Wisconsin, Madison 1979.
  2. R. W. R. Darling, Isotropic stochastic flows: a survey, Diffusion Processes and Related Problems in Analysis (M. Pinsky and M. Wihstutz, eds.) 2 75–94. Birkh a ¨   user, Boston, 1992.
  3. D. A. Dawson, Z. H. Li, H. Wang, Superprocesses with dependent spatial motion and general branching densities, Elect. J. Probab. 6 (2001) 25, 1–33.
  4. D. A. Dawson, Z. H. Li, Construction of immigration superprocesses with dependent spatial motion from one–dimensional excursions, Probab. Theory Related Fields, 127 (1) (2003), 37–61.
  5. D. A. Dawson, Z. H. Li, X. Zhou, Rescaled limit of a superprocess with dependent spatial motion, J. Theo. Probab. 17 (3) (2004), 673–692.
  6. S. N. Ethier, T. G. Kurtz, Markov Processes: Characterization and Convergence, Wiley, New York, 1986.
  7. J. Jacod, A.N. Shiryaev, Limit Theorems for Stochastic Processes, Springer-Verlag, Berlin, 1987.
  8. O. Kallenberg, Random Measures, Academic Press, New York, 1976.
  9. Z. Ma, K. Xiang, Superprocesses of stochastic flows, Ann. Probab. 29 (1) (2001) 317–343.
  10. E. Perkins, Dawson-Watanabe Superprocesses and Measure-valued Diffusions, Lectures on probability theory and statistics (Saint-Flour, 1999), 125–329, Lecture Notes in Math., 1781, Springer, Berlin, 2002.
  11. D. Revuz, M. Yor, Continuous Martingales and Brownian Motion, Springer, Berlin, 1991.
  12. G. Skoulakis, R. J. Adler, Superprocesses over a stochastic flow, Ann. Appl. Probab. 11 (2) (2002) 488–543.
  13. F. Soucaliuc, B. Tóth, W. Werner, Reflection and coalescence between independent one-dimensional Brownian paths, Ann. Inst. H. Poincaré Probab. Statist. 36 (4) (2000) 509–545.
  14. B. Tóth, W. Werner, The True self–repelling motion, Probab. Theory Related Fields, 111 (3) (1997) 375–452.
  15. R. Tribe, The behavior of superprocesses near extinction, Ann. Probab. 20 (1) (1992) 286–311.
  16. J. Xiong, X. Zhou, On the duality between coalescing Brownian motions, Canad. J. Math. 57 (1) (2005) 204–224.