0
$\begingroup$

I'm just curious if it's possible to use forward dynamics and work out the pricing formulas. Does anyone know if there's a reference (paper/url) I can look at?

$\endgroup$
  • $\begingroup$ Is your question how to price European vanilla options on a forward when it follows a Heston process? $\endgroup$ – LocalVolatility Jun 15 '17 at 16:19
  • $\begingroup$ yes. I'm asking whether it's possible to switch from spot dynamics to forward dynamics similar to black-scholes <=>black formula $\endgroup$ – Danny Jun 16 '17 at 13:20
1
$\begingroup$

We define a process $X$ to start at $X_0 = 0$ have the Heston-like $\mathbb{Q}$-dynamics

\begin{eqnarray} \mathrm{d}X_t & = & -\frac{1}{2} V_t \mathrm{d}t + \sqrt{V_t} \mathrm{d}W_t^{(1)},\\ \mathrm{d}V_t & = & \kappa \left( \theta - V_t \right) \mathrm{d}t + \xi \sqrt{V_t} \mathrm{d}W_t^{(2)}. \end{eqnarray}

with $\mathrm{d} \langle W^{(1)}, W^{(2)} \rangle_t = \rho \mathrm{d}t$ and such that $e^{X_t}$ is a $\mathbb{Q}$-martingale. The characteristic function of $X_t$ under $\mathbb{Q}$ is given by

\begin{equation} \phi_{X_t}(\omega) = \exp \left\{ C(\omega, t) + D(\omega, t) V_0 \right\}, \end{equation}

where

\begin{eqnarray} C(\omega, t) & = & \frac{\kappa \theta}{\xi^2} \left[ \left( \kappa - \rho \xi \mathrm{i} \omega - d(\omega) \right) t - 2 \ln \left( \frac{1 - c(\omega) e^{-d(\omega) t}}{1 - c(\omega)} \right) \right],\\ D(\omega, t) & = & \frac{\kappa - \rho \xi \mathrm{i} \omega - d(\omega)}{\xi^2} \left( \frac{1 - e^{-d(\omega) t}}{1 - c(\omega) e^{-d(\omega) t}} \right),\\ c(\omega) & = & \frac{\kappa - \rho \xi \mathrm{i} \omega - d(\omega)}{\kappa - \rho \xi \mathrm{i} \omega + d(\omega)},\\ d(\omega) & = & \sqrt{(\rho\xi \mathrm{i} \omega - \kappa)^2 + \xi^2 \left( \mathrm{i} \omega + \omega^2 \right)}. \end{eqnarray}

This representation of the characteristic function is based on the numerically stable version provided in Albrecher et al. (2007).

Using $X$ you can now construct the underlying assets of interest and immediately obtain the relevant characteristic function. In case of a stock with dividend yield $q$, you set

\begin{equation} S_t = S_0 e^{X_t + (r - q) t} \qquad \Rightarrow \qquad \phi_{\ln \left( S_t \right)}(\omega) = \phi_{X_t}(\omega) e^{\mathrm{i} \omega \left( \ln \left( S_0 \right) + (r - q) t \right)}. \end{equation}

For a forward, you get

\begin{equation} F_t = F_0 e^{X_t} \qquad \Rightarrow \qquad \phi_{\ln \left( F_t \right)}(\omega) = \phi_{X_t}(\omega) e^{\mathrm{i} \omega \ln \left( F_0 \right)} \end{equation}

You can then directly employ your favourite characteristic function-based European vanilla option pricer - e.g. Fang and Oosterlee (2008).

References

Albrecher, Hansjoerg, Philipp Mayer, Wim Shoutens and Jurgen Tistaert (2007) "The Little Heston Trap," Wilmott Magazine, Mar.

Fang, Fang and Cornelis W. Oosterlee (2008) "A Novel Pricing Method for European Options Based on Fourier-Cosine Series Expansions," Siam Journal of Scientific Computing, Vol. 31, No. 2, pp. 826-848

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.