I have posted a new research paper to the arXiv titled “Exact Solutions for a GBM-type Stochastic Volatility Model having a Stationary Distribution”. The article may be found here.

Let me explain the title and motivation. The most popular stochastic volatility model is certainly the Heston ’93 model. Its virtues are that it is exactly solvable and it has some very nice features: a mean-reverting, stochastic volatility process that can be negatively correlated with equity returns. This process has the well-known “square-root model” form

\displaystyle Heston: \quad d v_t = (\omega - \theta v_t) \, dt + \xi \,\sqrt{ v_t} \, dW_t,

where v_t is the instantaneous variance rate of the asset in question. We could also write this evolution in terms of the instantaneous volatility \sigma_t = \sqrt{v_t}. However, time series analysis of the volatility of broad-based indices, such as the S&P500 index, suggests that a better specification would have d \sigma_t \sim \xi \sigma_t dW_t, which is a geometric Brownian motion (GBM)-type volatility. Taking d \sigma_t = \xi \sigma_t dW_t is the well-known SABR model. The problem with that one is that the volatility is not mean-reverting and does not admit a stationary density \psi(\sigma). To get a stationary density, you need a mean-reverting drift term. Also, it would be very nice to also have a model that has exact solutions. The XGBM model may be the first to combine all these properties. It is a standard bivariate model for the pair (S_t,\sigma_t), but here I’m just going to write the volatility evolution, which is

\displaystyle XGBM: \quad d \sigma_t = \sigma_t (\omega - \theta \sigma_t) \, dt + \xi \, \sigma_t \, dW_t.

To see how this model is closer to the “real-world” than the Heston ’93 model, take a look at the figures at the end, which show Maximum Likelihood fits to a proxy series for \{\sigma_t\}. The proxy is the daily (annualized) volatility for the S&P500 taken from the Oxford-Man Institutes “realized library”. They maintain a number of estimators — I am using a basic one (rv5) which is simply the daily (intraday) volatility using 5-minute log-return observations. I am using all the data available at the time of this study, which is January 3, 2000 – September 28, 2018 (4705 volatility observations). The first figure shows this volatility time series. You can see there is a maximum of around 140% which is the annualized volatility at the height of the 2007-2008 Financial Crisis.

The next figures show the stationary volatility fits. For these histograms, I am using the annualized volatility in decimal (not percent), so the time series maximum is a histogram entry at \sigma \approx 1.4. That small bump is not really visible, but you can see it is accounted for by the axis extension out to that value.

As one can see, the visual fit is better for XGBM vs. Heston, with corresponding log-likelihoods: 5356 (Heston) vs 5920 (XGBM). In fact, the fit is even better (LL = 6055) for the figure labeled “GARCH -3/2 model”, which is the stationary density corresponding to both the GARCH diffusion model and the 3/2-model. In any event, the point of the exercise is to motivate my new paper and interested readers may find it at the link above.

**Update**. This research paper is now published. The cite is “Lewis, Alan L. (2019), “Exact solutions for a GBM-type stochastic volatility model having a stationary distribution”, Wilmott mag, Vol. 2019, Issue 101, May, 20-41. The editors were kind enough to make it a

cover story. As part of that publication, I wrote an additional Introduction, which is available here: XGBMIntro.WMag.Final.