Barrier Option Pricing under SABR Model Using Monte Carlo Methods
by
Junling Hu
A Project Report
Submitted to the Faculty
of the
WORCESTER POLYTECHNIC INSTITUTE
In partial fulfillment of the requirements for the
Degree of Master of Science
in
Financial Mathematics
May 2013
APPROVED:
Professor Stephan Sturm
Professor Bogdan Vernescu, Head of Department ACKNOWLEDGMENTS
This project could not have been possible without the support of my professors and my friends. In particular, I would like to express my gratitude towards Prof. Stephan Sturm for his guidance throughout the study.
1 Abstract
The project investigates the prices of barrier options from the constant underly- ing volatility in the Black-Scholes model to stochastic volatility model in SABR framework. The constant volatility assumption in derivative pricing is not able to capture the dynamics of volatility. In order to resolve the shortcomings of the Black-Scholes model, it becomes necessary to find a model that reproduces the smile effect of the volatility. To model the volatility more accurately, we look into the recently developed SABR model which is widely used by practitioners in the financial industry.
Pricing a barrier option whose payoff to be path dependent intrigued us to find a proper numerical method to approximate its price. We discuss the basic sam- pling methods of Monte Carlo and several popular variance reduction techniques. Then, we apply Monte Carlo methods to simulate the price of the down-and-out put barrier options under the Black-Scholes model and the SABR model as well as compare the features of these two models. Contents
1 Introduction of Options 1 1.1 Definition of Options ...... 1 1.2 The Development of Option Trading and Option Pricing ...... 2 1.3 Summary of Vanilla Option Pricing under the Black Scholes .... 3
2 Barrier Option Pricing under the Black Scholes 7 2.1 Valuation Methods ...... 8 2.2 Analytical Solution under the Black-Scholes Model ...... 8 2.2.1 A Continuity Correction to Discrete Barrier Option ..... 14 2.3 Monte Carlo Methods ...... 15 2.3.1 The Main Framework of Monte Carlo ...... 16 2.3.2 Variance Reduction Techniques ...... 18 2.3.3 Path Generation of Barrier Options ...... 29 2.3.4 Application of Variance Reduction to Barrier Options .... 31
3 Barrier Option Pricing under SABR Model 41 3.1 The Need for a Stochastic Volatility Model ...... 41 3.2 Building SABR Model Step by Step ...... 43 3.2.1 Black Model and Implied Volatility ...... 43 3.2.2 Local Volatility Models ...... 46 3.2.3 The SABR Model ...... 49 3.3 Estimating the SABR Parameters ...... 51 3.3.1 Choice of β ...... 53 3.3.2 Estimate α, ρ and ν directly ...... 54 3.3.3 α from σATM ...... 55 3.3.4 Algorithm for Calibration ...... 56 3.4 Dynamics of Parameters ...... 60 3.4.1 α ...... 60 3.4.2 β ...... 60 3.4.3 ρ ...... 61 3.4.4 ν ...... 62 3.4.5 f ...... 62 3.5 Monte Carlo Methods under SABR ...... 64 3.5.1 Path Generation of Barrier Options ...... 65
1 3.5.2 Result Discussion ...... 66
Appendices 69
A Mathematical Appendix 70 A.1 Probability Space ...... 70 A.2 Conditional Expectation ...... 72 A.3 Stochastic Processes ...... 72 A.3.1 Brownian Motions ...... 73 A.3.2 Martingales ...... 73 A.4 Convergence and Central Limit Theorem ...... 74 A.4.1 Convergence ...... 74 A.4.2 The Law of Large Numbers ...... 75 A.4.3 Central Limit Theorem ...... 75 A.5 Change of Measure ...... 75 A.5.1 Radon-Nikodym Derivative ...... 76 A.5.2 Girsanov’s Theorem ...... 76 A.5.3 Equivalent Martingale Measure ...... 77
2 List of Figures
2.1 Antithetic Path Example ...... 19 3.1 β =0.399 α =1.1649 ρ =0.1659 ν =1.2543 ...... 56 3.2 α from σATM approach ...... 57 3.3 α dynamics ...... 61 3.4 β dynamics ...... 61 3.5 ρ dynamics ...... 62 3.6 ν dynamics ...... 62 3.7 f dynamics ...... 63 3.8 One Simulation Path Example ...... 66 3.9 50 Simulation Paths Example ...... 67
3 Chapter 1
Introduction of Options
1.1 Definition of Options
Options are traded both on exchanges and in the over-the counter market. There are two types of options. A call option gives the holder the right to buy the under- lying asset by a certain date for a certain price. A put option gives the holder the right to sell the underlying asset by a certain date for certain price. The price in the contract is known as the exercise price or strike price; the date in the contract is known as the expiration date or maturity. [15]
Options can be either American or European, a distinction that has nothing to do with geographical location. American options can be exercised at any time up to the expiry date, whereas European options can be exercised only on the expiration date itself.
Options has been considered to be the most dynamic segment of the derivative markets since the inception of the Chicago Board Options Exchange (CBOE) in April 1973, with more than 1 million contracts per day, CBOE is the largest option exchange in the world. After that, several other option exchanges such as London
1 2
International Financial Futures and Options Exchange had been set up.
1.2 The Development of Option Trading and Option Pric-
ing
The history of stock options trading began with the 1973 establishment of the CBOE and the development of the Black-Scholes option pricing model. Over the last a few decades due to the famous work of Black and Scholes, the option val- uation problem has gained a lot of attention. In Black and Scholes 1973 seminal paper [4], the assumption of log-normality on stock price was applied. Moreover, in the same year, Robert Merton [19] extended the Black-Scholes model in several important ways. The application of the Black-Scholes-Merton option pricing mod- el for valuing various range of financial instruments and derivatives is considered essential.
Options form the foundation of innovative financial instruments, which are ex- tremely versatile securities that can be used in many different ways. Over the past decade, option has been utilized to speculation – leverage of pay-off on the underlying asset and hedging – reducing risk or providing portfolio insurance for financial institutions.
Following Black-Scholes, a number of other popular approaches was developed to pricing options with payoff depending on the value of the underlying asset at a single time horizon, including solving PDEs, binomial tree model, the nu- merical methods such as Monte-Carlo Simulation, Finite Differences Method and Replication. These work pushed option pricing into a very fascinating position in 3
quantitative finance.
The Black-Scholes formula is still around, even though it depends on several unrealistic assumptions. In special cases, though, we can improve the formula by making more realistic assumptions. However, we haven’t produced a formula that works better across a wide range of circumstances. One theory of the 1987 crash relies on incorrect beliefs, held before the crash, about the extent to which investors were using portfolio insurance, and about how changes in stock prices cause changes in expected returns.
The complexity of innovative derivatives came up with the complexity of option pricing analytical formulas. At the same time, the demand of speed in derivatives trading require fast ways to process these calculations. As a result, the develop- ment of computational methods for option pricing models can be a better solution.
1.3 Summary of Vanilla Option Pricing under the Black
Scholes
To price the option, we denote the value of the option C, on an underlying asset
St which pays a function f(ST ) at maturity time T. The interest rate, which is constant, to be r.
The payoff at maturity of an European call option with strike price K is defined by f(S ) = max(S K, 0) (1.1) T T − 4
The payoff at maturity of an European put option with strike price K is defined by f(S ) = max(K S , 0) (1.2) T − T
The Black-Scholes Assumption
In order to incorporate above factors as variables into one consistent model. Black and Scholes made some explicit assumptions on the market for a particular stock:
There is no arbitrage opportunity (no way to make riskless profits). • It is possible to borrow and lend cash at a known constant risk-free interest • rate r.
The stock price follows a geometric Brownian motion with constant drift • and volatility σ.
It is possible to buy and sell any amount, even fractional, of stock (this • includes short selling).
The above transactions do not incur any fees or costs (i.e., frictionless mar- • ket).
The underlying security does not pay a dividend. • Now we can define the process of the stock price to be a stochastic differential equations below:
dSt = uStdt + σStdWt (1.3)
Apply Itˆo’s lemma to function ln St where St is driven by the diffusion process above. Then ln St follows the SDE:
1 d ln S =(u σ2)dt + σdW (1.4) t − 2 t 5
Integrating from t to T , we have:
1 2 (u 2 σ )(T t)+σ(WT Wt)
ST = Ste − − − (1.5) É É denotes a new probability measure where there is a Brownian motion that
Q u r Wt = Wt + −σ t
Its dynamics under the É measure is:
Q dSt = rStdt + σStdWt (1.6)
˜ ˜ St Apply Itˆo’s lemma to St, the discounted stock price St = ert is driven by a simple SDE: ˜ ˜ Q dSt = σStdWt (1.7)
1 2 Q 2 σ t+σWt S˜t = S˜0e− (1.8)
Therefore, from the risk neutral valuation argument, the arbitrage-free price under
risk neutral measure É:
r(T t) Q
C (S , t)= e− − E [max (S K, 0) F ] (1.9) BS t T − | t 1 2 Q Q r(T t) Q (r 2 σ )(T t)+σ(WT Wt ) + C (S , t)= e− − E [S e − − − K] (1.10) BS t { t − }
r(T t) C (S , t)=Φ(d )S Φ(d )Ke− − (1.11) BS t 1 t − 2 2 St σ ln K +(r + 2 )(T t) d1 = − σ√T t − d = d σ√T t 2 1 − − (1.12) 6
Φ( ) is the cumulative distribution function of a standard normal distribution.
Similarly, the value of the put option is given by
r(T t) P = Φ( d )Ke− − Φ( d )S (1.13) BS − 2 − − 1 t
The Black-Scholes formula gives the price of European put and call options. This price is consistent with the Black-Scholes PDE: this follows since the formula can be obtained by solving the PDE for the corresponding terminal and boundary conditions.
∂V 1 ∂2V ∂V + σ2S2 + rS rV = 0 (1.14) ∂t 2 ∂S2 ∂S − Chapter 2
Barrier Option Pricing under the
Black Scholes
A barrier option is a type of exotic option, in which the payoff demands reaching or crossing of a barrier (predetermined price) by the underlying product. They include call options and put options, and are similar to common options in many aspects. Barrier options become active/inactive when the underlying product crosses the barrier.
Barrier options can be grouped into two as knock-in options and knock-out op- tions. Knock-in barrier options are inactive options at the beginning, but become active on reaching the barrier. On the other hand knock-out options starts as active but become inactive on reaching the barrier. Many barrier options carry rebates, which are paid off to the holder on reaching the barrier.
Barrier options are available in both European and American forms. In barrier options trading, premiums are paid in advance. Barrier options come in 4 types like up & out, up & in, down & out, and down & in. Of these four types 4 we can
7 8 take either call or put- giving us total 8 single barrier types. Barriers also come in various other forms including double barriers, parisians, and partial time barriers. Here, our discussion only target on single barrier types.
2.1 Valuation Methods
The valuation of barrier options can be tricky, because unlike other simpler op- tions they are path-dependent — that is, the value of the option at any time depends not just on the underlying at that point, but also on the path taken by the underlying since if the barrier is crossed the option is initiated or exterminat- ed. Although the classical Black–Scholes approach does not directly apply, several more complex methods can be used.
An approach is to study the law of the maximum (or minimum) of the underly- ing. This approach gives explicit prices to barrier options under the Black-Scholes framework. Also when an exact formula is difficult to obtain, barrier options can be priced with the Monte Carlo path simulation. However, computing the Greeks (sensitivities) using this approach is numerically unstable.
2.2 Analytical Solution under the Black-Scholes Model
Using the methods of equivalent martingale pricing (risk neutral valuation princi- ple), the price of a down-and-out put barrier option at time zero with respect to maturity T , strike K, barrier Sb, initial stock price S is given by:
rT Q
½ T Pdop(ST ,K,Sb)= e− E [(K ST ) K ST ,m Sb ] (2.1) − { ≥ 0 ≥ } 9
Note:
1. The realized minimum value of underlying stock from time zero to time t is
t defined as m0 = min Su. 0
2. The indicator function ½A(X) where:
1 if X A
½A(X)= ∈ 0 otherwise
is a useful notational device, in probability theory: if X is a probability space ½ with probability measure È and A is a measurable set, then A becomes a
random variable whose expected value is equal to the probability of A:
½ È È E(½A)= A(x)d = d = P (A) (2.2) X A
In view of linearity of the expectation, the above equation of down-and-out put
option price Pdop(ST ,K,Sb) can break into four terms:
Pdop(ST ,K,Sb)
rT Q ½ = e− (E [(K ST )½ ST − { } − − { } ½ ½ T T E[(K ST ) ST >Sb;m In the above equation, the first term is the price of a vanilla European put option under the Black-Scholes model. 10 The second term is similar to the Black Scholes put price. rT Q rT e− E [(K ST )½ ST 2 ln St +(r σ )(T t) Sb 2 d1,2(Sb)= ± − σ√T t − The value of the last two terms in price function Pdop(ST ,K,Sb) requires the T determination of the joint distribution function of ST and m0 . Therefore, we apply the following theorems to solve the problem. Theorem 1 (Reflection Principle) • 0 Let Wt denote the Brownian motion that starts at zero, with constant volatil- ity σ and zero drift rate. ξ = inf t : W 0 = h . { t } 0 Wt , t<ξ ˆ 0 Wt = 0 2h Wt , ξ Theorem 2 (Girsanov’s Theorem) • Girsanov’s Theorem begin with a Brownian motion under measure È and then construct a new measure É under which a “translated” process is a Brownian motion. P So assume that Wt t 0 defined on (Ω,F, È) is a Brownian motion starting { } ≥ from zero. Now define for an process γs s 0 adapted to the filtration Ft t 0 { } ≥ { } ≥ t P 1 t 2 R0 γs dWs 2 R0 γs ds Zt = e − (2.5) 11 E[Zt] = 1 and for each T > 0 define a new probability measure É on Ft t 0 { } ≥ by ½ É = E( Z ) for A F (2.6) A T ∈ T Define a process W Q = W P t γ ds. Then for each fixed T < the process t t − 0 s ∞ Q Wt is a Brownian motion on (Ω, FT , É) We illustrate how the reflection principle is applied to derive the joint law of the minimum value over [0,T ] and terminal value of a Brownian motion. Regarding u t to pricing a barrier option, we would like to find P (Wt > k,m0 < h) To figure out the last two terms, we construct three transformation equations to apply the reflection principle and the Girsanov’s Theorem. r 1 σ2 g = − 2 (2.7) σ 1 S h = ln( b ) (2.8) σ S0 1 S k = ln( b ) (2.9) σ S0 12 To calculate the last term of (2.4): rT Q ½ T e− E [(K ST ) ST >K;m Now we can re-plug in g,h,k, thus get the expression of the last term of (2.4) Then we put all derivations into a formula that represents Pdop to be the val- ue of this down-and-out put barrier option: rT p (S ,K,S )= Ke− Φ( d ) S Φ( d )+ S Φ( x ) dop 0 b − 2 − 0 − 1 0 − 1 rT Sb 2η Ke− Φ( x1 + σ√T ) S0( ) [Φ(y) Φ(y1)] − − − S0 − rt Sb 2η 2 + Ke− ( ) − [Φ(y σ√T ) Φ(y1 σ√T )] S0 − − − with S2 r + 1 σ2 ln( b ) η = 2 , y = S0K + ησ√T, σ2 σ√T ln( S0 ) ln( Sb ) Sb S0 x1 = + ησ√T, y1 = + ησ√T σ√T σ√T (2.10) 13 This analytical solution for continuously monitored down-and-out put barrier op- tion is also in accordance with the formula given in [13]: rT P (S ,K,S )= S e− (Φ(d ) Φ(d ) b[Φ(d ) Φ(d )]) dop o b − 0 3 − 1 − 8 − 6 rT + Ke− (Φ(d ) N(d ) a[Φ(d ) Φ(d )]) (2.11) 4 − 2 − 7 − 5 with 2r Sb 1+ 2 a =( )− σ S0 2 Sb 1+ r b =( ) σ2 S0 S0 1 2 ln K +(r + 2 σ )T d1 = σ√T S0 1 2 ln K +(r 2 σ )T d2 = − σ√T ln S0 +(r + 1 σ2)T Sb 2 d3 = σ√T ln S0 +(r 1 σ2)T Sb 2 d4 = − σ√T ln S0 (r 1 σ2)T Sb 2 d5 = − − σ√T ln S0 (r + 1 σ2)T Sb 2 d6 = − σ√T S0K 1 2 ln S2 (r 2 σ )T b − − d7 = σ√T S0K 1 2 ln S2 (r + 2 σ )T b − d8 = σ√T 14 2.2.1 A Continuity Correction to Discrete Barrier Option The barrier option pricing formulas presented so far assume continuous monitor- ing of the barrier. In practice, the barrier is normally monitored only at discrete points in time. An exception is the currency options market, where the barrier is frequently monitored almost continuously. For equity option in our case, the barrier is typically monitored against an periodical time point price. Discrete monitoring will naturally lower the probability of barrier hits compared with con- tinuous barrier monitoring. Broadie, Glasserman and Kou [7] showed theoretically and through examples that discrete barrier options can be priced with remarkable accuracy using the following simple correction. Sb is the continuous barrier options formulas with a discrete barrier level Sd e- qual to βσ√ t Sd = Sbe △ if the barrier is above the underlying security, and to βσ√ t Sd = Sbe− △ if the barrier is below the underlying security. At is the time between monitoring evens, and β 0.5826. ≈ 15 2.3 Monte Carlo Methods Monte Carlo methods are based on the analogy between probability and volume. The mathematics of measure formalizes the intuitive notion of probability, asso- ciating an event with a set of outcomes and defining the probability of the event to be its volume or measure relative to that of a universe of possible outcomes. Monte Carlo uses this identity in reverse, calculating the volume of a set by in- terpreting the volume as a probability. In the simplest case, this means sampling randomly from a universe of possible outcomes and taking the fraction of random draws that fall in a given set as an estimate of the set’s volume. The law of large numbers ensures that this estimate converges to the correct value as the number of draws increases. The central limit theorem provides information about the likely magnitude of the error in the estimate after a finite number of draws. In mathematical finance, a Monte Carlo option model uses Monte Carlo meth- ods to calculate the value of an option with multiple sources of uncertainty or with complicated features. The term ‘Monte Carlo method’ was coined by Stanislaw Ulam in the 1940s. The first application to option pricing was by Phelim Boyle in 1977 (for European options). In 1996, M. Broadie and P. Glasserman showed how to price Asian op- tions by Monte Carlo. In 2001 F. A. Longstaff and E. S. Schwartz developed a practical Monte Carlo method for pricing American-style options. We look at a simple example of European call option. Consider the option grant- ing the holder the right to buy the stock at a fixed price K at a fixed time T in 16 the future; the current time is t = 0. The payoff to option holder is (S K)+ = max[0,S K] (2.12) T − T − rT To get the present value of this payoff we multiply by a discount factor e− . We rT + denote the expected present value by E[e− (S K) ] T − As we know we apply the risk neutral measure into the calculation of the above Q expectation. The solution to the Black-Scholes SDE dSt/St = rdt + σdWt under risk neutral measure is 1 2 [(r 2 σ )T +σWT ] ST = S0e − (2.13) The random variable WT is normally distributed with mean 0 and variance T . We may therefore represent the terminal stock price as 1 2 [(r 2 σ )T +σ√TZ] ST = S0e − (2.14) where Z is a standard normal variable. Then we can simulate the path of stock price and apply the Monte Carlo sampling rT + methods to estimate price of the option by the expectation form E[e− (S K) ]. T − 2.3.1 The Main Framework of Monte Carlo Suppose we want to estimate some quantity θ = E[h(X)], where X = X ,X , ..., X { 1 2 n} n n is a random vector in R . h( ) is a function from R to R, and E[ h(X) ] < . | | ∞ Note that X could represent the values of a stochastic process at different points in time. For example, Xi might be the price of a particular stock at time i and 17 h( ) might be given by: X + X + ... + X h(X)= 1 2 n (2.15) n so then θ is the expected average value of the stock price. To estimate θ we use the following algorithm: Monte Carlo Algorithm for i=1 to n generate Xi set hi=h(Xi) h1+h2+...+hn set θn= n