The double-mean-reverting model, introduced by Gatheral [(2008). Consistent modeling of SPX and VIX options. In The Fifth World Congress of the Bachelier Finance Society London, July 18], is known to be a successful three-factor model that can be calibrated to both CBOE Volatility Index (VIX) and S&P 500 Index (SPX) options. However, the calibration of this model may be slow because there is no closed-form solution formula for European options. In this paper, we use a rescaled version of the model developed by Huh et al. [(2018). A scaled version of the double-mean-reverting model for VIX derivatives. Mathematics and Financial Economics 12: 495–515] and obtain explicitly a closed-form pricing formula for European option prices. Our formulas for the first and second-order approximations do not require any complicated calculation of integral. We demonstrate that a faster calibration result of the double-mean revering model is available and yet the practical implied volatility surface of SPX options can be produced. In particular, not only the usual convex behavior of the implied volatility surface but also the unusual concave down behavior as shown in the COVID-19 market can be captured by our formula.

We study the numerical Adomian decomposition method for the pricing of European options under the well-known Black–Scholes model. However, because of the nondifferentiability of the pay-off function for such options, applying the Adomian decomposition method to the Black–Scholes model is not straightforward. Previous works on this assume that the pay-off function is differentiable or is approximated by a continuous estimation. Upon showing that these approximations lead to incorrect results, we provide a proper approach, in which the singular point is relocated to infinity through a coordinate transformation. Further, we show that our technique can be extended to pricing digital options and European options under the Vasicek interest rate model, in both of which the pay-off functions are singular. Numerical results show that our approach overcomes the difficulty of directly dealing with the singularity within the Adomian decomposition method and gives very accurate results.

The stochastic volatility jump diffusion model with jumps in both return and volatility leads to a two-dimensional partial integro-differential equation (PIDE). We exploit a fast exponential time integration scheme to solve this PIDE. After spatial discretization and temporal integration, the solution of the PIDE can be formulated as the action of an exponential of a block Toeplitz matrix on a vector. The shift-invert Arnoldi method is employed to approximate this product. To reduce the computational cost, matrix splitting is combined with the multigrid method to deal with the shift-invert matrix-vector product in each inner iteration. Numerical results show that our proposed scheme is more robust and efficient than the existing high accurate implicit-explicit Euler-based extrapolation scheme.

We present two new stochastic volatility models in which option prices for European plain-vanilla options have closed-form expressions. The models are motivated by the well-known SABR model, but use modified dynamics of the underlying asset. The asset process is modelled as a product of functions of two independent stochastic processes: a Cox-Ingersoll-Ross process and a geometric Brownian motion. An application of the models to options written on foreign currencies is studied.