---
title: Double Heston Model Explained
url: https://www.emergentmind.com/topics/double-heston-model
type: topic
---

# Double Heston Model Explained

The Double Heston model is a two-factor extension of the well-known Heston stochastic volatility model, designed to capture more complex dynamics observed in equity and foreign exchange option markets, particularly those arising from stochastic spot/volatility correlation and volatility-of-volatility clustering. The model prescribes asset price dynamics governed by two independent Cox–Ingersoll–Ross (CIR) variance factors, each with distinct means, mean-reversion speeds, volatilities-of-volatility, and asset-variance correlation structures, producing richer implied volatility surfaces and more realistic exotic option pricing [2602.01376][2601.00815][2501.17100].

## 1. Model Specification and Mathematical Framework

The standard risk-neutral form of the Double Heston model for a non-dividend-paying asset $S_t$ with two stochastic variances $\nu_t^{(1)}, \nu_t^{(2)}$ is:

\[
\begin{aligned}
dS_t &= r S_t dt + \sqrt{\nu_t^{(1)}} S_t dW_t^{(1)} + \sqrt{\nu_t^{(2)}} S_t dW_t^{(2)}, \qquad S_0 = s>0, \\
d\nu_t^{(i)} &= \kappa_i (\theta_i - \nu_t^{(i)}) dt + \sigma_i \sqrt{\nu_t^{(i)}} dW_t^{(2+i)}, \qquad \nu_0^{(i)} = v_i>0, \quad i=1,2,
\end{aligned}
\]

with
\[
\langle W^{(1)}, W^{(2+i)}\rangle_t = \rho_i t, \quad \langle W^{(1)}, W^{(2)}\rangle_t = \langle W^{(2)}, W^{(3)}\rangle_t = \langle W^{(2)}, W^{(4)}\rangle_t = 0.
\]

Key parameters include the risk-free rate $r$, mean-reversions $\kappa_i > 0$, long-run variances $\theta_i > 0$, volatilities-of-volatility $\sigma_i > 0$, and asset-variance correlations $\rho_i \in [-1,1]$. Positivity is enforced by the Feller conditions $2\kappa_i\theta_i \geq \sigma_i^2$ [2601.00815].

A prominent variant focuses on the case where both variance factors share $\kappa$, $\sigma$, but have different $\theta_+, \theta_-$ and crucially, different spot–vol correlations $\rho_+, \rho_-$. The total variance is $v_t = v_t^+ + v_t^-$, and correlation parameters are often expressed as $\rho_\pm = \bar\rho \pm \eta$, thus enabling stochastic spot/volatility correlation effects [2602.01376].

## 2. Calibration and Statistical Estimation

Calibration of the Double Heston model to vanilla options employs joint characteristic function methods and least-squares fits to implied volatility surfaces. The parameter vector $(\kappa, \theta_+, \theta_-, \sigma, \bar\rho, \eta)$ is typically bounded as $\kappa \in [0.1, 20]$, $\theta_{\pm} \in [0.0001,1]$, $\sigma \in [0.01,5]$, $\bar\rho\in(-0.999,0.999)$, $\eta\in[0,0.8]$.

Fourier inversion approaches, such as Carr–Madan single-integral or Lewis’ $P_1,P_2$ representations, are employed. The affine structure yields a joint characteristic function:
\[
\phi(u; \tau) = \exp\left\{A(\tau;u) + B_+(\tau;u)\, v_t^+ + B_-(\tau;u)\, v_t^- + iu X_t\right\},
\]
with explicit formulas for $A,B_\pm$ via Riccati ODEs [2602.01376].

For statistical inference, global parameter estimation is performed via maximum likelihood (MLE) and conditional least squares (CLS) based on continuous-time observations. In the ergodic regime, the drift parameter vector $\tau$ admits closed-form estimators, with strong consistency and asymptotic normality following from martingale central limit theory:
\[
\hat\tau_T = \left(\int_0^T \Lambda R^{-1} \Lambda^\top ds\right)^{-1}\int_0^T \Lambda R^{-1} dZ,
\]
where $\Lambda(Z)$ is a design matrix, and $R(Z)$ is the diffusion covariance [2501.17100].

## 3. Simulation Schemes and Numerical Methods

Exact simulation of the CIR variance factors is achieved via non-central chi-square sampling:
\[
\nu^{(i)}_{t+\Delta t} \mid \nu^{(i)}_t \sim c_i\,\chi^2(\delta_i, \lambda_i),
\]
where $c_i = \frac{\sigma_i^2}{4\kappa_i} (1 - e^{-\kappa_i \Delta t})$, $\delta_i = \frac{4\kappa_i\theta_i}{\sigma_i^2}$, and $\lambda_i = \frac{4\kappa_i e^{-\kappa_i \Delta t}}{\sigma_i^2(1-e^{-\kappa_i \Delta t})}\nu^{(i)}_t$.

The Almost-Exact Simulation (AES) scheme couples this sampling with a partially corrected update for $\ln S$, replacing two of the four stochastic integrals in the discretization by exact CIR-integrals. The AES method substantially outperforms Euler–Maruyama in accuracy and time-step economy, especially when the number of simulation steps matches exercise opportunities in Bermudan and American path-dependent options. For typical scenarios, AES achieves pricing errors $<1.5\%$ for American puts ($M=12$, $N=10^6$) where Euler errors exceed $12\%$, with AES CPU time $150$ s versus $190$ s for Euler [2601.00815].

## 4. Stochastic Spot/Volatility Correlation and Model Dynamics

The Double Heston construction enables stochastic spot/volatility correlation:
\[
\rho_t = \frac{\rho_+ v_t^+ + \rho_- v_t^-}{v_t^+ + v_t^-} = \bar\rho + \eta\, \frac{v_t^+ - v_t^-}{v_t^+ + v_t^-},
\]
so that $\rho_t$ fluctuates within $[\bar\rho - \eta, \bar\rho + \eta]$ as the relative weights of the two variance factors evolve [2602.01376].

Instantaneous correlation between spot and the spot–volatility correlation,
\[
\rho_{c\rho} = \mathrm{Corr}(dX_t, d\rho_t) = \sqrt{\eta^2 - (\rho_t - \bar\rho)^2} \approx \eta\quad \text{(for $\rho_t\approx \bar\rho$)},
\]
translates into a positive “skew–spot beta,” reproducing phenomena such as the risk-reversal beta observed in FX markets.

Exotic derivative pricing is affected: the additivity and stochasticity of correlation generate path-dependent implied skew and realized correlations, which modify one-touch and barrier option pricing as well as volatility swap fair strikes. Numerical experiments show price uplifts for one-touch options of order $1$–$2\%$ for realistic parameters; volatility-swap strikes experience upward shifts of ${\sim}10$ bps as $\eta$ increases from $0$ to $0.6$ [2602.01376].

## 5. Stationarity and Ergodicity

Under “subcritical” conditions ($a_1>0$, $a_2 > \tfrac{\sigma_{12}^2}{2}$, $b_{11}, b_{22} > 0$, $\theta > 0$), the Double Heston process admits a unique invariant distribution. The system exhibits exponential ergodicity; averages of functionals converge rapidly to their stationary means. The stationary joint law can be expressed in terms of affine transforms and is characterized via solutions to two-dimensional Riccati ODEs.

A direct implication is the strong law of large numbers for long-time empirical averages, enabling robust risk management (e.g., long-run realized variance targeting) and statistical inference for model parameters [2501.17100].

## 6. Applications, Performance, and Practical Guidance

The Double Heston model is calibrated to market data using least-squares fits to implied volatilities and global optimizers (e.g., Differential Evolution followed by Levenberg-Marquardt). In single-expiry calibration, $\kappa$ is often fixed for robustness.

AES simulation excels for Bermudan options when the number of time steps matches the number of exercise dates, rendering Monte Carlo option prices virtually bias-free. For American options, adopting a reasonable exercise grid (e.g., 12 or 50 per year) preserves high accuracy. Typical runs use $10^5$–$10^6$ paths.

The model is particularly effective for capturing risk-reversal beta term structure, stochastic skew effects, and providing realistic pricing for barrier products and volatility swaps in foreign exchange and equity option markets [2602.01376][2601.00815].

## 7. Limitations and Best-Practice Recommendations

The AES scheme retains Euler-type bias for two log-price integrals, so in the presence of extreme “vol-of-vol” or correlations, smaller time steps ($\Delta t$) are advised. The scheme generalizes to jump processes or further variance factors, with complexity scaling with the number of CIR draws. For out-of-the-money options with few time steps, relative error may be larger and can be mitigated by moderately increasing $M$.

In summary, the Double Heston model, when paired with advanced simulation and calibration methods, enables a tractable, expressive framework for modeling and pricing a broad class of derivatives with features closely matching observed market phenomena in volatility surfaces and exotics pricing [2602.01376][2601.00815][2501.17100].

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