Abstract
This paper develops an eigenfunction expansion approach to solve discretely monitored first passage time problems for a rich class of Markov processes, including diffusions and subordinate diffusions with jumps, whose transition or Feynman–Kac semigroups possess eigenfunction expansions in \(L^{2}\)-spaces. Many processes important in finance are in this class, including OU, CIR, (JD)CEV diffusions and their subordinate versions with jumps. The method represents the solution to a discretely monitored first passage problem in the form of an eigenfunction expansion with expansion coefficients satisfying an explicitly given recursion. A range of financial applications is given, drawn from across equity, credit, commodity, and interest rate markets. Numerical examples demonstrate that even in the case of frequent barrier monitoring, such as daily, approximating discrete first passage time problems with continuous solutions may result in unacceptably large errors in financial applications. This highlights the relevance of the method to financial applications.
Similar content being viewed by others
References
Abramowitz, M., Stegun, I.A.: Handbook of Mathematical Functions. Dover, New York (1972)
Ahn, D.H., Gao, B.: A parametric nonlinear model of term structure dynamics. Rev. Financ. Stud. 12, 721–762 (1999)
Albanese, C., Kuznetsov, A.: Unifying the three volatility models. Risk 17(3), 94–98 (2004)
Alili, L., Kyprianou, A.E.: Some remarks on first passage of Lévy processes, the American put and pasting principles. Ann. Appl. Probab. 15, 2062–2080 (2005)
Alili, L., Patie, P., Pedersen, J.L.: Representations of the first hitting time density of an Ornstein–Uhlenbeck process. Stoch. Models 21, 967–980 (2004)
Barndorff-Nielsen, O.E.: Processes of normal inverse Gaussian type. Finance Stoch. 2, 41–68 (1998)
Barndorff-Nielsen, O.E., Levendorskiĭ, S.: Feller processes of normal inverse Gaussian type. Quant. Finance 1, 318–331 (2001)
Beaglehole, D.R., Tenney, M.: A nonlinear equilibrium model of the term structure of interest rates: corrections and additions. J. Financ. Econ. 32, 345–353 (1992)
Borodin, A., Salminen, P.: Handbook of Brownian Motion—Facts and Formulae. Birkhäuser, Basel (2002)
Boyarchenko, M., Levendorskiĭ, S.: Valuation of continuously monitored double barrier options and related securities. Math. Finance 22, 419–444 (2012)
Boyarchenko, N., Levendorskiĭ, S.: The eigenfunction expansion method in multifactor quadratic term structure models. Math. Finance 17, 503–539 (2007)
Boyarchenko, S., Levendorskiĭ, S.: Barrier options and touch-and-out options under regular Lévy processes of exponential type. Ann. Appl. Probab. 12, 1261–1298 (2002)
Brigo, D., Mercurio, F.: Interest Rate Models—Theory and Practice, 2nd edn. Springer, Berlin (2006)
Broadie, M., Glasserman, P., Kou, S.G.: A continuity correction for discrete barrier options. Math. Finance 7, 325–348 (1997)
Broadie, M., Glasserman, P., Kou, S.G.: Connecting discrete and continuous path-dependent options. Finance Stoch. 3, 55–82 (1999)
Broadie, M., Yamamoto, Y.: A double-exponential fast Gauss transform algorithm for pricing discrete path-dependent options. Oper. Res. 53, 764–779 (2005)
Carr, P., Linetsky, V.: A jump to default extended CEV model: an application of Bessel processes. Finance Stoch. 10, 303–330 (2006)
Chen, A., Feng, L., Song, R.: On the monitoring error of the supremum of a normal jump-diffusion process. J. Appl. Probab. 48, 1021–1034 (2011)
Collin-Dufresne, P., Goldstein, R.S.: Do credit spreads reflect stationary leverage ratios? J. Finance 56, 1929–1957 (2001)
Cox, J.C.: Notes on option pricing I: constant elasticity of variance diffusions. Working paper, Stanford University (1975) (Reprinted in J. Portf. Manag. 22, 15–17 (1996))
Cox, J.C., Ingersoll, J.E., Ross, S.A.: A theory of the term structure of interest rates. Econometrica 53, 385–407 (1985)
Davies, E.B.: Linear Operators and Their Spectra. Cambridge University Press, Cambridge (2007)
Davydov, D., Linetsky, V.: Pricing and hedging path-dependent options under the CEV process. Manag. Sci. 47, 949–965 (2001)
Davydov, D., Linetsky, V.: Pricing options on scalar diffusions: an eigenfunction expansion approach. Oper. Res. 51, 185–209 (2003)
De Innocentis, M., Levendorskiĭ, S.: Pricing discrete barrier options and credit default swaps under Lévy processes. Quant. Finance 14, 1337–1365 (2014)
Dia, E.H.A., Lamberton, D.: Connecting discrete and continuous lookback or hindsight options in exponential Lévy models. Adv. Appl. Probab. 43, 1136–1165 (2011)
Dia, E.H.A., Lamberton, D.: Continuity correction for barrier options in jump-diffusion models. SIAM J. Financ. Math. 2, 866–900 (2011)
Feng, L., Linetsky, V.: Pricing discretely monitored barrier options and defaultable bonds in Lévy process models: a fast Hilbert transform approach. Math. Finance 18, 337–384 (2008)
Feng, L., Linetsky, V.: Computing exponential moments of the discrete maximum of a Lévy process and lookback options. Finance Stoch. 13, 501–529 (2009)
Fukushima, M., Oshima, Y., Takeda, M.: Dirichlet Forms and Symmetric Markov Processes, 2nd edn. de Gruyter, Berlin (2011)
Fulton, C., Pruess, S.: Eigenvalue and eigenfunction asymptotics for regular Sturm–Liouville problems. J. Math. Anal. Appl. 188, 297–340 (1994)
Fusai, G., Abrahams, I.D., Sgarra, C.: An exact analytical solution for discrete barrier options. Finance Stoch. 10, 1–26 (2006)
Fusai, G., Recchioni, M.C.: Analysis of quadrature methods for pricing discrete barrier options. J. Econ. Dyn. Control 31, 826–860 (2007)
Geman, H.: Commodities and Commodity Derivatives: Modeling and Pricing for Agriculturals, Metals and Energy. Wiley, Hoboken (2005)
Göing-Jaeschke, A., Yor, M.: A clarification about hitting time densities for Ornstein–Uhlenbeck processes. Finance Stoch. 7, 413–415 (2003)
Gorovoi, V., Linetsky, V.: Black’s model of interest rates as options, eigenfunction expansion and Japanese interest rates. Math. Finance 14, 49–78 (2004)
Hörfelt, P.: Extension of the corrected barrier approximation by Broadie, Glasserman and Kou. Finance Stoch. 7, 231–243 (2003)
Howison, S., Steinberg, M.: A matched asymptotic expansions approach to continuity corrections for discretely sampled options. Part 1: barrier options. Appl. Math. Finance 14, 63–89 (2007)
Karlin, S., Taylor, H.M.: A Second Course in Stochastic Processes. Academic Press, New York (1981)
Keilson, J., Ross, H.F.: Passage time distributions for Gaussian Markov (Ornstein–Uhlenbeck) statistical processes. In: Selected Tables in Mathematical Statistics, vol. III, pp. 233–259. Am. Math. Soc., Providence (1975)
Koekoek, R., Lesky, P.A., Swarttouw, R.F.: Hypergeometric Orthogonal Polynomials and Their \(q\)-Analogues. Springer, Berlin (2010)
Kou, S.G.: On pricing of discrete barrier options. Stat. Sin. 13, 955–964 (2003)
Kou, S.G., Wang, H.: First passage times of a jump diffusion process. Adv. Appl. Probab. 35, 504–531 (2003)
Kuan, C.H., Webber, N.: Pricing barrier options with one-factor interest rate models. J. Deriv. 10(4), 33–50 (2003)
Kuznetsov, A., Kyprianou, A.E., Pardo, J.C.: Meromorphic Lévy processes and their fluctuation identities. Ann. Appl. Probab. 22, 1101–1135 (2012)
Leblanc, B., Renault, O., Scaillet, O.: A correction note on the first passage time of an Ornstein–Uhlenbeck process to a boundary. Finance Stoch. 4, 109–111 (2000)
Leblanc, B., Scaillet, O.: Path-dependent options on yield in the affine term structure model. Finance Stoch. 2, 349–367 (1998)
Leippold, M., Wu, L.: Asset pricing under the quadratic class. J. Financ. Quant. Anal. 37, 271–295 (2002)
Li, J., Li, L., Mendoza-Arriaga, R.: Additive subordination and its applications in finance (2014), preprint, available at SSRN, http://ssrn.com/abstract=2470859
Li, L., Linetsky, V.: Optimal stopping and early exercise: an eigenfunction expansion approach. Oper. Res. 61, 625–643 (2013)
Li, L., Linetsky, V.: Time-changed Ornstein–Uhlenbeck processes and their applications in commodity derivative models. Math. Finance 24, 289–330 (2014)
Li, L., Mendoza-Arriaga, R.: Ornstein–Uhlenbeck processes time changed with additive subordinators and their applications in commodity derivative models. Oper. Res. Lett. 41, 521–525 (2013)
Lim, D., Li, L., Linetsky, V.: Evaluating callable and putable bonds: an eigenfunction expansion approach. J. Econ. Dyn. Control 36, 1888–1908 (2012)
Linetsky, V.: Computing hitting time densities for CIR and OU diffusions: applications to mean reverting models. J. Comput. Finance 7, 1–22 (2004)
Linetsky, V.: The spectral decomposition of the option value. Int. J. Theor. Appl. Finance 7, 337–384 (2004)
Linetsky, V.: Pricing equity derivatives subject to bankruptcy. Math. Finance 16, 255–282 (2006)
Linetsky, V.: Spectral methods in derivatives pricing. In: Birge, J.R., Linetsky, V. (eds.) Handbook of Financial Engineering, Handbooks in Operations Research and Management Sciences, pp. 213–289. Elsevier, Amsterdam (2008). Chap. 6
McKean, H.: Elementary solutions for certain parabolic partial differential equations. Trans. Am. Math. Soc. 82, 519–548 (1956)
Mendoza-Arriaga, R., Carr, P., Linetsky, V.: Time changed Markov processes in unified credit-equity modeling. Math. Finance 20, 527–569 (2010)
Mendoza-Arriaga, R., Linetsky, V.: Pricing equity default swaps under the jump-to-default extended CEV model. Finance Stoch. 15, 513–540 (2010)
Mendoza-Arriaga, R., Linetsky, V.: Multivariate subordination of Markov processes with financial applications. Math. Finance (2014). doi:10.1111/mafi.12061
Mendoza-Arriaga, R., Linetsky, V.: Time-changed CIR default intensities with two-sided mean-reverting jumps. Ann. Appl. Probab. 24, 811–856 (2014)
Mijatović, A., Pistorius, M.: Continuously monitored barrier options under Markov processes. Math. Finance 23, 1–38 (2013)
Nikiforov, A.F., Uvarov, V.B.: Special Functions of Mathematical Physics: A Unified Introduction with Applications. Birkhäuser, Basel (1988)
Petrella, G., Kou, S.G.: Numerical pricing of discrete barrier and lookback options via Laplace transforms. J. Comput. Finance 8, 1–37 (2004)
Prudnikov, A.P., Brychkov, Y.A., Marichev, O.I.: Integrals and Series, vol. 2, Gordon & Breach, New York (1986)
Reed, M., Simon, B.: Methods of Modern Mathematical Physics, Vol. 1: Functional Analysis. Academic Press, San Diego (1980)
Ricciardi, L.M., Sato, S.: First-passage-time density and moments of the Ornstein–Uhlenbeck process. J. Appl. Probab. 25, 43–57 (1988)
Sato, K.: Lévy Processes and Infinitely Divisible Distributions. Cambridge University Press, Cambridge (1999)
Schilling, R., Song, R., Vondraček, Z.: Bernstein Functions: Theory and Applications. de Gruyter, Berlin (2010)
Vasiček, O.: An equilibrium characterisation of the term structure. J. Financ. Econ. 5, 177–188 (1977)
Wong, E.: The construction of a class of stationary Markov processes. In: Bellman, R. (ed.) Sixteenth Symposium in Applied Mathematics—Stochastic Processes in Mathematical Physics and Engineering, pp. 264–276. Am. Math. Soc., Providence (1964)
Yi, C.: On the first passage time distribution of an Ornstein–Uhlenbeck process. Quant. Finance 10, 957–960 (2010)
Author information
Authors and Affiliations
Corresponding author
Additional information
The research of the first author was supported by the Chinese University of Hong Kong Direct Grant for Research (project code 4055005). The research of the second author was supported by the National Science Foundation under grant DMS-1109506.
Appendix: Supplementary proofs
Appendix: Supplementary proofs
Proof of Proposition 2.4
It is easy to see that we have \(\mathcal{P}^{\phi}_{t}\varphi_{n}(x)=e^{-\phi (\lambda_{n}^{d})t}\varphi_{n}(x)\) after subordination. Hence
We also note that
Define \(p^{\phi}_{t}(x,y):=\int_{(0,\infty)}p^{d}_{s}(x,y)\pi_{t}(ds)\). It is clear that \(p^{\phi}_{t}(x,y)\) is symmetric. We want to show it is also jointly continuous in \(x\) and \(y\).
(i) \(\gamma>0\). Since \(\phi(\lambda)=\gamma\lambda+\int_{(0,\infty )}(1-e^{-\lambda s})\nu(ds)\), we have \(\phi(\lambda_{n}^{d})\geq\lambda _{n}^{d}\gamma t\). Hence
Applying Theorem 2.3 to \((\mathcal{P}^{d}_{t})_{t\geq0}\), the eigenfunctions satisfy
(A.1) implies
From the expression of \(\phi\), \(\sum_{n=1}^{\infty}e^{-\phi(\lambda ^{d}_{n})t}e^{2\gamma t\lambda^{d}_{n}/3}<\infty\). This shows that
So we can apply the dominated convergence theorem and get
(A.2) also implies that \(\sum_{n=1}^{\infty}e^{-\phi (\lambda^{d}_{n})t}\varphi_{n}(x)\varphi_{n}(y)\) converges uniformly on compact sets to \(p^{\phi}_{t}(x,y)\). Since \(\varphi_{n}(x)\) is continuous by applying Theorem 2.3 to \((\mathcal{P}^{d}_{t})_{t\geq0}\), \(p^{\phi}_{t}(x,y)\) is jointly continuous in \(x\) and \(y\).
(ii) The condition in (ii) implies that \(\mathcal{P}^{\phi}_{t}\) is trace-class, and for \(x\) and \(y\) on compacts, \(\sum_{n=1}^{\infty}e^{-\phi(\lambda^{d}_{n})t}|\varphi_{n}(x)\varphi_{n}(y)|<\infty\). This allows us to apply the dominated convergence theorem, and the rest of the proof is similar to the case \(\gamma>0\). □
Proof of Theorem 2.7
Denote by \(\|f\|_{\infty}\) the \(L^{\infty}\)-norm of \(f\). For a given \(t\), from (2.6) in Theorem 2.3, \(|\varphi _{n}(x)|\leq e^{\lambda_{n}t'/2}\sqrt{p_{t'}(x,x)}\) for all \(t'\) such that \(0< t'< t\). Hence
Together with the trace-class condition and our assumption, this implies that the quantity \(\int_{I}|f(y)|\sum_{n=1}^{\infty}e^{-\lambda_{n}t}|\varphi_{n}(x)\varphi_{n}(y)|m(dy)\) is finite. So by the dominated convergence theorem,
We note that
which is finite by (A.3) and the trace-class condition. Hence
i.e. \(\mathcal{P}_{t}f\in L^{2}(I,m)\). Since
is finite and \(p_{t'}(x,x)\) is continuous in \(x\), the convergence in (2.7) is uniform on compact sets and \(\mathcal {P}_{t}f(x)\) is continuous due to the continuity of \(\varphi_{n}(x)\). We also note that \(\|\sum_{n=N+1}^{\infty}f_{n}e^{-\lambda_{n}t}\varphi_{n}(x)\|=\sum _{n=N+1}^{\infty}f^{2}_{n}e^{-2\lambda_{n}t}\rightarrow0 \text{ as } N\rightarrow\infty\). Therefore, the convergence is also in \(L^{2}\). □
Proof of Proposition 4.1
Note that \(\varphi^{2}_{n}(x)=A^{\nu}|\mu+b|x^{2}\frac{n!}{\varGamma(\nu +n+1)}(L_{n}^{\nu}(x))^{2}\). From [1, 22.14.13], we have \(|L^{\nu}_{n}(x)|\leq\frac{\varGamma(\nu+n+1)}{n!\varGamma(\nu+1)}e^{\frac {x}{2}}\) for \(x\geq0\). Hence
Therefore, we have
Since \(\frac{\varGamma(\nu+n+1)}{n!}\sim n^{\nu}\) as \(n\rightarrow\infty \), condition (4.4) implies \(\sum_{n=0}^{\infty}e^{-\lambda_{n}t}\frac{\varGamma(\nu+n+1)}{n!}<\infty\) in the above. To prove the conclusion, note that when \(\mu +b<0\), we obtain that \(\int_{0}^{\infty}xe^{Ax^{-2\beta }/2}m(dx)<\infty\). When \(\mu+b>0\), for any \(S\) bounded above, we have \(\int _{S}xe^{-Ax^{-2\beta}/2}m(dx)<\infty\). □
Proof of Proposition 4.3
The calculation is based on the following identities and then simplification:
They can be obtained from [66, Eq. (1.14.3.3) and Eq. (1.14.3.7)] by a change of variable. In the simplification procedure, we use \(_{2}F_{2}\left(a_{1},a;b_{1},a;x\right )= {}_{1}F_{1}\left(a_{1};b_{1};x\right)\) and \(_{1}F_{1}(-n;\alpha+1;x)=\frac {\varGamma(n+1)}{(\alpha+1)_{n}}L_{n}^{(\alpha)}(x)\). □
Proof of Proposition 4.4
Define \(S(x,h)=\mathbb{P}_{x}[\zeta>h]\). Then
The first term is already calculated in Proposition 4.3. For the second term, we note that when \(\mu+b<0\),
From this we can easily derive the expression for the second term by interchanging integration and summation, using the continuity of inner products. For \(\mu+b>0\),
Again the second term can be calculated by interchanging integration and summation since the above series converges uniformly on \((0,u)\). The only thing left to compute is \(d_{m,n}(0,u)\). By change of variable,
We denote the above integral by \(I\). To calculate \(I\), we use the following identities:
For \(n=0\), by (A.5),
For \(n\geq1\), if \(\frac{c}{|\beta|}=m-n-1\), \(_{1}F_{1}(1-m+\frac{c}{|\beta |};\nu+1;y)=\frac{n!\varGamma(\nu+1)}{\varGamma(\nu+1+n)}L_{n}^{(\nu )}(y)\). Hence
If \(\frac{c}{|\beta|}\neq m-n-1\), by (A.5), integration by parts and (A.6), we get
By (A.7), integration by parts and (A.4), we can also obtain
Multiplying (A.8) by \(\frac{1-m+c/|\beta|}{n(\nu+1)}\), adding it to (A.9) and simplifying gives the expression for \(I\). □
Proof of Proposition 4.5
The interchange of summation and integration in the second equality is justified by the continuity of inner products. The formula for \(g^{i}_{n}((\ell,u)^{c})\) is proved similarly using the identity \(\int _{0}^{\infty}e^{-sy}y^{\alpha}L_{n}^{(\alpha)}(y)dy=\frac{\varGamma(\alpha +n+1)(s-1)^{n}}{n!s^{\alpha+n+1}}\). The calculation of \(\psi_{n}^{(\nu )}(s,x)\) can be done similarly as in [53]. □
Rights and permissions
About this article
Cite this article
Li, L., Linetsky, V. Discretely monitored first passage problems and barrier options: an eigenfunction expansion approach. Finance Stoch 19, 941–977 (2015). https://doi.org/10.1007/s00780-015-0271-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00780-015-0271-1