Branching Brownian motion seen from the tip

Size: px
Start display at page:

Download "Branching Brownian motion seen from the tip"

Transcription

1 Branching Brownian motion seen from the tip J. Berestycki 1 1 Laboratoire de Probabilité et Modèles Aléatoires, UPMC, Paris 09/02/2011 Joint work with Elie Aidekon, Eric Brunet and Zhan Shi J Berestycki (Paris) ANR MANEGE 09/02/ / 35

2 Outline 1 BBM and FKPP J Berestycki (Paris) ANR MANEGE 09/02/ / 35

3 Outline 1 BBM and FKPP 2 BBM seen from the tip J Berestycki (Paris) ANR MANEGE 09/02/ / 35

4 Outline 1 BBM and FKPP 2 BBM seen from the tip J Berestycki (Paris) ANR MANEGE 09/02/ / 35

5 BBM The model = Particles X 1 (t),..., X N(t) (t) on R Start with one particle at 0. Movement = independent Brownian motions Branching = at rate 1 into two new particles (more general possible). J Berestycki (Paris) ANR MANEGE 09/02/ / 35

6 Rightmost particle Define M(t) = max i=1,...,n(t) X i (t). Theorem (Rightmost particle M(t)) c = 2. lim t M(t) t = c. c t M(t) a.s. Proof by martingale techniques J Berestycki (Paris) ANR MANEGE 09/02/ / 35

7 Rightmost particle Define M(t) = max i=1,...,n(t) X i (t). Define u(t, x) := P(M(t) < x). Theorem (Rightmost particle M(t)) c = 2. lim t M(t) t = c. c t M(t) a.s. Proof by martingale techniques J Berestycki (Paris) ANR MANEGE 09/02/ / 35

8 Rightmost particle Define M(t) = max i=1,...,n(t) X i (t). Theorem (Rightmost particle M(t)) c = 2. lim t M(t) t = c. c t M(t) a.s. Proof by martingale techniques Define u(t, x) := P(M(t) < x). The map u(t, x) = P(M(t) < x) solves t u = x u + u(u 1) ( avec 1 u(x, 0) = ) J Berestycki (Paris) ANR MANEGE 09/02/ / 35

9 Rightmost particle Define M(t) = max i=1,...,n(t) X i (t). Theorem (Rightmost particle M(t)) c = 2. lim t M(t) t = c. c t M(t) a.s. Proof by martingale techniques Define u(t, x) := P(M(t) < x). The map u(t, x) = P(M(t) < x) solves t u = x u + u(u 1) ( avec 1 u(x, 0) = ) Idea : Initial particle in 0. After dt - with proba (1 dt) no split but diffuse : u(t + dt, x) u(t, x ξ t ) - with proba dt branch u(t + dt, x) u(t, x) 2 J Berestycki (Paris) ANR MANEGE 09/02/ / 35

10 Bramson s result Define m t by u(t, m t ) = 1/2. i.e. m t is the median of M(t) (results still valid if we take the expectation). KPP 37 u(t, m t + x) w(x) unif. in x as t where m t = 2t + a(t) and a(t). (McKean shows that a(t) 2 3/2 log t) and w solution to 1 2 w + 2w + (w 2 w) = 0. Purely analytical method. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

11 Bramson s result Define m t by u(t, m t ) = 1/2. i.e. m t is the median of M(t) (results still valid if we take the expectation). KPP 37 u(t, m t + x) w(x) unif. in x as t where m t = 2t + a(t) and a(t). (McKean shows that a(t) 2 3/2 log t) and w solution to 1 2 w + 2w + (w 2 w) = 0. Purely analytical method. Bramson 83 m t = 2t 3 2 3/2 log t + c + o t (1) J Berestycki (Paris) ANR MANEGE 09/02/ / 35

12 Bramson s result Define m t by u(t, m t ) = 1/2. i.e. m t is the median of M(t) (results still valid if we take the expectation). KPP 37 u(t, m t + x) w(x) unif. in x as t where m t = 2t + a(t) and a(t). (McKean shows that a(t) 2 3/2 log t) and w solution to 1 2 w + 2w + (w 2 w) = 0. Purely analytical method. Bramson 83 m t = 2t 3 log t + c + o 2 3/2 t (1) ( ) ( ) 1 1 If change initial condition from, to Then only the constant c changes J Berestycki (Paris) ANR MANEGE 09/02/ / 35

13 Bramson s result Define m t by u(t, m t ) = 1/2. i.e. m t is the median of M(t) (results still valid if we take the expectation). KPP 37 u(t, m t + x) w(x) unif. in x as t where m t = 2t + a(t) and a(t). (McKean shows that a(t) 2 3/2 log t) and w solution to 1 2 w + 2w + (w 2 w) = 0. Purely analytical method. Bramson 83 m t = 2t 3 log t + c + o 2 3/2 t (1) ( ) ( ) 1 1 If change initial condition from, to Then only the constant c changes Hence C B R such that with P(W x) = w(x). M t ( 2t log t + C B) dist W 0 J Berestycki (Paris) ANR MANEGE 09/02/ /

14 Lalley-Sellke KPP, Bramson P(M t m t < x) w(x) so M t m t converges in law to a variable with dist. w. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

15 Lalley-Sellke KPP, Bramson P(M t m t < x) w(x) so M t m t converges in law to a variable with dist. w. Can t be ergodic i.e. 1 t 1 Ms<m t s+xds w(x) a.s.. 0 J Berestycki (Paris) ANR MANEGE 09/02/ / 35

16 Lalley-Sellke KPP, Bramson P(M t m t < x) w(x) so M t m t converges in law to a variable with dist. w. Can t be ergodic i.e. 1 t 1 Ms<m t s+xds w(x) a.s.. 0 Suppose it holds. Then t 1 t 0 1 M x s <m s+xds w(0) a.s.. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

17 Lalley-Sellke KPP, Bramson P(M t m t < x) w(x) so M t m t converges in law to a variable with dist. w. Can t be ergodic i.e. 1 t 1 Ms<m t s+xds w(x) a.s.. 0 Suppose it holds. Then t 1 t 0 1 M x s <m s+xds w(0) a.s.. Start two independent BBM, one at x and one at 0. Positive probability that they meet before any branching successfull coupling so w(0) = w(x). Contradiction. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

18 The derivative martingale The derivative martingale Z(t) = u N(t) ( 2t X u (t))e 2Xu(t) 2t (additive martingale W 2 (t) = u N(t) e 2X u(t) 2t 0) Z := lim Z(t) exists, is finite and positive with proba. 1. Not UI. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

19 The derivative martingale The derivative martingale Z(t) = u N(t) ( 2t X u (t))e 2Xu(t) 2t (additive martingale W 2 (t) = u N(t) e 2X u(t) 2t 0) Z := lim Z(t) exists, is finite and positive with proba. 1. Not UI. Theorem C > 0 s.t. x R lim lim P (M(t + s) m(t + s) x F s) = exp{ CZe 2x }, a.s. s t J Berestycki (Paris) ANR MANEGE 09/02/ / 35

20 The derivative martingale The derivative martingale Z(t) = u N(t) ( 2t X u (t))e 2Xu(t) 2t (additive martingale W 2 (t) = u N(t) e 2X u(t) 2t 0) Z := lim Z(t) exists, is finite and positive with proba. 1. Not UI. Theorem C > 0 s.t. x R lim lim P (M(t + s) m(t + s) x F s) = exp{ CZe 2x }, a.s. s t Thus the TW has representation [ ] w(x) = E exp{ CZe 2x } J Berestycki (Paris) ANR MANEGE 09/02/ / 35

21 The derivative martingale The derivative martingale Z(t) = u N(t) ( 2t X u (t))e 2Xu(t) 2t (additive martingale W 2 (t) = u N(t) e 2X u(t) 2t 0) Z := lim Z(t) exists, is finite and positive with proba. 1. Not UI. Theorem C > 0 s.t. x R lim lim P (M(t + s) m(t + s) x F s) = exp{ CZe 2x }, a.s. s t Thus the TW has representation 1 w(x) Cxe 2x. [ ] w(x) = E exp{ CZe 2x } J Berestycki (Paris) ANR MANEGE 09/02/ / 35

22 Structure of W Suggests that, once we "know" Z, P (M(t) m(t) x) exp{ e 2x+log CZ }, a.s. so that P( 2(M(t) m(t)) log(cz) x) exp( e x ) and hence M(t) m(t) 2 1/2 log(cz) dist 2 1/2 G where G is a Gumbel variable. W = dist (G + log(cz))/ 2 J Berestycki (Paris) ANR MANEGE 09/02/ / 35

23 Structure of W Suggests that, once we "know" Z, P (M(t) m(t) x) exp{ e 2x+log CZ }, a.s. so that P( 2(M(t) m(t)) log(cz) x) exp( e x ) and hence M(t) m(t) 2 1/2 log(cz) dist 2 1/2 G where G is a Gumbel variable. W = dist (G + log(cz))/ 2 Purple line is m(t) J Berestycki (Paris) ANR MANEGE 09/02/ / 35

24 Structure of W Suggests that, once we "know" Z, P (M(t) m(t) x) exp{ e 2x+log CZ }, a.s. so that P( 2(M(t) m(t)) log(cz) x) exp( e x ) and hence M(t) m(t) 2 1/2 log(cz) dist 2 1/2 G where G is a Gumbel variable. W = dist (G + log(cz))/ 2 Purple line is m(t) Because of early fluctuation a random delay is created (2 1/2 log(cz)) J Berestycki (Paris) ANR MANEGE 09/02/ / 35

25 Structure of W Suggests that, once we "know" Z, P (M(t) m(t) x) exp{ e 2x+log CZ }, a.s. so that P( 2(M(t) m(t)) log(cz) x) exp( e x ) and hence M(t) m(t) 2 1/2 log(cz) dist 2 1/2 G where G is a Gumbel variable. W = dist (G + log(cz))/ 2 Purple line is m(t) Because of early fluctuation a random delay is created (2 1/2 log(cz)) Fluctuations around m(t) + 2 1/2 log(cz) are Gumbel. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

26 Outline 1 BBM and FKPP 2 BBM seen from the tip J Berestycki (Paris) ANR MANEGE 09/02/ / 35

27 Conjecture Let X 1 (t) > X 2 (t) >... be the particles ordered by position at time t. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

28 Conjecture Let X 1 (t) > X 2 (t) >... be the particles ordered by position at time t. Lalley and Sellke conjecture that the point process of particles converges to an equilibrium state. X i (t) m(t) + c log(cz(t))/ 2 J Berestycki (Paris) ANR MANEGE 09/02/ / 35

29 Conjecture Let X 1 (t) > X 2 (t) >... be the particles ordered by position at time t. Lalley and Sellke conjecture that the point process of particles X i (t) m(t) + c log(cz(t))/ 2 converges to an equilibrium state. First we show a simple argument due to Brunet Derrida to show that the PP X i (t) m(t) converges. Uses Bramson and McKean representation. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

30 McKean Representation Recall : u(t, x) := P(M(t) < x) solves t u = x u + u(u 1) with ( ) 1 1 u(x, 0) =. More generally, let g : R [0, 1] then Theorem (McKean, 1975) 0 If u : R R + [0, 1] solves the FKPP equation 0 t u = x u + u(u 1) with initial condition u(0, x) = g(x), then u(t, x) = E[ u N(t) g(x u (t) + x)]. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

31 McKean Representation Recall : u(t, x) := P(M(t) < x) solves t u = x u + u(u 1) with ( ) 1 1 u(x, 0) =. More generally, let g : R [0, 1] then Theorem (McKean, 1975) 0 If u : R R + [0, 1] solves the FKPP equation 0 t u = x u + u(u 1) with initial condition u(0, x) = g(x), then u(t, x) = E[ u N(t) In general t u = x u + β(f (u) u). g(x u (t) + x)]. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

32 Brunet Derrida argument (seen from m t ) 1 0 exp( k) H φ (x, t) = E[ φ(x X u (t))], H φ solves KPP with H φ (x, 0) = φ(x). If φ H φ (x, t) = P(M(t) < x) Bramson : H φ (m(t) + z, t) w(z + δ(φ)). If φ H φ (x, t) = E[e kn(x,t) ] with N(x, t) = #u : X u (t) > x. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

33 Brunet Derrida argument (seen from m t ) 1 0 exp( k) H φ (x, t) = E[ φ(x X u (t))], H φ solves KPP with H φ (x, 0) = φ(x). If φ H φ (x, t) = P(M(t) < x) If φ H φ (x, t) = E[e kn(x,t) ] with N(x, t) = #u : X u (t) > x. Bramson : H φ (m(t) + z, t) w(z + δ(φ)). For any Borel set A R, the Laplace transform of #{u : X u (t) m(t) + A} converges. Theorem (Brunet Derrida 2010) The point process of the particles seen from m(t) converges in distribution as t. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

34 Not too hard to show : the limit point process (X i, i = 1, 2,...) has the superposition property, i.e., α, β s.t. e α + e β = 1 : X α + X β = d X. Clearly the PPP (e x ) has this property. Who else? J Berestycki (Paris) ANR MANEGE 09/02/ / 35

35 Not too hard to show : the limit point process (X i, i = 1, 2,...) has the superposition property, i.e., α, β s.t. e α + e β = 1 : X α + X β = d X. Clearly the PPP (e x ) has this property. Who else? each atom is replaced by an iid copy of a point process Maillard (2011) show those are the only superposable PP. What is the decoration of BBM? J Berestycki (Paris) ANR MANEGE 09/02/ / 35

36 Not too hard to show : the limit point process (X i, i = 1, 2,...) has the superposition property, i.e., α, β s.t. e α + e β = 1 : X α + X β = d X. Clearly the PPP (e x ) has this property. Who else? each atom is replaced by an iid copy of a point process Maillard (2011) show those are the only superposable PP. What is the decoration of BBM? Independently, Bovier, Arguin and Kestler 2010 obtain results that are similar to part of what follows. Seem to use techniques. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

37 Normalization We do the following change of coordinates : we suppose that the Brownian motions have diffusion σ 2 = 2 and drift ρ = 2. Instead of rightmost focus on leftmost. Tilts the cone in which the BBM lives. Left-most particle m(t) = 3 2 log t + C B. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

38 Normalization We do the following change of coordinates : we suppose that the Brownian motions have diffusion σ 2 = 2 and drift ρ = 2. Instead of rightmost focus on leftmost. Tilts the cone in which the BBM lives. Left-most particle m(t) = 3 2 log t + C B. In this framework Z(t) := u N(t) X u (t)e Xu(t) is the derivative martingale. Recall its limit exists and is positive a.s. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

39 BBM seen from the tip Y i (t) := X i (t) m t + log(cz), 1 i N(t). Theorem (Aidekon, B., Brunet, Shi 11) As t the point process (Y i (t), 1 i N(t)) converges in distribution to the point process L obtained as follows (i) Define P a Poisson point process on R +, with intensity measure 2 σ e 2x/σ dx where a is a constant. (ii) For each atom x of P, we attach a point process x + Q (x) where Q (x) are i.i.d. copies of a certain point process Q. (iii) L is then the superposition of all the point processes x + Q (x) : L := {x + y : x P {0}, y Q (x) } Bovier Arguin and Kistler (2010) obtain a very similar result (and much more). J Berestycki (Paris) ANR MANEGE 09/02/ / 35

40 Structure of Q τ X 1,t 2 τ X 1,t 1 t X 1,t N τ 1 t,x 1,t N τ 2 t,x 1,t J Berestycki (Paris) ANR MANEGE 09/02/ / 35

41 Structure of Q τ X 1,t 2 τ X 1,t 1 Y (s) := X 1,t(t s) X 1,t(t) is the red path. Q is the superposition of the BBM that branch on this path. t X 1,t N τ 1 t,x 1,t N τ 2 t,x 1,t J Berestycki (Paris) ANR MANEGE 09/02/ / 35

42 Structure of Q τ X 1,t 2 τ X 1,t 1 Y (s) := X 1,t(t s) X 1,t(t) is the red path. Q is the superposition of the BBM that branch on this path. t X 1,t N τ 1 t,x 1,t N τ 2 t,x 1,t N x (t) = BBM at time t started from one ptc at x. Nx (t) = BBM at time t conditioned to min Nx (t) > 0 started from one ptc at x. G t (x) := P{min N 0 (t) x}, so that G t (x + m t ) P(σW x) by Bramson. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

43 Law of Y. U (b) v := B v, if v [0, T b ], b R v Tb, if v T b. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

44 Law of Y. U (b) v := B v, if v [0, T b ], b R v Tb, if v T b. Take b random : P(b dx) = f (x)/c 1 where f (x) := E [e 2λ 0 ] G v (σu v (x) )dv (1) and c 1 = constant. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

45 Law of Y. U (b) v := B v, if v [0, T b ], b R v Tb, if v T b. Take b random : P(b dx) = f (x)/c 1 where f (x) := E [e 2λ 0 ] G v (σu v (x) )dv (1) and c 1 = constant. Conditionally on b, Y has a density /U (b) given by 1 f (b) e 2λ 0 G v (σu v (b) )dv (2) J Berestycki (Paris) ANR MANEGE 09/02/ / 35

46 Structure of Q Theorem (Aidekon, B., Brunet, Shi 11) Let (Y (t), t 0) be as above. Start independent BBMs N Y (t) (t) conditioned to finish to the right of X 1,t along the path Y at rate 2λ(1 P(min N Y (t) (t) 0))dt. Then t π N Y (t) (t) is distributed as Q. t X 1,t τ X 1,t 1 τ X 1,t 2 N τ 1 t,x 1,t N τ 2 t,x 1,t D(ζ, t) := τ X N (i) 1,t i t ζ t,x 1,t particles born off X 1,t less than ζ unit of time ago, then we have the following joint convergence in distribution lim lim {(X 1,t(t s) X 1,t (t), s 0), D(ζ, t)} = {(Y (s), s 0), Q}. ζ t J Berestycki (Paris) ANR MANEGE 09/02/ / 35

47 Structure of Q { ( I ζ (t) = E exp i 1 )} nj=1 t τ i <ζ α j #[N τ i t,x 1,t (X 1 (t) + A j )] J Berestycki (Paris) ANR MANEGE 09/02/ / 35

48 Structure of Q { ( I ζ (t) = E exp i 1 )} nj=1 t τ i <ζ α j #[N τ i t,x 1,t (X 1 (t) + A j )] Theorem (Aidekon, B., Brunet, Shi 11) α j 0 (for 1 j n), E {e n j=1 α jq(a j ) } = lim lim I 0 ζ(t) = ζ t 0 E(e 2λ 0 E(e 2λ 0 G v (σu(b) v )dv )db G v (σu v (b) )dv )db, where G v (x) := P{min N (v) x}, Gv (x) := 1 E [e n j=1 α j#[n (v) (x+a j )] 1{min N (v) x} ]. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

49 Define { I(t) := E F (X 1,t (s), s [0, t]) ( exp i f (t τ X 1,t i ) n j=1 )} α j #[N τ i t,x 1,t (X 1 (t) + A j )], τ X 1,t 2 τ X 1,t 1 t X 1,t N τ 1 t,x 1,t N τ 2 t,x 1,t J Berestycki (Paris) ANR MANEGE 09/02/ / 35

50 Define { I(t) := E F (X 1,t (s), s [0, t]) ( exp i f (t τ X 1,t i ) n j=1 )} α j #[N τ i t,x 1,t (X 1 (t) + A j )], τ X 1,t 2 r 0, x R define τ X 1,t 1 G r (x) := P{min N (r) x}, n ] G (f ) r (x) := E [e f (r) α j #[N (r) (x+a j )] j=1 1{min N (r) x}. Hence, when f 0 we have G (f ) r (x) = 1 G r (x)). t X 1,t N τ 1 t,x 1,t N τ 2 t,x 1,t J Berestycki (Paris) ANR MANEGE 09/02/ / 35

51 Define { I(t) := E F (X 1,t (s), s [0, t]) ( exp i f (t τ X 1,t i ) n j=1 )} α j #[N τ i t,x 1,t (X 1 (t) + A j )], τ X 1,t 2 r 0, x R define τ X 1,t 1 G r (x) := P{min N (r) x}, n ] G (f ) r (x) := E [e f (r) α j #[N (r) (x+a j )] j=1 1{min N (r) x}. Hence, when f 0 we have G (f ) r (x) = 1 G r (x)). t X 1,t N τ 1 t,x 1,t N τ 2 t,x 1,t Theorem (Aidekon, B., Brunet, Shi 11) We have I(t) = E [e σbt F (σb s, s [0, t]) e 2λ t where B is BM. 0 [1 G(f t s(σb ) t σb s)]ṣ ], J Berestycki (Paris) ANR MANEGE 09/02/ / 35

52 Theorem (Aidekon, B., Brunet, Shi 11) In particular, the path (s X 1,t (s), 0 s t) is a standard Brownian motion in a potential : [ ] E F (X 1,t (s), s [0, t]) = E [e σbt F (σb s, s [0, t]) e 2λ t ] G 0 t s(σb t σb s)ds. (3) J Berestycki (Paris) ANR MANEGE 09/02/ / 35

53 path of the leftmost particle Theorem (Aidekon, B., Brunet, Shi 11) ɛ > 0 x > 0 s.t. P(A t (x)) > 1 ɛ for all t large enough. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

54 path of the leftmost particle Theorem (Aidekon, B., Brunet, Shi 11) ɛ > 0 x > 0 s.t. P(A t (x)) > 1 ɛ for all t large enough. Again see also Arguin Bovier Kistler. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

55 Theorem (Aidekon, B., Brunet, Shi 11) As t the point process (Y i (t), 1 i N(t)) converges in distribution to the point process L obtained as follows (i) Define P a Poisson point process on R +, with intensity measure 2 σ e 2x/σ dx where a is a constant. (ii) For each atom x of P {0}, we attach a point process x + Q (x) where Q (x) are i.i.d. copies of a certain point process Q. (iii) L is then the superposition of all the point processes x + Q (x) : L := {x + y : x P {0}, y Q (x) } With extreme value theory. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

56 Idea : Fix K > 0 large and stop particles when they hit k for the first time. (no escape). H k = # particles stoped at k k Z = lim k 2 1/2 ke k H k, (4) exists a.s., is in (0, ). ξ k,u u J Berestycki (Paris) ANR MANEGE 09/02/ / 35

57 Idea : Fix K > 0 large and stop particles when they hit k for the first time. (no escape). H k = # particles stoped at k k Z = lim k 2 1/2 ke k H k, (4) ξ k,u u exists a.s., is in (0, ). Start iid BBM so that u H k, X u 1 (t) = d k + X 1 (t ξ k,u ), where ξ k,u = time when u reaches k. By Bramson k 1, u H k, X u 1 (t) m t law k + Wu, t. (see also Bovier Arguin and Kistler (2010) : extremal particles in BBM either branch at the very beginning or at the end) J Berestycki (Paris) ANR MANEGE 09/02/ / 35

58 Define Pk, := δ k+w (u)+log Z. u H k Recall that P = PPP with intensity ae x dx. Proposition Pk, P In the sense of convergence in distribution. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

59 Take (X i, i N) a sequence of i.i.d. r.v. such that P(X i x) Cxe x, as x. Call M n = max i=1,...,n X i the record. Define b n = log n + log log n. Then P(M n b n y) = (P(X i y + b n )) n = (1 (1 + o(1))c(y + b n )e (y+bn) ) n ( ) 1 exp nc(y + b n ) n log n e y exp( Ce y ) J Berestycki (Paris) ANR MANEGE 09/02/ / 35

60 Take (X i, i N) a sequence of i.i.d. r.v. such that P(X i x) Cxe x, as x. Call M n = max i=1,...,n X i the record. Define b n = log n + log log n. Then P(M n b n y) = (P(X i y + b n )) n = (1 (1 + o(1))c(y + b n )e (y+bn) ) n ( ) 1 exp nc(y + b n ) n log n e y exp( Ce y ) Therefore P(M n (b n + log C) y) exp( e y ). By classical results the point process n ζ n := δ Xi b n log C i=1 converges in distribution to a Poisson point process on R with intensity e x dx. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

61 Recall that P(W x) Cxe x so apply to u H k δ W (u)+(log Hk +log log H k )+log C which converges in dist. to a PPP on R with intensity e x dx. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

62 Recall that P(W x) Cxe x so apply to u H k δ W (u)+(log Hk +log log H k )+log C which converges in dist. to a PPP on R with intensity e x dx. Use k2 1/2 e k H k Z to obtain H k σk 1 e k Z and therefore and hence log H k + log log H k k + log Z + c P k, = u H k δ k+w (u)+log Z also converges (as k ) towards a PPP on R with intensity 2 σ e 2x/σ dx. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

63 many-to-one W t := u N (t) e X(u), t 0, is the additive martingale. Because critical not UI and 0. Q F t = W t P F t. Under Q spine = BM drift (0), branch at rate 2 into two particles. Q{Ξ t = u F t } = e X(u) W t. For each u N (t), G u = a r.v. in F t. E P [ u N (t) G u ] = E Q 1 W t u N (t) [ ] = E Q e X(Ξt) G Ξ. G u Suppose that we want to check if a path X u (s) with some property A P( u = t : (X u (s), s [0, t]) A) = P(e σbt ; (σb s, s [0, t]) A). J Berestycki (Paris) ANR MANEGE 09/02/ / 35

64 path of the leftmost particle Theorem (Aidekon, B., Brunet, Shi 11) ɛ > 0 x > 0 s.t. P(A t (x)) > 1 ɛ for all t large enough. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

65 A t (x) := E 1 (x) E 2 (x) E 3 (x) E 4 (x) where the events E i are defined by { X1,t E 1 (x) := 3 } 2 log t x, { E 2 (x) := E 3 (x) := E 4 (x) := min X 1,t(s) x, min X 1,t(s) 3 [0,t] [t/2,t] 2 log t x { s [x, t/2], X 1,t (s) s 1/3} { } s [t/2, t x], X 1,t (s) X 1,t [(t s) 1/3, (t s) 2/3 ]. } J Berestycki (Paris) ANR MANEGE 09/02/ / 35

66 the event E 3 Suppose we have E 1(z) and E 2(z) for z large enough. By the many-to-one lemma, we get P(E 3(x) c, E 1(z), E 2(z)) { e z t 3/2 P s [x, t/2] : B s s 1/3, min B s z, [0,t/2] min B s 3 t/2,t 2 log t z, Bt 3 }. 2 log t + z Applying the Markov property at time t/2, it yields that { P s [x, t/2] : B s s 1/3, min B s z, min B s 3 [0,t/2] t/2,t 2 log t z, Bt 3 } 2 log t + z = E [ 1 { s [x,t/2]:bs s 1/3 } 1 {min [0,t/2] B s z} { P Bt/2 min B s 3 2 log t z, B t/2 3 } ] 2 log t + z s [0,t/2] c2zt 3/2 E [1 { s [x,t/2]:bs s 1/3} 1 {min[0,t/2] Bs z}(b t/2 3 ] 2 log t + z) c2zt 3/2 E [ 1 { s [x,t/2]:bs s 1/3 } 1 {min [0,t/2] B s z}(b t/2 + z) ] where { the second inequality comes from } bound on P min s [0,t] B s x, B t x + y. We recognize the h-transform of the Bessel. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

67 We end up with P(E 3 (x) c, E 1 (z), E 2 (z)) e z c2zp z ( s [x, t/2] : R s s 1/3 ) e z c2zp z ( s x : R s s 1/3 ) which is less than ε for x large enough. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

68 proof Step 1 : The process V x (t) := u N(t) w( 2t X i (t) + x) is a F t -martingale J Berestycki (Paris) ANR MANEGE 09/02/ / 35

69 proof Step 1 : The process V x (t) := u N(t) w( 2t X i (t) + x) is a F t -martingale Step 2 : u N(t) e 2X u(t) 2t is a positive martingale, converges to a finite value, so min u ( 2t X u (t)) = + a.s. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

70 proof Step 1 : The process V x (t) := u N(t) w( 2t X i (t) + x) is a F t -martingale Step 2 : u N(t) e 2X u(t) 2t is a positive martingale, converges to a finite value, so min u ( 2t Xu x (t)) = + a.s. Step 3 : 1 w(y) = Cye 2y. log V x (t) = u N(t) u N(t) log w( 2t X i (t) + x) C( 2t X i (t) + x)e 2t+ 2X i (t) 2x CZ(t)e 2x CY (t)xe 2x with Y (t) = u N(t) e 2X i (t) 2t. Clearly lim Y /Z = 0. lim Y (t) = Y 0 exists a.s. so Z(t) a.s. on the event Y > 0, this V x = 0. But since E(V x ) = w(x) 1 when x, P(Y > 0) = 0. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

71 proof 2 Step 4 : Thus lim Z(t) = e 2x C 1 log V x. We conclude that lim Z(t) exists and > 0 a.s. [ { }] w(x) = E[V x ( )] = E exp CZe 2x. J Berestycki (Paris) ANR MANEGE 09/02/ / 35

72 proof 2 Step 4 : Thus lim Z(t) = e 2x C 1 log V x. We conclude that lim Z(t) exists and > 0 a.s. [ { }] w(x) = E[V x ( )] = E exp CZe 2x. Step 5 : P(M(t +s) m(t +s)+x F s ) = u N(t) u(t, x +m(t +s) X u(s)). J Berestycki (Paris) ANR MANEGE 09/02/ / 35

73 proof 2 Step 4 : Thus lim Z(t) = e 2x C 1 log V x. We conclude that lim Z(t) exists and > 0 a.s. [ { }] w(x) = E[V x ( )] = E exp CZe 2x. Step 5 : P(M(t +s) m(t +s)+x F s ) = u N(t) u(t, x +m(t +s) X u(s)).recall that u(t, x + m(t)) = P(M(t) m(t) + x) w(x) and that lim t (m(t + s) m(t) 2s) = 0 J Berestycki (Paris) ANR MANEGE 09/02/ / 35

74 proof 2 Step 4 : Thus lim Z(t) = e 2x C 1 log V x. We conclude that lim Z(t) exists and > 0 a.s. [ { }] w(x) = E[V x ( )] = E exp CZe 2x. Step 5 : P(M(t +s) m(t +s)+x F s ) = u N(t) u(t, x +m(t +s) X u(s)).recall that u(t, x + m(t)) = P(M(t) m(t) + x) w(x) and that lim t (m(t + s) m(t) 2s) = 0 so that lim t P(M(t + s) m(t + s) + x F s ) = u N(t) = u N(t) := V x (s) w(x + m(t + s) X u (s) m(t) w(x + 2s X u (s)) J Berestycki (Paris) ANR MANEGE 09/02/ / 35

A. Bovier () Branching Brownian motion: extremal process and ergodic theorems

A. Bovier () Branching Brownian motion: extremal process and ergodic theorems Branching Brownian motion: extremal process and ergodic theorems Anton Bovier with Louis-Pierre Arguin and Nicola Kistler RCS&SM, Venezia, 06.05.2013 Plan 1 BBM 2 Maximum of BBM 3 The Lalley-Sellke conjecture

More information

Critical branching Brownian motion with absorption. by Jason Schweinsberg University of California at San Diego

Critical branching Brownian motion with absorption. by Jason Schweinsberg University of California at San Diego Critical branching Brownian motion with absorption by Jason Schweinsberg University of California at San Diego (joint work with Julien Berestycki and Nathanaël Berestycki) Outline 1. Definition and motivation

More information

Extremal Particles in Branching Brownian Motion - Part 1

Extremal Particles in Branching Brownian Motion - Part 1 Extremal Particles in Branching Brownian Motion - Part 1 Sandra Kliem (Part 1) und Kumarjit Saha (Part 2) Singapore - Genealogies of Interacting Particle Systems - 2017 Overview: Part 1 Part 1: Extremal

More information

Recent results on branching Brownian motion on the positive real axis. Pascal Maillard (Université Paris-Sud (soon Paris-Saclay))

Recent results on branching Brownian motion on the positive real axis. Pascal Maillard (Université Paris-Sud (soon Paris-Saclay)) Recent results on branching Brownian motion on the positive real axis Pascal Maillard (Université Paris-Sud (soon Paris-Saclay)) CMAP, Ecole Polytechnique, May 18 2017 Pascal Maillard Branching Brownian

More information

The genealogy of branching Brownian motion with absorption. by Jason Schweinsberg University of California at San Diego

The genealogy of branching Brownian motion with absorption. by Jason Schweinsberg University of California at San Diego The genealogy of branching Brownian motion with absorption by Jason Schweinsberg University of California at San Diego (joint work with Julien Berestycki and Nathanaël Berestycki) Outline 1. A population

More information

The extremal process of branching Brownian motion

The extremal process of branching Brownian motion Probab. Theory Relat. Fields DOI 1.17/s44-1-464-x The extremal process of branching Brownian motion Louis-Pierre Arguin Anton Bovier Nicola Kistler Received: 11 June 1 / Revised: 1 December 1 Springer-Verlag

More information

Yaglom-type limit theorems for branching Brownian motion with absorption. by Jason Schweinsberg University of California San Diego

Yaglom-type limit theorems for branching Brownian motion with absorption. by Jason Schweinsberg University of California San Diego Yaglom-type limit theorems for branching Brownian motion with absorption by Jason Schweinsberg University of California San Diego (with Julien Berestycki, Nathanaël Berestycki, Pascal Maillard) Outline

More information

arxiv: v4 [math.pr] 24 Sep 2012

arxiv: v4 [math.pr] 24 Sep 2012 The Annals of Applied Probability 2012, Vol. 22, No. 4, 1693 1711 DOI: 10.1214/11-AAP809 c Institute of Mathematical Statistics, 2012 arxiv:1010.2376v4 [math.pr] 24 Sep 2012 POISSONIAN STATISTICS IN THE

More information

Slowdown in branching Brownian motion with inhomogeneous variance

Slowdown in branching Brownian motion with inhomogeneous variance Slowdown in branching Brownian motion with inhomogeneous variance Pascal Maillard Ofer Zeitouni February 2, 24 Abstract We consider the distribution of the maximum M T of branching Brownian motion with

More information

Lecture 17 Brownian motion as a Markov process

Lecture 17 Brownian motion as a Markov process Lecture 17: Brownian motion as a Markov process 1 of 14 Course: Theory of Probability II Term: Spring 2015 Instructor: Gordan Zitkovic Lecture 17 Brownian motion as a Markov process Brownian motion is

More information

Mouvement brownien branchant avec sélection

Mouvement brownien branchant avec sélection Mouvement brownien branchant avec sélection Soutenance de thèse de Pascal MAILLARD effectuée sous la direction de Zhan SHI Jury Brigitte CHAUVIN, Francis COMETS, Bernard DERRIDA, Yueyun HU, Andreas KYPRIANOU,

More information

Variable Speed Branching Brownian Motion 1. Extremal Processes in the Weak Correlation Regime

Variable Speed Branching Brownian Motion 1. Extremal Processes in the Weak Correlation Regime ALEA, Lat. Am. J. Probab. Math. Stat. 1 1, 61 91 15 Variable Speed Branching Brownian Motion 1. Extremal Processes in the Weak Correlation Regime Anton Bovier and Lisa Hartung Institut für Angewandte Mathematik,

More information

Introduction to Random Diffusions

Introduction to Random Diffusions Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales

More information

1. Stochastic Processes and filtrations

1. Stochastic Processes and filtrations 1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S

More information

On a class of stochastic differential equations in a financial network model

On a class of stochastic differential equations in a financial network model 1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University

More information

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3 Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................

More information

arxiv: v3 [math.pr] 21 Sep 2016

arxiv: v3 [math.pr] 21 Sep 2016 Submitted to the Annals of Applied Probability arxiv: arxiv:1412.5975 EXTENDED CONVERGENCE OF THE EXTREMAL PROCESS OF BRANCHING BROWNIAN MOTION arxiv:1412.5975v3 [math.pr] 21 Sep 2016 By Anton Bovier,

More information

Introduction to self-similar growth-fragmentations

Introduction to self-similar growth-fragmentations Introduction to self-similar growth-fragmentations Quan Shi CIMAT, 11-15 December, 2017 Quan Shi Growth-Fragmentations CIMAT, 11-15 December, 2017 1 / 34 Literature Jean Bertoin, Compensated fragmentation

More information

Hydrodynamics of the N-BBM process

Hydrodynamics of the N-BBM process Hydrodynamics of the N-BBM process Anna De Masi, Pablo A. Ferrari, Errico Presutti, Nahuel Soprano-Loto Illustration by Eric Brunet Institut Henri Poincaré, June 2017 1 Brunet and Derrida N branching particles

More information

Convergence to a single wave in the Fisher-KPP equation

Convergence to a single wave in the Fisher-KPP equation Convergence to a single wave in the Fisher-KPP equation James Nolen Jean-Michel Roquejoffre Lenya Ryzhik Dedicated to H. Brezis, with admiration and respect Abstract We study the large time asymptotics

More information

Almost sure asymptotics for the random binary search tree

Almost sure asymptotics for the random binary search tree AofA 10 DMTCS proc. AM, 2010, 565 576 Almost sure asymptotics for the rom binary search tree Matthew I. Roberts Laboratoire de Probabilités et Modèles Aléatoires, Université Paris VI Case courrier 188,

More information

Branching Processes II: Convergence of critical branching to Feller s CSB

Branching Processes II: Convergence of critical branching to Feller s CSB Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied

More information

LAN property for ergodic jump-diffusion processes with discrete observations

LAN property for ergodic jump-diffusion processes with discrete observations LAN property for ergodic jump-diffusion processes with discrete observations Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Arturo Kohatsu-Higa (Ritsumeikan University, Japan) &

More information

Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )

Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( ) Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical

More information

I forgot to mention last time: in the Ito formula for two standard processes, putting

I forgot to mention last time: in the Ito formula for two standard processes, putting I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy

More information

Crump Mode Jagers processes with neutral Poissonian mutations

Crump Mode Jagers processes with neutral Poissonian mutations Crump Mode Jagers processes with neutral Poissonian mutations Nicolas Champagnat 1 Amaury Lambert 2 1 INRIA Nancy, équipe TOSCA 2 UPMC Univ Paris 06, Laboratoire de Probabilités et Modèles Aléatoires Paris,

More information

Real roots of random polynomials and zero crossing properties of diffusion equation

Real roots of random polynomials and zero crossing properties of diffusion equation Real roots of random polynomials and zero crossing properties of diffusion equation Grégory Schehr Laboratoire de Physique Théorique et Modèles Statistiques Orsay, Université Paris XI G. S., S. N. Majumdar,

More information

Dynamics of the evolving Bolthausen-Sznitman coalescent. by Jason Schweinsberg University of California at San Diego.

Dynamics of the evolving Bolthausen-Sznitman coalescent. by Jason Schweinsberg University of California at San Diego. Dynamics of the evolving Bolthausen-Sznitman coalescent by Jason Schweinsberg University of California at San Diego Outline of Talk 1. The Moran model and Kingman s coalescent 2. The evolving Kingman s

More information

Solution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have

Solution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have 362 Problem Hints and Solutions sup g n (ω, t) g(ω, t) sup g(ω, s) g(ω, t) µ n (ω). t T s,t: s t 1/n By the uniform continuity of t g(ω, t) on [, T], one has for each ω that µ n (ω) as n. Two applications

More information

Probability Background

Probability Background Probability Background Namrata Vaswani, Iowa State University August 24, 2015 Probability recap 1: EE 322 notes Quick test of concepts: Given random variables X 1, X 2,... X n. Compute the PDF of the second

More information

9.2 Branching random walk and branching Brownian motions

9.2 Branching random walk and branching Brownian motions 168 CHAPTER 9. SPATIALLY STRUCTURED MODELS 9.2 Branching random walk and branching Brownian motions Branching random walks and branching diffusions have a long history. A general theory of branching Markov

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

LogFeller et Ray Knight

LogFeller et Ray Knight LogFeller et Ray Knight Etienne Pardoux joint work with V. Le and A. Wakolbinger Etienne Pardoux (Marseille) MANEGE, 18/1/1 1 / 16 Feller s branching diffusion with logistic growth We consider the diffusion

More information

Stability of Stochastic Differential Equations

Stability of Stochastic Differential Equations Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010

More information

Uniformly Uniformly-ergodic Markov chains and BSDEs

Uniformly Uniformly-ergodic Markov chains and BSDEs Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,

More information

Exercises. T 2T. e ita φ(t)dt.

Exercises. T 2T. e ita φ(t)dt. Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.

More information

p 1 ( Y p dp) 1/p ( X p dp) 1 1 p

p 1 ( Y p dp) 1/p ( X p dp) 1 1 p Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p

More information

1 IEOR 6712: Notes on Brownian Motion I

1 IEOR 6712: Notes on Brownian Motion I Copyright c 005 by Karl Sigman IEOR 67: Notes on Brownian Motion I We present an introduction to Brownian motion, an important continuous-time stochastic process that serves as a continuous-time analog

More information

Brownian Motion and Stochastic Calculus

Brownian Motion and Stochastic Calculus ETHZ, Spring 17 D-MATH Prof Dr Martin Larsson Coordinator A Sepúlveda Brownian Motion and Stochastic Calculus Exercise sheet 6 Please hand in your solutions during exercise class or in your assistant s

More information

Selected Exercises on Expectations and Some Probability Inequalities

Selected Exercises on Expectations and Some Probability Inequalities Selected Exercises on Expectations and Some Probability Inequalities # If E(X 2 ) = and E X a > 0, then P( X λa) ( λ) 2 a 2 for 0 < λ

More information

Stochastic Process (ENPC) Monday, 22nd of January 2018 (2h30)

Stochastic Process (ENPC) Monday, 22nd of January 2018 (2h30) Stochastic Process (NPC) Monday, 22nd of January 208 (2h30) Vocabulary (english/français) : distribution distribution, loi ; positive strictement positif ; 0,) 0,. We write N Z,+ and N N {0}. We use the

More information

Theoretical Tutorial Session 2

Theoretical Tutorial Session 2 1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 7 9/25/2013

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 7 9/25/2013 MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 7 9/5/013 The Reflection Principle. The Distribution of the Maximum. Brownian motion with drift Content. 1. Quick intro to stopping times.

More information

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant

More information

Strong approximation for additive functionals of geometrically ergodic Markov chains

Strong approximation for additive functionals of geometrically ergodic Markov chains Strong approximation for additive functionals of geometrically ergodic Markov chains Florence Merlevède Joint work with E. Rio Université Paris-Est-Marne-La-Vallée (UPEM) Cincinnati Symposium on Probability

More information

AN EXTREME-VALUE ANALYSIS OF THE LIL FOR BROWNIAN MOTION 1. INTRODUCTION

AN EXTREME-VALUE ANALYSIS OF THE LIL FOR BROWNIAN MOTION 1. INTRODUCTION AN EXTREME-VALUE ANALYSIS OF THE LIL FOR BROWNIAN MOTION DAVAR KHOSHNEVISAN, DAVID A. LEVIN, AND ZHAN SHI ABSTRACT. We present an extreme-value analysis of the classical law of the iterated logarithm LIL

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

Partial Differential Equations, Winter 2015

Partial Differential Equations, Winter 2015 Partial Differential Equations, Winter 215 Homework #2 Due: Thursday, February 12th, 215 1. (Chapter 2.1) Solve u xx + u xt 2u tt =, u(x, ) = φ(x), u t (x, ) = ψ(x). Hint: Factor the operator as we did

More information

Stat 150 Practice Final Spring 2015

Stat 150 Practice Final Spring 2015 Stat 50 Practice Final Spring 205 Instructor: Allan Sly Name: SID: There are 8 questions. Attempt all questions and show your working - solutions without explanation will not receive full credit. Answer

More information

On Ergodic Impulse Control with Constraint

On Ergodic Impulse Control with Constraint On Ergodic Impulse Control with Constraint Maurice Robin Based on joint papers with J.L. Menaldi University Paris-Sanclay 9119 Saint-Aubin, France (e-mail: maurice.robin@polytechnique.edu) IMA, Minneapolis,

More information

MA8109 Stochastic Processes in Systems Theory Autumn 2013

MA8109 Stochastic Processes in Systems Theory Autumn 2013 Norwegian University of Science and Technology Department of Mathematical Sciences MA819 Stochastic Processes in Systems Theory Autumn 213 1 MA819 Exam 23, problem 3b This is a linear equation of the form

More information

Hybrid Atlas Model. of financial equity market. Tomoyuki Ichiba 1 Ioannis Karatzas 2,3 Adrian Banner 3 Vassilios Papathanakos 3 Robert Fernholz 3

Hybrid Atlas Model. of financial equity market. Tomoyuki Ichiba 1 Ioannis Karatzas 2,3 Adrian Banner 3 Vassilios Papathanakos 3 Robert Fernholz 3 Hybrid Atlas Model of financial equity market Tomoyuki Ichiba 1 Ioannis Karatzas 2,3 Adrian Banner 3 Vassilios Papathanakos 3 Robert Fernholz 3 1 University of California, Santa Barbara 2 Columbia University,

More information

Supplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION. September 2017

Supplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION. September 2017 Supplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION By Degui Li, Peter C. B. Phillips, and Jiti Gao September 017 COWLES FOUNDATION DISCUSSION PAPER NO.

More information

The Central Limit Theorem: More of the Story

The Central Limit Theorem: More of the Story The Central Limit Theorem: More of the Story Steven Janke November 2015 Steven Janke (Seminar) The Central Limit Theorem:More of the Story November 2015 1 / 33 Central Limit Theorem Theorem (Central Limit

More information

Lecture 21 Representations of Martingales

Lecture 21 Representations of Martingales Lecture 21: Representations of Martingales 1 of 11 Course: Theory of Probability II Term: Spring 215 Instructor: Gordan Zitkovic Lecture 21 Representations of Martingales Right-continuous inverses Let

More information

PDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation

PDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation PDEs, Homework #3 Solutions. Use Hölder s inequality to show that the solution of the heat equation u t = ku xx, u(x, = φ(x (HE goes to zero as t, if φ is continuous and bounded with φ L p for some p.

More information

Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet 1

Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet 1 Stochastic Calculus and Black-Scholes Theory MTH772P Exercises Sheet. For ξ, ξ 2, i.i.d. with P(ξ i = ± = /2 define the discrete-time random walk W =, W n = ξ +... + ξ n. (i Formulate and prove the property

More information

Hypothesis testing for Stochastic PDEs. Igor Cialenco

Hypothesis testing for Stochastic PDEs. Igor Cialenco Hypothesis testing for Stochastic PDEs Igor Cialenco Department of Applied Mathematics Illinois Institute of Technology igor@math.iit.edu Joint work with Liaosha Xu Research partially funded by NSF grants

More information

Lower Tail Probabilities and Normal Comparison Inequalities. In Memory of Wenbo V. Li s Contributions

Lower Tail Probabilities and Normal Comparison Inequalities. In Memory of Wenbo V. Li s Contributions Lower Tail Probabilities and Normal Comparison Inequalities In Memory of Wenbo V. Li s Contributions Qi-Man Shao The Chinese University of Hong Kong Lower Tail Probabilities and Normal Comparison Inequalities

More information

MS&E 321 Spring Stochastic Systems June 1, 2013 Prof. Peter W. Glynn Page 1 of 10

MS&E 321 Spring Stochastic Systems June 1, 2013 Prof. Peter W. Glynn Page 1 of 10 MS&E 321 Spring 12-13 Stochastic Systems June 1, 2013 Prof. Peter W. Glynn Page 1 of 10 Section 3: Regenerative Processes Contents 3.1 Regeneration: The Basic Idea............................... 1 3.2

More information

On the submartingale / supermartingale property of diffusions in natural scale

On the submartingale / supermartingale property of diffusions in natural scale On the submartingale / supermartingale property of diffusions in natural scale Alexander Gushchin Mikhail Urusov Mihail Zervos November 13, 214 Abstract Kotani 5 has characterised the martingale property

More information

Continuous-state branching processes, extremal processes and super-individuals

Continuous-state branching processes, extremal processes and super-individuals Continuous-state branching processes, extremal processes and super-individuals Clément Foucart Université Paris 13 with Chunhua Ma Nankai University Workshop Berlin-Paris Berlin 02/11/2016 Introduction

More information

(A n + B n + 1) A n + B n

(A n + B n + 1) A n + B n 344 Problem Hints and Solutions Solution for Problem 2.10. To calculate E(M n+1 F n ), first note that M n+1 is equal to (A n +1)/(A n +B n +1) with probability M n = A n /(A n +B n ) and M n+1 equals

More information

Preliminary Exam: Probability 9:00am 2:00pm, Friday, January 6, 2012

Preliminary Exam: Probability 9:00am 2:00pm, Friday, January 6, 2012 Preliminary Exam: Probability 9:00am 2:00pm, Friday, January 6, 202 The exam lasts from 9:00am until 2:00pm, with a walking break every hour. Your goal on this exam should be to demonstrate mastery of

More information

A Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1

A Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1 Chapter 3 A Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1 Abstract We establish a change of variable

More information

Brownian Motion. Chapter Definition of Brownian motion

Brownian Motion. Chapter Definition of Brownian motion Chapter 5 Brownian Motion Brownian motion originated as a model proposed by Robert Brown in 1828 for the phenomenon of continual swarming motion of pollen grains suspended in water. In 1900, Bachelier

More information

Research Article Existence and Uniqueness Theorem for Stochastic Differential Equations with Self-Exciting Switching

Research Article Existence and Uniqueness Theorem for Stochastic Differential Equations with Self-Exciting Switching Discrete Dynamics in Nature and Society Volume 211, Article ID 549651, 12 pages doi:1.1155/211/549651 Research Article Existence and Uniqueness Theorem for Stochastic Differential Equations with Self-Exciting

More information

ERRATA: Probabilistic Techniques in Analysis

ERRATA: Probabilistic Techniques in Analysis ERRATA: Probabilistic Techniques in Analysis ERRATA 1 Updated April 25, 26 Page 3, line 13. A 1,..., A n are independent if P(A i1 A ij ) = P(A 1 ) P(A ij ) for every subset {i 1,..., i j } of {1,...,

More information

Universal examples. Chapter The Bernoulli process

Universal examples. Chapter The Bernoulli process Chapter 1 Universal examples 1.1 The Bernoulli process First description: Bernoulli random variables Y i for i = 1, 2, 3,... independent with P [Y i = 1] = p and P [Y i = ] = 1 p. Second description: Binomial

More information

x ct x + t , and the characteristics for the associated transport equation would be given by the solution of the ode dx dt = 1 4. ξ = x + t 4.

x ct x + t , and the characteristics for the associated transport equation would be given by the solution of the ode dx dt = 1 4. ξ = x + t 4. . The solution is ( 2 e x+ct + e x ct) + 2c x+ct x ct sin(s)dx ( e x+ct + e x ct) + ( cos(x + ct) + cos(x ct)) 2 2c 2. To solve the PDE u xx 3u xt 4u tt =, you can first fact the differential operat to

More information

Erdős-Renyi random graphs basics

Erdős-Renyi random graphs basics Erdős-Renyi random graphs basics Nathanaël Berestycki U.B.C. - class on percolation We take n vertices and a number p = p(n) with < p < 1. Let G(n, p(n)) be the graph such that there is an edge between

More information

Extreme values of two-dimensional discrete Gaussian Free Field

Extreme values of two-dimensional discrete Gaussian Free Field Extreme values of two-dimensional discrete Gaussian Free Field Marek Biskup (UCLA) Joint work with Oren Louidor (Technion, Haifa) Midwest Probability Colloquium, Evanston, October 9-10, 2015 My collaborator

More information

The Chaotic Character of the Stochastic Heat Equation

The Chaotic Character of the Stochastic Heat Equation The Chaotic Character of the Stochastic Heat Equation March 11, 2011 Intermittency The Stochastic Heat Equation Blowup of the solution Intermittency-Example ξ j, j = 1, 2,, 10 i.i.d. random variables Taking

More information

Travelling-waves for the FKPP equation via probabilistic arguments

Travelling-waves for the FKPP equation via probabilistic arguments Travelling-waves for the FKPP equation via probabilistic arguments Simon C. Harris Department of Mathematical Sciences, University of Bath, Bath, BA2 7AY, United Kingdom. Abstract. We outline a completely

More information

Solutions to the Exercises in Stochastic Analysis

Solutions to the Exercises in Stochastic Analysis Solutions to the Exercises in Stochastic Analysis Lecturer: Xue-Mei Li 1 Problem Sheet 1 In these solution I avoid using conditional expectations. But do try to give alternative proofs once we learnt conditional

More information

Question 1. Question Solution to final homework!! Friday 5 May Please report any typos immediately for a timely response.

Question 1. Question Solution to final homework!! Friday 5 May Please report any typos immediately for a timely response. 1 Please report any typos immediately for a timely response. Call EW tn W tn 1 near u n 1,...,..., W t1 near u 1 ) LHS for short. Add the quantity W tn 1 W tn 1 into the expected value of the LHS. Apply

More information

Strong Markov property of determinatal processes

Strong Markov property of determinatal processes Strong Markov property of determinatal processes Hideki Tanemura Chiba university (Chiba, Japan) (August 2, 2013) Hideki Tanemura (Chiba univ.) () Markov process (August 2, 2013) 1 / 27 Introduction The

More information

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539 Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory

More information

Kolmogorov Equations and Markov Processes

Kolmogorov Equations and Markov Processes Kolmogorov Equations and Markov Processes May 3, 013 1 Transition measures and functions Consider a stochastic process {X(t)} t 0 whose state space is a product of intervals contained in R n. We define

More information

Exponential functionals of Lévy processes

Exponential functionals of Lévy processes Exponential functionals of Lévy processes Víctor Rivero Centro de Investigación en Matemáticas, México. 1/ 28 Outline of the talk Introduction Exponential functionals of spectrally positive Lévy processes

More information

Two viewpoints on measure valued processes

Two viewpoints on measure valued processes Two viewpoints on measure valued processes Olivier Hénard Université Paris-Est, Cermics Contents 1 The classical framework : from no particle to one particle 2 The lookdown framework : many particles.

More information

Chapter 1. Poisson processes. 1.1 Definitions

Chapter 1. Poisson processes. 1.1 Definitions Chapter 1 Poisson processes 1.1 Definitions Let (, F, P) be a probability space. A filtration is a collection of -fields F t contained in F such that F s F t whenever s

More information

Large Deviations Principles for McKean-Vlasov SDEs, Skeletons, Supports and the law of iterated logarithm

Large Deviations Principles for McKean-Vlasov SDEs, Skeletons, Supports and the law of iterated logarithm Large Deviations Principles for McKean-Vlasov SDEs, Skeletons, Supports and the law of iterated logarithm Gonçalo dos Reis University of Edinburgh (UK) & CMA/FCT/UNL (PT) jointly with: W. Salkeld, U. of

More information

1.1 Definition of BM and its finite-dimensional distributions

1.1 Definition of BM and its finite-dimensional distributions 1 Brownian motion Brownian motion as a physical phenomenon was discovered by botanist Robert Brown as he observed a chaotic motion of particles suspended in water. The rigorous mathematical model of BM

More information

Weak solutions of mean-field stochastic differential equations

Weak solutions of mean-field stochastic differential equations Weak solutions of mean-field stochastic differential equations Juan Li School of Mathematics and Statistics, Shandong University (Weihai), Weihai 26429, China. Email: juanli@sdu.edu.cn Based on joint works

More information

P (A G) dp G P (A G)

P (A G) dp G P (A G) First homework assignment. Due at 12:15 on 22 September 2016. Homework 1. We roll two dices. X is the result of one of them and Z the sum of the results. Find E [X Z. Homework 2. Let X be a r.v.. Assume

More information

{σ x >t}p x. (σ x >t)=e at.

{σ x >t}p x. (σ x >t)=e at. 3.11. EXERCISES 121 3.11 Exercises Exercise 3.1 Consider the Ornstein Uhlenbeck process in example 3.1.7(B). Show that the defined process is a Markov process which converges in distribution to an N(0,σ

More information

Other properties of M M 1

Other properties of M M 1 Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected

More information

Exact Simulation of Diffusions and Jump Diffusions

Exact Simulation of Diffusions and Jump Diffusions Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction

More information

A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin

A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin Mario Bukal A. Jüngel and D. Matthes ACROSS - Centre for Advanced Cooperative Systems Faculty of Electrical Engineering and Computing

More information

Mathematical Methods - Lecture 9

Mathematical Methods - Lecture 9 Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr

More information

On continuous time contract theory

On continuous time contract theory Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem

More information

Reflected Brownian Motion

Reflected Brownian Motion Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide

More information

Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao

Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley

More information

Scaling limits for random trees and graphs

Scaling limits for random trees and graphs YEP VII Probability, random trees and algorithms 8th-12th March 2010 Scaling limits for random trees and graphs Christina Goldschmidt INTRODUCTION A taste of what s to come We start with perhaps the simplest

More information

E. Vitali Lecture 6. Brownian motion

E. Vitali Lecture 6. Brownian motion Markov chains, presented in the previous chapter, belong to a wide class of mathematical objects, the stochastic processes, which will introduced in the present chapter. Stochastic processes form the basis

More information

On Stopping Times and Impulse Control with Constraint

On Stopping Times and Impulse Control with Constraint On Stopping Times and Impulse Control with Constraint Jose Luis Menaldi Based on joint papers with M. Robin (216, 217) Department of Mathematics Wayne State University Detroit, Michigan 4822, USA (e-mail:

More information

Quenched invariance principle for random walks on Poisson-Delaunay triangulations

Quenched invariance principle for random walks on Poisson-Delaunay triangulations Quenched invariance principle for random walks on Poisson-Delaunay triangulations Institut de Mathématiques de Bourgogne Journées de Probabilités Université de Rouen Jeudi 17 septembre 2015 Introduction

More information

Wiener Measure and Brownian Motion

Wiener Measure and Brownian Motion Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u

More information

Weak convergence and large deviation theory

Weak convergence and large deviation theory First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov

More information