We extend the approach of Carr, Itkin and Muravey, 2021 for getting semi-analytical prices of barrier options for the time-dependent Heston model with time-dependent barriers by applying it to the so-called $lambda$-SABR stochastic volatility model.
In doing so we modify the general integral transform method (see Itkin, Lipton, Muravey, Generalized integral transforms in mathematical finance, World Scientific, 2021) and deliver solution of this problem in the form of Fourier-Bessel series. The weights of this series solve a linear mixed Volterra-Fredholm equation (LMVF) of the second kind also derived in the paper. Numerical examples illustrate speed and accuracy of our method which are comparable with those of the finite-difference approach at small maturities and outperform them at high maturities even by using a simplistic implementation of the RBF method for solving the LMVF.
We continue a series of papers devoted to construction of semi-analytic solutions for barrier options. These options are written on underlying following some simple one-factor diffusion model, but all the parameters of the model as well as the barrie
rs are time-dependent. We managed to show that these solutions are systematically more efficient for pricing and calibration than, eg., the corresponding finite-difference solvers. In this paper we extend this technique to pricing double barrier options and present two approaches to solving it: the General Integral transform method and the Heat Potential method. Our results confirm that for double barrier options these semi-analytic techniques are also more efficient than the traditional numerical methods used to solve this type of problems.
We continue a series of papers where prices of the barrier options written on the underlying, which dynamics follows some one factor stochastic model with time-dependent coefficients and the barrier, are obtained in semi-closed form, see (Carr and It
kin, 2020, Itkin and Muravey, 2020). This paper extends this methodology to the CIR model for zero-coupon bonds, and to the CEV model for stocks which are used as the corresponding underlying for the barrier options. We describe two approaches. One is generalization of the method of heat potentials for the heat equation to the Bessel process, so we call it the method of Bessel potentials. We also propose a general scheme how to construct the potential method for any linear differential operator with time-independent coefficients. The second one is the method of generalized integral transform, which is also extended to the Bessel process. In all cases, a semi-closed solution means that first, we need to solve numerically a linear Volterra equation of the second kind, and then the option price is represented as a one-dimensional integral. We demonstrate that computationally our method is more efficient than both the backward and forward finite difference methods while providing better accuracy and stability. Also, it is shown that both method dont duplicate but rather compliment each other, as one provides very accurate results at small maturities, and the other one - at high maturities.
In this paper we derive semi-closed form prices of barrier (perhaps, time-dependent) options for the Hull-White model, ie., where the underlying follows a time-dependent OU process with a mean-reverting drift. Our approach is similar to that in (Carr
and Itkin, 2020) where the method of generalized integral transform is applied to pricing barrier options in the time-dependent OU model, but extends it to an infinite domain (which is an unsolved problem yet). Alternatively, we use the method of heat potentials for solving the same problems. By semi-closed solution we mean that first, we need to solve numerically a linear Volterra equation of the first kind, and then the option price is represented as a one-dimensional integral. Our analysis shows that computationally our method is more efficient than the backward and even forward finite difference methods (if one uses them to solve those problems), while providing better accuracy and stability.
In this paper we modify the model of Itkin, Shcherbakov and Veygman, (2019) (ISV2019), proposed for pricing Quanto Credit Default Swaps (CDS) and risky bonds, in several ways. First, it is known since the Lehman Brothers bankruptcy that the recovery
rate could significantly vary right before or at default, therefore, in this paper we consider it to be stochastic. Second, to reduce complexity of the model, we treat the domestic interest rate as deterministic, because, as shown in ISV2019, volatility of the domestic interest rate does not contribute much to the value of the Quanto CDS spread. Finally, to solve the corresponding systems of 4D partial differential equations we use a different flavor of the Radial Basis Function (RBF) method which is a combination of localized RBF and finite-difference methods, and is known in the literature as RBF-FD. Results of our numerical experiments presented in the paper demonstrate that the influence of volatility of the recovery rate is significant if the correlation between the recovery rate and the log-intensity of the default is non-zero. Also, the impact of the recovery mean-reversion rate on the Quanto CDS spread could be comparable with the impact due to jump-at-default in the FX rate.
We derive a backward and forward nonlinear PDEs that govern the implied volatility of a contingent claim whenever the latter is well-defined. This would include at least any contingent claim written on a positive stock price whose payoff at a possibl
y random time is convex. We also discuss suitable initial and boundary conditions for those PDEs. Finally, we demonstrate how to solve them numerically by using an iterative finite-difference approach.
In this paper we apply Markovian approximation of the fractional Brownian motion (BM), known as the Dobric-Ojeda (DO) process, to the fractional stochastic volatility model where the instantaneous variance is modelled by a lognormal process with drif
t and fractional diffusion. Since the DO process is a semi-martingale, it can be represented as an Ito diffusion. It turns out that in this framework the process for the spot price $S_t$ is a geometric BM with stochastic instantaneous volatility $sigma_t$, the process for $sigma_t$ is also a geometric BM with stochastic speed of mean reversion and time-dependent colatility of volatility, and the supplementary process $calV_t$ is the Ornstein-Uhlenbeck process with time-dependent coefficients, and is also a function of the Hurst exponent. We also introduce an adjusted DO process which provides a uniformly good approximation of the fractional BM for all Hurst exponents $H in [0,1]$ but requires a complex measure. Finally, the characteristic function (CF) of $log S_t$ in our model can be found in closed form by using asymptotic expansion. Therefore, pricing options and variance swaps (by using a forward CF) can be done via FFT, which is much easier than in rough volatility models.
This paper proposes a numerical method for pricing foreign exchange (FX) options in a model which deals with stochastic interest rates and stochastic volatility of the FX rate. The model considers four stochastic drivers, each represented by an It^{o
}s diffusion with time--dependent drift, and with a full matrix of correlations. It is known that prices of FX options in this model can be found by solving an associated backward partial differential equation (PDE). However, it contains non--affine terms, which makes its difficult to solve it analytically. Also, a standard approach of solving it numerically by using traditional finite--difference (FD) or finite elements (FE) methods suffers from the high computational burden. Therefore, in this paper a flavor of a localized radial basis functions (RBFs) method, RBF--FD, is developed which allows for a good accuracy at a relatively low computational cost. Results of numerical simulations are presented which demonstrate efficiency of such an approach in terms of both performance and accuracy for pricing FX options and computation of the associated Greeks.
This paper describes another extension of the Local Variance Gamma model originally proposed by P. Carr in 2008, and then further elaborated on by Carr and Nadtochiy, 2017 (CN2017), and Carr and Itkin, 2018 (CI2018). As compared with the latest versi
on of the model developed in CI2018 and called the ELVG (the Expanded Local Variance Gamma model), here we provide two innovations. First, in all previous papers the model was constructed based on a Gamma time-changed {it arithmetic} Brownian motion: with no drift in CI2017, and with drift in CI2018, and the local variance to be a function of the spot level only. In contrast, here we develop a {it geometric} version of this model with drift. Second, in CN2017 the model was calibrated to option smiles assuming the local variance is a piecewise constant function of strike, while in CI2018 the local variance is a piecewise linear} function of strike. In this paper we consider 3 piecewise linear models: the local variance as a function of strike, the local variance as function of log-strike, and the local volatility as a function of strike (so, the local variance is a piecewise quadratic function of strike). We show that for all these new constructions it is still possible to derive an ordinary differential equation for the option price, which plays a role of Dupires equation for the standard local volatility model, and, moreover, it can be solved in closed form. Finally, similar to CI2018, we show that given multiple smiles the whole local variance/volatility surface can be recovered which does not require solving any optimization problem. Instead, it can be done term-by-term by solving a system of non-linear algebraic equations for each maturity which is fast.
