Willow tree algorithms for pricing Guaranteed Minimum Withdrawal Benefits under jump-diffusion and CEV models

This paper presents the willow tree algorithms for pricing variable annuities with Guaranteed Minimum Withdrawal Benefits (GMWB), where the underlying fund dynamics evolve under the Merton jump-diffusion process or constant-elasticity-of-variance (CEV) process. The GMWB rider gives the policyholder the right to make periodic withdrawals from his policy account throughout the life of the contract. The dynamic nature of the withdrawal policy allows the policyholder to decide how much to withdraw on each withdrawal date, or even to surrender the contract. For numerical valuation of the GMWB rider, we use willow tree algorithms that adopt more effective placement of the lattice nodes based on better fitting of the underlying fund price distribution. When compared with other numerical algorithms, like the finite difference method and fast Fourier transform method, the willow tree algorithms compute GMWB prices with significantly less computational time to achieve a similar level of numerical accuracy. The design of our pricing algorithm also includes an efficient search method for the optimal dynamic withdrawal policies. We perform sensitivity analysis of various model parameters on the prices and fair participating fees of the GMWB riders. We also examine the effectiveness of delta hedging when the fund dynamics exhibit various jump levels.     

Pricing and hedging guaranteed minimum withdrawal benefits under a general Lévy framework using the COS method

This paper extends the Fourier-cosine (COS) method to the pricing and hedging of variable annuities embedded with guaranteed minimum withdrawal benefit (GMWB) riders. The COS method facilitates efficient computation of prices and hedge ratios of the GMWB riders when the underlying fund dynamics evolve under the influence of the general class of Lévy processes. Formulae are derived to value the contract at each withdrawal date using a backward recursive dynamic programming algorithm. Numerical comparisons are performed with results presented in Bacinello et al. [Scand. Actuar. J., 2014, 1–20], and Luo and Shevchenko [Int. J. Financ. Eng., 2014, 2, 1–24], to confirm the accuracy of the method. The efficiency of the proposed method is assessed by making comparisons with the approach presented in Bacinello et al. [op. cit.]. We find that the COS method presents highly accurate results with notably fast computational times. The valuation framework forms the basis for GMWB hedging. A local risk minimisation approach to hedging intra-withdrawal date risks is developed. A variety of risk measures are considered for minimisation in the general Lévy framework. While the second moment and variance have been considered in existing literature, we show that the Value-at-Risk (VaR) may also be of interest as a risk measure to minimise risk in variable annuities portfolios.     

Portfolio optimisation with jumps: Illustration with a pension accumulation scheme

In this paper, we address portfolio optimisation when stock prices follow general Lévy processes in the context of a pension accumulation scheme. The optimal portfolio weights are obtained in quasi-closed form and the optimal consumption in closed form. To solve the optimisation problem, we show how to switch back and forth between the stochastic differential and standard exponentials of the Lévy processes. We apply this procedure to both the Variance Gamma process and a Lévy process whose arrival rate of jumps exponentially decreases with size. We show through a numerical example that when jumps, and therefore asymmetry and leptokurtosis, are suitably taken into account, then the optimal portfolio share of the risky asset is around half that obtained in the Gaussian framework.     

Risk based capital for guaranteed minimum withdrawal benefit

The guaranteed minimum withdrawal benefit (GMWB), which is sold as a rider to variable annuity contracts, guarantees the return of total purchase payment regardless of the performance of the underlying investment funds. The valuation of GMWB has been extensively covered in the previous literature, but a more challenging problem is the computation of the risk based capital for risk management and regulatory reasons. One needs to find the tail distribution of the profit–loss function, which differs from its expected payoff required for pricing the GMWB contract. GMWB has embedded two option-like features: Management fees are proportional to the current value of the policyholder’s account which results in an average price of the account. Thus the contract resembles an Asian option. However, the fees are charged only up to the time of the account hitting zero which resembles a barrier option payoff. Thus the GMWB is mathematically more complicated than Asian or barrier options traded on the financial markets. To the authors’ best knowledge, this is the first paper in the literature to formulate and analyse profit–loss distribution using PDE methods of such a product with intricate option-like features. Our approach is much more efficient than the current market practice of rather intensive and expensive Monte Carlo simulations due to the lack of samples for extreme cases.     

A multivariate Lévy process model with linear correlation

In this paper, we develop a multivariate risk-neutral Lévy process model and discuss its applicability in the context of the volatility smile of multiple assets. Our formulation is based upon a linear combination of independent univariate Lévy processes and can easily be calibrated to a set of one-dimensional marginal distributions and a given linear correlation matrix. We derive conditions for our formulation and the associated calibration procedure to be well-defined and provide some examples associated with particular Lévy processes permitting a closed-form characteristic function. Numerical results of the option premiums on three currencies are presented to illustrate the effectiveness of our formulation with different linear correlation structures.     

Pure jump Lévy processes for asset price modelling

The goal of the paper is to show that some types of Lévy processes such as the hyperbolic motion and the CGMY are particularly suitable for asset price modelling and option pricing. We wish to review some fundamental mathematic properties of Lévy distributions, such as the one of infinite divisibility, and how they translate observed features of asset price returns. We explain how these processes are related to Brownian motion, the central process in finance, through stochastic time changes which can in turn be interpreted as a measure of the economic activity. Lastly, we focus on two particular classes of pure jump Lévy processes, the generalized hyperbolic model and the CGMY models, and report on the goodness of fit obtained both on stock prices and option prices.     

Leveraged Lévy processes as models for stock prices

Adopting a constant elasticity of variance formulation in the context of a general Lévy process as the driving uncertainty we show that the presence of the leverage effect? ?One explanation of the documented negative relation between market volatilities and the level of asset prices (the ‘smile’ or ‘skew’), we term the ‘leverage effect’, argues that this negative relation reflects greater risk taking by the management, induced by a fall in the asset price, with a view of maximizing the option value of equity shareholders. in this form has the implication that asset price processes satisfy a scaling hypothesis. We develop forward partial integro-differential equations under a general Markovian setup, and show in two examples (both continuous and pure-jump Lévy) how to use them for option pricing when stock prices follow our leveraged Lévy processes. Using calibrated models we then show an example of simulation-based pricing and report on the adequacy of using leveraged Lévy models to value equity structured products.     

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.     

Computing the finite-time expected discounted penalty function for a family of Lévy risk processes

Ever since the first introduction of the expected discounted penalty function (EDPF), it has been widely acknowledged that it contains information that is relevant from a risk management perspective. Expressions for the EDPF are now available for a wide range of models, in particular for a general class of Lévy risk processes. Yet, in order to capitalize on this potential for applications, these expressions must be computationally tractable enough as to allow for the evaluation of associated risk measures such as Value at Risk (VaR) or Conditional Value at Risk (CVaR). Most of the models studied so far offer few interesting examples for which computation of the associated EDPF can be carried out to the last instances where evaluation of risk measures is possible. Another drawback of existing examples is that the expressions are available for an infinite-time horizon EDPF only. Yet, realistic applications would require the computation of an EDPF over a finite-time horizon. In this paper we address these two issues by studying examples of risk processes for which numerical evaluation of the EDPF can be readily implemented. These examples are based on the recently introduced meromorphic processes, including the beta and theta families of Lévy processes, whose construction is tailor-made for computational ease. We provide expressions for the EDPF associated with these processes and we discuss in detail how a finite-time horizon EDPF can be computed for these families. We also provide numerical examples for different choices of parameters in order to illustrate how ruin-based risk measures can be computed for these families of Lévy risk processes.     

VaR limits for pension funds: an evaluation

In a series of papers during the last ten years an interest rate theory with models which are driven by Lévy or more general processes has been developed. In this paper we derive explicit formulas for the correlations of interest rates as well as zero coupon bonds with different maturities. The models considered in this general setting are the forward rate (HJM), the forward process and the LIBOR model as well as the multicurrency extension of the latter. Specific subclasses of the class of generalized hyperbolic Lévy motions are studied as driving processes. Based on a data set of parametrized yield curves derived from German government bond prices we estimate correlations. In a second step the empirical correlations are used to calibrate the Lévy forward rate model. The superior performance of the Lévy driven models becomes obvious from the graphs.     

A transform approach to compute prices and Greeks of barrier options driven by a class of Lévy processes

In this paper we propose a transform method to compute the prices and Greeks of barrier options driven by a class of Lévy processes. We derive analytical expressions for the Laplace transforms in time of the prices and sensitivities of single barrier options in an exponential Lévy model with hyper-exponential jumps. Inversion of these single Laplace transforms yields rapid, accurate results. These results are employed to construct an approximation of the prices and sensitivities of barrier options in exponential generalized hyper-exponential Lévy models. The latter class includes many of the Lévy models employed in quantitative finance such as the variance gamma (VG), KoBoL, generalized hyperbolic, and the normal inverse Gaussian (NIG) models. Convergence of the approximating prices and sensitivities is proved. To provide a numerical illustration, this transform approach is compared with Monte Carlo simulation in cases where the driving process is a VG and a NIG Lévy process. Parameters are calibrated to Stoxx50E call options.     

The Minimal Entropy Martingale Measure (MEMM) for a Markov-Modulated Exponential Lévy Model

This paper deals with the characterization problem of the minimal entropy martingale measure (MEMM) for a Markov-modulated exponential Lévy model. This model is characterized by the presence of a background process modulating the risky asset price movements between different regimes or market environments. This allows to stress the strong dependence of financial assets price with structural changes in the market conditions. Our main results are obtained from the key idea of working conditionally on the modulator-factor process. This reduces the problem to studying the simpler case of processes with independent increments. Our work generalizes some previous works in the literature dealing with either the exponential Lévy case or the exponential-additive case.     

Minimal <Emphasis Type="Italic">q</Emphasis>-entropy martingale measures for exponential time-changed Lévy processes

We consider structure preserving measure transforms for time-changed Lévy processes. Within this class of transforms preserving the time-changed Lévy structure, we derive equivalent martingale measures minimizing relative q-entropy. They combine the corresponding transform for the Lévy process with an Esscher transform on the time change. Structure preservation is found to be an inherent property of minimal q-entropy martingale measures under continuous time changes, whereas it imposes an additional restriction for discontinuous time changes.     

Monte Carlo Option Pricing for Tempered Stable (CGMY) Processes

Lévy processes are popular models for stock price behavior since they allow to take into account jump risk and reproduce the implied volatility smile. In this paper, we focus on the tempered stable (also known as CGMY) processes, which form a flexible 6-parameter family of Lévy processes with infinite jump intensity. It is shown that under an appropriate equivalent probability measure a tempered stable process becomes a stable process whose increments can be simulated exactly. This provides a fast Monte Carlo algorithm for computing the expectation of any functional of tempered stable process. We use our method to price European options and compare the results to a recent approximate simulation method for tempered stable process by Madan and Yor (CGMY and Meixner Subordinators are absolutely continuous with respect to one sided stable subordinators, 2005).     

GMWB for life an analysis of lifelong withdrawal guarantees

In this paper, we analyze the latest guarantee feature in the variable annuities market: guaranteed minimum withdrawal benefits for life (GMWB for life) which are also called guaranteed lifelong withdrawal benefits (GLWB). This option gives the client the right to deduct a certain amount annually from the policy??s account value until death??even if a unit-linked account value drops to zero. We show how such products can be analyzed within a general framework presented in Bauer et al. (ASTIN Bull. 38(2):621?C651, 2008). We price the embedded guarantee for different product designs and parameters under both, deterministic and optimal client behavior.     

Marginal consistent dependence modelling using weak subordination for Brownian motions

We present an approach for modelling dependencies in exponential Lévy market models with arbitrary margins originated from time changed Brownian motions. Using weak subordination of Buchmann et al. [Bernoulli, 2017], we face a new layer of dependencies, superior to traditional approaches based on pathwise subordination, since weakly subordinated processes are not required to have independent components considering multivariate stochastic time changes. We apply a subordinator being able to incorporate any joint or idiosyncratic information arrivals. We emphasize multivariate variance gamma and normal inverse Gaussian processes and state explicit formulae for the Lévy characteristics. Using maximum likelihood, we estimate multivariate variance gamma models on various market data and show that these models are highly preferable to traditional approaches. Consistent values of basket-options under given marginal pricing models are achieved using the Esscher transform, generating a non-flat implied correlation surface.     

Regression-based Monte Carlo methods for stochastic control models: variable annuities with lifelong guarantees

We present regression-based Monte Carlo simulation algorithm for solving the stochastic control models associated with pricing and hedging of the guaranteed lifelong withdrawal benefit (GLWB) in variable annuities, where the dynamics of the underlying fund value is assumed to evolve according to the stochastic volatility model. The GLWB offers a lifelong withdrawal benefit, even when the policy account value becomes zero, while the policyholder remains alive. Upon death, the remaining account value will be paid to the beneficiary as a death benefit. The bang-bang control strategy analysed under the assumption of maximization of the policyholder’s expected cash flow reduces the strategy space of optimal withdrawal policies to three choices: zero withdrawal, withdrawal at the contractual amount or complete surrender. The impact on the GLWB value under various withdrawal behaviours of the policyholder is examined. We also analyse the pricing properties of GLWB subject to different model parameter values and structural features.     

Pricing of guaranteed minimum withdrawal benefits in variable annuities under stochastic volatility,stochastic interest rates and stochastic mortality via the componentwise splitting method

This paper values guaranteed minimum withdrawal benefit (GMWB) riders embedded in variable annuities assuming that the underlying fund dynamics evolve under the influence of stochastic interest rates, stochastic volatility, stochastic mortality and equity risk. The valuation problem is formulated as a partial differential equation (PDE) which is solved numerically by employing the operator splitting method. Sensitivity analysis of the fair guarantee fee is performed with respect to various model parameters. We find that (i) the fair insurance fee charged by the product provider is an increasing function of the withdrawal rate; (ii) the GMWB price is higher when stochastic interest rates and volatility are incorporated in the model, compared to the case of static interest rates and volatility; (iii) the GMWB price behaves non-monotonically with changing volatility of variance parameter; (iv) the fair fee increases with increasing volatility of interest rates parameter, and increasing correlation between the underlying fund and the interest rates; (v) the fair fee increases when the speed of mean-reversion of stochastic volatility or the average long-term volatility increases; (vi) the GMWB fee decreases when the speed of mean-reversion of stochastic interest rates or the average long-term interest rates increase. We investigate both static and dynamic (optimal) policyholder's withdrawal behaviours; we present the optimal withdrawal schedule as a function of the withdrawal account and the investment account for varying volatility and interest rates. When incorporating stochastic mortality, we find that its impact on the fair guarantee fee is rather small. Our results demonstrate the importance of correct quantification of risks embedded in GMWBs and provide guidance to product providers on optimal hedging of various risks associated with the contract.     

COS method for option pricing under a regime-switching model with time-changed Lévy processes

We extend the regime-switching model to the rich class of time-changed Lévy processes and use the Fourier cosine expansion (COS) method to price several options under the resulting models. The extension of the COS method to price under the regime-switching model is not straightforward because it requires the evaluation of the characteristic function which is based on a matrix exponentiation which is not an easy task. For a two-state economy, we give an analytical expression for computing this matrix exponential, and for more than two states, we use the Carathéodory–Fejér approximation to find the option prices efficiently. In the new framework developed here, it is possible to allow switches not only in the model parameters as is commonly done in literature, but we can also completely switch among various popular financial models under different regimes without any additional computational cost. Calibration of the different regime-switching models with real market data shows that the best models are the regime-switching time-changed Lévy models. As expected by the error analysis, the COS method converges exponentially and thus outperforms all other numerical methods that have been proposed so far.