arxiv: v1 [q-fin.mf] 19 Jan 2019

Transcription

1 Opton Prcng n Illqud Markets wth Jumps arxv: v1 [q-fn.mf 19 Jan 219 José M. T. S. Cruz a and Danel Ševčovčb a ISEG, Unversty of Lsbon, ua de Quelhas 6, Lsbon, Portugal; b Comenus Unversty n Bratslava, Mlynská dolna, Bratslava, Slovaka ATICLE HISTOY Compled January 23, 219 ABSTACT The classcal lnear Black Scholes model for prcng dervatve securtes s a popular model n fnancal ndustry. It reles on several restrctve assumptons such as completeness, and frctonless of the market as well as the assumpton on the underlyng asset prce dynamcs followng a geometrc Brownan moton. The man purpose of ths paper s to generalze the classcal Black Scholes model for prcng dervatve securtes by takng nto account feedback effects due to an nfluence of a large trader on the underlyng asset prce dynamcs exhbtng random jumps. The assumpton that an nvestor can trade large amounts of assets wthout affectng the underlyng asset prce tself s usually not satsfed, especally n llqud markets. We generalze the Frey Stremme nonlnear opton prcng model for the case the underlyng asset follows a Lévy stochastc process wth jumps. We derve and analyze a fully nonlnear parabolc partal-ntegro dfferental equaton for the prce of the opton contract. We propose a sem-mplct numercal dscretzaton scheme and perform varous numercal experments showng nfluence of a large trader and ntensty of jumps on the opton prce. KEYWODS Nonlnear partal ntegro-dfferental equaton; Lévy measure; Fnte dfference approxmaton 1. Introducton Over recent decades, the Black Scholes model and ts generalzatons become wdely used n fnancal markets because of ts smplcty and exstence of the analytc formula for prcng European style optons. Accordng to the classcal theory developed by Black, Scholes and Merton the prce V (t, S of an opton n a stylzed fnancal market at tme t [, T and the underlyng asset prce S can be computed as a soluton to the lnear Black Scholes parabolc equaton: t (t, S σ2 S 2 2 V (t, S + rs (t, S rv (t, S =, t [, T, S >. (1 2 Here σ > s the hstorcal volatlty of the underlyng asset drven by the geometrc Brownan moton, r > s the rsk-free nterest rate of zero-coupon bond. A soluton s subject to the termnal pay-off condton V (T, S = Φ(S at maturty t = T. Evdence CONTACT J. Cruz. Emal: jmcruz@seg.ulsboa.pt to appear n: Appled Mathematcal Fnance

2 from stock markets observatons ndcates that ths model s not the most realstc one, snce t assumes that the market s lqud, complete, frctonless and wthout transacton costs. We also recall that the lnear Black Scholes equaton provdes a soluton correspondng to a perfectly replcated portfolo whch need not be a desrable property. In the last two decades some of these assumptons have been relaxed n order to model, for nstance, the presence of transacton costs (see e.g. Kwok [1 and Avellaneda and Paras [2, feedback and llqud market effects due to large traders choosng gven stock-tradng strateges (Schönbucher and Wllmott [13, Frey and Pate [7, Frey and Stremme [6, rsk from the unprotected portfolo (Jandačka and Ševčovč [8. In all aforementoned generalzatons of the lnear Black Scholes equaton (1 the constant volatlty σ s replaced by a nonlnear functon σ(s S 2 V dependng on the second dervatve 2 V of the opton prce tself. In the class of generalzed Black Scholes equaton wth such a nonlnear dffuson functon, an mportant role s played by the nonlnear Black Scholes model derved by Frey and Stremme n [8 (see also [7,[5. In ths model the asset dynamcs takes nto account the presence of feedback effects due to a large trader choosng hs/her stock-tradng strategy (see also [13. The dffuson coeffcent s agan non-constant: σ(s 2 SV 2 = σ 2 ( 1 ϱs 2 SV 2, (2 where σ, ϱ > are constants. Another mportant drecton n generalzng the orgnal Black Scholes equaton arses from the fact that the sample paths of a Brownan moton are contnuous, but the realzed stock prce of a typcal company exhbts random jumps over the ntraday scale, makng the prce trajectores dscontnuous. In the classcal Black Scholes model the underlyng asset prce process s assumed to follow a geometrc Brownan moton. However the emprcal dstrbuton of stock returns exhbts fat tals. Several alternatves have been proposed n the lterature for the generalzaton of ths model. The models wth jumps and dffuson can, at least n part, solve the problems nherent to the lnear Black Scholes model and they have also an mportant role n the optons market. Whle n the Black Scholes model the market s complete, mplyng that every pay-off can be perfectly replcated, n jump dffuson models there s no perfect hedge and ths way the optons are not redundant. It turns out that the opton prce can be computed from the soluton V (t, S of the followng partal ntegrodfferental (PIDE Black Scholes equaton: t (t, S σ2 S 2 2 V (t, S + rs (t, S rv (t, S 2 + V (t, S + H(z, S V (t, S H(z, S (t, Sν(dz =, (3 where H(z, S = S(e z 1 and ν s the so-called Lévy measure characterzng the underlyng asset process wth random jumps n tme and space. Note that, f ν = then (3 reduces to the classcal lnear Black Scholes equaton (1. The novelty and man purpose of ths paper s to take nto account both drectons of generalzatons of the Black Scholes equaton. The assumpton that an nvestor can trade large amounts of the underlyng asset wthout affectng ts prce s no longer true, especally n llqud markets. Therefore, we wll derve, analyze, and perform numercal computaton of the model. We relax the assumpton of lqud market followng the Frey Stremme model under the assumpton that the underlyng asset prce follows 2

3 a Lévy stochastc process wth jumps. We wll show that the correspondng PIDE nonlnear equaton has the form: t (1 ϱs S φ 2 S2 2 V + rs 2 rv + V (t, S + H(t, z, S V (t, S H(t, z, S ν(dz =, (4 σ 2 where the functon H(t, z, S may depend e.g. on the large trader strategy functon φ = φ(t, S. Ths functon may depend on the delta S V of the prce V, f ϱ >. The paper s organzed as follows. In the next Secton 2 we recall known facts regardng exponental Lévy models. We also recall mportant classes of Lévy measures wth fnte and nfnte actvty. Secton 3 s devoted to dervaton of the novel opton prcng model takng nto account feedback effects of a large trader on the underlyng asset followng jump-dffuson Lévy process. We show that the prce of an opton can be computed from a soluton to a fully nonlnear partal ntegro-dfferental equaton (PIDE (4. We also derve a formula for the tradng strategy functon φ whch mnmzes the varance of the trackng error. Next, n Secton 4 we present a sem-mplct numercal dscretzaton scheme for solvng the resultng nonlnear PIDE. The scheme s based on fnte dfference approxmaton. In Secton 5 we present numercal results when consderng the Varance Gamma process. We also perform senstvty numercal analyss of a soluton wth respect to model parameters ρ and ν. 2. Prelmnares, defntons and motvaton We consder a stylzed economy wth two traded assets, a rsky asset, usually a stock wth a prce S t, and a rsk-less asset, typcally a bond wth a prce B t whch s taken as numerare. The bond market s assumed to be perfectly elastc as bonds are assumed to be more lqud when compared to stocks. In ths economy there are two type of traders: the reference traders and the program traders. The program traders are also referred to as portfolo nsurers snce they use dynamc hedgng strateges to hedge portfolo aganst jumps n stock prces. They are ether sngle traders or a group of traders actng together. It s assumed that ther trades nfluence the stock prce equlbrum. The reference traders can be consdered as representatve traders of many small agents. We assume they act as prce takers. Typcally, t s assumed that D(t, Y t, S t s the reference trader demand functon whch depends on the ncome process Y t or some other fundamental state varable nfluencng the reference trader demand. The aggregate demand of program traders s denoted by ϕ(t, S t = ξφ(t, S t, where ξ s the number of wrtten dentcal securtes that the program traders are tryng to hedge and φ(t, S t s the demand per unt of the securty beng hedged. For smplcty, we assume that ξ s the same for every program trader. The general case where dfferent securtes are consdered can be found e.g. n [14. Assume the supply of a stock wth the prce S s constant. Let D(t, Y, S = D(t,Y,S denote the quantty S demanded by a reference trader per unt of supply. Then the total demand relatve to the supply at tme t s gven by G(t, Y, S = D(t, Y, S + ρφ(t, S, where ρ = ξ S and ρφ(t, S s the proporton of the total supply of the stock that s beng traded by program traders. In order to obtan the market equlbrum the varables Y and S should satsfy G(t, Y, S = 1. Assume that the functon G s monotone wth respect to 3

4 Y and S varables, and t s suffcently smooth. Then we can solve the mplct equaton G(t, Y t, S t = 1 to obtan S t = ψ(t, Y t where ψ s a suffcently smooth functon. Followng [14, we assume that the stochastc process Y t has the followng dynamcs: dy t = µ(t, Y t dt + η(t, Y t dw t. Then, by usng Itô s lemma for the process S t = ψ(t, Y t we obtan ds t = ( t ψ + µ y ψ + η2 2 2 yψ dt + η y ψdw t b(t, S t S t dt + v(t, S t S t dw t. (5 It means that S t follows a geometrc Brownan moton wth a nonconstant volatlty functon v(t, S = η(t, Y Y ψ(t, Y /ψ(t, Y where Y = ψ 1 (t, S. Followng the argument used n the dervaton of the orgnal Black Scholes equaton we obtan a generalzaton of the Black Scholes partal dfferental equaton wth a nonconstant volatlty functon σ = v(t, S. In ths paper we follow Frey and Stremme s approach (cf. [7, 5. The dea s to prescrbe a dynamcs for the underlyng stock prce nstead of dervng t by usng the market equlbrum and dynamcs for the ncome process Y t as t s done e.g. n [14. Ths way Frey and Stremme derved the same stock prce dynamcs as n [14 correspondng to a stuaton where the demand functon s of logarthmc type, D(Y, S = ln( Y γ S, where γ = σ η, and the ncome process Y t follows a geometrc Brownan moton,.e. Y D(Y, S = γ 1 Y, SD(Y, S = 1 S, dy t = µ Y t dt + η Y t dw t, (6 v(t, S = η(t, Y Y ψ(t, Y ψ(t, Y = η Y S γ 1 Y 1 S + ρ φ = σ 1 ρs φ. Assumng the delta hedgng strategy wth φ(t, S = S V (t, S and nsertng the volatlty functon v(t, S nto (5 we obtan the generalzed Black Scholes equaton wth the nonlnear dffuson functon of the form (2. Our man goal s to extend the Frey Stremme model to an underlyng asset followng a Lévy process. Next, we recall basc propertes of Lévy jump-dffuson processes Exponental Lévy models Let X t, t, be a stochastc process. The measure ν(a of a Borel set A B( defned by ν(a = E [J X ([, 1 A where J X ([, t A = # {s [, t : X s A} s the Posson random measure. It gves the mean number, per unt of tme, of jumps whose ampltude belongs to the set A. ecall that the Lévy-Itô decomposton provdes a representaton of X t whch can be nterpreted as a combnaton of a Brownan moton wth a drft ω and an nfnte sum of ndependent compensated Posson processes wth varable jump szes (see [4,.e. dx t = ωdt + σdw t + xj X (dt, dx + x JX (dt, dx, x 1 x <1 where JX ([, t A = J X ([, t A tν (A s the compensaton of J X. 4

5 Any Lévy process s a strong Markov process, the assocated semgroup s a convoluton semgroup. Its nfntesmal generator L : u L[u s a nonlocal partal ntegro-dfferental operator gven by (see [1: E [u (x + X h u (x L[u(x = lm h + h = σ2 2 u 2 x 2 + γ u [ x + u u (x + y u (x y1 y 1 x (x ν(dy, (7 whch s well defned for any compactly supported functon u C 2 (. Let S t, t, be a stochastc process representng an underlyng asset process under a fltered probablty space (Ω, F, {F t }, P. The fltraton {F t } represents the prce hstory up to the tme t. If the market s arbtrage-free then there s an equvalent measure Q under whch dscounted prces of all traded fnancal assets are Q martngales. Ths result s known as the fundamental theorem of asset prcng (see [4. The measure Q s also known as the rsk neutral measure. We consder the exponental Lévy model n whch the rsk-neutral prce process S t under Q s gven by S t = e rt+xt, where X t s a Lévy process under Q wth the characterstc trplet (σ, γ, ν. Then the arbtrage-free market hypothess mposes that Ŝ t = S t e rt = e Xt s a martngale, whch s equvalent to the followng condtons mposed on the trplet (σ, γ, ν: y 1 e y ν(dy <, γ, γ = σ2 2 + The rsk-neutral dynamcs of S t under Q s gven by ( e y 1 y1 y 1 ν(dy. (8 ds t = rs t dt + σs t dw t + (e y 1 S t JX (dt, dy. (9 The exponental prce process e Xt, t, s also a Markov process wth the state space (, and the nfntesmal generator: L S [V (S = Xh E[V (Se V (S lm h h + = rs + σ2 [ V (Se y V (S S(e y 1 2 S2 2 V 2 (1 ν(dy (11 (see [4. ecall that a Lévy process s called the Lévy type stochastc ntegral f t has the followng representaton: dx t = ωdt + σdw t + K(t, xj X (dt, dx + H(t, x J X (dt, dx. x 1 x <1 An mportant result that wll be needed later s the followng varant of Itô s lemma. Theorem 2.1. [1 Let f C 1,2 ([, T and H, K C([, T. Let X t, t, 5

6 be a Lévy stochastc process. Then df(t, X t = f f dt + t x dx t f 2 x 2 d[x t, X t + f(t, X t + K(t, x f(t, X t J X (dt, dx ( x 1 x <1 x <1 f(t, X t + H(t, x f(t, X t J X (dt, dx f(t, X t + H(t, x f(t, X t H(t, x f x (t, X tν(dxdt. A classcal example of a Lévy process s the jump-dffuson model ntroduced by Merton n [12. It has the followng dynamcs: ( dx t = b + x <1 xν(dx Then, by applyng Itô s lemma to S t = e Xt dt + xj X (dt, dx + σdw t + x J X (dt, dx. x 1 x <1 we obtan: ds t = (b + σ 2 /2S t dt + σs t dw t + S t (e x 1 J X (dt, dx. In fnancal applcatons, exponental Lévy models are of several types. In ths paper we are concerned wth the so-called jump-dffuson models n whch we represent the log-prce as a Lévy process wth a non-zero dffuson part (σ > and a jump process wth ether fnte actvty wth (ν( < or nfnte actvty (ν( =. In the context of fnancal modellng a jump-dffuson model was proposed by Merton n [12. The random jump varables are normally dstrbuted wth the mean m and varance δ 2. Its Lévy densty s gven by: ν(dx = λ 1 δ (x m 2 2π e 2δ 2 dx. (13 Another popular and frequently used model s the so-called double exponental model whch was ntroduced by Kou n [9. In ths model the dstrbuton of jumps have a Lévy measure of the form: ν(dx = λ (θλ + e λ+x 1 x> + (1 θλ e λ x 1 x< dx, (14 where λ s the ntensty of jumps, θ s the probablty of havng a postve jump and λ ± > correspond to the level of decay of the dstrbuton of postve and negatve jumps. Ths mples that the dstrbuton of jumps s asymmetrc and the tals of the dstrbuton of returns are sem-heavy. Among examples of nfnte actvty Lévy processes used n the fnancal modellng there are e.g. the Varance Gamma (see [11, Normal Inverse Gaussan (NIG (see [3, or CGMY processes. The Varance Gamma process s a process of nfnte actvty 6

7 ν( = and fnte varaton x 1 x ν(dx <. Its Lévy measure s gven by: ν(dx = 1 κ x eax B x dx wth A = θ/σ 2 and B = σ 2 θ 2 + 2σ 2 /κ, (15 where the parameters σ and θ are related to volatlty and drft of the Brownan moton and κ s a parameter related to the varance of the subordnator, n ths case the Gamma process (see [4. All aforementoned examples of jump-dffuson models have Lévy measures belongng to the class of the so-called admssble actvty Lévy measures. Defnton 2.2. A Lévy measure ν s called an admssble actvty Lévy measure f ν(dz dz h(z C z α ( e D z 1 z + e D+z 1 z< e µz2, (16 for any z and shape parameters α, D ± and µ. emark 1. Note that the addtonal condtons mn(z2, 1ν(dz < and z >1 ez ν(dz < are satsfed provded that ν s an admssble Lévy measure wth shape parameters α < 3, and, ether µ >, D ±, or µ = and D + 1 < < D +. For the Merton model we have α =, D ± = and µ = 1/(2δ 2 >. In the Kou model α = µ =, D + = λ, D = λ +. As for the Varance Gamma process we have α = 1, µ =, D ± = A ± B. 3. Feedback effects under jump-dffuson underlyng asset prce dynamcs Let us suppose that a large trader uses a stock-holdng strategy α t and S t s a cadlag process (rght contnuous wth lmts to the left. Henceforth, we shall dentfy S t wth S t. We assume S t has the followng dynamcs: ds t = µs t dt + σs t dw t + ρs t dα t + S t (e x 1J X (dt, dx. (17 It can be vewed as a perturbaton of the classcal jump-dffuson model. Indeed, f a large trader does not trade then α t = or the market lqudty parameter ρ s set to zero then the stock prce S t follows the classcal jump-dffuson model. In what follows, we wll assume the followng structural hypothess: Assumpton 3.1. Assume the tradng strategy α t = φ(t, S t and the parameter ρ satsfy ρl < 1, where L = sup S> S φ. Next we show an explct formula for the dynamcs of S t satsfyng (17 under certan regularty assumptons made on the stock-holdng functon φ(t, S. Proposton 3.2. Suppose that the stock-holdng strategy α t = φ(t, S t satsfes Assumpton 3.1 where φ C 1,2 ([, T +. If the process S t, t, satsfes the mplct stochastc equaton (17 then the process S t satsfes the followng SDE: ds t = b(t, S t S t dt + v(t, S t S t dw t + H(t, x, S t J X (dt, dx, (18 7

8 where 1 b(t, S = 1 ρs φ (t, S σ v(t, S = ( ( φ µ + ρ t v(t, S2 S 2 2 φ 2, (19 1 ρs φ (2 (t, S, H(t, x, S = S(e x 1 + ρs [φ(t, S + H(t, x, S φ(t, S. (21 Proof. We can rewrte the SDE (18 for S t, n the followng way: ds t = ( b(t, S t S t + x <1 H(t, x, S t ν(dx + H(t, x, S t J X (dt, dx + x 1 x <1 dt + v(t, S t S t dw t H(t, x, S t J X (dt, dx. Snce φ(t, S s assumed to be a smooth functon then, by applyng Itô formula (12 to the process φ(t, S t, we obtan dα t = ( φ t v(t, S t 2 St φ 2 dt + φ ds t (22 φ(t, S t + H(t, x, S t φ(t, S t H(t, x, S t φ (t, S tj X (dt, dx. Now, nsertng the dfferental dα t nto (17, we obtan ds t = µs t dt + σs t dw t + ( φ +ρs t t v(t, S t 2 St 2 +ρs t earrangng terms n (23 we conclude S t (e x φ 1J X (dt, dx + ρs t ds t 2 φ 2 dt (23 φ(t, S t + H(t, x, S t φ(t, S t H(t, x, S t φ (t, S tj X (dt, dx. φ (1 ρs t (t, S tds t = (µs t + ρs t ( φ t v(t, S t 2 St 2 2 φ 2 dt +σs t dw t ρs t H(t, x, S t φ (t, S tj X (dt, dx (24 + S t (e x 1 + ρs t (φ(t, S t + H(t, x, S t φ(t, S t J X (dt, dx. Comparng terms n (18 and (24 we end up wth expressons (19, (2, and the 8

9 mplct equaton for the functon H: H(t, x, S = 1 1 ρs φ (t, S (S(ex 1 + ρs (φ(t, S + H(t, x, S φ(t, S 1 1 ρs φ φ (t, SH(t, x, S. (25 (t, SρS Smplfyng ths expresson for H we conclude (21, as clamed. The functon H s gven mplctly by equaton (21. If we expand ts soluton H n terms of a small parameter ρ,.e. H(t, x, S = H (t, x, S + ρh 1 (t, x, S + O(ρ 2 as ρ, we conclude the followng proposton: Proposton 3.3. Assume ρ s small. Then the frst order approxmaton of the functon H(t, x, S reads as follows: H(t, x, S = S(e x 1 + ρs (φ(t, Se x φ(t, S + O(ρ 2 as ρ. (26 Proposton 3.4. Assume that the asset prce process S t = e Xt+rt fulflls SDE (18 where the Lévy measure ν s such that x 1 e2x ν (dx <. Denote by V (t, S the prce of a dervatve securty gven by [ [ V (t, S = E e r(t t Φ(S T S t = S = e r(t t E Φ(Se r(t t+xt t. (27 Assume that the pay-off functon Φ s a Lpschtz contnuous functon and the functon φ has a bounded dervatve. Then V (t, S s a soluton to the PIDE: t v(t, S2 S 2 2 V + rs 2 rv + V (t, S + H(t, x, S V (t, S H(t, x, S (t, Sν(dx =, (28 where v(t, S and H(t, x, S are gven by (2 and (21, respectvely. Proof. The asset prce dynamcs of S t under the Q measure s gven by ds t = rs t dt + v(t, S t S t dw t + H(t, x, S t J X (dt, dx. (29 If we apply Itô s lemma to V (t, S t we obtan d(v (t, S t e rt = a(tdt + dm t where a(t = t v(t, S t 2 St 2 2 V 2 + rs t rv + V (t, S t + H(t, x, S t V (t, S t H(t, x, S t (t, S tν(dx, dm t = e rt S t v(t, S t dw t + e rt V (t, S t + H(t, x, S t V (t, S t J X (dt, dx. Our goal s to show that M t s a martngale. Consequently, we have a a.s., and V s a soluton to (28 (see Proposton 8.9 of [4. To prove the term T e rt V (t, S t + 9

10 H(t, x, S t V (t, S t J S (dt, dy s a martngale t s suffcent to show that [ T ( E e 2rt V (t, S t + H(t, x, S t V (t, S t ν(dx 2 dt <. (3 Snce sup t T E [ e XT t < and the pay-off functon Φ s Lpschtz contnuous, V (t, S s Lpschtz contnuous as well wth some Lpschtz constant C >. As the functon φ(t, S has bounded dervatves we obtan S φ(t, S + H(t, x, S φ(t, S S φ H(t, x, S L H(t, x, S (see Assumpton 3.1. Snce H(t, x, S = S(e x 1 + ρs(φ(t, S + H(t, x, S φ(t, S we obtan H(t, x, S 2 S 2 (e x 1 2 /(1 ρl 2. As V s Lpschtz contnuous wth the Lpschtz constant C > we have [ T ( E e 2rt C 2 [ T (1 ρl 2 E V (t, S t + H(t, x, S t V (t, S t ν(dx 2 dt e 2rt S t 2 (e x 1 2 ν(dxdt <, because sup t [,T E [ St 2 <. Here C = (ex 1 2 ν(dx < due to the assumptons made on the measure ν. It remans to prove that T e rt S t v(t, S t (t, S tdw t s a martngale. Snce S φ (t, S s assumed to be bounded we obtan < v(t, S = σ 1 ρs φ σ (t, S 1 ρl C 1 <. Therefore E[ T e 2rt ( (t, S tv(t, S t S t 2 dt C 2 C1 2 T e 2rt E[St 2 dt < because S t s a martngale. Hence M t s a martngale as well. As a consequence, a and so V s a soluton to PIDE (28, as clamed. emark 2. If ρ = then H(t, x, S = S(e x 1 and equaton (28 reduces to: t + σ2 2 S2 2 V 2 + rs rv + V (t, Se x V (t, S S(e x 1 (t, Sν(dx =, (31 whch s the well-known classcal PIDE. If there are no jumps (ν = and a trader follows the delta hedgng strategy,.e. φ(t, S = S V (t, S, then equaton (28 reduces to the Frey Stremme opton prcng model: t + 1 ( 2 1 ϱs 2 S V 2 S2 2 V + rs 2 σ 2 rv = (32 (cf. [5. Fnally, f ρ = and ν = equaton (28 reduces to the classcal lnear Black Scholes equaton. 1

11 For smplcty, we assume the nterest rate s zero, r =. Then the functon V (t, S s a soluton to the PIDE: t v(t, S2 S 2 2 V 2 + V (t, S + H(t, x, S V (t, S H(t, x, S (t, Sν(dx =. (33 Let us defne the trackng error of a tradng strategy α t = φ(t, S t as follows: e M T := Φ(S T V = V (T, S T V T α tds t. By applyng Itô s formula to V (t, S t and usng (33 we obtan V (T, S T V = V (T, S T V (, S = = = = T T T T T dv (t, S t ds t + t v(t, S t 2 St 2 2 V 2 S dt T + V (t, S t + H(t, x, S t V (t, S t H(t, x, S t J X(dt, dx T ds t V (t, S t + H(t, x, S t V (t, S t H(t, x, S t ν(dxdt T + V (t, S t + H(t, x, S t V (t, S t H(t, x, S t J X(dt, dx T ds t + V (t, S t + H(t, x, S t V (t, S t H(t, x, S t J X (dt, dx. Usng expresson (29 for the dynamcs of the asset prce S t (wth r =, the trackng error e M T can be expressed as follows: T e M T = V (T, S T V = T + T T + α t ds t = T ( (t, S t α t ds t V (t, S t + H(t, x, S t V (t, S t H(t, x, S t J X (dt, dx ( v(t, S t S t α t dw t (34 V (t, S t + H(t, x, S t V (t, S t α t H(t, x, S t J X (dt, dx. emark 3. For the delta hedgng strategy α t = φ(t, S t = (t, S t the trackng error functon e M T can be expressed as follows: e M T = T V (t, S t + H(t, x, S t V (t, S t H(t, x, S t (t, S t J X (dt, dx. Clearly, the trackng error for the delta hedgng strategy need not be zero for ν. Next, we propose a crteron that can be used to fnd the optmal hedgng strategy. 11

12 Proposton 3.5. The tradng strategy α t = φ(t, S t of a large trader mnmzng the varance E [ (ɛ M T 2 of the trackng error s gven by the mplct equaton: φ(t, S t = β ρ (t, S t [ v(t, S t 2 St 2 (t, S t + (V (t, S t + H(t, x, S t V (t, S t H(t, x, S t ν(dx, (35 where β ρ (t, S t = 1/[v(t, S t 2 S 2 t + H(t, x, S t 2 ν(dx and H(t, x, S = S(e x 1 + ρs[φ(t, S + H(t, x, S φ(t, S. Proof. Usng expresson (34 for the trackng error ɛ M T ( 2 (t, S t α t dt E [ [ (ɛ M T 2 T = E v(t, S t 2 St 2 [ T +E and Itô s sometry we obtan (V (t, S t + H(t, x, S t V (t, S t α t H(t, x, S t 2 ν(dxdt. The mnmzer α t of the above convex quadratc mnmzaton problem satsfes the frst order necessary condtons d(e [ ɛ 2 T, αt =, that s, [ T = 2E ( ( v(t, S t 2 St 2 (t, S t α t + H(t, x, S t ( V (t, S t + H(t, x, S t V (t, S t α t H(t, x, S t ν(dx ω t dt for any varaton ω t. Thus the trackng error mnmzng strategy α t s gven by (35. emark 4. The optmal tradng strategy mnmzng the varance of the trackng error need not satsfy the structural Assumpton 3.1. For nstance, f ν = then the trackng error mnmzer s just the delta hedgng strategy φ = S V. In the case of a call or put opton ts gamma,.e. 2 V (t, S becomes nfnte as t T and S = K. Gven a level L > we can however mnmze the trackng error E [ ɛ 2 T under the addtonal constrant sup S> S φ (t, S L. That s we can solve the followng convex constraned nonlnear optmzaton problem mn φ E [ ɛ 2 T s.t. S φ L nstead of the unconstraned mnmzaton problem proposed n Proposton 3.5. emark 5. Notce that, f ν = and ρ, the tradng strategy α t reduces to the Black Scholes delta hedgng strategy,.e. α t = (t, S t. If ν and ρ =, then the optmal tradng strategy becomes α t = φ (t, S t where ( φ (t, S t = β (t, S t σ 2 St 2 (t, S t + S t (e x 1 (V (t, S t e x V (t, S t ν(dx, where β (t, S t = 1/[σ 2 S 2 t + S2 t (e x 1 2 ν(dx. 12

13 We conclude ths secton by the followng proposton provdng the frst order approxmaton of the trackng error mnmzng tradng strategy for the case when the parameter ρ 1 s small. In what follows, we derve the frst order approxmaton of φ ρ (t, S t n the form φ ρ (t, S t = φ (t, S t + ρφ 1 (t, S t + O(ρ 2 as ρ. Clearly, the frst order Taylor expanson for the volatlty functon v(t, S has the form: v(t, S 2 = σ 2 (1 ρs S φ 2 = σ2 + 2ρσ 2 S φ (t, S + O(ρ2, as ρ. Wth regard to Proposton 3.3 (see (26 we have H(t, x, S = H (t, x, S + ρh 1 (t, x, S + O(ρ 2, where H (t, x, S = S(e x 1, H 1 (t, x, S = S[φ (t, Se x φ (t, S. (36 The functon β ρ can be expanded as follows: β ρ (t, S = β (t, S + ρβ (1 (t, S + O(ρ 2, β (t, S = 1/[σ 2 S 2 + S 2 (e x 1 2 ν(dx, (37 [ β (1 (t, S = (β (t, S 2 2σ 2 S 3 φ (t, S + 2S2 (e x 1[φ (t, Se x φ (t, Sν(dx. Usng the frst order expansons of the functons v 2, β ρ and H we obtan the followng results. Proposton 3.6. For small values of the parameter ρ 1, the trackng error varance mnmzng strategy α t = φ ρ (t, S t s gven by where φ (1 (t, S = β (t, S + φ ρ (t, S t = φ (t, S t + ρφ (1 (t, S t + O(ρ 2, as ρ, (38 [ 2σ 2 S 3 (t, S φ (t, S ( V (t, Se x V (t, S + +β (1 (t, S [ σ 2 S 2 (t, S + (t, Sex H (t, x, S H 1 (t, x, Sν(dx (V (t, Se x V (t, S H (t, x, Sν(dx and the functons H, H 1, β and β (1 are defned as n (36 and ( Implct-explct numercal dscretzaton scheme The am of ths secton s to propose a full tme-space dscretzaton scheme for solvng the nonlnear PIDE (28. The method of dscretzaton s based on a fnte dfference approxmaton of all dervatves occurrng n (28 and approxmaton of the ntegral term by means of the trapezodal ntegraton rule on a truncated doman. 13

14 In order to solve (28 we transform t nto a nonlnear parabolc PIDE defned on. Indeed, usng the standard transformatons V (t, S = e rτ u(τ, x, φ(t, S = ψ(τ, x where τ = T t, x = ln( S K we conclude that V (t, S s a soluton to (28 f and only f the functon u(τ, x solves the followng nonlnear parabolc equaton: u τ σ 2 2 u (1 ρ ψ x 2 2 x + ( σ 2 = 1 r (1 ρ ψ (39 x 2 + u(τ, x + ξ(τ, z, x u(τ, x H(T τ, z, Ke x 1 u e x (τ, xν(dz, K x u(, x = h(x Φ(Ke x, (τ, x [, T, (4 u x and H(t, z, S = S(e z 1 + ρs[φ(t, S + H(t, z, S φ(t, S, (41 ξ(τ, z, x = ln(1 + K 1 e x H(T τ, z, Ke x. ( A numercal scheme for solvng nonlnear PIDEs wth fnte actvty Lévy measures We frst consder the case when the Lévy measure ν has a fnte actvty,.e. ν( <. Let us denote λ = ν(dz, and ω(τ, x = H(T τ, z, S e x 1 e x ν(dz. S We have λ <. Observe that (39 s equvalent to u τ = 1 2 σ 2 2 u (1 ρ ψ x 2 2 x + ( r 1 2 σ 2 (1 ρ ψ x 2 ω u x λu+ u(τ, x+ξ(τ, z, xν(dz. (43 We proceed to solve (43 by means of the sem-mplct fnte dfference scheme proposed n [15. The dea s to separate the rght-hand sde nto two parts: the dfferental part and the ntegral part. Let u j = u(τ j, x, τ j = j τ, x = z = x for = N + 1,, N 1 and j = 1,, M. We approxmate the dfferental part mplctly except of ψ(τ, x ( u x ( 2 u x 2 ( u τ j j j { u j+1 +1 u j+1 x, f (σj 2 2 r + ω j <, u j+1 u j+1 1 x, f (σj 2 2 r + ω j, σ j = σ 1 ρdψ j, uj uj+1 + u j+1 1 ( x 2, uj+1 u j. t ( ψ x j ψj +1 ψj x = Dψ j. As for the ntegral operator, frst we have to truncate the ntegraton doman to a bounded nterval [B l, B r. We approxmate ths ntegral by choosng ntegers K l and 14

15 K r such that [B l, B r [(K l 1/2 x, (K l + 1/2 x. Then Br B l u(τ j, x + ξ(τ j, z, x ν(dz K r k=k l u(τ j, x + ξ(τ j, z k, x ν k, (44 where ν k = 1 2 ( ν(zk+1/2 + ν(z k 1/2 x. Analogously, ω j e x K K r k=k l H(T τ j, z k, Ke x ν k, and λ K r k=k l ν k, where ξ(τ, z, x s gven as n (42. Insertng the fnte dfference approxmatons of dervatves of u nto (43 we obtan u j+1 u j t = 1 uj+1 2 (σj +1 2uj+1 + u j ( x 2 λu j+1 (45 +(r 1 2 (σj 2 ω j uj+1 +1 uj+1 x + K r k=k l u(τ j, x + ξ(τ j, z k, x ν k, provded that 1 2 (σj 2 r + ω j <. Smlarly, we can derve a dfference equaton for the case 1 2 (σj 2 r + ω j. If we defne coeffcents βj ±, and βj as follows: β j ± = τ 2( x 2 (σj 2 τ ( r 1 ± x 2 (σj 2 ω j, (46 β j = 1 + τλ (β j + βj +, (47 where (a + = max(a,, (a = mn(a,. Then the trdagonal system of lnear equatons for the soluton u j = (u j N+1,, uj N 1 T, j =,, M, reads as follows: where β j + uj βj uj+1 + β j uj+1 u = h(x, for = N + 1,, N 1, u j+1 = g(τ j+1, x, for = N + 1,, N/2 1, Kr 1 = u j + τ k=k l u(τ j, x + ξ(τ j, z k, x ν k, (48 for = N/2 + 1,, N/2 1, u j+1 = g(τ j+1, x, for = N/2,, N 1, ξ(τ j, z k, x = ln(1 + K 1 e x H(T τ j, z k, Ke x, and g s a functon of ponts x lyng outsde the localzaton nterval. Followng Proposton n [15 the recommended choce s g(τ, x = h(x + rτ = Φ(Ke rτ+x. The 15

16 term u(τ j, x + ξ(τ j, z k, x enterng the sum n the rght-hand sde of (48 s approxmated by means of the frst order Taylor seres expanson: u(τ j, x + ξ(τ j, z k, x u j + uj +1 uj x ξ(τ j, z k, x Numercal scheme for solvng nonlnear PIDEs wth nfnte actvty Lévy measures Next we consder the case when the Lévy measure has an nfnte actvty, e.g. the Varance Gamma process where ts Lévy densty explodes at zero and ν( =. Equaton (39 s equvalent to u τ ( = 1 σ 2 2 u 2 (1 ρ ψ x 2 2 x + r 1 σ 2 2 (1 ρ ψ ω u x 2 x + u(τ, x + ξ(τ, z, x u(τ, xν(dz. (49 Equaton (49 dffers from (43 as the term u(τ, x s contaned n the ntegral part because λ = ν(dz =. Proceedng smlarly as for dscretzaton of (43 we can solve (49 numercally by means of a sem-mplct fnte dfference scheme. If the coeffcents β j ± are defned as n (46 and βj = 1 (βj +βj +, then the soluton vector u j = (u j N+1,, uj N 1 T, j =,, M, s a soluton to the followng trdagonal system of lnear equatons: β j + uj βj uj+1 + β j uj+1 u = h(x, for = N + 1,, N 1, u j+1 = g(τ j+1, x, for = N + 1,, N/2 1, (5 Kr 1 = u j + τ k=k l (u(τ j, x + ξ(τ j, z k, x u(τ j, x ν k, for = N/2 + 1,, N/2 1, u j+1 = g(τ j+1, x, for = N/2,, N 1. The term u(τ j, x + ξ(τ j, z k, x u(τ j, x enterng the sum n the rght-hand sde of (5 s agan approxmated by means of the frst order Taylor seres expanson,.e. u(τ j, x + ξ(τ j, z k, x u(τ j, x uj +1 uj x ξ(τ j, z k, x. 5. Numercal results In ths secton we present results of numercal experments usng the fnte dfference scheme descrbed n Secton 4 for the case of a European put opton,.e. Φ(S = (K S +. As for the Lévy process we consder the Varance Gamma process wth parameters θ =.33, σ =.12, κ =.16, and other opton prcng model parameters: r =, K = 1, T = 1. Numercal dscretzaton parameters were chosen as follows: 16

17 PIDE 25 PIDEFrey 2 BS 2 Frey S ClassPIDE S Fgure 1. Comparson of European put opton prces between the classcal PIDE and the lnear Black Scholes model (left. Comparson between the classcal PIDE and the Frey Stremme PIDE model (rght. Impled Volatlty PIDEFEY Frey PIDE Fgure 2. Comparson of mpled volatltes between the Frey Stremme model, classcal PIDE and Frey Stremme PIDE generalzaton. x =.1, t =.5. Snce the Varance Gamma process has nfnte actvty, we employ numercal dscretzaton scheme descrbed n Secton 4.2. In what follows, we present varous opton prces computed by means of the fnte-dfference numercal scheme descrbed n Secton 4 for the lnear Black Scholes (ρ = and the Frey Stremme model (ρ > and ther jump-dffuson PIDE generalzatons. In Fg. 1 we show comparson of European put opton prces between the classcal PIDE and the lnear Black Scholes model, and comparson between the classcal PIDE and the Frey Stremme PIDE model for the case when the large trader s nfluence s small, ρ =.1. In Fg. 2 we plot dependence of the mpled volatltes as decreasng functons of the strke prce K for the Frey Stremme model and ts PIDE generalzatons. We can observe that the mpled volatltes for the Frey Stremme PIDE model s always hgher when varyng the strke prce of the European Put opton. Numercal values of opton prces for varous models and parameter settngs are summarzed n Tables 1 and 2. Numercal results confrm our expectaton that assumng rsk arsng from sudden jumps n the underlyng asset process yelds a hgher opton prce when comparng to the Frey Stremme model opton prce. In Fg. 3 (left we compare European put opton prces V (, S computed by means of the Black Scholes and Frey-Stremme models dependng on the parameter ρ measurng nfluence of a large trader. We can observe that the prce of the European put opton ncreases wth respect to ρ, as expected. Furthermore, the prce computed from 17

18 Table 1. European put opton prces V (, S for the Black-Scholes and Frey Stremme models wth ρ =.1 and ther PIDE generalzatons. B-S F-S B-S PIDE F-S PIDE S ν =, ρ = ν =, ρ ν, ρ = ν, ρ Table 2. European put opton prces V (, S for the Frey Stremme and Frey-Stremme PIDE models for varous values of ρ. F-S F-S PIDE F-S F-S PIDE F-S F-S PIDE S ρ =.1 ρ =.1 ρ =.2 ρ =.2 ρ =.3 ρ = Frey Ρ.3 15 Frey Ρ.2 Frey Ρ.1 BS PIDE PIDEFrey Ρ.1 PIDEFrey Ρ.2 PIDEFrey Ρ.3 Fgure 3. Comparson of European put opton prces for the Black Scholes and the Frey Stremme models (left and Frey-Stremme PIDE model for varous ρ. 18

19 the Frey Stremme PIDE model s larger than the one obtaned from the lnear Black Scholes equaton. Moreover, the prce computed from Frey Stremme PIDE model s hgher than the one computed by means of the nonlnear Frey Stremme model. Ths s due to the fact that the jump part of the underlyng asset process enhances rsk, and, consequently ncreases the opton prce. Fg. 3 (rght shows comparson of the opton prces for the Black Scholes and Frey Stremme PIDE model for varous values of ρ. 6. Concluson In ths paper we nvestgated a novel nonlnear opton prcng model generalzng the Frey Stremme model under the assumpton that the underlyng asset prce follows a Lévy stochastc process. We derved the fully-nonlnear PIDE for prcng optons under nfluence of a large trader. We also proposed the hedgng strategy mnmzng the varance of the trackng error. We derved a sem-mplct fnte dfference numercal approxmaton scheme for solvng the nonlnear PIDE. We presented varous numercal experments llustratng the nfluence of the large trader under the Lévy process wth jumps. 7. Acknowledgements The research was supported by the project CEMAPE MULTI/491 fnanced by FCT/MEC through natonal funds and the Slovak esearch Project VEGA 1/62/18. eferences [1 Davd Applebaum, Lévy processes and stochastc calculus, second ed., Cambrdge Studes n Advanced Mathematcs, vol. 116, Cambrdge Unversty Press, Cambrdge, 29. [2 Marco Avellaneda, Avron Levy, and Antono Paras, Prcng and hedgng dervatve securtes n markets wth uncertan volatltes, Appled Mathematcal Fnance 2 (1995, [3 Ole E. Barndorff-Nelsen and Serge Z. Levendorskĭ, Feller processes of normal nverse Gaussan type, Quant. Fnance 1 (21, [4 ama Cont and Peter Tankov, Fnancal Modelng wth Jump Processes, Chapman & Hall/CC, Boca aton, FL, 24. [5 üdger Frey, Perfect opton hedgng for a large trader, Fnance Stochast 2 (1998, [6 üdger Frey and Perre Pate, sk management for dervatves n llqud markets: a smulaton study, Advances n fnance and stochastcs, Sprnger, Berln, 22, pp M [7 üdger Frey and Alexander Stremme, Market volatlty and feedback effects from dynamc hedgng, Math. Fnance 7 (1997, no. 4, M [8 Martn Jandačka and Danel Ševčovč, On the rsk-adjusted prcng-methodologybased valuaton of vanlla optons and explanaton of the volatlty smle, J. Appl. Math. (25, no. 3, M

20 [9 Steven Kou, A jump-dffuson model for opton prcng, Management Scence 48 (22, [1 Yue-Kuen Kwok, Mathematcal models of fnancal dervatves, second ed., Sprnger Fnance, Sprnger, Berln, 28. M [11 Dlp B. Madan, Peter Carr, and Erc C. Chang, The Varance Gamma Process and Opton Prcng, European Fnance evew 2 (1998, [12 obert C. Merton, Opton prcng when underlyng stock returns are dscontnuous., J. Fnancal Economcs 3 (1976, [13 Phlpp J. Schönbucher and Paul Wlmott, The feedback effect of hedgng n llqud markets, SIAM J. Appl. Math. 61 (2, no. 1, M [14 K. onne Srcar and George Papancolaou, General Black-Scholes models accountng for ncreased market volatlty from hedgng strateges, Appled Mathematcal Fnance 5 (1998, no. 1, [15 Ekaterna Voltchkova, Equatons ntégro-dfférentelles d évoluton: méthodes numérques et applcatons en fnance, Ph.D. thess, École Polytechnque, 25. 2