PRICING COUPON-BOND OPTIONS AND SWAPTIONS IN AFFINE TERM STRUCTURE MODELS

This paper provides a numerically accurate and computationally fast approximation to the prices of European options on coupon-bearing instruments that is applicable to the entire family of affine term structure models. Exploiting the typical shapes of the conditional distributions of the risk factors in affine diffusions, we show that one can reliably compute the relevant probabilities needed for pricing options on coupon-bearing instruments by the same Fourier inversion methods used in the pricing of options on zero-coupon bonds. We apply our theoretical results to the pricing of options on coupon bonds and swaptions, and the calculation of "expected exposures" on swap books. As an empirical illustration, we compute the expected exposures implied by several affine term structure models fit to historical swap yields.     

SWAPTION PRICING IN AFFINE AND OTHER MODELS

This paper shows that Singleton and Umantsev's method for swaption pricing in affine models can be simplified and extended to other models. Two alternative methods for approximating the option exercise boundary are introduced: one based on the multivariate Taylor series expansion, and the other based on duration‐matched zero‐coupon bond approximation. Applied to affine models and quadratic‐Gaussian models, these methods are found to give accurate swaption prices.     

FAST SWAPTION PRICING IN GAUSSIAN TERM STRUCTURE MODELS

We propose a fast and accurate numerical method for pricing European swaptions in multifactor Gaussian term structure models. Our method can be used to accelerate the calibration of such models to the volatility surface. The pricing of an interest rate option in such a model involves evaluating a multidimensional integral of the payoff of the claim on a domain where the payoff is positive. In our method, we approximate the exercise boundary of the state space by a hyperplane tangent to the maximum probability point on the boundary and simplify the multidimensional integration into an analytical form. The maximum probability point can be determined using the gradient descent method. We demonstrate that our method is superior to previous methods by comparing the results to the price obtained by numerical integration.     

PRICING OPTIONS ON VARIANCE IN AFFINE STOCHASTIC VOLATILITY MODELS

We consider the pricing of options written on the quadratic variation of a given stock price process. Using the Laplace transform approach, we determine semi‐explicit formulas in general affine models allowing for jumps, stochastic volatility, and the leverage effect. Moreover, we show that the joint dynamics of the underlying stock and a corresponding variance swap again are of affine form. Finally, we present a numerical example for the Barndorff‐Nielsen and Shephard model with leverage. In particular, we study the effect of approximating the quadratic variation with its predictable compensator.     

SOLUTION OF THE EXTENDED CIR TERM STRUCTURE AND BOND OPTION VALUATION

The extended Cox-Ingersoll-Ross (ECIR) models of interest rates allow for time-dependent parameters in the CIR square-root model. This article presents closed-form pathwise unique solutions of these unsolved stochastic differential equations (s.d.e.s) in terms of functionals of their driving Brownian motion and parameters. It is shown that quadratics in solution of linear s.d.e.s solve the ECIR model if and only if the dimension of the model is a positive integer and that this solution can be achieved by construction of a pathwise unique generalized Ornstein-Uhlenbeck process from the ECIR Brownian motion. For real valued dimensions an extension of the time-change theorem of Dubins and Schwarz (1965) is presented and applied to show that a lognormal process solves the model through a stochastic time change. Pathwise equivalence to a rescaled time-changed Bessel square process is also established. These novel results are applied to characterize zero-hitting time and to produce transition density and zero-hitting conditions for the ECIR spot rate. the CIR term structure is then extended to ECIR under no arbitrage, and its solutions and the transition density are represented under a new ECIR martingale measure. the findings are employed to derive a closed-form ECIR bond option valuation formula which generalizes that obtained by CIR (1985).     

QUADRATIC TERM STRUCTURE MODELS FOR RISK-FREE AND DEFAULTABLE RATES

In this paper, quadratic term structure models (QTSMs) are analyzed and characterized in a general Markovian setting. The primary motivation for this work is to find a useful extension of the traditional QTSM, which is based on an Ornstein–Uhlenbeck (OU) state process, while maintaining the analytical tractability of the model. To accomplish this, the class of quadratic processes, consisting of those Markov state processes that yield QTSM, is introduced. The main result states that OU processes are the only conservative quadratic processes. In general, however, a quadratic potential can be added to allow QTSMs to model default risk. It is further shown that the exponent functions that are inherent in the definition of the quadratic property can be determined by a system of Riccati equations with a unique admissible parameter set. The implications of these results for modeling the term structure of risk-free and defaultable rates are discussed.     

ANAYTICAL SOLUTIONS FOR THE PRICING OF AMERICAN BOND AND YIELD OPTIONS1

In this paper we use the Cox, Ingersoll, and Ross (1985b) single-factor, term structure model and extend it to the pricing of American default-free bond puts. We provide a quasi-analytical formula for these option prices based on recently established mathematical results for Bessel bridges, coupled with the optimal stopping time method. We extend our results to another interest rate contingent claim and provide a quasi-analytical solution for American yield option prices which illustrates the flexibility of our framework.     

EFFICIENT PRICING OF BARRIER OPTIONS AND CREDIT DEFAULT SWAPS IN LÉVY MODELS WITH STOCHASTIC INTEREST RATE

Recently, advantages of conformal deformations of the contours of integration in pricing formulas for European options have been demonstrated in the context of wide classes of Lévy models, the Heston model, and other affine models. Similar deformations were used in one‐factor Lévy models to price options with barrier and lookback features and credit default swaps (CDSs). In the present paper, we generalize this approach to models, where the dynamics of the assets is modeled as , where X is a Lévy process, and the interest rate is stochastic. Assuming that X and r are independent, and , the infinitesimal generator of the pricing semigroup in the model for the short rate, satisfies weak regularity conditions, which hold for popular models of the short rate, we develop a variation of the pricing procedure for Lévy models which is almost as fast as in the case of the constant interest rate. Numerical examples show that about 0.15 second suffices to calculate prices of 8 options of same maturity in a two‐factor model with the error tolerance and less; in a three‐factor model, accuracy of order 0.001–0.005 is achieved in about 0.2 second. Similar results are obtained for quanto CDS, where an additional stochastic factor is the exchange rate. We suggest a class of Lévy models with the stochastic interest rate driven by 1–3 factors, which allows for fast calculations. This class can satisfy the current regulatory requirements for banks mandating sufficiently sophisticated credit risk models.     

COMPLEX LOGARITHMS IN HESTON‐LIKE MODELS

The characteristic functions of many affine jump‐diffusion models, such as Heston's stochastic volatility model and all of its extensions, involve multivalued functions such as the complex logarithm. If we restrict the logarithm to its principal branch, as is done in most software packages, the characteristic function can become discontinuous, leading to completely wrong option prices if options are priced by Fourier inversion. In this paper, we prove without any restrictions that there is a formulation of the characteristic function in which the principal branch is the correct one. Because this formulation is easier to implement and numerically more stable than the so‐called rotation count algorithm of Kahl and Jäckel, we solely focus on its stability in this paper. This paper shows how complex discontinuities can be avoided in the Variance Gamma and Schöbel–Zhu models, as well as in the exact simulation algorithm of the Heston model, recently proposed by Broadie and Kaya.