We consider an asset whose risk-neutral dynamics are described by a general class of local-stochastic volatility models and derive a family of asymptotic expansions for European-style option prices and implied volatilities. We also establish rigorous error estimates for these quantities. Our implied volatility expansions are explicit; they do not require any special functions nor do they require numerical integration. To illustrate the accuracy and versatility of our method, we implement it under four different model dynamics: constant elasticity of variance local volatility, Heston stochastic volatility, three-halves stochastic volatility, and SABR local-stochastic volatility.

We propose a multi-scale stochastic volatility model in which a fast mean-reverting factor of volatility is built on top of the Heston stochastic volatility model. A singular pertubative expansion is then used to obtain an approximation for European option prices. The resulting pricing formulas are semi-analytic, in the sense that they can be expressed as integrals. Difficulties associated with the numerical evaluation of these integrals are discussed, and techniques for avoiding these difficulties are provided. Overall, it is shown that computational complexity for our model is comparable to the case of a pure Heston model, but our correction brings significant flexibility in terms of fitting to the implied volatility surface. This is illustrated numerically and with option data.

We consider the Cauchy problem associated with a general parabolic partial differential equation in d dimensions. We find a family of closed-form asymptotic approximations for the unique classical solution of this equation as well as rigorous short-time error estimates. Using a boot-strapping technique, we also provide convergence results for arbitrarily large time intervals.