The Heston (1993) model is a continuous-time stochastic-volatility option-pricing model in which the instantaneous variance Vt follows a mean-reverting square-root (Cox-Ingersoll-Ross) diffusion while the asset price follows a geometric diffusion driven by Vt. Its two defining features are (i) a diffusion term on variance proportional to Vt, which keeps variance non-negative and gives the process an affine structure, and (ii) a non-zero correlation ρ between the return and variance innovations, which captures the empirical leverage effect. Because the log-price's conditional characteristic function is exponential-affine and available in closed form, European options admit a semi-closed-form price via Fourier inversion — making Heston the tractable workhorse of the affine stochastic-volatility class.
Key Ideas
Square-root variance dynamics.dVt=κ(θ−Vt)dt+σvVtdWtV: mean-reversion speed κ, long-run variance θ, and vol-of-vol σv. The Vt diffusion is what distinguishes Heston from the GARCH-diffusion limit, whose vol-of-vol is ∝V (see GARCH and BEKK-GARCH).
Leverage via correlated Brownian motions.corr(dWV,dWreturn)=ρ, typically ρ<0, generating the asymmetric put-skew that zero-correlation models (e.g. the arithmetic Ornstein-Uhlenbeck model of Stein-Stein) cannot produce.
Affine / characteristic-function pricing. The conditional characteristic function of the log-price is exponential-affine in Vt, so option prices follow from a one-dimensional Fourier inversion rather than simulation.
Non-central chi-square transition law. The square-root process has a known non-central χ2 transition density; the moment generating function of average variance can be read directly off the CIR bond-pricing formula (Stochastic Volatility Option Pricing).
Under the physical measure the joint dynamics are
dSt/St=μdt+VtdWtS,dVt=κ(θ−Vt)dt+σvVtdWtV,corr(dWS,dWV)=ρ.
The variance is non-negative and, under the Feller condition 2κθ≥σv2, strictly positive. The crucial analytical property is that the log-price xt=lnSt has a conditional characteristic function of exponential-affine form,
E[eiuxT∣xt,Vt]=exp(A(u,τ)+B(u,τ)Vt+iuxt),τ=T−t,
where A and B solve a system of Riccati ordinary differential equations with closed-form solutions. A European call is then recovered by Fourier inversion as C=StΠ1−Ke−rτΠ2, with the risk-neutral exercise probabilities Π1,Π2 obtained from integrals of the characteristic function — the "semi-closed-form" solution. In the uncorrelated special case (ρ=0) the same prices follow from the moment generating function of average variance AV=T−1∫0TVsds, which for the square-root model equals exp(N∗(T)+M∗(T)V0) read off the CIR bond price (Stochastic Volatility Option Pricing); the general correlated case requires Heston's characteristic-function route because ST∣AV is no longer lognormal. The convexity of the Black-Scholes Option Pricing price in variance, combined with the dispersion of AV, produces the U-shaped volatility smile analytically.
Why It Matters
Heston is the reference stochastic-volatility model in derivatives pricing because it delivers non-negative variance, a leverage-driven skew, and near-closed-form option prices simultaneously. Empirically, incorporating stochastic volatility is the first-order improvement over Black-Scholes in pricing and hedging S&P 500 options, reducing absolute pricing errors by 25-70%; adding stochastic interest rates or jumps yields only second-order gains (Bakshi-Cao-Chen (1997) Empirical Performance of Alternative Option Pricing Models). The model's two structural parameters — the leverage ρ<0 and the volatility risk premium λv<0 — jointly and analytically explain a cluster of return-volatility "puzzles": the ambiguous sign of the volatility-feedback slope, the stronger leverage asymmetry of implied than Realized Volatility, and the downward bias of implied-volatility forecasts (Bollerslev-Zhou (2006) Volatility Puzzles: A Simple Framework for Gauging Return-Volatility Regressions). As the canonical single-factor affine model it also anchors the broader affine jump-diffusion pricing framework (Sundaresan (2000) Continuous-Time Methods in Finance: A Review and an Assessment).
Feller condition. Whether 2κθ≥σv2 holds at empirically calibrated parameters (variance touching zero) affects simulation schemes and boundary behaviour.
One-factor limitation. A single variance factor cannot jointly fit the level, slope, and curvature of the implied-vol surface, motivating multi-factor and rough-volatility extensions.