An Asymptotic Expansion Approach to Currency Options with a Market Model of Interest Rates under Stochastic Volatility Processes of Spot Exchange Rates

This paper proposes an asymptotic expansion scheme of currency options with a libor market model of interest rates and stochastic volatility models of spot exchange rates. In particular, we derive closed-form approximation formulas for the density functions of the underlying assets and for pricing currency options based on a third order asymptotic expansion scheme; we do not model a foreign exchange rate’s variance such as in Heston [(1993) The Review of Financial studies, 6, 327–343], but its volatility that follows a general time-inhomogeneous Markovian process. Further, the correlations among all the factors such as domestic and foreign interest rates, a spot foreign exchange rate and its volatility, are allowed. Finally, numerical examples are provided and the pricing formula are applied to the calibration of volatility surfaces in the JPY/USD option market.     

Pricing options on discrete realized variance with partially exact and bounded approximations

We derive efficient and accurate analytic approximation formulas for pricing options on discrete realized variance (DRV) under affine stochastic volatility models with jumps using the partially exact and bounded (PEB) approximations. The PEB method is an enhanced extension of the conditioning variable approach commonly used in deriving analytic approximation formulas for pricing discrete Asian style options. By adopting either the conditional normal or gamma distribution approximation based on some asymptotic behaviour of the DRV of the underlying asset price process, we manage to obtain PEB approximation formulas that achieve a high level of numerical accuracy in option values even for short-maturity options on DRV.     

Executive Stock Options

Abstract

As a part of the compensation package many companies provide executives with executive stock options, which are call options with additional restrictions. They provide some financial advantages to the executives and help the company retain the service of the executives who improve the company’s earnings and management.

Until recently the values of the executive stock options were not required to be disclosed in the company?s financial reports. But recent statements from the Financial Accounting Standards Board (FASB) have made it necessary to value these executive stock options. The valuation of executive stock options is also required for investors and financial practitioners. This paper considers the award of performance-based executive stock options when the stock price at the time of stock option award exceeds a given preassigned value. It is assumed that the stock price follows a geometric Brownian motion, and that the number of stock options awarded at any time depends on the stock price at that time.

A valuation formula is derived using the method of Esscher transforms for a multiyear award plan. The closed-form formula derived is similar to the Black-Scholes formula for options and utilizes the standard bivariate normal distribution function, which is available in statistical software. In this paper the number of stock options awarded is assumed to be in a specific form, but the theory presented can be modified to suit other forms of award structure. Moreover, by suitable choice of parameters, a valuation formula is also presented for the award of fixed-value executive stock options grants; this formula is also in a closed form and involves cumulative distribution values of the standard normal random variable. Numerical illustrations of the use of the valuation formulas are presented.     

On the valuation of fader and discrete barrier options in Heston's stochastic volatility model

We focus on closed-form option pricing in Heston's stochastic volatility model, where closed-form formulas exist only for a few option types. Most of these closed-form solutions are constructed from characteristic functions. We follow this closed-form approach and derive multivariate characteristic functions depending on at least two spot values for different points in time. The derived characteristic functions are used as building blocks to set up (semi-) analytical pricing formulas for exotic options with payoffs depending on finitely many spot values such as fader options and discretely monitored barrier options. We compare our result with different numerical methods and examine the computational accuracy.     

Optimal investment and risk control for an insurer with partial information in an anticipating environment

Abstract

This paper considers an optimal investment and risk control problem under the criterion of logarithm utility maximization. The risky asset process and the insurance risk process are described by stochastic differential equations with jumps and anticipating coefficients. The insurer invests in the financial assets and controls the number of policies based on some partial information about the financial market and the insurance claims. The forward integral and Malliavin calculus for Lévy processes are used to obtain a characterization of the optimal strategy. Some special cases are discussed and the closed-form expressions for the optimal strategies are derived.     

Generic pricing of FX,inflation and stock options under stochastic interest rates and stochastic volatility

We consider the pricing of FX, inflation and stock options under stochastic interest rates and stochastic volatility, for which we use a generic multi-currency framework. We allow for a general correlation structure between the drivers of the volatility, the inflation index, the domestic (nominal) and the foreign (real) rates. Having the flexibility to correlate the underlying FX/inflation/stock index with both stochastic volatility and stochastic interest rates yields a realistic model that is of practical importance for the pricing and hedging of options with a long-term exposure. We derive explicit valuation formulas for various securities, such as vanilla call/put options, forward starting options, inflation-indexed swaps and inflation caps/floors. These vanilla derivatives can be valued in closed form under Schöbel and Zhu [Eur. Finance Rev., 1999, 4, 23–46] stochastic volatility, whereas we devise an (Monte Carlo) approximation in the form of a very effective control variate for the general Heston [Rev. Financial Stud., 1993, 6, 327–343] model. Finally, we investigate the quality of this approximation numerically and consider a calibration example to FX and inflation market data.     

On the Malliavin approach to Monte Carlo approximation of conditional expectations

Given a multi-dimensional Markov diffusion X, the Malliavin integration by parts formula provides a family of representations of the conditional expectation E[g(X ₂)|X₁]. The different representations are determined by some localizing functions. We discuss the problem of variance reduction within this family. We characterize an exponential function as the unique integrated mean-square-error minimizer among the class of separable localizing functions. For general localizing functions, we prove existence and uniqueness of the optimal localizing function in a suitable Sobolev space. We also provide a PDE characterization of the optimal solution which allows to draw the following observation : the separable exponential function does not minimize the integrated mean square error, except for the trivial one-dimensional case. We provide an application to a portfolio allocation problem, by use of the dynamic programming principle.Mathematics Subject Classification: 60H07, 65C05, 49-00JEL Classification: G10, C10The authors gratefully acknowledge for the comments raised by an anonymous referee, which helped understanding the existence result of Sect. [4.2] of this paper.     

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.     

Pricing variance and volatility swaps in a stochastic volatility model with regime switching: discrete observations case

This study presents a set of closed-form exact solutions for pricing discretely sampled variance swaps and volatility swaps, based on the Heston stochastic volatility model with regime switching. In comparison with all the previous studies in the literature, this research, which obtains closed-form exact solutions for variance and volatility swaps with discrete sampling times, serves several purposes. (1) It verifies the degree of validity of Elliott et al.'s [Appl. Math. Finance, 2007, 14(1), 41–62] continuous-sampling-time approximation for variance and volatility swaps of relatively short sampling periods. (2) It examines the effect of ignoring regime switching on pricing variance and volatility swaps. (3) It contributes to bridging the gap between Zhu and Lian's [Math. Finance, 2011, 21(2), 233–256] approach and Elliott et al.'s framework. (4) Finally, it presents a semi-Monte-Carlo simulation for the pricing of other important realized variance based derivatives.     

Computations of Greeks in a market with jumps via the Malliavin calculus

Using the Malliavin calculus on Poisson space we compute Greeks in a market driven by a discontinuous process with Poisson jump times and random jump sizes, following a method initiated on the Wiener space in [5]. European options do not satisfy the regularity conditions required in our approach, however we show that Asian options can be considered due to a smoothing effect of the integral over time. Numerical simulations are presented for the Delta and Gamma of Asian options, and confirm the efficiency of this approach over classical finite difference Monte-Carlo approximations of derivatives.Received: July 2003, Mathematics Subject Classification (1991): 90A09, 90A12, 90A60, 60H07JEL Classification: C15, G12We thank M. Coutaud for contributions to the simulations.     

On an efficient multiple time step Monte Carlo simulation of the SABR model

In this paper, we will present a multiple time step Monte Carlo simulation technique for pricing options under the Stochastic Alpha Beta Rho model. The proposed method is an extension of the one time step Monte Carlo method that we proposed in an accompanying paper Leitao et al. [Appl. Math. Comput. 2017, 293, 461–479], for pricing European options in the context of the model calibration. A highly efficient method results, with many very interesting and nontrivial components, like Fourier inversion for the sum of log-normals, stochastic collocation, Gumbel copula, correlation approximation, that are not yet seen in combination within a Monte Carlo simulation. The present multiple time step Monte Carlo method is especially useful for long-term options and for exotic options.     

Pricing American options under stochastic volatility and stochastic interest rates

We introduce a new analytical approach to price American options. Using an explicit and intuitive proxy for the exercise rule, we derive tractable pricing formulas using a short-maturity asymptotic expansion. Depending on model parameters, this method can accurately price options with time-to-maturity up to several years. The main advantage of our approach over existing methods lies in its straightforward extension to models with stochastic volatility and stochastic interest rates. We exploit this advantage by providing an analysis of the impact of volatility mean-reversion, volatility of volatility, and correlations on the American put price.     

Summation-formulas of Lubbock's type

Abstract

If it is required to calculate a sum consisting of a great many terms, it is natural to ask oneself whether an approximation might not be obtained by adding up every mth term and multiplying the result by m. If the approximation obtained in this way is not considered sufficient, certain supplementary terms must be added to the result, and these may either be expressed as differences or as differential coefficients of the function under consideration. In the former case we have formulas of Lubbock's, in the latter of Woolhouse's type. We here confine our attention to the formulas of Lubbock's type.     

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.     

A constraint-free approach to optimal reinsurance

Reinsurance is available for a reinsurance premium that is determined according to a convex premium principle H. The first insurer selects the reinsurance coverage that maximizes its expected utility. No conditions are imposed on the reinsurer's payment. The optimality condition involves the gradient of H. For several combinations of H and the first insurer's utility function, closed-form formulas for the optimal reinsurance are given. If H is a zero utility principle (for example, an exponential principle or an expectile principle), it is shown, by means of Borch's Theorem, that the optimal reinsurer's payment is a function of the total claim amount and that this function satisfies the so-called 1-Lipschitz condition. Frequently, authors impose these two conclusions as hypotheses at the outset.     

Optimal portfolio choice in the bond market

We consider the Merton problem of optimal portfolio choice when the traded instruments are the set of zero-coupon bonds. Working within a Markovian Heath–Jarrow–Morton model of the interest rate term structure driven by an infinite-dimensional Wiener process, we give sufficient conditions for the existence and uniqueness of an optimal trading strategy. When there is uniqueness, we provide a characterization of the optimal portfolio as a sum of mutual funds. Furthermore, we show that a Gauss–Markov random field model proposed by Kennedy [Math. Financ. 4, 247–258(1994)] can be treated in this framework, and explicitly calculate the optimal portfolio. We show that the optimal portfolio in this case can be identified with the discontinuities of a certain function of the market parameters.     

The survival probability of the SABR model: asymptotics and application

The stochastic-alpha-beta-rho (SABR) model is widely used by practitioners in interest rate and foreign exchange markets. The probability of hitting zero sheds light on the arbitrage-free small strike implied volatility of the SABR model (see, e.g. De Marco et al. [SIAM J. Financ. Math., 2017, 8(1), 709–737], Gulisashvili [Int. J. Theor. Appl. Financ., 2015, 18, 1550013], Gulisashvili et al. [Mass at zero in the uncorrelated SABR modeland implied volatility asymptotics, 2016b]), and the survival probability is also closely related to binary knock-out options. Besides, the study of the survival probability is mathematically challenging. This paper provides novel asymptotic formulas for the survival probability of the SABR model as well as error estimates. The formulas give the probability that the forward price does not hit a nonnegative lower boundary before a fixed time horizon.     

On accurate and provably efficient GARCH option pricing algorithms

The GARCH model has been very successful in capturing the serial correlation of asset return volatilities. As a result, applying the model to options pricing attracts a lot of attention. However, previous tree-based GARCH option pricing algorithms suffer from exponential running time, a cut-off maturity, inaccuracy, or some combination thereof. Specifically, this paper proves that the popular trinomial-tree option pricing algorithms of Ritchken and Trevor (Ritchken, P. and Trevor, R., Pricing options under generalized GARCH and stochastic volatility processes. J. Finance, , 54(1), 377–402.) and Cakici and Topyan (Cakici, N. and Topyan, K., The GARCH option pricing model: a lattice approach. J. Comput. Finance, , 3(4), 71–85.) explode exponentially when the number of partitions per day, n, exceeds a threshold determined by the GARCH parameters. Furthermore, when explosion happens, the tree cannot grow beyond a certain maturity date, making it unable to price derivatives with a longer maturity. As a result, the algorithms must be limited to using small n, which may have accuracy problems. The paper presents an alternative trinomial-tree GARCH option pricing algorithm. This algorithm provably does not have the short-maturity problem. Furthermore, the tree-size growth is guaranteed to be quadratic if n is less than a threshold easily determined by the model parameters. This level of efficiency makes the proposed algorithm practical. The surprising finding for the first time places a tree-based GARCH option pricing algorithm in the same complexity class as binomial trees under the Black–Scholes model. Extensive numerical evaluation is conducted to confirm the analytical results and the numerical accuracy of the proposed algorithm. Of independent interest is a simple and efficient technique to calculate the transition probabilities of a multinomial tree using generating functions.     

“Relative Importance of Risk Sources in Insurance Systems”, Edward W. Frees,April 1998

Abstract

This paper considers a Sparre Andersen collective risk model in which the distribution of the interclaim time is that of a sum of n independent exponential random variables; thus, the Erlang(n) model is a special case. The analysis is focused on the function φ(u), the expected discounted penalty at ruin, with u being the initial surplus. The penalty may depend on the deficit at ruin and possibly also on the surplus immediately before ruin. It is shown that the function φ(u) satisfies a certain integro-differential equation and that this equation can be solved in terms of Laplace transforms, extending a result found in Lin (2003). As a consequence, a closed-form expression is obtained for the discounted joint probability density of the deficit at ruin and the surplus just before ruin, if the initial surplus is zero. For this formula and other results, the roots of Lundberg’s fundamental equation in the right half of the complex plane play a central role. Also, it is shown that φ(u) satisfies Li’s (2003) renewal equation. Under the assumption that the penalty depends only on the deficit at ruin and that the individual claim amount density is a combination of exponential densities, a closed-form expression for φ(u) is derived. In this context, known results of the Cauchy matrix are useful. Surprisingly, certain results are best expressed in terms of divided differences, a topic deleted from the actuarial examinations at the end of last century.