Application of Homotopy Analysis Method to Option Pricing Under Lévy Processes

Option pricing under the Lévy process has been considered an important research direction in the field of financial engineering, where a closed-form expression for the standard European option is available due to the existence of analytically tractable characteristic function according to the Lévy–Khinchin representation. However, this approach cannot be applied to exotic derivatives (such as barrier options) directly, although a large volume of exotic derivatives are actively traded in the current options market. An alternative approach is to solve the corresponding partial integro-differential equation (PIDE) numerically, which is, in fact, time-consuming and is not computationally tractable in general. In this paper, we apply the so-called homotopy analysis method (HAM) to solve the corresponding PIDE in a semi analytic form, being obtained from the following three steps: (1) Apply the Fourier transform to convert the PIDE to an ordinal differential equitation (ODE), and construct a differential system of ODEs. (2) Solve the system of ODEs, where each differential equation is shown to have an analytical solution. (3) Express the option price using the sum of infinite series, where each term may be expressed analytically and derived by applying Steps (1) and (2) recursively. To illustrate our technique more precisely, we take the variance gamma model as an example and provide the semi-analytic form. Numerical examples demonstrate a fast convergence of our proposed method to the prices of European and down-and-out call options with a few number of terms. Note that this method is easy to implement and can be applied to other types of options under general Lévy processes.     

Valuation of asset and volatility derivatives using decoupled time-changed Lévy processes

In this paper we propose a general derivative pricing framework that employs decoupled time-changed (DTC) Lévy processes to model the underlying assets of contingent claims. A DTC Lévy process is a generalized time-changed Lévy process whose continuous and pure jump parts are allowed to follow separate random time scalings; we devise the martingale structure for a DTC Lévy-driven asset and revisit many popular models which fall under this framework. Postulating different time changes for the underlying Lévy decomposition allows the introduction of asset price models consistent with the assumption of a correlated pair of continuous and jump market activity rates; we study one illustrative DTC model of this kind based on the so-called Wishart process. The theory we develop is applied to the problem of pricing not only claims that depend on the price or the volatility of an underlying asset, but also more sophisticated derivatives whose payoffs rely on the joint performance of these two financial variables, such as the target volatility option. We solve the pricing problem through a Fourier-inversion method. Numerical analyses validating our techniques are provided. In particular, we present some evidence that correlating the activity rates could be beneficial for modeling the volatility skew dynamics.     

Time consistency of Lévy models

Time consistency of the models used is an important ingredient to improve risk management. The empirical investigation in this article gives evidence for some models driven by Lévy processes to be highly consistent. This means that they provide a good statistical fit of empirical distributions of returns not only on the timescale used for calibration but on various other timescales as well. As a result these models produce more reliable risk numbers and derivative prices.     

On The Expected Discounted Penalty function for Lévy Risk Processes

Abstract

Dufresne et al. (1991) introduced a general risk model defined as the limit of compound Poisson processes. Such a model is either a compound Poisson process itself or a process with an infinite number of small jumps. Later, in a series of now classical papers, the joint distribution of the time of ruin, the surplus before ruin, and the deficit at ruin was studied (Gerber and Shiu 1997, 1998a, 1998b; Gerber and Landry 1998). These works use the classical and the perturbed risk models and hint that the results can be extended to gamma and inverse Gaussian risk processes.

In this paper we work out this extension to a generalized risk model driven by a nondecreasing Lévy process. Unlike the classical case that models the individual claim size distribution and obtains from it the aggregate claims distribution, here the aggregate claims distribution is known in closed form. It is simply the one-dimensional distribution of a subordinator. Embedded in this wide family of risk models we find the gamma, inverse Gaussian, and generalized inverse Gaussian processes. Expressions for the Gerber-Shiu function are given in some of these special cases, and numerical illustrations are provided.     

Existence of Lévy term structure models

Lévy driven term structure models have become an important subject in the mathematical finance literature. This paper provides a comprehensive analysis of the Lévy driven Heath–Jarrow–Morton type term structure equation. This includes a full proof of existence and uniqueness in particular, which seems to have been lacking in the finance literature so far.      

Spectral calibration of exponential Lévy models

We investigate the problem of calibrating an exponential Lévy model based on market prices of vanilla options. We show that this inverse problem is in general severely ill-posed and we derive exact minimax rates of convergence. The estimation procedure we propose is based on the explicit inversion of the option price formula in the spectral domain and a cut-off scheme for high frequencies as regularisation.     

Implied Lévy volatility

We study the cause of large fluctuations in prices on the London Stock Exchange. This is done at the microscopic level of individual events, where an event is the placement or cancellation of an order to buy or sell. We show that price fluctuations caused by individual market orders are essentially independent of the volume of orders. Instead, large price fluctuations are driven by liquidity fluctuations, variations in the market's ability to absorb new orders. Even for the most liquid stocks there can be substantial gaps in the order book, corresponding to a block of adjacent price levels containing no quotes. When such a gap exists next to the best price, a new order can remove the best quote, triggering a large midpoint price change. Thus, the distribution of large price changes merely reflects the distribution of gaps in the limit order book. This is a finite size effect, caused by the granularity of order flow: in a market where participants place many small orders uniformly across prices, such large price fluctuations would not happen. We show that this also explains price fluctuations on longer timescales. In addition, we present results suggesting that the risk profile varies from stock to stock, and is not universal: lightly traded stocks tend to have more extreme risks.     

Evaluation of the MEMM, Parameter Estimation and Option Pricing for Geometric Lévy Processes

In this paper, we shall propose a useful approach to evaluate concretely the MEMM (minimal entropy martingale measure) for the typical geometric Lévy processes such as compound Poisson, stable, VG (Variance Gamma), CGMY (Carr-Geman-Madan-Yor), NIG (Normal Inverse Gaussian), etc. In addition, we shall estimate the parameters of geometric Lévy processes and value the European call option and Asian call option using the Nikkei financial data.     

Pricing Surrender Risk in Ratchet Equity-Index Annuities under Regime-Switching Lévy Processes

This article presents a numerical method of pricing the surrender risk in Ratchet equity-index annuities (EIAs). We assume that log-returns of the underlying fund belong to a class of regime-switching models where the parameters are allowed to change randomly according to a hidden Markov chain. The defining feature of these models is the fact that in each regime the characteristic function of log-returns is assumed to have an analytical form. The presented method provides an unified pricing framework within this class and includes the recently developed COS method as a particular case. This aspect of the method is particularly useful when pricing Ratchet options embedded in EIAs, for which the COS method exhibits a low rate of convergence. Our numerical results confirm that for models considered in this article the proposed approach improves convergence of the COS method without increasing the computational burden.     

Tangent Lévy market models

In this paper, we introduce a new class of models for the time evolution of the prices of call options of all strikes and maturities. We capture the information contained in the option prices in the density of some time-inhomogeneous Lévy measure (an alternative to the implied volatility surface), and we set this static code-book in motion by means of stochastic dynamics of It?’s type in a function space, creating what we call a tangent Lévy model. We then provide the consistency conditions, namely, we show that the call prices produced by a given dynamic code-book (dynamic Lévy density) coincide with the conditional expectations of the respective payoffs if and only if certain restrictions on the dynamics of the code-book are satisfied (including a drift condition à la HJM). We then provide an existence result, which allows us to construct a large class of tangent Lévy models, and describe a specific example for the sake of illustration.     

A cross-currency Lévy market model

The Lévy Libor or market model which was introduced in Eberlein and Özkan (The Lévy Libor model. Financ. Stochast., 2005, 9, 327–348) is extended to a multi-currency setting. As an application we derive closed form pricing formulas for cross-currency derivatives. Foreign caps and floors and cross-currency swaps are studied in detail. Numerically efficient pricing algorithms based on bilateral Laplace transforms are derived. A calibration example is given for a two-currency setting (EUR, USD).     

Numerical methods for Lévy processes

We survey the use and limitations of some numerical methods for pricing derivative contracts in multidimensional geometric Lévy models.      

Skewness premium with Lévy processes

We study the skewness premium (SK) introduced by Bates [J. Finance, 1991, 46(3), 1009–1044] in a general context using Lévy processes. Under a symmetry condition, Fajardo and Mordecki [Quant. Finance, 2006, 6(3), 219–227] obtained that SK is given by Bates' x% rule. In this paper, we study SK in the absence of that symmetry condition. More exactly, we derive sufficient conditions for the excess of SK to be positive or negative, in terms of the characteristic triplet of the Lévy process under a risk-neutral measure.     

Confidence sets in nonparametric calibration of exponential Lévy models

Confidence intervals and joint confidence sets are constructed for the nonparametric calibration of exponential Lévy models based on prices of European options. To this end, we show joint asymptotic normality in the spectral calibration method for the estimators of the volatility, the drift, the jump intensity and the Lévy density at finitely many points.     

Asian options and meromorphic Lévy processes

One method to compute the price of an arithmetic Asian option in a Lévy driven model is based on an exponential functional of the underlying Lévy process: If we know the distribution of the exponential functional, we can calculate the price of the Asian option via the inverse Laplace transform. In this paper, we consider pricing Asian options in a model driven by a general meromorphic Lévy process. We prove that the exponential functional is equal in distribution to an infinite product of independent beta random variables, and its Mellin transform can be expressed as an infinite product of gamma functions. We show that these results lead to an efficient algorithm for computing the price of the Asian option via the inverse Mellin–Laplace transform, and we compare this method with some other techniques.