# On the Resolution of a Remarkable Bond Pricing Model from Financial Mathematics: Application of the Deductive Group Theoretical Technique

Authors: Taha Aziz
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The classical Cox–Ingersoll–Ross (CIR) bond-pricing model is based on the evolution space-time dependent partial differential equation (PDE) which represents the standard European interest rate derivatives. In general, such class of evolution partial differential equations (PDEs) has generally been resolved by classical methods of PDEs and by ansatz-based techniques which have been previously applied in a similar context. The author here shows the application of an invariant approach, a systematic method based on deductive group-theoretical analysis. The invariant technique reduces the scalar linear space-time dependent parabolic PDE to one of the four classical Lie canonical forms. This method leads us to exactly solve the scalar linear space-time dependent parabolic PDE representing the CIR model. It was found that CIR PDE is transformed into the first canonical form, which is the heat equation. Under the proper choice of emerging parameters of the model, the CIR equation is also reduced to the second Lie canonical form. The equivalence transformations which map the CIR PDE into the different canonical forms are deduced. With the use of these equivalence transformations, the invariant solutions of the underlying model are found by using some well-known results of the heat equation and the second Lie canonical form. Furthermore, the Cauchy initial-value model of the CIR problem along with the terminal condition is discussed and closed-form solutions are deduced. Finally, the conservation laws associated with the CIR equation are derived by using the general conservation theorem.


---

## Body

## 1. Theory and Background

The bond-pricing theory began in the 19th century when Bachelier [1] gave the assumption that the pricing of bond depends on stock price which follows a Brownian motion. The pioneering work of Merton [2] and later by Black and Scholes [3] also gave a new dimension to mathematical formulation of finance problems. The actual mathematical formulation of these models was initially given in the form of stochastic differential equations (DEs), but after incorporating certain assumptions, these models were derived in terms of parabolic PDEs with variable coefficients. The fundamental model of bond-pricing theory was first proposed by Vasicek [4], which was modified by Cox et al. in [5], now famous in literature as the CIR bond-pricing model.For the Vasicek model, the instantaneous ratezt is represented by an Ornstein–Uhlenbeck process, which is written in the form [4, 5]:(1)dzt=βdWt+μα−ztdt,with Wt is the standard Brownian motion and μα−zt is the drift term. For the CIR model, the instantaneous rate followed a stochastic process of the form(2)dzt=βztdWt+μα−ztdt.The mean-reversion feature of these bond-pricing models is consistent with economic phenomenon. However, the main disadvantage of the Vasicek model is that it allows the negative interest rates. On the contrary, the CIR model provides improvement in the Vasicek model by ensuring that the interest rate can never be negative. The CIR model, therefore, incorporates the number of important characteristics of interest rate models.For Vasicek model, the price bondFz,t satisfies the following parabolic equation [6, 7]:(3)Ft+12β2Fzz+μα−zFz−zF=0.In the CIR model, the price of a zero-couponFz,t satisfies the parabolic 1+1 PDE [6, 7], which is written as(4)Ft+12β2zFzz+μα−zFz−zF=0.For both the above PDEs, the valueFz,t of the zero-coupon bond matures at time t=T with a value of unity, i.e.,(5)Fz,T=1,which is known as terminal condition, where T denotes the expiration time and gives the payoff of the derivative on maturity.A mapping has been constructed in [8] which embeds the Vasicek model into the CIR model such that, under this mapping, any solution of the CIR equation can be transformed into a solution of the Vasicek equation. The solutions of the CIR equation are, therefore, rather more general and provide deeper insights into interest rate models.It is important to note that the general PDE in which both Vasicek equation (3) and CIR equation (4) are embedded is given by(6)Ft+12β2Fzz+F−αβFz−zF=0.The above equation models the price of a zero-coupon bond when the spot rate is determined by the general stochastic process:(7)dzt=βx,tdWt+Fz,tdt.The functionαz,t in equation (7) represents the market price of risk. The particular cases of PDE (10) have been considered in [9], in which the authors performed the group classification of the general equation by using some simplifying assumptions. The group invariant solutions of the CIR equation (4) were first found in [10, 11]. The complete Lie symmetry group admitted by (4) was found for all choices of the parameters. It was found that, besides trivial symmetries, some additional symmetries were obtained. The self-similar solutions of PDE (4) subject to the terminal condition (5) were deduced for special parametric choices (refer to [10, 11]).The applications of Lie groups to problems in mathematical finance are not new, and it is certainly young enough, with a possibility to expand. Generally, the nonlinear DEs that arise in mathematical finance are usually tackled with either approximate or numerical methods. Whereas, the Lie symmetry approach allows us to find complete local one-parameter transformation groups that can be used to find exact solutions or to reduce nonlinear DEs to linear DEs. Concerning the relationship between Lie symmetry methods and financial mathematics and economic models, we would like to note the presence of several recent studies [12–22], which represent the applications and utility of well-known concepts of Lie theory to several problems from mathematical finance.Lie [23] was the first to introduce the classification of DEs in terms of continuous groups. Lie classified the scalar second-order (1 + 1) PDEs and developed a method for their integration. Lie provided all the canonical forms of the scalar linear (1 + 1) PDE for which the equation admits six-, four-, two-, and one-dimensional nontrivial symmetry algebras. Ovsiannikov [24] further investigated the parabolic linear (1 + 1) PDE by using the reduction to the fourth Lie canonical form. Bluman in [25] studied the symmetry properties of the parabolic linear (1 + 1) PDE and developed a mapping algorithm. Johnpillai and Mahomed [26] were the first to propose the invariant approach to Fokker–Planck equations. They proved that the linear (1 + 1) parabolic equation was reducible to the heat equation via a local equivalence transformation. The necessary and sufficient conditions were discussed in [26], for a scalar linear (1 + 1) parabolic equation which can be reduced to a one-dimensional heat equation through equivalence transformations. A few years later, Mahomed [27] discussed the complete invariant characterization of scalar linear (1 + 1) parabolic PDE into four canonical forms. The refined invariant criteria, for the reduction of (1 + 1) parabolic equation into different canonical forms, in terms of the coefficients of the parabolic equation were derived in [27]. Recently, the usefulness of the invariant criterion has been thoroughly demonstrated in studies [28–30].From the definition point of view, a general scalar1+1 PDE is written as(8)Ft=At,zFzz+Bt,zFz+Ct,zF,with At,z, Bt,z, and Ct,z are continuously differentiable functions. The equivalence transformations of this PDE are an infinite local group of linear transformations of the dependent variables, given as(9)F¯=Dt,zF.Together with equivalence transformations of independent variables,(10)t¯=ft,z¯=Wt,z,ft≠0,Wx≠0,where ft, Dt,z, and Wt,z are arbitrary functions. Two parabolic equations of the form given in equation (8) are considered to be equivalent only if they can be mapped to one another by the suitable combinations of equivalence transformations. Lie was the first to deduce that every linear scalar (1 + 1) parabolic PDE can be mapped on the following four canonical forms:(11)Ft=Fzz,Ft=Fzz+Az2F,a≠0,Ft=Fzz+CzF,c≠0,Ft=Fzz+Ct,zF,c≠0.The first form is the classical heat PDE which contains six nontrivial symmetry generators. The second canonical equation has four nontrivial symmetry generators, the third generally has two symmetry generators, and the fourth has one nontrivial symmetry generator. The following theorems give the algorithm on the invariant characterization of scalar1+1 parabolic PDE into different canonical forms given in equation (11).Theorem 1.
The scalar1+1 parabolic PDE (8) can be reduced to the heat equation by a transformation if and only if the coefficients in equation (8) satisfy the condition(12)2lz+2mz−nz=0,where(13)l=A1/2A1/2Πzz,m=A1/2A1/2∂tB/2Az,n=A1/2∂t21/A1/2,with Π given by(14)Π=C−Bz2+BAz2A+Azz4−316Az2a−At2A−B24A.Theorem 2.
The scalar linear1+1 parabolic PDE (8) is reducible to the second Lie canonical equation if and only if the coefficients in equation (8) satisfy the following invariant equation:(15)20lz+20mz−10nz+10A1/2mzz∫dzA1/2−5A1/2nzz∫dzA1/2+10A1/2lzz∫dzA1/2+A1/2A1/2lzzz∫dzA1/22+A1/2A1/2mzzz∫dzA1/22−12A1/2A1/2nzzz∫dzA1/22=0,with l, m, n, and Π are defined in equations (13) and (14).Theorem 3.
The parabolic PDE (8) reduces to the classical heat PDE by the transformation relations:(16)t¯=χt,z¯=±∫χ˙At,z−11/2dz+σt,F¯=ftAt,z−1/4Fexp∫Bt,z2At,zdz−18χ¨χ˙∫dzAt,z1/22−12∫1At,z1/2∂∂t∫dzAt,z1/2dz±12σ˙σ˙1/2∫dzAt,z1/2,where χ, σ, and f satisfy the relation(17)χ˙c¯=Π+∂∂t∫Bt,z2At,zdz−12∫1At,z1/2∂2∂t2∫dzAt,z1/2dz+ut∫dzAt,z1/22+vt∫dzAt,z1/2+wt,where Π is defined in equation (14) and(18)ut=116χ¨2χ˙2−18χ¨χ˙t,vt=±14χ¨χ˙σ˙χ˙1/2±12σ˙χ˙1/2t,wt=14χ¨χ˙+14σ˙2χ˙+f˙f.
The functionsu, v, and w are constrained by relation (18).
The derivation of conservation laws plays an important role in understanding the mathematical models described by DEs. For several mathematical and physical models, the conservation laws are used to study the existence and uniqueness of solutions of nonlinear DEs. The conservation laws are applied in the case of PDEs to search for possible potential symmetries and related nonclassical symmetries. Consequently, it is necessary to derive the conservation laws of the CIR equation, arguably one of the most popular mathematical finance models. There are many approaches available to the construction of conservation laws for PDEs [31–37]. The classical Noether theorem provides a sophisticated way to construct conservation laws once the Lagrangian associated with Noether symmetries are found for the corresponding Euler–Lagrange equation. Recently, a general result has been established on the conservation laws for arbitrary DEs, which is based on the concept of an adjoint equation for associated DE [38]. According to this result, it is possible to find a conservation law with any group of Lie symmetries and to derive conservation laws without seeking classical Lagrangians.
Motivated by the facts stated above, the purpose of this paper is to describe the systematic application of invariant criterion to precisely solve the traditional CIR model from mathematical finance. The majority of financial interest rate problems are modeled using linear parabolic (1 + 1) PDEs. Because many of the DEs in the literature addressing financial models are highly rich in symmetry, the Lie group theory of transformation groups is widely useful. As a result, the remainder of the study demonstrates the algebraic features of the CIR model in the context of an invariant criterion based on Lie groups. The closed-form solutions will be derived for the model problem under a suitable choice of terminals condition. By not using the ansatz method or relying on adhoc approach, this will be achieved through an algorithmic technique: the invariant criterion.


## 2. Fundamental Solutions of Cox–Ingersoll–Ross Bond-Pricing Equation

The fundamental solutions of CIR bond-pricing equation under the implication of invariant criterion are constructed in this section. Two separate cases will be considered for the reduction of CIR equation (4) into heat equation and second canonical equation. Thereafter, the invariant solutions of the model are deduced.


### 2.1. Case 1: Reduction to Heat Equation—Ft=Fzz

First, find the equivalence transformations which reduce equation (4) into different canonical forms. With the comparison of PDE (4) with scalar linear parabolic PDE (8), the coefficients functions are given by(19)At,z=−12β2z,Bt,z=−μα−z,Ct,z=z.With the use of Theorem1, the parametric values are deduced such that the PDE (4) is mapped to the heat equation. The value of Π given in equation (14) for the CIR equation is found as(20)Π=μ2z−α22β2z+μz−α2z+332β2z−μ2+z.It is clear thatm and n are equal to zero in the condition defined in Theorem 1 because the coefficients of PDE (4) are independent of time t. Under the value of Π, condition (12) takes the form(21)lz=z3−38β4+2μαβ2−2μ2α2.From the above equation, one obtains the following parametric conditions:(22)β1=−4μα3,β2=4μα3,β3=−4μα,β4=4μα.For specific values of the parameterβ, equation (22) gives the nontrivial cases such that invariant condition (12) is holds, which reduces the CIR PDE (4) to the heat equation.Now, the transformation relations will be computed which give reduction of the CIR equation (4) to the heat equation. These transformations are further used to find the fundamental solutions of the PDE (4).With the use of the coefficients as defined in Theorem3, one obtains the following results:(23)ut=−18β2+μ22,vt=0,wt=μ2αβ2,χt=k12KtanK2t,σt=k2∫dtlncosK/2t,ft=exp12lncosK/2t−k324∫dtlncosK/2t2+t+k4,where K=8ft and kjj=1,2,3,4 are the integration constants and β can take any value as given in equation (22). Substituting the values from equation (23) to equation (16), the equivalence transformations which give the reduction of CIR equation (4) into the heat equation (after performing some lengthy and tedious calculation) are given as(24)t¯=k12/Ktan2/Kt,z¯=±4k1cos2/Ktz+k2∫tTdtlncos2/Kt+k3,F¯=12β2z−1/4exp12lncos2/Kt−k424∫tTdtlncos2/Kt2+k5×Fz,texpμβ2αlnz−z±2k3lncos2/Ktz,with kjj=1,…,5 are integration constants. The terminal condition (5) for the CIR equation is also transformed to(25)Fz,T=12β2z−1/4exp12lncosK/2t+μ2αβ2T+k5×expμβ2αlnz−z±2k3lncosK/2tz.In order to seek for the fundamental solutions of the PDE (4), it is seen that there exist equivalence transformations (24) which gives reduction of PDE (4) to the heat equation. This result is employed to find the fundamental solution of the CIR equation.It is well known that the fundamental solution of the heat is given in barred coordinates [39]:(26)F¯z¯,t¯=12t¯πexp−z¯24t¯.To look for the fundamental solution of CIR equation (4), solution (26) is transformed by means of (24). With the use of equations (24) and (26), the solution Fz,t is given by(27)Fz,t=12β2z1/412t¯πexp−z24t¯×exp−12lncosK/2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz.Finally, using the values oft¯ and z¯ into equation (27), the above solution takes the form(28)Fz,t=12β2z1/412k1K/2tanK/2tπ×exp−±4k1/cosK/2tz+k2∫tTdt/lncosK/2t+k324k1K/2tanK/2t×exp−12lncosK2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz.Under the same set of equivalence transformations for reduction into the heat equation, one can find another closed-form solution for PDE (4). The fundamental solution of the heat equation is also given as series solution of the form [37](29)F¯z¯,t¯=∑r=1∞prexp−rπL2qt¯sinrπz¯L,which converges for all 0≤x≤L. Again, one needs to transform the solution given in equation (29) subject to the transformations given in equation (24). This solution is written as(30)Fz,t=12β2z1/4∑r=1∞prexp−rπL2qt¯sinrπz¯L×exp−12lncosK2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz.Substituting the values ofz¯ and t¯, the series solution Fz,t takes the form(31)Fz,t=12β2z1/4∑r=1∞prexp−rπL2qk12KtanK2t×∑r=1∞sinrπ±4k1/cosK/2tz+k2∫tTdt/lncosK/2t+k3L×exp−12lncosK2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz.


### 2.2. Case 2: Reduction to Second Lie Canonical Equation—Ft=Fzz+A/z2Fz,t

We look for the case when the CIR PDE (4) is transformed to the second Lie canonical equation. In this situation, one deduce following parametric constraint(32)β2≠−8μα±4μα−3⇒lz≠0.Under the same procedure used previously, the equivalence transformations that transform equation (4) into the second Lie canonical form are given by(33)t¯=k12KtanK2t,z¯=±8k1cosK/2tz,F¯=cosK/2texpαt+k2expαlnz−z−2Kβ2tanK/2tzF,with k1 and k2 are the integration constants.Now, to deduce the closed-form solution such that the CIR equation is transformed into the second Lie canonical equation subject to the equivalence transformations given in equation (33), first, rewrite the second canonical equation:(34)F¯t¯−F¯z¯z¯−Pz¯2F¯=0.In order to find the exact solution of CIR equation, one first needs to seek the solution of PDE (34) and then employ transformations (33). To obtain the solution for PDE (34), the separation of variables technique is used.Let us assume that the PDE (34) has the separable solution of the form(35)F¯z¯,t¯=Uz¯Vt¯.Using equation (35) in equation (34) and separating the variables, one obtains(36)Vt¯+γV=0,(37)Uz¯z¯+Pz¯2+γU=0,with γ≥0 is the separation constant. The following two independent cases are considered here.Case I.
γ=0
The solutions of ODEs (36) and (37) are written as(38)Vt¯=c1,(39)Uz¯=c2z¯1/21−1−4P+c3x¯−1/21−1−4P,Subject to equations (38) and (39), and the exact solution of equation (34) is written as(40)F¯z¯,t¯=c1c2z¯1/21−1−4P+c3x¯−1/21−1−4P.
In order to obtain the solution of PDE (4), solution (40) is transformed subject to the transformation relations given in equation (33). One can write the solution Fz,t as(41)Fz,t=c1c2z¯1/21−1−4P+c3x¯−1/21−1−4P,×1cosK/2texp−αt+k2.
Substituting the values ofz¯ and t¯ from equation (33) into equation (41), the solution for Fz,t takes the form(42)Fz,t=c1c2±8k1cosK/2tz1/21−1−4P+c3±8k1cosK/2tz−1/21−1−4P×exp−αlnz−z−2Kβ2tanK2tz×1cosK/2texp−αt+k2.Case II.
γ>0
Here, the solutions of the ODEs (36) and (37) are found as(43)Vt¯=c1exp−γt¯,Uz¯=c2z¯J1121−4P;γz¯+c3z¯J2121−4P;γz¯,with J1 and J2 as the first and the second kind Bessel functions [25], respectively. Therefore, the solution of equation (34) is given by(44)F¯z¯,t¯=c1exp−γt¯c2z¯J1121−4P;γz¯+c3z¯J2121−4P;γz¯.
The solution of CIR PDE (4), subject to equation (45), is given as(45)Fz,t=c1exp−γt¯c2z¯J1121−4P;γz¯+c3z¯J2121−4P;γz¯×1cosK/2texp−αt+k2×exp−αlnz−z−2Kβ2tanK2tz.
Incorporating the values ofz¯ and t¯ from equation (33) in equation (45), the solution Fz,t in equation (45) takes the form(46)Fz,t=c1exp−γk12KtanK2t×c2±8zk1cosK/2tJ1121−4P;γ±8zk1cosK/2t+c3±8zk1cosK/2tJ2121−4P;γ±8zk1cosK/2t×1cosK/2texp−αt+k2×exp−αlnz−z−2Kβ2tanK2tz.


## 2.1. Case 1: Reduction to Heat Equation—Ft=Fzz

First, find the equivalence transformations which reduce equation (4) into different canonical forms. With the comparison of PDE (4) with scalar linear parabolic PDE (8), the coefficients functions are given by(19)At,z=−12β2z,Bt,z=−μα−z,Ct,z=z.With the use of Theorem1, the parametric values are deduced such that the PDE (4) is mapped to the heat equation. The value of Π given in equation (14) for the CIR equation is found as(20)Π=μ2z−α22β2z+μz−α2z+332β2z−μ2+z.It is clear thatm and n are equal to zero in the condition defined in Theorem 1 because the coefficients of PDE (4) are independent of time t. Under the value of Π, condition (12) takes the form(21)lz=z3−38β4+2μαβ2−2μ2α2.From the above equation, one obtains the following parametric conditions:(22)β1=−4μα3,β2=4μα3,β3=−4μα,β4=4μα.For specific values of the parameterβ, equation (22) gives the nontrivial cases such that invariant condition (12) is holds, which reduces the CIR PDE (4) to the heat equation.Now, the transformation relations will be computed which give reduction of the CIR equation (4) to the heat equation. These transformations are further used to find the fundamental solutions of the PDE (4).With the use of the coefficients as defined in Theorem3, one obtains the following results:(23)ut=−18β2+μ22,vt=0,wt=μ2αβ2,χt=k12KtanK2t,σt=k2∫dtlncosK/2t,ft=exp12lncosK/2t−k324∫dtlncosK/2t2+t+k4,where K=8ft and kjj=1,2,3,4 are the integration constants and β can take any value as given in equation (22). Substituting the values from equation (23) to equation (16), the equivalence transformations which give the reduction of CIR equation (4) into the heat equation (after performing some lengthy and tedious calculation) are given as(24)t¯=k12/Ktan2/Kt,z¯=±4k1cos2/Ktz+k2∫tTdtlncos2/Kt+k3,F¯=12β2z−1/4exp12lncos2/Kt−k424∫tTdtlncos2/Kt2+k5×Fz,texpμβ2αlnz−z±2k3lncos2/Ktz,with kjj=1,…,5 are integration constants. The terminal condition (5) for the CIR equation is also transformed to(25)Fz,T=12β2z−1/4exp12lncosK/2t+μ2αβ2T+k5×expμβ2αlnz−z±2k3lncosK/2tz.In order to seek for the fundamental solutions of the PDE (4), it is seen that there exist equivalence transformations (24) which gives reduction of PDE (4) to the heat equation. This result is employed to find the fundamental solution of the CIR equation.It is well known that the fundamental solution of the heat is given in barred coordinates [39]:(26)F¯z¯,t¯=12t¯πexp−z¯24t¯.To look for the fundamental solution of CIR equation (4), solution (26) is transformed by means of (24). With the use of equations (24) and (26), the solution Fz,t is given by(27)Fz,t=12β2z1/412t¯πexp−z24t¯×exp−12lncosK/2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz.Finally, using the values oft¯ and z¯ into equation (27), the above solution takes the form(28)Fz,t=12β2z1/412k1K/2tanK/2tπ×exp−±4k1/cosK/2tz+k2∫tTdt/lncosK/2t+k324k1K/2tanK/2t×exp−12lncosK2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz.Under the same set of equivalence transformations for reduction into the heat equation, one can find another closed-form solution for PDE (4). The fundamental solution of the heat equation is also given as series solution of the form [37](29)F¯z¯,t¯=∑r=1∞prexp−rπL2qt¯sinrπz¯L,which converges for all 0≤x≤L. Again, one needs to transform the solution given in equation (29) subject to the transformations given in equation (24). This solution is written as(30)Fz,t=12β2z1/4∑r=1∞prexp−rπL2qt¯sinrπz¯L×exp−12lncosK2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz.Substituting the values ofz¯ and t¯, the series solution Fz,t takes the form(31)Fz,t=12β2z1/4∑r=1∞prexp−rπL2qk12KtanK2t×∑r=1∞sinrπ±4k1/cosK/2tz+k2∫tTdt/lncosK/2t+k3L×exp−12lncosK2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz.


## 2.2. Case 2: Reduction to Second Lie Canonical Equation—Ft=Fzz+A/z2Fz,t

We look for the case when the CIR PDE (4) is transformed to the second Lie canonical equation. In this situation, one deduce following parametric constraint(32)β2≠−8μα±4μα−3⇒lz≠0.Under the same procedure used previously, the equivalence transformations that transform equation (4) into the second Lie canonical form are given by(33)t¯=k12KtanK2t,z¯=±8k1cosK/2tz,F¯=cosK/2texpαt+k2expαlnz−z−2Kβ2tanK/2tzF,with k1 and k2 are the integration constants.Now, to deduce the closed-form solution such that the CIR equation is transformed into the second Lie canonical equation subject to the equivalence transformations given in equation (33), first, rewrite the second canonical equation:(34)F¯t¯−F¯z¯z¯−Pz¯2F¯=0.In order to find the exact solution of CIR equation, one first needs to seek the solution of PDE (34) and then employ transformations (33). To obtain the solution for PDE (34), the separation of variables technique is used.Let us assume that the PDE (34) has the separable solution of the form(35)F¯z¯,t¯=Uz¯Vt¯.Using equation (35) in equation (34) and separating the variables, one obtains(36)Vt¯+γV=0,(37)Uz¯z¯+Pz¯2+γU=0,with γ≥0 is the separation constant. The following two independent cases are considered here.Case I.
γ=0
The solutions of ODEs (36) and (37) are written as(38)Vt¯=c1,(39)Uz¯=c2z¯1/21−1−4P+c3x¯−1/21−1−4P,Subject to equations (38) and (39), and the exact solution of equation (34) is written as(40)F¯z¯,t¯=c1c2z¯1/21−1−4P+c3x¯−1/21−1−4P.
In order to obtain the solution of PDE (4), solution (40) is transformed subject to the transformation relations given in equation (33). One can write the solution Fz,t as(41)Fz,t=c1c2z¯1/21−1−4P+c3x¯−1/21−1−4P,×1cosK/2texp−αt+k2.
Substituting the values ofz¯ and t¯ from equation (33) into equation (41), the solution for Fz,t takes the form(42)Fz,t=c1c2±8k1cosK/2tz1/21−1−4P+c3±8k1cosK/2tz−1/21−1−4P×exp−αlnz−z−2Kβ2tanK2tz×1cosK/2texp−αt+k2.Case II.
γ>0
Here, the solutions of the ODEs (36) and (37) are found as(43)Vt¯=c1exp−γt¯,Uz¯=c2z¯J1121−4P;γz¯+c3z¯J2121−4P;γz¯,with J1 and J2 as the first and the second kind Bessel functions [25], respectively. Therefore, the solution of equation (34) is given by(44)F¯z¯,t¯=c1exp−γt¯c2z¯J1121−4P;γz¯+c3z¯J2121−4P;γz¯.
The solution of CIR PDE (4), subject to equation (45), is given as(45)Fz,t=c1exp−γt¯c2z¯J1121−4P;γz¯+c3z¯J2121−4P;γz¯×1cosK/2texp−αt+k2×exp−αlnz−z−2Kβ2tanK2tz.
Incorporating the values ofz¯ and t¯ from equation (33) in equation (45), the solution Fz,t in equation (45) takes the form(46)Fz,t=c1exp−γk12KtanK2t×c2±8zk1cosK/2tJ1121−4P;γ±8zk1cosK/2t+c3±8zk1cosK/2tJ2121−4P;γ±8zk1cosK/2t×1cosK/2texp−αt+k2×exp−αlnz−z−2Kβ2tanK2tz.


## 3. Fundamental Solution of Cauchy Initial-Value Problem: A Case Study

This section is devoted to discuss the fundamental solution of the Cauchy initial-value problem for the CIR model subject to the terminal condition. Since the transformation relations are already deduced for PDE (4) for reduction to heat equation, these transformations are used to discuss the fundamental solutions of the corresponding initial-value problem.The standard Cauchy initial-value problem for the heat equation is given by [40](47)F¯t¯=F¯z¯z¯,F¯z¯,T=Ωz¯,with some function Ωz¯. The solution of the above initial-value problem is given by [38](48)F¯z¯,t¯=12t¯π∫−∞+∞Ωηexp−z¯−η24t¯dη.One needs to transform solution (48) into the solution of CIR equation (4). Inserting F¯z¯,t¯ from equation (48) to equation (33), we deduce(49)Fz,t=12β2z−1/4exp−12lncosK2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz×12t¯π∫−∞+∞Ωηexp−z¯−η24t¯dη.Keeping in mind condition (25), the function Ωη in equation (49) is written as(50)Ωη=12β2z−1/4exp12lncosK2t+μ2αβ2T+k5×expμβ2αlnz−z±2k3lncosK/2tz.Finally, by substituting of the values ofz¯ and t¯ from equation (24) to equation (49), one can write(51)Fz,t=12β2z−1/4exp−12lncosK2t−k424∫tTdtlncosK/2t2+k5×exp−μβ2αlnz−z±2k3lncosK/2tz×12πk1K/2tanK/2t×∫−∞+∞Ωηexp−±4k1/cosK/2tz+k2∫tTdt/lncosK/2t+k3−η24k1K/2tanK/2tdη.The value of the integral in the above solution could be found in terms of the Kummer confluent hypergeometric functions [41] or can be computed numerically.


## 4. Conservation Laws of Cox–Ingersoll–Ross Equation

The conservation laws for the CIR equation (4) are derived in this section. These conservation laws are constructed for the first time by using the conservation theorem of Ibragimov [38]. The CIR equation admits six plus infinite set of Lie symmetries [11]. Every symmetry generator of the CIR equation gives rise to a conservation law.The formal second-order Lagrangian for the CIR equation (1) is written as(52)ℒ=GFt+12β2zFzz+μα−zFz−zF.The Lagrangian defined in (52) is used to derive the adjoint equation of the CIR equation, which is given as(53)Gt+2z+β22G−12zβ2zGzz−2μα−zGz=0.Every Lie symmetry generator of the CIR equation of the form(54)Γ=η1∂t+η2∂z+η3∂F,gives rise to a conservation law(55)DtWt+DzWz=0,with(56)Wt=η1ℒ+CℒFt,Wz=η2ℒ+CℒFz−DzℒFzz+DzCℒFzz,where(57)C=η3−η1Ft+η2Fz.The symmetries of the CIR equation (4) and the associated adjoint equation can be computed easily. Firstly, we consider the time-translation symmetry Γ1=∂t of CIR equation (4). By using the Ibragimov theorem [38], the components of the conserved vector associated with operator Γ1 are given by(58)W1t=FGt,W1z=z2−β2zFGt+12β2z2FzGt−FGtz.Similarly, if one uses the symmetry operatorΓ2=z∂z, one can deduce the associated conservation law with components:(59)W2t=−zFzG,W2z=zFt−μαz−z2F+12β2z2FzG+12β2z3Gz.


## 5. Final Remarks

The purpose of this work is to demonstrate the utility of an algorithmic method, an invariant approach based on the theory of Lie groups. Utilizing this invariant technique, the valuation parabolic PDE of interest rate derivative in the CIR model framework is studied in detail. The classification of the CIR equation was based on the classical results of Lie’s classification theory for the parabolic PDEs. Hence, the invariant theory helped in the computation of equivalence transformations which transform the underlying model to different Lie canonical equations. It was found that, for proper choice of parametric values, the CIR equation was transformed to the first and the second Lie canonical equations. For each case, the transformation relations are derived which reduced the CIR equation into the first and the second Lie canonical form. These transformations were further employed to find the invariant solutions of the CIR equation by using some well-known fundamental solutions of the first and the second Lie canonical equations. Moreover, the Cauchy initial-value problem of the CIR model, with a suitable choice of a terminal condition, is also studied in detail. The solution technique was based on translating the Cauchy initial-value problem for the CIR model into the classical heat equation Cauchy problem and then implementing the well-known Cauchy solutions of the heat equation.The conservation laws corresponding to the two Lie point symmetries of the CIR equation are constructed by using a general theorem on conservation laws. The notion of conservation laws and associated Lie point symmetries can be used to formulate exact solutions of such type of complicated equations arising in the study of financial mathematics models.Lie’s group theory provides a systematic way to the resolution of the various DEs which arise in financial mathematics. It should come as no surprise that many evolution PDEs in the field of mathematical finance have a rich symmetry group structure. These symmetry structures help to explain the ease with which ansatz-based techniques can yield solutions. However, such nonalgorithmic adhoc procedures have no place, especially when dealing with real-world problems that are described in terms of complicated PDEs. It is important to note that the invariant analysis performed on the CIR equation can be carried out on any other class of parabolic models which transform not only to the heat equation but also to other canonical equations.



---

*Source: 9974073-2021-06-15.xml*