Journal of Applied Mathematics

Volume 2013 (2013), Article ID 276238, 9 pages

http://dx.doi.org/10.1155/2013/276238

## Lie-Algebraic Approach for Pricing Zero-Coupon Bonds in Single-Factor Interest Rate Models

Institute of Theoretical Physics and Department of Physics, The Chinese University of Hong Kong, Shatin, New Territories, Hong Kong

Received 17 December 2012; Revised 9 April 2013; Accepted 11 April 2013

Academic Editor: Alvaro Valencia

Copyright © 2013 C. F. Lo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The Lie-algebraic approach has been applied to solve the bond pricing problem in single-factor interest rate models. Four of the popular single-factor models, namely, the Vasicek model, Cox-Ingersoll-Ross model, double square-root model, and Ahn-Gao model, are investigated. By exploiting the dynamical symmetry of their bond pricing equations, analytical closed-form pricing formulae can be derived in a straightfoward manner. Time-varying model parameters could also be incorporated into the derivation of the bond price formulae, and this has the added advantage of allowing yield curves to be fitted. Furthermore, the Lie-algebraic approach can be easily extended to formulate new analytically tractable single-factor interest rate models.

#### 1. Introduction

In this paper we apply the Lie-algebraic method to tackle the bond pricing problem in single-factor interest rate models. In particular, we investigate the bond pricing equation of the form where is a real parameter, is a real function of both spot interest rate and time-to-maturity , is the time-varying volatility, and denotes the price of a zero-coupon bond of duration with a value of unity at maturity; that is, . By exploiting the dynamical symmetry of the bond pricing equation, we derive analytically tractable single-factor interest rate models in a unified manner and obtain their closed-form bond pricing formulae. It is found that not only four of the popular single-factor models, namely the Vasicek model [1], Cox-Ingersoll-Ross model [2], double square-root model [3], and Ahn-Gao model [4], can be derived in a straightforward manner, but also new analytically tractable models can be formulated systematically. Moreover, time-varying model parameters can be incorporated into the derivation of the bond pricing formulae without difficulty. This has the added advantage of allowing yield curves to be fitted, and thus a “no-arbitrage” yield curve model can be developed to match the current market data.

The Lie-algebraic method was introduced by Lo and Hui [5–7] to the field of finance for the pricing of financial derivatives with time-dependent model parameters. This new method is very simple and consists of two basic ingredients: identifying the dynamical symmetries of the given pricing partial differential equations and applying the Wei-Norman theorem [8] to solve the equations and obtain analytical closed-form pricing formulae. For demonstration, the Lie-algebraic approach has already been applied to price European options for the constant elasticity of variance processes and corporate discount bonds with default risk, multiasset financial derivatives, and so forth. It should be noted that the Lie-algebraic method is different from the Lie group analysis which was introduced by Ibragimov and his coauthors [9, 10] to tackle partial differential equations occurring in financial problems. The Lie group analysis is a mathematical theory developed by Sophus Lie and classifies partial differential equations in terms of their symmetry groups, thereby identifying the set of equations which could be integrated or reduced to low-order equations by group theoretic algorithms. (Details of the Lie group analysis and its application to partial differential equations can be found in, for example, Hydon (2000) and Bluman et al. (2010) [11, 12].) Further applications of the Lie group analysis in mathematical finance were subsequently explored by a number of papers, for example, Goard [13], Pooe et al. [14], Goard et al. [15], Leach et al. [16], and Sinkala et al. [17]. A recent review of the applications of the Lie theory to problems in mathematical finance and economics can be found in [18].

#### 2. Lie-Algebraic Approach

We consider a possible set of differential operators realizing the Lie algebra [19]: where , , and are real adjustable parameters (see Appendix A). Then we try to look for appropriate linear combinations of these operators, namely, , where the coefficients are arbitrary scalar functions of only, to produce the bond pricing equation given in (1). Here are some illustrative examples.

For , , , , , , and being arbitrary, we recover the bond pricing equation of the *Vasicek model:*
By applying the Wei-Norman theorem, we can derive the bond price as (see Appendix B)
where
As a result, the bond price can be expressed as
In the special case of constant model parameters, that is, , , and , the and can be analytically determined as
and the bond price is reduced to the well-known closed-form expression [1]:

For , , , , , , , and being arbitrary, the bond pricing equation of the *double square-root model* is reproduced:
As in the Vasicek model, we apply the Wei-Norman theorem to determine the bond price as (see Appendix B)
where
Accordingly, the bond price is given by
In the special case of constant model parameters, that is, , , and , the and can be analytically determined as
with and . Consequently, we are able to recover the well-known closed-form expression of the bond price [3]:
where
with

For , , , , , and , a bond pricing equation with a special time-dependent nonlinear drift term can be obtained: which can be straightforwardly solved as in the Vasicek model. As a result, the bond price can be expressed as (see Appendix B) where

For , , , , , and , we can derive the bond pricing equation: which has both the dependence of volatility and a time-dependent nonlinear drift term. By performing the same analysis as in the other three cases, the bond price is found to be given by (see Appendix B) where

Next we apply the same analysis to derive the Cox-Ingersoll-Ross model and Ahn-Gao model from an alternative set of differential operators realizing the subalgebra of the Lie algebra : where , , and are real adjustable parameters (see Appendix A). The subalgebra is actually the reductive Lie algebra .

By choosing , , and , the bond pricing equation of the *Cox-Ingersoll-Ross model* with constant model parameters can be expressed in terms of the differential operators realizing the subalgebra as
(Strictly speaking, the model parameters , , and could be time-varying with the constraint that is independent of time .)

By the Wei-Norman theorem, the bond price can be easily determined as (see Appendix B) where The Riccati equation with constant coefficients in (26) can be straightforwardly solved to yield with . Once the has been found, we are also able to obtain via analytical integrations. Beyond question, our finding is in agreement with the well-known closed-form result [2]:

By setting , , and , we can cast the bond pricing equation of the *Ahn-Gao model,* which includes nonlinearity in the drift term and a realistic dependence in the volatility (not only the nonlinear drift term of the Ahn-Gao model is consistent with the empirical findings of Aït-Sahalia [20], but also the chosen dependence of volatility is the best fit power law for volatility [21, 22]), in terms of the differential operators realizing the subalgebra into the form
where is a real function of . Then applying the Wei-Norman theorem allows us to represent the bond price by (see Appendix B)
where and
Without loss of generality, we suppose that , where ,
and . It is not difficult to show that the bond price is given by
where and
The function is the Bessel function of the first kind of order . Here we have made use of the fact that is an eigenfunction of the operator with the eigenvalue . The integral over can be analytically evaluated to give [23]
for , , , and . The function is the modified Bessel function of the first kind of order . As a result, is found to be given by
Since , we can readily derive the bond price as follows:
where , denotes the Gamma function and is the standard confluent hypergeometric function [23, 24]. Furthermore, (38) will reproduce the well-known closed-form result if the model parameter is independent of time [4].

#### 3. Conclusion

In this paper the Lie-algebraic method has been applied to solve the bond pricing problem in single-factor interest rate models. Four of the popular single-factor models, namely the Vasicek model, Cox-Ingersoll-Ross model, double square-root model, and Ahn-Gao model, are investigated, and analytical closed-form pricing formulae are derived. Since all the four bond pricing equations exhibit the dynamical symmetry or its subgroup, their solutions can be derived in a unified manner and have very similar mathematical structures. This interesting feature helps shed new light upon the systematic formulation of new analytically tractable single-factor interest rate models, as demonstrated in Section 2. Time-varying model parameters could also be incorporated into the derivation of the bond price formulae without difficulty. This has the added advantage of allowing yield curves to be fitted, and thus a “no-arbitrage” yield curve model can be developed to match the current market data. Hence, we believe that the Lie-algebraic method will provide an easy-to-use analytical tool for the bond pricing problem. Furthermore, the Lie-algebraic approach can be easily extended to the pricing of other standard European interest rate derivatives for they differ from the zero-coupon bonds in the final payoff conditions only [25].

#### Appendices

#### A. Generators of the Lie Algebra and Its Subalgebras

The generators of the Heisenberg-Weyl Lie algebra obey the set of commutation relations [19]: By direct substitution, a possible set of differential operators realizing the Lie algebra can be identified as where , , and are real adjustable parameters. Then, in terms of these generators one can construct the generators of the Lie algebra as follows: which satisfy the set of commutation relations [19]: These six generators in turn form the Lie algebra , which is defined by the following set of commutation relations [19]: In addition to the subalgebras and , we may also form another subalgebra in terms of the four generators satisfying the commutation relations: The subalgebra is actually the reductive Lie algebra . Moreover, it is not difficult to show that an alternative realization of the set of generators of the Lie algebra is given by where , , and are real adjustable parameters.

#### B. Wei-Norman Theorem

Consider the linear operator differential equation of the first order where and are both time-dependent linear operators in a Banach space or a finite-dimensional space. According to the Wei-Norman theorem [8], if the operator can be expressed as where ’s are scalar functions of time and are the generators of an -dimensional solvable Lie algebra or a real split 3-dimensional simple Lie algebra, then the operator can assume the following form: Here the ’s are time-dependent scalar functions to be determined. To find the ’s, we simply substitute (B.2) and (B.3) into (B.1) and compare the two sides term by term to obtain a set of coupled nonlinear differential equations where are nonlinear functions of ’s. Thus, we have transformed the linear operator differential equation in (B.1) to a set of coupled nonlinear differential equations of scalar functions in (B.4).

Moreover, according to *Levi’s Theorem*, “If is a finite-dimensional Lie algebra with radical which is the maximal solvable ideal of the Lie algebra, then there exists a semisimple subalgebra of such that is the semidirect sum ,” in the equation , where generates , the decomposition gives rise to the corresponding decomposition , where and . Then it is easy to verify that where and satisfy
Since is an ideal in , we can easily see that is in . The fact that is solvable makes it easy to find once has been found. More details can be found in [8].

For illustration, we apply the Wei-Norman theorem to the following cases.

##### B.1. Heisenberg-Weyl Lie Algebra

The Heisenberg-Weyl Lie algebra is defined by the set of commutation relations [17]: of its generators . Given that the Wei-Norman theorem states that can be expressed as where the time-dependent functions ’s satisfy a set of three coupled nonlinear differential equations: It is obvious that the set of differential equations can be easily solved by quadrature: As a result, the operator is thus determined.

##### B.2. Lie Algebra

We consider the evolution equation of the operator : where the operators form the Lie algebra defined by the commutation relations [19]: According to the Wei-Norman theorem, the operator can be expressed in the product form where the time-dependent functions ’s satisfy a set of three coupled nonlinear differential equations: After further simplification, these three differential equations become Hence, once the is found by solving the Riccati equation, the and can be readily determined by quadrature.

##### B.3. Lie Algebra

If is a linear combination of the six generators of the Lie algebra , then, according to Levi’s theorem, we may decompose the into two parts and , and the operator assumes the product form where and satisfy (B.5) and (B.6), respectively. It is obvious that the operator is given by where Next, in order to determine , we need to evaluate . Using the explicit form of the operator , we can apply the Baker-Hausdorff formula [26] to derive the operator : where Then, the operator can be easily found to be given by where

#### References

- O. Vasicek, “An equilibrium characterization of the term structure,”
*Journal of Financial Economics*, vol. 5, no. 2, pp. 177–188, 1977. View at Google Scholar · View at Scopus - J. C. Cox, J. E. Ingersoll,, and S. A. Ross, “An intertemporal general equilibrium model of asset prices,”
*Econometrica*, vol. 53, no. 2, pp. 363–384, 1985. View at Publisher · View at Google Scholar · View at MathSciNet - F. A. Longstaff, “A nonlinear general equilibrium model of the term structure of interest rates,”
*Journal of Financial Economics*, vol. 23, no. 2, pp. 195–224, 1989. View at Google Scholar · View at Scopus - D. H. Ahn and B. Gao, “A parametric nonlinear model of term structure dynamics,”
*Review of Financial Studies*, vol. 12, no. 4, pp. 721–762, 1999. View at Google Scholar · View at Scopus - C. F. Lo and C. H. Hui, “Valuation of financial derivatives with time-dependent parameters: lie-algebraic approach,”
*Quantitative Finance*, vol. 1, no. 1, pp. 73–78, 2001. View at Publisher · View at Google Scholar · View at MathSciNet - C. F. Lo and C. H. Hui, “Pricing multi-asset financial derivatives with time-dependent parameters—lie algebraic approach,”
*International Journal of Mathematics and Mathematical Sciences*, vol. 32, no. 7, pp. 401–410, 2002. View at Publisher · View at Google Scholar · View at Scopus - C. F. Lo and C. H. Hui, “Lie-algebraic approach for pricing moving barrier options with time-dependent parameters,”
*Journal of Mathematical Analysis and Applications*, vol. 323, no. 2, pp. 1455–1464, 2006. View at Publisher · View at Google Scholar · View at Scopus - J. Wei and E. Norman, “Lie algebraic solution of linear differential equations,”
*Journal of Mathematical Physics*, vol. 4, no. 4, pp. 575–581, 1963. View at Google Scholar · View at Scopus - N. H. Ibragimov and C. W. Soh, “Solution of the Cauchy problem for the Black-Scholes equation using its symmetries,” in
*Proceedings of the International Conference on Modern Group Analysis*, N. H. Ibragimov, K. R. Naqvi, and E. Straume, Eds., MARS, Nordfjordeid, Norway, 1997. - R. K. Gazizov and N. H. Ibragimov, “Lie symmetry analysis of differential equations in finance,”
*Nonlinear Dynamics*, vol. 17, no. 4, pp. 387–407, 1998. View at Google Scholar · View at Scopus - P. E. Hydon,
*Symmetry Methods for Differential Equations: A Beginner's Guide*, Cambridge University, New York, NY, USA, 2000. View at Publisher · View at Google Scholar · View at MathSciNet - G. W. Bluman, A. F. Cheviakov, and S. C. Anco,
*Applications of Symmetry Methods to Partial Differential Equations*, Springer, New York, NY, USA, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - J. Goard, “New solutions to the bond-pricing equation via Lie's classical method,”
*Mathematical and Computer Modelling*, vol. 32, no. 3-4, pp. 299–313, 2000. View at Publisher · View at Google Scholar · View at Scopus - C. A. Pooe, F. M. Mahomed, and C. W. Soh, “Fundamental solutions for zero-coupon bond pricing models,”
*Nonlinear Dynamics*, vol. 36, no. 1, pp. 69–76, 2004. View at Publisher · View at Google Scholar · View at Scopus - J. Goard, P. Broadbridge, and G. Raina, “Tractable forms of the bond pricing equation,”
*Mathematical and Computer Modelling*, vol. 40, no. 1-2, pp. 151–172, 2004. View at Publisher · View at Google Scholar · View at Scopus - P. G. L. Leach, J. G. O'Hara, and W. Sinkala, “Symmetry-based solution of a model for a combination of a risky investment and a riskless investment,”
*Journal of Mathematical Analysis and Applications*, vol. 334, no. 1, pp. 368–381, 2007. View at Publisher · View at Google Scholar · View at Scopus - W. Sinkala, P. G. L. Leach, and J. G. O'Hara, “Zero-coupon bond prices in the Vasicek and CIR models: their computation as group-invariant solutions,”
*Mathematical Methods in the Applied Sciences*, vol. 31, no. 6, pp. 665–678, 2008. View at Publisher · View at Google Scholar · View at Scopus - I. Hernández, C. Mateos, J. Núñez, and Á. F. Tenorio, “Lie theory: applications to problems in mathematical finance and economics,”
*Applied Mathematics and Computation*, vol. 208, no. 2, pp. 446–452, 2009. View at Publisher · View at Google Scholar · View at Scopus - A. Perelomov,
*Generalized Coherent States and Their Applications*, Springer, New York, NY, USA, 1986. View at MathSciNet - Y. Aït-Sahalia, “Testing continuous-time models of the spot interest rate,”
*Review of Financial Studies*, vol. 9, no. 2, pp. 385–426, 1996. View at Google Scholar · View at Scopus - K. C. Chan, G. A. Karolyi, F. Longstaff, and A. B. Sanders, “An empirical comparison of alternative models of the short-term interest rate,”
*Journal of Finance*, vol. 47, pp. 1209–1227, 1992. View at Google Scholar - J. Y. Campbell, A. W. Lo, and A. C. MacKinlay,
*The Econometrics of Financial Markets*, Princeton University, Princeton, NJ, USA, 1996. - I. S. Gradshteyn and I. M. Ryzhik,
*Table of Integrals, Series, and Products*, Academic Press, Boston, Mass, USA, 5th edition, 1994. View at MathSciNet - N. N. Lebedev,
*Special Functions and Their Applications*, Dover Publications, New York, NY, USA, 1972. View at MathSciNet - P. Wilmott, S. Howison, and J. Dewynne,
*The Mathematics of Financial Derivatives: A Student Introduction*, Cambridge University, New York, NY, USA, 1989. - W. H. Louisell,
*Quantum Statistical Properties of Radiation*, John Wiley & Sons, New York, NY, USA, 1990. View at MathSciNet