Short dated option pricing under rough volatility

Size: px
Start display at page:

Download "Short dated option pricing under rough volatility"

Transcription

1 Department of Mathematics, Imperial College London Mathematics of Quantitative Finance Oberwolfach, February 2017 Based on joint works with P. Harms, A. Jacquier, C. Lacombe and with C. Bayer, P. Friz, A. Gulisashvili and B. Stemper

2 Implied volatility Asset price process: S t = e Xt ) t 0, with X 0 = 0. Black-Scholes-Merton BSM) framework: C BS τ, k, σ) := E 0 e Xτ e k) = N d+) + ek N d ), d ± := k σ τ ± 1 2 σ τ. Spot implied volatility σ τ k): the unique non-negative) solution to C observed τ, k) = C BS τ, k, σ τ k)). Implied volatility: unit-free measure of option prices.

3 Implied volatility Asset price process: S t = e Xt ) t 0, with X 0 = 0. Black-Scholes-Merton BSM) framework: C BS τ, k, σ) := E 0 e Xτ e k) = N d+) + ek N d ), d ± := k σ τ ± 1 2 σ τ. Spot implied volatility σ τ k): the unique non-negative) solution to C observed τ, k) = C BS τ, k, σ τ k)). Implied volatility: unit-free measure of option prices. Implied volatility is not available in closed form generally. Its asymptotic behaviour is available via small/large k, τ) approximations.

4 Literature Implied volatility σ τ k)) asymptotics as k, τ 0 or τ : Hagan-Kumar-Lesniewski-Woodward 2003/2015), Ob lój 2008): small-maturity for the SABR model. Berestycki-Busca-Florent 2004): small-τ using PDE methods for diffusions. Henry-Labordère 2009): small-τ asymptotics using differential geometry. Forde et al.2012), Jacquier et al.2012): small/ large τ using large deviations. Lee 2003), Benaim-Friz 2009), Gulisashvili ), Caravenna-Corbetta 2016), De Marco-Jacquier-Hillairet 2013): k. Laurence-Gatheral-Hsu-Ouyang-Wang 2012): small-τ in local volatility models. Mijatović-Tankov 2012): small-τ for jump models. Bompis-Gobet 2015): asymptotic expansions in the presence of both local and stochastic volatility using Malliavin calculus. Fouque et al ): perturbation techniques for slow and fast mean-reverting stochastic volatility models.

5 Related works: Deuschel-Friz-Jacquier-Violante CPAM 2014), De Marco-Friz 2014): small-noise expansions using Laplace method on Wiener space Ben Arous-Bismut approach). Baudoin-Ouyang 2015): small-noise expansions in a fully) fractional setting Gatheral-Jaisson-Rosenbaum 2014), and Bayer-Gatheral-Friz 2015), Forde-Zhang 2015): large deviations in a fractional stochastic volatility setting Fukasawa 2011,2015), Alós-León-Vives 2007) small-time fractional) skew Guennoun-Jacquier-Roome 2015), El Euch-Rosenbaum 2016) fractional Heston.

6 Classical case: H i = 1 ) case: Deuschel-Friz-Jacquier-Violante 2011) 2 dx ε t = bε, X ε t )dt + ε m i=1 σ i X ε t )dw i t, X ε 0 = xε 0 Rd

7 Classical case: H i = 1 ) case: Deuschel-Friz-Jacquier-Violante 2011) 2 dx ε t = bε, X ε t )dt + ε m i=1 σ i X ε t )dw i t, X ε 0 = xε 0 Rd Fractional case: H i = H 1, 1)) case: Baudoin-Ouyang 2015) 4 m dx ε t = bε, Xε t )dt + ε σ i X ε t )dw H ) i t, i=1 Xε 0 = xε 0 Rd

8 Classical case: H i = 1 ) case: Deuschel-Friz-Jacquier-Violante 2011) 2 dx ε t = bε, X ε t )dt + ε m i=1 σ i X ε t )dw i t, X ε 0 = xε 0 Rd Fractional case: H i = H 1, 1)) case: Baudoin-Ouyang 2015) 4 m dx ε t = bε, Xε t )dt + ε σ i X ε t )dw H ) i t, i=1 Xε 0 = xε 0 Rd Our main interest: H 1 = 1 2, H

9 Rough volatility models Short-term data suggests a time decay of the ATM skew proportional to τ θ, with θ 0, 1/2) while classical stochastic volatility models generate a constant short-maturity skew. Gatheral-Jaisson-Rosenbaum and Bayer-Gatheral-Friz 2014,1015) proposed a fractional volatility model: ds t = S tσ tdz t + µ tdt), σ t = exp Y t), 1) where dy t = µdw H t by t m)dt, for µ, b > 0, m R for a Bm Z and a fbm motion W H with Hurst parameter H. Time series of the Oxford-Man SPX realised variance as well as implied volatility smiles of the SPX suggest that H 0, 1/2): short-memory volatility. Main drawback: loss of Markovianity H 1/2) rules out PDE techniques, and Monte Carlo is computationally intensive. One way out is an efficient Hybrid scheme of Bennedsen, Lunde and Pakkanen 2015).

10 Rough volatility models: our setting dx t = b 1 Y t)dt + σ 1 Y t)dw t dy t = b 2 Y t)dt + σ 2 dw H t. 2) σ 2 > 0 extendible to bounded, elliptic) and H 0, 1), particular interest in H < 1/2. Forde-Zhang 16: b 1 b 2 0, and σ 1 C α, α 0, 1] second part of the talk); Fractional Stein-Stein: b 1 y) y 2 /2, σ 1 y) y, b 2 y) a by first part). Gatheral-Jaisson-Rosenbaum 14: b 1 y) e 2y /2, σ 1 y) e y, b 2 y) a by.

11 Rough volatility models: our setting dx t = b 1 Y t)dt + σ 1 Y t)dw t dy t = b 2 Y t)dt + σ 2 dw H t. 2) σ 2 > 0 extendible to bounded, elliptic) and H 0, 1), particular interest in H < 1/2. Forde-Zhang 16: b 1 b 2 0, and σ 1 C α, α 0, 1] second part of the talk); Fractional Stein-Stein: b 1 y) y 2 /2, σ 1 y) y, b 2 y) a by first part). Gatheral-Jaisson-Rosenbaum 14: b 1 y) e 2y /2, σ 1 y) e y, b 2 y) a by. To introduce correlation we consider B and B independent, and set W = ρ B + ρb and t Wt H = Kt, s)db s 0 where K the Volterra kernel of the standard) fbm W H.

12 Rough volatility models: our setting dx t = b 1 Y t)dt + σ 1 Y t)dw t dy t = b 2 Y t)dt + σ 2 dw H t. 2) σ 2 > 0 extendible to bounded, elliptic) and H 0, 1), particular interest in H < 1/2. Forde-Zhang 16: b 1 b 2 0, and σ 1 C α, α 0, 1] second part of the talk); Fractional Stein-Stein: b 1 y) y 2 /2, σ 1 y) y, b 2 y) a by first part). Gatheral-Jaisson-Rosenbaum 14: b 1 y) e 2y /2, σ 1 y) e y, b 2 y) a by. To introduce correlation we consider B and B independent, and set W = ρ B + ρb and t Wt H = Kt, s)db s 0 where K the Volterra kernel of the standard) fbm W H. An intuitive remark: Mandelbrot-van Ness representation for fbm: At time zero, the volatility process in 2) has already accumulated some randomness. Hot-start of the process Y with A. Jacquier and C. Lacombe).

13 1 Introduction 2 3 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics

14 Varadhan-type asymptotics Recall the Black-Scholes density expansion: heat-kernel asymptotics f BS t, x) t 1/2 exp 1 x ) ) 2, as t 0. 2t σ In the homogenous fractional) case: H 1 = H 2 =... = H d = H Baudoin-Ouyang) dx t = bx t)dt + Extended) Varadhan formula m σ i X t)dw H ) i t, X 0 = x 0 R d 3) i=1 ) f X t, x) cst t H exp d2 x 0, x) 2t 2H, as t 0.

15 Varadhan-type asymptotics Recall the Black-Scholes density expansion: heat-kernel asymptotics f BS t, x) t 1/2 exp 1 x ) ) 2, as t 0. 2t σ In the homogenous fractional) case: H 1 = H 2 =... = H d = H Baudoin-Ouyang) dx t = bx t)dt + Extended) Varadhan formula m σ i X t)dw H ) i t, X 0 = x 0 R d 3) i=1 ) f X t, x) cst t H exp d2 x 0, x) 2t 2H, as t 0. What if Hurst parameters are different H 1 H 2?

16 Rescalings Key to asymptotic expansions for H 1 H 2 : Rescalings Recall the considered processes dx t = b 1 Y t)dt + σ 1 Y t)dw t dy t = b 2 Y t)dt + σ 2 dw H t.

17 Rescalings Key to asymptotic expansions for H 1 H 2 : Rescalings Recall the considered processes Define appropriate rescalings dx ε t dy ε t dx t = b 1 Y t)dt + σ 1 Y t)dw t dy t = b 2 Y t)dt + σ 2 dw H t. = b 1 ε κ 1, Y ε t )dt + εβ σ 1 Y ε t )dwt = b 2 ε κ 2, Y ε t )dt + ε β σ 2 dw H t.

18 Rescalings dx ε t = b 1 ε κ 1, Y ε t )dt + ε β σ 1 Y ε t )dw t dy ε t = b 2 ε κ 2, Y ε t )dt + ε β σ 2 dw H t. 4) Fractional Stein-Stein: Consider X 0, Y 0 ) = 0, y 0 ) and dx t = Y t 2 H dt + YtdWt, dyt = a byt)dt + cdwt 2. 5)

19 Rescalings dx ε t = b 1 ε κ 1, Y ε t )dt + ε β σ 1 Y ε t )dw t dy ε t = b 2 ε κ 2, Y ε t )dt + ε β σ 2 dw H t. 4) Fractional Stein-Stein: Consider X 0, Y 0 ) = 0, y 0 ) and dx t = Y t 2 H dt + YtdWt, dyt = a byt)dt + cdwt 2. 5) Rescaling 1 short-time): Xt ε, Y t ε) := ε2h 1 X ε 2 t, Y ε 2 t ) 4) with κ 1 = 2H + 1, κ 2 = 2 and β = 2H, x0 ε, y 0 ε) = 0, y 0): dx ε t 2H+1 Y ε t = ε )2 dt + ε 2H Yt ε 2 dwt, dy ε t = ε 2 a by ε t )dt + ε2h cdw H t. 6)

20 Rescalings dx ε t = b 1 ε κ 1, Y ε t )dt + ε β σ 1 Y ε t )dw t dy ε t = b 2 ε κ 2, Y ε t )dt + ε β σ 2 dw H t. 4) Fractional Stein-Stein: Consider X 0, Y 0 ) = 0, y 0 ) and dx t = Y t 2 H dt + YtdWt, dyt = a byt)dt + cdwt 2. 5) Rescaling 1 short-time): Xt ε, Y t ε) := ε2h 1 X ε 2 t, Y ε 2 t ) 4) with κ 1 = 2H + 1, κ 2 = 2 and β = 2H, x0 ε, y 0 ε) = 0, y 0): dx ε t 2H+1 Y ε t = ε )2 dt + ε 2H Yt ε 2 dwt, dy ε t = ε 2 a by ε t )dt + ε2h cdw H t. Rescaling 2 tails): X ε t, Y ε t ) := ε2h X t, ε H Y t) 4) with κ 1 = 0, κ 2 = β = H, x ε 0, y ε 0 ) = 0, εh y 0 ): dx ε t = Y t ε)2 dt + ε H Yt ε 2 dwt, dy ε t 6) = aε H + by ε t )dt + εh cdw H t. 7)

21 Theorem Harms-H-Jacquier) Consider an SDE of the form 4). Then the density of XT ε admits an expansion f εt, x) = exp Λx) ε 2β + Λ x) X ) ) T ε β ε minκ 1,β) c 0 + Oε δκ1,β) ), as ε 0, where and d X ) t = [ x b 1 0, φ k 0 t + x σ 1 { } 1 Λx) = inf 2 k 2 H H, k K x x0 0,y 0 0 = 1 2 k 0 2 H H, φ h 0 t ) k ] 0 t) Xtdt+ ε β b 1 0, φ k 0 t ) dt, X0 = ε β x ε 0 ε=0, where φ k 0 denotes the ODE solution of the same SDE 4) replacing ε β dw by k 0 and x0 ε by x0 0.

22 Notations H: absolutely continuous paths [0, T ] R 2 starting at 0 such that ḣ 2 <. H H H := K H H and k := K H h, where K H denotes the Volterra kernel. For fixed x 0, y 0 ) R 2, φ k is the unique) ODE solution to ) ) ) φ k t = b 1 0, φ k t dt + σ 1 φ h t dk 1 t + σ 2 φ h t dk 2 t, φk 0 = x0 0, y 0 0 ). Denote ψ k := Π 1 φ k its projection on to the first coordinate X. K a := { k H H : ψt k = a R} by Hörmander condition ). { } Λa) := inf 1 2 k 2 H : k Ka.

23 Corollary short-time asymptotics in Stein-Stein) Corollary: Varadhan-type asymptotics dy t = a by t)dt + cdw H t In the fractional Stein-Stein model X t, Y t) with X 0 = 0, Y 0 = y 0 > 0 the density of X t satisfies in a neighbourhood of x 0, y 0 ) the following asymptotic expansion as t 0 f X t, x) = exp Λx) ) ) 1 t 2H t H 2π + OtδH,H+1/2) ) where { } 1 Λx) = inf 2 k 2 H H, k Kx x 0,y 0.

24 Corollary short-time asymptotics in Stein-Stein) Corollary: Varadhan-type asymptotics dy t = a by t)dt + cdw H t In the fractional Stein-Stein model X t, Y t) with X 0 = 0, Y 0 = y 0 > 0 the density of X t satisfies in a neighbourhood of x 0, y 0 ) the following asymptotic expansion as t 0 f X t, x) = exp Λx) ) ) 1 t 2H t H 2π + OtδH,H+1/2) ) where { } 1 Λx) = inf 2 k 2 H H, k Kx x 0,y 0. Proof: Take T = 1, ε 2 = t and consider Xt ε, Y t ε) := ε2h 1 X ε 2 t, Y ε 2 t ) with X0 ε = 0, Y 0 ε = y 0 > 0. Short-time scaling: dx ε t 2H+1 Y ε t = ε )2 dt + ε 2H Yt ε 2 dwt, dy ε t = ε 2 a + by ε t )dt + ε2h cdw H t, 6) Note that the drift vanishes in the limit ε 0 and x ε 0 = x 0 = 0. X t, Ŷt) 0, so that there is no 1/εβ = 1/t β/2 term in the exponential.

25 Corollary tail expansion in Stein-Stein) Corollary: tail asymptotics dy t = a + by t)dt + cdw H t Consider the fractional Stein-Stein model with X 0 = 0, Y 0 = y 0 > 0. Then as x, f X T, x) = exp c 1 x + c 2 x 1/2) 1 x 1/2 c 0 + O x 1/2)) where c 1 := Λ1), c 2 := X T Λ 1). Note that the expression on the RHS is independent of the Hurst-parameter!

26 Corollary tail expansion in Stein-Stein) Corollary: tail asymptotics dy t = a + by t)dt + cdw H t Consider the fractional Stein-Stein model with X 0 = 0, Y 0 = y 0 > 0. Then as x, f X T, x) = exp c 1 x + c 2 x 1/2) 1 x 1/2 c 0 + O x 1/2)) where c 1 := Λ1), c 2 := X T Λ 1). Note that the expression on the RHS is independent of the Hurst-parameter! Proof: Consider X ε T, Y ε T ) := ε2h X T, ε H Y T ) with X ε 0 = ε2h X 0 and Y ε 0 = εh Y 0. dx ε t = Y t ε)2 dt + ε H Yt ε 2 dwt, dy ε t = aε H + by ε t )dt + εh cdw H t, 7) Note that XT ε = ε 2H X T. PXT ε y) = PX T y/ε 2H ), f X T, y/ε 2H ) = ε 2H f εt, y). ) Take y = 1, that is x := ε 2H. By the theorem, f εt, 1) exp Λ1) ε H ) ε, hence H f X T, x) exp Λ1) ε +... ε H = exp Λ1)x +...) 2H 1 x 1/2.

27 From density to implied volatility: small-time Recall the Black-Scholes density expansion: f BS t, x) t 1/2 exp 1 x ) ) 2, as t 0, for any x R. 2t σ The corollary Varadhan-type asymptotics) implies that in the fractional Stein-Stein model ) f X t, x) cst t H exp d2 x 0, y 0 ; x) 2t 2H, as t 0.

28 From density to implied volatility: small-time Recall the Black-Scholes density expansion: f BS t, x) t 1/2 exp 1 x ) ) 2, as t 0, for any x R. 2t σ The corollary Varadhan-type asymptotics) implies that in the fractional Stein-Stein model ) f X t, x) cst t H exp d2 x 0, y 0 ; x) 2t 2H, as t 0. Matching the leading-orders gives σ BS t, x) x dx 0, y 0 ; x) th 1/2 as t 0. Skew explodes with rate H 1/2 in the short end whenever H < 1/2.

29 From density to implied volatility: tails Recall the Black-Scholes density expansion: f BS t, x) exp x2 2σ 2 t x ) 4 as x, for any t > 0. Our theorem corollary) says that in the fractional Stein-Stein model 5), we have f X t, x) cst x 1/2 exp c 1 x + c 2 x ), as x.

30 From density to implied volatility: tails Recall the Black-Scholes density expansion: f BS t, x) exp x2 2σ 2 t x ) 4 as x, for any t > 0. Our theorem corollary) says that in the fractional Stein-Stein model 5), we have f X t, x) cst x 1/2 exp c 1 x + c 2 x ), as x. Matching the leading-orders gives x 2 c 1 x + c 2 x 2σ 2 t x 4, and we recover Roger Lee s formula independently of the Hurst exponent in 5).

31 Proof of Theorem 1 dx t = ε 2H Y 2 t dt + ε 2H Y tdw t, dy t = ε 2H dw H t, with the same initial condition X 0 = Y 0 = 0. [ Density: f εt, x) = exp Λx) ε 4H + Λ x) X ] ) T ε 2H ε 2H c 0 + Oε 2H ).

32 Proof of Theorem 1 dx t = ε 2H Y 2 t dt + ε 2H Y tdw t, dy t = ε 2H dw H t, with the same initial condition X 0 = Y 0 = 0. [ Density: f εt, x) = exp Λx) ε 4H + Λ x) X ] ) T ε 2H ε 2H c 0 + Oε 2H ). Proof: Take x R and a C -bounded function F such that F x) = 0. f εt, x)e F x)/ε4h = 1 { [ 2πε 2H E exp iζ, 0) R X ε T x, 0) ε 2H ) F X T ε ) ]} ε 4H dζ. Choose F such that F ) + Λ x0 ) has a non-degenerate minimum at z. This implies that k F φ k T x 0, y 0 )) k 2 H H has a non-degenerate minimum at k 0 H H. For instance F z) = λ z x 2 [Λ x0,y 0 z) Λ x0,y 0 x)] with λ > 0).

33 Proof of Theorem 1 Replace ε 2H db B := W, W H )) in the SDE by ε 2H dw + k 0. Call the corresponding Girsanov-transformed process Z t ε = X t ε, Ỹ t ε): d X t ε = ε 2H+1 1 Ỹt 2 dt + 2 Ỹ t ε ε 2H dw t + k 0 ) 1 ), Girsanov factor dỹt = ε2h dw H t + k 0 ) 2. G = exp 1 T ε 2H ψk 0 ) tdb t 1 ) 0 2ε 4H k 0 2 H H. Therefore f x, T )e F x)/4ε4h = 1 2πε 2H R = 1 2πε 2H R where [ E e ε2h iζ X T x) ε 4H F X ] T ) G dζ [ E e )] dζ T ) = ε 2H iζ X T x) ε 4H F X T ) ε 2H ψγ) tdb t ε 4H γ 2 H H.

34 Proof of Theorem 1 Replace ε 2H db B := W, W H )) in the SDE by ε 2H dw + k 0. Call the corresponding Girsanov-transformed process Z t ε = X t ε, Ỹ t ε): d X t ε = ε 2H+1 1 Ỹt 2 dt + 2 Ỹ t ε ε 2H dw t + k 0 ) 1 ), Girsanov factor dỹt = ε2h dw H t + k 0 ) 2. G = exp 1 T ε 2H ψk 0 ) tdb t 1 ) 0 2ε 4H k 0 2 H H. By a stochastic Taylor expansion of Z t ε ) exp F X ε t ε 4H = exp = X ε t, Ỹ ε t ) for ε2h 0, [ 1 T )] F ε 4H x) ε 2H ψk 0 ) tdb t ε 2H XT Λ x 0 x) + Oε 4H ) 0

35 Proof of Theorem 1 Replace ε 2H db B := W, W H )) in the SDE by ε 2H dw + k 0. Call the corresponding Girsanov-transformed process Z t ε = X t ε, Ỹ t ε): d X t ε = ε 2H+1 1 Ỹt 2 dt + 2 Ỹ t ε ε 2H dw t + k 0 ) 1 ), Girsanov factor dỹt = ε2h dw H t + k 0 ) 2. G = exp 1 T ε 2H ψk 0 ) tdb t 1 ) 0 2ε 4H k 0 2 H H. By a stochastic Taylor expansion of Z t ε ) exp F X ε t ε 4H = exp = X ε t, Ỹ ε t ) for ε2h 0, [ F x) ε 4H + 1 T ε 2H ψk 0 ) tdb t + 1 ] X 0 ε 2H T Λ x 0 x) + O1).

36 Proof of Theorem 1 Replace ε 2H db B := W, W H )) in the SDE by ε 2H dw + k 0. Call the corresponding Girsanov-transformed process Z t ε = X t ε, Ỹ t ε): where d X t ε = ε 2H+1 1 Ỹt 2 dt + 2 Ỹ t ε ε 2H dw t + k 0 ) 1 ), dỹt = ε2h dw H t + k 0 ) 2. G = exp 1 T ε 2H ψk 0 ) tdb t 1 ) 0 2ε 4H k 0 2 H H. f x, T )e F x)/4ε4h = 1 [ 2πε 2H E gz ε )e ε 4H F X ] T ) G dζ R = 1 [ 2πε 2H E gz ε )e ε2h R e 1 ε 4H F x)+ 1 ] 2 k 0 HH ) 1 exp ε 2H X T Λ x x) 0 dζ R ) = ε 4H F x) ε 2H T T ε 2H ψγ) tdb t ε 4H γ 2 H H. 0 ) ψk 0 ) tdb t ε 2H XT Λ x 0 x)

37 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics 1 Introduction 2 3 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics

38 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics Direct call price expansion: For x 0 the option price satisfies cε 1 2H x, t) : = E [expx t) expε 1 2H x)) +] = E [expz t) expε 1 2H x)) + G ] I x)g 1 e ε 2H where Z is the controlled process around the optimal path k and G I x) = e ε 4H is the Girsanov factor for the optimal path, and g 1 a Gaussian random variable. Then cε 1 2H x, t) = exp [ Jε, x) := E exp I x) ε 4H ) exp xε 1 2H) Jε, x), where Û := Ẑ 1 ε x and I x) ε 2H ) Û expε 1 2H Û) 1) exp ] I ) x)r 2 1Û 0. Expansion uniformly in x) via epansion of the energy I directly Osajima).

39 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics Moderate regimes Moderate Regimes in the sense of Friz-Gerhold-Pinter ) 16) interpolate between out-of-the-money calls with fixed strike log K = k > 0 and at-the-money S0 k = 0 calls: Now k t = ct θ MOTM for 0 < θ < 1 ) and AATM for larger θ) 2 Reflects market reality: options closer expiry strikes closer to the money first observed by Mijatović-Tankov on FX markets The moderate regime MOTM) permits explicit computations for the rate function Λk) in terms of the model parameters Moderate deviations Advantage over large deviations OTM) case where the Λk) often related to geodesic distance problems and not explicitly available. MOTM expansions naturally involve quantities very familiar to practitioners, notably spot implied) volatility, implied volatility skew... In some cases fractional volatility models) the scaling θ permits a fine-tuning to understand the behavior and derivatives of the energy function.

40 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics Moderate regimes for rough volatility Rescalings = We tacitly agreed to consider P X t t 1/2 H x ). Now it is only a small step to consider instead for some suitable small θ > 0) ) P X t t 1/2 H+θ x.

41 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics Moderate regimes for rough volatility Rescalings = We tacitly agreed to consider P X t t 1/2 H x ). Now it is only a small step to consider instead for some suitable small θ > 0) ) P X t t 1/2 H+θ x. Theorem Bayer-Friz-Gulisashvili-H-Stemper) Consider a moderately out-of-the-money call k t = xt 1/2 H+θ ; θ 0, H) resp. θ 0, 2H ). Then as t 0, the following holds 3 log ck t, t) 1 x 2 2 Λ 0) t 2H 2θ + 1 x 3 6 Λ 0) t 2H 3θ, where we have explicit expressions: Λ 0) = 1 σ 2 0 Here K denotes the Volterra kernel and K, 1 := 1 0 and Λ 0) = ρ 6σ 0 σ0 4 K, 1. t 0 Kt, s)dsdt.

42 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics Moderate regimes for rough volatility Corollary MOTM Implied volatility skew) In the moderately out-of-the-money case k t = xt 1/2 H+θ ; θ 0, H); θ 0, 2H )) the 3 implied volatility satisfies the expansion σ impl k t, t) = σ 0 ρ σ 0 σ 0 1 t 0 0 Kt, s) dsdt k t t H 1/2 ) 1 + o1).

43 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics Moderate regimes for rough volatility Corollary MOTM Implied volatility skew) In the moderately out-of-the-money case k t = xt 1/2 H+θ ; θ 0, H); θ 0, 2H )) the 3 implied volatility satisfies the expansion σ impl k t, t) = σ 0 ρ σ 0 σ 0 1 t 0 0 Kt, s) dsdt k t t H 1/2 ) 1 + o1). Proof: The statement follows from the Theorem via matching components in the asymptotic expansions and by see Gao-Lee 2014)) using tσimpl 2 kt 2 kt, t) 2 log ck t, t).

44 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics Moderate regimes for rough volatility Corollary MOTM Implied volatility skew) In the moderately out-of-the-money case k t = xt 1/2 H+θ ; θ 0, H); θ 0, 2H )) the 3 implied volatility satisfies the expansion σ impl k t, t) = σ 0 ρ σ 0 σ 0 1 t 0 0 Kt, s) dsdt k t t H 1/2 ) 1 + o1). Proof: The statement follows from the Theorem via matching components in the asymptotic expansions and by see Gao-Lee 2014)) using tσimpl 2 kt 2 kt, t) 2 log ck t, t). This formula for the skew is in accord with ones previously derived by Alòs-León-Vives 2007) and Fukasawa 2011, 2016) in different settings.

45 Refined expansions and moderate regimes A non-markovian extention of Osajima s energy expansion Implied volatility asymptotics Thank you!

Asymptotic expansions for fractional stochastic volatility models

Asymptotic expansions for fractional stochastic volatility models Inroducion Main resul and moivaion Corollaries and oulook Asympoic expansions for fracional sochasic volailiy models Blanka Horvah Deparmen of Mahemaics, Imperial College London 3rd London-Paris Bachelier

More information

On density asymptotics for diffusions, old and new

On density asymptotics for diffusions, old and new On density asymptotics for diffusions, old and new Stefano De Marco (Ecole Polytechnique) Paris, January 2014 S. De Marco (Ecole Polytechnique) 7/01/14 0 / 15 d dx = b(x)dt + σ j (X)dW j, X 0 = x 0, t

More information

Malliavin Calculus in Finance

Malliavin Calculus in Finance Malliavin Calculus in Finance Peter K. Friz 1 Greeks and the logarithmic derivative trick Model an underlying assent by a Markov process with values in R m with dynamics described by the SDE dx t = b(x

More information

Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility

Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,

More information

Stochastic Volatility and Correction to the Heat Equation

Stochastic Volatility and Correction to the Heat Equation Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century

More information

The Uses of Differential Geometry in Finance

The Uses of Differential Geometry in Finance The Uses of Differential Geometry in Finance Andrew Lesniewski Bloomberg, November 21 2005 The Uses of Differential Geometry in Finance p. 1 Overview Joint with P. Hagan and D. Woodward Motivation: Varadhan

More information

Rough paths methods 1: Introduction

Rough paths methods 1: Introduction Rough paths methods 1: Introduction Samy Tindel Purdue University University of Aarhus 216 Samy T. (Purdue) Rough Paths 1 Aarhus 216 1 / 16 Outline 1 Motivations for rough paths techniques 2 Summary of

More information

Introduction to numerical simulations for Stochastic ODEs

Introduction to numerical simulations for Stochastic ODEs Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical

More information

A Class of Fractional Stochastic Differential Equations

A Class of Fractional Stochastic Differential Equations Vietnam Journal of Mathematics 36:38) 71 79 Vietnam Journal of MATHEMATICS VAST 8 A Class of Fractional Stochastic Differential Equations Nguyen Tien Dung Department of Mathematics, Vietnam National University,

More information

Dyson series for the PDEs arising in Mathematical Finance I

Dyson series for the PDEs arising in Mathematical Finance I for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by

More information

LAN property for sde s with additive fractional noise and continuous time observation

LAN property for sde s with additive fractional noise and continuous time observation LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,

More information

2012 NCTS Workshop on Dynamical Systems

2012 NCTS Workshop on Dynamical Systems Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,

More information

Harnack Inequalities and Applications for Stochastic Equations

Harnack Inequalities and Applications for Stochastic Equations p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline

More information

Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models

Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models Jean-Pierre Fouque North Carolina State University SAMSI November 3, 5 1 References: Variance Reduction for Monte

More information

Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises

Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,

More information

A new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009

A new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009 A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance

More information

Bridging the Gap between Center and Tail for Multiscale Processes

Bridging the Gap between Center and Tail for Multiscale Processes Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations

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

Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations.

Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint work with Michael Salins 2, Konstantinos Spiliopoulos 3.) 1. Department of Mathematics

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

Accurate approximation of stochastic differential equations

Accurate approximation of stochastic differential equations Accurate approximation of stochastic differential equations Simon J.A. Malham and Anke Wiese (Heriot Watt University, Edinburgh) Birmingham: 6th February 29 Stochastic differential equations dy t = V (y

More information

Gillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde

Gillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation

More information

Metastability for the Ginzburg Landau equation with space time white noise

Metastability for the Ginzburg Landau equation with space time white noise Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz Metastability for the Ginzburg Landau equation with space time white noise Barbara Gentz University of Bielefeld, Germany

More information

ON THE PROBABILITY OF HITTING THE BOUNDARY OF A BROWNIAN MOTION ON THE SABR PLANE

ON THE PROBABILITY OF HITTING THE BOUNDARY OF A BROWNIAN MOTION ON THE SABR PLANE ON THE PROBABILITY OF HITTING THE BOUNDARY OF A BROWNIAN MOTION ON THE SABR PLANE ARCHIL GULISASHVILI, BLANKA HORVATH, AND ANTOINE JACQUIER Abstract. Starting from the hyperbolic Brownian motion as a time-changed

More information

Geometric projection of stochastic differential equations

Geometric projection of stochastic differential equations Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically

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

Mean-field SDE driven by a fractional BM. A related stochastic control problem

Mean-field SDE driven by a fractional BM. A related stochastic control problem Mean-field SDE driven by a fractional BM. A related stochastic control problem Rainer Buckdahn, Université de Bretagne Occidentale, Brest Durham Symposium on Stochastic Analysis, July 1th to July 2th,

More information

Brownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion

Brownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of

More information

Lecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University

Lecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University Lecture 4: Ito s Stochastic Calculus and SDE Seung Yeal Ha Dept of Mathematical Sciences Seoul National University 1 Preliminaries What is Calculus? Integral, Differentiation. Differentiation 2 Integral

More information

A MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT

A MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields

More information

Convoluted Brownian motions: a class of remarkable Gaussian processes

Convoluted Brownian motions: a class of remarkable Gaussian processes Convoluted Brownian motions: a class of remarkable Gaussian processes Sylvie Roelly Random models with applications in the natural sciences Bogotá, December 11-15, 217 S. Roelly (Universität Potsdam) 1

More information

Ernesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010

Ernesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010 Optimal stopping for Hunt and Lévy processes Ernesto Mordecki 1 Lecture III. PASI - Guanajuato - June 2010 1Joint work with Paavo Salminen (Åbo, Finland) 1 Plan of the talk 1. Motivation: from Finance

More information

Paracontrolled KPZ equation

Paracontrolled KPZ equation Paracontrolled KPZ equation Nicolas Perkowski Humboldt Universität zu Berlin November 6th, 2015 Eighth Workshop on RDS Bielefeld Joint work with Massimiliano Gubinelli Nicolas Perkowski Paracontrolled

More information

Strong Solutions and a Bismut-Elworthy-Li Formula for Mean-F. Formula for Mean-Field SDE s with Irregular Drift

Strong Solutions and a Bismut-Elworthy-Li Formula for Mean-F. Formula for Mean-Field SDE s with Irregular Drift Strong Solutions and a Bismut-Elworthy-Li Formula for Mean-Field SDE s with Irregular Drift Thilo Meyer-Brandis University of Munich joint with M. Bauer, University of Munich Conference for the 1th Anniversary

More information

Stochastic Differential Equations

Stochastic Differential Equations Chapter 5 Stochastic Differential Equations We would like to introduce stochastic ODE s without going first through the machinery of stochastic integrals. 5.1 Itô Integrals and Itô Differential Equations

More information

arxiv: v1 [q-fin.mf] 8 Jan 2017

arxiv: v1 [q-fin.mf] 8 Jan 2017 FUNCTIONAL ANALYTIC (IR-)REGULARITY PROPERTIES OF SABR-TYPE PROCESSES LEIF DÖRING, BLANKA HORVATH, AND JOSEF TEICHMANN arxiv:1701.02015v1 [q-fin.mf] 8 Jan 2017 Abstract. The SABR model is a benchmark stochastic

More information

Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models

Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models Scott Robertson October 19, 29 Abstract Sample path Large Deviation Principles LDP) of the Freidlin-Wentzell

More information

The concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.

The concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt. The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes

More information

Fast-slow systems with chaotic noise

Fast-slow systems with chaotic noise Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David

More information

Approximating diffusions by piecewise constant parameters

Approximating diffusions by piecewise constant parameters Approximating diffusions by piecewise constant parameters Lothar Breuer Institute of Mathematics Statistics, University of Kent, Canterbury CT2 7NF, UK Abstract We approximate the resolvent of a one-dimensional

More information

Affine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009

Affine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009 Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine

More information

Numerical Integration of SDEs: A Short Tutorial

Numerical Integration of SDEs: A Short Tutorial Numerical Integration of SDEs: A Short Tutorial Thomas Schaffter January 19, 010 1 Introduction 1.1 Itô and Stratonovich SDEs 1-dimensional stochastic differentiable equation (SDE) is given by [6, 7] dx

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

Stochastic Integration and Stochastic Differential Equations: a gentle introduction

Stochastic Integration and Stochastic Differential Equations: a gentle introduction Stochastic Integration and Stochastic Differential Equations: a gentle introduction Oleg Makhnin New Mexico Tech Dept. of Mathematics October 26, 27 Intro: why Stochastic? Brownian Motion/ Wiener process

More information

Convergence to equilibrium for rough differential equations

Convergence to equilibrium for rough differential equations Convergence to equilibrium for rough differential equations Samy Tindel Purdue University Barcelona GSE Summer Forum 2017 Joint work with Aurélien Deya (Nancy) and Fabien Panloup (Angers) Samy T. (Purdue)

More information

Stochastic optimal control with rough paths

Stochastic optimal control with rough paths Stochastic optimal control with rough paths Paul Gassiat TU Berlin Stochastic processes and their statistics in Finance, Okinawa, October 28, 2013 Joint work with Joscha Diehl and Peter Friz Introduction

More information

Nonparametric Drift Estimation for Stochastic Differential Equations

Nonparametric Drift Estimation for Stochastic Differential Equations Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,

More information

Lecture 4: Introduction to stochastic processes and stochastic calculus

Lecture 4: Introduction to stochastic processes and stochastic calculus Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London

More information

Dynamic Bertrand and Cournot Competition

Dynamic Bertrand and Cournot Competition Dynamic Bertrand and Cournot Competition Effects of Product Differentiation Andrew Ledvina Department of Operations Research and Financial Engineering Princeton University Joint with Ronnie Sircar Princeton-Lausanne

More information

Hypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations.

Hypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations. Hypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint works with Michael Salins 2 and Konstantinos Spiliopoulos 3.) 1. Department of Mathematics

More information

HJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011

HJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011 Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance

More information

Orthonormal polynomial expansions and lognormal sum densities

Orthonormal polynomial expansions and lognormal sum densities 1/28 Orthonormal polynomial expansions and lognormal sum densities Pierre-Olivier Goffard Université Libre de Bruxelles pierre-olivier.goffard@ulb.ac.be February 22, 2016 2/28 Introduction Motivations

More information

1 Brownian Local Time

1 Brownian Local Time 1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =

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

Effective dynamics for the (overdamped) Langevin equation

Effective dynamics for the (overdamped) Langevin equation Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath

More information

Stochastic contraction BACS Workshop Chamonix, January 14, 2008

Stochastic contraction BACS Workshop Chamonix, January 14, 2008 Stochastic contraction BACS Workshop Chamonix, January 14, 2008 Q.-C. Pham N. Tabareau J.-J. Slotine Q.-C. Pham, N. Tabareau, J.-J. Slotine () Stochastic contraction 1 / 19 Why stochastic contraction?

More information

Approximation of BSDEs using least-squares regression and Malliavin weights

Approximation of BSDEs using least-squares regression and Malliavin weights Approximation of BSDEs using least-squares regression and Malliavin weights Plamen Turkedjiev (turkedji@math.hu-berlin.de) 3rd July, 2012 Joint work with Prof. Emmanuel Gobet (E cole Polytechnique) Plamen

More information

Multi-Factor Finite Differences

Multi-Factor Finite Differences February 17, 2017 Aims and outline Finite differences for more than one direction The θ-method, explicit, implicit, Crank-Nicolson Iterative solution of discretised equations Alternating directions implicit

More information

Integration by parts and Quasi-invariance of Horizontal Wiener Measure on a foliated compact manifold

Integration by parts and Quasi-invariance of Horizontal Wiener Measure on a foliated compact manifold Integration by parts and Quasi-invariance of Horizontal Wiener Measure on a foliated compact manifold Qi Feng (Based on joint works with Fabrice Baudoin and Maria Gordina) University of Connecticut Mathematical

More information

Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model

Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model Xiaowei Chen International Business School, Nankai University, Tianjin 371, China School of Finance, Nankai

More information

Obstacle problems for nonlocal operators

Obstacle problems for nonlocal operators Obstacle problems for nonlocal operators Camelia Pop School of Mathematics, University of Minnesota Fractional PDEs: Theory, Algorithms and Applications ICERM June 19, 2018 Outline Motivation Optimal regularity

More information

Perturbed Gaussian Copula

Perturbed Gaussian Copula Perturbed Gaussian Copula Jean-Pierre Fouque Xianwen Zhou August 8, 6 Abstract Gaussian copula is by far the most popular copula used in the financial industry in default dependency modeling. However,

More information

Small ball probabilities and metric entropy

Small ball probabilities and metric entropy Small ball probabilities and metric entropy Frank Aurzada, TU Berlin Sydney, February 2012 MCQMC Outline 1 Small ball probabilities vs. metric entropy 2 Connection to other questions 3 Recent results for

More information

Large Deviations for Small-Noise Stochastic Differential Equations

Large Deviations for Small-Noise Stochastic Differential Equations Chapter 21 Large Deviations for Small-Noise Stochastic Differential Equations This lecture is at once the end of our main consideration of diffusions and stochastic calculus, and a first taste of large

More information

Particle models for Wasserstein type diffusion

Particle models for Wasserstein type diffusion Particle models for Wasserstein type diffusion Vitalii Konarovskyi University of Leipzig Bonn, 216 Joint work with Max von Renesse konarovskyi@gmail.com Wasserstein diffusion and Varadhan s formula (von

More information

SMALL-TIME EXPANSIONS FOR LOCAL JUMP-DIFFUSIONS MODELS WITH INFINITE JUMP ACTIVITY. 1. Introduction

SMALL-TIME EXPANSIONS FOR LOCAL JUMP-DIFFUSIONS MODELS WITH INFINITE JUMP ACTIVITY. 1. Introduction SMALL-TIME EXPANSIONS FOR LOCAL JUMP-DIFFUSIONS MODELS WITH INFINITE JUMP ACTIVITY JOSÉ E. FIGUEROA-LÓPEZ AND CHENG OUYANG Abstract. We consider a Markov process {X (x) t } t with initial condition X (x)

More information

Vast Volatility Matrix Estimation for High Frequency Data

Vast Volatility Matrix Estimation for High Frequency Data Vast Volatility Matrix Estimation for High Frequency Data Yazhen Wang National Science Foundation Yale Workshop, May 14-17, 2009 Disclaimer: My opinion, not the views of NSF Y. Wang (at NSF) 1 / 36 Outline

More information

ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER

ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose

More information

Utility Maximization in Hidden Regime-Switching Markets with Default Risk

Utility Maximization in Hidden Regime-Switching Markets with Default Risk Utility Maximization in Hidden Regime-Switching Markets with Default Risk José E. Figueroa-López Department of Mathematics and Statistics Washington University in St. Louis figueroa-lopez@wustl.edu pages.wustl.edu/figueroa

More information

Pathwise volatility in a long-memory pricing model: estimation and asymptotic behavior

Pathwise volatility in a long-memory pricing model: estimation and asymptotic behavior Pathwise volatility in a long-memory pricing model: estimation and asymptotic behavior Ehsan Azmoodeh University of Vaasa Finland 7th General AMaMeF and Swissquote Conference September 7 1, 215 Outline

More information

Optimal exit strategies for investment projects. 7th AMaMeF and Swissquote Conference

Optimal exit strategies for investment projects. 7th AMaMeF and Swissquote Conference Optimal exit strategies for investment projects Simone Scotti Université Paris Diderot Laboratoire de Probabilité et Modèles Aléatories Joint work with : Etienne Chevalier, Université d Evry Vathana Ly

More information

Thomas Knispel Leibniz Universität Hannover

Thomas Knispel Leibniz Universität Hannover Optimal long term investment under model ambiguity Optimal long term investment under model ambiguity homas Knispel Leibniz Universität Hannover knispel@stochastik.uni-hannover.de AnStAp0 Vienna, July

More information

SDE Coefficients. March 4, 2008

SDE Coefficients. March 4, 2008 SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)

More information

Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias

Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias José E. Figueroa-López 1 1 Department of Statistics Purdue University Seoul National University & Ajou University

More information

Theory and applications of random Poincaré maps

Theory and applications of random Poincaré maps Nils Berglund nils.berglund@univ-orleans.fr http://www.univ-orleans.fr/mapmo/membres/berglund/ Random Dynamical Systems Theory and applications of random Poincaré maps Nils Berglund MAPMO, Université d'orléans,

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

A class of globally solvable systems of BSDE and applications

A class of globally solvable systems of BSDE and applications A class of globally solvable systems of BSDE and applications Gordan Žitković Department of Mathematics The University of Texas at Austin Thera Stochastics - Wednesday, May 31st, 2017 joint work with Hao

More information

Some Tools From Stochastic Analysis

Some Tools From Stochastic Analysis W H I T E Some Tools From Stochastic Analysis J. Potthoff Lehrstuhl für Mathematik V Universität Mannheim email: potthoff@math.uni-mannheim.de url: http://ls5.math.uni-mannheim.de To close the file, click

More information

Probability Distribution in the SABR Model of Stochastic Volatility

Probability Distribution in the SABR Model of Stochastic Volatility Probability Distribution in the SABR Model of Stochastic Volatility Patrick Hagan Bloomberg LP 499 Park Avenue New York, NY 10022 Andrew Lesniewski Ellington Management Group 53 Forest Avenue Old Greenwich,

More information

Optimal portfolio strategies under partial information with expert opinions

Optimal portfolio strategies under partial information with expert opinions 1 / 35 Optimal portfolio strategies under partial information with expert opinions Ralf Wunderlich Brandenburg University of Technology Cottbus, Germany Joint work with Rüdiger Frey Research Seminar WU

More information

Precise Asymptotics for Infinite Dimensional Itô-Lyons Maps of Brownian Rough Paths

Precise Asymptotics for Infinite Dimensional Itô-Lyons Maps of Brownian Rough Paths Precise Asymptotics for Infinite Dimensional Itô-Lyons Maps of Brownian Rough Paths Jointwork with Yuzuru INAHAMA (Nagoya Univ.) Hiroshi KAWABI (Department of Mathematics, Faculty of Science Okayama University,

More information

Modelling Interactions Between Weather and Climate

Modelling Interactions Between Weather and Climate Modelling Interactions Between Weather and Climate p. 1/33 Modelling Interactions Between Weather and Climate Adam Monahan & Joel Culina monahana@uvic.ca & culinaj@uvic.ca School of Earth and Ocean Sciences,

More information

Bernstein-gamma functions and exponential functionals of Lévy processes

Bernstein-gamma functions and exponential functionals of Lévy processes Bernstein-gamma functions and exponential functionals of Lévy processes M. Savov 1 joint work with P. Patie 2 FCPNLO 216, Bilbao November 216 1 Marie Sklodowska Curie Individual Fellowship at IMI, BAS,

More information

Introduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0

Introduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0 Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods

More information

Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process

Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process Asymptotic properties of maximum likelihood estimator for the growth rate for a jump-type CIR process Mátyás Barczy University of Debrecen, Hungary The 3 rd Workshop on Branching Processes and Related

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

Skew Brownian Motion and Applications in Fluid Dispersion

Skew Brownian Motion and Applications in Fluid Dispersion Skew Brownian Motion and Applications in Fluid Dispersion Ed Waymire Department of Mathematics Oregon State University Corvallis, OR 97331 *Based on joint work with Thilanka Appuhamillage, Vrushali Bokil,

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

An adaptive numerical scheme for fractional differential equations with explosions

An adaptive numerical scheme for fractional differential equations with explosions An adaptive numerical scheme for fractional differential equations with explosions Johanna Garzón Departamento de Matemáticas, Universidad Nacional de Colombia Seminario de procesos estocásticos Jointly

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart

More information

Generalised Fractional-Black-Scholes Equation: pricing and hedging

Generalised Fractional-Black-Scholes Equation: pricing and hedging Generalised Fractional-Black-Scholes Equation: pricing and hedging Álvaro Cartea Birkbeck College, University of London April 2004 Outline Lévy processes Fractional calculus Fractional-Black-Scholes 1

More information

ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME

ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME SAUL D. JACKA AND ALEKSANDAR MIJATOVIĆ Abstract. We develop a general approach to the Policy Improvement Algorithm (PIA) for stochastic control problems

More information

Pseudo-marginal MCMC methods for inference in latent variable models

Pseudo-marginal MCMC methods for inference in latent variable models Pseudo-marginal MCMC methods for inference in latent variable models Arnaud Doucet Department of Statistics, Oxford University Joint work with George Deligiannidis (Oxford) & Mike Pitt (Kings) MCQMC, 19/08/2016

More information

Stochastic Calculus. Kevin Sinclair. August 2, 2016

Stochastic Calculus. Kevin Sinclair. August 2, 2016 Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed

More information

Large Deviations for Small-Noise Stochastic Differential Equations

Large Deviations for Small-Noise Stochastic Differential Equations Chapter 22 Large Deviations for Small-Noise Stochastic Differential Equations This lecture is at once the end of our main consideration of diffusions and stochastic calculus, and a first taste of large

More information

Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches

Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes

More information

Problems 5: Continuous Markov process and the diffusion equation

Problems 5: Continuous Markov process and the diffusion equation Problems 5: Continuous Markov process and the diffusion equation Roman Belavkin Middlesex University Question Give a definition of Markov stochastic process. What is a continuous Markov process? Answer:

More information

DUBLIN CITY UNIVERSITY

DUBLIN CITY UNIVERSITY DUBLIN CITY UNIVERSITY SAMPLE EXAMINATIONS 2017/2018 MODULE: QUALIFICATIONS: Simulation for Finance MS455 B.Sc. Actuarial Mathematics ACM B.Sc. Financial Mathematics FIM YEAR OF STUDY: 4 EXAMINERS: Mr

More information

Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise

Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Robert J. Elliott 1 Samuel N. Cohen 2 1 Department of Commerce, University of South Australia 2 Mathematical Insitute, University

More information

Ambiguity and Information Processing in a Model of Intermediary Asset Pricing

Ambiguity and Information Processing in a Model of Intermediary Asset Pricing Ambiguity and Information Processing in a Model of Intermediary Asset Pricing Leyla Jianyu Han 1 Kenneth Kasa 2 Yulei Luo 1 1 The University of Hong Kong 2 Simon Fraser University December 15, 218 1 /

More information