This study conducts an optimal surrender analysis of reverse mortgage (RM) loans offered to elderly homeowners as a financing option. Recent market evidence on borrower early surrenders has raised concerns about the marketability of RM products and their impact on the program viability. In this article, we derive the borrower optimal surrender strategy as a function of the underlying value of the home used as collateral for RM contracts with tenure payment option. Using a probabilistic approach to American option pricing, we present a decomposition result for the value of the contract as the sum of its European counterpart without the surrendering provision and an early exercise premium. The methodology allows policymakers to assess the financial incentive of their policy design, from which we explain the existing market evidence about borrower rational lapse by means of the resulting surrender boundary and reference probabilities.

We consider the super-replication problem for a class of exotic options known as life-contingent options within the framework of the Black–Scholes market model. The option is allowed to be exercised if the death of the option holder occurs before the expiry date, otherwise there is a compensation payoff at the expiry date. We show that there exists a minimal super-replication portfolio and determine the associated initial investment. We then give a characterisation of when replication of the option is possible. Finally, we give an example of an explicit super-replicating hedge for a simple life-contingent option.

This paper extends the standard double-exponential jump-diffusion (DEJD) model to allow for successive jumps to bring about different effects on the asset price process. The double-exponentially distributed jump sizes are no longer assumed to have the same parameters; instead, we assume that these parameters may take a series of different values to reflect growing or diminishing effects from these jumps. The mathematical analysis of the stock price requires an introduction of a number of distributions that are extended from the hypoexponential (HE) distribution. Under such a generalized setting, the European option price is derived in closed-form which ensures its computational convenience. Through our numerical examples, we examine the effects on the return distributions from the growing and diminishing severity of the upcoming jumps expected in the near future, and investigate how the option prices and the shapes of the implied volatility smiles are influenced by the varying severity of jumps. These results demonstrate the benefits of the modeling flexibility provided by our extension.

We propose a Legendre–Laguerre spectral approximation to price the European and double barrier options in the time-fractional framework. By choosing an appropriate basis function, the spectral discretization is used for the approximation of the spatial derivatives of the time-fractional Black–Scholes equation. For the time discretization, we consider the popular $L1$ finite difference approximation, which converges with order $\mathcal {O}((\Delta \tau )^{2-\alpha })$ for functions which are twice continuously differentiable. However, when using the $L1$ scheme for problems with nonsmooth initial data, only the first-order accuracy in time is achieved. This low-order accuracy is also observed when solving the time-fractional Black–Scholes European and barrier option pricing problems for which the payoffs are all nonsmooth. To increase the temporal convergence rate, we therefore consider a Richardson extrapolation method, which when combined with the spectral approximation in space, exhibits higher order convergence such that high accuracies over the whole discretization grid are obtained. Compared with the traditional finite difference scheme, numerical examples clearly indicate that the spectral approximation converges exponentially over a small number of grid points. Also, as demonstrated, such high accuracies can be achieved in much fewer time steps using the extrapolation approach.

We present an analytical option pricing formula for the European options, in which the price dynamics of a risky asset follows a mean-reverting process with a time-dependent parameter. The process can be adapted to describe a seasonal variation in price such as in agricultural commodity markets. An analytical solution is derived based on the solution of a partial differential equation, which shows that a European option price can be decomposed into two terms: the payoff of the option at the initial time and the time-integral over the lifetime of the option driven by a time-dependent parameter. Finally, results obtained from the formula have been compared with Monte Carlo simulations and a Black–Scholes-type formula under various kinds of long-run mean functions, and some examples of option price behaviours have been provided.

We investigate the European call option pricing problem under the fractional stochastic volatility model. The stochastic volatility model is driven by both fractional Brownian motion and standard Brownian motion. We obtain an analytical solution of the European option price via the Itô’s formula for fractional Brownian motion, Malliavin calculus, derivative replication and the fundamental solution method. Some numerical simulations are given to illustrate the impact of parameters on option prices, and the results of comparison with other models are presented.

In this study, we consider option pricing under a Markov regime-switching GARCH-jump (RS-GARCH-jump) model. More specifically, we derive the risk neutral dynamics and propose a lattice algorithm to price European and American options in this framework. We also provide a method of parameter estimation in our RS-GARCH-jump setting using historical data on the underlying time series. To measure the pricing performance of the proposed algorithm, we investigate the convergence of the tree-based results to the true option values and show that this algorithm exhibits good convergence. By comparing the pricing results of RS-GARCH-jump model with regime-switching GARCH (RS-GARCH) model, GARCH-jump model, GARCH model, Black–Scholes (BS) model, and Regime-Switching (RS) model, we show that accommodating jump effect and regime switching substantially changes the option prices. The empirical results also show that the RS-GARCH-jump model performs well in explaining option prices and confirm the importance of allowing for both jump components and regime switching.

Nowadays many financial derivatives, such as American or Bermudan options, are of early exercise type. Often the pricing of early exercise options gives rise to high-dimensional optimal stopping problems, since the dimension corresponds to the number of underlying assets. High-dimensional optimal stopping problems are, however, notoriously difficult to solve due to the well-known curse of dimensionality. In this work, we propose an algorithm for solving such problems, which is based on deep learning and computes, in the context of early exercise option pricing, both approximations of an optimal exercise strategy and the price of the considered option. The proposed algorithm can also be applied to optimal stopping problems that arise in other areas where the underlying stochastic process can be efficiently simulated. We present numerical results for a large number of example problems, which include the pricing of many high-dimensional American and Bermudan options, such as Bermudan max-call options in up to 5000 dimensions. Most of the obtained results are compared to reference values computed by exploiting the specific problem design or, where available, to reference values from the literature. These numerical results suggest that the proposed algorithm is highly effective in the case of many underlyings, in terms of both accuracy and speed.

In this paper we derive several explicit results on one special sticky diffusion process which is constructed as a time-changed version of a diffusion with no sticky points. A theorem concerning the process-related Green operators defined on some nonnegative piecewise continuous functions is provided. Then, based on this theorem, we explore the distributional properties of the sticky diffusion. A financial application is presented where we compute the value of the European vanilla call option written on the underlying with sticky price dynamics.

We consider the pricing of European options under a modified Black–Scholes equation having fractional derivatives in the “spatial” (price) variable. To be specific, the underlying price is assumed to follow a geometric Koponen–Boyarchenko–Levendorski process. This pure jump Lévy process could better capture the real behaviour of market data. Despite many difficulties caused by the “globalness” of the fractional derivatives, we derive an explicit closed-form analytical solution by solving the fractional partial differential equation analytically, using the Fourier transform technique. Based on the newly derived formula, we also examine, in theory, many basic properties of the option price under the current model. On the other hand, for practical purposes, we impose a reliable implementation method for the current formula so that it can be easily used in the trading market. With the numerical results, the impact of different parameters on the option price are also investigated.

A numerical comparison of the Monte Carlo (MC) simulation and the finite-difference method for pricing European options under a regime-switching framework is presented in this paper. We consider pricing options on stocks having two to four volatility regimes. Numerical results show that the MC simulation outperforms the Crank–Nicolson (CN) finite-difference method in both the low-frequency case and the high-frequency case. Even though both methods have linear growth, as the number of regimes increases, the computational time of CN grows much faster than that of MC. In addition, for the two-state case, we propose a much faster simulation algorithm whose computational time is almost independent of the switching frequency. We also investigate the performances of two variance-reduction techniques: antithetic variates and control variates, to further improve the efficiency of the simulation.

While a classic result by Merton (1973, Bell J. Econ. Manage. Sci., 141–183) is that one should never exercise an American call option just before expiration if the underlying stock pays no dividends, the conclusion of a very recent empirical study conducted by Jensen and Pedersen (2016, J. Financ. Econ.121(2), 278–299) suggests that one should ‘never say never’. This paper complements Jensen and Pedersen's empirical study by presenting a theoretical study on how to price American call options under a hard-to-borrow stock model proposed by Avellaneda and Lipkin (2009, Risk22(6), 92–97). Our study confirms that it is the lending fee that results in the early exercise of American call options and we shall also demonstrate to what extent lending fees have affected the early exercise decision.

A class of finite volume methods is developed for pricing either European or American options under jump-diffusion models based on a linear finite element space. An easy to implement linear interpolation technique is derived to evaluate the integral term involved, and numerical analyses show that the full discrete system matrices are M-matrices. For European option pricing, the resulting dense linear systems are solved by the generalised minimal residual (GMRES) method; while for American options the resulting linear complementarity problems (LCP) are solved using the modulus-based successive overrelaxation (MSOR) method, where the H+-matrix property of the system matrix guarantees convergence. Numerical results are presented to demonstrate the accuracy, efficiency and robustness of these methods.

The notional defined contribution pension scheme combines pay-as-you-go financing and a defined contribution pension formula. The return on contributions is based on an index set by law, such as the growth rate of GDP, average wages or contribution payments. The volatility of this return compromises the system's pension adequacy and therefore guarantees may be needed. Here, we provide a minimum return guarantee to the pension contributions. The price is calculated in a utility indifference framework. We obtain a closed-form solution for a general dependence structure with exponential preferences and in presence of stochastic short interest rates.

In this paper, we study carbon emission trading whose market is gaining in popularity as a policy instrument for global climate change. The mathematical model is presented for pricing options on CO2 emission allowance futures with jump diffusion processes, and a so-called fitted finite volume method is proposed to solve the pricing model for the spatial discretization, in which the Crank-Nicolson is employed for time stepping. In addition, the stability and the convergence of the fully discrete scheme are given, and some numerical results, which are compared with the closed form solution and the Monte Carlo simulation solution, are provided to demonstrate the rates of convergence and the robustness of the numerical method.

The objective of the research presented in this paper is the development of a stochastic adoption threshold. The option pricing approach for modeling investment under uncertainty is extended for the case of comparing two stochastic input prices associated with inputs that are perfect substitutes in a production process. Based on this methodology, a threshold decision rule influenced by the drift and volatility of these two input prices is developed. Theoretical results establish an empirical link for measuring the tradeoff of a relatively more expensive input (biodiesel) with lower price drift and volatility compared with a lower but more volatile priced input (petroleum diesel).

The pricing of option contracts is one of the classical problems in Mathematical Finance. While useful exact solution formulas exist for simple contracts, typically numerical simulations are mandated due to the fact that standard features, such as early-exercise, preclude the existence of such solutions. In this paper we consider derivatives which generalize the classical Black-Scholes setting by not only admitting the early-exercise feature, but also considering assets which evolve by the Constant Elasticity of Variance (CEV) process (which includes the Geometric Brownian Motion of Black-Scholes as a special case). In this paper we investigate a Discontinuous Galerkin method for valuing European and American options on assets evolving under the CEV process which has a number of advantages over existing approaches including adaptability, accuracy, and ease of parallelization.

Transform-based algorithms have wide applications in applied probability, but rarely provide computable error bounds to guarantee the accuracy. We propose an inversion algorithm for two-sided Laplace transforms with computable error bounds. The algorithm involves a discretization parameter C and a truncation parameter N. By choosing C and N using the error bounds, the algorithm can achieve any desired accuracy. In many cases, the bounds decay exponentially, leading to fast computation. Therefore, the algorithm is especially suitable to provide benchmarks. Examples from financial engineering, including valuation of cumulative distribution functions of asset returns and pricing of European and exotic options, show that our algorithm is fast and easy to implement.