Derivatives and Volatility · Derivatives
11SABR and Smile Dynamics
In 2002 four authors working on an interest-rate desk published a short paper with an uncomfortable finding. The model many desks used to fit swaption smiles, local volatility, predicted that when the forward rate rose the smile would move the other way, towards lower strikes; the market’s smile moved with the forward. A model that gets the direction wrong gets the delta wrong, and the authors showed that its hedges could be worse than hedging with no smile at all. Their replacement had four parameters, a closed-form approximation for the implied volatility that a spreadsheet could evaluate, and the right dynamics. SABR became the market standard for interest-rate options and a common one for currency options. This chapter states the model and its formula, explains the backbone that governs its dynamics, derives the delta that its dynamics imply, measures by simulation what that delta is worth, and shows where the formula breaks.
11.1 The SABR model
Definition 11.1 (SABR model)
The SABR model (stochastic alpha, beta, rho) describes a forward , driftless under its forward measure, and its volatility :
with an exponent , a correlation and a volatility of volatility ; the initial value is the fourth parameter.
It is a stochastic volatility model in the sense of chapter 10, with two differences from Heston: the volatility is lognormal and does not mean-revert, which suits a model calibrated one expiry at a time, and the forward has a constant-elasticity diffusion , which interpolates between a normal model (, the Bachelier model of One Quant Book 2, chapter 13) and a lognormal one (). Interest-rate desks calibrate a separate SABR to each expiry and tenor of the volatility cube; the normal and shifted versions they use when rates approach zero are One Quant Book 6, chapter 5.
11.2 The asymptotic formula
SABR has no closed-form price. Its authors derived, by singular perturbation in the time to expiry and the log-moneyness, an explicit approximation of the Black implied volatility.
Definition 11.2 (Hagan formula)
The Hagan formula for the SABR Black volatility at strike and expiry is
with , and .
At the money and the formula reduces to , a cubic in . Desks therefore parametrise SABR by the at-the-money volatility they observe, solve the cubic for , fix by convention or from the backbone (next section), and fit and to the smile. The normal (Bachelier) volatility the rates market quotes is the same price read with the normal model.
Example 11.3 (A three-percent forward)
A one-year option on a forward rate of 3%, at-the-money Black volatility 20%, , , : the cubic gives . The Black smile runs from 28.8% at a 2% strike to 15.8% at 4%; the same prices read as normal volatilities run from 70.8 to 55.0 basis points a year, 59.9 at the money. The normal smile is far less skewed: much of the Black skew is the lognormal model’s own scaling.
Example 11.4 (SABR on an equity smile)
With (a lognormal backbone), SABR fitted to the one-year smile of chapter 9’s surface returns , , , with a root-mean-square error of 0.07 volatility point and at most 0.15 point, at the 70 strike (Figure 11.1). One expiry, three free parameters, a good fit: SABR is a smile interpolator as much as a model.
11.3 Backbone and smile dynamics
Definition 11.5 (Backbone)
The backbone of a smile model is the curve traced by the at-the-money implied volatility as the forward moves, all other state variables fixed. In SABR it is : flat for , falling like for .
The backbone is how SABR encodes dynamics. When the forward moves, the whole smile moves with it: the at-the-money point slides along the backbone and the smile’s shape travels with the new forward (Figure 11.2). That is the behaviour Hagan and his coauthors observed in rates markets and the opposite of local volatility (chapter 9), whose smile moves against the forward. Two parameters produce skew: gives a skew along the backbone, a skew relative to it. They are nearly interchangeable for fitting one smile, and distinguishable only by dynamics: is chosen from how at-the-money volatility has moved with the forward historically, or by convention.
11.4 Hedging under a smile model
A smile model’s delta is the change of the option’s value when the underlying moves, with the model’s own view of what the other state variables do meanwhile. Three answers are in use.
- The Black delta freezes the option’s implied volatility: .
- Hagan’s delta moves the volatility along the smile as moves, with fixed: .
- The third takes into account that itself tends to move when moves, since the two Brownian motions are correlated.
Definition 11.6 (Minimum-variance delta)
The minimum-variance delta of an option in a stochastic volatility model is the position in the underlying that minimises the variance of the hedged position’s instantaneous P&L. In SABR (Bartlett’s delta) it is
Proposition 11.7 (Why that delta)
The option’s value changes by to first order. The regression of on has slope , so the position removes the part of correlated with and leaves only the uncorrelated part, which no position in the underlying can hedge.
Proof. and , so . Substituting, the P&L of is ; its variance is minimal when the first bracket vanishes. With , and . ∎
On the example of the chapter, the one-year at-the-money call has a Black delta of 0.540, a Hagan delta of 0.580 and a minimum-variance delta of 0.461 (Figure 11.3). With a negative correlation, a rise in the forward tends to bring a fall in and so in the call’s value: the right hedge is smaller than Black’s. Hagan’s delta moves the other way because it ignores the correlation and sees only the backbone.
How much does it matter? The tutorial simulates one day of SABR dynamics on 200 000 paths and reprices the call with the formula at the new forward and . The standard deviation of the hedged one-day P&L is with the Black delta, with Hagan’s and with the minimum-variance delta: 19% less variance than Black, and 25% more with Hagan’s. What remains is the part of the volatility’s move that is independent of the forward, hedged only with vega (another option).
11.5 Limits: wings and long expiries
Hagan’s formula is an expansion in and in . It is accurate for the expiries and strikes of the liquid swaption market, and it fails in two known ways. At long expiries with a large volatility of volatility it produces call prices that are not convex at low strikes: a negative density (Figure 11.4). And its zero-order term was shown to be slightly wrong (Obłój’s correction), which matters in the wings. The fixes used in practice are to solve an effective one-dimensional forward equation for the density with the same asymptotic accuracy (“arbitrage-free SABR”), or to shift the forward (One Quant Book 6, chapter 5), or to price the wings with a different extrapolation.
11.6 Tutorial: SABR from formula to hedge
Goal. Implement Hagan’s formula, calibrate SABR with tied to the at-the-money volatility, compare three deltas, and measure their hedging errors by simulation. End state: the four figures of the chapter and the numbers of the weekend problem.
The formula, with the at-the-money limit of handled explicitly:
def hagan_lognormal(f: float, k: float, t: float, alpha: float, beta: float, rho: float, nu: float) -> float: """Black implied volatility of the SABR model (Hagan et al. 2002, 2.17).""" fk = f * k lfk = math.log(f / k) omb = 1.0 - beta pre = alpha / (fk ** (omb / 2) * (1 + omb ** 2 / 24 * lfk ** 2 + omb ** 4 / 1920 * lfk ** 4)) z = nu / alpha * fk ** (omb / 2) * lfk if abs(z) < 1e-8: zx = 1.0 - 0.5 * rho * z else: x = math.log((math.sqrt(1 - 2 * rho * z + z * z) + z - rho) / (1 - rho)) zx = z / x corr = 1 + (omb ** 2 / 24 * alpha ** 2 / fk ** omb + 0.25 * rho * beta * nu * alpha / fk ** (omb / 2) + (2 - 3 * rho * rho) / 24 * nu * nu) * t return pre * zx * corrListing 11.1. Hagan’s Black-volatility formula. code/firm/sabr/firm_sabr.py The three deltas, the correlation correction coming from :
def deltas(f: float, k: float, t: float, alpha: float, beta: float, rho: float, nu: float, h: float = 1e-6) -> dict: """Black delta (volatility frozen), Hagan's delta (volatility moves along the smile with F, alpha fixed) and Bartlett's delta (alpha also moves with its expected co-move rho nu / F^beta dF).""" def vol(ff: float, aa: float) -> float: return hagan_lognormal(ff, k, t, aa, beta, rho, nu) sig = vol(f, alpha) g = greeks(f, k, t, 0.0, 0.0, sig, "C") # on the forward, zero rates dsig_df = (vol(f + h, alpha) - vol(f - h, alpha)) / (2 * h) ha = 1e-6 * alpha dsig_da = (vol(f, alpha + ha) - vol(f, alpha - ha)) / (2 * ha) black_d = g["delta"] hagan_d = black_d + g["vega"] * dsig_df bartlett_d = hagan_d + g["vega"] * dsig_da * rho * nu / f ** beta return {"black": black_d, "hagan": hagan_d, "bartlett": bartlett_d, "vega": g["vega"], "sigma": sig}Listing 11.2. Black, Hagan and minimum-variance deltas. code/firm/sabr/firm_sabr.py - Run
dv_sabr.hedge_experiment(),equity_fit(),wing_density()andfig_sabr.py.
What to change next. Set and compare the three deltas’ hedging errors (the minimum-variance delta then nearly equals Black’s); then fit the equity smile with and compare and with the fit.
11.7 Build: the SABR smile
Purpose. The miniature firm’s smile for rates and currency options, and the parametrisation that One Quant Book 6 extends to the volatility cube.
Interface. hagan_lognormal(f, k, t, alpha, beta, rho, nu); alpha_from_atm(f, t, sigma_atm, beta, rho, nu); calibrate(f, t, strikes, vols, beta, atm_vol) -> ((alpha, rho, nu), rmse); normal_vol; deltas(…) -> black, hagan, bartlett; one_day_hedge_errors.
Rules. is an input, not a fitted parameter; is always solved from the at-the-money volatility (the smallest positive root of the cubic); report the density of the wings, and do not use the formula where it is negative.
Acceptance tests. code/firm/sabr/tests/: at and the smile is flat at ; the at-the-money volatility is reproduced exactly; the fit recovers known parameters; the minimum-variance delta has the lowest simulated hedging error for negative correlations.
Stretch. The arbitrage-free density of Hagan and his coauthors (2014), solved on a grid, and Obłój’s corrected leading term.
Sources and further reading
- P. S. Hagan, D. Kumar, A. S. Lesniewski and D. E. Woodward, “Managing smile risk”, Wilmott (2002).
- B. Bartlett, “Hedging under SABR model”, Wilmott (July/August 2006) 2–4.
- J. Obłój, “Fine-tune your smile: correction to Hagan et al.”, Wilmott (May 2008) 102–109.
- P. S. Hagan, D. Kumar, A. S. Lesniewski and D. E. Woodward, “Arbitrage-free SABR”, Wilmott 69 (2014) 60–75.
11.8 Exercises
Exercise 11.1 ★
With and , what is the SABR smile? With and ?
Solution
Solution of Exercise 11.1.
, : a lognormal forward with constant volatility , a flat Black smile at . , : the Bachelier model, a flat normal smile at , which in Black volatility is a skew falling roughly like .
Exercise 11.2 ★
Using the leading term with set for 20% at a 3% forward, give the at-the-money volatility at a 2.5% forward for , 0.5 and 1.
Solution
Solution of Exercise 11.2.
and : 24.0% for , 21.9% for , 20.0% for .
Exercise 11.3 ★
Why does SABR not need mean reversion when it is calibrated expiry by expiry?
Solution
Solution of Exercise 11.3.
Each expiry is fitted with its own parameters, so the model never has to describe how volatility behaves across expiries: mean reversion, which shapes a term structure, has nothing to do.
Exercise 11.4 ★★
Derive the correction term of the minimum-variance delta from the regression of on .
Solution
Solution of Exercise 11.4.
and with correlation : and , so the slope is . The option’s sensitivity to times that slope is the part of its -risk hedgeable with the forward: .
Exercise 11.5 ★★
and both produce skew. How would you choose for a currency pair, and what does the wrong choice cost?
Solution
Solution of Exercise 11.5.
From the historical co-movement of at-the-money volatility and the forward: regress the change of on the change of ; the slope estimates . A wrong is compensated by in the fit, so prices are unchanged, but the deltas are wrong by the difference in backbone slope times vega.
Exercise 11.6 ★★
Give the normal volatility of the 2%, 3% and 4% strikes in Example 11.3 and explain why the normal smile is flatter than the Black smile.
Solution
Solution of Exercise 11.6.
70.8, 59.9 and 55.0 basis points a year, against 28.8%, 20.0% and 15.8% in Black volatility. A normal volatility is roughly the Black volatility times ; much of the fall of the Black smile with the strike is the lognormal model’s own scaling, which the normal quote removes.
Exercise 11.7 ★★★
Coding. Repeat the one-day hedging experiment with , : which delta is best, and by how much?
Solution
Solution of Exercise 11.7.
The minimum-variance and Black deltas give almost the same error (variance ratio 0.999); Hagan’s is 10% worse. With a positive correlation the backbone and the volatility’s co-move cancel nearly exactly at the money.
Exercise 11.8 ★★★
Find the flaw. “Our SABR fit is perfect at every expiry, so its ten-year 1% strike price is reliable.”
Solution
Solution of Exercise 11.8.
The fit is to quoted strikes and expiries; a ten-year 1% strike is outside them, where Hagan’s formula can produce a negative density (below about 1.6% in the chapter’s ten-year example), so its price there is not even arbitrage-free. Use the arbitrage-free density or a shifted model, and a reserve for the extrapolation.
11.9 Problem: Which Delta
Problem 11.1
Weekend problem — the delta a smile model should use
A rates desk is long one-year at-the-money calls on a 3% forward (payer swaptions in normalised units), priced with SABR: , , , 20% at-the-money Black volatility, zero discounting. It must choose the delta it hedges with.
Part I — The model.
- Solve for .
- Give the Black volatility at the 2%, 3% and 4% strikes.
- Give the at-the-money normal volatility in basis points.
- What does the backbone predict for the at-the-money volatility if the forward rises to 3.5%?
- How does the smile move when the forward rises?
Part II — Three deltas.
- Give the Black, Hagan and minimum-variance deltas of the at-the-money call.
- Why is Hagan’s delta larger than Black’s?
- Why is the minimum-variance delta smaller than both?
- Give the regression slope of on .
- What risk remains after the minimum-variance hedge, and how is it hedged?
Part III — The experiment.
- Give the standard deviations of the one-day hedged P&L with each delta.
- Give the variance of each relative to Black’s.
- Why does Hagan’s delta do worse than Black’s here?
- What would the result be with ?
- If daily hedging errors are independent, by how much does a 19% lower daily variance reduce the standard deviation of a year’s hedging P&L?
Part IV — Judgement.
- Which delta would you report to risk, and which would you trade?
- How does the choice of change the deltas?
- What if the market’s smile dynamics differ from SABR’s?
- State the named result: the reduction in the variance of the one-day hedging error from the minimum-variance delta relative to the Black delta.
- In one sentence: what does a smile model’s delta depend on that its prices do not?
Solution
Solution of Problem 11.1.
1. . 2. 28.8%, 20.0%, 15.8%. 3. 59.9 basis points a year. 4. 18.5%. 5. With the forward: its minimum and shape travel to the right while its level slides down the backbone. 6. 0.540, 0.580 and 0.461. 7. Hagan’s delta adds vega times the slope of the implied volatility in , positive here because the at-the-money point moves up the smile’s right side as rises (the volatility of a fixed strike rises with ). 8. With a rise in comes with a fall in , which lowers the call’s value; the hedge should anticipate it. 9. . 10. The part of uncorrelated with : vega risk, hedged with another option. 11. (Black), (Hagan), (minimum-variance). 12. 1.25 and 0.81 of Black’s. 13. It moves the volatility along the smile with fixed, while in the model tends to fall when rises: it corrects in the wrong direction. 14. Hagan’s and the minimum-variance delta coincide, 1% less variance than Black’s in this example. 15. By , about 10%. 16. Report the model’s delta used for hedging, and trade the minimum-variance delta if the model’s dynamics are believed; show the others as sensitivities to the dynamics. 17. A lower steepens the backbone, raising Hagan’s delta for a call; refits the smile and changes the correction term; the deltas move by vega times the change in slope. 18. The minimum-variance delta is only as good as the model’s co-movement of and ; estimate the co-movement from data and compare. 19. 19% less variance than with the Black delta (ratio 0.81). 20. On the model’s dynamics: how volatility moves with the underlying.
11.10 Interview questions
Interview question 11.1 ★ trader, researcher
Write the SABR model. What does each parameter do?
Interview question 11.2 ★ trader
What is the backbone, and how would you estimate from market data?
Solution
Solution of Interview question 11.2.
The at-the-money volatility as a function of the forward, . Regress changes of on changes of over a period of stable regime; the slope is .
What the interviewer is looking for: the definition and an estimation procedure.
Interview question 11.3 ★★ researcher
Derive the minimum-variance delta in a stochastic volatility model.
Solution
Solution of Interview question 11.3.
; ; regress on : slope ; hedge slope. The remaining term is uncorrelated with .
What the interviewer is looking for: the regression argument.
Interview question 11.4 ★★ researcher, risk
Why did local volatility give bad hedges in the rates market, according to the authors of SABR?
Solution
Solution of Interview question 11.4.
Local volatility predicts that the smile moves opposite to the forward, while the rates market’s smile moved with it; the model’s delta therefore had the wrong correction and its hedges could be worse than Black’s.
What the interviewer is looking for: dynamics in the wrong direction.
Interview question 11.5 ★★ developer
Your SABR implementation returns negative butterfly prices for long-dated low strikes. Why, and what do you do?
Solution
Solution of Interview question 11.5.
Hagan’s formula is an asymptotic expansion; far from the money and at long expiries it implies non-convex call prices. Detect it by computing the density, and price there with the arbitrage-free SABR density (a one-dimensional forward equation), a shifted model, or a separate wing extrapolation.
What the interviewer is looking for: the expansion’s domain and a fix.
Interview question 11.6 ★★★ trader, researcher
and fit the same smile equally well. Which trades distinguish them, and how would you hedge a book in which they are uncertain?
Solution
Solution of Interview question 11.6.
Trades whose value depends on how the smile moves with the forward: delta-hedged options (through the delta), forward-start options and spread options between tenors. Hedge by computing risks under several pairs fitted to the same smile and holding a reserve for the spread of their P&Ls, or by hedging vega by strike (bucketed) so that the model’s dynamics matter less.
What the interviewer is looking for: dynamics-sensitive products and model-uncertainty management.