Passport options with stochastic volatility

A passport option is a call option on the profits of a trading account. In this article, the robustness of passport option pricing is investigated by incorporating stochastic volatility. The key feature of a passport option is the holders' optimal strategy. It is known that in the case of exponential Brownian motion the strategy is to be long if the trading account is below zero and short if the account is above zero. Here this result is extended to models with stochastic volatility where the volatility is defined via an autonomous SDE. It is shown that if the Brownian motions driving the underlying asset and the volatility are independent then the form of the optimal strategy remains unchanged. This means that the strategy is robust to misspecification of the underlying model. A second aim of this article is to investigate some of the biases which become apparent in a stochastic volatility regime. Using an analytic approximation, comparisons are obtained for passport option prices using the exponential Brownian motion model and some well-known stochastic volatility models. This is illustrated with numerical examples. One conclusion is that if volatility and price are uncorrelated, then prices are sometimes lower in a model with stochastic volatility than in a model with constant volatility.     

Multifactor Heston's stochastic volatility model for European option pricing

In this study, we extend the multiscale stochastic volatility model of [Fouque J‐P, Lorig MJ, SIAM J Financial Math. 2011;2(1):221‐254] by incorporating a slow varying factor of volatility. The resulting model can be viewed as a multifactor extension of the Heston model with two additional factors driving the volatility levels. An asymptotic analysis consisting of singular and regular perturbation expansions is developed to obtain an approximation to European option prices. We also find explicit expressions for some essential functions that are available only in integral formulas in the work of [Fouque J‐P, Lorig MJ, SIAM J Financial Math. 2011;2(1):221‐254]. This finding basically leads to considerable reduction in computational time for numerical calculation as well as calibration problems. An accuracy result of the asymptotic approximation is also provided. For numerical illustration, the multifactor Heston model is calibrated to index options on the market, and we find that the resulting implied volatility surfaces fit the market data better than those produced by the multiscale stochastic volatility model of [Fouque J‐P, Lorig MJ, SIAM J Financial Math. 2011;2(1):221‐254], particularly for long‐maturity call options.     

Pricing VXX option with default risk and positive volatility skew

This paper proposes and makes a comparative study of alternative models for VXX option pricing. Factors such as mean-reversion, jumps, default risk and positive volatility skew are taken into consideration. In particular, default risk is characterized by jump-to-default framework and the “positive volatility skew” issue is addressed by stochastic volatility of volatility and jumps. Daily calibration is conducted and comparative study of the models is performed to check whether they properly fit market prices and generate reasonable positive volatility skews and deltas. Overall, jump-to-default extended LRJ model with positive correlated stochastic volatility (called JDLRJSV in the paper) serves as the best model in all the required aspects.     

Pricing VIX options in a stochastic vol‐of‐vol model

This paper proposes and makes a study of a new model for volatility index option pricing. Factors such as mean‐reversion, jumps, and stochastic volatility are taken into consideration. In particular, the positive volatility skew is addressed by the jump and the stochastic volatility of volatility. Daily calibration is used to check whether the model fits market prices and generates positive volatility skews. Overall, the results show that the mean‐reverting logarithmic jump and stochastic volatility model (called MRLRJSV in the paper) serves as the best model in all the required aspects. Copyright © 2015 John Wiley & Sons, Ltd.     

High-order computational methods for option valuation under multifactor models

Many of the different numerical techniques in the partial differential equations framework for solving option pricing problems have employed only standard second-order discretization schemes. A higher-order discretization has the advantage of producing low size matrix systems for computing sufficiently accurate option prices and this paper proposes new computational schemes yielding high-order convergence rates for the solution of multi-factor option problems. These new schemes employ Galerkin finite element discretizations with quadratic basis functions for the approximation of the spatial derivatives in the pricing equations for stochastic volatility and two-asset option problems and time integration of the resulting semi-discrete systems requires the computation of a single matrix exponential. The computations indicate that this combination of high-order finite elements and exponential time integration leads to efficient algorithms for multi-factor problems. Highly accurate European prices are obtained with relatively coarse meshes and high-order convergence rates are also observed for options with the American early exercise feature. Various numerical examples are provided for illustrating the accuracy of the option prices for Heston’s and Bates stochastic volatility models and for two-asset problems under Merton’s jump-diffusion model.     

A semi-analytic pricing formula for lookback options under a general stochastic volatility model

In general, the pricing problems of exotic options in finance do not have analytic solutions under stochastic volatility and so it is hard to compute the option prices or at least it requires much of time to compute them. This paper investigates a semi-analytic pricing method for lookback options in a general stochastic volatility framework. The resultant formula is well connected to the Black–Scholes price that is the first term of a series expansion, which makes computing the option prices relatively efficient. Further, a convergence condition for the expansion is provided with an error bound.     

Optimizing bounds on security prices in incomplete markets. Does stochastic volatility specification matter?

We extend and generalize some results on bounding security prices under two stochastic volatility models that provide closed-form expressions for option prices. In detail, we compute analytical expressions for benchmark and standard good-deal bounds. For both models, our findings show that our benchmark results generate much tighter bounds. A deep analysis of the properties of option prices and bounds involving a sensitivity analysis and analytical derivation of Greeks for both option prices and bounds is also presented. These results provide strong practical applications taking into account the relevance of pricing and hedging strategies for traders, financial institutions, and risk managers.     

Numerical Studies on Asymptotics of European Option Under Multiscale Stochastic Volatility

Multiscale stochastic volatilities models relax the constant volatility assumption from Black-Scholes option pricing model. Such models can capture the smile and skew of volatilities and therefore describe more accurately the movements of the trading prices. Christoffersen et al. Manag Sci 55(2):1914–1932 (2009) presented a model where the underlying price is governed by two volatility components, one changing fast and another changing slowly. Chiarella and Ziveyi Appl Math Comput 224:283–310 (2013) transformed Christoffersen’s model and computed an approximate formula for pricing American options. They used Duhamel’s principle to derive an integral form solution of the boundary value problem associated to the option price. Using method of characteristics, Fourier and Laplace transforms, they obtained with good accuracy the American option prices. In a previous research of the authors (Canhanga et al. 2014), a particular case of Chiarella and Ziveyi Appl Math Comput 224:283–310 (2013) model is used for pricing of European options. The novelty of this earlier work is to present an asymptotic expansion for the option price. The present paper provides experimental and numerical studies on investigating the accuracy of the approximation formulae given by this asymptotic expansion. We present also a procedure for calibrating the parameters produced by our first-order asymptotic approximation formulae. Our approximated option prices will be compared to the approximation obtained by Chiarella and Ziveyi Appl Math Comput 224:283–310 (2013).     

Static super-replicating strategies for a class of exotic options

In this paper, we investigate static super-replicating strategies for European-type call options written on a weighted sum of asset prices. This class of exotic options includes Asian options and basket options among others. We assume that there exists a market where the plain vanilla options on the different assets are traded and hence their prices can be observed in the market. Both the infinite market case (where prices of the plain vanilla options are available for all strikes) and the finite market case (where only a finite number of plain vanilla option prices are observed) are considered. We prove that the finite market case converges to the infinite market case when the number of observed plain vanilla option prices tends to infinity.We show how to construct a portfolio consisting of the plain vanilla options on the different assets, whose pay-off super-replicates the pay-off of the exotic option. As a consequence, the price of the super-replicating portfolio is an upper bound for the price of the exotic option. The super-hedging strategy is model-free in the sense that it is expressed in terms of the observed option prices on the individual assets, which can be e.g. dividend paying stocks with no explicit dividend process known. This paper is a generalization of the work of Simon et al. [Simon, S., Goovaerts, M., Dhaene, J., 2000. An easy computable upper bound for the price of an arithmetic Asian option. Insurance Math. Econom. 26 (2–3), 175–184] who considered this problem for Asian options in the infinite market case. Laurence and Wang [Laurence, P., Wang, T.H., 2004. What’s a basket worth? Risk Mag. 17, 73–77] and Hobson et al. [Hobson, D., Laurence, P., Wang, T.H., 2005. Static-arbitrage upper bounds for the prices of basket options. Quant. Fin. 5 (4), 329–342] considered this problem for basket options, in the infinite as well as in the finite market case.As opposed to Hobson et al. [Hobson, D., Laurence, P., Wang, T.H., 2005. Static-arbitrage upper bounds for the prices of basket options. Quant. Fin. 5 (4), 329–342] who use Lagrange optimization techniques, the proofs in this paper are based on the theory of integral stochastic orders and on the theory of comonotonic risks.     

Arithmetic Asian Options under Stochastic Delay Models

Abstract

Motivated by the increasing interest in past-dependent asset pricing models, shown in recent years by market practitioners and prominent authors such as Hobson and Rogers (1998 Hobson, D. and Rogers, L. C. G. 1998. Complete models with stochastic volatility. Mathematical Finance, 8(1): 27–48.  [Google Scholar], Complete models with stochastic volatility, Mathematical Finance, 8(1), pp. 27–48), we explore option pricing techniques for arithmetic Asian options under a stochastic delay differential equation approach. We obtain explicit closed-form expressions for a number of lower and upper bounds and compare their accuracy numerically.     

Market calibration under a long memory stochastic volatility model

In this article, we study a long memory stochastic volatility model (LSV), under which stock prices follow a jump-diffusion stochastic process and its stochastic volatility is driven by a continuous-time fractional process that attains a long memory. LSV model should take into account most of the observed market aspects and unlike many other approaches, the volatility clustering phenomenon is captured explicitly by the long memory parameter. Moreover, this property has been reported in realized volatility time-series across different asset classes and time periods. In the first part of the article, we derive an alternative formula for pricing European securities. The formula enables us to effectively price European options and to calibrate the model to a given option market. In the second part of the article, we provide an empirical review of the model calibration. For this purpose, a set of traded FTSE 100 index call options is used and the long memory volatility model is compared to a popular pricing approach – the Heston model. To test stability of calibrated parameters and to verify calibration results from previous data set, we utilize multiple data sets from NYSE option market on Apple Inc. stock.     

Generalized Lipschitz Spaces on Vilenkin Groups

It is shown that, for certain choices of the defining indices, the generalized Lipschitz spaces on Vilenkin groups are included in certain Figà-Talamanca spaces A_p, and that the Fourier series of functions in the latter spaces converge uniformly. This result includes an extension of the classical Dini test, and endpoint versions of Bernstein's theorem on absolute convergence of Fourier series.     

Option prices under stochastic volatility

The well known Heston model for stochastic volatility captures the reality of the motion of stock prices in our financial market. However, the solution of this model is expressed as integrals in the complex plane and has difficulties in numerical evaluation. Here, we present closed-form solutions for option prices and implied volatilities in terms of series expansions. We show that our theoretical predictions are in remarkably good agreement with numerical solutions of the Heston model of stochastic volatility.     

Option pricing under joint dynamics of interest rates, dividends, and stock prices

This paper proposes a unified framework for option pricing, which integrates the stochastic dynamics of interest rates, dividends, and stock prices under the transversality condition. Using the Vasicek model for the spot rate dynamics, I compare the framework with two existing option pricing models. The main implication is that the stochastic spot rate affects options not only directly but also via an endogenously determined dividend yield and return volatility; consequently, call prices can be decreasing with respect to interest rates.     

A method of constructing conditional bases of power functions in certain Banach spaces of analytic functions

In this note, a, method of constructing conditional exponential bases in some Banach spaces of analytic functions is presented. The method of constructing basic families due to Figà-Talamanca is generalized. Concrete applications and examples of conditional basic families in the space of multipliers of power series with the sequence of Taylor coefficients from l^p, in the space of multipliers of Cauchy-type integrals, and in other spaces are given. Bibliography: 12 titles. Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 247, 1997, pp. 170–199. Translated by S. V. Kislyakov.     

Option pricing and hedging in incomplete market driven by Normal Tempered Stable process with stochastic volatility

This paper develops a distribution class, termed Normal Tempered Stable, by subordinating a drifted Brownian motion through a strictly increasing Tempered Stable process that generalizes the Variance Gamma and the Normal Inverse Gaussian and is used to model the logarithm asset returns. The newly added parameter is to create subclasses for all the distributions discovered in financial market. The empirical test suggests that time series of Technology stock returns in US market reject both the Variance Gamma distribution and the Normal Inverse Gaussian distribution and admit instead another subclass of the Normal Tempered Stable distribution. Furthermore, we introduce stochastic volatilities into the Normal Tempered Stable process and derive explicit formulae for option pricing and hedging by means of the characteristic function based methods. To answer the question of how well different models work in practice, we investigate four models adopting data on daily equity option prices and obtain several findings from the numerical results. To sum up, the Normal Tempered Stable process with stochastic volatility is able to adequately capture implied volatility dynamics and seen as a superior model relative to the jump-diffusion stochastic volatility model, based on the construction methodology that incorporates more sophisticated and flexible jump structure and the systematic and realistic treatment of volatility dynamics. The Normal Tempered Stable model turns out to have the competitive performance in an efficient manner given that it only requires three parameters.     

The calibration of volatility for option pricing models with jump diffusion processes

This paper is devoted to calibrate smooth local volatility surface under jump-diffusion processes. This calibration problem is posed as an inverse problem: given a finite set of observed European option prices, find a local volatility function such that the theoretical option prices matches the observed ones optimally with respect to a prescribed performance criterion. Firstly, we obtain an Euler-Lagrange equation for the calibration problem using Tikhonov regularization method. Then we solve the Euler–Lagrange equation using an iterative algorithm and obtain the volatility. Finally, numerical experiments show the effectiveness of the proposed method.     

Option pricing with mean reversion and stochastic volatility

Many underlying assets of option contracts, such as currencies, commodities, energy, temperature and even some stocks, exhibit both mean reversion and stochastic volatility. This paper investigates the valuation of options when the underlying asset follows a mean-reverting lognormal process with stochastic volatility. A closed-form solution is derived for European options by means of Fourier transform. The proposed model allows the option pricing formula to capture both the term structure of futures prices and the market implied volatility smile within a unified framework. A bivariate trinomial lattice approach is introduced to value path-dependent options with the proposed model. Numerical examples using European options, American options and barrier options demonstrate the use of the model and the quality of the numerical scheme.     

A spectral-collocation method for pricing perpetual American puts with stochastic volatility

Based on the Legendre pseudospectral method, we propose a numerical treatment for pricing perpetual American put option with stochastic volatility. In this simple approach, a nonlinear algebraic equation system is first derived, and then solved by the Gauss-Newton algorithm. The convergence of the current scheme is ensured by constructing a test example similar to the original problem, and comparing the numerical option prices with those produced by the classical Projected SOR (PSOR) method. The results of our numerical experiments suggest that the proposed scheme is both accurate and efficient, since the spectral accuracy can be easily achieved within a small number of iterations. Moreover, based on the numerical results, we also discuss the impact of stochastic volatility term on the prices of perpetual American puts.     

Computation and analysis for a constrained entropy optimization problem in finance

In [T. Coleman, C. He, Y. Li, Calibrating volatility function bounds for an uncertain volatility model, Journal of Computational Finance (2006) (submitted for publication)], an entropy minimization formulation has been proposed to calibrate an uncertain volatility option pricing model (UVM) from market bid and ask prices. To avoid potential infeasibility due to numerical error, a quadratic penalty function approach is applied. In this paper, we show that the solution to the quadratic penalty problem can be obtained by minimizing an objective function which can be evaluated via solving a Hamilton–Jacobian–Bellman (HJB) equation. We prove that the implicit finite difference solution of this HJB equation converges to its viscosity solution. In addition, we provide computational examples illustrating accuracy of calibration.