Skip to content

Stochastic volatility

Stochastic volatility makes future variance uncertain. Read the Heston model, separate expected variance from expected volatility, and see why positive theory does not guarantee positive Euler steps.

Stochastic volatility means that an asset's instantaneous volatility evolves randomly through time. The model treats the future intensity of price fluctuations as uncertain, alongside the future price itself. For an options desk, this adds a risk that a constant-volatility model suppresses: tomorrow's volatility state can change the value of a position even when spot barely moves.
Start with volatility and Brownian motion. A stochastic differential equation then describes how today's state responds to drift and new random shocks.

Random price, random variance

In the constant-volatility Black–Scholes model, the stock price is random but its volatility parameter is fixed. Replacing that constant by a known calendar schedule σ(t)\sigma(t) changes the amount of uncertainty over time, but the schedule itself remains deterministic.
A stochastic-volatility model gives volatility, or often variance, its own random dynamics. In the Heston model, the pricing equations with constant risk-free rate rr and continuous dividend yield qq are
dStSt=(r−q) dt+vt dWtS,dvt=κ(θ−vt) dt+ηvt dWtv,dWtS dWtv=ρ dt.\begin{aligned} \frac{dS_t}{S_t} &= (r-q)\,dt+\sqrt{v_t}\,dW_t^S,\\ dv_t &= \kappa(\theta-v_t)\,dt+\eta\sqrt{v_t}\,dW_t^v,\\ dW_t^S\,dW_t^v &= \rho\,dt. \end{aligned}
These equations are under the risk-neutral measure Q\mathbb{Q}. The two Brownian shocks have instantaneous correlation ρ\rho. Gatheral develops this square-root variance specification in Chapter 2 of The Volatility Surface; Heston's original paper explains the role of correlation in return skewness and option pricing.

The distinction between vtv_t and vt\sqrt{v_t} matters:

QuantityWhat it means
vtv_tInstantaneous variance rate: decimal return squared per year when time is in years
vt\sqrt{v_t}Instantaneous volatility: decimal return per square root of a year
θ\thetaLong-run mean variance, not long-run mean volatility
κ>0\kappa>0Speed of variance mean reversion, per year
η>0\eta>0Scale of the variance diffusion, conventionally called “vol of vol”
ρ∈[−1,1]\rho\in[-1,1]Correlation between spot-return and variance shocks

For example, vt=0.04v_t=0.04 corresponds to 20% annualised instantaneous volatility, not 4%. Nor is vtv_t the variance of the price in currency units: the stock's diffusion coefficient is StvtS_t\sqrt{v_t}.

With ρ<0\rho<0, a negative spot shock tends to accompany a positive variance shock. This is a probabilistic relationship between shocks. It does not mean every fall in the stock must raise volatility; the variance drift and the other component of its random shock still matter.

A six-month variance forecast inside the model

Suppose a hypothetical pricing model uses

v0=0.04,κ=2,θ=0.09,η=0.30,ρ=−0.70,v_0=0.04,\qquad \kappa=2,\qquad \theta=0.09,\qquad \eta=0.30,\qquad \rho=-0.70,

with time in years and returns in decimals. These are illustrative parameters, not estimates from a traded asset.

Mean reversion pulls variance toward 0.090.09, but individual paths wander around that pull. Taking the conditional expectation of the variance equation gives the moment equation m′(t)=κ(θ−m(t))m'(t)=\kappa(\theta-m(t)), with m(0)=v0m(0)=v_0. Hence
EQ[vT∣v0]=θ+(v0−θ)e−κT.\mathbb{E}^{\mathbb{Q}}[v_T\mid v_0] =\theta+(v_0-\theta)e^{-\kappa T}.

At T=0.5T=0.5 years,

EQ[v0.5∣v0=0.04]=0.09−0.05e−1≈0.071606.\mathbb{E}^{\mathbb{Q}}[v_{0.5}\mid v_0=0.04] =0.09-0.05e^{-1} \approx0.071606.

Variance starts below its long-run mean, so its conditional mean rises. That statement does not predict a monotone rise on a particular path. The initial state is fixed here; the process is not started in its stationary distribution.

Three synthetic Heston variance paths start at 0.04 and fluctuate around a rising conditional mean. The mean reaches 0.071606 at half a year and tends toward the long-run mean variance 0.09. The vertical axis is decimal return squared per year; the horizontal axis is years.
Three synthetic Heston variance paths start at 0.04 and fluctuate around a rising conditional mean. The mean reaches 0.071606 at half a year and tends toward the long-run mean variance 0.09. The vertical axis is decimal return squared per year; the horizontal axis is years.
Original synthetic example. Three seeded paths use exact CIR transition draws at 104 steps per year; joining the observations does not resolve motion between them. The solid blue curve is the analytic conditional mean, not an average of these three paths. Only the variance process is simulated.

Taking the square root gives

EQ[v0.5∣v0]≈26.76%.\sqrt{\mathbb{E}^{\mathbb{Q}}[v_{0.5}\mid v_0]} \approx26.76\%.
This is a root-mean-variance quantity. It is not the expected instantaneous volatility. The square-root function is concave, so Jensen's inequality gives
EQ[vT∣v0]≤EQ[vT∣v0].\mathbb{E}^{\mathbb{Q}}[\sqrt{v_T}\mid v_0] \leq\sqrt{\mathbb{E}^{\mathbb{Q}}[v_T\mid v_0]}.
For non-degenerate future variance, the inequality is strict. Likewise, θ=30%\sqrt{\theta}=30\% is the square root of the long-run mean variance; it is not the stationary mean volatility. Neither number is automatically the implied volatility of a six-month option: that comes from an option price, with its strike and full model dynamics.

What this changes for pricing and hedging

An option's value now depends on both spot and the variance state. A hedge that offsets a small spot move can still lose when variance changes. With two independent sources of Brownian uncertainty and only stock and cash to trade, spot delta alone generally cannot remove the extra variance risk.

The volatility smile provides information about pricing dynamics, but it does not hand over a unique physical forecast. Parameters fitted to option prices belong to the chosen pricing measure. Parameters estimated from historical observations describe dynamics under the physical measure P\mathbb{P}. A volatility-risk premium can make the variance drift differ between them. Simply replacing a historical stock drift by r−qr-q does not justify carrying the historical variance drift across unchanged.

Also keep three modelling choices separate:

  • A deterministic time schedule σ(t)\sigma(t) is known in advance.
  • Local volatility σ(St,t)\sigma(S_t,t) is a specified function of spot and time. Its realised value is random because spot is random, but it has no separate variance state.
  • Stochastic volatility introduces a random volatility or variance factor, as Heston does above. Heston is one member of this family, not its definition.

Using the blue conditional-mean curve as a deterministic variance schedule would preserve that curve, but discard the uncertainty around it. Equal expected variance curves do not imply equal option prices, forward-smile behaviour or hedge performance.

Positive theory does not make Euler steps positive

For the exact square-root variance process, positive parameters and v0>0v_0>0 give non-negative variance. The Feller condition

2κθ≥η22\kappa\theta\geq\eta^2

makes zero unattainable. Here 2(2)(0.09)=0.362(2)(0.09)=0.36, whereas 0.302=0.090.30^2=0.09, so the condition holds. If it fails, the exact process can reach zero; failure does not authorize negative variance.

An ordinary Euler discretisation has a different issue:

vn+1=vn+κ(θ−vn)Δt+ηvnΔt Z,Z∼N(0,1).v_{n+1}=v_n+\kappa(\theta-v_n)\Delta t +\eta\sqrt{v_n\Delta t}\,Z,\qquad Z\sim\mathcal{N}(0,1).

Because the normal shock is unbounded, an Euler step can be negative even when the Feller condition holds. With the same initial state, a deliberately coarse Δt=0.5\Delta t=0.5 and Z=−3Z=-3, the update is about −0.037279-0.037279. That is an artefact of the scheme, not a valid state of the exact model.

Truncation schemes, other specialised discretisations and exact transition sampling address this in different ways. Their treatment of the variance boundary and their pricing bias require checks. Exact sampling of the variance at observation dates alone does not provide an exact joint simulation of the stock and its integrated variance. Gatheral's Chapter 2 makes this theoretical-versus-numerical distinction explicit.

For the next step, Stochastic Volatility and the Heston Model develops parameter interpretation, pricing routes and simulation. The public page is a preview; the full lesson requires membership.

Sources

Last updated
  • volatility
  • stochastic processes
  • heston
  • variance
  • option pricing
← Back to the glossaryBack to top ↑