CALIBRATION TO SWAPTIONS IN THE LIBOR
MARKET MODEL
PIERRE BERET
NATIONAL UNIVERSITY OF SINGAPORE
2007 CALIBRATION TO SWAPTIONS IN THE LIBOR MARKET MODEL
PIERRE BERET (Ingenieur, Ecole Centrale Paris)
A THESIS SUBMITTED FOR THE DEGREE OF MASTER OF SCIENCE DEPARTMENT OF MATHEMATICS NATIONAL UNIVERSITY OF SINGAPORE 2007 i
Name : Pierre Beret Degree : Master of Science Supervisor : Dr Oliver Chen Department : Department of Mathematics Thesis Title : Calibration to swaptions in the Libor Market Model
Abstract
In this dissertation, the Libor Market Model is presented and its calibration process is derived. We assume the Forward Libor Rates follow log-normal stochastic processes with a d-dimensional Brownian motion and build an in- terest rates model able to price interest rate derivatives. We emphasize how dierent it is from the usual short-term interest rates models (Hull-White). Nevertheless, this pricing model only makes sense if vanilla products, namely caps and European swaptions, can be well priced with respect to their market value. To check this, we propose dierent parametric forms of instantaneous volatilities σi(t) and correlations ρij to obtain the best results. Then, we show a method to reduce the dimensionality of the Libor Market model compared to the number of Forward rates involved by using Rebonato An- gles and Frobenius norm. Finally, we derive approximations formula for European swaptions and show we can avoid Monte-Carlo simulations for the calculations of the swaptions during the calibration. Some numerical results are given on a 3 factors model. We discuss then dierent issues raised and current developments, more specif- ically the SABR skew form and cross-asset products.
Keywords : Interest Rate Derivatives, Libor Market Model, Cali- bration, Rank reduction methods, Swaption Approxi- mations. ii
Acknowledgment
I consider myself extremely fortunate to have been given the opportunity and privilege of doing this research work at the National University of Singapore. I would like to thank all the people who have helped me during my Master's degree program. All my gratitude to Doctor Oliver Chen who accepted to be my supervisor and provided warm and constant guidance throughout progress of this work. My warmest thanks to the Royal Bank of Scotland who welcomed me in its Exotic Rates Structuring Team for 6 months. This experience was very rich and I learned a lot with Serge Pomonti. I am happy to continue this collaboration in January. I am also thankful for the graduate research scholarship oered to me by the National University of Singapore without which this Master's degree program would not have been possible. Finally, I would like to express my deep aection for my family and my friends in Singapore who have encouraged me throughout this work and for Camille who supported me everyday.
March 3, 2007 Contents
1 Interest Rates Models 1 1.1 Important concepts ...... 2 1.1.1 Zero coupon bonds ...... 2 1.1.2 Short-Term interest rate ...... 3 1.1.3 The Arbitrage free assumption ...... 3 1.1.4 Forward Interest rates ...... 5 1.1.5 LIBOR interest rate and swaps ...... 6 1.1.6 Stochastic tools ...... 9 1.2 Interest Rates Models ...... 12 1.2.1 Short term interest rates ...... 13 1.2.2 Heath Jarrow and Morton Framework ...... 15 1.2.3 The Libor Market Model ...... 17 1.2.4 Libor Market model summary ...... 24 1.3 Pricing Vanilla Derivatives ...... 24 1.3.1 Interest rate options: cap and oor ...... 24 1.3.2 Swaptions ...... 29
2 Calibration of the Libor Market Model 32 2.1 The settings: Main purpose of the Calibration ...... 32 2.2 Structure of the instantaneous volatility ...... 34 2.2.1 Total parameterized volatility structure ...... 34 2.2.2 General Piecewise-Constant Parameterization ..... 35 iv CONTENTS
2.2.3 Laguerre function linear combination type volatility . 36 2.3 Structure of the correlation among the Forward Rates .... 40 2.3.1 Historic correlation vs parametric correlation ..... 41 2.3.2 Rank Reduction methods ...... 50 2.4 Swaption Approximation formulas ...... 62 2.4.1 Rebonato Formula ...... 62 2.4.2 Hull and White Formula ...... 64 2.4.3 Andersen and Andereasen Formula ...... 65 2.5 Monte Carlo Simulation and Results on 3 Factors BGM ... 66 2.5.1 Monte Carlo Method ...... 66 2.5.2 Numerical Results ...... 67
3 Perspectives and issues 71 3.1 Stochastic volatility models applied to Libor Market Model . 71
3.1.1 Stochastic α β ρ model - SABR ...... 71 3.2 Hybrids Products ...... 74 3.3 Issues raised ...... 76 3.3.1 Choice between Historical and Implied volatility ... 76 3.3.2 Interest-rates skew ...... 76 3.3.3 Approximation formula ...... 77 3.3.4 Market liquidity ...... 77
4 General Methodology proposed for calibration 78 4.1 Assumptions ...... 78 4.2 Modeling choices ...... 78 4.3 Market data ...... 79 4.4 Calibration process ...... 79 4.5 Conclusion ...... 80 List of Figures
1.1 Zero-coupon bond mechanism ...... 2 1.2 Swap mechanism ...... 7
2.1 Laguerre-type volatility ...... 39 2.2 Historical correlation among Forward rates ...... 42 2.3 Simple Exponential Parameterized correlation ...... 45 2.4 Modied Exponential Parameterized correlation ...... 47 2.5 Schoenmakers Coey correlation ...... 49 2.6 Eigenvectors comparison ...... 57 2.7 2Y Forward Libor Rate Correlation ...... 59 2.8 5Y Forward Libor Rate Correlation ...... 60 2.9 10Y Forward Libor Rate Correlation ...... 61 List of Tables
2.1 General volatility structure ...... 35 2.2 Piecewise-constant volatility structure ...... 36 2.3 Laguerre type volatility structure ...... 36 2.4 Eigenvalues of the correlation matrix ...... 55 2.5 Swaption approximation accuracy ...... 69 Chapter 1
Interest Rates Models
At the end of the 70's, after Black and Scholes breakthrough with their formula to value a European option, Black also proposed the alter ego of this formula in the world of interest rates. This was the beginning of the interest rates derivatives. Since 1976 and Black's formula [2], a lot has been proposed on the interest rates topic. First were presented models that tried to adapt the frameworks coming from the equity world : those used a stochastic equation to describe a short-term rate as it was done for a stock. From this basic idea dierent evolutions rose by changing the form of this stochastic dierential equation to t the economic behavior of the interest rates generally observed - for instance the mean reversion phenomenon. Finally in 1997, Brace, Gatarek and Musiela proposed a new concept where observable rates were modeled using the work of Heath, Jarrow and Morton in 1992. This completely redened the vision of pricing and everything needs to be done in this eld. The purpose of this model is undoubtedly to be able to t the market. Hence, we call calibration the choice of the dierent assumptions and inputs so that we obtain the best t to the market. Calibration is always a huge issue for market operators as they may face severe misprices if the model they use is not well calibrated and I will be 2 Interest Rates Models
presenting how this can be handled in the second part; before explaining what are the main issues and how some are managed (skew/smile, liquidity..) and what are the next challenges faced by the Libor Market Model (Cross-asset hybrid products). In this rst chapter the main denitions and the models currently used in the world of interest rates are dened and explained.
1.1 Important concepts
1.1.1 Zero coupon bonds
The rst concept we have to dene when discussing interest rates products is the Zero coupon bond (Z.C.). In this thesis, the underlying assets are not stocks like in Black-Scholes original framework in 1973 in [1] but bonds. Several bonds can be dened, paying various coupons, depending on some conditions...1Hence, it is necessary to dene a simplest underlying: this one is the set of discount factors for dierent maturities. We will denote them by B(t, T ). This bond represents at time t the price of 1 paid at time T , the maturity of the bond. See Figure 1.1 for a more visual explanation.
Figure 1.1: Zero-coupon bond mechanism
1For instance, a daily range accrual coupon: I pay n where is the number of X% N n days 3-months LIBOR rate stays below 6.5% and N the number of days in the accrual period. 3 Interest Rates Models
One can observe that at any date t, those prices are not all quoted on the market but can be obtained from other zero coupons bonds. This bond does not pay any coupon, that is why we generally call the discount factors
B(t, T ) the Zero coupon bonds (Z.C.). We introduce very generally the log-normal dynamic for a Zero Coupon bond as:
B dB(t, T ) = m(t, T )tB(t, T )dt + σ B(t, T )dWt,B(T,T ) = 1 (1.1)
With m(t, T ), the drift, equal to the short term interest rate rt in a risk- neutral world, σB, the volatility eventually stochastic or time-dependent and
Wt a Brownian motion.
1.1.2 Short-Term interest rate
We just mentioned the short term interest rate in the previous section. Tra- ditional stochastic interest rates models are based on the exogenous speci- cation of a short-term interest rate and its dynamic. We will denote by rt the instantaneous interest rate or short-term interest rate the rate one can borrow in a risk free loan beginning at t over the innitesimal period dt.
In general, we assume that rt is an adapted process on a ltered proba- bility space. The important thing about short term interest rate is that by consideration over the absence of arbitrage in the market we can create links between rt and B(t, T ).
1.1.3 The Arbitrage free assumption
This classic assumption introduces constraints on the payo of derivatives. Here when we study rate issues, this assumption is made on the Zero coupon bonds as we can link long maturities (more than 1 year) bonds with coupons with Zero coupon bonds by considering the Arbitrage free assumption. 4 Interest Rates Models
The price of an asset delivering xed cash-ows in the future is given by the sum of its cash-ows weighted by the price of the Zero coupon bonds of the settlement dates.
We make the usual mathematical assumption: all processes are dened on a probability space (Ω, {Ft; t ≥ 0}, Q0). The probability measure Q0 is any risk neutral probability measure whose existence is given by the no-arbitrage as- sumption (See The Girsanov transformation in section 1.1.6). The ltration
2 {Ft; t ≥ 0} is the ltration generated in Q0 by a d-dimensional Brownian motion W Q0 = {W Q0 (t); t ≥ 0}.
Now, we infer that one can invest in a savings account continuously compounded with the stochastic short rate rs prevailing at time s over the time [s; s + ds]. The value of 1 invested at time t at time T is βT :
Z T βT = exp rsds t
Therefore, if we invest B(t, T ) in a Z.C. of maturity T and the same amount in our saving account, the fundamental theorem of asset pricing (this will be detailled in 1.1.6) ensures that they produce on average over all the paths the same amount namely 1. This equality at time t can be written:
Z T Q0 B(t, T ) = Et exp −rsds |Ft t
In the case of a deterministic rate rs, as B(T,T ) = 1:
Z T B(t, T ) = exp −rsds t
2 In a nancial point of view, the ltration{Ft; t ≥ 0} represents the structure of all the information known by every market agent. 5 Interest Rates Models
And in the case of a constant deterministic rate r compound n-times per year: 1 B(t, T ) = (1.2) r (T −t) (1 + n )
1.1.4 Forward Interest rates
We can dene Forward Interest Rates for all the previous rates we saw:
Bt(T,T + δ) is the forward value at t of a Z.C. invested at T which will pay 1 at T + δ. By arbitrage we know it is worth:
B(t, T + δ) B (T,T + δ) = t B(t, T )
The equivalent rate simply compounded to this Zero Coupon Bond can be computed writing:
1 B(t, T ) F (t, T ) = − 1 (1.3) δ δ B(t, T + δ)
This rate is named the Forward Rate and is the constant rate simply com- pounded to be paid if you want to borrow money at time t for a future time period between T and T + δ. We can also dene f(t, T ) the instantaneous forward interest rate, the for- ward version of rt. Formally, f(t, T ) is the forward rate at t one can borrow in a risk free loan beginning at T over the innitesimal period dt. This con- cept is rather a mathematical idealization as it can not be observed in the market but is useful to describe bond price models. One can write:
Z T B(t, T ) = exp − f(t, u)du , ∀t ∈ [0,T ] (1.4) t 6 Interest Rates Models
1.1.5 LIBOR interest rate and swaps
Libor interest rate
During the 80's, Libor (which stands for London Inter Bank Oered Rates) interest rates have become more and more traded. This rate is declined for dierent short maturities (inferior to one year) and is a benchmark of the main banks of their loan rate for those maturities. It is xed everyday at 11h00 am, London Time. It is considered in general as the risk-free interest rate by the investors: even credit default swaps values are given with respect to the LIBOR curve. However, this is not true, those nancial institutions have a probability of default and hence this default risk is quantied. In the markets, the risk free does not really exist but it can be assumed that the main central banks (More specically: US Fed, ECB, CBE) have an almost nil probability of default as they can literally print their money and hence the bonds they issue called treasuries have almost no probability of default3. The spread between the LIBOR and the treasury rate represents this risk to default. For the USD Market, LIBOR rates trade around 50 basis points above treasury rates.
We call Lδ(t, t), the LIBOR Interest rate at time t for a maturity of δ:
1 = B(t, t + δ) (1.5) 1 + δLδ(t, t) with δ is three or six months usually. Using the arbitrage free rule and applying the previous section about Forward
Interest rates to Libor Interest Rates and their Forwards Lδ(t, T ) the Libor rate at time t at which one can borrow money at time T for a maturity of δ we can write: 1 B(t, T + δ) = 1 + δLδ(t, T ) B(t, T ) 3It should be emphasized that the sovereign risk is real: in July 1998, Russia defaulted on its bonds causing the fall of the famous hedge-fund LTCM. 7 Interest Rates Models
That is, B(t, T ) − B(t, T + δ) L (t, T ) = (1.6) δ δB(t, T + δ)
We will skip the index δ when there will be no ambiguities about the matu- rity.
Swap rate
The rst swap contracts were also negotiated in the early 1980s. Since, it has shown an amazing growth becoming more and more important in the exotic derivatives market. A swap is a contract between two companies to exchange a predened cash ow in the future. The schedule of the cash ows and the way they are calculated is specied in this agreement. At the beginning, swaps were tailored for companies who wanted to hedge their loans exposure and lock in a good level of interest rate. Hence one can decide to enter a swap where he will exchange his semi-annual
xed rates cash-ows at x% against a oating rate, for instance the value of the 6-months LIBOR rate with xing date at the beginning of the 6-months period (Fixing in advance 4) The following Figure 1.2 explains how is built the exchange of cash-ows from the customer point of view. This type of
Figure 1.2: Exchange of cash-ows for a Payer Swap
4Several issues are not mentioned here about the xing dates and the convexity ad- justment that are necessary when pricing non perfectly scheduled structure or in arrears xing structures, for instance see [3] 8 Interest Rates Models
swap is called is called a payer swap. The symmetric version is called receiver swap.
As a matter of fact, from this denition appears the swap rate Sp,n(t)dened as the rate which gives a net present value of 0 at time t to the swap which exchange this swap rate against a oating one (δ-months Libor Lδ(t, Ti)) on a schedule Ti, i = p, . . . , n. We can compute this swap rate Sp,n(t) by arbitrage considerations and, it is worth noticing it, independently of any model assumption. The xed leg is:
n−p X F ixedp,n(t) = Sp,n(t)δB(t, Tp+i) i=0
And the oating leg is:
n−p X F loatingp,n(t) = B(t, Ti+p)δL(t, Ti−1+p) i=1 n−p X B(t, Ti−1+p) = B(t, T ) − 1 i+p B(t, T ) i=1 i+p n−p X = B(t, Ti−1+p) − B(t, Ti+p) i=1
= B(t, Tp) − B(t, Tn)
The swap rate is by denition the one that equalize both legs:
F ixedp,n(t) = F loatingp,n(t) B(t, T ) − B(t, T ) S (t) = p n p,n Pn−p i=0 δB(t, Tp+i)
This swap was more precisely a forward start interest rate swap which rst settlement date is Tp. Once this product was well understood by every one on the markets, it naturally gave rise to its rst most natural derivative: 9 Interest Rates Models
the European swaption 5. A European swaption is a one-time option on a swap rate. From now, we will always refer to European swaptions when we describe swaptions. When one is long a swaption strike Sp,n, he owns the right and not the obligation to enter a swap of tenor Tn at maturity Tp. A swaption can be computed through dierent methods but the market in general quotes the implied volatility of the swaption with the generalization of the Black formula (See section 1.3.2). On the mathematical side this arise issues as one can show that swap rates and forward rates can not be log normal at the same time. We will discuss later this point in section 2.4.
1.1.6 Stochastic tools
This subsection is going to present a few stochastic tools we need to describe the basics of the Libor Market Model. This subsection does not seek to be exhaustive and totally rigorous in stochastic calculus but just to give a general idea about the tools we will be using in the construction of the models in the next section. For further details about stochastic calculus please refer to the excellent [5].
Numeraire
A Numeraire is a price process (A(t))T (a process is a sequence of random variables), which is strictly positive for all t ∈ [O,T ]. Numeraires are used to express prices in order to have relative prices. The application of this rather abstract concept can be seen in what follows.
Change of numeraire
6 Let P and Q be equivalent measures with respect to the numeraires A(T ) and B(t). The Radon-Nikodym derivative that changes the equivalent mea-
5American and Bermudean swaption also exist but are not as liquid and as vanilla than European 6 P and Q are equivalent if and only if : P(M) = 0 ↔ Q(M) = 0, ∀M ∈ F 10 Interest Rates Models
sure P in Q is given by:
d A(T )B(t) R = P = (1.7) dQ A(t)B(T )
This derivative is very useful: due to the no arbitrage rule the price of an asset X should be independent from the choice of the measure and numeraire:
X(T ) X(T ) A(t) P |F = B(t) Q |F E A(T ) t E B(T ) t
If one introduces: X(T ) and doing some simple manipulation on the G(T ) = A(T ) previous equation:
A(T )B(t) P (G(T )|F ) = Q G(T ) |F E t E A(t)B(T ) t
Q = E (G(T )R|Ft)
We can see that we can change the probability measure just by multiplying the martingale by its Radon-Nikodym derivative.
Girsanov theorem
For any adapted stochastic process k(t) which satises the following con- dition: