On the smoothness of the partition function for multiple Schramm-Loewner evolutions

Size: px
Start display at page:

Download "On the smoothness of the partition function for multiple Schramm-Loewner evolutions"

Transcription

1 On the smoothness of the partition function for multiple Schramm-Loewner evolutions Mohmammad Jahangoshahi Gregory F. Lawler August Abstract We consider the measure on multiple chordal Schramm-Loewner evolution (SLE κ ). We estabalish a derivative estimate and use it to give a direct proof that the partition function is C 2 if κ < 4. 1 Introduction The (chordal) Schramm-Loewner evolution with parameter κ > 0 (SLE κ ) is a measure on curves connecting two distinct boundary points z w of a simply connected domain D. As originally defined by Schramm 6 this is a probability measure on paths. If κ 4 the measure is supported on simple curves that do not touch the boundary. Although Schramm 6 originally defined SLE κ as a probability measure if κ 4 and z w are locally analytic boundary points it is natural to consider SLE κ as a finite measure µ(z w) with partition function that is as a measure with total mass Ψ D (z w) = H D (z w) b. Here H denotes the boundary Poisson kernel normalized so that H H (0 x) = x 2 and b = (6 κ)/2κ is the boundary scaling exponent. If f : D f(d) is a conformal transformation then H D (z w) = f (z) f (w) H f(d) (f(z) f(w)). There are several reasons for considering SLE κ as a measure with a total mass. First SLE κ is known to be the scaling limit of various two-dimensional discrete models that are considered as measures with partition functions and hence it is natural to consider the (appropriately normalized) partition function in the scaling limit. Second the restriction property or boundary perturbation can be described more naturally for the nonprobability measures; see (1) below. This description leads to one way to define SLE κ in multiply connected domains or as is important for this paper for multiple SLE κ in a simply connected domain. See 3 9 for more information. We write µ # D (z w) = µ D(z w)/ψ D (z w) for the probability measure which is well defined even for rough boundaries. The definition of the measure on multiple SLE κ paths immediately gives a partition function defined as the total mass of the measure. The measure on multiple SLE κ paths γ = (γ 1... γ n ) Research supported by National Science Foundation grant DMS

2 has been constructed in Even though the definition in 2 is given for the so-called rainbow arrangement of the boundary points it can be easily extended to the other arrangements One can see that unlike SLE κ measure on single curves conformal invariance and domain Markov property do not uniquely specify the measure when n 2. This definition makes it unique by requiring the measure to satisfy the restriction property which is explained in Section 2. Study of the multiple SLE κ measure involves characterizing the partition function. For n = 2 the partition function is explicitly given in terms of the hypergeometric function. For n 3 the goal is to characterize the partition function by a particular second-order PDE. However it does not directly follow from the definition that the partition function is C 2. There are two main approaches to address this problem. One approach is to show that the PDE system has a solution and use it to describe the partition function. In 1 it is shown that a family of integrals taken on a specific set of cycles satisfy the required PDE system. In 8 conformal field theory and partial differential equation techniques such as Hörmander s theorem are used to show that the partition function satisfies the PDE system. The other approach which is the one we take in this work is to directly prove that the partition function is C 2. Then Itô s formula can be used next to show that the partition function satisfies the PDEs. The basic idea of our proof is to interchange derivatives and expectations in expressions for the partition function. This interchange needs justification and we prove an estimate about SLE κ to justfiy this. Here we summarize the paper. We finish this introduction by reviewing examples of partition functions for SLE κ. Definitions and properties of multiple SLE κ and the outline of the proof are given in Section 2. Section 3 includes an estimate for SLE κ using techniques similar to the ones in 4. Proof of Lemma 2.3 which explains estimates for derivatives of the Poisson kernel is given in Section Examples SLE κ in a subset of H. Let κ 4 and suppose D H is a simply connected domain such that K = H \ D is bounded and dist(0 K) > 0. Also assume that γ is parameterized with half-plane capacity. By the restriction property we have dµ D (0 ) { c } dµ H (0 ) (γ) = 1{γ D = } exp 2 m H(γ K) (1) where m H (γ D) denotes the Brownian loop measure of the loops that intersect both γ and K and (6 κ)(3κ 8) c = 2κ is the central charge. We normalize the partition functions so that Ψ H (0 ) = 1. For an initial segment of the curve γ t let g t : H \ γ t H be the unique conformal transformation with g t (z) = z + o(1) as z. Then t g t (z) = a g t (z) U t 2

3 where a = 2/κ and U t is a standard Brownian motion. Suppose γ t K = and let D t = g t (D \ γ t ). One can see that m H (γ t K) = a 6 t 0 SΦ s (U s )ds where S denotes the Schwarzian derivative and Φ s (U s ) = H Ds (U s ). It follows from conditioning on γ t that { c } M t = exp 2 m H(γ t K) Ψ Dt (U t ) is a martingale. We assume the function V (t x) = Ψ Dt (x ) is C 2 for a moment. Therefore we can apply Itô s formula and we get ac 12 V (t U t) SΦ t (U t ) + t V (t U t ) xxv (t U t ) = 0. Straightforward calculation shows that V (t x) = H Dt (x ) b is C 2 and satisfies this PDE. Here b is the boundary scaling exponent b = 6 κ 2κ. Other examples. Similar ideas were used in 2 to describe the partition function of two SLE κ curves with a PDE. Differentiability of the partition function was justified using the explicit form of the solution in terms of the hypergeometric function. The PDE system in 9 characterizes the partition function of the annulus SLE κ. That PDE is more complicated and one cannot find an explicit form for the solution. In fact it is not easy to even show that the PDE has a solution. Instead it was directly proved that the partition function is C 2 and Itô s formula was used to derive the PDE. 2 Definitions and Preliminaries We will consider the multiple SLE κ measure only for κ 4 on simply connected domains D and distinct locally analytic boundary points x = (x 1... x n ) y = (y 1... y n ). The measure is supported on n-tuples of curves γ = (γ 1... γ n ) where γ j is a curve connecting x j to y j in D. If n = 1 then µ D (x 1 y 1 ) is SLE κ from x 1 to y 1 in D with total mass H D (x 1 y 1 ) b whose corresponding probability measure µ # D (x 1 y 1 ) = µ D (x 1 y 1 )/H D (x 1 y 1 ) b is (a time change of) SLE κ from x 1 to y 1 as defined by Schramm. Definition If κ 4 and n 1 then µ D (x y) is the measure absolutely continuous with respect to µ D (x 1 y 1 ) µ D (x n y n ) with Radon-Nikodym derivative c n Y (γ) := I(γ) exp m K j (γ) 2. j=2 3

4 Here c = (6 κ)(3κ 8)/2κ is the central charge I(γ) is the indicator function of the event {γ j γ k = 1 j < k n} and m K j (γ) denotes the Brownian loop measure of loops that intersect at least j of the paths γ 1... γ n. Brownian loop measure is a measure on (continuous) curves η : 0 t γ C with η(0) = η(t η ). Let ν # (0 0; 1) be the law of the Brownian bridge starting from 0 and returning to 0 at time 1. Brownian loop measure can be considered as the measure ( ) 1 m = area 2πt 2 dt ν # (0 0; 1) on the triplets (z t η η) where η(t) = tη 1/2 η(t/t η ) for t 0 1. For a domain D C we denote the restriction of m to the loops η D by m D. One important property of m D is conformal invariance. More precisely if f : D f(d) is a conformal transformation then f m D = m f(d) where f m D is the pushforward measure. Note that if σ is a permutation of {1... n} and γ σ = (γ σ(1)... γ σ(n) ) then Y (γ) = Y (γ σ ). The partition function is the total mass of this measure Ψ D (x y) = µ D (x y). We also write which can also be written as Ψ D (x y) = Ψ D (x y) n j=1 H D(x j y j ) b Ψ D (x y) = EY where the expectation is with respect to the probability measure µ # D (x 1 y 1 ) µ # D (x n y n ). Note that Ψ D (x y) is a conformal invariant f Ψ D (x y) = Ψ f(d) (f(x) f(y)) and hence is well defined even if the boundaries are rough. Since SLE κ is reversible 7 interchanging x j and y j does not change the value. To compute the partition function we use an alternative description of the measure µ D (x y). We will give a recursive definition. For n = 1 µ D (x 1 y 1 ) is the usual SLE κ measure with total mass H D (x 1 y 1 ) b. Suppose the measure has been defined for all n-tuples of paths. Suppose x = (x x n+1 ) y = (y y n+1 ) are given and write an (n + 1)-tuple of paths as γ = (γ γ (n+1) ). The marginal measure on γ induced by µ D (x y) is absolutely continuous with respect to µ D (x y ) with Radon-Nikodym derivative H D(x n+1 y n+1 ) b. Here D is the component of D \ γ containing x n+1 y n+1 on its boundary. (If there is no such component then we set H D(x n+1 y n+1 ) = 0 and µ D (x y) is the zero measure.) 4

5 Given γ the curve γ n+1 is chosen using the probability distribution µ # D(z n+1 y n+1 ). One could try to use this description of the measure as the definition but it is not obvious that it is consistent. However one can see that the first definition satisfies this property using the following lemma. Lemma 2.1. Let γ denote (n + 1)-tuple of paths as γ = (γ γ (n+1) ) and let D be the connected component of D \ γ containing the end points of γ (n+1) on its boundary. Then n+1 m K j (γ) = j=2 n m K j (γ ) + m D (γ (n+1) D \ D). j=2 Proof. Let Kj 1(γ) denote the set of loops in K j(γ) that intersect γ (n+1) and let Kj 2 (γ) denote the set of loops that do not intersect γ (n+1). Then m K 1 2(γ) = m D (γ (n+1) D \ D). (2) Note that K 1 j (γ) is equivalent to the set of loops in D that intersect γ(n+1) and at least j 1 paths from γ. Moreover K 2 j (γ) is equivalent to the set of loops that intersect least j paths from γ but do not intersect γ (n+1). Therefore K j (γ ) = K 1 j+1(γ) K 2 j (γ). Now the result follows from this the fact that Kn+1 2 (γ) = and (2). We can also take the marginals in a different order. For example we could have defined the recursive step above as follows. The marginal measure on γ n+1 induced by µ D (x y) is absolutely continuous with respect to µ D (x n+1 y n+1 ) with Radon-Nikodym derivative Ψ D(x y ) where D = D \ γ. (It is possible that D has two separate components in which case we multiply the partition functions on the two components.) We will consider boundary points on the real line. We write just H Ψ Ψ µ µ # for H H Ψ H Ψ H µ H µ # H ; and note that Ψ(x y) = E Y = Ψ(x y) n y j x j 2b where the expectation is with respect to the probability measure If n = 1 then Y 1 and Ψ(x y) = 1. j=1 µ # (x 1 y 1 ) µ # (x n y n ). 5

6 For n = 2 and γ = (γ 1 γ 2 ) then HD\γ EY γ 1 1(x 2 y 2 ) b =. H D (x 2 y 2 ) The right-hand side is well defined even for non smooth boundaries provided that γ 1 stays a positive distance from x 2 y 2. In particular EY = E E(Y γ 1 ) (HD\γ ) 1(x 2 y 2 ) b = E 1. H D (x 2 y 2 ) If 8/3 < κ 4 then c > 0 and Y > 1 on the event I(γ) so the inequality EY 1 is not obvious. More generally if γ = (γ γ n+1 ) Using this we see that Ψ D (x y) 1. EY γ = Y (γ HD\γ (x n+1 y n+1 ) b ) Y (γ ). H D (x n+1 y n+1 ) For n = 2 if x 1 = 0 y 1 = y 2 = 1 and x 2 = x with 0 < x < 1 we have (see for example 2 (3.7)) Γ(2a) Γ(6a 1) Ψ(x y) = φ(x) := Γ(4a) Γ(4a 1) xa F (2a 1 2a 4a; x) (3) where F = 2 F 1 denotes the hypergeometric function and a = 2/κ. finding E H H\γ 1(x 1) b. In fact this calculation is valid for κ < 8 if it is interpreted as E H H\γ 1(x 1) b ; H H\γ 1(x 1) > 0. This is computed by It will be useful to write the conformal invariant (3) in a different way. If V 1 V 2 are two arcs of a domain D let E D (V 1 V 2 ) = H D (z w) dz dw. V 1 V 2 This is π times the usual excursion measure between V 1 and V 2 ; the factor of π comes from our choice of Poisson kernel. Note that 1 0 dr ds 1 E H (( 0 x 1) = x (s r) 2 = dr x r = log(1/x) Hence we can write (3) as φ (exp { E H (( 0 x 1)}). More generally if x 1 < y 1 < x 2 < y 2 ( }) Ψ(x y) = φ (exp { E H (x 1 y 1 x 2 y 2 )}) = φ exp 6 { y1 y2 x 1 x 2 dr ds (s r) 2

7 and if D is a simply connected subdomain of H containing x 1 y 2 x 2 y 2 on its boundary then ( { y1 y2 }) Ψ D (x y) = φ (exp { E D (x 1 y 1 x 2 y 2 )}) = φ exp H D (r s) dr ds. (4) This expression is a little bulky but it allows for easy differentiation with respect to x 1 x 2 y 1 y 2. At this point we can state the main proposition. Proposition 2.2. Ψ and Ψ are C 2 functions. It clearly suffices to prove this for Ψ. The idea is simple we will write the partition function as an expectation and differentiate the expectation by interchanging the expectation and the derivatives. This interchange requires justification and this is the main work of this paper. We will use the following fact which is an analogue of derivative estimates for positive harmonic functions. The proof is straightforward but we delay it to Section 4. Lemma 2.3. There exists c < such that for every x 1 < y 1 < x 2 < y 2 the following holds. Suppose D H is a simply connected domain whose boundary contains an open real neighborhood of x 1 y 1 and suppose that Then if z 1 z 2 {x 1 y 1 } δ := min { x 1 y 1 dist {x 1 y 1 } H \ D} > 0. x 1 z1 H D (x 1 y 1 ) c δ 1 H D (x 1 y 1 ). z1 z 2 H D (x 1 y 1 ) c δ 2 H D (x 1 y 1 ). Suppose D H is a simply connected domain whose boundary contains open real neighborhoods of x 1 y 1 and x 2 y 2 and suppose that δ := min {{ w 1 w 2 ; w 1 w 2 and w 1 w 2 {x 1 x 2 y 1 y 2 }} dist {x 1 y 1 x 2 y 2 } H \ D}. Then if z 1 {x 1 y 1 } z 2 {x 2 y 2 } z1 z 2 ΨD (x y) c δ 2 ΨD (x y). Moreover the constant can be chosen uniformly in neighborhoods of x 1 y 1 x 2 y 2. We will also need to show that expectations do not blow up when paths get close to starting points. We prove this lemma in Section 3. Let jk (γ) = dist { {x k y k } γ j} (γ) = min j k jk(γ). x 2 7

8 Lemma 2.4. If κ < 4 then for every n and every (x y) there exists c < such that for all ɛ > 0 and all j k E Y ; ɛ c ɛ 12 κ 1. In particular E Y 2 m= Proof. It suffices to show that for each j k 2 2m E Y ; 2 m < 2 m+1 <. E Y ; jk ɛ c ɛ 12 κ 1 and by symmetry we may assume j = 1 k = 2. If we write γ = (γ 1 γ 2 γ ) then the event { 12 ɛ} is measurable with respect to (γ 1 γ 2 ) and EY γ 1 γ 2 Y (γ 1 γ 2 ). Hence it suffices to prove the result when n = 2. This will be done in Section 3; in that section we consider κ < 8. For n = 1 2 it is clear that Ψ is C from the exact expression so we will assume that n 3. By invariance under permutation of indices it suffices to consider second order derivatives involving only x 1 x 2 y 1 y 2. We will assume x j < y j for j = 1 2 and x 1 < x 2 (otherwise we just relabel the vertices). The configuration x 1 < x 2 < y 1 < y 2 is impossible for topological reasons. If x 1 < x 2 < y 2 < y 1 we can find a Móbius transformation taking a point y (y 2 y 1 ) to and then the images would satisfy y 1 < x 1 < x 2 < y 2 and this reduces to above. So we may assume that x 1 < y 1 < x 2 < y 2. Case I: Derivatives involving only x j y j for some j. We assume j = 1. We will write x = (x x ) y = (y y ) γ = (γ 1 γ ) and let D be the connected component of H \ γ containing x y on the boundary. Then EY γ = Y (γ HD (x y) b ) = Y (γ ) Q D (x y) b H(x y) where Q D (x y) is the probability that a (Brownian) excursion in H from x to y stays in D. Hence Ψ(x y) = E Y (γ ) Q D (x y) b. Let δ = δ(γ ) = dist{{x y} γ }. Using Lemma 2.3 we see that x Q D (x y) b c δ 1 Q D (x 1 y 1 ) b xy Q D (x y) b + xx Q D (x y) b c δ 2 Q D (x 1 y 1 ) b. (Here c may depend on x y but not on D). Hence E Y (γ ) x Q D (x y) b c E Y (γ ) δ(γ ) 1 Q D (x y) b 8

9 and if z = x or y E Y (γ ) xz Q D (x y) b c E Y (γ ) δ(γ ) 2 Q D (x y) b. Since E Y (γ ) δ(γ ) 2 Q D (x y) b = E E ( Y δ 2 γ ) = EY δ 2 EY 2 < the interchange of expectation and derivative is valid x Ψ(x y) = E Y (γ ) x Q D (x y) b xz Ψ(x y) = E Y (γ ) xz Q D (x y) b. Case 2: The partial z1 z 2 where z 1 {x j y j } z 2 {x k y k } with j k. We assume j = 1 k = 2. We will write x = (x 1 x 2 x ) y = (y 1 y 2 y ) γ = (γ 1 γ 2 γ ). We will write D = D \ γ and let D 1 D 2 be the connected components of D containing {x 1 y 1 } and {x 2 y 2 } on the boundary. It is possible that D 1 = D 2 or D 1 D 2. If D 1 D 2 then If D 1 = D 2 = D then EY γ = Y (γ ) Q D1 (x 1 y 1 ) b Q D2 (x 2 y 2 ) b. EY γ = Y (γ ) Q D1 (x 1 y 1 ) b Q D2 (x 2 y 2 ) b ΨD ((x 1 x 2 ) (y 1 y 2 )) where Ψ D is defined as in (4). In either case we have written EY γ = Y (γ ) Φ(z; γ ) where z = (x 1 y 1 x 2 y 2 ) and we can use Lemma 2.3 to see that z1 z 2 Φ(z; γ ) c (γ z) 2 Φ(z γ ) (γ z) = dist{γ {x 1 y 1 x 2 y 2 }}. As in the previous case we can now interchange the derivatives and the expectation. 3 Estimate In this section we will derive an estimate for SLE κ κ < 8. While the estimate is valid for all κ < 8 the result is only strong enough to prove our main result for κ < 4. We follow the ideas in 4 where careful analysis was made of the boundary exponent for SLE. Let g t denote the usual conformal transformation associated to the SLE κ path γ from 0 to parametrized so that t g t (z) = a g t (z) U t (5) where a = 2/κ and U t = W t is a standard Brownian motion. Throughout we assume that κ < 8 so that D = D = H \ γ is a nonempty set. If 0 < x < y < we let Φ = Φ(x y) = H D(x y) H H (x y) 9

10 where H denotes the boundary Poisson kernel. If x and y are on the boundary of different components of D (which can only happen for 4 < κ < 8) then H D (x y) = 0. As usual we let b = 6 κ 2κ = 3a 1. 2 As a slight abuse of notation we will write Φ b for Φ b 1{Φ > 0} even if b 0. Proposition 3.1. For every κ < 8 and δ > 0 there exists 0 < c < such that for all δ x < y 1/δ and all 0 < ɛ < (y x)/10 E Φ b ; dist({x y} γ) < ɛ c ɛ 6a 1. It is already known that P {dist({x y} γ) < ɛ} ɛ 4a 1 and hence we can view this as the estimate E Φ b dist({x y} γ) < ɛ c ɛ 2a. Using reversibility 5 7 and scaling of SLE κ we can see that to prove the proposition it suffices to show that for every δ > 0 there exists c = c δ such that if δ x < 1 E Φ b ; dist(1 γ) < ɛ c ɛ 6a 1. This is the result we will prove. Proposition 3.2. If κ < 8 there exists c < such that if γ is an SLE κ curve from 0 to 0 < x < 1 Φ = Φ(x 1) 0 < ɛ 1/2 E Φ b ; dist(γ 1) < ɛ (1 x) c x a (1 x) 4a 1 ɛ 6a 1. We will relate the distance to the curve to a conformal radius. In order to do this we will need 1 to be an interior point of the domain. Let D t be the unbounded component of K t = H \ ( x γ t { z : z γ t } and let T = T 1 = inf{t : 1 Dt }. Then for t < T the distance from 1 to Dt is the minimum of 1 x and dist(1 γ t ). In particular if t < T and ɛ < 1 x then dist(γ t 1) ɛ if and only if dist(1 Dt ) < ɛ. We define Υ t to be 4(1 x) 1 times the conformal radius of 1 with respect to Dt and Υ = Υ. Note that Υ 0 = 1 and if dist(1 Dt ) ɛ(1 x) then Υ ɛ. It suffices for us to show that E Φ b ; Υ < ɛ c ɛ 6a 1. We set up some notation. We fix 0 < x < 1 and assume that g t satisfies (5). Let X t = g t (1) U t Z t = g t (x) U t Y t = X t Z t K t = Z t X t 10

11 and note that the scaling rule for conformal radius implies that The Loewner equation implies that Υ t = Y t (1 x) g t (1). dx t = a X t dt + db t dz t = a Z t dt + db t t g t(1) = a g t(1) X 2 t t Υ t = Υ t a X 2 t t g t(x) = ag t(x) a X t Z t Let D t be the unbounded component of H \ γ t and let Z 2 t = aυ t 1 X 2 t Φ t = H D t (x 1) H D0 (x 1) = x2 g t(x) g t(1) Yt 2 t Y t = a Y t X t Z t. 1 K t K t. where we set Φ t = 0 if x is not on the boundary of D t that is if x has been swallowed by the path (this is relevant only for 4 < κ < 8). Note that Φ = Φ and Itô s formula implies that t Φ b t = Φ b t ab Xt 2 ab Zt 2 + 2ab ( ) = ab Φb t 1 2 Kt X t Z t Xt 2 K t { t ( ) Φ b Ks t = exp ab ds}. d 1 X t = 1 X 2 t dx t + 1 X 3 t 0 X 2 s K s d X t = 1 X t 1 a X 2 t and the product rule gives 1 a d1 K t = 1 K t Xt 2 dt a dt 1 dw t = 1 K t X t Z t X t Xt 2 dt 1 X t dw t (1 a) akt dt 1 K t X t dw t. which can be written as dk t = 1 K t X 2 t a + a 1 dt + 1 K t dw t. K t X t As in 4 we consider the local martingale Mt = (1 x) 1 4a Xt 1 4a g t(1) 4a 1 = (1 x) 1 4a (1 K t ) 4a 1 Υ 1 4a t which satisfies dm t = 1 4a X t M t dw t M 0 = 1 11

12 If we use Girsanov and tilt by the local martingale we see that dk t = 1 K t a Xt 2 3a dt + 1 K t dwt. K t X t where Wt is a standard Brownian motion in the new measure P. We reparametrize so that log Υ t decays linearly. More precisely we let σ(t) = inf{t : Υ t = e at } and define ˆX t = X σ(t) Ŷt = Y σ(t) etc. Since ˆΥ t := Υ σ(t) = e at and a ˆΥ t = t ˆΥt = a ˆΥ t 1 ˆX 2 t 1 ˆK t ˆK t σ(t) we see that Therefore σ(t) = ˆX t 2 ˆK t 1 ˆK t t ˆΦ b t := Φ b σ(t) { ab = exp 1 ˆK } { t s ds = e 0 ˆK abt exp ab s 0 } 1 ds ˆK s d ˆK t = a 3a ˆK t ( = ˆK t aˆkt dt + ˆK t (1 ˆK t ) dbt ) 3a dt + 1 ˆK t ˆK t db t. for a standard Brownian motion B t (in the measure P ). Let λ = 2a 2 and { } { a(7a 1) N t = e λt a(3a 1) t ˆΦb t ˆKa t = exp t exp } 1 ds ˆK s ˆK a t. Itô s formula shows that N t is a local P -martingale satisfying 1 dn t = N t a ˆK t dbt ˆK N 0 = x a t One can show it is a martingale by using Girsanov to see that d ˆK t = 2a 4a ˆK t dt + ˆK t (1 ˆK t ) d B t where B t is a Brownian motion in the new measure P. By comparison with a Bessel process we see that the solution exists for all time. Equivalently we can say that ˆM t := ˆM t N t 12

13 is a P-martingale with ˆM 0 = x a. (Although Mt is only a local martingale the time-changed version ˆM t := Mσ(t) is a martingale.) Using (3) we see that E Φ b γ σ(t) c ˆKa t ˆΦb t. If ɛ = e at then E Φ b ; σ(t) < = c E E(Φ b 1{σ(t) < } γ σ(t) ) c E ˆKa t ˆΦb t ; σ(t) < = c e λt e (1 4a)at (1 x) 4a 1 1 ˆM 0 E ˆMt (1 ˆK t ) 1 4a ; σ(t) < = c e a(1 6a)t x a (1 x) 4a 1 Ẽ (1 ˆK t ) 1 4a = c ɛ 6a 1 x a (1 x) 4a 1 Ẽ (1 ˆK t ) 1 4a. So the result follows once we show that Ẽ (1 ˆK t ) 1 4a < is uniformly bounded for t t 0. The argument for this proceeds as in 4. If we do the change of variables ˆK t = 1 cos Θ t /2 then Itô s formula shows that ( dθ t = 4a 1 ) cot Θ t dt + db t. 2 This is a radial Bessel process that never reaches the boundary. It is known that the invariant distribution is proportional to sin 8a 1 θ and that it approaches the invariant distribution exponentially fast. One then computes that the invariant distribution for ˆK t is proportional to x 4a 1 (1 x) 4a 1. In particular (1 ˆK t ) 1 4a is integrable with respect to the invariant distribution. 4 Proof of Lemma 2.3 We prove the first part of lemma 2.3 for x 1 = 0 y 1 = 1. Other cases follow from this and a Möbius transformation sending x 1 y 1 to 0 1. Lemma 4.1. There exists c < such that if D is a simply connected subdomain of H containing 0 1 on its boundary then where δ = dist({0 1} D H). x H D (0 1) + y H D (0 1) c δ 1 H D (0 1) xx H D (0 1) + xy H D (0 1) + yy H D (0 1) c δ 2 H D (0 1) Proof. Let g : D H be a conformal transformation with g(0) = 0 g(1) = 1 g (0) = 1. Then if x < δ y 1 < δ H D (x y) = g (x) g (y) g(y) g(x) 2. (6) 13

14 In particular g (0) g (1) = H D (0 1) H H (0 1) = 1 and hence g (1) 1. Using Schwartz reflection we can extend g to be a conformal transformations of disks of radius δ about 0 and 1. By the distortion estimates (the fact that a 2 2 a 3 3 for schlicht functions) we have g (0) 4 δ 1 g (0) 4 δ 1 g (0) 18 δ 2 g (0) 18 δ 2 and similarly g (1) 4 δ 1 g (1) and g (1) 18 δ 1 g (1). right-hand side of (6) we get the result. By direct differentiation of the Lemma 4.2. There exists c < such that if x 1 < y 1 0 < 1 x 2 < y 2 Ψ D (x y) is as in (4) and z 1 {x 1 y 1 } z 2 {x 2 y 2 } then z1 ΨD (x y) + z2 ΨD (x y) c δ 1 Ψ D (x y) z1 z 2 ΨD (x y) c δ 2 ΨD (x y) where δ := min {{ w 1 w 2 ; w 1 w 2 and w 1 w 2 {x 1 x 2 y 1 y 2 }} dist {x 1 y 1 x 2 y 2 } H \ D}. Proof. Let Ψ D (x y) = φ (u D (x y)). where u D (x y) = e E D(xy) E D (x y) = y1 y2 Using the Harnack inequality we can see that for j = 1 2 x 1 x 2 H D (r s) dr ds. H D (x s) H D (x j s) H D (r y) H D (r y j ) if x x j δ/2 y y j δ/2. From this we see that y2 x 2 H D (z 1 s) ds + y1 Let z 1 be x 1 or y 1 and let z 2 be x 2 or y 2. Then x 1 H D (r z 2 ) ds c δ 1 E D (x y) H D (z 1 z 2 ) c δ 2 E D (x y). z1 ΨD (x y) = φ (u D (x y)) z1 u D (x y) z1 z 2 ΨD (x y) = φ (u D (x y)) z1 u D (x y) z2 u D (x y) + φ (u D (x y)) z1 z 2 u D (x y). y2 z1 u D (x y) = ± H D (z 1 s) ds u D (x y). z2 u D (x y) = ± x 2 y1 x 1 H D (r z 2 ) ds u D (x y). 14

15 This gives y1 y2 z2 z 1 u D (x y) = ± H D (r z 2 ) dr H D (z 1 s) ds ± H D (z 1 z 2 ) u D (x y) x 1 x 2 The result will follow if we show that z1 u D (x y) + z2 u D (x y) c δ 1 E D (x y) u D (x y) z2 z 1 u D (x y) c δ 2 E D (x y) u D (x y). z2 z 1 u D (x y) u D (x y) (1 x) φ (1 x) φ(x) c δ 2 E D (x y). φ (x) φ(x) are uniformly bounded for x > x 0. Recall that u(x) = c x a F (x) where F (x) = 2 F 1 (2a 1 2a 4a; x). We recall that F is analytic in the unit disk with power series expansion where the coefficients b j satisfy F (x) = 1 + b n x n n=1 b n = C n 4a O(n 1 ). We therefore get asymptotic expansions for the coefficients of the derivatives of F. The important thing for us is that if κ < 8 then 4a 1 > 1 and we have as x 1 In other words the quantities F (1 x) = O(1) F (x) = o(x 1 ) F (x)o(x 2 ). F (x) (1 x) F (x) F (x) (1 x) 2 F (x) F (x) are uniformly bounded for 0 x < 1. If g(x) = x a F (x) then a g (x) = g(x) x + F (x) F (x) ( a g (x) = g(x) x + F ) (x) 2 a F (x) x 2 + F (x) F (x) F (x) 2 F (x) 2. Therefore for every x 0 > 0 the quantities φ(x) (1 x) φ (x) φ(x) (1 x) 2 φ (x) φ(x) are uniformly bounded for x 0 < x < 1. 15

16 References 1 J. Dubédat (2006). Euler integrals for commuting SLEs J. Stat. Phys M. Kozdron G. Lawler (2007). The configurational measure pn mutuallly avoiding SLE paths in Universality and Renormalization: From Stochastic Evolution to Renormalization of Quantum Fields I. Binder D Kreimer ed. Amer. Math. Soc G. Lawler (2009). Schramm-Loewner Evolution in Statistical Mechanics S. Sheffield T Spenceer ed. IAS/Park City Mathematics Series G. Lawler (2015). Minkowski content of the intersection of a Schramm-Loewner evolution (SLE) curve with the real line J. Math Soc. Japan J. Miller S. Sheffield (2016). Imaginary geometry III: reversibilty of SLE κ for κ (4 8) Annals of Math O. Schramm (2000). Scaling limits for loop-erased random walks and uniform spanning trees Israel J. Math D. Zhan (2008). Reversibility of chordal SLE Annals of Prob Eveliina Peltola Hao Wu (2017). Global Multiple SLE κ for κ 4 and Connection Probabilities for Level Lines of GFF Preprint in arxiv: G. Lawler (2011). Defining SLE in multiply connected domains with the Brownian loop measure arxiv: G. Lawler (2009). Partition Functions Loop Measure and Versions of SLE Journal of Statistical Physics

Research Statement. Dapeng Zhan

Research Statement. Dapeng Zhan Research Statement Dapeng Zhan The Schramm-Loewner evolution (SLE), first introduced by Oded Schramm ([12]), is a oneparameter (κ (0, )) family of random non-self-crossing curves, which has received a

More information

SLE and CFT. Mitsuhiro QFT2005

SLE and CFT. Mitsuhiro QFT2005 SLE and CFT Mitsuhiro Kato @ QFT2005 1. Introduction Critical phenomena Conformal Field Theory (CFT) Algebraic approach, field theory, BPZ(1984) Stochastic Loewner Evolution (SLE) Geometrical approach,

More information

Conformal Invariance of the Exploration Path in 2D Critical Bond Percolation in the Square Lattice

Conformal Invariance of the Exploration Path in 2D Critical Bond Percolation in the Square Lattice Conformal Invariance of the Exploration Path in 2D Critical Bond Percolation in the Square Lattice Chinese University of Hong Kong, STAT December 12, 2012 (Joint work with Jonathan TSAI (HKU) and Wang

More information

M5MF6. Advanced Methods in Derivatives Pricing

M5MF6. Advanced Methods in Derivatives Pricing Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................

More information

The stochastic calculus

The stochastic calculus Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations

More information

Risk Neutral Measures

Risk Neutral Measures CHPTER 4 Risk Neutral Measures Our aim in this section is to show how risk neutral measures can be used to price derivative securities. The key advantage is that under a risk neutral measure the discounted

More information

Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs.

Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs. Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs Andrea Cosso LPMA, Université Paris Diderot joint work with Francesco Russo ENSTA,

More information

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce

More information

Rohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)

Rohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,

More information

Basic Arbitrage Theory KTH Tomas Björk

Basic Arbitrage Theory KTH Tomas Björk Basic Arbitrage Theory KTH 2010 Tomas Björk Tomas Björk, 2010 Contents 1. Mathematics recap. (Ch 10-12) 2. Recap of the martingale approach. (Ch 10-12) 3. Change of numeraire. (Ch 26) Björk,T. Arbitrage

More information

Hedging under Arbitrage

Hedging under Arbitrage Hedging under Arbitrage Johannes Ruf Columbia University, Department of Statistics Modeling and Managing Financial Risks January 12, 2011 Motivation Given: a frictionless market of stocks with continuous

More information

Change of Measure (Cameron-Martin-Girsanov Theorem)

Change of Measure (Cameron-Martin-Girsanov Theorem) Change of Measure Cameron-Martin-Girsanov Theorem Radon-Nikodym derivative: Taking again our intuition from the discrete world, we know that, in the context of option pricing, we need to price the claim

More information

Stochastic Dynamical Systems and SDE s. An Informal Introduction

Stochastic Dynamical Systems and SDE s. An Informal Introduction Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x

More information

PAPER 27 STOCHASTIC CALCULUS AND APPLICATIONS

PAPER 27 STOCHASTIC CALCULUS AND APPLICATIONS MATHEMATICAL TRIPOS Part III Thursday, 5 June, 214 1:3 pm to 4:3 pm PAPER 27 STOCHASTIC CALCULUS AND APPLICATIONS Attempt no more than FOUR questions. There are SIX questions in total. The questions carry

More information

arxiv:math/ v1 [math.pr] 24 Apr 2003

arxiv:math/ v1 [math.pr] 24 Apr 2003 ICM 2002 Vol. III 1 3 arxiv:math/0304369v1 [math.pr] 24 Apr 2003 Conformal Invariance, Universality, and the Dimension of the Brownian Frontier G. Lawler Abstract This paper describes joint work with Oded

More information

Lecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.

Lecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r. Lecture 7 Overture to continuous models Before rigorously deriving the acclaimed Black-Scholes pricing formula for the value of a European option, we developed a substantial body of material, in continuous

More information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Optimal stopping problems for a Brownian motion with a disorder on a finite interval

Optimal stopping problems for a Brownian motion with a disorder on a finite interval Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal

More information

Valuation of derivative assets Lecture 8

Valuation of derivative assets Lecture 8 Valuation of derivative assets Lecture 8 Magnus Wiktorsson September 27, 2018 Magnus Wiktorsson L8 September 27, 2018 1 / 14 The risk neutral valuation formula Let X be contingent claim with maturity T.

More information

Lecture 4. Finite difference and finite element methods

Lecture 4. Finite difference and finite element methods Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation

More information

SPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin

SPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin SPDE and portfolio choice (joint work with M. Musiela) Princeton University November 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 Performance measurement of investment strategies 2 Market

More information

Local vs Non-local Forward Equations for Option Pricing

Local vs Non-local Forward Equations for Option Pricing Local vs Non-local Forward Equations for Option Pricing Rama Cont Yu Gu Abstract When the underlying asset is a continuous martingale, call option prices solve the Dupire equation, a forward parabolic

More information

Brownian Motion. Richard Lockhart. Simon Fraser University. STAT 870 Summer 2011

Brownian Motion. Richard Lockhart. Simon Fraser University. STAT 870 Summer 2011 Brownian Motion Richard Lockhart Simon Fraser University STAT 870 Summer 2011 Richard Lockhart (Simon Fraser University) Brownian Motion STAT 870 Summer 2011 1 / 33 Purposes of Today s Lecture Describe

More information

Economathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t

Economathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t Economathematics Problem Sheet 1 Zbigniew Palmowski 1. Calculate Ee X where X is a gaussian random variable with mean µ and volatility σ >.. Verify that where W is a Wiener process. Ws dw s = 1 3 W t 3

More information

The ruin probabilities of a multidimensional perturbed risk model

The ruin probabilities of a multidimensional perturbed risk model MATHEMATICAL COMMUNICATIONS 231 Math. Commun. 18(2013, 231 239 The ruin probabilities of a multidimensional perturbed risk model Tatjana Slijepčević-Manger 1, 1 Faculty of Civil Engineering, University

More information

RMSC 4005 Stochastic Calculus for Finance and Risk. 1 Exercises. (c) Let X = {X n } n=0 be a {F n }-supermartingale. Show that.

RMSC 4005 Stochastic Calculus for Finance and Risk. 1 Exercises. (c) Let X = {X n } n=0 be a {F n }-supermartingale. Show that. 1. EXERCISES RMSC 45 Stochastic Calculus for Finance and Risk Exercises 1 Exercises 1. (a) Let X = {X n } n= be a {F n }-martingale. Show that E(X n ) = E(X ) n N (b) Let X = {X n } n= be a {F n }-submartingale.

More information

On the pricing equations in local / stochastic volatility models

On the pricing equations in local / stochastic volatility models On the pricing equations in local / stochastic volatility models Hao Xing Fields Institute/Boston University joint work with Erhan Bayraktar, University of Michigan Kostas Kardaras, Boston University Probability

More information

Multiname and Multiscale Default Modeling

Multiname and Multiscale Default Modeling Multiname and Multiscale Default Modeling Jean-Pierre Fouque University of California Santa Barbara Joint work with R. Sircar (Princeton) and K. Sølna (UC Irvine) Special Semester on Stochastics with Emphasis

More information

EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS

EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society

More information

Hints on Some of the Exercises

Hints on Some of the Exercises Hints on Some of the Exercises of the book R. Seydel: Tools for Computational Finance. Springer, 00/004/006/009/01. Preparatory Remarks: Some of the hints suggest ideas that may simplify solving the exercises

More information

Equivalence between Semimartingales and Itô Processes

Equivalence between Semimartingales and Itô Processes International Journal of Mathematical Analysis Vol. 9, 215, no. 16, 787-791 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.215.411358 Equivalence between Semimartingales and Itô Processes

More information

Chapter 3: Black-Scholes Equation and Its Numerical Evaluation

Chapter 3: Black-Scholes Equation and Its Numerical Evaluation Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random

More information

AMH4 - ADVANCED OPTION PRICING. Contents

AMH4 - ADVANCED OPTION PRICING. Contents AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5

More information

Risk Minimization Control for Beating the Market Strategies

Risk Minimization Control for Beating the Market Strategies Risk Minimization Control for Beating the Market Strategies Jan Večeř, Columbia University, Department of Statistics, Mingxin Xu, Carnegie Mellon University, Department of Mathematical Sciences, Olympia

More information

Generalized Affine Transform Formulae and Exact Simulation of the WMSV Model

Generalized Affine Transform Formulae and Exact Simulation of the WMSV Model On of Affine Processes on S + d Generalized Affine and Exact Simulation of the WMSV Model Department of Mathematical Science, KAIST, Republic of Korea 2012 SIAM Financial Math and Engineering joint work

More information

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.

More information

STOCHASTIC INTEGRALS

STOCHASTIC INTEGRALS Stat 391/FinMath 346 Lecture 8 STOCHASTIC INTEGRALS X t = CONTINUOUS PROCESS θ t = PORTFOLIO: #X t HELD AT t { St : STOCK PRICE M t : MG W t : BROWNIAN MOTION DISCRETE TIME: = t < t 1

More information

Characterization of the Optimum

Characterization of the Optimum ECO 317 Economics of Uncertainty Fall Term 2009 Notes for lectures 5. Portfolio Allocation with One Riskless, One Risky Asset Characterization of the Optimum Consider a risk-averse, expected-utility-maximizing

More information

The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations

The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.

More information

Market interest-rate models

Market interest-rate models Market interest-rate models Marco Marchioro www.marchioro.org November 24 th, 2012 Market interest-rate models 1 Lecture Summary No-arbitrage models Detailed example: Hull-White Monte Carlo simulations

More information

SHORT-TERM RELATIVE ARBITRAGE IN VOLATILITY-STABILIZED MARKETS

SHORT-TERM RELATIVE ARBITRAGE IN VOLATILITY-STABILIZED MARKETS SHORT-TERM RELATIVE ARBITRAGE IN VOLATILITY-STABILIZED MARKETS ADRIAN D. BANNER INTECH One Palmer Square Princeton, NJ 8542, USA adrian@enhanced.com DANIEL FERNHOLZ Department of Computer Sciences University

More information

Calculating Implied Volatility

Calculating Implied Volatility Statistical Laboratory University of Cambridge University of Cambridge Mathematics and Big Data Showcase 20 April 2016 How much is an option worth? A call option is the right, but not the obligation, to

More information

CONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES

CONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES CONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES D. S. SILVESTROV, H. JÖNSSON, AND F. STENBERG Abstract. A general price process represented by a two-component

More information

Chapter 3 Common Families of Distributions. Definition 3.4.1: A family of pmfs or pdfs is called exponential family if it can be expressed as

Chapter 3 Common Families of Distributions. Definition 3.4.1: A family of pmfs or pdfs is called exponential family if it can be expressed as Lecture 0 on BST 63: Statistical Theory I Kui Zhang, 09/9/008 Review for the previous lecture Definition: Several continuous distributions, including uniform, gamma, normal, Beta, Cauchy, double exponential

More information

Saddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance

Saddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is

More information

Constructing Markov models for barrier options

Constructing Markov models for barrier options Constructing Markov models for barrier options Gerard Brunick joint work with Steven Shreve Department of Mathematics University of Texas at Austin Nov. 14 th, 2009 3 rd Western Conference on Mathematical

More information

D MATH Departement of Mathematics Finite dimensional realizations for the CNKK-volatility surface model

D MATH Departement of Mathematics Finite dimensional realizations for the CNKK-volatility surface model Finite dimensional realizations for the CNKK-volatility surface model Josef Teichmann Outline 1 Introduction 2 The (generalized) CNKK-approach 3 Affine processes as generic example for the CNNK-approach

More information

1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components:

1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components: 1 Mathematics in a Pill The purpose of this chapter is to give a brief outline of the probability theory underlying the mathematics inside the book, and to introduce necessary notation and conventions

More information

On Complexity of Multistage Stochastic Programs

On Complexity of Multistage Stochastic Programs On Complexity of Multistage Stochastic Programs Alexander Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0205, USA e-mail: ashapiro@isye.gatech.edu

More information

Lecture Notes for Chapter 6. 1 Prototype model: a one-step binomial tree

Lecture Notes for Chapter 6. 1 Prototype model: a one-step binomial tree Lecture Notes for Chapter 6 This is the chapter that brings together the mathematical tools (Brownian motion, Itô calculus) and the financial justifications (no-arbitrage pricing) to produce the derivative

More information

1 The continuous time limit

1 The continuous time limit Derivative Securities, Courant Institute, Fall 2008 http://www.math.nyu.edu/faculty/goodman/teaching/derivsec08/index.html Jonathan Goodman and Keith Lewis Supplementary notes and comments, Section 3 1

More information

Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes

Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes Fabio Trojani Department of Economics, University of St. Gallen, Switzerland Correspondence address: Fabio Trojani,

More information

BACHELIER FINANCE SOCIETY. 4 th World Congress Tokyo, Investments and forward utilities. Thaleia Zariphopoulou The University of Texas at Austin

BACHELIER FINANCE SOCIETY. 4 th World Congress Tokyo, Investments and forward utilities. Thaleia Zariphopoulou The University of Texas at Austin BACHELIER FINANCE SOCIETY 4 th World Congress Tokyo, 26 Investments and forward utilities Thaleia Zariphopoulou The University of Texas at Austin 1 Topics Utility-based measurement of performance Utilities

More information

25857 Interest Rate Modelling

25857 Interest Rate Modelling 25857 UTS Business School University of Technology Sydney Chapter 20. Change of Numeraire May 15, 2014 1/36 Chapter 20. Change of Numeraire 1 The Radon-Nikodym Derivative 2 Option Pricing under Stochastic

More information

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets Chapter 5: Jump Processes and Incomplete Markets Jumps as One Explanation of Incomplete Markets It is easy to argue that Brownian motion paths cannot model actual stock price movements properly in reality,

More information

Large Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models

Large Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied

More information

Optimal trading strategies under arbitrage

Optimal trading strategies under arbitrage Optimal trading strategies under arbitrage Johannes Ruf Columbia University, Department of Statistics The Third Western Conference in Mathematical Finance November 14, 2009 How should an investor trade

More information

Drawdowns Preceding Rallies in the Brownian Motion Model

Drawdowns Preceding Rallies in the Brownian Motion Model Drawdowns receding Rallies in the Brownian Motion Model Olympia Hadjiliadis rinceton University Department of Electrical Engineering. Jan Večeř Columbia University Department of Statistics. This version:

More information

Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model

Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model 1(23) Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility

More information

MSc Financial Engineering CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL. To be handed in by monday January 28, 2013

MSc Financial Engineering CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL. To be handed in by monday January 28, 2013 MSc Financial Engineering 2012-13 CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL To be handed in by monday January 28, 2013 Department EMS, Birkbeck Introduction The assignment consists of Reading

More information

1.1 Basic Financial Derivatives: Forward Contracts and Options

1.1 Basic Financial Derivatives: Forward Contracts and Options Chapter 1 Preliminaries 1.1 Basic Financial Derivatives: Forward Contracts and Options A derivative is a financial instrument whose value depends on the values of other, more basic underlying variables

More information

Logarithmic derivatives of densities for jump processes

Logarithmic derivatives of densities for jump processes Logarithmic derivatives of densities for jump processes Atsushi AKEUCHI Osaka City University (JAPAN) June 3, 29 City University of Hong Kong Workshop on Stochastic Analysis and Finance (June 29 - July

More information

Exam Quantitative Finance (35V5A1)

Exam Quantitative Finance (35V5A1) Exam Quantitative Finance (35V5A1) Part I: Discrete-time finance Exercise 1 (20 points) a. Provide the definition of the pricing kernel k q. Relate this pricing kernel to the set of discount factors D

More information

Optimal Allocation of Policy Limits and Deductibles

Optimal Allocation of Policy Limits and Deductibles Optimal Allocation of Policy Limits and Deductibles Ka Chun Cheung Email: kccheung@math.ucalgary.ca Tel: +1-403-2108697 Fax: +1-403-2825150 Department of Mathematics and Statistics, University of Calgary,

More information

Variance Reduction for Monte Carlo Simulation in a Stochastic Volatility Environment

Variance Reduction for Monte Carlo Simulation in a Stochastic Volatility Environment Variance Reduction for Monte Carlo Simulation in a Stochastic Volatility Environment Jean-Pierre Fouque Tracey Andrew Tullie December 11, 21 Abstract We propose a variance reduction method for Monte Carlo

More information

Optimal asset allocation under forward performance criteria Oberwolfach, February 2007

Optimal asset allocation under forward performance criteria Oberwolfach, February 2007 Optimal asset allocation under forward performance criteria Oberwolfach, February 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 References Indifference valuation in binomial models (with

More information

An Introduction to Point Processes. from a. Martingale Point of View

An Introduction to Point Processes. from a. Martingale Point of View An Introduction to Point Processes from a Martingale Point of View Tomas Björk KTH, 211 Preliminary, incomplete, and probably with lots of typos 2 Contents I The Mathematics of Counting Processes 5 1 Counting

More information

GPD-POT and GEV block maxima

GPD-POT and GEV block maxima Chapter 3 GPD-POT and GEV block maxima This chapter is devoted to the relation between POT models and Block Maxima (BM). We only consider the classical frameworks where POT excesses are assumed to be GPD,

More information

Replication under Price Impact and Martingale Representation Property

Replication under Price Impact and Martingale Representation Property Replication under Price Impact and Martingale Representation Property Dmitry Kramkov joint work with Sergio Pulido (Évry, Paris) Carnegie Mellon University Workshop on Equilibrium Theory, Carnegie Mellon,

More information

Laurence Boxer and Ismet KARACA

Laurence Boxer and Ismet KARACA SOME PROPERTIES OF DIGITAL COVERING SPACES Laurence Boxer and Ismet KARACA Abstract. In this paper we study digital versions of some properties of covering spaces from algebraic topology. We correct and

More information

1 Implied Volatility from Local Volatility

1 Implied Volatility from Local Volatility Abstract We try to understand the Berestycki, Busca, and Florent () (BBF) result in the context of the work presented in Lectures and. Implied Volatility from Local Volatility. Current Plan as of March

More information

A NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS. Burhaneddin İZGİ

A NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS. Burhaneddin İZGİ A NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS Burhaneddin İZGİ Department of Mathematics, Istanbul Technical University, Istanbul, Turkey

More information

On Using Shadow Prices in Portfolio optimization with Transaction Costs

On Using Shadow Prices in Portfolio optimization with Transaction Costs On Using Shadow Prices in Portfolio optimization with Transaction Costs Johannes Muhle-Karbe Universität Wien Joint work with Jan Kallsen Universidad de Murcia 12.03.2010 Outline The Merton problem The

More information

Monte Carlo Methods for Uncertainty Quantification

Monte Carlo Methods for Uncertainty Quantification Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)

More information

Log-linear Dynamics and Local Potential

Log-linear Dynamics and Local Potential Log-linear Dynamics and Local Potential Daijiro Okada and Olivier Tercieux [This version: November 28, 2008] Abstract We show that local potential maximizer ([15]) with constant weights is stochastically

More information

KØBENHAVNS UNIVERSITET (Blok 2, 2011/2012) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours

KØBENHAVNS UNIVERSITET (Blok 2, 2011/2012) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours This question paper consists of 3 printed pages FinKont KØBENHAVNS UNIVERSITET (Blok 2, 211/212) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours This exam paper

More information

Continuous Time Finance. Tomas Björk

Continuous Time Finance. Tomas Björk Continuous Time Finance Tomas Björk 1 II Stochastic Calculus Tomas Björk 2 Typical Setup Take as given the market price process, S(t), of some underlying asset. S(t) = price, at t, per unit of underlying

More information

Collinear Triple Hypergraphs and the Finite Plane Kakeya Problem

Collinear Triple Hypergraphs and the Finite Plane Kakeya Problem Collinear Triple Hypergraphs and the Finite Plane Kakeya Problem Joshua Cooper August 14, 006 Abstract We show that the problem of counting collinear points in a permutation (previously considered by the

More information

Are the Azéma-Yor processes truly remarkable?

Are the Azéma-Yor processes truly remarkable? Are the Azéma-Yor processes truly remarkable? Jan Obłój j.obloj@imperial.ac.uk based on joint works with L. Carraro, N. El Karoui, A. Meziou and M. Yor Welsh Probability Seminar, 17 Jan 28 Are the Azéma-Yor

More information

Dr. Maddah ENMG 625 Financial Eng g II 10/16/06

Dr. Maddah ENMG 625 Financial Eng g II 10/16/06 Dr. Maddah ENMG 65 Financial Eng g II 10/16/06 Chapter 11 Models of Asset Dynamics () Random Walk A random process, z, is an additive process defined over times t 0, t 1,, t k, t k+1,, such that z( t )

More information

MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.

MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS. MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS May/June 2006 Time allowed: 2 HOURS. Examiner: Dr N.P. Byott This is a CLOSED

More information

Optimal robust bounds for variance options and asymptotically extreme models

Optimal robust bounds for variance options and asymptotically extreme models Optimal robust bounds for variance options and asymptotically extreme models Alexander Cox 1 Jiajie Wang 2 1 University of Bath 2 Università di Roma La Sapienza Advances in Financial Mathematics, 9th January,

More information

Computational Finance

Computational Finance Path Dependent Options Computational Finance School of Mathematics 2018 The Random Walk One of the main assumption of the Black-Scholes framework is that the underlying stock price follows a random walk

More information

Risk Neutral Valuation

Risk Neutral Valuation copyright 2012 Christian Fries 1 / 51 Risk Neutral Valuation Christian Fries Version 2.2 http://www.christian-fries.de/finmath April 19-20, 2012 copyright 2012 Christian Fries 2 / 51 Outline Notation Differential

More information

Asian Options under Multiscale Stochastic Volatility

Asian Options under Multiscale Stochastic Volatility Contemporary Mathematics Asian Options under Multiscale Stochastic Volatility Jean-Pierre Fouque and Chuan-Hsiang Han Abstract. We study the problem of pricing arithmetic Asian options when the underlying

More information

Hedging under arbitrage

Hedging under arbitrage Hedging under arbitrage Johannes Ruf Columbia University, Department of Statistics AnStAp10 August 12, 2010 Motivation Usually, there are several trading strategies at one s disposal to obtain a given

More information

Local Volatility Dynamic Models

Local Volatility Dynamic Models René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Columbia November 9, 27 Contents Joint work with Sergey Nadtochyi Motivation 1 Understanding

More information

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question

More information

A note on the existence of unique equivalent martingale measures in a Markovian setting

A note on the existence of unique equivalent martingale measures in a Markovian setting Finance Stochast. 1, 251 257 1997 c Springer-Verlag 1997 A note on the existence of unique equivalent martingale measures in a Markovian setting Tina Hviid Rydberg University of Aarhus, Department of Theoretical

More information

Risk, Return, and Ross Recovery

Risk, Return, and Ross Recovery Risk, Return, and Ross Recovery Peter Carr and Jiming Yu Courant Institute, New York University September 13, 2012 Carr/Yu (NYU Courant) Risk, Return, and Ross Recovery September 13, 2012 1 / 30 P, Q,

More information

USC Math. Finance April 22, Path-dependent Option Valuation under Jump-diffusion Processes. Alan L. Lewis

USC Math. Finance April 22, Path-dependent Option Valuation under Jump-diffusion Processes. Alan L. Lewis USC Math. Finance April 22, 23 Path-dependent Option Valuation under Jump-diffusion Processes Alan L. Lewis These overheads to be posted at www.optioncity.net (Publications) Topics Why jump-diffusion models?

More information

IEOR 3106: Introduction to Operations Research: Stochastic Models SOLUTIONS to Final Exam, Sunday, December 16, 2012

IEOR 3106: Introduction to Operations Research: Stochastic Models SOLUTIONS to Final Exam, Sunday, December 16, 2012 IEOR 306: Introduction to Operations Research: Stochastic Models SOLUTIONS to Final Exam, Sunday, December 6, 202 Four problems, each with multiple parts. Maximum score 00 (+3 bonus) = 3. You need to show

More information

M.I.T Fall Practice Problems

M.I.T Fall Practice Problems M.I.T. 15.450-Fall 2010 Sloan School of Management Professor Leonid Kogan Practice Problems 1. Consider a 3-period model with t = 0, 1, 2, 3. There are a stock and a risk-free asset. The initial stock

More information

Stochastic Calculus - An Introduction

Stochastic Calculus - An Introduction Stochastic Calculus - An Introduction M. Kazim Khan Kent State University. UET, Taxila August 15-16, 17 Outline 1 From R.W. to B.M. B.M. 3 Stochastic Integration 4 Ito s Formula 5 Recap Random Walk Consider

More information

PDE Approach to Credit Derivatives

PDE Approach to Credit Derivatives PDE Approach to Credit Derivatives Marek Rutkowski School of Mathematics and Statistics University of New South Wales Joint work with T. Bielecki, M. Jeanblanc and K. Yousiph Seminar 26 September, 2007

More information

A utility maximization proof of Strassen s theorem

A utility maximization proof of Strassen s theorem Introduction CMAP, Ecole Polytechnique Paris Advances in Financial Mathematics, Paris January, 2014 Outline Introduction Notations Strassen s theorem 1 Introduction Notations Strassen s theorem 2 General

More information

arxiv: v1 [q-fin.pm] 13 Mar 2014

arxiv: v1 [q-fin.pm] 13 Mar 2014 MERTON PORTFOLIO PROBLEM WITH ONE INDIVISIBLE ASSET JAKUB TRYBU LA arxiv:143.3223v1 [q-fin.pm] 13 Mar 214 Abstract. In this paper we consider a modification of the classical Merton portfolio optimization

More information

Girsanov s Theorem. Bernardo D Auria web: July 5, 2017 ICMAT / UC3M

Girsanov s Theorem. Bernardo D Auria   web:   July 5, 2017 ICMAT / UC3M Girsanov s Theorem Bernardo D Auria email: bernardo.dauria@uc3m.es web: www.est.uc3m.es/bdauria July 5, 2017 ICMAT / UC3M Girsanov s Theorem Decomposition of P-Martingales as Q-semi-martingales Theorem

More information

The Lognormal Interest Rate Model and Eurodollar Futures

The Lognormal Interest Rate Model and Eurodollar Futures GLOBAL RESEARCH ANALYTICS The Lognormal Interest Rate Model and Eurodollar Futures CITICORP SECURITIES,INC. 399 Park Avenue New York, NY 143 Keith Weintraub Director, Analytics 1-559-97 Michael Hogan Ex

More information

Disaster risk and its implications for asset pricing Online appendix

Disaster risk and its implications for asset pricing Online appendix Disaster risk and its implications for asset pricing Online appendix Jerry Tsai University of Oxford Jessica A. Wachter University of Pennsylvania December 12, 2014 and NBER A The iid model This section

More information