

Singular exotic perturbation
A solution based on local volatility and sensitivities is proposed to calculate exotics' prices
Florian Monciaud and Adil Reghai combine singular perturbation techniques with a price-adjustment argument to analyse the impact of the smile dynamics, ie, the price difference between local stochastic volatility and local volatility on exotic products. They obtain an elegant formula that is exact for vanilla options and they propose a set of well-chosen scenarios to compute the impact efficiently
The main driver when selecting a model for pricing and risk management derivatives products is its capacity to explain the profit and loss (PnL) evolution (Reghai 2015; Reghai & Kettani 2020; Bergomi 2018). Choosing the right model to successfully price and hedge financial instruments is based on a careful study of the financial structure to be considered and the market in which it evolves. The quantitative finance literature initially promoted the local volatility (LV) and then the pure stochastic volatility (SV) models as a means of explaining the observed market smile. However, when we consider the dynamic hedging of exotic products such as autocalls, we rapidly conclude that matching the smile is not enough; we also need to control the way in which the latter evolves when the spot moves. Neither the LV model nor the SV model can describe the smile and its evolution properly. However, a fine-tuned mix between the two gives the flexibility to fit both the vanilla options and the way they evolve when the spot moves. These kinds of models are known in the literature as the local stochastic volatility (LSV) models.
The literature on this topic is vast and covers a diversity of approaches to the definition of the models or the way to calibrate them. Lipton (2002) introduces a universal diffusion model presentation with applications to foreign exchange derivatives. Lipton et al (2014) survey LSV models applied to a variety of first-generation exotics. Many papers cover the different approaches to calibration, as this is one of the most important building blocks of the computation. These include Monte Carlo based approaches (Henry-Labordère 2009) and McKean’s particle method (Guyon & Henry-Labordère 2011). Fouque et al (2011) take a perturbation approach to SV models with one or two factors. However, the calibration of the vanilla is not considered therein. Reghai et al (2012) introduce a mixing weight to control the correlation and the volatility of volatility of the process. The LSV impact is computed for exotics. However, this method does not span all possible stochastic volatility parameters and only works for mild parameters. Hagan et al (2018) perform a singular perturbation analysis on a term-structure stochastic alpha-beta-rho model with fast-varying parameters. However, this is a pure stochastic volatility model.
The objective of this paper is to apply a singular perturbation approach in the case of the LSV model. For that purpose, we use a singular perturbation approach without focusing on the vanilla calibration as described in Fouque et al (2011). We then recover the vanilla fit using a price adjustment, as described in Hagan (2005). The obtained formulas are then computed effectively using a well-chosen scenario, as is done in the computation of the Exotic Theta by Bergomi (2018).
This is the main motivation of this work. We indeed propose an extremely fast algorithm that prices the LSV impact at a much lower computational cost than traditional LSV implementations. It is based on LV prices calculated on a well-chosen volatility scenario. This is not only a game changer for real-time risk management but also a powerful way to infer the stochastic volatility parameters in the presence of exotic prices. One last property of the proposed technique is that it reprices vanilla options perfectly, removing the known burden of LSV calibration. Ultimately, this formula offers a rapid, robust and easy implementation of an essential model in real-time management. Finally, this technique is quite general and could open the door for other industrial applications, which will make it possible to enhance all those existing perturbation formulas developed over the years that did not find industrial applications due to the lack of their equivalent in the presence of the smile.
Problem formulation and main result
In this section, for clarity we first formulate the problem, and then go on to state the main result of the paper.
Formulation of the problem
Assume that the dynamics of the stock are given by the following LSV process:
(1) |
where and satisfies an Orstein-Uhlenbeck process:
(2) |
with .
In order to preserve the vanilla calibration, we make an adjustment, , of the Dupire local volatility , which we define as (Henry-Labordère 2009):
(3) |
The particular choice to represent the mean reversion and the volatility of volatility of the dynamic as and , respectively, highlights the fact that these parameters are large in practice, in order to fit the anticipated breakeven values. At this stage we note that is homogeneous to . This choice is dictated by the fact that has an inverse time dimension, whereas has an inverse of the square root of time dimension.
Statement of the main result
The main contributions of the paper are, first, to provide a methodology that combines perturbation techniques with calibration and, second, to suggest a computation strategy based on exotic Greeks that performs well both theoretically and numerically.
The paper shows a detailed application to the stochastic volatility model. To efficiently compute the LSV impact , we derive the following formula:
(4) |
where:
- •
is the price under the LSV process (1), which can be decomposed as follows:
(5) - •
represents the variance of the invariant distribution of when ;
- •
is defined as pricing under the LV process, which depends on the initial conditions :
(6) - •
denotes a particular exotic Greek.
We recall that any exotic Greek of a payoff is formed from the standard Greek adjusted from the vanilla contribution. More precisely, if we define:
(7) |
represents the quantity of vanilla to be retained in order to hedge volatility surface movement:
We also show below how to compute the exotic Greek through a well-chosen scenario at a very low computational cost. We conclude by showing numerical results on the most-traded instruments in equity derivatives, ie, autocalls.
PnL explain
The most important feature of a model is its ability to explain the PnL evolution on a daily basis. For this exercise, the model is based on the following three pillars:
- •
Option: , which denotes the exposure.
- •
Market: , , which denote the breakeven values of the stock volatility.
- •
Model: , , , , etc, which detail the intrinsic property of the model.
A delta hedged position under the Black-Scholes model gives the following PnL explanatory formula:
(Gamma) |
The fair price is obtained by putting the model parameter to its corresponding realised value , or at least its anticipated level.
Likewise, if we use an advanced model such as an LV or LSV model and perform only a delta hedge strategy, we obtain the following PnL explanatory formula:
(Vega) | ||||
(Gamma) | ||||
(Vanna) | ||||
(Volga) |
Before tackling the fair value pricing through matching the breakeven value, we need to cancel the Vega risk, which is of first-order magnitude and is in practice 20 times more important than second-order terms. By doing so, the PnL explanatory formula changes not only with the first-order term, which disappears, but also with exposures, which are now adjusted by Vega hedging. The new exposures are called exotic exposures as they are nil for vanilla options or any product that is replicable with vanillas.
More precisely, to cancel Vega risk, the trader needs to sell vanilla options (see (7)).
The trader’s new PnL equation is then given by:
(Exotic Gamma) | ||||
(Exotic Vanna) | ||||
(Exotic Volga) |

As an illustration, we can compute the exotic Greeks of an autocall product (figure 1).
The Exotic Vanna for the autocall product (long position) remains positive regardless of the spot scenarios. This has important consequences for the hedging with the LV model, for which the model breakeven point is exactly , and therefore generates a negative carry for the seller of the autocall.

At this stage, we can conclude that hedging an autocall with the LV model will generate a systematic loss on the Exotic Vanna term. A way to compensate for this and retain a fair pricing is to move to a local stochastic volatility one-factor (LSV1F) model that fits the correlation breakeven value (figure 2). This is done via extreme values of stochastic volatility parameters: and close to .
Price-adjustment technique
We start with the following modified Newton lemma.
Lemma 1 (Modified Newton).
Suppose we have a function of a variable . Suppose that is set such that , ie, satisfies a given constraint. The value of on is given by the following formula:
Typically, this calibration can be denoted by . In many dimensions, the adjustment takes the following form:
(8) |
Proof of lemma 1.
First we follow the lines of Newton’s approach by searching for as a perturbation of , ie, . As we can expand as follows:
Then, . Now, we expand to obtain the final result. ∎
This result is the basis of price adjustment in order to fit a given set of constraints. Fitting is another word for calibration. Indeed, the constraint is usually written in finance as follows:
where plays the role of model parameters.
We can adjust our LV model price, in order to match vanilla prices exactly in more classical financial notation using the following vector formula:
where is defined in (7).
The objective of this paper is to combine this idea with the design of a well-chosen scenario that will permit the precise computation of the LSV impact at a very low computational cost, by interpreting the adjustment as an exotic Greek.
Singular exotic perturbation
In this section, we present the method for a singular exotic perturbation. The objective is to solve the pricing dynamic as a function of and see how it converges when for (1).
We denote by the conditional expectation in such a way that we fit the vanilla. It satisfies (3).
Making the expansion in the presence of is hard. Instead, we use the following method, which we name the ‘singular exotic perturbation’:
- •
We perform the singular perturbation without calibrating the vanilla, ie, . The zeroth order gives the Dupire local volatility model. We identify the higher orders as Volga and Vanna contributions.
- •
We apply the modified Newton lemma in order to adjust the expansion and recover an exact calibration.
- •
We explicitly compute zeroth-, first- and second-order order adjustments due to the singular perturbation and correct them in order to maintain the vanilla fit.
- •
Finally, we design well-chosen scenarios in order to simplify the above computations.
Singular perturbation on the non-calibrated process
The non-calibrated dynamics have the following form:
(9)
Let be the price of the derivative. It satisfies the following partial differential equation (PDE):
where . We search for as follows:
We then apply the Poisson formula in (centring condition):
where denotes the integration over the invariant distribution of .
This becomes:
where is the Dupire operator:
Therefore, by applying the Feynman-Kac formula, the first-order equation above can be computed using only the local volatility model and its derivatives:
(I) | ||||
(II) |
The first term in the last equality, (I), shows the full Vanna of the product, summing up all contributions for a comovement of . The second term, (II), has a contribution coming only from vanilla options weighted by the Vega of the product. Therefore, when the product is Vega KT hedged, the first term, which is a Vanna, becomes an Exotic Vanna, and the second term, which is a combination of European contributions, just vanishes. More precisely:111 We use the following equation: .
Probabilistic interpretation
In the computation of the second-order impact, the stochastic volatility is important. We start as if there is no recalibration:
(10) |
where is the invariant distribution of the previous process. Let be as defined in (9). Note that:222 .
Then:
Adjusting the prices to recover the vanilla fit
At this stage, we can apply the modified Newton lemma in order to compensate for the non-calibration. Let us denote by the price obtained using which is non-calibrated. Let us denote by the calibrated (to the vanilla) price. We see that this pricing does not match vanilla options due to the extra term:
We can construct the adjusted price by compensating at first order:
where we have used .
We introduce the ‘Exotic Volga’ Greek:
The Exotic Volga appears naturally for perturbation Greeks once we adjust for the calibration of the vanilla.
Similarly, we adjust the term by introducing the effect of the calibration on the vanilla options and calculate :
We obtain the final result described in (4).
We have obtained the desired result. In particular we have identified the functions that intervene in the expansion, and their exotic nature yields an important property of the formula: it ensures that the vanilla has exactly zero impact on the LSV. This formula of great interest for understanding the effect of stochastic volatility on top of the local volatility. However, a brute force implementation will be needed to provide all the . We can, for example, use an offline computation, as proposed in Hagan (2005).
We could also use algorithmic automatic differentiation (AAD) to compute all the at a cost that does not exceed four price computations. We instead propose the design of particular scenarios in order to compute these exotic Greeks based on two building blocks: the implied Black-Scholes calculator and local volatility pricer. We approach the problem as a computation of the Exotic Theta as presented in Bergomi (2018).
Designing exotic scenarios
In order to compute an exotic version of a sensitivity we:
- •
compute the sensitivity using a classical scenarios approach; and
- •
adjust the scenarios in such a way that the effect on vanilla options disappears.
We use the perturbed local volatility (9).
We therefore obtain an exact definition of the asymmetric bump that makes the vanilla Volga disappear:
(11) |

The ‘Exotic Volga’ Greek is then computed as follows:333 Note that is the price of vanilla options obtained with the local vol process. This is done using a forward PDE sweep.
This Greek is the result of three local volatility prices. One of them has already been computed as it is the central price with no deformation of the volatility surface. The two others are computed by first generating a scenario of volatility deformation using a bump with a value of and then creating an implied volatility bump constructed point by point by inferring the volatility from (11). This construction guarantees that the exotic term is mechanically zero on vanilla options. This property gives it its name: exotic.
It is non-zero if and only if the product is a non-vanilla option, ie, not replicable using a vanilla option.
At this stage we have shown an analytic formula for the LSV impact and given details of the low-complexity scenarios that permit the computation of the exotic in order to implement the formula.
In the next section we shall apply the previous formula to autocalls and show how precise it is for these types of product.
Numerical examples
Twenty-eight different structures are exchanged via the Totem service from IHS Markit:
- •
Seven barrier levels: .
- •
Four maturities: years, years, years, years.
- •
Quarterly coupon equal to 1.25%.
Totem provides the running cost of these structures (figure 3):

We compute the exotic Greeks for different payoffs and then simulate hundreds of random stochastic volatility model parameters . We compare the full LSV impact computed using a full implementation with the formula for the different barrier levels. We show the results in figures 4 and 5.

These exotic Greeks depend only on the products and the volatility surface. They do not depend on the stochastic volatility parameters.
We then use our formula for randomly simulated stochastic volatility parameters and compare the full LSV with the proxy formula. Needless to say, the formula is instantaneous, whereas the full LSV takes a non-negligible time to run. We note that the results are extremely good with no bias (average error 0 basis points) and a standard deviation of around 2bp to be compared with the 200bp we ought to match.
Conclusion
In this paper we have introduced a new methodology: singular exotic perturbation. This is an efficient approach to compute the impact of smile dynamics without running a costly LSV model. It instead builds on a well-chosen scenario priced completely under the simpler LV model. Our proposed formula ensures zero impact on vanilla and performs very well on more complex products such as autocalls. The methodology proposed can be used in different contexts, as it combines three different building blocks: singular perturbation, first-order price adjustment and computability through the introduction of exotic Greeks. We suggest that this methodology be used in other cases such as those involving stochastic rates, multiple assets or correlation skew. Further work is needed to detail these practical applications.
Adil Reghai is head of the quantitative research team for equities, commodities and hybrid business at Natixis Investment Bank’s Groupe BPCE. Florian Monciaud is a quantitative analyst in the same team. Both are based in Paris.
Email: florian.monciaud@gmail.com, adil.reghai@polytechnique.org.
References
- Bergomi L, 2018
Exotic theta
Conference presentation at ‘Quant minds’ - Bergomi L, 2016
Stochastic Volatility Modelling
Wiley - Lipton A, 2002
The vol smile problem
Risk February, pages 81–85, http://www.risk.net/1530435 - Lipton A, A Gal and A Lasis, 2014
Pricing of vanilla and first-generation exotic options in the local stochastic volatility framework: survey and new results
Quantitative Finance 14(11), pages 1,899–1,922 - Fouque JP, G Papanicolaou, R Sircar and K Sølna, 2011
Multiscale Stochastic Volatility for Equity, Interest Rate, and Credit Derivatives
Cambridge University Press - Hagan PS, 2005
Adjusters: turning good prices into great prices
Best of Wilmott , page 297 - Hagan PS, A Lesniewski and D Woodward, 2018
Managing vol surfaces
Wilmott 2018(93), pages 24–43 - Henry-Labordère P, 2009
Calibration of local stochastic volatility models to market smiles: a Monte-Carlo approach
Risk September, pages 112–117 - Guyon J and P Henry-Labordère, 2011
The smile calibration problem solved
SSRN preprint, https://ssrn.com/abstract=1885032 - Reghai A and O Kettani, 2020
Financial Models in Production
Springer - Reghai A, O Kettani, V Klaeyle and A Boukhaffa, 2012
LSV models with a mixing weight
SSRN preprint, https://ssrn.com/abstract=2008207 - Reghai A, 2015
Quantitative Finance: Back to Basics
Palgrave MacMillan
Only users who have a paid subscription or are part of a corporate subscription are able to print or copy content.
To access these options, along with all other subscription benefits, please contact info@risk.net or view our subscription options here: http://subscriptions.risk.net/subscribe
You are currently unable to print this content. Please contact info@risk.net to find out more.
You are currently unable to copy this content. Please contact info@risk.net to find out more.
Copyright Infopro Digital Limited. All rights reserved.
You may share this content using our article tools. Printing this content is for the sole use of the Authorised User (named subscriber), as outlined in our terms and conditions - https://www.infopro-insight.com/terms-conditions/insight-subscriptions/
If you would like to purchase additional rights please email info@risk.net
Copyright Infopro Digital Limited. All rights reserved.
You may share this content using our article tools. Copying this content is for the sole use of the Authorised User (named subscriber), as outlined in our terms and conditions - https://www.infopro-insight.com/terms-conditions/insight-subscriptions/
If you would like to purchase additional rights please email info@risk.net
More on Cutting Edge
Leveraged wrong-way risk
A model to assess the exposure to leveraged and collateralised counterparties is presented
Neural joint S&P 500/VIX smile calibration
A one-factor stochastic local volatility model can solve the joint calibration problem
Pricing the transition of Scope 3 emissions
A framework to measure banks’ costs associated with carbon emissions is proposed
Optimal allocation to cryptocurrencies in diversified portfolios
Asset allocation methods assign positive weights to cryptos in diversified portfolios
Getting more for less: better A / B testing via causal regularisation
A causal machine learning algorithm is used to estimate trades’ price impact
Analytic risk-free rates option pricing with smile and skew
An arbitrage-free short-rate model for backward-looking compounded rates is presented
Smile-consistent basket skew
An analytic approximation for the implied volatility surface of basket options is introduced
Podcast: Artur Sepp on rates volatility and decentralised finance
Quant says high volatility requires pricing and risk management models to be revisited