

Local stochastic volatility: shaken, not stirred
Dominique Bang introduces a novel LSV approach to term distribution modelling
CLICK HERE TO LISTEN TO DOMINIQUE IN CONVERSATION WITH RISK.NET
For the past two decades, local stochastic volatility (LSV) models have been enjoying a great deal of popularity among practitioners, not only as term structure models for the pricing and risk management of complex products (see, for example, Andersen & Piterbarg 2010; Blacher 2001; Lipton 2002), but also for the modelling of vanilla options (see, for example, Berestycki et al (2004) and the widely used stochastic alpha, beta, rho (SABR) formula from Hagan et al (2002)). Non-LSV models have also attracted some interest, eg, the Lévy process class (see Ticot & Charvet (2013) for an application to inflation modelling), but are used far less in practical applications.
In the world of vanilla interest rate derivatives (a particular focus of this article), standard LSV frameworks have recently been put under increasing pressure from changing market environments, be it low or negative rates in major currencies or regulatory stress tests scenarios. To compound these difficulties, there are ever-growing demands on these models’ ability to calibrate to a large number of swaptions, while simultaneously capturing the convexity embedded in constant maturity swap (CMS) products accurately. As a result, the near-standard SABR model volatility marking approach, based on the small-time asymptotics in Hagan et al (2002), is now problematic, as the model is both insufficiently rich to match the market and also plagued by well-known arbitrage problems for high and low option strikes. A number of authors have proposed different approaches to tackle the current challenges. For example, Andreasen & Huge (2013) use an extended SABR model based on a one-step partial differential equation (PDE) method, producing an arbitrage-free interpolator, but with a term distribution that deviates from the original dynamics as time to maturity increases. Balland & Tran (2013) propose a semi-analytical approach that is not guaranteed to be arbitrage free (though it is stated that it will be in practical applications). Antonov et al (2017) present a mixture approach.
While the improvements to the original SABR framework represent meaningful progress, none is, in our opinion, completely satisfactory from a theoretical or a practical perspective, and many involve non-intuitive parametrisation, a possibility of arbitrage, a lack of flexibility or challenges in the implementation. Here, we present a new approach that allows us to mix a generic local volatility (LV) function, parametric or non-parametric, with a pure stochastic volatility (SV) process in an intuitive, tractable, easy-to-implement and arbitrage-free way, in the sense that there are no negative butterflies for a given maturity. One design strength of the approach is that the choice of the LV component can be guided by criteria based on relevance rather than numerical and analytical feasibility.
The rest of this article is organised as follows. We first describe the LSV dynamics, present the fundamental idea of our method and show how to handle the LV component (via a one-dimensional integral), provided that the pure SV process is tractable, either via its cumulative density function (CDF) or, equivalently, via its European options prices. Then, we illustrate the method, pairing a general LV with a normal SABR process. The normal SABR dynamics are somewhat tractable, allowing for option prices through a one-dimensional integral of special functions (see Antonov et al 2017); unfortunately, the existing formulas become computationally too expensive when embedded in our method and combined with an LV process. Thus, we develop new analytical formulas for option prices and probability densities in the normal SABR model. These approximations are arbitrage free by construction, computationally very simple (only a few calls to the Gaussian CDF are required) and accurate, even for maturities as long as 50 years. Next, we make a concrete choice for the LV function and eventually show that the model (which supports negative rates) can be successfully calibrated to a large strip of options and to CMS forward levels. Last, we consider, for variety, an (arbitrary) piecewise linear LV combined with a normal Heston, describe the calibration strategy and provide option-pricing formulas, as a discrete sum of options under the Heston model. Finally, we conclude the paper.
Local stochastic volatility
Notation and dynamics
Our primary objective in this article is to compute European option pricing on some forward process , in the context of a general LSV model of the form:
(1) | ||||
(2) |
In (1) and (2), the LV function , assumed positive, controls the support of the distribution and the dynamics of the model (also known as ‘backbone’). The pure SV component, , is driven by a Brownian motion , with some generic volatility process that will typically depend on a second source of randomness. In most model configurations, the skewness of implied volatilities as a function of strike is jointly determined by the form of and the correlation between and , whereas the convexity of implied volatilities primarily originates with the variance of and, in some cases, with the convexity of . Popular choices for include the lognormal dynamics of SABR, and the mean-reverting Cox-Ingersoll-Ross (CIR) dynamics of the Heston process.
We are not focused on using (1) and (2) in a fully dynamic sense, as needed, for instance, for an arbitrage-free model for path-dependent interest rate options. Instead, we wish to use the model in the same way that the SABR model is used in practice: as a model for marking European option prices at single maturities, likely with different model parameters for each maturity. This will ultimately affect how we treat certain drift terms that arise in theory; our focus is on preserving only time 0 forwards, rather than on imposing hard intertemporal martingale conditions.
Local volatility: a reduced-form approach
In this section, we show how to reduce the pricing of the LSV model into independent components via a replication argument, and we present approximate options pricing formulas (4) and (5) in terms of one-dimensional integrals (using either call and put prices on or the CDF of ).
Lamperti’s harmonic transform and drift approximation.
For a positive and sufficiently regular LV function , we consider the Lamperti transform:
and the associated process . An application of Ito’s lemma to results in:
Intuitively, the Lamperti transform is a tool of choice to disentangle local from pure stochastic volatility: the effect of the local component is then captured via a straight change in co-ordinates, giving rise to a stochastic drift, , that ensures the martingality of . As explained earlier, we are solely concerned with the term distribution of for a given horizon , so we make the pragmatic assumption that , where is an ad hoc constant, fixed at a value that ensures preservation of the forward, in the sense that:
where is the expectation operator under the pricing measure . Accordingly, we adopt the following definition for the term distribution of from now on:
(3) |
The simplifying assumption made on the stochastic drift impacts only mildly the distribution of the underlying. This is certainly the case as time to maturity decreases, where in (3) will converge to the exact solution of the stochastic differential equation (SDE), which is consistent with Berestycki et al (2004). For longer expiries, one can expect the model to exhibit larger differences with the true LSV SDE, but qualitatively the LV parametrisation will have a similar effect on the implied volatility smile. Note that, regardless of the maturity, the model is exact for normal/lognormal dynamics and for pure SV, providing increased confidence for the general case. As per construction, (3) will produce a valid (ie, non-negative) density for and is not restrictive in terms of the resulting probability distribution, since the LV , being generic, can absorb any differences when calibrated.
Option pricing
For any given random variable observed at , we adopt the following notation throughout: (respectively, ) for a call value (respectively, a put value) and for its CDF.
Now, since , it follows that and can be written in terms of call/put options on , via a replication argument:
(4) |
We may also write these equations in terms of , the CDF of :
(5) |
Proof (sketch).
SABR LV model
In this section, we use the normal SABR process as the SV building block in our model. As mentioned earlier, fast closed-form expressions for options on are critical to the practical implementation of (4) or (5), so we present efficient, arbitrage-free approximation formulas for options prices and the CDF for normal SABR dynamics.
Normal SABR
We recall the normal SABR model dynamics sets , where:
Here, , , and and are two independent Brownian motions. From Gyongy’s lemma, we know that there exists a unique, pure LV process with the same marginal distributions as . In lemma 1, we describe explicitly the dynamics of .
Lemma 1.
Consider the process , defined by:
Then, for any , has the same marginal distributions as .
Proof (sketch).
Given the result of Gyongy, we know the squared volatility term in the process is , ensuring consistent marginal distributions between and . An explicit calculation of this expression can be performed by a measure change, along the lines of Balland & Tran (2013). ∎
Lemma 2.
For a given strike , let and set:
We then have:
(6) | ||||
Proof.
Integrating explicitly the relation:
and inverting the resulting expression leads to the reconstruction formula:
(7) |
From the definitions of and , we have:
Thus, subtracting from and reverting back to the initial variable , we see:
which shows is distributed as a sum of two increasing functions in , vanishing simultaneously at . Thus, we clearly have:
(8) | ||||
(9) |
Applying the expectation operator then proves the lemma. ∎
We remark that the elegant technique used to prove (8)–(9) is known as Jamshidian’s trick, which was originally used in the context of swaption valuation in a one-factor model. Our next goal is to work out accurate and non-arbitrageable closed-form approximations for , and . We start by using Ito’s lemma to note that:
(10) |
with:
The dynamics of in (10) resemble a mean-reverting drifted Brownian motion, so both its mean and variance will be affected by the (locally) stochastic drift term. Given is defined as an exponential of , one is tempted to use a simple lognormal approximation for , with moment-matched mean and variance. Indeed, we find that this works quite well but nonetheless present an alternative in proposition 1, which boasts even better precision. The idea behind the proposition is to apply a measure change to capture the subtle dynamics of more accurately, without sacrificing tractability.
Proposition 1.
The following (arbitrage-free) expressions can be used as very accurate closed-form solutions for option prices and the CDF for the process :
(11) | ||||
(12) | ||||
(13) |
with:
(14) | ||||
(15) | ||||
(16) | ||||
(17) |
and where the auxiliary functions , and have been defined below via the Gaussian kernel and re-expressed, using elementary integral manipulations, in terms of the normal CDF, :
(18) |
with:
(19) |
Proof.
We consider the exponential martingale in the measure :
The associated measure is defined via:
By construction, is a driftless Brownian motion under . Denoting by the Dirac distribution at :
where denotes expectation in measure . Defining:
and applying Ito’s lemma, we find . Moreover, using a second-order Taylor expansion in , we see:
Thus, we can write:
The constant is introduced to ensure:
Finally, we can write , where the density of under is:
(20) |
can be found by enforcing , which leads, after standard integral manipulations, to the expression in (15).
The parameter , which has been introduced to account for the small measure change approximations made on the way, can be found using the martingale condition in (7), that is:
where:
can be conveniently expressed as in (17). Thus, is the solution of a second-order equation with a unique admissible solution:
Note (20) is a valid density by construction and completely tractable. As a matter of fact, since , and from the definition of in (14), we have:
which can be expressed as in (11) using the function defined in (18). The proof for (12) is similar. For (13), we note that:
which yields the result. ∎
Bringing together results from lemma 2 and proposition 1, we arrive at explicit and fast approximations for options on . We emphasise that these expressions are arbitrage-free by construction, for any set of parameters and for any maturity. Empirically, the approximations perform surprisingly well, even for large maturities, as is shown by the results in figure 1; this compares our approximation to both the Hagan expansion in Hagan et al (2002) and to exact values in Antonov et al (2017) (computed by an expensive integral-based closed-form solution). The new formulas only require a few calls to the Gaussian CDF (with a cost similar to that of the Hagan formula).

Local volatility specification
Now the problem of determining fast and accurate approximations for the normal SABR model is out of the way, we can freely design the LV function in (1) to provide trading desks, using a model that can be fine-tuned to accommodate their specific requirements, for instance:
- •
adequate support, eg, the distribution for should allow for negative rates;
- •
flexibility, the parametrisation needs to be generic enough to be calibrated to a reasonably large number of liquid option strikes (eg, five);
- •
control over the dynamics between level and volatility of the underlying for risk computation purposes (backbone); this should comply with trader requirements and empirical evidence; and
- •
control over model behaviour at high and low strikes and, ultimately, over convex products such as CMS.
While many different LV specifications could be used to target these requirements, let us give a reasonable example for practical purposes. To address the first point, we choose an LV function that remains positive for negative values of . To this aim, for , we introduce the function , which regularises . To control the backbone, we consider a constant elasticity of variance (power) specification, with level-dependent power:
Note varies smoothly from to , with a speed controlled by the constant . So equipped, we are now able to explicitly define the LV function:
In the function , the parameter controls the smile for low strikes and typically the probability mass below 0. If , then remains above 0 (if it started above 0). is an arbitrary strike outside the region of observable strikes, but beyond which the smile still affects CMS forwards. The parameter controls the volatility beyond , and the parameter is used as a smoothing parameter. To avoid explosive behaviour from in the wings, we finally choose to extrapolate the LV flat outside some interval . One benefit of this extrapolation is the integrals in (4) are truncated over finite support, making them more amenable to fast numerical integration. For interest rate applications, a possible parameterisation for , and is shown in figure 2. The final expression for the LV reads:
(21) |
In figure 1 we demonstrate the impact of the parameters controlling the wings of the LV.

Model settings and numerical results
We consider the LV function defined in (21) via a set of LV parameters and the associated Lamperti transform . can be computed and cached, once and for all, using, for instance, a trapezoid quadrature rule on a fine grid. Options on struck at are then expressed in (4), where is computed using lemma 2 and proposition 1; put options may be computed using put-call parity.
For some fixed maturity , our first task is calibrating and to simultaneously reprice the forward (or, equivalently, ensuring ) and the at-the-money (ATM) option, , where is the normal (Bachelier) ATM implied volatility. This is achieved using ATM call and put values, derived from setting in (4). Note that, once calibrated, the method allows us to price a strip of options simultaneously by adequately discretising the integrals in (4). In the next paragraph, the model is (successfully) calibrated to market quotes (a strip of options and a CMS forward).
Calibration to market instruments
In this section, our model is used with a full set of LV and SV parameters, calibrated to best-fit option prices and CMS forwards seen in market quotes. In figure 3, we display the probability density function (PDF) and the normal (Bachelier) volatility implied by the model, using the normal SABR formula (6) and the full LV function (21). Evidently, a very accurate calibration to these instruments can be achieved, with the resulting implied PDF being smooth and well behaved.

Heston model with linear piecewise local volatility
As mentioned earlier, the proposed method for building LSV term volatility models is generic and can accommodate different dynamics. In this section, we briefly illustrate the breadth of the model, using it to equip the Heston model with a general piecewise linear LV. For this purpose, we consider, as a building block, a normal Heston dynamics for the process :
where , and control the ATM level, the skew and the convexity of the implied volatility of the model, respectively.
The characteristic function of , , is known in closed form; thus, the Lewis-Lipton formula can be employed to recover option prices on or on via an inverse Fourier integral (see Lipton 2002). The general pricing equation (4) would not be computationally workable, as it would involve double integrals. Nevertheless, in the next section we show how the assumption of piecewise linearity on the LV allows us to write options prices as a discrete sum of options in the standard Heston model.
For later use, let us denote:
with the convention , (respectively, ) for a call (respectively, for a put).
Replication formula using calls and puts
Now, assuming a piecewise linear LV (extrapolated flat, for concreteness), we consider a strip of strikes ; without loss of generality, we can assume the forward , with . We also define the associated volatilities . Let:
be the LV slopes, and then define the LV function:
We also denote as the (Lamperti) transformed variables on the strike grid. The following identities hold:
(22) |
Calibration
The calibration problem involves the constraints of repricing the forward and reproducing the ATM volatility . The two parameters available for this task are the drift and the initial value of the volatility, :
Thus, the calibration problem can be written as the system:
which can be solved in and by any standard two-dimensional root search algorithm.
Option pricing
Consider a call struck at . We have:
Note has already been computed at the calibration stage. Thus, it is straightforward to price an option on struck at . Note that in the degenerate case where the normal Heston process is Gaussian (), the model can be used as an arbitrage-free LV interpolator (see Schlenkrich (2018) for details).
Conclusion
We have presented a new method for term distribution modelling that allows the user to mix an arbitrary LV with an arbitrary SV in a tractable, intuitive and arbitrage-free way. We have applied the method to the SABR models family, pairing a flexible LV with a normal SABR SV model. In the process, we developed new, efficient closed-form option pricing formulas for the normal SABR SV model. The resulting LSV model is smooth and arbitrage free, and it successfully calibrates to market quotes for swaptions and CMS forwards. Finally, we have also shown how the method can be used in the context of the Heston model by providing closed-form pricing solutions when the LV is a general linear piecewise function and the SV is a normal Heston model.
Dominique Bang is head of interest rate vanilla analytics at Bank of America Merrill Lynch in London. He is indebted to Leif Andersen and Guillaume Blacher for insightful suggestions, comments and support, and he is very grateful to Mark Lake for an efficient implementation of the model. He also thanks Guillaume Royer, Julien Pantz, Elias Daboussi and an anonymous referee for constructive feedback.
Email: dominique.bang@baml.com.
References
- Andersen L and V Piterbarg, 2010
Interest Rate Modeling, Volume 2
Atlantic Financial Press - Andreasen J and B Huge, 2013
Expanded forward volatility
Risk January, 101–107 - Antonov A, M Konikov and M Spector, 2017
Mixing SABR models for negative rates
Risk September, 86–91 - Balland P and Q Tran, 2013
SABR goes normal
Risk May, 76–81 - Bang D, 2018
Local-stochastic volatility models for vanilla modeling: a tractable and arbitrage free approach to option pricing
SSRN preprint, available at https://papers.ssrn.com/sol3/papers.cfm?abstract_id=3171877 - Berestycki H, J Busca and I Florent, 2004
Computing the implied volatility in stochastic volatility models
Communications in Pure and Applied Mathematics 57, pages 1352–1373 - Blacher G, 2001
A new approach for designing and calibrating stochastic volatility models for optimal delta vega hedging of exotic options
Working Paper, Global Derivatives Conference, Juan-les-Pins, France - Hagan P, D Kumar, A Lesniewski and D Woodward, 2002
Managing smile risk
Wilmott January, 84–108 - Lipton A, 2002
The vol smile problem
Risk February, 81–85 - Schlenkrich S, 2018
Approximate local volatility model for vanilla rates options
SSRN preprint, available at https://ssrn.com/abstract=3150689 - Ticot Y and X Charvet, 2013
Inflation derivatives – LPI swaps with a smile
Risk June, 66–71
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 Derivatives
One strike and they’re out: traders threaten liquidity stoppage
PTFs vow to withdraw from Treasuries market in protest at SEC registration plans
Collateral markets in need of rewiring
New data suggests a tech upgrade is needed to avoid a large central bank footprint in markets
Dutch pensions have extra year to restructure hedges
January 2028 implementation date allows more time for long-dated swaps to roll off
Fast LPs accuse rivals of maxing out last look response times
Firms with sub-10ms checks complain of losing volumes to slower rivals, prompting one to ditch ECNs
Swap Connect shines light on US client clearing hurdles
New scheme may intensify calls for CFTC to reassess its exempt DCO limitations
US life insurer index options market hits $1trn mark
Counterparty Radar: Lincoln Financial emerges as top player in Q4 with $43 billion portfolio increase
Long-end euro swap pricing anomaly remains largely untapped
Deviation in swap curve attracts limited interest because of regulatory and pension reform barriers
SG1 growth slower than expected, say LPs
Despite sluggish take-up of Singapore FX matching engines, some hope a new NDF venue will offer a boost