A simple approach for pricing barrier options with time-dependent parameters

Abstract

In this paper we present a simple and easy-to-use method for computing accurate estimates (in closed form) of Black-Scholes barrier option prices with time-dependent parameters. This new approach is also able to provide tight upper and lower bounds (in closed form) for the exact barrier option prices.     

A general closed-form spread option pricing formula

We propose a new accurate method for pricing European spread options by extending the lower bound approximation of Bjerksund and Stensland (2011) beyond the classical Black–Scholes framework. This is possible via a procedure requiring a univariate Fourier inversion. In addition, we are also able to obtain a new tight upper bound. Our method provides also an exact closed form solution via Fourier inversion of the exchange option price, generalizing the Margrabe (1978) formula. The method is applicable to models in which the joint characteristic function of the underlying assets forming the spread is known analytically. We test the performance of these new pricing algorithms performing numerical experiments on different stochastic dynamic models.     

Explaining the volatility smile: non-parametric versus parametric option models

We employ a “non-parametric” pricing approach of European options to explain the volatility smile. In contrast to “parametric” models that assume that the underlying state variable(s) follows a stochastic process that adheres to a strict functional form, “non-parametric” models directly fit the end distribution of the underlying state variable(s) with statistical distributions that are not represented by parametric functions. We derive an approximation formula which prices S&P 500 index options in closed form which corresponds to the lower bound recently proposed by Lin et al. (Rev Quant Financ Account 38(1):109–129, 2012). Our model yields option prices that are more consistent with the data than the option prices that are generated by several widely used models. Although a quantitative comparison with other non-parametric models is more difficult, there are indications that our model is also more consistent with the data than these models.     

Orthogonal expansions for VIX options under affine jump diffusions

In this work we derive new closed-form pricing formulas for VIX options in the jump-diffusion SVJJ model proposed by Duffie et al. [Econometrica, 2000, 68, 1343–1376]. Our approach is based on the classic methodology of approximating a density function with an orthogonal expansion of polynomials weighted by a kernel. Orthogonal expansions based on the Gaussian distribution, such as Edgeworth or Gram–Charlier expansions, have been successfully employed by a number of authors in the context of equity options. However, these expansions are not quite suitable for volatility or variance densities as they inherently assign positive mass to the negative real line. Here we approximate option prices via expansions that instead are based on kernels defined on the positive real line. Specifically, we consider a flexible family of distributions, which generalizes the gamma kernel associated with the classic Laguerre expansions. The method can be employed whenever the moments of the underlying variance distribution are known. It provides fast and accurate price computations, and therefore it represents a valid and possibly more robust alternative to pricing techniques based on Fourier transform inversions.     

A numerical method to price exotic path-dependent options on an underlying described by the Heston stochastic volatility model

We consider the problem of pricing European exotic path-dependent derivatives on an underlying described by the Heston stochastic volatility model. Lipton has found a closed form integral representation of the joint transition probability density function of underlying price and variance in the Heston model. We give a convenient numerical approximation of this formula and we use the obtained approximated transition probability density function to price discrete path-dependent options as discounted expectations. The expected value of the payoff is calculated evaluating an integral with the Monte Carlo method using a variance reduction technique based on a suitable approximation of the transition probability density function of the Heston model. As a test case, we evaluate the price of a discrete arithmetic average Asian option, when the average over n = 12 prices is considered, that is when the integral to evaluate is a 2n = 24 dimensional integral. We show that the method proposed is computationally efficient and gives accurate results.     

The pricing of basket-spread options

Since the pioneering paper of Black and Scholes was published in 1973, enormous research effort has been spent on finding a multi-asset variant of their closed-form option pricing formula. In this paper, we generalize the Kirk [Managing Energy Price Risk, 1995] approximate formula for pricing a two-asset spread option to the case of a multi-asset basket-spread option. All the advantageous properties of being simple, accurate and efficient are preserved. As the final formula retains the same functional form as the Black–Scholes formula, all the basket-spread option Greeks are also derived in closed form. Numerical examples demonstrate that the pricing and hedging errors are in general less than 1% relative to the benchmark results obtained by numerical integration or Monte Carlo simulation with 10 million paths. An implicit correction method is further applied to reduce the pricing errors by factors of up to 100. The correction is governed by an unknown parameter, whose optimal value is found by solving a non-linear equation. Owing to its simplicity, the computing time for simultaneous pricing and hedging of basket-spread option with 10 underlying assets or less is kept below 1 ms. When compared against the existing approximation methods, the proposed basket-spread option formula coupled with the implicit correction turns out to be one of the most robust and accurate methods.     

Options on realized variance by transform methods: a non-affine stochastic volatility model

In this paper we study the pricing and hedging of options on realized variance in the 3/2 non-affine stochastic volatility model by developing efficient transform-based pricing methods. This non-affine model gives prices of options on realized variance that allow upward-sloping implied volatility of variance smiles. Heston's model [Rev. Financial Stud., 1993, 6, 327–343], the benchmark affine stochastic volatility model, leads to downward-sloping volatility of variance smiles—in disagreement with variance markets in practice. Using control variates, we propose a robust method to express the Laplace transform of the variance call function in terms of the Laplace transform of the realized variance. The proposed method works in any model where the Laplace transform of realized variance is available in closed form. Additionally, we apply a new numerical Laplace inversion algorithm that gives fast and accurate prices for options on realized variance, simultaneously at a sequence of variance strikes. The method is also used to derive hedge ratios for options on variance with respect to variance swaps.     

Analysis of the Nonlinear Option Pricing Model Under Variable Transaction Costs

In this paper we analyze a nonlinear Black–Scholes model for option pricing under variable transaction costs. The diffusion coefficient of the nonlinear parabolic equation for the price V is assumed to be a function of the underlying asset price and the Gamma of the option. We show that the generalizations of the classical Black–Scholes model can be analyzed by means of transformation of the fully nonlinear parabolic equation into a quasilinear parabolic equation for the second derivative of the option price. We show existence of a classical smooth solution and prove useful bounds on the option prices. Furthermore, we construct an effective numerical scheme for approximation of the solution. The solutions are obtained by means of the efficient numerical discretization scheme of the Gamma equation. Several computational examples are presented.     

American option pricing under the double Heston model based on asymptotic expansion

This paper focuses on pricing American put options under the double Heston model proposed by Christoffersen et al. By introducing an explicit exercise rule, we obtain the asymptotic expansion of the solution to the partial differential equation for pricing American put options. We calculate American option price by the sum of the European option price and the early exercise premium. The early exercise premium is calculated by the difference between the American and European option prices based on asymptotic expansions. The European option price is obtained by the efficient COS method. Based on the obtained American option price, the double Heston model is calibrated by minimizing the distance between model and market prices, which yields an optimization problem that is solved by a differential evolution algorithm combined with the Matlab function fmincon.m. Numerical results show that the pricing approach is fast and accurate. Empirical results show that the double Heston model has better performance in pricing short-maturity American put options and capturing the volatility term structure of American put options than the Heston model.     

The multivariate Variance Gamma model: basket option pricing and calibration

In this paper, we propose a methodology for pricing basket options in the multivariate Variance Gamma model introduced in Luciano and Schoutens [Quant. Finance 6(5), 385–402]. The stock prices composing the basket are modelled by time-changed geometric Brownian motions with a common Gamma subordinator. Using the additivity property of comonotonic stop-loss premiums together with Gauss-Laguerre polynomials, we express the basket option price as a linear combination of Black & Scholes prices. Furthermore, our new basket option pricing formula enables us to calibrate the multivariate VG model in a fast way. As an illustration, we show that even in the constrained situation where the pairwise correlations between the Brownian motions are assumed to be equal, the multivariate VG model can closely match the observed Dow Jones index options.     

A new well-posed algorithm to recover implied local volatility

Abstract

This paper presents a new algorithm to calibrate the option pricing model, i.e. the algorithm that recovers the implied local volatility function from market option prices in the optimal control framework. A unique optimal control is shown to exist. Our algorithm is well-posed. Our numerical experiments show that, with the help of the techniques developed in the field of optimal control, the local volatility function is recovered very well.     

Target volatility option pricing in the lognormal fractional SABR model

We examine in this article the pricing of target volatility options in the lognormal fractional SABR model. A decomposition formula of Itô's calculus yields an approximation formula for the price of a target volatility option in small time by the technique of freezing the coefficient. A decomposition formula in terms of Malliavin derivatives is also provided. Alternatively, we also derive closed form expressions for a small volatility of volatility expansion of the price of a target volatility option. Numerical experiments show the accuracy of the approximations over a reasonably wide range of parameters.     

Pricing currency options in the presence of time-varying volatility and non-normalities

A new framework is developed for pricing currency options in the case where the distribution of exchange rate returns exhibits time-varying volatility and non-normalities. A forward-looking volatility structure is adopted whereby volatility is expressed as a function of currency returns over the life of the contract. Time to maturity and moneyness effects in volatility are also modelled. An analytical solution for the option price is obtained up to a one-dimensional integral in the real plane, enabling option prices to be computed efficiently and accurately. The proposed modelling framework is applied to European currency call options for the UK pound written on the US dollar, over the period October 1997–June 1998. The results show that pricing higher order moments improves both within-sample fit and out-of-sample prediction of observed option prices, as well as having important implications for constructing hedged portfolios and managing risk.     

A chaos expansion approach under hybrid volatility models

In this paper, we propose an approximation method based on the Wiener–Ito chaos expansion for the pricing of European contingent claims. Our method is applicable to widely used option pricing models such as local volatility models, stochastic volatility models and their combinations. This method is useful in practice since the resulting approximation formula is not computationally expensive, hence it is suitable for calibration purposes. We will show through some numerical examples that our approximation remains quite good even for the long maturity and/or the high volatility cases, which is a desired feature. As an example, we propose a hybrid volatility model and apply our approximation formula to the JPY/USD currency option market obtaining very accurate results.     

GARCH vs. stochastic volatility: Option pricing and risk management

In this paper we compare the out-of-sample performance of two common extensions of the Black–Scholes option pricing model, namely GARCH and stochastic volatility (SV). We calibrate the three models to intraday FTSE 100 option prices and apply two sets of performance criteria, namely out-of-sample valuation errors and Value-at-Risk (VaR) oriented measures. When we analyze the fit to observed prices, GARCH clearly dominates both SV and the benchmark Black–Scholes model. However, the predictions of the market risk from hypothetical derivative positions show sizable errors. The fit to the realized profits and losses is poor and there are no notable differences between the models. Overall, we therefore observe that the more complex option pricing models can improve on the Black–Scholes methodology only for the purpose of pricing, but not for the VaR forecasts.     

Covariance dependent kernels,a Q-affine GARCH for multi-asset option pricing

This paper introduces a class of multivariate GARCH models that extends the existing literature by explicitly modeling correlation dependent pricing kernels. A large subclass admits closed-form recursive solutions for the moment generating function under the risk-neutral measure, which permits efficient pricing of multi-asset options. We perform a full calibration to three bivariate series of index returns and their corresponding volatility indexes in a joint maximum likelihood estimation. The results empirically confirm the presence of correlation dependance in addition to the well known variance dependance in the pricing kernel. The model improves both the overall likelihood and the VIX-implied likelihoods, with a better fitting of marginal distributions, e.g., 15% less error on one-asset option prices. The new degree of freedom is also shown to significantly impact the shape of marginal and joint pricing kernels, and leads to up to 53% differences for out-of-the-money two-asset correlation option prices.     

A Comparison of Option Prices Under Different Pricing Measures in a Stochastic Volatility Model with Correlation

This paper investigates option prices in an incomplete stochastic volatility model with correlation. In a general setting, we prove an ordering result which says that prices for European options with convex payoffs are decreasing in the market price of volatility risk.As an example, and as our main motivation, we investigate option pricing under the class of q-optimal pricing measures. The q-optimal pricing measure is related to the marginal utility indifference price of an agent with constant relative risk aversion. Using the ordering result, we prove comparison theorems between option prices under the minimal martingale, minimal entropy and variance-optimal pricing measures. If the Sharpe ratio is deterministic, the comparison collapses to the well known result that option prices computed under these three pricing measures are the same.As a concrete example, we specialize to a variant of the Hull-White or Heston model for which the Sharpe ratio is increasing in volatility. For this example we are able to deduce option prices are decreasing in the parameter q. Numerical solution of the pricing pde corroborates the theory and shows the magnitude of the differences in option price due to varying q.JEL Classification: D52, G13     

Computing the Gerber–Shiu function by frame duality projection

Inspired by some works of Kirkby, J. L. [2015. Efficient option pricing by frame duality with the fast Fourier transform. SIAM Journal on Financial Mathematics 6(1), 713–747; 2016. An efficient transform method for Asian option pricing. SIAM Journal on Financial Mathematics 7(1), 845–892], we present a systematic study on effectively computing the Gerber–Shiu function in the Lévy risk model, where the frame duality projection is used for approximation. By introducing an auxiliary function, we provide a smooth extension of the Gerber–Shiu function, which has closed-form Fourier transform and is differentiable over the whole real line under some conditions. The objective function is approximated by its frame duality projection onto a Riesz basis, and the projection coefficients are readily computed by the fast Fourier transform algorithm. Error analysis is made and the effectiveness of our results will be further illustrated in the numerical experiments.     

Pricing bounds and approximations for discrete arithmetic Asian options under time-changed Lévy processes

We derive efficient and accurate analytical pricing bounds and approximations for discrete arithmetic Asian options under time-changed Lévy processes. By extending the conditioning variable approach, we derive the lower bound on the Asian option price and construct an upper bound based on the sharp lower bound. We also consider the general partially exact and bounded (PEB) approximations, which include the sharp lower bound and partially conditional moment matching approximation as special cases. The PEB approximations are known to lie between a sharp lower bound and an upper bound. Our numerical tests show that the PEB approximations to discrete arithmetic Asian option prices can produce highly accurate approximations when compared to other approximation methods. Our proposed approximation methods can be readily applied to pricing Asian options under most common types of underlying asset price processes, like the Heston stochastic volatility model nested in the class of time-changed Lévy processes with the leverage effect.     

Fast and realistic European ARCH option pricing and hedging

The behavior of the implied volatility surface for European options was analysed in detail by Zumbach and Fernandez for prices computed with a new option pricing scheme based on the construction of the risk-neutral measure for realistic processes with a finite time increment. The resulting dynamics of the surface is static in the moneyness direction, and given by a volatility forecast in the time-to-maturity direction. This difference is the basis of a cross-product approximation of the surface. The subsequent speed-up for option pricing is large, allowing the computation of Greeks and the delta replication strategy in simulations with the cost of replication and the replication risk. The corresponding premia are added to the option arbitrage price in order to compute realistic implied volatility surfaces. Finally, the cross-product approximation for realistic prices can be used to analyse European options on the SP500 in depth. The cross-product approximation is used to compute a mean quotient implied volatility, which can be compared with the full theoretical computation. The comparison shows that the cost of hedging and the replication risk premium have contributions to the implied volatility smile that are of similar magnitude to the contribution from the process for the underlying asset.