The above model can be extended for variable (but deterministic) rates and volatilities. The model may also be used to value European options on instruments paying dividends. In this case, closed-form solutions are available if the dividend is a known proportion of the stock price.
American options and options on stocks paying a known cash dividend (in the short term, more realistic than a proportional dividend) are more difficult to value, and a choice of solution techniques is available (for example
lattices and
grids).
Instruments paying continuous yield dividends For options on indices, it is reasonable to make the simplifying assumption that dividends are paid continuously, and that the dividend amount is proportional to the level of the index. The dividend payment paid over the time period [t, t + dt] is then modelled as: :qS_t\,dt for some constant q (the
dividend yield). Under this formulation the arbitrage-free price implied by the Black–Scholes model can be shown to be: :C(S_t, t) = e^{-r(T - t)}[FN(d_1) - KN(d_2)]\, and :P(S_t, t) = e^{-r(T - t)}[KN(-d_2) - FN(-d_1)]\, where now :F = S_t e^{(r - q)(T - t)}\, is the modified forward price that occurs in the terms d_1, d_2: :d_1 = \frac{1}{\sigma\sqrt{T - t}}\left[\ln\left(\frac{S_t}{K}\right) + \left(r - q + \frac{1}{2}\sigma^2\right)(T - t)\right] and :d_2 = d_1 - \sigma\sqrt{T - t} = \frac{1}{\sigma\sqrt{T - t}}\left[\ln\left(\frac{S_t}{K}\right) + \left(r - q - \frac{1}{2}\sigma^2\right)(T - t)\right].
Instruments paying discrete proportional dividends It is also possible to extend the Black–Scholes framework to options on instruments paying discrete proportional dividends. This is useful when the option is struck on a single stock. A typical model is to assume that a proportion \delta of the stock price is paid out at pre-determined times t_1, t_2, \ldots, t_n . The price of the stock is then modelled as: :S_t = S_0(1 - \delta)^{n(t)}e^{ut + \sigma W_t} where n(t) is the number of dividends that have been paid by time t. The price of a call option on such a stock is again: :C(S_0, T) = e^{-rT}[FN(d_1) - KN(d_2)]\, where now :F = S_{0}(1 - \delta)^{n(T)}e^{rT}\, is the forward price for the dividend paying stock.
American options The problem of finding the price of an
American option is related to the
optimal stopping problem of finding the time to execute the option. Since the American option can be exercised at any time before the expiration date, the Black–Scholes equation becomes a variational inequality of the form: :\frac{\partial V}{\partial t} + \frac{1}{2}\sigma^2 S^2 \frac{\partial^2 V}{\partial S^2} + rS\frac{\partial V}{\partial S} - rV \leq 0 together with V(S, t) \geq H(S) where H(S) denotes the payoff at stock price S and the terminal condition: V(S, T) = H(S). In general this inequality does not have a closed form solution, though an American call with no dividends is equal to a European call and the Roll–Geske–Whaley method provides a solution for an American call with one dividend; see also
Black's approximation. Barone-Adesi and Whaley is a further approximation formula. Here, the stochastic differential equation (which is valid for the value of any derivative) is split into two components: the European option value and the early exercise premium. With some assumptions, a
quadratic equation that approximates the solution for the latter is then obtained. This solution involves
finding the critical value, s*, such that one is indifferent between early exercise and holding to maturity. Bjerksund and Stensland provide an approximation based on an exercise strategy corresponding to a trigger price. Here, if the underlying asset price is greater than or equal to the trigger price it is optimal to exercise, and the value must equal S - X, otherwise the option "boils down to: (i) a European
up-and-out call option... and (ii) a rebate that is received at the knock-out date if the option is knocked out prior to the maturity date". The formula is readily modified for the valuation of a put option, using
put–call parity. This approximation is computationally inexpensive and the method is fast, with evidence indicating that the approximation may be more accurate in pricing long dated options than Barone-Adesi and Whaley.
Perpetual put Despite the lack of a general analytical solution for American put options, it is possible to derive such a formula for the case of a perpetual option – meaning that the option never expires (i.e., T\rightarrow \infty). In this case, the time decay of the option is equal to zero, which leads to the Black–Scholes PDE becoming an ODE:{1\over{2}}\sigma^{2}S^{2}{d^{2}V\over{dS^{2}}} + (r-q)S{dV\over{dS}} - rV = 0Let S_{-} denote the lower exercise boundary, below which it is optimal to exercise the option. The boundary conditions are:V(S_{-}) = K-S_{-}, \quad {dV\over{dS}}(S_{-}) = -1, \quad V(S) \leq KThe solutions to the ODE are a linear combination of any two linearly independent solutions:V(S) = A_{1}S^{\lambda_{1}} + A_{2}S^{\lambda_{2}}For S_{-} \leq S, substitution of this solution into the ODE for i = {1,2} yields:\left[ {1\over{2}}\sigma^{2}\lambda_{i}(\lambda_{i}-1) + (r-q)\lambda_{i} - r \right]S^{\lambda_{i}} = 0Rearranging the terms gives:{1\over{2}}\sigma^{2}\lambda_{i}^{2} + \left(r-q - {1\over{2}} \sigma^{2}\right)\lambda_{i} - r = 0Using the
quadratic formula, the solutions for \lambda_{i} are:\begin{aligned} \lambda_{1} &= {-\left(r-q-{1\over{2}}\sigma^{2} \right ) + \sqrt{\left(r-q-{1\over{2}}\sigma^{2} \right )^{2} + 2\sigma^{2}r}\over{\sigma^{2}}} \\ \lambda_{2} &= {-\left(r-q-{1\over{2}}\sigma^{2} \right ) - \sqrt{\left(r-q-{1\over{2}}\sigma^{2} \right )^{2} + 2\sigma^{2}r}\over{\sigma^{2}}} \end{aligned}In order to have a finite solution for the perpetual put, since the boundary conditions imply upper and lower finite bounds on the value of the put, it is necessary to set A_{1} = 0, leading to the solution V(S) = A_{2}S^{\lambda_{2}}. From the first boundary condition, it is known that:V(S_{-}) = A_{2}(S_{-})^{\lambda_{2}} = K-S_{-} \implies A_{2} = {K-S_{-}\over{(S_{-})^{\lambda_{2}}}}Therefore, the value of the perpetual put becomes:V(S) = (K-S_{-})\left( {S\over{S_{-}}} \right)^{\lambda_{2}}The second boundary condition yields the location of the lower exercise boundary:{dV\over{dS}}(S_{-}) = \lambda_{2}{K-S_{-}\over{S_{-}}} = -1 \implies S_{-} = {\lambda_{2}K\over{\lambda_{2}-1}}To conclude, for S \geq S_{-} = {\lambda_{2}K\over{\lambda_{2}-1}}, the perpetual American put option is worth:V(S) = {K\over{1-\lambda_{2}}} \left( {\lambda_{2}-1\over{\lambda_{2}}}\right)^{\lambda_{2}} \left( {S\over{K}} \right)^{\lambda_{2}}
Binary options By solving the Black–Scholes differential equation with the
Heaviside function as a boundary condition, one ends up with the pricing of options that pay one unit above some predefined strike price and nothing below. In fact, the Black–Scholes formula for the price of a vanilla call option (or put option) can be interpreted by decomposing a call option into an asset-or-nothing call option minus a cash-or-nothing call option, and similarly for a put—the binary options are easier to analyze, and correspond to the two terms in the Black–Scholes formula.
Cash-or-nothing call This pays out one unit of cash if the spot is above the strike at maturity. Its value is given by: : C =e^{-r (T-t)}N(d_2). \,
Cash-or-nothing put This pays out one unit of cash if the spot is below the strike at maturity. Its value is given by: : P = e^{-r (T-t)}N(-d_2). \,
Asset-or-nothing call This pays out one unit of asset if the spot is above the strike at maturity. Its value is given by: : C = Se^{-q (T-t)}N(d_1). \,
Asset-or-nothing put This pays out one unit of asset if the spot is below the strike at maturity. Its value is given by: : P = Se^{-q (T-t)}N(-d_1),
Foreign Exchange (FX) Denoting by
S the FOR/DOM exchange rate (i.e., 1 unit of foreign currency is worth S units of domestic currency) one can observe that paying out 1 unit of the domestic currency if the spot at maturity is above or below the strike is exactly like a cash-or nothing call and put respectively. Similarly, paying out 1 unit of the foreign currency if the spot at maturity is above or below the strike is exactly like an asset-or nothing call and put respectively. Hence by taking r_{f}, the foreign interest rate, r_{d}, the domestic interest rate, and the rest as above, the following results can be obtained: In the case of a digital call (this is a call FOR/put DOM) paying out one unit of the domestic currency gotten as present value: : C = e^{-r_{d} T}N(d_2) \, In the case of a digital put (this is a put FOR/call DOM) paying out one unit of the domestic currency gotten as present value: : P = e^{-r_{d}T}N(-d_2) \, In the case of a digital call (this is a call FOR/put DOM) paying out one unit of the foreign currency gotten as present value: : C = Se^{-r_{f} T}N(d_1) \, In the case of a digital put (this is a put FOR/call DOM) paying out one unit of the foreign currency gotten as present value: : P = Se^{-r_{f}T}N(-d_1) \,
Skew In the standard Black–Scholes model, one can interpret the premium of the binary option in the risk-neutral world as the expected value = probability of being in-the-money * unit, discounted to the present value. The Black–Scholes model relies on symmetry of distribution and ignores the
skewness of the distribution of the asset. Market makers adjust for such skewness by, instead of using a single standard deviation for the underlying asset \sigma across all strikes, incorporating a variable one \sigma(K) where volatility depends on strike price, thus incorporating the
volatility skew into account. The skew matters because it affects the binary considerably more than the regular options. A binary call option is, at long expirations, similar to a tight call spread using two vanilla options. One can model the value of a binary cash-or-nothing option,
C, at strike
K, as an infinitesimally tight spread, where C_v is a vanilla European call: : C = \lim_{\epsilon \to 0} \frac{C_v(K-\epsilon) - C_v(K)}{\epsilon} Thus, the value of a binary call is the negative of the
derivative of the price of a vanilla call with respect to strike price: : C = -\frac{dC_v}{dK} When one takes volatility skew into account, \sigma is a function of K: : C = -\frac{dC_v(K,\sigma(K))}{dK} = -\frac{\partial C_v}{\partial K} - \frac{\partial C_v}{\partial \sigma} \frac{\partial \sigma}{\partial K} The first term is equal to the premium of the binary option ignoring skew: : -\frac{\partial C_v}{\partial K} = -\frac{\partial (S N(d_1) - Ke^{-r(T-t)} N(d_2))}{\partial K} = e^{-r (T-t)} N(d_2) = C_\text{no skew} \frac{\partial C_v}{\partial \sigma} is the
Vega of the vanilla call; \frac{\partial \sigma}{\partial K} is sometimes called the "skew slope" or just "skew". If the skew is typically negative, the value of a binary call will be higher when taking skew into account. : C = C_\text{no skew} - \text{Vega}_v \cdot \text{Skew}
Relationship to vanilla options' Greeks Since a binary call is a mathematical derivative of a vanilla call with respect to strike, the price of a binary call has the same shape as the delta of a vanilla call, and the delta of a binary call has the same shape as the gamma of a vanilla call. ==Black–Scholes in practice==