Contact Symmetries of a Model in Optimal Investment Theory

: It is generally known that Lie symmetries of differential equations can lead to a reduction of the governing equation(s), lead to exact solutions of these equations and, in the best case scenario, lead to a linearization of the original equation. In this paper, we consider a model from optimal investment theory where we show the governing equation possesses an extensive contact symmetry and, through this, we show it is linearizable. Several exact solutions are provided including a solution to a particular terminal value problem.


Introduction
Nonlinear partial differential equations (NLPDEs) play an integral part in describing the world around us. They can be found in the fields of nonlinear diffusion, wave propagation, Mathematical Biology, ray optics, solid mechanics, and financial mathematics to name just a few (see, for example, refs. [1][2][3][4][5] and the references within). However, obtaining exact solutions of these equations is usually a difficult task and techniques for obtaining solutions is a current area of research. One popular technique are symmetry methods probably due to the fact that the method is rather algorithmic and thus computer algebra systems such as Maple and Mathematica can be used. Symmetry methods have been extensively used in a number of fields, and we refer the reader to the books by Arrigo [6], Bluman et al. [7,8], Bordag [9], Cantwell [10], Cherniha et al. [11], and Olver [12] .
In this paper, we are interested in a model from optimal investment theory. Consider an investment portfolio consisting of n + 1 assets. Let the first asset be a bond and the next n assets be stocks, all of which are traded continuously. In the simplest case where n = 1, the value of the portfolio, u(t, x), for time t and investment amount x, one model presented by see Yong [13] is the NPDE where θ = b−r 2σ and the variables r, b, σ represent the interest rate, appreciation rate, and volatility, respectively, and are assumed constant with σ > 0 and b − r > 0.
A classical symmetry analysis was performed by Yang and Xu [14] who were able to show that (2) admitted the symmetry generator where T = c 1 , where c i , i = 1...5 are arbitrary constants. Through their analysis, they were able to exploit several of these symmetries to obtain a number of reductions and, in some cases, construct exact solutions. It is natural to ask whether (2) admits symmetries that are more general than Lie point symmetries. In this paper, we consider contact symmetries of (2), and we will show that, in fact, (2) admits a rather large contact symmetry which leads to its linearization. In addition to recovering known solutions, we will obtain new exact solutions. We also solve a particular terminal value problem.

Contact Symmetries
In this section, we construct contact symmetries of (2). If we denote this original NLPDE by ∆ so then contact symmetries are given by where the infinitesimal generator Γ is where U = U(t, x, u, u t , u x ); the first and second extensions of the generator are where the operators D t and D x are This leads to the set of determining equations Although somewhat a laborious calculation, we find the solution of (10) to be where c i , i = 1...6 are arbitrary constants and the function F(t, u x ) satisfies Equation (12) is linear and possesses an infinite number of solutions, which means that there are an infinite number of symmetries to (2). Furthermore, since there is a particular function F in (11) that satisfies a linear PDE, this suggests that (2) can be transformed to a linear PDE (Bluman and Kumei [15]). As both t and u x are independent variables in (12), this indicates that maybe we should use these as new independent variables.

A Linearization
Since the symmetry obtained in the last section contains an arbitrary function that satisfies a linear PDE, this suggests that the original PDE is linearizable. Introducing the new variables where U = U(T, X), derivatives transform as and (2) becomes which is exactly (12). Interestingly enough, (15) looks remarkably similar to the Black-Scholes-Merton equation [16,17], which is known to be mappable to the linear heat equation ( [16,18]), so it should come as no surprise that the same is true for (15). However, it would make the contact transformation (13) more complicated and, thus, we will not pursue this line any further.
In the next section, we obtain simple exact solutions of (2) in addition to exploiting Lie symmetries of the linearized Equation (15) to obtain additional solutions.

Exact Solutions
We have shown that the nonlinear PDE can be transformed to the linear PDE via the transformation t = T, x = U X , u = XU X − U.
We are now in a position to obtain a number of exact solutions to (16). For example, (17) admits separable solutions of the form where F and G satisfies where k is a separation constant. Equation (20a) is easily solved giving for some arbitrary constant F 0 . Equation (20b) possesses solutions of the form X m , where m is a solution of For example, if k = sm (where s is some constant), then the solution of (22) is leading to the exact solution Passing (24) through the transformation (18) (resetting the constants c 1 and c 2 ) gives which recovers the exact solution presented by Yang and Xu [14] by choosing c 1 = 0, c 2 = 1 and s = 1.
As a second example, if k = r, then, from (22), we obtain m = 1, r/θ ( = 1) and we obtain the solution to (20b) as G = g 1 X + g 2 X r/θ (26) and the general solution to (17) Passing (27) of (16), which we believe to be new. In the cases where r = θ, the solution of (20b) is and the general solution to (17) Passing (30) through the transformation (18) leads to the exact solution which was given in [14] in the case of r = θ. Of course, other choices of m that satisfy (22) would lead to exact solutions of (16) which we will not pursue here.
We also note that symmetries of (17) can be used to generate new solutions of (17), which could lead to new solutions of (16). It is well known (Broadbridge and Arrigo [19]) that, if (17) possesses symmetries with the generator where Σ, Φ, and Ξ have some particular forms and Q satisfies the original PDE (17), then, if one has one seed solution, say U = U 0 (T, X), then additional solutions can be obtained from For example, (17) admits the symmetry generator (32) where Σ, Φ, and Ξ are given by where c i , i = 1...6 are arbitrary constants.
One particularly simple solution of (17) is From (33), we obtain the solution Passing (36) through the transformation (18) leads to a solution that is parametric in its nature (which we do not list here); however, setting c 3 = 0, we obtain (omitting translational constants) which we believe is new. Of course, other seed solutions could lead to an abundance of exact solutions to (17) which, in turn, would lead to exact solutions to (16).

A Particular Terminal Value Problem
A particular problem of interest is one that is given in Koleva and Vulkow [20], which is to solve (2) subject to the terminal condition Here, we introduce a slight variation of (18) Under this transformation, the PDE (2) still transforms to (17); however, the terminal condition (38) turns into the initial condition U(X, 0) = X − X ln X. (40) At this point, we exploit the symmetries obtained in the previous section. With these symmetries, we associate an invariant surface condition where Σ, Φ and Ξ are given in (34). From (40) and (17), we obtain the initial conditions Requiring that (40) and (42) satisfy (41) on the boundary T = 0 gives where c 1 and c 2 are arbitrary. Here, we choose c 1 = 1, c 2 = 0, leading to the invariant surface condition This is easily solved giving Imposing the initial condition (40) on the solution (45) gives and, if we let λ = ln X, we obtain and, from (45), we obtain U = X(1 − ln X + (θ − r)T)e rT (48) and one can indeed verify that (48) does satisfy (17) and (42). As the final step, we pass (48) through the transformation (39), leading to which satisfies the original PDE (2) and the terminal condition (38).

Conclusions
It is well known that classical Lie symmetries can lead to a reduction of a given PDE and sometimes lead to exact solutions of the equation. The best case scenario, albeit rare, indicates that the original equation is linearizable. In this paper, we constructed the contact symmetries of a model from an optimal investment theory, which led to a linearization of the given PDE. Several exact solutions were obtained. The symmetries of the linearized equation were also considered where an additional exact solution was obtained in addition to solving a particular terminal value problem.