Misplaced Pages

Brownian snake

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
Stochastic Markov process

A Brownian snake is a stochastic Markov process on the space of stopped paths. It has been extensively studied., and was in particular successfully used as a representation of superprocesses.

Informally, superprocesses are the scaling limit of branching processes, except each particle splits and dies at infinite rates. The Brownian snake is a stochastic object that enables the representation of the genealogy of a superprocess, providing a link between super-Brownian motion and Brownian trees. In other words, even though infinitely many particles are constantly born, we can still keep track of individual trajectories in space, or of when two given present-day particles have split from a common ancestor in the past.

History

The Brownian snake approach was originally developed by Jean-François Le Gall. It has since been applied in fragmentation theory, partial differential equation or planar map

The simplest setting

Let D ( R + , R ) {\displaystyle D(\mathbb {R} _{+},\mathbb {R} )} be the space of càdlàg functions from R + {\displaystyle \mathbb {R} _{+}} to R {\displaystyle \mathbb {R} } , equipped with a metric d {\displaystyle d} compatible with the Skorokhod topology. We define a stopped path as a couple ( w , z ) {\displaystyle (w,z)} where w D ( R + , R ) {\displaystyle w\in D(\mathbb {R} _{+},\mathbb {R} )} and z R + {\displaystyle z\in \mathbb {R} _{+}} are such that w ( t ) = w ( t z ) {\displaystyle w(t)=w(t\wedge z)} . In other words, w {\displaystyle w} is constant after z {\displaystyle z} .

Now, we consider a jump process ( J s N ) s 0 {\displaystyle (J_{s}^{N})_{s\geq 0}} with states { + 1 , 1 } {\displaystyle \{+1,-1\}} and jump rate N {\displaystyle N} , such that J 0 N = + 1 {\displaystyle J_{0}^{N}=+1} . We set: β ^ s N := 0 s J s N d s {\displaystyle {\hat {\beta }}_{s}^{N}:=\int _{0}^{s}J_{s'}^{N}ds'} and then β s N := | β ^ s N | {\displaystyle \beta _{s}^{N}:=|{\hat {\beta }}_{s}^{N}|} to be the process reflected on 0.

In words, β s N {\displaystyle \beta _{s}^{N}} increases with speed 1, until J s N {\displaystyle J_{s}^{N}} jumps, in which case it decreases with speed 1, and so on. We define the stopping time σ N {\displaystyle \sigma _{N}} to be the N {\displaystyle N} -th hitting time of 0 by β N {\displaystyle \beta ^{N}} . We now define a stochastic process ( η s N , β s N ) s R + {\displaystyle (\eta _{s}^{N},\beta _{s}^{N})_{s\in \mathbb {R} _{+}}} on the set of stopped paths as follows:

  • η 0 N = 0 {\displaystyle \eta _{0}^{N}=0}
  • if J s N = + 1 {\displaystyle J_{s}^{N}=+1} for s [ s 1 , s 2 ] {\displaystyle s\in } then:
    • η s 1 ( t ) = η s 2 ( t ) {\displaystyle \eta _{s_{1}}(t)=\eta _{s_{2}}(t)} for t β s 1 {\displaystyle t\leq \beta _{s_{1}}}
    • ( η s 2 ( t β s 1 ) η s 1 ( β s 1 ) ) 0 t β s 2 β s 1 {\displaystyle {\Big (}\eta _{s_{2}}(t-\beta _{s_{1}})-\eta _{s_{1}}(\beta _{s_{1}}){\Big )}_{0\leq t\leq \beta _{s_{2}}-\beta _{s_{1}}}} is distributed as a Brownian motion independent from η s 1 {\displaystyle \eta _{s_{1}}}
  • if J s N = 1 {\displaystyle J_{s}^{N}=-1} for s [ s 1 , s 2 ] {\displaystyle s\in } then η s 1 ( t ) = η s 2 ( t ) {\displaystyle \eta _{s_{1}}(t)=\eta _{s_{2}}(t)} for t β s 2 {\displaystyle t\leq \beta _{s_{2}}}

See animation for an illustration. We call this process a snake and β s N {\displaystyle \beta _{s}^{N}} the head of the snake. This process is not yet the Brownian snake, but a good introduction. The path is erased when the snake head moves backwards, and is created anew when it moves forward.

Show/hide animation
The left panel is the "health bar", which goes to 0 as the number of times the snake hits 0 increases. The panel below just shows how the head of the snake moves. The large panel represents the non-constant part of the snake in black, with the head as a red dot.

Duality with a branching Brownian motion

We now consider a measure-valued branching process ( X t N ) t 0 {\displaystyle (X_{t}^{N})_{t\geq 0}} starting with N {\displaystyle N} particles, such that each particle dies with rate N {\displaystyle N} , and upon its death gives birth to two offspring with probability 1 / 2 {\displaystyle 1/2} .

On the other hand, we may define from our process ( η s N , β s N ) 0 s σ N {\displaystyle (\eta _{s}^{N},\beta _{s}^{N})_{0\leq s\leq \sigma ^{N}}} a measure-valued random process X ^ t {\displaystyle {\hat {X}}_{t}} as follows: note that for any t R + {\displaystyle t\in \mathbb {R} _{+}} , there will almost surely be finitely many times s 1 , s 2 , , s n [ 0 , σ N ] {\displaystyle s_{1},s_{2},\dots ,s_{n}\in } such that β s i = t {\displaystyle \beta _{s_{i}}=t} . We then set for any measurable function f {\displaystyle f} :

X ^ t N ( f ) := i = 1 n f ( η s N ( t ) ) {\displaystyle {\hat {X}}_{t}^{N}(f):=\sum \limits _{i=1}^{n}f(\eta _{s}^{N}(t))}

Then X {\displaystyle X} and X ^ {\displaystyle {\hat {X}}} are equal in distribution.

The Brownian snake

We take the limit of the previous system as N {\displaystyle N\to \infty } . In this setting, the head of the snake keeps jittering. In fact, the process β s N {\displaystyle \beta _{s}^{N}} tends towards a reflected Brownian motion β s {\displaystyle \beta _{s}} . The definitions are no longer valid for a number of reasons, in particular because β s {\displaystyle \beta _{s}} is almost surely never monotonous on any interval.

However, we may define a probability R a , b ( ( u , y ) , d ( w , z ) ) {\displaystyle R_{a,b}((u,y),d(w,z))} on stopped paths such that:

  • R a , b {\displaystyle R_{a,b}} -almost surely z = b {\displaystyle z=b} and w ( t ) = u ( t ) {\displaystyle w(t)=u(t)} for 0 t a {\displaystyle 0\leq t\leq a}
  • The law of ( w ( a + t ) ) 0 t b a {\displaystyle (w(a+t))_{0\leq t\leq b-a}} is the law of a standard Brownian motion.

We may also define γ s y ( d a , d b ) {\displaystyle \gamma _{s}^{y}(da,db)} to be the distribution of ( inf 0 r s β r , β s ) {\displaystyle (\inf _{0\leq r\leq s}\beta _{r},\beta _{s})} if β 0 = y {\displaystyle \beta _{0}=y} . Finally, define the transition semigroup on the set of stopped paths:

Q s ( ( u , y ) , d ( w , z ) ) = γ s y ( d a , d b ) R a , b ( ( u , y ) , d ( w , z ) ) {\displaystyle Q_{s}((u,y),d(w,z))=\int \gamma _{s}^{y}(da,db)R_{a,b}((u,y),d(w,z))}

A stochastic process with this semigroup is called a Brownian snake.

We may again find a duality between this process and a branching process. Here the branching process will be a super-Brownian motion ( X t ) t R + {\displaystyle (X_{t})_{t\in \mathbb {R} _{+}}} with branching mechanism ϕ ( z ) = z 2 {\displaystyle \phi (z)=z^{2}} , started on a Dirac in 0.

However, unlike the previous case, we must be more careful in the definition of the process X ^ {\displaystyle {\hat {X}}} . Indeed, for t R + {\displaystyle t\in \mathbb {R} _{+}} we cannot just list the times s 1 , s 2 , {\displaystyle s_{1},s_{2},\dots } such that β s = t {\displaystyle \beta _{s}=t} . Instead we use the local time l s ( t ) {\displaystyle l_{s}(t)} associated with β s {\displaystyle \beta _{s}} : we first define the stopping time σ = inf { s 0 , l s ( 0 ) u } {\displaystyle \sigma =\inf\{s\geq 0,l_{s}(0)\geq u\}} . Then we define for any measurable f {\displaystyle f} : X ^ t ( f ) := 0 σ f ( η s ( t ) ) d l s ( t ) {\displaystyle {\hat {X}}_{t}(f):=\int _{0}^{\sigma }f(\eta _{s}(t))dl_{s}(t)} Then, as before, we obtain that X {\displaystyle X} and X ^ {\displaystyle {\hat {X}}} are equal in distribution. See the animation for the construction of the branching process from the Brownian snake.

Animation for the branching process associated with the Brownian snake
The left panel shows the "health bar" of the snake, which decreases with the local time the head spends on 0. The panel below shows the movement of the snake head according to a Brownian motion reflected on 0. The central panel shows: in red the head of the current snake, in black the current snake, in green the past snakes. The branching superprocess X {\displaystyle X} is obtained once the health bar reaches 0, by taking all of the green paths.

Generalisation

The previous example can be generalized in many ways:

  • We may consider D ( R + , E ) {\displaystyle D(\mathbb {R} _{+},E)} where ( E , d ) {\displaystyle (E,d)} is a complete separable metric space.
  • Instead of a Brownian motion, the underlying movement of the snake can be very general class of Markov processes (see Superprocess).

Link with genealogy and the Brownian tree

The Brownian snake can be seen as a way to represent the genealogy of a superprocess, the same way a Galton-Watson tree may encode the hidden genealogy of a Galton–Watson process. Indeed, for two points of the Brownian snake, their common ancestor will be the infimum of the snake's head position between them.

If we take a Brownian snake and construct a real tree from it, we obtain a Brownian tree.

References

  1. Li, Zenghu (2011), Li, Zenghu (ed.), "Measure-Valued Branching Processes", Measure-Valued Branching Markov Processes, Probability and Its Applications, Berlin, Heidelberg: Springer, pp. 29–56, doi:10.1007/978-3-642-15004-3_2, ISBN 978-3-642-15004-3, retrieved 2022-12-20
  2. ^ Le Gall, Jean-Francois (1999-07-01). Spatial Branching Processes, Random Snakes and Partial Differential Equations. Springer Science & Business Media. ISBN 978-3-7643-6126-6.
  3. Le Gall, Jean-Francois (1991). "Brownian Excursions, Trees and Measure-Valued Branching Processes". The Annals of Probability. 19 (4): 1399–1439. ISSN 0091-1798. JSTOR 2244522.
  4. Abraham, Romain; Serlet, Laurent (2002-07-01). "Poisson Snake and Fragmentation". Electronic Journal of Probability. 7 (none). doi:10.1214/EJP.v7-116. ISSN 1083-6489. S2CID 12003800.
  5. Abraham, Romain (2000-10-01). "Reflecting Brownian snake and a Neumann–Dirichlet problem". Stochastic Processes and Their Applications. 89 (2): 239–260. doi:10.1016/S0304-4149(00)00027-2. ISSN 0304-4149.
  6. Le Gall, Jean-François (2019-09-01). "Brownian geometry". Japanese Journal of Mathematics. 14 (2): 135–174. doi:10.1007/s11537-019-1821-7. ISSN 1861-3624. S2CID 255314865.
  7. Miermont, Grégory (2013). "The Brownian map is the scaling limit of uniform random plane quadrangulations". Acta Mathematica. 210 (2): 319–401. arXiv:1104.1606. doi:10.1007/s11511-013-0096-8. ISSN 0001-5962. S2CID 119140342.
Category: