We know that correlation between empirical volatility and price returns in equities is negative with say but averaging at -0.02. We also know with confidence that volatility has long memory. Heston 1993 model is therefore not appropriate. But it has a nice feature that ready mathematical theory tells us that Heston model’s volatility process (the Cox-Ingersoll-Ross square root process)
is positive when with probability 1. The natural analogue is the fractional stochastic differential equation
and here we have no ready theory telling us positivity. So here we employ a comparison principle. Suppose so that we know that will be positive almost surely. Now consider ; since is positive for positive time, is positive. Now we consider
Now and an example of the graph is the following:
This is a function that is extremely negative for small times quickly converging to zero. Now since is positive and has a derivative that is much smaller than that of the latter must stay positive. This implies that with the same condition as the Heston model, i.e. the fractional Heston model will also stay positive.
By Lemma 3.1 of Wong-Heyde 2006 (Wong-Heyde-2006) for the Heston model we know that
for and .
Working backward from Wong-Heyde’s argument, we first note that
is a local martingale. Then we want to consider a function that solves the Cauchy problem
with boundary condition . The idea is that the solution can be given as an exponential function and therefore is nonnegative, one concludes the local martingale is a supermartingale and hence a true martingale.
The Heston model has an exact solution because in that case the Riccati equation that results from the PDE can be solved exactly. In our case we have can consider leading to
This second equation is of Riccati type
which reduces with the standard substitution to
which simplifies a bit to
for which it is still unclear to me whether there is a closed form solution. A quick numerical solution of the Riccati equation tells us that the solution with initial values and and parameters and is increasing, so we can assume this holds for the boundary condition will also be increasing and hence bounded by . So this gives a rough non-rigorous indication that the Wong-Heyde argument can be repeated to show the existence of an equivalent martingale measure for the fractional Heston model.
from scipy import integrate import matplotlib.pyplot as plt def A(t,kappa,alpha): v1 = mlfn(-kappa*t**alpha,alpha,alpha,1) v2 = mlf(-kappa*t**alpha,alpha,alpha) a = (alpha-1) - kappa*t**alpha * v1/v2 return a def B(t,kappa,alpha): return t**alpha * mlf(-kappa*t**alpha,alpha,alpha) def solvr(Y,t,kappa,alpha): return [Y,-B(t,kappa,alpha)*Y-A(t,kappa,alpha)*Y] kappa=1.2 alpha=0.8 a_t = np.arange(0.01,25,0.01) asol = integrate.odeint(solvr,[0.01,0.01],a_t,args=(kappa,alpha)) plt.plot( a_t, asol ) plt.show()