Vasicek model

The model specifies that the instantaneous interest rate follows the stochastic differential equation:

${\displaystyle dr_{t}=a(b-r_{t})\,dt+\sigma \,dW_{t}}$ 

where W_t is a Wiener process under the risk neutral framework modelling the random market risk factor, in that it models the continuous inflow of randomness into the system. The standard deviation parameter, ${\displaystyle \sigma }$ , determines the volatility of the interest rate and in a way characterizes the amplitude of the instantaneous randomness inflow. The typical parameters ${\displaystyle b,a}$  and ${\displaystyle \sigma }$ , together with the initial condition ${\displaystyle r_{0}}$ , completely characterize the dynamics, and can be quickly characterized as follows, assuming ${\displaystyle a}$  to be non-negative:

• ${\displaystyle b}$ : "long term mean level". All future trajectories of ${\displaystyle r}$  will evolve around a mean level b in the long run;
• ${\displaystyle a}$ : "speed of reversion". ${\displaystyle a}$  characterizes the velocity at which such trajectories will regroup around ${\displaystyle b}$  in time;
• ${\displaystyle \sigma }$ : "instantaneous volatility", measures instant by instant the amplitude of randomness entering the system. Higher ${\displaystyle \sigma }$  implies more randomness

The following derived quantity is also of interest,

• ${\displaystyle {\sigma ^{2}}/(2a)}$ : "long term variance". All future trajectories of ${\displaystyle r}$  will regroup around the long term mean with such variance after a long time.

${\displaystyle a}$  and ${\displaystyle \sigma }$  tend to oppose each other: increasing ${\displaystyle \sigma }$  increases the amount of randomness entering the system, but at the same time increasing ${\displaystyle a}$  amounts to increasing the speed at which the system will stabilize statistically around the long term mean ${\displaystyle b}$  with a corridor of variance determined also by ${\displaystyle a}$ . This is clear when looking at the long term variance,

${\displaystyle {\frac {\sigma ^{2}}{2a}}}$ 

which increases with ${\displaystyle \sigma }$  but decreases with ${\displaystyle a}$ .

This model is an Ornstein–Uhlenbeck stochastic process.

Vasicek's model was the first one to capture mean reversion, an essential characteristic of the interest rate that sets it apart from other financial prices. Thus, as opposed to stock prices for instance, interest rates cannot rise indefinitely. This is because at very high levels they would hamper economic activity, prompting a decrease in interest rates. Similarly, interest rates do not usually decrease much below 0. As a result, interest rates move in a limited range, showing a tendency to revert to a long run value.

The drift factor ${\displaystyle a(b-r_{t})}$  represents the expected instantaneous change in the interest rate at time t. The parameter b represents the long-run equilibrium value towards which the interest rate reverts. Indeed, in the absence of shocks (${\displaystyle dW_{t}=0}$ ), the interest rate remains constant when r_t = b. The parameter a, governing the speed of adjustment, needs to be positive to ensure stability around the long term value. For example, when r_t is below b, the drift term ${\displaystyle a(b-r_{t})}$  becomes positive for positive a, generating a tendency for the interest rate to move upwards (toward equilibrium).

The main disadvantage is that, under Vasicek's model, it is theoretically possible for the interest rate to become negative, an undesirable feature under pre-crisis assumptions. This shortcoming was fixed in the Cox–Ingersoll–Ross model, exponential Vasicek model, Black–Derman–Toy model and Black–Karasinski model, among many others. The Vasicek model was further extended in the Hull–White model. The Vasicek model is also a canonical example of the affine term structure model, along with the Cox–Ingersoll–Ross model. In recent research both models were used for data partitioning and forecasting.^[2]

We solve the stochastic differential equation to obtain

${\displaystyle r_{t}=r_{0}e^{-at}+b\left(1-e^{-at}\right)+\sigma e^{-at}\int _{0}^{t}e^{as}\,dW_{s}.\,\!}$ 

Using similar techniques as applied to the Ornstein–Uhlenbeck stochastic process we get that state variable is distributed normally with mean

${\displaystyle \mathrm {E} [r_{t}]=r_{0}e^{-at}+b(1-e^{-at})}$ 

and variance

${\displaystyle \mathrm {Var} [r_{t}]={\frac {\sigma ^{2}}{2a}}(1-e^{-2at}).}$ 

Consequently, we have

${\displaystyle \lim _{t\to \infty }\mathrm {E} [r_{t}]=b}$ 

and

${\displaystyle \lim _{t\to \infty }\mathrm {Var} [r_{t}]={\frac {\sigma ^{2}}{2a}}.}$ 

Under the no-arbitrage assumption, a discount bond may be priced in the Vasicek model. The time ${\displaystyle t}$  value of a discount bond with maturity date ${\displaystyle T}$  is exponential affine in the interest rate:

${\displaystyle P(t,T)=A(t,T)e^{-B(t,T)r(t)}}$ 

where

${\displaystyle B(t,T)={\frac {1-e^{-a(T-t)}}{a}}}$ 

${\displaystyle A(t,T)=\exp \left\{\left(b-{\frac {\sigma ^{2}}{2a^{2}}}\right)\left[B(t,T)-(T-t)\right]-{\frac {\sigma ^{2}}{4a}}B^{2}(t,T)\right\}}$ 

• Ornstein–Uhlenbeck process.
• Hull–White model
• Cox–Ingersoll–Ross model

• Hull, John C. (2003). Options, Futures and Other Derivatives. Upper Saddle River, NJ: Prentice Hall. ISBN 978-0-13-009056-0.
• Damiano Brigo, Fabio Mercurio (2001). Interest Rate Models – Theory and Practice with Smile, Inflation and Credit (2nd ed. 2006 ed.). Springer Verlag. ISBN 978-3-540-22149-4.
• Jessica James, Nick Webber (2000). Interest Rate Modelling. Wiley. ISBN 978-0-471-97523-6.

• The Vasicek Model, Bjørn Eraker, Wisconsin School of Business
• Yield Curve Estimation and Prediction with the Vasicek Model, D. Bayazit, Middle East Technical University