ANALYTICAL COMPARISONS OF OPTION PRICES IN STOCHASTIC VOLATILITY MODELS

This paper gives an ordering on option prices under various well-known martingale measures in an incomplete stochastic volatility model. Our central result is a comparison theorem that proves convex option prices are decreasing in the market price of volatility risk, the parameter governing the choice of pricing measure. The theorem is applied to order option prices under q -optimal pricing measures. In doing so, we correct orderings demonstrated numerically in Heath, Platen, and Schweizer ( Mathematical Finance , 11(4), 2001) in the special case of the Heston model.     

CONVERGENCE OF BARRIER OPTION PRICES IN THE BINOMIAL MODEL

In this paper, we study the rate of convergence of the European barrier call option price given by the CRR binomial model to the Black–Scholes price as the number of periods n tends to infinity. In general the error is of order and we give explicit formulas for the coefficients of and 1/n in the asymptotic expansion of the error. These coefficients depend on the positions of the barrier and strike in the binomial lattice and enable us to give a rigorous explanation of the observed fact that the error is of order 1/n when n is chosen in an appropriate way.     

ON PROPERTIES OF ANALYTICALLY SOLVABLE FAMILIES OF LOCAL VOLATILITY DIFFUSION MODELS

We present some further developments in the construction and classification of new solvable one‐dimensional diffusion models having transition densities, and other quantities that are fundamental to derivatives pricing, representable in analytically closed form. Our approach is based on so‐called diffusion canonical transformations that produce a large class of multiparameter nonlinear local volatility diffusion models that are mapped onto various simpler diffusions. Using an asymptotic analysis, we arrive at a rigorous boundary classification as well as a characterization with respect to probability conservation and the martingale property of the newly constructed diffusions. Specifically, we analyze and classify in detail four main families of driftless regular diffusion models that arise from the underlying squared Bessel process (the Bessel family), Cox–Ingersoll–Ross process (the confluent hypergeometric family), the Ornstein‐Uhlenbeck diffusion (the OU family), and the Jacobi diffusion (the hypergeometric family). We show that the Bessel family is a superset of the constant elasticity of variance model without drift. The Bessel family, in turn, is nested by the confluent hypergeometric family. For these two families we find further subfamilies of conservative strict supermartingales and nonconservative martingales with an exit boundary. For the new classes of nonconservative regular diffusions we also derive analytically exact first exit time densities that are given in terms of generalized inverse Gaussians and extensions. As for the two other new models, we show that the OU family of processes are conservative strict martingales, whereas the Jacobi family are nonconservative nonmartingales. Considered as asset price diffusion models, we also show that these models demonstrate a wide range of local volatility shapes and option implied volatility surfaces that include various pronounced skew and smile patterns.     

ASYMPTOTIC BEHAVIOR OF DISTRIBUTION DENSITIES IN MODELS WITH STOCHASTIC VOLATILITY. I

We study the asymptotic behavior of distribution densities arising in stock price models with stochastic volatility. The main objects of our interest in the present paper are the density of time averages of a geometric Brownian motion and the density of the stock price process in the Hull–White model. We find explicit formulas for leading terms in asymptotic expansions of these densities and give error estimates. As an application of our results, sharp asymptotic formulas for the price of an Asian option are obtained.     

PRICING ASIAN OPTIONS FOR JUMP DIFFUSION

We construct a sequence of functions that uniformly converge (on compact sets) to the price of an Asian option, which is written on a stock whose dynamics follow a jump diffusion. The convergence is exponentially fast. We show that each element in this sequence is the unique classical solution of a parabolic partial differential equation (not an integro‐differential equation). As a result we obtain a fast numerical approximation scheme whose accuracy versus speed characteristics can be controlled. We analyze the performance of our numerical algorithm on several examples.     

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).     

EXPLICIT SOLUTIONS OF CONSUMPTION-INVESTMENT PROBLEMS IN FINANCIAL MARKETS WITH REGIME SWITCHING

We consider a consumption and investment problem where the market presents different regimes. An investor taking decisions continuously in time selects a consumption–investment policy to maximize his expected total discounted utility of consumption. The market coefficients and the investor's utility of consumption are dependent on the regime of the financial market, which is modeled by an observable finite-state continuous-time Markov chain. We obtain explicit optimal consumption and investment policies for specific HARA utility functions. We show that the optimal policy depends on the regime. We also make an economic analysis of the solutions, and show that for every investor the optimal proportion to allocate in the risky asset is greater in a bull market than in a bear market. This behavior is not affected by the investor's risk preferences. On the other hand, the optimal consumption to wealth ratio depends not only on the regime, but also on the investor's risk tolerance: high risk-averse investors will consume relatively more in a bull market than in a bear market, and the opposite is true for low risk-averse investors.     

NONCONVEXITY OF THE OPTIMAL EXERCISE BOUNDARY FOR AN AMERICAN PUT OPTION ON A DIVIDEND‐PAYING ASSET

We prove that when the dividend rate of the underlying asset following a geometric Brownian motion is slightly larger than the risk‐free interest rate, the optimal exercise boundary of the American put option is not convex.     

TUG‐OF‐WAR,MARKET MANIPULATION,AND OPTION PRICING

We develop an option pricing model based on a tug‐of‐war game. This two‐player zero‐sum stochastic differential game is formulated in the context of a multidimensional financial market. The issuer and the holder try to manipulate asset price processes in order to minimize and maximize the expected discounted reward. We prove that the game has a value and that the value function is the unique viscosity solution to a terminal value problem for a parabolic partial differential equation involving the nonlinear and completely degenerate infinity Laplace operator.     

NO MARGINAL ARBITRAGE OF THE SECOND KIND FOR HIGH PRODUCTION REGIMES IN DISCRETE TIME PRODUCTION–INVESTMENT MODELS WITH PROPORTIONAL TRANSACTION COSTS

We consider a class of production–investment models in discrete time with proportional transaction costs. For linear production functions, we study a natural extension of the no‐arbitrage of the second kind condition introduced by Rásonyi. We show that this condition implies the closedness of the set of attainable claims and is equivalent to the existence of a strictly consistent price system under which the evaluation of future production profits is strictly negative. This allows us to discuss the closedness of the set of terminal wealth in models with nonlinear production, functions which may admit arbitrages of the second kind for low production regimes but not marginally for high production regimes.     

CONVEXITY OF THE EXERCISE BOUNDARY OF THE AMERICAN PUT OPTION ON A ZERO DIVIDEND ASSET

We show that the optimal exercise boundary for the American put option with non-dividend-paying asset is convex. With this convexity result, we then give a simple rigorous argument providing an accurate asymptotic behavior for the exercise boundary near expiry.     

ROBUST UTILITY MAXIMIZATION IN NONDOMINATED MODELS WITH 2BSDE: THE UNCERTAIN VOLATILITY MODEL

The problem of robust utility maximization in an incomplete market with volatility uncertainty is considered, in the sense that the volatility of the market is only assumed to lie between two given bounds. The set of all possible models (probability measures) considered here is nondominated. We propose studying this problem in the framework of second‐order backward stochastic differential equations (2BSDEs for short) with quadratic growth generators. We show for exponential, power, and logarithmic utilities that the value function of the problem can be written as the initial value of a particular 2BSDE and prove existence of an optimal strategy. Finally, several examples which shed more light on the problem and its links with the classical utility maximization one are provided. In particular, we show that in some cases, the upper bound of the volatility interval plays a central role, exactly as in the option pricing problem with uncertain volatility models.     

A MODEL OF OPTIMAL CONSUMPTION UNDER LIQUIDITY RISK WITH RANDOM TRADING TIMES

We consider a portfolio/consumption choice problem in a market model with liquidity risk. The main feature is that the investor can trade and observe stock prices only at exogenous Poisson arrival times. He may also consume continuously from his cash holdings, and his goal is to maximize his expected utility from consumption. This is a mixed discrete/continuous stochastic control problem, non‐standard in the literature. The dynamic programming principle leads to a coupled system of Integro‐Differential Equations (IDE), and we provide a convergent numerical algorithm for the resolution to this coupled system of IDE. Several numerical experiments illustrate the impact of the restricted liquidity trading opportunities, and we measure in particular the utility loss with respect to the classical Merton consumption problem.     

THE EIGENFUNCTION EXPANSION METHOD IN MULTI-FACTOR QUADRATIC TERM STRUCTURE MODELS

We propose the eigenfunction expansion method for pricing options in quadratic term structure models. The eigenvalues, eigenfunctions, and adjoint functions are calculated using elements of the representation theory of Lie algebras not only in the self-adjoint case, but in non-self-adjoint case as well; the eigenfunctions and adjoint functions are expressed in terms of Hermite polynomials. We demonstrate that the method is efficient for pricing caps, floors, and swaptions, if time to maturity is 1 year or more. We also consider subordination of the same class of models, and show that in the framework of the eigenfunction expansion approach, the subordinated models are (almost) as simple as pure Gaussian models. We study the dependence of Black implied volatilities and option prices on the type of non-Gaussian innovations.     

ACHIEVING HIGHER ORDER CONVERGENCE FOR THE PRICES OF EUROPEAN OPTIONS IN BINOMIAL TREES

A new family of binomial trees as approximations to the Black–Scholes model is introduced. For this class of trees, the existence of complete asymptotic expansions for the prices of vanilla European options is demonstrated and the first three terms are explicitly computed. As special cases, a tree with third-order convergence is constructed and the conjecture of Leisen and Reimer that their tree has second-order convergence is proven.     

OPTIMAL SELLING RULES FOR MONETARY INVARIANT CRITERIA: TRACKING THE MAXIMUM OF A PORTFOLIO WITH NEGATIVE DRIFT

Considering a positive portfolio diffusion X with negative drift, we investigate optimal stopping problems of the form where f is a nonincreasing function, τ is the next random time where the portfolio X crosses zero and θ is any stopping time smaller than τ. Hereby, our motivation is the obtention of an optimal selling strategy minimizing the relative distance between the liquidation value of the portfolio and its highest possible value before it reaches zero. This paper unifies optimal selling rules observed for the quadratic absolute distance criteria in this stationary framework with bang–bang type ones observed for monetary invariant criteria but in finite horizon. More precisely, we provide a verification result for the general stopping problem of interest and derive the exact solution for two classical criteria f of the literature. For the power utility criterion with , instantaneous selling is always optimal, which is consistent with previous observations for the Black‐Scholes model in finite observation. On the contrary, for a relative quadratic error criterion, , selling is optimal as soon as the process X crosses a specified function φ of its running maximum . These results reinforce the idea that optimal stopping problems of similar type lead easily to selling rules of very different nature. Nevertheless, our numerical experiments suggest that the practical optimal selling rule for the relative quadratic error criterion is in fact very close to immediate selling.     

EXISTENCE OF AN EQUILIBRIUM WITH DISCONTINUOUS PRICES, ASYMMETRIC INFORMATION, AND NONTRIVIAL INITIAL σ-FIELDS

In this paper, we solve the problems of optimization and equilibrium on a continuous-time financial market with discontinuous prices, in which agents have different random endowments and different information on the structure and future behavior of the prices. Our purpose is to go over and to extend the work of Pikovsky and Karatzas (1996) by using the theory developed by Amendinger (2000) about martingale representation theorems for initially enlarged filtrations, and to generalize the results in the case of discontinuous prices.     

VALUATION OF FLOATING RANGE NOTES IN LÉVY TERM-STRUCTURE MODELS

Turnbull (1995) as well as Navatte and Quittard-Pinon (1999) derived explicit pricing formulae for digital options and range notes in a one-factor Gaussian Heath–Jarrow–Morton (henceforth HJM) model. Nunes (2004) extended their results to a multifactor Gaussian HJM framework. In this paper, we generalize these results by providing explicit pricing solutions for digital options and range notes in the multivariate Lévy term-structure model of Eberlein and Raible (1999) , that is, an HJM-type model driven by a d -dimensional (possibly nonhomogeneous) Lévy process. As a byproduct, we obtain a pricing formula for floating range notes in the special case of a multifactor Gaussian HJM model that is simpler than the one provided by Nunes (2004) .     

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.     

The Market Model of Interest Rate Dynamics

A class of term structure models with volatility of lognormal type is analyzed in the general HJM framework. The corresponding market forward rates do not explode, and are positive and mean reverting. Pricing of caps and floors is consistent with the Black formulas used in the market. Swaptions are priced with closed formulas that reduce (with an extra assumption) to exactly the Black swaption formulas when yield and volatility are flat. A two–factor version of the model is calibrated to the U.K. market price of caps and swaptions and to the historically estimated correlation between the forward rates.