Next Article in Journal
Natural Compounds as DPP-4 Inhibitors: 3D-Similarity Search, ADME Toxicity, and Molecular Docking Approaches
Next Article in Special Issue
Some Explicit Properties of Frobenius–Euler–Genocchi Polynomials with Applications in Computer Modeling
Previous Article in Journal
Attribute Network Representation Learning with Dual Autoencoders
Previous Article in Special Issue
Beyond the Beta Integral Method: Transformation Formulas for Hypergeometric Functions via Meijer’s G Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Symmetry Analysis of a Model of Option Pricing and Hedging

by
Sergey M. Sitnik
1,*,
Khristofor V. Yadrikhinskiy
2 and
Vladimir E. Fedorov
2,3
1
Department of Applied Mathematics and Computer Modeling, Institute of Engineering and Digital Technologies, Belgorod State National Research University (BelGU), Pobedy Street, 85, 308015 Belgorod, Russia
2
Yakut Branch of Far Eastern Center for Mathematical Research, North Eastern Federal University, 58, Belinskiy St., 677000 Yakutsk, Russia
3
Mathematical Analysis Department, Mathematics Faculty, Chelyabinsk State University, 129, Kashirin Brothers St., 454001 Chelyabinsk, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(9), 1841; https://doi.org/10.3390/sym14091841
Submission received: 15 August 2022 / Revised: 29 August 2022 / Accepted: 1 September 2022 / Published: 5 September 2022

Abstract

:
The Guéant and Pu model of option pricing and hedging, which takes into account transaction costs, and the impact of operations on the market is studied by group analysis methods. The infinite-dimensional continuous group of equivalence transforms of the model is found. It is applied to get the group classification of the model under consideration. In addition to the general case, the classification contains three specifications of a free element in the equation, which correspond to models with groups of symmetries of a special kind. Optimal systems of subalgebras for some concrete models from the obtained classification are derived and used for the calculation of according invariant submodels.

1. Introduction

The classical Black–Scholes model [1,2] of option pricing dynamics is based on the perfect market hypothesis. Under this hypothesis, there are no execution costs and market participants use only the prevailing market prices and cannot influence the prices by their operations. The Black–Scholes model gives useful results when the underlying asset is liquid and the transaction amount is not too large for the market. However, the perfect market hypothesis contradicts to the market practice in many aspects, this fact makes the classical model too limited in application.
Last decades many researchers actively studied improvements of the classical Black–Scholes model, which would take into account the market illiquidity and the impact of transactions on prices. One may see works of Magill and Constantinides [3], Kyle [4], Leland [5], Cvitanić and Karatzas [6], Barles and Soner [7], Grossman [8], Platen and Schweizer [9], Sircar and Papanicolaou [10], Schönbucher and Wilmott [11], Bank and Baum [12], Çetin, Jarrow and Protter ([13] Section 4), Çetin and Rogers ([14] Section 6), Rogers and Singh [15].
New models proposed in these works have been investigated by many researchers both numerically and analytically. The work of Ibragimov and Gazizov [16] contains the first analytical investigation of the Black–Scholes equation by the group analysis methods [17,18].
Note that differential equations play a major role in solving problems of modeling processes and phenomena of the surrounding world and their importance does not decrease over the years (see, e. g., recent papers [19,20,21]). The group analysis (or symmetry analysis) of differential equations is one of the few theories that provide methods for finding exact solutions to nonlinear differential equations and systems of equations of wide classes. Since the middle of the 20th century, a huge number of results have been obtained on the group structure and exact solutions of many equations and systems of equations encountered in describing the dynamics of various physical processes, especially in gas dynamics, elasticity theory, etc. (see [17,18,22,23] and the bibliographies therein).
Let us determine that, following the classical works [17,18], we call a symmetry group of a differential equation a one-parameter group of transformations of independent and dependent variables. After such transformations, the differential equation under study does not change its form in the new variables. Each such group uniquely corresponds to a first-order differential operator, and these operators form a Lie algebra of the differential equation under consideration. The commutator is the Lie multiplication in such algebra. Therefore, when talking about the symmetry groups of a differential equation or about its group structure, at the same time we are talking about the Lie algebra of operators of the mentioned kind for this equation. Strict definitions of the listed objects are very cumbersome, they can be found by the reader in monographs [17,18] or in any other monograph on group analysis of differential equations.
In the last decade, in the works of Bordag [24,25], of Dyshaev and Fedorov [26,27,28,29,30,31,32] group properties of various nonlinear Black–Scholes type models were studied, and their invariant solutions and submodels were calculated. In the papers of Dyshaev and Fedorov, group classifications for various classes of nonlinear Black–Scholes type models were obtained.
Guéant and Pu in [33,34] carried out an analysis of options pricing taking into account transaction costs and the impact of operations on the market under the next assumptions:
(1)
the risk-free rate r, the absolute risk aversion parameter γ and the volatility σ are constant;
(2)
the process of market trading volume V t is deterministic, non-negative, and bounded;
(3)
there exists a maximum degree of participation ρ m , i.e., processes ν are such that | ν t | ρ m V t almost everywhere;
(4)
the number of shares in the hedged portfolio is q t = q 0 + 0 t ν s d s ;
(5)
the price process is modeled by the stochastic differential equation d S t = μ d t + σ d W t , where μ is the expected return of the underlying asset;
(6)
to model execution costs, a continuous, non-negative, even, strictly convex function L : R R + is used, which is increasing on R + , L ( 0 ) = 0 , and coercive, i.e., lim ρ + L ( ρ ) / ρ = + ;
(7)
the dynamics of the account X is described by the equation d X t = r X t d t ν t S t d t V t L ( ν t / V t ) d t .
As result, Guéant and Pu derived a differential equation
θ t = r θ + ( μ r S ) q μ θ S 1 2 σ 2 θ S S 1 2 γ σ 2 e r ( T t ) ( θ S q ) 2 + V t H ( θ q ) ,
where H ( p ) = sup | ρ | ρ m [ p ρ L ( ρ ) ] . It is a model of the dynamics of the indifference price θ ( t , S , q ) for a call option. This price depends on time t, the price of the underlying asset S, and the number of shares in the hedged portfolio q.
In [35], the group classification of the Guéant—Pu model (1) with a constant market trading volume V t is obtained, and for all specifications of the free element H from the classification optimal systems of subalgebras of the Lie algebra is found, invariant solutions and submodels for subalgebras from the optimal systems are derived.
In the present paper, the Guéant—Pu model
θ t = r θ + ( μ r S ) q μ θ S 1 2 σ 2 θ S S 1 2 γ σ 2 e r ( T t ) ( θ S q ) 2 + F ( t , θ q )
is investigated. Here a free element F depends on t and θ q , i.e., the market trading volume V t may depend on t, in contrast to the model, which is considered in [35]. In the Section 2, the continuous group of equivalence transformations of Equation (2) is calculated. It correponds to an infinite-dimensional Lie algebra of the equation with three-dimensional finite part and with two basis operators, which coefficients are defined by two arbitrary functions of t and by their derivatives. In the Section 3 the search of the symmetry groups for general Equation (2) started. The equivalence transformations are used in the Section 4 for the search of the specifications of the free element F, such that F θ q θ q 0 , which corresponds to equations of form (2) with different Lie algebras. The obtained theorem on group classification is formulated in the Section 5. In the Section 6, optimal systems of subalgebras are found for the Lie algebra of model (2) with a general function F ( t , θ q ) and with a specification F = e r t Φ ( θ q ) , which were obtained in the group classification. For every subalgebra from the optimal system, the invariant submodel of the Guéant—Pu model is calculated, if it exists.

2. Continuous Groups of Equivalence Transformations

Consider the Gueant—Pu Equation (2), where θ = θ ( t , S , q ) , F ( t , θ q ) is a free element. Assume that r γ σ μ 0 , T > 0 . For the search of continouos equivalence transformations groups of Equation (2) we will consider the function F and all its derivatives as additional variables. Generators of such groups have a form
Y = τ t + ξ S + α q + η θ + ζ F ,
where τ , ξ , α , η depend on t , S , q , θ , and ζ depends on t , S , q , θ , F , θ t , θ S , θ q . Hereafter β : = β is the partial derivative with respect to a variable β . Equation (2) with a new variable F we will consider in the system with additional equations
θ t = r θ + ( μ r S ) q μ θ S 1 2 σ 2 θ S S 1 2 γ σ 2 e r ( T t ) ( θ S q ) 2 + F ,
F S = 0 , F q = 0 , F θ = 0 , F θ t = 0 , F θ S = 0 ,
which show the dependence of F on t and θ q only. System (3), (4) is considered as a manifold M in the expanded space of the corresponding variables. Let us act by the prolongated operator [17,18]
Y 2 = Y + η t θ t + η S θ S + η q θ q + η S S θ S S + ζ t F t + ζ S F S + ζ q F q + ζ θ F θ + + ζ θ t F θ t + ζ θ S F θ S + ζ θ q F θ q
on the both sides of Equation (3). In order to use the geometric invariance criterion [17], we restrict the result on the manifold M and obtain the equation
η t r η ( μ r S ) α + r q ξ + μ η S + γ σ 2 e r ( T t ) ( θ S q ) η S α r 2 ( θ S q ) τ ζ + 1 2 σ 2 η S S | M = η t r η + r S α + r q ξ + ( μ + γ σ 2 e r ( T t ) ( θ S q ) ) ( η S α ) r 2 γ σ 2 e r ( T t ) ( θ S q ) 2 τ ζ + 1 2 σ 2 η S S | M = 0 .
The coefficients of the prolongated operator Y 2 are calculated using the total derivatives operators
D t = t + θ t θ + , D S = S + θ S θ + θ S S θ S + , D q = q + θ q θ + , D ˜ t = t + F t F + , D ˜ S = S + F S F + , D ˜ q = q + F q F + , D ˜ θ = θ + F θ F + , D ˜ θ t = θ t + F θ t F + , D ˜ θ S = θ S + F θ S F +
and the prolongation formulas (see detail in [17,18])
η t = D t η θ t D t τ θ S D t ξ θ q D t α , η S = D S η θ t D S τ θ S D S ξ θ q D S α , η q = D q η θ t D q τ θ S D q ξ θ q D q α , η S S = D S η S θ S t D S τ θ S S D S ξ θ S q D S α , ζ S = D ˜ S ζ F t D ˜ S τ F S D ˜ S ξ F q D ˜ S α F θ D ˜ S η F θ t D ˜ S η t F θ S D ˜ S η S F θ q D ˜ S η q , ζ q = D ˜ q ζ F t D ˜ q τ F S D ˜ q ξ F q D ˜ q α F θ D ˜ q η F θ t D ˜ q η t F θ S D ˜ q η S F θ q D ˜ q η q , ζ θ = D ˜ θ ζ F t D ˜ θ τ F S D ˜ θ ξ F q D ˜ θ α F θ D ˜ θ η F θ t D ˜ θ η t F θ S D ˜ θ η S F θ q D ˜ θ η q , ζ θ t = D ˜ θ t ζ F t D ˜ θ t τ F S D ˜ θ t ξ F q D ˜ θ t α F θ D ˜ θ t η F θ t D ˜ θ t η t F θ S D ˜ θ t η S F θ q D ˜ θ t η q , ζ θ S = D ˜ θ S ζ F t D ˜ θ S τ F S D ˜ θ S ξ F q D ˜ θ S α F θ D ˜ θ S η F θ t D ˜ θ S η t F θ S D ˜ θ S η S F θ q D ˜ θ S η q .
The result of the action of Y 2 on Equation (4) after restricting on the manifold M gives
ζ S | M = ζ S F t τ S F θ q η S q | M = ζ S F t τ S F θ q ( η S q + θ q η S θ θ t ( τ S q + θ q τ S θ ) θ S ( ξ S q + θ q ξ S θ ) θ q ( α S q + θ q α S θ ) ) | M = 0 , ζ q | M = ζ q F t τ q F θ q η q q | M = ζ q F t τ q F θ q ( η q q + θ q η q θ θ t ( τ q q + θ q τ q θ ) θ S ( ξ q q + θ q ξ q θ ) θ q ( α q q + θ q α q θ ) ) | M = 0 , ζ θ | M = ζ θ F t τ θ F θ q η θ q | M = ζ θ F t τ θ F θ q ( η q θ + θ q η θ θ θ t ( τ q θ + θ q τ θ θ ) θ S ( ξ q θ + θ q ξ θ θ ) θ q ( α q θ + θ q α θ θ ) ) | M = 0 ζ θ t | M = ζ θ t F θ q η θ t q | M = ζ θ t + F θ q ( τ q + θ q τ θ ) | M = 0 , ζ θ S | M = ζ θ S F θ q η θ S q | M = ζ θ S + F θ q ( ξ q + θ q ξ θ ) | M = 0 .
The transition on the manifold M means the substitution for θ t the right-hand side of (3) and vanishing of variables F S , F q , F θ , F θ t , F θ S . It does not change the form of the last two equations in (6). Therefore, the separation of variables F θ q and θ q gives ζ θ t = 0 , ζ θ S = 0 , τ q = 0 , τ θ = 0 , ξ q = 0 , ξ θ = 0 .
We substitute the prolongation formulas into Equation (5) and after the transition on M equate the result to zero:
η t θ S ξ t θ q α t r η + r S α + r q ξ r 2 γ σ 2 e r ( T t ) ( θ S q ) 2 τ ζ + + ( μ + γ σ 2 e r ( T t ) ( θ S q ) ) ( η S + θ S η θ θ S ξ S θ q ( α S + θ S α θ ) α ) + + 1 2 σ 2 ( η S S + 2 θ S η S θ + θ S 2 η θ θ 2 θ S t τ S + θ S S ( η θ θ q α θ 2 ξ S ) 2 θ S q ( α S + θ S α θ ) θ S ξ S S θ q ( α S S + 2 θ S α S θ + θ S 2 α θ θ ) ) + + r θ + ( μ r S ) q μ θ S 1 2 σ 2 θ S S 1 2 γ σ 2 e r ( T t ) ( θ S q ) 2 + F × × η θ τ t θ q α θ ( μ + γ σ 2 e r ( T t ) ( θ S q ) ) τ S σ 2 2 τ S S = 0 .
Since all the functions in (7) do not depend on θ S t , θ S q , equate to zero the coefficients in Equation (7) at θ S t , θ S q and obtain the equalities τ S = 0 , α S + θ S α θ = 0 ; since α and its derivatives does not depend on θ S , we have α S = α θ = 0 . Thus,
τ S = 0 , τ q = 0 , τ θ = 0 , ξ q = 0 , ξ θ = 0 , α S = 0 , α θ = 0 , ζ θ t = 0 , ζ θ S = 0 .
The first 3 equations in (6) now have the form
ζ S | M = ζ S F θ q ( η S q + θ q η S θ ) = 0 , ζ q | M = ζ q F θ q ( η q q + θ q η q θ θ q α q q ) = 0 , ζ θ | M = ζ θ F θ q ( η q θ + θ q η θ θ ) = 0 .
The separation of the variables F θ q and θ q here gives
η S θ = 0 , η S q = 0 , α q q = η q θ = 0 , η q q = 0 , η θ θ = 0 , ζ S = 0 , ζ q = 0 , ζ θ = 0 , ζ θ t = 0 , ζ θ S = 0 .
By substitution into Equation (7) equalities (8) and (9) we get
η t θ S ξ t θ q α t r η + r S α + r q ξ + ( μ + γ σ 2 e r ( T t ) ( θ S q ) ) ( η S + θ S η θ θ S ξ S α ) r 2 γ σ 2 e r ( T t ) ( θ S q ) 2 τ ζ + 1 2 σ 2 ( η S S + θ S S ( η θ 2 ξ S ) θ S ξ S S ) + + ( η θ τ t ) r θ + ( μ r S ) q μ θ S 1 2 σ 2 θ S S 1 2 γ σ 2 e r ( T t ) ( θ S q ) 2 + F = 0 .
We separate this equation by the variables θ S S , θ S , since ζ does not depend on them in view of (9), and after a reduction we obtain
θ S S : 2 ξ S = τ t , θ S 2 : η θ 2 ξ S r τ + τ t = 0 ,
θ S : γ σ 2 e r ( T t ) ( r q τ + q ξ S α + η S q τ t ) μ ξ S + μ τ t ξ t σ 2 2 ξ S S = 0 ,
1 : η t θ q α t r η + r S α + r q ξ + ( μ γ σ 2 e r ( T t ) q ) ( η S α ) r 2 γ σ 2 e r ( T t ) q 2 τ ζ + σ 2 2 η S S + ( η θ τ t ) r θ + ( μ r S ) q 1 2 γ σ 2 e r ( T t ) q 2 + F = 0 .
From 2 ξ S = τ t in (10) in view of τ S = 0 due to (8) we get ξ S S = 0 . Substitution 2 ξ S = τ t from the first equation to the second one in (10) yields η θ = r τ . Now we differentiate (11) by q and using (8), (9) we get r τ + ξ S α q τ t = 0 . Substitute ξ S = τ t / 2 from (10), then α q = r τ τ t / 2 . Next, by differentiating (11) by S and using (8) and (9), we obtain γ σ 2 e r ( T t ) η S S ξ t S = 0 , or γ σ 2 e r ( T t ) η S S = τ t t / 2 . Therefore, η S S S = 0 . Thus,
ξ S S = 0 , α q = r τ τ t 2 , η θ = r τ , η S S = e r ( t T ) 2 γ σ 2 τ t t , η S S S = 0 .
The differentiation of Equation (12) by S twice with the substitution of the vanishing functions from (8), (9) and (13) gives η t S S r η S S = 0 . The substitution of η S S from (13) here leads to τ t t t = 0 . By the differentiation of Equation (12) by θ q we obtain α t + ζ θ q = 0 . The differentiation by q gives α t q = 0 , i.e., due to (13) r τ t τ t t / 2 = 0 . Since r 0 , this differential equation implies the equality τ t t = 0 , hence τ t = 0 . Then due to (10), (13)
τ t = 0 , α q = r τ , ξ S = 0 , η S S = 0 .
From (8), (9), (14) it follows that τ is a constant, ξ = ξ ( t ) , α = r τ q + A ( t ) . Substituting these equalities into (11) we get γ σ 2 e r ( T t ) ( η S A ( t ) ) ξ t = 0 . Therefore, due to (9), (13), (14) η = r τ θ + B ( t ) q + ( A ( t ) + e r ( t T ) ξ ( t ) / ( γ σ 2 ) ) S + C ( t ) . So,
ξ = ξ ( t ) , α = r τ q + A ( t ) , η = r τ θ + B ( t ) q + A ( t ) + e r ( t T ) ξ ( t ) γ σ 2 S + C ( t ) .
Note that arbitrary functions used to represent the solution of the considered system of equations here and further we denote with capital letters A, B, C, etc. Substituting equalities (15) into (12) and shortening we obtain
( B r B + r ξ ξ ) q + C r C + μ γ σ 2 e r ( t T ) ξ + + A + e r ( t T ) γ σ 2 ξ S A θ q ζ + r τ F = 0 .
The differentiation by q of Equation (16) implies that B r B + r ξ ξ = 0 , hence B ( t ) = ξ ( t ) + D e r t . Next, differentiate by S Equation (16) and obtain
B = ξ + D e r t , A = e r T γ σ 2 t 0 t e r s ξ ( s ) d s + c .
We substitute (17) into (16) and (15) and get
ξ = ξ ( t ) , α = r τ q e r T γ σ 2 t 0 t e r s ξ ( s ) d s + c , η = r τ θ + ( ξ ( t ) + D e r t ) q + e r T γ σ 2 e r t ξ ( t ) t 0 t e r s ξ ( s ) d s S + c S + C ( t ) , ζ = C ( t ) r C ( t ) + μ γ σ 2 e r ( t T ) ξ ( t ) + e r ( t T ) γ σ 2 ξ ( t ) θ q + r τ F .
We formulate this result in the form of a theorem. To do this, we fix alternately one of the integration constants or an arbitrary function, equating the others to zero, and obtain the corresponding basic operator of the resulting Lie algebra. By Y k we denote here these basis operators. If coefficients of an operator are defined by a function and its derivatives, we use this function denotation as the lower index for such operator.
Theorem 1.
The Lie algebra of continuous equivalence transformations for Equation (2) is generated by the operators
Y 1 : = e r t q θ , Y 2 : = q + S θ , Y 3 : = t + r q q + r θ θ + r F F , Y ϕ : = ϕ ( t ) θ + ( ϕ ( t ) r ϕ ( t ) ) F , Y ψ : = ψ ( t ) S e r T γ σ 2 t 0 t e r s ψ ( s ) d s q + + ψ ( t ) q + e r T γ σ 2 e r t ψ ( t ) t 0 t e r s ψ ( s ) d s S θ + + μ γ σ 2 e r ( t T ) ψ ( t ) + e r ( t T ) γ σ 2 ψ ( t ) θ q F .
Solving the Lie equations for the obtained Lie algebras and taking the projections on the variables t, θ q , F we get
Y 1 : θ ¯ q = θ q + a 1 e r t ; Y 3 : t ¯ = t + a 3 , F ¯ = e r a 1 F ; Y ϕ : F ¯ = F r ϕ ( t ) + ϕ ( t ) ; Y ψ : θ ¯ q = θ q + ψ ( t ) , F ¯ = F + μ γ σ 2 e r ( t T ) ψ ( t ) + e r ( t T ) 2 γ σ 2 ψ ( t ) ψ ( t ) + e r ( t T ) γ σ 2 ψ ( t ) θ q .
Remark 1.
A Lie algebra is called principal [17] for Equation (2), if it is admissible for (2) with any specification of F. From (18) it follows that the principal Lie algebra of Equation (2) is generated by Y 2 and by Y ϕ at ϕ ( t ) = e r t . Indeed, for such ϕ the group of transformations, which is generated by Y ϕ , does not change t, θ q and F.
Remark 2.
We see that the Lie algebra of continuous equivalence transformations for Equation (2) is infinite-dimensional, since its operators depend on arbitrary functions ϕ and ψ. Note that such equation with a function F depending on θ q only has a 5-dimensional Lie algebra of continuous equivalence transformations (see [35]). It generated by Y 2 , Y 3 , Y ϕ for ϕ ( t ) 1 , Y ϕ for ϕ ( t ) = e r t and Y ψ at ψ ( t ) 1 .

3. Calculation of the Symmetry Groups in General Case

Our purpose is to obtain the so-called group classification [17] for equation
θ t = r θ + ( μ r S ) q μ θ S σ 2 2 θ S S 1 2 γ σ 2 e r ( T t ) ( θ S q ) 2 + F ( t , θ q ) .
For this aim, firstly we will search generators of the symmetry groups for the equation under general assumptions.
On Equation (19) we act by the second prolongation X 2 = X + η q θ q + η S θ S + η t θ t + η S S θ S S for a generator X = τ t + ξ S + α q + η θ of a continuous group of transformations, where functions τ , ξ , α , η depend on t, S, q, θ . So,
η t r η + r S α + r q ξ + ( μ + γ σ 2 e r ( T t ) ( θ S q ) ) ( η S α ) r 2 γ σ 2 e r ( T t ) ( θ S q ) 2 τ + σ 2 2 η S S F t τ F θ q η q | M = 0 .
After the substitution into (20) of the prolongation formulas and the restriction on the manifold M , using the Equation (19) for θ t , we obtain
r θ + ( μ r S ) q μ θ S σ 2 2 θ S S 1 2 γ σ 2 e r ( T t ) ( θ S q ) 2 + F × × η θ τ t θ S ξ θ θ q α θ ( μ + γ σ 2 e r ( T t ) ( θ S q ) ) ( τ S + θ S τ θ ) + + F θ q ( τ q + θ q τ θ ) σ 2 2 ( θ S S τ θ + θ S τ S S + 2 θ S τ S θ + θ S 2 τ θ θ ) r θ + ( μ r S ) q μ θ S σ 2 2 θ S S 1 2 γ σ 2 e r ( T t ) ( θ S q ) 2 + F 2 τ θ + + η t θ S ξ t θ q α t r η + r S α + r q ξ + ( μ + γ σ 2 e r ( T t ) ( θ S q ) ) ( η S + θ S η θ θ S ( ξ S + θ S ξ θ ) θ q ( α S + θ S α θ ) α ) r 2 γ σ 2 e r ( T t ) ( θ S q ) 2 τ F t τ F θ q ( η q + θ q η θ θ S ( ξ q + θ q ξ θ ) θ q ( α q + θ q α θ ) ) + + σ 2 2 ( η S S + 2 θ S η S θ + θ S 2 η θ θ 2 θ t S ( τ S + θ S τ θ ) + θ S S ( η θ θ q α θ 2 ξ S 3 θ S ξ θ ) 2 θ S q ( α S + θ S α θ ) θ S ( ξ S S + 2 θ S ξ S θ + θ S 2 ξ θ θ ) θ q ( α S S + 2 θ S α S θ + θ S 2 α θ θ ) ) = 0 .
The differentiation of this equation by the variables θ S q and θ t S leads to the equations τ S = 0 , τ θ = 0 , α S = 0 , α θ = 0 .
Equating the coefficient at θ S S in (21) to zero, obtain ξ θ = 0 , τ t 2 ξ S F θ q τ q = 0 , and using the equality τ S = 0 we get ξ S S = 0 . Therefore,
τ S = 0 , τ θ = 0 , α S = 0 , α θ = 0 , ξ θ = 0 , ξ S S = 0 , τ t 2 ξ S F θ q τ q = 0 .
Applying these equalities in (21) we obtain the equality
r θ + ( μ r S ) q μ θ S 1 2 γ σ 2 e r ( T t ) ( θ S q ) 2 + F ) ( η θ τ t + F θ q τ q + η t θ S ξ t θ q α t r η + ( μ + γ σ 2 e r ( T t ) ( θ S q ) ) ( η S + θ S η θ θ S ξ S α ) + + r S α + r q ξ r 2 γ σ 2 e r ( T t ) ( θ S q ) 2 τ F t τ F θ q ( η q + θ q η θ θ S ξ q θ q α q ) + + 1 2 σ 2 ( η S S + 2 θ S η S θ + θ S 2 η θ θ ) = 0 .
We separate this equation by the variable θ S taking into account the last equation from (22) and get the equations
θ S 2 : γ e r ( T t ) ( η θ r τ ) + η θ θ = 0 ,
θ S : γ σ 2 e r ( T t ) ( q ξ S + η S α + r q τ ) + F θ q ξ q + σ 2 η S θ ξ t + 2 μ ξ S = 0 ,
1 : r θ + ( μ r S ) q γ σ 2 2 e r ( T t ) q 2 + F ( η θ τ t + F θ q τ q ) + η t θ q α t r η + + r S α + r q ξ + ( μ γ σ 2 e r ( T t ) q ) ( η S α ) r 2 γ σ 2 e r ( T t ) q 2 τ F t τ F θ q ( η q + θ q η θ θ q α q ) + σ 2 2 η S S = 0 .
From (22) it follows that ξ = A ( t , q ) S + B ( t , q ) , Equation (23) implies that η θ = r τ + C 0 ( t , S , q ) e γ e r ( T t ) θ . Therefore,
ξ = A ( t , q ) S + B ( t , q ) , η = r θ τ + C ( t , S , q ) e γ e r ( T t ) θ + D ( t , S , q ) .
Substitute these equalities into (24), (25) and get
γ σ 2 e r ( T t ) ( q A + D S α + r q τ ) A t S B t + F θ q ( A q S + B q ) + 2 μ A = 0 ,
r θ + ( μ r S ) q γ σ 2 2 e r ( T t ) q 2 + F ( r τ γ e r ( T t ) C e γ e r ( T t ) θ 2 A ) + + r θ τ t + C t e γ e r ( T t ) θ + r γ e r ( T t ) θ C e γ e r ( T t ) θ + D t r 2 θ τ r C e γ e r ( T t ) θ r D θ q α t + r S α + r q ( A S + B ) + ( μ γ σ 2 e r ( T t ) q ) ( C S e γ e r ( T t ) θ + D S α ) r 2 γ σ 2 e r ( T t ) q 2 τ F t τ + σ 2 2 ( C S S e γ e r ( T t ) θ + D S S ) F θ q ( r θ τ q + C q e γ e r ( T t ) θ + D q + θ q ( r τ γ e r ( T t ) C e γ e r ( T t ) θ ) θ q α q ) = 0 .
In the last equation the variable θ is present explicitly, after the reduction of similar terms the equation has a form a + b e q θ = 0 , q 0 . Hence a = b = 0 and we have the equations
a = ( μ r S ) q γ σ 2 2 e r ( T t ) q 2 + F ( r τ 2 A ) + D t r D θ q α t + r S α + r q ( A S + B ) + ( μ γ σ 2 e r ( T t ) q ) ( D S α ) r 2 γ σ 2 e r ( T t ) q 2 τ F t τ + σ 2 2 D S S F θ q ( D q + r θ q τ θ q α q ) = 0 ,
b = ( μ r S ) q γ σ 2 2 e r ( T t ) q 2 + F γ e r ( T t ) C + C t r C + + ( μ γ σ 2 e r ( T t ) q ) C S + σ 2 2 C S S F θ q ( C q γ e r ( T t ) C θ q ) = 0 .

4. Calculation of the Group Classification in the Case F θ q θ q 0

Let us continue the calculations using the assumption F θ q θ q 0 . Differentiating the last equation in (22), (27) and (29) by θ q , we obtain that τ q = 0 , A q = 0 , B q = 0 , C = 0 . Taking into account form (26) of ξ , we get A = τ t / 2 . Hence
τ q = 0 , A q = 0 , A = τ t 2 , B q = 0 , C = 0 .
Differentiate (27) by S and due to (30) obtain the equality γ σ 2 e r ( T t ) D S S = τ t t / 2 , hence D S S S = 0 and D S S q = 0 . Therefore, the differentiation of Equation (28) twice by S gives r D S S + D t S S = 0 and substituting the expression for D S S we get τ t t t = 0 .
Next, differentiating Equation (27) by q and Equation (28) by θ q and S we obtain
τ t 2 + D S q α q + r τ = 0 , D S q = 0 , α q = r τ τ t 2 .
Differentiate (28) by S and q and get
r ( r τ τ t ) + r α q + r A γ σ 2 e r ( T t ) D S S = r τ t τ t t / 2 = 0 .
From this equation and the equality τ t t t = 0 it follows that τ t = 0 .
Therefore, τ is a constant, α q = r τ , α = r q τ + E ( t ) . Substitute it in Equation (27) and obtain γ σ 2 e r ( T t ) ( D S E ) B t = 0 . Hence,
ξ = B ( t ) , α = r q τ + E ( t ) , D = G ( t , q ) + E ( t ) S + e r ( t T ) γ σ 2 B ( t ) S .
Substituting (31) into (28) and reducing we get
r τ F + G t + E S + e r ( t T ) γ σ 2 B S E θ q r G + r B q + μ e r ( t T ) γ σ 2 B B q τ F t G q F θ q = 0 .
Differentiate (32) by θ q and q and get G q q = 0 . Then G = H ( t ) q + J ( t ) and the separation of Equation (32) by q and S gives
G = H ( t ) q + J ( t ) , E + e r ( t T ) γ σ 2 B = 0 , H r H + r B B = 0 , r τ F + J E θ q r J + μ e r ( t T ) γ σ 2 B τ F t H F θ q = 0 .
The third equation in (33) implies that B = H + K e r t . Substitute this equality into the second equation in (33), then
B ( t ) = H ( t ) + K e r t , E ( t ) = t 0 t e r ( s T ) γ σ 2 ( H ( s ) + r 2 K e r s ) d s + L .
Now equalities (31) implies that
ξ = H ( t ) + K e r t , α = r q τ t 0 t e r ( s T ) γ σ 2 ( H ( s ) + r 2 K e r s ) d s + L , η = r θ τ + t 0 t e r ( s T ) γ σ 2 ( H ( s ) + r 2 K e r s ) d s + L S + + e r ( t T ) γ σ 2 ( H ( t ) + r K e r t ) S + H ( t ) q + J ( t ) .
Substituting (34) into the last equation in (33) we get
r τ F τ F t H F θ q + J r J + e r ( t T ) γ σ 2 ( H + r 2 K e r t ) θ q + μ e r ( t T ) γ σ 2 ( H + r K e r t ) = 0 .
This equation has the form r τ F τ F t H ( t ) F θ q + u ( t ) θ q + v ( t ) = 0 . Consider possible situations.

4.1. The Case τ = 0 , H 0

If τ = 0 , H 0 , then K = 0 , J r J = 0 , J = J 0 e r t . Due to (35) we get the generators of symmetry groups X 1 = e r t θ , X 2 = q + S θ for arbitrary F, such that F θ q θ q 0 .

4.2. The Case τ 0 , H 0

If τ 0 , H 0 , then F = Φ 1 ( θ q ) e r t + b ( t ) θ q + c ( t ) . Using the equivalence transformation of the group, which is generated by Y ϕ (18) with ϕ , such that ϕ r ϕ + c = 0 , we obtain F = Φ 1 ( θ q ) e r t + b ( t ) θ q . Since F θ q θ q 0 , then Φ 1 0 . Substitute F in (36), then
r τ b θ q τ b θ q + J r J + e r ( t T ) γ σ 2 r 2 K e r t θ q + μ e r ( t T ) γ σ 2 r K e r t = 0 .
Separating by the variable θ q , obtain
b ( t ) = b 0 e r t + r K e r ( 2 t T ) τ γ σ 2 , J ( t ) = J 0 e r t μ K e r ( 2 t T ) γ σ 2 .
Denote Φ ( θ ) : = Φ 1 ( θ q ) + b 0 θ q , then Φ = Φ 1 0 . Thus,
F = Φ ( θ q ) e r t + 2 r b e 2 r t θ q , Φ 0 , b R , τ = τ 0 , ξ = 2 τ γ σ 2 b e r ( t + T ) , α = r q τ r τ b e 2 r t + L , η = r θ τ + r τ b e 2 r t S + L S + J 0 e r t 2 μ τ b e 2 r t .
Therefore, we obtained the specialization and the symmetry group, which is generated by operators
X 1 = e r t θ , X 2 = q + S θ , X 3 = t + 2 γ σ 2 b 1 e r ( t + T ) S + ( r q r b 1 e 2 r t ) q + ( r θ + r b 1 e 2 r t S 2 μ b 1 e 2 r t ) θ .

4.3. The Case τ = 0 , H 0

If τ = 0 , H 0 , then u 0 , otherwise, F θ q θ q 0 . Therefore, F = a ( t ) θ q 2 + b ( t ) θ q + c ( t ) , a ¬ 0 . We use the equivalence transformation of the group with the generator Y ψ (18), where ψ is a solution of the equation
ψ ( t ) + γ σ 2 e r ( T t ) ( 2 a ( t ) ψ ( t ) + b ( t ) ) = 0 ,
and get F = a ( t ) θ q 2 + c ( t ) , then by a transformation with a generator Y ϕ we obtain the equivalent function F = a ( t ) θ q 2 . Then (36) implies the equation
J + e r ( t T ) γ σ 2 ( H + r 2 K e r t ) θ q r J + μ e r ( t T ) γ σ 2 ( H + r K e r t ) 2 H a θ q = 0 .
Therefore,
H = 2 γ σ 2 e r ( T t ) a ( t ) H r 2 K e r t , J r J + μ e r ( t T ) γ σ 2 ( H + r K e r t ) = 0 .
Solving the second equation in (37) we get
J = J 0 e r t μ e r ( t T ) γ σ 2 ( H + K e r t ) .
Then (35) has the form
τ = 0 , ξ = H ( t ) + K e r t , α = 2 t 0 t a ( s ) H ( s ) d s + L , η = 2 t 0 t a ( s ) H ( s ) d s + L S + e r ( t T ) γ σ 2 ( H ( t ) + r K e r t ) S + + H ( t ) q + J 0 e r t μ e r ( t T ) γ σ 2 ( H ( t ) + K e r t ) .
Let Ψ ( t ) is a partial solution of the first equation in (37) for K = 1 , then a general solution of the equation is H ( t ) = c 1 φ 1 ( t ) + c 2 φ 2 ( t ) + K Ψ ( t ) , where φ 1 and φ 2 are two linearly independent solutions of the homogeneous equation H = 2 γ σ 2 e r ( T t ) a ( t ) H . Therefore, (38) implies that
X 1 = e r t θ , X 2 = q + S θ , X 3 = φ 1 ( t ) S 2 t 0 t a ( s ) φ 1 ( s ) d s q + + e r ( t T ) φ 1 ( t ) γ σ 2 S 2 S t 0 t a ( s ) φ 1 ( s ) d s + φ 1 ( t ) q μ e r ( t T ) γ σ 2 φ 1 ( t ) θ , X 4 = φ 2 ( t ) S 2 t 0 t a ( s ) φ 2 ( s ) d s q + + e r ( t T ) φ 2 ( t ) γ σ 2 S 2 S t 0 t a ( s ) φ 2 ( s ) d s + φ 2 ( t ) q μ e r ( t T ) γ σ 2 φ 2 ( t ) θ , X 5 = ( Ψ ( t ) + e r t ) S 2 t 0 t a ( s ) Ψ ( s ) d s q + + e r ( t T ) γ σ 2 ( Ψ ( t ) + r e r t ) S 2 S t 0 t a ( s ) Ψ ( s ) d s + Ψ ( t ) q μ e r ( t T ) γ σ 2 ( Ψ ( t ) + e r t ) θ .

4.4. The Case τ 0 , H 0

For the case τ 0 , H 0 make a replacement
F = e r t Φ ( t , θ q ) + 1 τ e r t t 0 t e r s u ( s ) d s θ q + 1 τ e r t t 0 t e r s v ( s ) d s
and obtain the equation τ Φ t + H ( t ) Φ θ q = g ( t ) for some g. Therefore, we have F = e r t Φ ( θ q H ( t ) d t / τ ) + b 1 ( t ) θ q + c 1 ( t ) . After using the equivalence transformation of the group for Y ψ (18) with ψ = H ( t ) d t / τ we obtain F = e r t Φ ( θ q ) + b ( t ) θ q + c ( t ) , Φ 0 . Substitute the result in (36), then
r τ b θ q + r τ c τ b θ q τ c e r t H Φ H b + J r J + + e r ( t T ) γ σ 2 ( H + r 2 K e r t ) θ q + μ e r ( t T ) γ σ 2 ( H + r K e r t ) = 0 .
Hence Φ ( θ q ) = a 0 + a 1 θ q + a θ q 2 , by the equivalence transformation for X ψ , where
ψ ( t ) + γ σ 2 e r ( T t ) ( 2 a e r t ψ ( t ) + a 1 + b ( t ) ) = 0 ,
then by an equivalence transformation for a group with a generator X ϕ obtain F = a e r t θ q 2 with a constant a 0 . So, we obtain a partial case to the previous one, but with a nonzero τ , which gives additional symmetry. Thus,
X 1 = e r t θ , X 2 = q + S θ , X 3 = t + r q q + r θ θ , X 4 = φ 1 ( t ) S 2 a t 0 t e r s φ 1 ( s ) d s q + + e r ( t T ) φ 1 ( t ) γ σ 2 S 2 a S t 0 t e r s φ 1 ( s ) d s + φ 1 ( t ) q μ e r ( t T ) γ σ 2 φ 1 ( t ) θ , X 5 = φ 2 ( t ) S 2 a t 0 t e r s φ 2 ( s ) d s q + + e r ( t T ) φ 2 ( t ) γ σ 2 S 2 a S t 0 t e r s φ 2 ( s ) d s + φ 2 ( t ) q μ e r ( t T ) γ σ 2 φ 2 ( t ) θ , X 6 = ( Ψ ( t ) + e r t ) S 2 a t 0 t e r s Ψ ( s ) d s q + + e r ( t T ) γ σ 2 ( Ψ ( t ) + r e r t ) S 2 a S t 0 t e r s Ψ ( s ) d s + Ψ ( t ) q μ e r ( t T ) γ σ 2 ( Ψ ( t ) + e r t ) θ .
Instead of the first equation in (37), we have the equation with constant coefficients
H 2 a γ σ 2 e r T H + r 2 K e r t = 0 .
Therefore, we can calculate a solution of this equation analitically. If a γ > 0 , a r 2 / 2 γ σ 2 e r T , then
φ 1 ( t ) = e 2 a γ σ 2 e r T t , φ 2 ( t ) = e 2 a γ σ 2 e r T t , Ψ ( t ) = r 2 K e r t r 2 2 a γ σ 2 e r T .
For a γ > 0 , a = r 2 / 2 γ σ 2 e r T we have
φ 1 ( t ) = e 2 a γ σ 2 e r T t , φ 2 ( t ) = e 2 a γ σ 2 e r T t , Ψ ( t ) = r K t e r t 2 .
Finally, if a γ < 0 , then
φ 1 ( t ) = sin 2 a γ σ 2 e r T t , φ 2 ( t ) = cos 2 a γ σ 2 e r T t , Ψ ( t ) = r 2 K e r t r 2 2 a γ σ 2 e r T .
The equality a = r 2 / 2 γ σ 2 e r T in this case is not possible.

5. Theorem on Group Classification

Let us formulate the results of calculations in the previous section as the following theorem on group classification. We denote here by X k obtained basis operators in the corresponding Lie algebras.
Theorem 2.
Let r , γ , σ , μ , T R .
1. The Lie algebra for the equation
θ t = r θ + ( μ r S ) q μ θ S σ 2 2 θ S S γ σ 2 2 e r ( T t ) ( θ S q ) 2 + F ( t , θ q ) ,
where F is not equivalent to a ( t ) θ q 2 or e r t Φ ( θ q ) + b 0 e r t θ q + b 1 e 2 r t θ q , F θ q θ q 0 , is generated by the operators
X 1 : = e r t θ , X 2 : = q + S θ .
2. The Lie algebra for the equation
θ t = r θ + ( μ r S ) q μ θ S σ 2 2 θ S S γ σ 2 2 e r ( T t ) ( θ S q ) 2 + e r t Φ ( θ q ) + b e 2 r t θ q ,
where b R , Φ is a nonlinear function, which is not equivalent to a θ q 2 , is generated by the operators
X 1 : = e r t θ , X 2 : = q + S θ , X 3 : = t + 2 γ σ 2 b e r ( t + T ) S + ( r q r b e 2 r t ) q + ( r θ + r b e 2 r t S 2 μ b e 2 r t ) θ .
3. The Lie algebra for the equation
θ t = r θ + ( μ r S ) q μ θ S σ 2 2 θ S S γ σ 2 2 e r ( T t ) ( θ S q ) 2 + a ( t ) θ q 2 ,
where a ( t ) is a nonzero function, which is not equivalent to a 0 e r t , is generated by the operators
X 1 : = e r t θ , X 2 : = q + S θ , X 3 : = φ 1 ( t ) S 2 t 0 t a ( s ) φ 1 ( s ) d s q + + e r ( t T ) φ 1 ( t ) γ σ 2 S 2 S t 0 t a ( s ) φ 1 ( s ) d s + φ 1 ( t ) q μ e r ( t T ) γ σ 2 φ 1 ( t ) θ , X 4 : = φ 2 ( t ) S 2 t 0 t a ( s ) φ 2 ( s ) d s q + + e r ( t T ) φ 2 ( t ) γ σ 2 S 2 S t 0 t a ( s ) φ 2 ( s ) d s + φ 2 ( t ) q μ e r ( t T ) γ σ 2 φ 2 ( t ) θ , X 5 : = ( Ψ ( t ) + e r t ) S 2 t 0 t a ( s ) Ψ ( s ) d s q + + e r ( t T ) γ σ 2 ( Ψ ( t ) + r e r t ) S 2 S t 0 t a ( s ) Ψ ( s ) d s + Ψ ( t ) q μ e r ( t T ) γ σ 2 ( Ψ ( t ) + e r t ) θ .
Here φ 1 , φ 2 are linearly independent solutions of the equation H ( t ) = 2 γ σ 2 e r ( T t ) a ( t ) H ( t ) , Ψ is a partial solution of the equation H ( t ) = 2 γ σ 2 e r ( T t ) a ( t ) H ( t ) r 2 e r t .
4. The Lie algebra for the equation
θ t = r θ + ( μ r S ) q μ θ S σ 2 2 θ S S γ σ 2 2 e r ( T t ) ( θ S q ) 2 + a e r t θ q 2 ,
where a is a nonzero constant, is generated by the operators
X 1 : = e r t θ , X 2 : = q + S θ , X 3 : = t + r q q + r θ θ , X 4 : = φ 1 ( t ) S 2 a t 0 t e r s φ 1 ( s ) d s q + + e r ( t T ) φ 1 ( t ) γ σ 2 S 2 a S t 0 t e r s φ 1 ( s ) d s + φ 1 ( t ) q μ e r ( t T ) γ σ 2 φ 1 ( t ) θ , X 5 : = φ 2 ( t ) S 2 a t 0 t e r s φ 2 ( s ) d s q + + e r ( t T ) φ 2 ( t ) γ σ 2 S 2 a S t 0 t e r s φ 2 ( s ) d s + φ 2 ( t ) q μ e r ( t T ) γ σ 2 φ 2 ( t ) θ , X 6 : = ( Ψ ( t ) + e r t ) S 2 a t 0 t e r s Ψ ( s ) d s q + + e r ( t T ) γ σ 2 ( Ψ ( t ) + r e r t ) S 2 a S t 0 t e r s Ψ ( s ) d s + Ψ ( t ) q μ e r ( t T ) γ σ 2 ( Ψ ( t ) + e r t ) θ ,
where φ 1 , φ 2 , Ψ are from (39), 40), or (41), depending on the sign of a γ and the value of a.
Remark 3.
In the second part of this theorem at b = 0 and in the fourth one we have the market trading volume V t = a e r t with a constant a 0 , as multiplier at a function of θ q in an expression for F. If Φ 0 in the second part, then the market trading volume is V t = b e 2 r t . In the third part of the theorem V t = a ( t ) .
Remark 4.
A theorem on the group classification of Equation (2) with a free element F depending on θ q only is obtained in [35]. It contains the specifications F = e ν θ q and F = θ q 2 , which correspond to additional symmetries of the equation.

6. Application to the Search of Some Submodels

Using a symmetry group for a differential equation we can reduce the number of variables on which an unknown function depends by the dimension of the considered group. If the resulting equation can be integrated, we obtain an exact solution of the original equation, invariant with respect to the group of symmetries under consideration. If the resulting equation is not integrable, following L.V. Ovsyannikov [22], we will call such an equation an invariant submodel of the initial equation (initial model).
In order to find invariant solutions or submodels that are not translated into each other by transformations of variables, we must find the so-called optimal system of subalgebras of the Lie algebra of the equation under study. To do this, the internal automorphisms of this algebra are used, which can be found through nonzero structural constants of the algebra. Below we will do this for the two simplest Lie algebras of the symmetry groups generators obtained in Theorem 2 on the group classification. Invariant solutions or invariant submodels for different subalgebras of the optimal system are not equivalent, i.e., they cannot be obtained from each other by replacing variables.

6.1. Optimal System of Subalgebras and Submodels for the General Case

Lie algebra L 2 (43) is commutative, hence it has no continuous groups of internal automorphisms. Thus, its optimal system of one-dimensional subalgebras is Θ 1 = { X 2 , X 1 + c X 2 , c R } .
The subalgebra X 2 has the invariants J 1 = t , J 2 = S , J 3 = θ q S . Writing J 3 = w ( J 1 , J 2 ) we obtain the form of the corresponding invariant solution θ = w ( t , S ) + S q . Substitute it into Equation (42) and obtain the submodel
w t = r w μ w S σ 2 2 w S S γ σ 2 2 e r ( T t ) w S 2 + F ( t , S ) ,
which is invariant for X 2 . Analogously we get the invariants t, e r t c + S , θ e r t c q S q for the subalgebra X 1 + c X 2 , c 0 . The invariant submodel for it has the form
w t = r w μ w S σ 2 2 w S S γ σ 2 2 e r ( T t ) w S 2 + F t , e r t c + S ,
where θ = w ( t , S ) + e r t c q + S q . If c = 0 , then the subalgebra X 1 has no invariant submodels, since its invariants t , S , q do not depend on θ .

6.2. Optimal System of Subalgebras and Submodels for the Specification F = Φ ( θ q ) e r t

Nonzero structural constants for a Lie algebra with a basis { X 1 , X 2 , , X n } are coefficients c i j k in the decomposition of a commutator [ X i , X j ] [17] by the basis: [ X i , X j ] = k = 1 n c i j k X k . Generators of continuous groups of internal automorphisms can be calculated by the formula E i = j , k = 1 n c i j k e j e k , where e j are coefficients in the decomposition of an element of the Lie algebra by its basis, which depend on group parameters.
Consider the Lie algebra L 3 with basis (45). For L 3 we have c 23 1 = c 32 1 = 2 γ σ 2 b e r T , c 23 2 = c 32 2 = r . The integration of the Lie equations for the generators gives E 2 : e ¯ 1 = e 1 2 γ σ 2 b e r T e 3 a 2 , e ¯ 2 = e 2 + r e 3 a 2 ; E 3 : e ¯ 1 = e 1 2 γ σ 2 r b e r T e 2 ( 1 e r a 3 ) , e ¯ 2 = e r a 3 e 2 . Also, we add a mirror automorphism E : e ¯ 1 = e ¯ 1 , which does not change the commutators of the basis operators of this Lie algebra L 3 .
Let b 0 , for e 3 0 by the internal automorphism E 2 obtain e 2 = 0 , then we have ( e 1 , e 2 , e 3 ) = ( c , 0 , 1 ) after scaling, i.e., we get c X 1 + X 3 , c R . If e 3 = 0 , e 2 0 , then by E 3 get e 1 = 0 , ( e 1 , e 2 , e 3 ) = ( 0 , 1 , 0 ) , if we take
a 3 = 1 r ln 1 r e r T e 1 2 γ σ 2 b e 2
in the case r e r T e 1 2 γ σ 2 b e 2 < 1 . If r e r T e 1 2 γ σ 2 b e 2 1 , we will use E to go to the previous case. For e 2 = e 3 = 0 we have ( e 1 , e 2 , e 3 ) = ( 1 , 0 , 0 ) . Thus, Θ 1 1 = { X 1 , X 2 , c X 1 + X 3 , c R } .
Let us search for a system of two-dimensional subalgebras for L 3 with b 0 . For the basis vector X 1 of the one-dimensional subalgebra X 1 , consider the second basis vector in the form α X 2 + β X 3 , then the commutator has the form [ X 1 , α X 2 + β X 3 ] = 0 . We get subalgebras X 1 , X 2 for e 3 = 0 , X 1 , X 3 for e 3 0 , if we use E 2 .
For the basis vector X 2 , consider the second basis vector in the form α X 1 + β X 3 . Their commutator is [ X 2 , α X 1 + β X 3 ] = r β X 2 2 β γ σ 2 b e r T X 1 . Therefore, a subalgebra is formed at β = 0 , which is already found in the form X 1 , X 2 .
For c X 1 + X 3 , consider the second basis vector in the form α X 1 + β X 2 . Then we have [ c X 1 + X 3 , α X 1 + β X 2 ] = 2 β γ σ 2 b e r T X 1 r β X 2 and get the subalgebra c X 1 + X 3 , 2 γ σ 2 b e r T X 1 r X 2 . By E 3 and E reduce it to c X 1 + X 3 , X 2 .
Thus, we proved the optimal systems of one-dimensional and two-dimensional subalgebras of the Lie algebra under consideration.
Lemma 1.
Optimal systems of one-dimensional and two-dimensional subalgebras of Lie algebra L 3 (45) with b 0 are
Θ 1 1 = { X 1 , X 2 , c X 1 + X 3 , c R } , Θ 2 1 = { X 1 , X 2 , X 1 , X 3 , c X 1 + X 3 , X 2 , c R } .
Here a denotation Θ k j is used for the optimal system of k-dimensional subalgebras in the j-th case.
In the case b = 0 , for e 3 0 we obtain the vector ( c , 0 , 1 ) , c R , using E 2 . If e 3 = 0 , then using E 3 , E we get ( 1 , 1 , 0 ) , ( 1 , 0 , 0 ) , ( 0 , 1 , 0 ) . So, Θ 1 2 = { X 1 , X 2 , X 1 + X 2 , c X 1 + X 3 , c R } . In this case, we have two-dimensional subalgebras X 1 , X 2 , X 1 , X 3 also. Moreover, [ X 2 , α X 1 + β X 3 ] = r β X 2 , and we have the subalgebra X 2 , α X 1 + β X 3 for any α , β R . If β = 0 , it will be a partial case of X 1 , X 2 , for β 0 we obtain the subalgebra c X 1 + X 3 , X 2 . Since [ c X 1 + X 3 , α X 1 + β X 2 ] = r β X 2 , we obtain another subalgebra c X 1 + X 3 , X 1 .
Lemma 2.
The optimal system of one-dimensional and two-dimensional subalgebras of Lie algebra L 3 (45) with b 1 = 0 are Θ 1 2 = { X 1 , X 2 , X 1 + X 2 , c X 1 + X 3 , c R } and Θ 2 2 = { X 1 , X 2 , c X 1 + X 3 , X 1 , c X 1 + X 3 , X 2 , c R } .
The subalgebras X 1 , X 1 , X 2 , X 1 , X 3 do not have invariant submodels, since X 1 does not have invariants depending on θ .
Consider the case b = 0 , then the subalgebra X 1 + X 2 has invariants t, S, θ ( e r t + S ) q , therefore, an invariant solution has the form θ = w ( t , S ) + ( e r t + S ) q and the invariant submodel is
w t = r w μ w S σ 2 2 w S S γ σ 2 2 e r ( T t ) w S 2 + F ( t , e r t + S ) .
The subalgebra c X 1 + X 3 has invariants x : = q e r t , S, θ e r t c t , hence we will look for an invariant solution in the form θ = c t e r t + e r t w ( q e r t , S ) , where w is a function of two variables. Substitute it into (44) and obtain the invariant for c X 1 + X 3 submodel
Φ ( w x ) + r x w x = σ 2 2 w S S + μ w S + γ σ 2 2 e r T ( w S x ) 2 + ( r S μ ) x + c .
The subalgebra c X 1 + X 3 , X 1 has no invariants depending on θ and, therefore, invariant submodels. Let us find the invariant submodel with respect to c X 1 + X 3 , X 2 . Consider a function G = G ( x , S , y ) , where x : = q e r t , y : = θ e r t c t are invariants for the subalgebra c X 1 + X 3 . Then X 2 G = e r t G x + S e r t G y and invariants of the subalgebra c X 1 + X 3 , X 2 are S and y S x = ( θ S q ) e r t c t . Therefore, we will search an invariant solution for this subalgebra in the form θ = c t e r t + S q + e r t w ( S ) . The invariant submodel will have the form
w ( S ) + 2 μ σ 2 w ( S ) + γ e r T w ( S ) 2 2 σ 2 Φ ( S ) + 2 c σ 2 = 0 .

7. Conclusions

In this paper, we develop a theorem on group classification of the Guéant and Pu model of the option pricing taking into account transaction costs and the impact of operations on the market. For this aim, the Lie algebra of generators of continuous groups of equivalence transformations is calculated. For the general case, and for the case of the equation with the right-hand side F = e r t Φ ( θ q ) , optimal systems of subalgebras and corresponding invariant submodels are derived. The results of this work will be applied to the analogous research into the Guéant and Pu model with the specifications F = e r t Φ ( θ q ) + b e 2 r t , F = a ( t ) θ q 2 , F = a e r t θ q 2 , which is presented in the obtained theorem on group classification. The knowledge of the group structure of the studied models obtained in this way will make it possible to calculate their exact solutions and conservation laws.

Author Contributions

Conceptualization, S.M.S. and V.E.F.; methodology, V.E.F. and K.V.Y.; software, K.V.Y.; validation, V.E.F. and K.V.Y.; formal analysis, K.V.Y.; investigation, K.V.Y.; resources, S.M.S.; data curation, S.M.S.; writing—original draft preparation, V.E.F.; writing—review and editing, S.M.S. and V.E.F.; visualization, K.V.Y.; supervision, S.M.S. and V.E.F.; project administration, S.M.S.; funding acquisition, S.M.S. and V.E.F. All authors have read and agreed to the published version of the manuscript.

Funding

The work is supported by the Ministry of Science and Higher Education of the Russian Federation, agreement No. 075-02-2022-881, 2 February 2022. The work of V.E.F. is funded by the grant of the President of the Russian Federation to support leading scientific schools, project number NSh-2708.2022.1.1.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Black, F. The pricing of Commodity Contracts. J. Financ. Econ. 1976, 3, 167–179. [Google Scholar] [CrossRef]
  2. Black, F.; Scholes, M. The pricing of options and corporate liabilities. J. Political Econ. 1973, 81, 637–659. [Google Scholar] [CrossRef]
  3. Magill, M.J.P.; Constantinides, G.M. Portfolio selection with transactions costs. J. Econ. Theory 1976, 13, 245–263. [Google Scholar] [CrossRef]
  4. Kyle, A.S. Continuous auctions and insider trading. Econometrica 1985, 53, 1315–1335. [Google Scholar] [CrossRef]
  5. Leland, H.E. Option pricing and replication with transactions costs. J. Financ. 1985, 40, 1283–1301. [Google Scholar] [CrossRef]
  6. Cvitanić, J.; Karatzas, I. Hedging and portfolio optimization under transaction costs: A martingale approach. Math. Financ. 1996, 6, 133–165. [Google Scholar] [CrossRef]
  7. Barles, G.; Soner, H.M. Option pricing with transaction costs and a nonlinear Black—Scholes equation. Financ. Stochastics 1998, 2, 369–397. [Google Scholar] [CrossRef]
  8. Grossman, S. An analysis of the implications for stock and futures price volatility of program trading and dynamic hedging strategies. J. Bus. 1998, 61, 275–298. [Google Scholar] [CrossRef]
  9. Platen, E.; Schweizer, M. On feedback effects from hedging derivatives. Math. Financ. 1998, 8, 67–84. [Google Scholar] [CrossRef]
  10. Sircar, R.; Papanicolaou, G. Generalized Black—Scholes models accounting for increased market volatility from hedging strategies. Appl. Math. Financ. 1998, 5, 45–82. [Google Scholar] [CrossRef]
  11. Schönbucher, P.; Wilmott, P. The feedback-effect of hedging in illiquid markets. SIAM J. Appl. Math. 2000, 61, 232–272. [Google Scholar]
  12. Bank, P.; Baum, D. Hedging and portfolio optimization in financial markets with a large trader. Math. Financ. 2004, 14, 1–18. [Google Scholar] [CrossRef]
  13. Çetin, U.; Jarrow, R.; Protter, P. Liquidity risk and arbitrage pricing theory. Financ. Stoch. 2004, 8, 311–341. [Google Scholar] [CrossRef]
  14. Çetin, U.; Rogers, L.C. Modelling liquidity effects in discrete time. Math. Financ. 2007, 17, 15–29. [Google Scholar] [CrossRef]
  15. Rogers, L.C.; Singh, L.S. The cost of illiquidity and its effects on hedging. Math. Financ. 2010, 20, 597–615. [Google Scholar] [CrossRef]
  16. Gazizov, R.K.; Ibragimov, N.H. Lie symmetry analysis of differential equations in finance. Nonlinear Dyn. 1998, 17, 387–407. [Google Scholar] [CrossRef]
  17. Ovsiannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  18. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: New York, NY, USA; Berlin/Heidelberg, Germany; Tokyo, Japan, 1986. [Google Scholar]
  19. Turab, A.; Sintunavarat, W. A unique solution of the iterative boundary value problem for a second-order differential equation approached by fixed point results. Alex. Eng. J. 2021, 60, 5797–5802. [Google Scholar] [CrossRef]
  20. Li, W.-J.; Yu, Y.-H. Oscillatory behavior of third-order nonlinear differential equations with a sublinear neutral term. Acta Math. Appl. Sin. Engl. Ser. 2022, 38, 484–496. [Google Scholar] [CrossRef]
  21. Santra, S.S.; Scapellato, A. Some conditions for the oscillation of second-order differential equations with several mixed delays. J. Fixed Point Theory Appl. 2022, 24, 18. [Google Scholar] [CrossRef]
  22. Ovsyannikov, L.V. The “podmodeli” programm. Gas dynamics. J. Appl. Math. Mech. 1994, 58, 601–627. [Google Scholar] [CrossRef]
  23. Annin, B.D.; Bytev, V.O.; Senashev, C.I. Group Properties of the Equations of Elasticity and Plasticity; Nauka: Novosibirsk, Russia, 1985. (In Russian) [Google Scholar]
  24. Bordag, L.A. Geometrical Properties of Differential Equations: Applications of the Lie Group Analysis in Financial Mathematics; World Scientific Publishing Company: Singapore, 2015. [Google Scholar]
  25. Bordag, L.A.; Mikaelyan, A. Models of self-financing hedging strategies in illiquid markets: Symmetry reductions and exact solution. J. Lett. Math. Phys. 2011, 96, 191–207. [Google Scholar] [CrossRef] [Green Version]
  26. Dyshaev, M.M. Group analysis of a nonlinear generalization for Black—Scholes equation. Chelyabinsk Phys. Math. J. 2016, 1, 7–14. [Google Scholar]
  27. Dyshaev, M.M.; Fedorov, V.E. Symmetry analysis and exact solutions for a nonlinear model of the financial markets theory. Math. Notes NEFU 2016, 23, 28–45. [Google Scholar]
  28. Fedorov, V.E.; Dyshaev, M.M. Group classification for a general nonlinear model of option pricing. Ural Math. J. 2016, 2, 37–44. [Google Scholar] [CrossRef]
  29. Dyshaev, M.M.; Fedorov, V.E. Symmetries and exact solutions of a nonlinear pricing options equation. Ufa Math. J. 2017, 9, 29–40. [Google Scholar] [CrossRef]
  30. Fedorov, V.E.; Dyshaev, M.M. Invariant solutions for nonlinear models in illiquid markets. Math. Methods Appl. Sci. 2018, 41, 8963–8972. [Google Scholar] [CrossRef]
  31. Dyshaev, M.M.; Fedorov, V.E. Comparing of some sensitivities (Greeks) for nonlinear models of option pricing with market illiquidity. Math. Notes NEFU 2019, 26, 94–108. [Google Scholar]
  32. Fedorov, V.E.; Dyshaev, M.M. Group classification for a class of non-linear models of the RAPM type. Commun. Nonlinear Sci. Numer. Simul. 2021, 92, 105471. [Google Scholar] [CrossRef]
  33. Guéant, O.; Pu, J. Option pricing and hedging with execution costs and market impact. arXiv 2015, arXiv:1311.4342. [Google Scholar]
  34. Guéant, O. The Financial Mathematics of Market Liquidity: From Optimal Execution to Market Making; Chapman and Hall/CRC: Boca Raton, FL, USA, 2016. [Google Scholar]
  35. Yadrikhinskiy, K.V.; Fedorov, V.E.; Dyshaev, M.M. Group analysis of the Guéant and Pu model of option pricing and hedging. In Nonlinear Physical Science; Springer Nature: Singapore, 2021; pp. 173–204. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sitnik, S.M.; Yadrikhinskiy, K.V.; Fedorov, V.E. Symmetry Analysis of a Model of Option Pricing and Hedging. Symmetry 2022, 14, 1841. https://doi.org/10.3390/sym14091841

AMA Style

Sitnik SM, Yadrikhinskiy KV, Fedorov VE. Symmetry Analysis of a Model of Option Pricing and Hedging. Symmetry. 2022; 14(9):1841. https://doi.org/10.3390/sym14091841

Chicago/Turabian Style

Sitnik, Sergey M., Khristofor V. Yadrikhinskiy, and Vladimir E. Fedorov. 2022. "Symmetry Analysis of a Model of Option Pricing and Hedging" Symmetry 14, no. 9: 1841. https://doi.org/10.3390/sym14091841

APA Style

Sitnik, S. M., Yadrikhinskiy, K. V., & Fedorov, V. E. (2022). Symmetry Analysis of a Model of Option Pricing and Hedging. Symmetry, 14(9), 1841. https://doi.org/10.3390/sym14091841

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop