arxiv: v1 [math.ap] 10 Apr 2008

Transcription

1 Optimal systems of subalgebras and invariant solutions for a nonlinear Black-Scholes equation arxiv: v1 [math.ap] 10 Apr 2008 Maxim Bobrov Halmstad University, Box 823, Halmstad, Sweden Abstract The main object of our study is a four dimensional Lie algebra which describes the symmetry properties of a nonlinear Black-Scholes model. This model implements a feedback effect which is typical for an illiquid market. The structure of the Lie algebra depends on one parameter, i.e. we have to do with one parametric family of algebras. We provide a classification of these algebras using Patera Winternitz method. Optimal systems of one-, two- and three- dimensional subalgebras are described for the family of symmetry algebras of the nonlinear Black-Scholes equation. The optimal systems give the possibility to describe a complete set of invariant solutions to equation. Key words and phrases: Black - Scholes model, nonlinearity, optimal system AMS classification: 35K55 1 Introduction In [2] Frey and Patie examined the feedback effect of the option replication strategy of the large trader on the asset price process. They obtain a new model by introduction of a liquidity coefficient which depends on the current stock price. The feedback-effect described leads to a nonlinear version of the Black-Scholes partial differential equation, u t + σ2 S 2 2 u SS = 0, (1) (1 ρsλ(s)u SS ) 2 with S [0, ), t [0, T]. As usual, S denotes the price of the underlying asset and u(s, t) denotes the hedge-cost of the claim with a payoff h(s) which will be defined later. The hedge-cost is different from the price of the derivatives product in illiquid markets. In the sequel t is the time variable, σ defines the volatility of the underlying asset. The proof is based on methods of Lie point symmetries, i.e. the Lie symmetry algebras and groups to the corresponding equations were found; see [1] for details and [4, 5, 3] for a general introduction to the methodology. 1

2 Theorem 1.1. The differential equation (1) with an arbitrary function λ(s) possesses a trivial three dimensional Lie algebra Diff (M) spanned by infinitesimal generators V 1 = t, V 2 = S V 3 = u. Only for the special form of the function λ(s) ωs k, where ω, k R equation (1) admits a non-trivial four dimensional Lie algebra L spanned by generators V 1 = t, V 2 = S V 3 = V 4 = S S with commutator relations + (1 k)u u (2) [V 1, V 2 ] = [V 1, V 3 ] = [V 1, V 4 ] = [V 2, V 3 ] = 0, [V 2, V 4 ] = kv 2, [V 3, V 4 ] = (1 k)v 3. (3) In the paper [8] authors try to construct the optimal system of subalgebras for the symmetry algebra (2). The authors used the method suggested in the series of well known papers by P.Winternitz and Co [6, 7] where all three and four dimensional Lie algebras where classified. The investigation in [8] contains some misprints and mistakes which demand corrigendum. In the paper the structure of the optimal system of subalgebras do not contains some of one dimensional algebras. Other sides the classification do not depends on the parameter k from which the structure of the algebra deeply depends. To be able to construct correct families of invariant solutions we need the correct optimal system of subalgebras. In our paper we present the correct optimal system of one-, two-, three- dimensional systems of subalgebras. 2 Classification of the algebra L Let us consider the following four dimensional Lie algebra L v 1 = t, v 2 v 3 = v 4 + (1 k)u x u. with the following commutator table We want to classify this algebra. We use the classification which was introduce 2

3 v 1 v 2 v 3 v 4 v v kv 2 v (1 k)v 3 v 4 0 kv 2 (1 k)v 3 0 Table 1: The commutator table of the algebra L by J. Patera and P. Winternitz [6]. As we said before the structure of the algebra L depend on the parameter k. We can see that we have two special cases for this algebra k = 0 and k = 1. Let us consider case k = 0 in this case we obtain following generators v 1 = t, v 2 v 3 = v 4 x + u u. (4) with the following commutator table v 1 v 2 v 3 v 4 v v v v 3 v v 3 0 Table 2: The commutator table of the algebra L 4.1 Let us consider the algebra L 4.1 with operators e 1 = v 4, e 2 = v 3, e 3 = v 1, e 4 = v 2. (5) Then L 4.1 = A 2 2A 1 with commutator relations [e 1, e 2 ] = e 2. The case k = 1 leads us to the algebra L 4.2 with generators v 1 = t, v 2 v 3 = v 4 x. (6) which is isomorphic to the algebra L 4.2 spanned on generators e 1 = v 4, e 2 = v 2, e 3 = v 1, e 4 = v 3. (7) 3

4 In this case algebra L 4.2 has the same decomposition A 2 2A 1 with the same commutator relations [e 1, e 2 ] = e 2. In the case k > 1 2 the algebra L isomorphic to the algebra L 4.3 = Aα 3.5 A 1 with commutator relations [e 1, e 3 ] = e 1, [e 2, e 3 ] = αe 2 where and α = k 1 k. e 1 = v 2, e 2 = v 3, e 3 = 1 k v 4, e 4 = v 1, (8) The case k < 1 2 leads us to the algebra L 4.4 = A α 3.5 A 1 with commutator relations [e 1, e 3 ] = e 1, [e 2, e 3 ] = αe 2 where and α = k k 1. e 1 = v 3, e 2 = v 2, e 3 = 1 1 k v 4, e 4 = v 1, (9) The last case k = 1 2 where L with generators and commutator table v 1 = t, v 2 v 3 = v 4 x u (10) v 1 v 2 v 3 v 4 v v v v 3 v v 4 0 v 2 2 1v Table 3: The commutator table of the algebra L isomorphic to the algebra L 4.5 = A 3.4 A 1 with generators e 1 = v 3, e 2 = v 2, e 3 = 2v 4, e 4 = v 1 (11) and the commutator relations [e 1, e 3 ] = e 1, [e 2, e 3 ] = e 2 4

5 3 Optimal system of subalgebras The main goal of this paper is to find the correct optimal system of subalgebras for the Lie algebra L. The procedure was described by Pattera & Winternitz in [6]. In the paper all three and four-dimensional algebras were classified and the optimal systems for these algebras were listed. We repeat this algorithm for the case of the algebra L in case k 0, 1, 1. In this case the Lie algebra L is 2 isomorphic to the algebra A α 3.5 A 1 with following commutator table e 1 e 2 e 3 e 4 e e 1 0 e αe 2 0 e 3 e 1 αe e Table 4: The commutator table of the algebra L where 0 < α < 1. In case k 0, 1, 1 the Lie algebra L has one cental element and can be represent 2 as a direct sum of one- and three- dimensional Lie algebras L = {e 4 } L 3, (12) where e 4 is a central element of the algebra L and L 3 = L \ {e 4 }. The representation (12) simplify the procedure of construction of the optimal system of subalgebras. We start with construction corresponding system of subalgebras for both algebras in (12) and then complete the study with the investigation of non-splitting extensions. We follow the paper ([6]) and describe solution of this problem in step-by-step method. Step 1. We find all subalgebras of e 4 it is obviously only two subalgebras {0} and {e 4 } Step 2. We have to classify all subalgebras of L 3 under conjugation which defined by interior isomorphism of the algebra. This isomorphism can be presented by adjoint representation. Definition 3.1 (Olver [5]). Let G be a Lie group. For each g G, group conjugation K g (h) = ghg 1, h G, determines a diffeomorphism on G. Moreover, K g K g = K gg, K e = 1 G, so K g determines a global group action of G on itself, with each conjugacy map K g being a group homomorphism: K g (hh ) = K g (h)k g (h ) etc. The differential dk g : TG h TG Kg(h) is readily seen to preserve the right-in variance of vector fields, and hence determines a linear map on the Lie algebra of G, called the adjoint representation: Ad g(v) = dk g (v) (13) 5

6 The simplest way to represent the adjoint representation is Lie series Ad(exp(εv))w = w ε[v, w] + ε2 [v, [v, w]]... (14) 2! The adjoint representation table for the algebra L is following Ad e 1 e 2 e 3 e 4 e 1 e 1 e 2 e 3 εe 1 e 4 e 2 e 1 e 2 e 3 αεe 2 e 4 e 3 e ε e 1 e αε e 2 e 3 e 4 e 4 e 1 e 2 e 3 e 4 Table 5: The adjoint representation table of the algebra L = A α 3.5 A 1 with (i, j)-th entry indicating Ad(exp(εe i ))e j By using the adjoint representation we want to classify all subalgebras of L 3 under conjugacy. Firstly we consider one-dimensional subalgebra of the general type A = ae 1 + be 2 + ce 3. If c 0 then we can use first and second lines of the Table 5. Ad(exp(ξe 1 + ζe 2 ))A = (a cξ)e 1 + (b cαζ)e 2 + ce 3 (15) with ξ = a c, ζ = b cα we can see that A isomorphic to e 3. If c = 0 we have three cases. If a 0, b = 0 then A isomorphic to e 1, if a = 0, b 0 then A isomorphic to e 2 and the last one ab 0 then we can use the third line of the Table 5 Ad(exp(ξe 3 ))A = ae ξ e 1 + be αξ e 2 (16) by division on ae ξ and choose ξ = 1 log a α 1 b we can see that A isomorphic to e 1 ± e 2 Then the optimal system of one-dimensional subalgebras of L 3 is following {0}, {e 1 }, {e 2 }, {e 3 }, {e 1 ± e 2 }. (17) Let us consider now two dimensional subalgebras of L 3. Let B is one of subalgebras (17) and A = ae 1 + be 2 + ce 3. For a subalgebra M = {B, A} we want [A, B] M. Let B = e 1 then without loss of generality we can represent A in form ae 2 + be 3. Let b 0, by using second line of the adjoint representation table we can see that A isomorphic to e 3. If b = 0 then A = e 2. In this case we obtain to subalgebras {e 1, e 2 }, {e 1, e 3 } (18) 6

7 By the same way we obtain subalgebras {e 1, e 2 }, {e 2, e 3 } in case B = e 2. Let B = e 3 then without loss of generality we can represent A = ae 1 + be 2. Let us check the commutator relations [ae 1 + be 2, e 3 ] = ae 1 + αbe 2. We can see that the pair {A, B} is algebra just under condition ab = 0. Then we obtain two subalgebras {e 1, e 3 }, {e 2, e 3 } (19) In case B = e 1 ± e 2 we can say that A = e 3. Then [e 1 ± e 2, e 3 ] = e 1 ± αe 2 and it obviously is not an algebra. Then the optimal system of subalgebras of L 3 is following {0}, {e 1 }, {e 2 }, {e 3 }, {e 1 ± e 2 }, {e 1, e 2 }, {e 1, e 3 }, {e 2, e 3 }. (20) Step 3. We have to find all splitting extensions of the algebra {e 4 }. To do this we have to find all subalgebras N a of L 3 such that and classify all such subalgebras under Nor L e 4. [e 4, N a ] N a (21) Definition 3.2. Let L be a Lie algebra. The set of generators y of the algebra L such that [x, y] N (22) for any x N and y L called normalizer N in algebra L or Nor L N As soon as e 4 is a central element and Nor L e 4 = L any adjoint representation does not affect on e 4 and N a is any subalgebra of L 3. This step is trivial and we obtain the subalgebras of the type {e 4, S} where S running through subalgebras (20). Step 4. We have to find all subalgebras of type { e 4 + } a i e i, N a i (23) where N a is a subalgebra of L 3 with an normalizer not contained in L 3 and x i R are not all equal to zero. Generators e 4 + x i e i must not be conjugate to e 4. i Hence e 4 is an central element of L all of conditions are correct. Let N a going through the list of algebras (17) and let A = e 4 + a i e i. Let N a = {0}. This i case is trivial in cause e 4 is a central element and the procedure is similar to the second step. We obtain four subalgebras. {ae 1 + e 4 }, {ae 2 + e 4 }, {ae 3 + e 4 }, {a(e 1 ± e 2 ) + e 4 }, (24) 7

8 where a 0. We divided by a all of those subalgebras. Then we can see that in two first one we can use adjoint representation generated by e 1 and e 2. We obtain {e 1 + be ε e 4 }, {e 2 + be αε e 4 }, (25) where b = 1 a 0 By choosing ε = log 1 b in first the case and ε = 1 α log 1 b in the second one, we obtain finally following list of one dimensional non-splitting extensions {e 1 ± e 4 }, {e 2 ± e 4 }, {e 3 + ae 4 }, {e 1 ± e 2 + ae 4 }, (26) where a 0. Let us now consider a two dimensional non-splitting extensions. To simplify this procedure we use as N a subalgebras (20) and e 4 + a i e i we left i in the general type. We can do it as soon as by simplification under conjugacy subalgebras N a the general view of e 4 + a i e i does not change. Let N a be equal i to e 1 then without loss of generality we can represent A = e 4 + a 2 e 2 + a 3 e 3. If a 3 does not equal to zero we can use the second line of adjoint representation table and reduce A to e 4 + ae 3. And we obtain the following subalgebra {e 3 + ae 4, e 1 }, (27) where a 0. On the other case we rewrite A = ae 4 + e 2 and use third line of the adjoint representation table then we obtain or {e αε e 2 + ae 4, e ε e 1 } (28) {e 2 + ae αε e 4, e (1 α)ε e 1 } (29) where a 0. By choosing ε = 1 log a and divided second generator of the α algebra below on corresponding constant we obtain the following algebra {e 2 ± e 4, e 1 } (30) The same procedure for the case N a = e 2 lead us to subalgebras {e 1 ± e 4, e 2 }, {e 3 + ae 4, e 2 }, (31) where a 0. Let us consider the case N a = e 3, then we can choose A = e 4 + a 1 e 1 +a 2 e 2. Note that {A, N a } is algebra just in case a 1 a 2 = 0. Those subalgebras we consider in previous cases. Let N a = e 1 ±e 2 then {N a, e 4 +a 1 e 1 +a 2 e 2 +a 3 e 3 } is subalgebra only in case a 3 = 0. Without loss of generality we can represent A = ae 4 + e 1 then by using third line we can see that algebras {A, N a } are isomorphic to the following algebras {e 1 ± e 4, e 1 + ae 2 }, (32) 8

9 where a R. Note that case a = 0 we consider on the third step then finally we obtain the following algebra {e 1 ± e 4, ae 1 + e 2 }, (33) where a 0. Now we consider case N a = {e 1, e 2 } in this case we can represent A as ae 4 + e 3 and we obtain following three dimensional subalgebra {e 3 + ae 4, e 1, e 2 }, (34) where a 0. It is easy to see that other cases do not provide any other subalgebras. Then we obtain the following list of the optimal system of subalgebras of the algebra L. {0}, {e 1 }, {e 2 }, {e 3 }, {e 1 ± e 2 }, {e 1, e 2 }, {e 1, e 3 }, {e 2, e 3 }, {e 4 }, {e 1, e 4 }, {e 2, e 4 }, {e 3, e 4 }, {e 1 ± e 2, e 4 }, {e 1, e 2, e 4 }, {e 1, e 3, e 4 }, {e 2, e 3, e 4 }, {e 1 ± e 4 }, {e 2 ± e 4 }, {e 3 + ae 4 }, {e 1 ± e 2 + ae 4 }, {e 3 + ae 4, e 1 }, {e 2 ± e 4, e 1 }, {e 1 ± e 4, e 2 }, {e 3 + ae 4, e 2 }, {e 1 ± e 4, ae 1 + e 2 }, {e 3 + ae 4, e 1, e 2 }. We have to remark that some of those subalgebras are partial cases of other more general then after simplification of this list finally we obtain the following optimal system Dimension Subalgebras 1 {e 2 }, {e 4 }, {e 1 + ae 2 }, {e 1 + ǫe 4 }, {e 2 + ǫe 4 }, {e 3 + ae 4 }, {e 1 + ǫe 2 + ae 4 } 2 {e 1, e 2 }, {e 1, e 4 }, {e 2, e 4 }, {e 3, e 4 }, {e 1 + ǫe 2, e 4 } {e 2 + ǫe 4, e 1 }, {e 1 + ǫe 4, ae 1 + e 2 }, {e 3 + ae 4, e 1 }, {e 3 + ae 4, e 2 } 3 {e 1, e 2, e 4 }, {e 1, e 3, e 4 }, {e 2, e 3, e 4 }, {e 1, e 2, e 3 + ae 4 }, Table 6: The Optimal System of subalgebras of the algebra L with a R, ǫ = ±1 9

10 We remark here that in case k = 1 the structure of the algebra L is the same 2 as in the case above and the optimal system of subalgebras is the same. For k = 0 or k = 1 by the same procedure we obtain the following system of subalgebras Dimension Subalgebras 1 {e 2 }, {e 3 cos ϕ + e 4 sin ϕ}, {e 1 + a(e 3 cosϕ + e 4 sin ϕ)}, {e 2 + ǫ(e 3 cosϕ + e 4 sin ϕ)} 2 {e 1 + a(e 3 cosϕ + e 4 sin ϕ), e 2 }, {e 3, e 4 }, {e 1 + a(e 3 cos ϕ + e 4 sin ϕ), e 3 sin ϕ e 4 cosϕ}, {e 2 + ǫ(e 3 cos ϕ + e 4 sin ϕ), e 3 sin ϕ e 4 cos ϕ}, {e 2, e 3 sin ϕ e 4 cosϕ} 3 {e 1, e 3, e 4 }, {e 2, e 3, e 4 }, {e 1 + a(e 3 cosϕ + e 4 sin ϕ), e 3 sin ϕ e 4 cosϕ, e 2 }, Table 7: where a R, ǫ = ±1, 0 ϕ Π 4 Results In this chapter we return to our starting algebra and introduce optimal system in starting generators. 4.1 Algebra L with k = 0 We obtain the algebra spanned on generators v 1 = t, v 2 v 3 = v 4 x + u u. The optimal system of subalgebras for this algebra is following 10

11 Dimension Subalgebras 1 {v 3 }, {v 1 cos ϕ + v 2 sin ϕ}, {v 4 + a(v 1 cosϕ + v 2 sin ϕ)}, {v 3 + ǫ(v 1 cosϕ + v 2 sin ϕ)} 2 {v 1 + a(v 1 cosϕ + v 2 sin ϕ), v 3 }, {v 1, v 2 }, {v 4 + a(v 1 cos ϕ + v 2 sin ϕ), v 1 sin ϕ v 2 cosϕ}, {v 3 + ǫ(v 1 cos ϕ + v 2 sin ϕ), v 1 sin ϕ v 2 cosϕ}, {v 3, v 1 sin ϕ v 2 cosϕ} 3 {v 4, v 1, v 2 }, {v 3, v 1, v 2 }, {v 4 + a(v 1 cosϕ + v 2 sin ϕ), v 1 sin ϕ v 2 cosϕ, v 3 }, Table 8: where a R, ǫ = ±1, 0 ϕ Π 4.2 Algebra L with k = 1 We obtain the algebra spanned on generators v 1 = t, v 2 v 3 = v 4 x. The optimal system of subalgebras for this algebra is following Dimension Subalgebras 1 {v 2 }, {v 1 cos ϕ + v 3 sin ϕ}, {v 4 + a(v 1 cosϕ + v 3 sin ϕ)}, {v 2 + ǫ(v 1 cosϕ + v 3 sin ϕ)} 2 {v 4 + a(v 1 cosϕ + v 3 sin ϕ), v 2 }, {v 1, v 3 }, {v 4 + a(v 1 cos ϕ + v 3 sin ϕ), v 1 sin ϕ v 3 cosϕ}, {v 2 + ǫ(v 1 cos ϕ + v 3 sin ϕ), v 1 sin ϕ v 3 cosϕ}, {v 2, v 1 sin ϕ v 3 cosϕ} 3 {v 4, v 1, v 3 }, {v 2, v 1, v 3 }, {v 4 + a(v 1 cosϕ + v 3 sin ϕ), v 1 sin ϕ v 3 cosϕ, v 2 }, Table 9: where a R, ǫ = ±1, 0 ϕ Π 11

12 4.3 Algebra L with k < 1 2 and k 0 We obtain the algebra spanned on generators v 1 = t, v 2 v 3 = v 4 + (1 k)u x u. The optimal system of subalgebras for this algebra is following Dimension Subalgebras 1 {v 2 }, {v 1 }, {v 3 + av 2 }, {v 3 + ǫv 1 }, {v 2 + ǫv 1 }, {v 4 + av 1 }, {v 3 + ǫv 2 + av 1 } 2 {v 3, v 2 }, {v 3, v 1 }, {v 2, v 1 }, {v 4, v 1 }, {v 3 + ǫv 2, v 1 } {v 2 + ǫv 1, v 3 }, {v 3 + ǫv 1, av 3 + v 2 }, {v 4 + av 1, v 3 }, {v 4 + av 1, v 2 } 3 {v 3, v 2, v 1 }, {v 3, v 4, v 1 }, {v 2, v 4, v 1 }, {v 3, v 2, v 4 + av 1 }, Table 10: The Optimal System of subalgebras of the algebra L with a R, ǫ = ±1 4.4 Algebra L with k > 1 2 and k 1 We obtain the algebra spanned on generators v 1 = t, v 2 v 3 = v 4 + (1 k)u x u. The optimal system of subalgebras for this algebra is following 12

13 Dimension Subalgebras 1 {v 3 }, {v 1 }, {v 2 + av 3 }, {v 2 + ǫv 1 }, {v 3 + ǫv 1 }, {v 4 + av 1 }, {v 2 + ǫv 3 + av 1 } 2 {v 2, v 3 }, {v 2, v 1 }, {v 3, v 1 }, {v 4, v 1 }, {v 2 + ǫv 3, v 1 } {v 3 + ǫv 1, v 2 }, {v 2 + ǫv 1, av 2 + v 3 }, {v 4 + av 1, v 2 }, {v 4 + av 1, v 3 } 3 {v 2, v 3, v 1 }, {v 2, v 4, v 1 }, {v 3, v 4, v 1 }, {v 2, v 3, v 4 + av 1 }, Table 11: The Optimal System of subalgebras of the algebra L with a R, ǫ = ±1 4.5 Algebra L with k = 1 2 We obtain the algebra spanned on generators v 1 = t, v 2 v 3 = v 4 x u u. The optimal system of subalgebras for this algebra is following Dimension Subalgebras 1 {v 2 }, {v 1 }, {v 3 + av 2 }, {v 3 + ǫv 1 }, {v 2 + ǫv 1 }, {v 4 + av 1 }, {v 3 + ǫv 2 + av 1 } 2 {v 3, v 2 }, {v 3, v 1 }, {v 2, v 1 }, {v 4, v 1 }, {v 3 + ǫv 2, v 1 } {v 2 + ǫv 1, v 3 }, {v 3 + ǫv 1, av 3 + v 2 }, {v 4 + av 1, v 3 }, {v 4 + av 1, v 2 } 3 {v 3, v 2, v 1 }, {v 3, v 4, v 1 }, {v 2, v 4, v 1 }, {v 3, v 2, v 4 + av 1 }, Table 12: The Optimal System of subalgebras of the algebra L with a R, ǫ = ±1 13

14 5 Conclusion We classified the Lie algebra L under parameter k. For each of meanings of parameter k we found optimal system of subalgebras. We showed that structer of this algebra strongly depend on parameter and optimal systems of subalgebras in each case are different. In [8] was not considered depending References [1] L. A. Bordag, On option-valuation in illiquid markets: invariant solutions to a nonlinear model, Proceedings of Workshop on Mathematical Control Theory and Finance, Lisbon, April, (2007), pp. 1-18, WMCT&F2007, Springer, [2] R. Frey and P. Patie, Risk management for derivatives in illiquid markets: a simulation study, in Advances in Finance and Stochastics ; (K. Sandmann and P. Schönbucher, eds), Springer, Berlin, [3] Nail H. Ibragimov, Elementary Lie Group Analysis and Ordinary Differential Equations, John Wiley&Sons, Chichester, USA, [4] Sophus Lie, Vorlesungen über Differentialgleichungen mit bekannten infinitesimalen Transformationen, Teubner, Leipzig, 1912 (German). [5] P. J. Olver, Application of Lie groups to differential equations, Springer- Verlag, New York, USA, [6] J. Patera and P. Winternitz Subalgebras of real three- and four- dimencional Lie algebras, J. Math. Phys., Vol. 18, No 7 (1977) pp [7] J. Patera and P. Winternitz Continuous subgroups of the fundamental groups of physics. I. General method and Poincaré group, J. Math. Phys., Vol. 16, No 8 (1975) pp [8] Yang Xuan-Liu, Zhang Shun-Li and Qu Chang-Zheng Symmetry Breaking for Black-Scholes Equations, Theor. Phys. (Beijling, China), 47 (2007) pp