---
title: Heston Model Example
url: https://www.emergentmind.com/topics/heston-model-example
type: topic
---

# Heston Model Example

The Heston model is a stochastic volatility framework widely employed in quantitative finance for pricing options and calibrating volatility surfaces. Core to its appeal are semi-closed-form solutions for European claims, robust calibration methodologies, tractable simulation algorithms, and extensions enabling joint calibration and improved volatility smile reproduction. This article provides a detailed exposition of the Heston model through worked examples, highlighting practical implementation workflows, calibration, numerical schemes, and recent advancements.

## 1. Core Heston Model and Semi-Closed-Form Pricing

Under the risk-neutral measure \( \mathbb{Q} \), the dynamics of the spot price \( S_t \) and its instantaneous variance \( V_t \) are governed by:
\[
\begin{cases}
  dS_t = r S_t dt + \sqrt{V_t} S_t dW_{1,t}, \\
  dV_t = \kappa (\theta - V_t) dt + \eta \sqrt{V_t} dW_{2,t}, \\
  d\langle W_1, W_2 \rangle_t = \rho\, dt,
\end{cases}
\]
where \( r \) is the risk-free rate, \( \kappa \) the variance mean-reversion speed, \( \theta \) the long-run variance, \( \eta \) the volatility of volatility, and \( \rho \) the spot-variance correlation [2101.03626].

For a European call, the model admits a semi-closed-form price using characteristic functions:
\[
C(S, V, t; K, T) = S\, P_1 - K e^{-r(T-t)}\, P_2,
\]
where the probabilities \( P_{1,2} \) are given by one-dimensional Fourier inversions of exponential-affine transforms of the log-moneyness [2101.03626]. This tractability makes the Heston model particularly suitable for calibration to market data and fast computation of vanilla prices.

## 2. Maximum Likelihood Estimation and Empirical Calibration

Parameter estimation for the Heston SDEs can be performed using explicit closed-form maximum likelihood estimators (MLEs) based on discretized likelihoods. The calibrated parameters capture mean reversion (\( \kappa \)), long-run variance (\( \theta \)), volatility-of-volatility (\( \sigma \)), and correlation (\( \rho \)). For observations \( v_n \) at discretized times:
\[
\hat\kappa_N = -\,\frac{2\,b_N + c_N\,d_N}{\,T\bigl(d_N\,f_N -4\bigr)\,},\\
\hat\theta_N = \frac{\,b_N\,f_N + 2\,c_N\,}{\,2\,b_N + c_N\,d_N\,},\\
\hat\sigma^2_N= \frac{a_N}{T} \;-\;\frac{b_N^2\,f_N+4\,b_N\,c_N+c_N^2\,d_N}{2T(d_N\,f_N-4)},
\]
where \( a_N, b_N, c_N, d_N, f_N \) are computed from observed variance time series [1403.4893].

Applications to S&P 500 daily data or minute-level equity data yield parameter estimates that are robust under large-sample regimes, provided the canonical condition \( \kappa\theta/\sigma^2 > 1 \) (for Gaussian asymptotics) is satisfied.

## 3. Risk-Neutral Densities and Explicit Approximations

The risk-neutral density (RND) for \( S_T \) implied by the Heston model falls within a scale-family with scale parameter equal to the forward price. Five explicit one-parameter RND families can be used for efficient approximation:

| Family           | RND Parameterization  | Calibration Note                        |
|------------------|----------------------|-----------------------------------------|
| Log-Normal       | variance \( \nu^2 \) | Good fit, matches positive skew         |
| Gamma            | \(\alpha=1/\nu^2\)   | Slightly inferior to Log-Normal/IG      |
| Inverse-Gaussian | mean 1, var \( \nu^2 \) | Comparable to Log-Normal with positive skew |
| Weibull          | (\( \xi, \lambda \)) | Generally less accurate                 |
| Inverse-Weibull  | (\( \xi, \alpha \))  | Good for positive skew                  |

Numerical calibration on real equity options (e.g., AMD) demonstrates that Log-Normal and Inverse-Gaussian RNDs achieve mean squared errors near those of full Heston, and deliver significant computational savings [2101.03626].

## 4. Advanced Calibration: Joint Equity and Volatility Index Fitting

For robust parameter inference, especially the vol-of-vol, practitioners increasingly incorporate volatility index (VIX/VVIX) data. Four approaches for calibrating the vol-of-vol (\( \xi \)) are demonstrated:

- Closed-form: Based on variance moments; fast and accurate within 5% of full PDE solution.
- Log-contract and replication integrals: Use transition densities for exact or approximate estimation.
- PDE-based double-replication: Numerically solve the full Heston PDE for forward-starting option portfolios, extracting VVIX by discrete replication.

Calibration to both SPX and VVIX data stabilizes \( \xi \) estimates, with values for \( \xi \) remaining within a narrow band across methodologies (0.33–0.37), in contrast to the broad variability under SPX-only calibration [2512.19611, 1706.00873].

## 5. Efficient Simulation and Numerical Schemes

Simulation of Heston paths for pricing or calibration relies on accurate schemes due to the nonlinearity of the variance process:

- Adaptive Simulation: The “ADAPT” framework recasts the CIR variance in terms of a bridge of squared-Bessel processes, adaptively refining quadrature intervals to meet a prescribed variance tolerance while maintaining positivity [1111.6067].
- Hybrid Tree-Finite Difference: Combines a recombining Markov chain for \( V_t \) with 1D finite-difference solves for the log-price, allowing fast, stable weak convergence for European and American options [1307.7178].
- Artificial Boundary Methods: Enhanced boundary condition schemes (ApABC, MApABC1/2) facilitate high-accuracy PDE solutions on truncated domains, significantly reducing boundary-induced bias (<0.5% relative error versus standard Heston BCs) [1912.00691].

Table: Numerical accuracy with various boundary conditions [1912.00691]

| Boundary Condition  | Relative Error (%) |
|---------------------|-------------------|
| Heston Standard     | 2.5               |
| MApABC1             | 0.6               |
| MApABC2             | 0.3               |

## 6. Extensions: Stationary Heston and Multiscale Vol-of-Vol

Two notable Heston extensions improve calibration to short-term smile/skew or enable joint VIX/SPX fitting:

- **Stationary Heston**: Samples the initial variance from the model's invariant Gamma law, correcting for insufficient steepness of short-maturity implied volatility skews, and is calibrated on the same market data as the standard model, yielding significantly improved near-expiry fits [2001.03101].
- **Heston Stochastic Vol-of-Vol (SVV)**: Introduces fast and slow stochastic factors into \( \eta \). The first-order perturbative pricing expansion reduces to Heston’s formula plus Fourier-corrective terms constructed from ODE systems in the Fourier variable. Joint calibration to SPX and VIX option data yields MSE reductions of ~50% compared to classic Heston, at minimal incremental computational cost [1706.00873].

## 7. Practical Implementation and Robustness

Robust calibration approaches include:

- Space Mapping: Use a PDE-calibrated surrogate to accelerate the calibration of SDE-based fine models (e.g., Asian options), with convergence to within 90% of residual error in <4 outer iterations [2501.14521].
- Iterative Splitting for PDEs: Decompose the option pricing PDE into tractable Black-Scholes and correction problems, solved with rapid 1D tridiagonal solvers, yielding reduced computational complexity and improved accuracy of Greeks [2003.12934].

Implementation in MATLAB/Python is facilitated by direct translation of the PDE discretization, ODE solvers for Fourier corrections, and adjoint-based gradient descent for PDE-constrained calibration.

---

The Heston model ecosystem thus supports semi-closed-form vanilla pricing, stable calibration to rich data sources, robust simulation schemes, and further admits powerful extensions for joint smile/skew and volatility-index calibration, making it a cornerstone of contemporary quantitative finance practice [2101.03626, 1706.00873, 2512.19611, 1307.7178].

Source: https://www.emergentmind.com/topics/heston-model-example