Heston Characteristic Function FFT Pricing Versus SVI Smile Interpolation

Steven Heston's 1993 square-root stochastic volatility model inspired Jim Gatheral's original SVI hyperbola—in fact, Gatheral proved in 2004 that in the large-maturity limit T → ∞, Heston implied total variance converges exact-asymptotically to Gatheral raw SVI five-parameter total variance! Yet calibrating pure 5-parameter Heston (v0, κ, θ, ξ, ρ) directly to listed option chains requires two numerical breakthroughs covered in our SVI & SABR Volatility Surface Calibration & No-Arbitrage Conditions guide: eliminating complex logarithm branch-cut discontinuities via Albrecher's "Little Heston Trap" formulation, and replacing slow Carr-Madan damped FFT grids with Fang-Oosterlee Fourier-Cosine (COS) or Lewis contour integration.

Eliminating Branch-Cut Explosions via the Little Heston Trap#

Heston's original 1993 paper wrote the log-spot characteristic function φT(u) = E[exp(i·u·ln ST)] with complex term g1 = (b − ρξiu + d) / (b − ρξiu − d) and exp(+d·T) in the numerator. Along real Fourier integration frequencies u > 15 on maturities T > 2 years, principal-value complex square roots flip sign across the negative real axis in C, introducing violent jump discontinuities in the integrand that corrupt call prices by several dollars! Albrecher, Mayer, Schoutens, and Tistaert (2007) inverted the ratio to g2 = 1 / g1 = (b − ρξiu − d) / (b − ρξiu + d) with exp(−d·T), guaranteeing |g2 · exp(−d·T)| < 1 and rock-solid branch continuity for all (u, T)!

Branch-Safe Heston Characteristic Function & COS Pricing Engine#

Carr-Madan FFT vs Lewis Contour vs Fang-Oosterlee COS#

Peter Carr and Dilip Madan (1999) pioneered damped Fast Fourier Transform pricing using αdamp ∈ [0.75, 1.50] across N = 4,096 uniform log-strike nodes. In modern calibration loops, however, Carr-Madan FFT wastes 92% of its grid outside traded OPRA strikes and requires careful tuning of αdamp inside Roger Lee's critical moment strip so E[ST^(1 + αdamp)] stays finite! Fang and Oosterlee's (2008) COS method and Alan Lewis's (2001) Im(u) = −1/2 strip contour evaluate only the exact listed market strikes with exponential O(exp(−c·N)) convergence:

Traders can cross-check individual European strike prices and Greeks from calibrated Heston implied volatilities in our Black-Scholes Option Calculator and Implied Volatility Calculator.

Why Pure Heston Needs Jumps (Bates) or Roughness Below 14 DTE#

When calibrated jointly across 7-DTE, 30-DTE, 90-DTE, and 365-DTE Cboe SPX chains, standard 5-parameter Heston hits a structural wall: because Heston is a continuous Brownian diffusion in both spot and variance, short-maturity ATM skew ∂σBS/∂k|{k=0} converges to a finite constant ρξ / (4√v0) as T → 0, whereas empirical SPX skew explodes as T^(−0.45)! Forcing Heston to fit 7-DTE skew pushes vol-of-vol ξ above 1.40, massively violating Feller's condition 2κθ > ξ2 and ruining 6-month smiles. Desks resolve this either by adding Merton lognormal jumps (Bates 1996), switching to Rough Bergomi fractional volatility, or layering SLV Gyongy leverage multipliers on top of moderate-ξ Heston.

Exact Analytical Variance Swap Fair Strike Under Heston#

One enduring advantage of keeping a calibrated Heston curve alongside SVI on desks regulated by the CFTC and ESMA is exact analytical pricing of OTC SPX variance swaps inside SVI & SABR Volatility Surface Calibration & No-Arbitrage Conditions: expected integrated variance equals E[∫0^T V_t dt] / T = θ + (v0 − θ)·(1 − e^(−κT))/(κT), providing an instant sanity check against discrete log-strip replication.