No Arabic abstract
Efficient sampling for the conditional time integrated variance process in the Heston stochastic volatility model is key to the simulation of the stock price based on its exact distribution. We construct a new series expansion for this integral in terms of double infinite weighted sums of particular independent random variables through a change of measure and the decomposition of squared Bessel bridges. When approximated by series truncations, this representation has exponentially decaying truncation errors. We propose feasible strategies to largely reduce the implementation of the new series to simulations of simple random variables that are independent of any model parameters. We further develop direct inversion algorithms to generate samples for such random variables based on Chebyshev polynomial approximations for their inverse distribution functions. These approximations can be used under any market conditions. Thus, we establish a strong, efficient and almost exact sampling scheme for the Heston model.
We present an option pricing formula for European options in a stochastic volatility model. In particular, the volatility process is defined using a fractional integral of a diffusion process and both the stock price and the volatility processes have jumps in order to capture the market effect known as leverage effect. We show how to compute a martingale representation for the volatility process. Finally, using It^o calculus for processes with discontinuous trajectories, we develop a first order approximation formula for option prices. There are two main advantages in the usage of such approximating formulas to traditional pricing methods. First, to improve computational effciency, and second, to have a deeper understanding of the option price changes in terms of changes in the model parameters.
How to reconcile the classical Heston model with its rough counterpart? We introduce a lifted version of the Heston model with n multi-factors, sharing the same Brownian motion but mean reverting at different speeds. Our model nests as extreme cases the classical Heston model (when n = 1), and the rough Heston model (when n goes to infinity). We show that the lifted model enjoys the best of both worlds: Markovianity and satisfactory fits of implied volatility smiles for short maturities with very few parameters. Further, our approach speeds up the calibration time and opens the door to time-efficient simulation schemes.
We consider a stochastic volatility model with Levy jumps for a log-return process $Z=(Z_{t})_{tgeq 0}$ of the form $Z=U+X$, where $U=(U_{t})_{tgeq 0}$ is a classical stochastic volatility process and $X=(X_{t})_{tgeq 0}$ is an independent Levy process with absolutely continuous Levy measure $ u$. Small-time expansions, of arbitrary polynomial order, in time-$t$, are obtained for the tails $bbp(Z_{t}geq z)$, $z>0$, and for the call-option prices $bbe(e^{z+Z_{t}}-1)_{+}$, $z eq 0$, assuming smoothness conditions on the {PaleGrey density of $ u$} away from the origin and a small-time large deviation principle on $U$. Our approach allows for a unified treatment of general payoff functions of the form $phi(x){bf 1}_{xgeq{}z}$ for smooth functions $phi$ and $z>0$. As a consequence of our tail expansions, the polynomial expansions in $t$ of the transition densities $f_{t}$ are also {Green obtained} under mild conditions.
The Heston stochastic volatility model is a standard model for valuing financial derivatives, since it can be calibrated using semi-analytical formulas and captures the most basic structure of the market for financial derivatives with simple structure in time-direction. However, extending the model to the case of time-dependent parameters, which would allow for a parametrization of the market at multiple timepoints, proves more challenging. We present a simple and numerically efficient approach to the calibration of the Heston stochastic volatility model with piecewise constant parameters. We show that semi-analytical formulas can also be derived in this more complex case and combine them with recent advances in computational techniques for the Heston model. Our numerical scheme is based on the calculation of the characteristic function using Gauss-Kronrod quadrature with an additional control variate that stabilizes the numerical integrals. We use our method to calibrate the Heston model with piecewise constant parameters to the foreign exchange (FX) options market. Finally, we demonstrate improvements of the Heston model with piecewise constant parameters upon the standard Heston model in selected cases.
In a discrete-time financial market, a generalized duality is established for model-free superhedging, given marginal distributions of the underlying asset. Contrary to prior studies, we do not require contingent claims to be upper semicontinuous, allowing for upper semi-analytic ones. The generalized duality stipulates an extended version of risk-neutral pricing. To compute the model-free superhedging price, one needs to find the supremum of expected values of a contingent claim, evaluated not directly under martingale (risk-neutral) measures, but along sequences of measures that converge, in an appropriate sense, to martingale ones. To derive the main result, we first establish a portfolio-constrained duality for upper semi-analytic contingent claims, relying on Choquets capacitability theorem. As we gradually fade out the portfolio constraint, the generalized duality emerges through delicate probabilistic estimations.