Small-noise limit of the quasi-Gaussian log-normal HJM model

Quasi-Gaussian HJM models are a popular approach for modeling the dynamics of the yield curve. This is due to their low dimensional Markovian representation, which greatly simplifies their numerical implementation. We present a qualitative study of the solutions of the quasi-Gaussian log-normal HJM model. Using a small-noise deterministic limit we show that the short rate may explode to infinity in finite time. This implies the explosion of the Eurodollar futures prices in this model. We derive explicit explosion criteria under mild assumptions on the shape of the yield curve.     

Lognormality of rates and term structure models

A term structure model with lognormal type volatility structure is proposed. The Heath, Jarrow and Morton (HJM) framework, coupled with the theory of stochastic evolution equations in infinite dimensions, is used to show that the resulting instantaneous rates are well defined (they do not explode) and remain positive, contrary to those derived in [2]. They are also bounded from below and above by lognormal processes. The model can be used to price and hedge caps, swaptions and other interest rate and currency derivatives including the Eurodollar futures contract, which requires integrability of one over zero coupon bond. This extends results obtained by Sandmann and Sondermann in [22] and [23] for Markovian lognormal short rates to (non-Markovian) lognormal forward rates. We show also existence of invariant measures for the proposed term structure dynamics     

Arbitrary Initial Term Structure within the CIR Model: A Perturbative Solution

Single‐factor interest rate models with constant coefficients are not consistent with arbitrary initial term structures. An extension which allows both arbitrary initial term structure and analytical tractability has been provided only in the Gaussian case. In this paper, within the context of the HJM methodology, an extension of the CIR model is provided which admits arbitrary initial term structure. It is shown how to calculate bond prices via a perturbative approach, and closed formulas are provided at every order. Since the parameter selected for the expansion is typically estimated to be small, the perturbative approach turns out to be adequate to our purpose. Using results on affine models, the extended CIR model is estimated via maximum likelihood on a time series of daily interest rate yields. Results show that the CIR model has to be rejected with respect to the proposed extension, and it is pointed out that the extended CIR model provides a more flexible characterization of the link between risk neutral and natural probability.     

Calibrating the Black-Derman-Toy model: some theoretical results

The Black-Derman-Toy (BDT) model is a popular one-factor interest rate model that is widely used by practitioners. One of its advantages is that the model can be calibrated to both the current market term structure of interest rate and the current term structure of volatilities. The input term structure of volatility can be either the short term volatility or the yield volatility. Sandmann and Sondermann derived conditions for the calibration to be feasible when the conditional short rate volatility is used. In this paper conditions are investigated under which calibration to the yield volatility is feasible. Mathematical conditions for this to happen are derived. The restrictions in this case are more complicated than when the short rate volatilities are used since the calibration at each time step now involves the solution of two non-linear equations. The theoretical results are illustrated by showing numerically that in certain situations the calibration based on the yield volatility breaks down for apparently plausible inputs. In implementing the calibration from period n to period n + 1, the corresponding yield volatility has to lie within certain bounds. Under certain circumstances these bounds become very tight. For yield volatilities that violate these bounds, the computed short rates for the period (n, n + 1) either become negative or else explode and this feature corresponds to the economic intuition behind the breakdown.     

Rare Shock,Two-Factor Stochastic Volatility and Currency Option Pricing

Abstract

In this paper, we develop an option valuation model where the dynamics of the spot foreign exchange rate is governed by a two-factor Markov-modulated jump-diffusion process. The short-term fluctuation of stochastic volatility is driven by a Cox–Ingersoll–Ross (CIR) process and the long-term variation of stochastic volatility is driven by a continuous-time Markov chain which can be interpreted as economy states. Rare events are governed by a compound Poisson process with log-normal jump amplitude and stochastic jump intensity is modulated by a common continuous-time Markov chain. Since the market is incomplete under regime-switching assumptions, we determine a risk-neutral martingale measure via the Esscher transform and then give a pricing formula of currency options. Numerical results are presented for investigating the impact of the long-term volatility and the annual jump intensity on option prices.     

Energy futures prices: term structure models with Kalman filter estimation

We present a class of multi-factor stochastic models for energy futures prices, similar to the interest rate futures models recently formulated by Heath. We do not postulate directly the risk-neutral processes followed by futures prices, but define energy futures prices in terms of a spot price, not directly observable, driven by several stochastic factors. Our formulation leads to an expression for futures prices which is well suited to the application of Kalman filtering techniques together with maximum likelihood estimation methods. Based on these techniques, we perform an empirical study of a one- and a two-factor model for futures prices for natural gas.     

GENERALIZED STOCHASTIC DURATION INMARKOVIAN HEATH-JARROW-MORTONFRAMEWORK

1 IntroductionIll moderll financial industry, risk managen1ent is a major task tliat fiuancial institutionsnlust deal witli in every tradiug day alld it is becouilng urore and more important for mailltaining tl1eir access to clieap capital and meeting risk-based caPital requirements. Meanwhile,sonle fiuaucial iustitutions (FI) sucl1 as co1nuercial ba1lks, iusurallce companies and securitiescolllpanies. etc., hold a large proportioll of fixed income security which price is sensitive totlle mar…     

An extended Heath–Jarrow–Morton risk-neutral drift

Using a finite dimensional Hilbert space framework, this work proposes a new derivation of the HJM [D. Heath, R. Jarrow, A. Morton, Bond pricing and the term structure of interest rates: A new methodology for contingent claims valuation, Econometrica 60 (1992) 77–105] risk-neutral drift that takes into account nonzero instantaneous correlations between factors. The results obtained generalize the original HJM risk-neutral drift and provide an approach by which interest rate derivatives can be priced using functions of directly observable factors.     

Pricing Fixed-Income Securities in an Information-Based Framework

Abstract

The purpose of this article is to introduce a class of information-based models for the pricing of fixed-income securities. We consider a set of continuous-time processes that describe the flow of information concerning market factors in a monetary economy. The nominal pricing kernel is assumed to be given at any specified time by a function of the values of information processes at that time. Using a change-of-measure technique, we derive explicit expressions for the prices of nominal discount bonds and deduce the associated dynamics of the short rate of interest and the market price of risk. The interest rate positivity condition is expressed as a differential inequality. An example that shows how the model can be calibrated to an arbitrary initial yield curve is presented. We proceed to model the price level, which is also taken at any specified time to be given by a function of the values of the information processes at that time. A simple model for a stochastic monetary economy is introduced in which the prices of the nominal discount bonds and inflation-linked notes can be expressed in terms of aggregate consumption and the liquidity benefit generated by the money supply.     

Convexity Bias in Eurodollar Futures Prices: A Dimension-Free HJM Criterion

In the theory of interest rate futures, the difference between the futures rate and forward rate is called the “convexity bias,” and there are several widely offered reasons why the convexity bias should be positive. Nevertheless, it is not infrequent that the empirical the bias is observed to be negative. Moreover, in its most general form, the benchmark Heath–Jarrow–Morton (HJM) term structure model is agnostic on the question of the sign of the bias; it allows for models where the convexity bias can be positive or negative. In partial support of the practitioner’s arguments, we develop a simple scalar condition within the HJM framework that suffices to guarantee that the convexity bias is positive. Moreover, when we check this condition on the LIBOR futures data, we find strong empirical support for the new condition. The empirical validity of the sufficient condition and the periodic observation of negative bias, therefore leads one to a paradoxical situation where either (1) there are arbitrage possibilities or (2) a large subclass of HJM models provide interest rate dynamics that fail to capture a fundamental feature of LIBOR futures.     

Three Ways to Solve for Bond Prices in the Vasicek Model

Three approaches in obtaining the closed-form solution of the Vasicek bond pricing problem are discussed in this exposition. A derivation based solely on the distribution of the short rate process is reviewed. Solving the bond price partial differential equation (PDE) is another method. In this paper, this PDE is derived via a martingale approach and the bond price is determined by integrating ordinary differential equations. The bond pricing problem is further considered within the Heath-Jarrow-Morton (HJM) framework in which the analytic solution follows directly from the short rate dynamics under the forward measure.     

A multiple-curve HJM model of interbank risk

In the aftermath of the 2007?C2009 financial crisis, a variety of spreads have developed between quantities that had been essentially the same until then, notably LIBOR?COIS spreads, LIBOR?COIS swap spreads, and basis swap spreads. By the end of 2011, with the sovereign credit crisis, these spreads were again significant. In this paper we study the valuation of LIBOR interest rate derivatives in a multiple-curve setup, which accounts for the spreads between a risk-free discount curve and LIBOR curves. Towards this end we resort to a defaultable HJM methodology, in which these spreads are explained by an implied default intensity of the LIBOR contributing banks, possibly in conjunction with an additional liquidity factor. Markovian short rate specifications are given in the form of an extended CIR and a Lévy Hull?CWhite model for a risk-free short rate and a LIBOR short spread. The use of Lévy drivers leads to the more parsimonious specification. Numerical values of the FRA spreads and the basis swap spreads computed with the latter largely cover the ranges of values observed even at the peak of the 2007?C2009 crisis.     

The American Foreign Exchange Option in Time-Dependent One-Dimensional Diffusion Model for Exchange Rate

The classical Garman-Kohlhagen model for the currency exchange assumes that the domestic and foreign currency risk-free interest rates are constant and the exchange rate follows a log-normal diffusion process. In this paper we consider the general case, when exchange rate evolves according to arbitrary one-dimensional diffusion process with local volatility that is the function of time and the current exchange rate and where the domestic and foreign currency risk-free interest rates may be arbitrary continuous functions of time. First non-trivial problem we encounter in time-dependent case is the continuity in time argument of the value function of the American put option and the regularity properties of the optimal exercise boundary. We establish these properties based on systematic use of the monotonicity in volatility for the value functions of the American as well as European options with convex payoffs together with the Dynamic Programming Principle and we obtain certain type of comparison result for the value functions and corresponding exercise boundaries for the American puts with different strikes, maturities and volatilities. Starting from the latter fact that the optimal exercise boundary curve is left continuous with right-hand limits we give a mathematically rigorous and transparent derivation of the significant early exercise premium representation for the value function of the American foreign exchange put option as the sum of the European put option value function and the early exercise premium. The proof essentially relies on the particular property of the stochastic integral with respect to arbitrary continuous semimartingale over the predictable subsets of its zeros. We derive from the latter the nonlinear integral equation for the optimal exercise boundary which can be studied by numerical methods.     

Limiting Connection between Discrete and Continuous Time Forward Interest Rate Curve Models

We study discrete time Heath–Jarrow–Morton (HJM) type of interest rate curve models, where the forward interest rates – in contrast to the classical HJM models – are driven by a random field. Our main aim is to investigate the relationship between the discrete time forward interest rate curve model and its continuous time counterpart. We derive a general result on the convergence of discrete time models and we give special focus on the nearly unit root spatial autoregression model.     

Monte Carlo Euler approximations of HJM term structure financial models

We present Monte Carlo-Euler methods for a weak approximation problem related to the Heath-Jarrow-Morton (HJM) term structure model, based on Itô stochastic differential equations in infinite dimensional spaces, and prove strong and weak error convergence estimates. The weak error estimates are based on stochastic flows and discrete dual backward problems, and they can be used to identify different error contributions arising from time and maturity discretization as well as the classical statistical error due to finite sampling. Explicit formulas for efficient computation of sharp error approximation are included. Due to the structure of the HJM models considered here, the computational effort devoted to the error estimates is low compared to the work to compute Monte Carlo solutions to the HJM model. Numerical examples with known exact solution are included in order to show the behavior of the estimates.     

Pricing caps with HJM models: The benefits of humped volatility

In this paper we compare different multifactor HJM models with humped volatility structures, to each other and to models with strictly decreasing volatility. All the models are estimated on Euribor and swap rates panel data maximizing the quasi-likelihood function obtained from the Kalman filter. We develop the analysis in two steps: first we study the in-sample properties of the estimated models, then we test the pricing performance on caps. We find the humped volatility specification to greatly improve the model estimation and to provide sufficiently accurate cap prices, although the models has been calibrated on interest rates data and not on cap prices. Moreover, we find the two-factor humped volatility model to outperform the three-factor models in pricing caps.     

A Wiener Chaos Approach to Hyperbolic SPDEs

We construct generalized weighted Wiener chaos solutions for hyperbolic linear SPDEs driven by a cylindrical Brownian motion. Explicit conditions for the existence, uniqueness, and regularity of generalized (Wiener Chaos) solutions are established in Sobolev spaces. An equivalence relation between the Wiener Chaos solution and the traditional one is established. The Heath–Jarrow–Morton (HJM) forward rate model is used as an example to illustrate the general construction.     

A rheological model for arbitrary symmetric distortion of the yield surface

A phenomenological model of metal viscoplasticity, which takes combined isotropic, kinematic, and distortional hardening into account, is motivated by a new rheological model. The distinctive advantage of the material model is that any smooth convex saturated form of the yield surface which is symmetric with respect to the recent loading direction can be captured. In particular, an arbitrary sharpening of the saturated yield locus in the loading direction combined with a flattening on the opposite side can be covered. Moreover, the yield locus evolves smoothly and its convexity is guaranteed at each hardening stage. The underlying two-dimensional rheological analogy can be used to provide insight into the main constitutive assumptions. This rheological model is utilized as a guideline for the construction of phenomenological constitutive relations. The distortion of the yield surface is described with the help of a so-called distortional backstress. Thus, 2nd rank tensors are utilized only. The resulting material model is thermodynamically consistent. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Valuation of Two-Factor Interest Rate Contingent Claims Using Green's Theorem

Abstract

Over the years a number of two-factor interest rate models have been proposed that have formed the basis for the valuation of interest rate contingent claims. This valuation equation often takes the form of a partial differential equation that is solved using the finite difference approach. In the case of two-factor models this has resulted in solving two second-order partial derivatives leading to boundary errors, as well as numerous first-order derivatives. In this article we demonstrate that using Green's theorem, second-order derivatives can be reduced to first-order derivatives that can be easily discretized; consequently, two-factor partial differential equations are easier to discretize than one-factor partial differential equations. We illustrate our approach by applying it to value contingent claims based on the two-factor CIR model. We provide numerical examples that illustrate that our approach shows excellent agreement with analytical prices and the popular Crank–Nicolson method.     

Two-boundary lattice paths and parking functions

We describe an involution on a set of sequences associated with lattice paths with north or east steps constrained to lie between two arbitrary boundaries. This involution yields recursions (from which determinantal formulas can be derived) for the number and area enumerator of such paths. An analogous involution can be defined for parking functions with arbitrary lower and upper bounds. From this involution, we obtained determinantal formulas for the number and sum enumerator of such parking functions. For parking functions, there is an alternate combinatorial inclusion–exclusion approach. The recursions also yield Appell relations. In certain special cases, these Appell relations can be converted into rational or algebraic generating functions.