On the Generation of Inﬁnitely Many Conservation Laws of the Black-Scholes Equation

: Construction of conservation laws of differential equations is an essential part of the mathematical study of differential equations. In this paper we derive, using two approaches, general formulas for ﬁnding conservation laws of the Black-Scholes equation. In one approach, we exploit nonlinear self-adjointness and Lie point symmetries of the equation, while in the other approach we use the multiplier method. We present illustrative examples and also show how every solution of the Black-Scholes equation leads to a conservation law of the same equation.


Introduction
An important study of mathematical models described by differential equations concerns construction of the inherent conservation laws of equations. This is because of the many uses of conservation laws, which include characterisation of the conserved physical quantities of the phenomenon being modelled. In some cases, conservation laws are used to investigate integrability, existence, uniqueness, and stability of solutions of differential equations. In the case of partial differential equations (PDEs), conservation laws can also be used to search for potential symmetries, which in turn lead to new solutions of the equations via admitted nonlocal symmetries. Therefore, construction of conservation laws of the Black-Scholes equation [1], arguably the most famous equation in financial mathematics, represents an important aspect of the study of the Black-Scholes market.
Hashemi [24] also uses Lie point symmetries of the Black-Scholes equation to compute conservation laws of the equation via Ibragimov's new conservation theorem [4,14].
In this paper, we augment the work by Edelstein and Govinder [23] and Hashemi [24]. We use two methods to construct general formulas for finding conservation laws of the Black-Scholes equation. In the first method, we employ the general conservation theorem by Ibragimov [4,14] by means of which conservation laws for a system of equations consisting of the given system and its adjoint are obtained. The second method used is the direct method proposed by Anco and Bluman in 1996 [3,5]. This method essentially reduces the construction of conservation laws to solving a system of linear determining equations similar to that for finding Lie point symmetries. An explicit formula is then derived which yields a conservation law for each solution of the determining system. Using this method we characterise conservation laws of the Black-Scholes equation in terms of solutions of the associated adjoint equation. Furthermore, we construct a mapping between the Black-Scholes equation and the associated adjoint equation so that every solution of the Black-Scholes equation yields a conservation law of the equation. Mathematica [25] is used to perform all the calculations reported in this paper.
In its simplest form, the Black-Scholes equation is a (1 + 1) linear parabolic equation, where u = u(x, t) is the fair option price depending on the current value of the underlying asset x and time t. The parameters σ and r are the market volatility of the underlying asset price and the interest rate, respectively. The paper is organised as follows. In Section 2, we present relevant preliminaries. In Section 3, we exploit nonlinear self-adjointness of the Black-Scholes equation to derive a general formula for constructing conservation laws of the quation. In Section 4, we derive a general formula for constructing conservation laws of the Black-Scholes via the direct method. We provide illustrative examples in Section 5. An equivalence transformation between the Black-Scholes equation and its adjoint equation is derived in Section 6. Finally, we give concluding remarks in Section 7.

Preliminaries
Consider a system of m PDEs of r-th order where x = (x 1 , . . . , x n ) is an independent variable set and u = (u 1 , . . . , u m ) is a dependent variable set, with u (i) denoting all i-th x derivatives of u. The summation convention for repeated indices is assumed unless otherwise stated. The formal Lagrangian, introduced in [4], associated to the system of Equations (2), is given by the expression where The system of adjoint equations to (2) is defined by where δ/δu α is the Euler operator and D i is the total derivative operator with respect to x i defined by Definition 1. The system of m differential Equations (2) is said to be nonlinearly self-adjoint if the adjoint Equations (4) are satisfied for all solutions u of the original system (2) upon a substitution We are now able to obtain the explicit formulas for conservation laws of any nonlinearly self-adjoint Equation (2) that admits symmetries. We shall, however, omit a discussion of Lie symmetry analysis as this is well-documented in many standard books [16][17][18][19][20][21][22].

Theorem 1. Every infinitesimal Lie point, Lie-Bäcklund and non-local symmetry
of the system (2) leads to a conservation law that is constructed by the formula where W α = η α − ξ i u α j and L is the formal Lagrangian (3), written in the symmetric form in the mixed derivatives.
An alternate method for constructing conservation laws, that circumvents Noether's theorem is the direct/multiplier method [5]. In this method we exploit the fact that every admitted nontrivial conservation law arises from multipliers on the Equation (2). Multipliers for the PDE system (2) are a set of functions {Λ β (x, u, u (1) , . . .)} satisfying.
where (11) holds identically for arbitrary function u(x 1 , x 2 , . . . , x n ). It follows therefore that the conservation law D i C i = 0 holds for all solutions of the system of PDEs (2).
The set of Equation (12) are the linear determining equations for multipliers associated with conservation laws admitted by the given system (2). Once the multipliers are computed, the conserved vectors are derived systematically using (11).

Conservation Laws of the Black-Scholes Equation via Nonlinear Self-Adjointness
The Black-Scholes equation admits 6 + ∞ Lie point symmetries [26] and, being a linear equation, is nonlinearly self-adjoint [15]. The equation is therefore amenable to Ibragimov's method for constructing conservation laws. Every symmetry of the Black-Scholes equation gives rise to a conservation law of the equation. We follow the theory outlined in Section 2 to construct the adjoint equation of the Black-Scholes Equation (1). According to (3) the formal Lagrangian for the Equation (1) is from which we derive the adjoint equation of the Black-Scholes equation as defined in (4): where The Black-Scholes equation, being a linear equation, is nonlinearly self-adjoint [15].
where ϕ is any solution of the adjoint Equation (14) obviously solves the adjoint equation for any solution of the Black-Scholes equation. According to Theorem 1 this leads to the following result: Proposition 1. Every infinitesimal symmetry generator of the Black-Scholes equation gives rise to a conservation law, where and v is any solution of the adjoint Equation (14).

Conservation Laws of the Black-Scholes Equation via the Direct Method
We will exploit the well-known fact that for any linear PDE system, each solution of its adjoint system yields a conservation law of the system [21]. According to (12) a multiplier Λ of a conservation law of the Black-Scholes equation satisfies the equation If we take Λ to be a first-order differential function Λ(x, t, u, u x , u t ), it is easy to show that Λ is a multiplier of the Black-Scholes Equation (1)

if and only if
where v is any solution of the adjoint Equation (14). This is arrived at by expanding Equation (19) and solving the resulting set of determining equations for the multiplier.
Furthermore, it is not hard to show that if Λ = v(x, t) is a solution of the adjoint Equation (14), then where f = f (t, x) and g = g(t, x) are any functions satisfying The leads to the following result: is a conserved vector of the Black-Scholes equation provided that v is a solution of the adjoint Equation (14).
Therefore, any solution of the adjoint Equation (14) leads to a conservation law of the Black-Scholes equations.

Illustrative Examples
In this section we provide examples of conservation laws of the Black-Scholes equation constructed via Propositions 1 and 2. Before we do this, however, we present results of basic Lie symmetry analysis of the Black-Scholes equation and the adjoint equation. We present the admitted Lie point symmetries and associated invariant solutions.

Lie Point Symmetries and Invariant Solutions of the Equations (1) and (14)
Lie point symmetries of the Black-Scholes Equation (1) and the associated adjoint equation are easily obtained. Using Program Lie [27], for example, we determine that symmetries of the Black-Scholes equation are where φ is any solution of the Black-Scholes equation. Similarly, Lie point symmetries of the associated adjoint Equation (14) are where ϕ is any solution of the adjoint Equation (14). Using Lie point symmetries of the Black-Scholes equation and those of the associated adjoint equation, we construct invariant solutions of these two equations via the usual routine [16,17,19,22]. These solutions are given in Tables 1 and 2, respectively.  (1).  (14).

Construction of Conservation Laws of the Black-Scholes Equation via Proposition 1
Consider symmetries (24) of the Black-Scholes equation.

Construction of Conservation Laws of the Black-Scholes Equation via Proposition 2
Consider invariant solutions of the adjoint equations given in Table 2.

Every Solution of the Black-Scholes Equation Gives Rise to a Conservation Law of the Equation
In both Propositions 1 and 2, the constructed conservation laws involve solutions of the adjoint equation. We can avoid making reference to the adjoint equation in the constructed conservation laws if we exploit the equivalence between the Black-Scholes equation and the associated adjoint equation, which is done via an equivalence transformation. , Ω = D + σ 2 2 2 σ 2 (41) is a solution of the adjoint Equation (14).

Concluding Remarks
Two approaches have been employed in this paper to establish general characterizations of conservation laws of the Black-Scholes equation. In one approach the self-adjointness of the Black-Scholes equation was exploited following a method due to Ibragimov [4], while in the other approach the direct method, proposed by Anco and Bluman [3], was used. We have provided illustrations of how infinitely many conservation laws of the Black-Scholes equation may be determined easily from the derived general characterizations of conservation laws of the equation. Furthermore, we have constructed an equivalence transformation between the Black-Scholes equation and its adjoint equation, which provides a correspondence between every solution of the Black-Scholes equation and a conservation law of the equation.