Black-Scholes modeli, matematiksel finansta bir opsiyon fiyatlama modelidir. İsmini, bu modeli 1973 yılında yayınlayan^[1] Fischer Black ve Myron Scholes'tan almıştır. Bu opsiyon modelinin sonucunda, halen opsiyon fiyatlamada piyasa katılımcılarınca yoğun olarak kullanılmakta olan Black-Scholes formülü elde edilmiştir. Black-Scholes modeli, aslında rassal hareketler izleyen sıvı moleküllerini ortaya koyan Brown hareketinin hisse fiyatlarına ve finansal hareketlere uyarlanması sonucu ortaya çıkmıştır. Daha önce bu uyarlamanın öncüsü sayılabilecek varsayımı Louis Bachelier 1900'de "Théorie de la spéculation" başlığıyla yazdığı doktora tezinde^[2] yapmıştır. Yine, benzer uyarlamalar Paul Samuelson, Sheen Kassouf, Edward O. Thorp and Case Sprenkle tarafından da yapılmıştır. Ancak, Black ve Scholes'un zamandaşlarının önüne geçtiği nokta opsiyon fiyatlarına ihtiyaç duyan opsiyon piyasa katılımcılarına piyasada gözlemlenen veri ve değişkenlerle pratik bir şekilde hesaplanabilen analitik bir formül ortaya koymalarıdır.

Robert Merton'un modelde çözülemeyen bir bölümü çözmesinden sonra, model, Black-Scholes-Merton modeli olarak anılmaya da başlamıştır. Bu çalışmaları sayesinde, Merton ve Scholes, 1997de Ekonomi alanında Nobel Ödülü almışlardır.^[3]

Black-Scholes modeli ve bunun sonucunda elde edilen Black-Scholes formülü şu varsayımlara dayanmaktadır:

• Söz konusu dayanak varlığın (Black-Scholes özelinde hisse senedinin) fiyatının hareketleri (S_t) geometrik Brown hareketini izlemektedir. Yani, sabit bir sapma (${\displaystyle \mu }$) ve volatilite (${\displaystyle \sigma }$) olmak üzere;

${\displaystyle dS_{t}=\mu S_{t}\,dt+\sigma S_{t}\,dW_{t}\,}$

• Söz konusu hissede açığa satış (short sell) yapılması mümkündür.
• Arbitraj imkânı yoktur.
• Hisselerde el değiştirme süreklidir.
• Alımsatım maliyeti veya vergi yoktur.
• Bütün yatırım araçları kesirli bir şekilde alınıp satılabilmelidir; örneğin, bir dayanak varlığın yüzde birini almak mümkün olmalıdır.
• Risksiz faiz ile borç alınabilmelidir.
• Hisse temettü dağıtmamalıdır; bu kural sadece basit Black-Scholes modeli için geçerlidir.

Black-Scholes formülü 1973 yılında Fischer Black ve Myron Scholes tarafından yazılan makalede^[1] ilk defa bahsedilen ve Black-Scholes modeline bağlı olarak elde edilmiş bir opsiyon fiyatlama formülüdür. Avrupa tipi ödenişleri olan alım ve satım opsiyonlarının fiyatlanmasında piyasa katılımcılarınca yoğun olarak kullanılmaktadır.

Black–Scholes modelinin varsayımları altında, Avrupa tipi alım opsiyonu (European call option) için,

• opsiyon kullanma fiyatı K
• hissenin şu andaki fiyatı S, (yani opsiyonun verdiği hak ile T zaman sonra, hisseyi K fiyatından alma imkânımız var),
• sabit faiz r ve sabit volatilite ${\displaystyle \sigma }$

olmak üzere, opsiyonun bugünkü fiyatı ${\displaystyle C(S,T)}$ şu şekilde verilir.

${\displaystyle C(S,T)=SN(d_{1})-Ke^{-rT}N(d_{2})\,}$

Burada;

${\displaystyle d_{1}={\frac {\ln(S/K)+(r+\sigma ^{2}/2)T}{\sigma {\sqrt {T))))}$

${\displaystyle d_{2}={\frac {\ln(S/K)+(r-\sigma ^{2}/2)T}{\sigma {\sqrt {T))))=d_{1}-\sigma {\sqrt {T)).}$

Bu formülde ${\displaystyle N}$ standart normal dağılımın kümülatif dağılım fonksiyonudur.

Bir satım opsiyonunun fiyatı ${\displaystyle P(S,T)}$, yukarıda verilen ${\displaystyle C(S,T)}$ formülü ve alım-satım paritesi (put-call parity) kullanılarak hesaplanabilir ve aşağıdaki şekilde düzenlenebilir:

${\displaystyle P(S,T)=Ke^{-rT}N(-d_{2})-SN(-d_{1}).\,}$

Black-Scholes formülünün kanıtı bugün Black-Scholes kısmi diferansiyel denklemi olarak bilinen diferansiyel denklemlerin çözümünden geçmektedir.^[1] Burada esas ilk fikir, opsiyon ve opsiyon dayanak varlığından oluşan bir portföy yaratmak ve bu portföyü küçük zaman aralıklarında dayanak varlığın piyasa fiyatına duyarsız hale getirmektir. Sonucunda, Black-Scholes kısmi diferansiyel denklemi elde edilir. İkinci esas fikir ise bu diferansiyel denklemi, değişik dönüşümler ve yerine koymalar vasıtasıyla ısı denklemine dönüştürmektir.

Bu yöntemde Girsanov teoremi aracılığıyla riske duyarsız ölçü olan ${\displaystyle \mathbb {Q} }$'ya geçilir ve iskontolu hisse fiyatı sürecinin martingal olduğu elde edilir. Bu sayede, opsiyonun vadesindeki ödenişin iskontolu halinin beklenen değeri kolaylıkla hesaplanabilir. Diyelim ki bir olasılık ölçüsü ${\displaystyle \mathbb {P} }$'de

${\displaystyle dS_{t}=\mu S_{t}dt+\sigma S_{t}dW_{t}^{\mathbb {P} ))$

verilmiş olsun. Girsanov teoreminde ${\displaystyle \theta ={\frac {\mu -r}{\sigma ))}$ alırsak, o zaman ${\displaystyle W_{t}^{\mathbb {Q} }:=W_{t}^{\mathbb {P} }+\int _{0}^{t}\theta ds}$ yeni ölçüde de bir Brown hareketi olur ve ${\displaystyle dW_{t}=dW_{t}^{\mathbb {Q} }-{\frac {\mu -r}{\sigma ))dt}$ sağlanır. Bu halde, herhangi bir ${\displaystyle \tau \geq 0}$ için

${\displaystyle {\begin{aligned}d(e^{-r\tau }S_{t})&=d(e^{-r\tau })S_{t}+e^{-r\tau }d(S_{t})\\&=-re^{-r\tau }S_{t}dt+e^{-r\tau }dS_{t}\\&=-re^{-r\tau }S_{t}dt+e^{-r\tau }S_{t}(\mu S_{t}dt+\sigma S_{t}dW_{t}^{\mathbb {P} })\\&=\left(e^{-r\tau }S_{t}\right)(\mu -r)dt+\sigma \left(e^{-r\tau }S_{t}\right)dW_{t}^{\mathbb {P} }\\&=(\mu -r)\left(e^{-r\tau }S_{t}\right)dt+\sigma \left(e^{-r\tau }S_{t}\right)(dW_{t}^{\mathbb {Q} }-{\frac {\mu -r}{\sigma ))dt)\\&=\sigma \left(e^{-r\tau }S_{t}\right)dW_{t}^{\mathbb {Q} }\end{aligned))}$

olur. Elimizde sadece difüzyon terimleri kaldığı için artık iskontolu dayanak varlık spot fiyatı sürecinin riske duyarsız ölçüde martingal olduğu açıktır. Aynı zamanda,

${\displaystyle {\begin{aligned}dS_{t}&=\mu S_{t}dt+\sigma S_{t}dW_{t}=\mu S_{t}dt+\sigma S_{t}\left(dW_{t}^{\mathbb {Q} }-{\frac {\mu -r}{\sigma ))dt\right)\\&=\mu S_{t}dt+\sigma S_{t}dW_{t}^{\mathbb {Q} }-\mu S_{t}dt+rS_{t}dt\\&=rS_{t}dt+\sigma S_{t}dW_{t}^{\mathbb {Q} }.\end{aligned))}$

olur. Diğer deyişle, riske duyarsız ölçüde, dayanak varlık spot fiyatı sürecinin ${\displaystyle \mathbb {P} }$'deki deterministik terimi ${\displaystyle \mu }$, riske duyarsız ölçüde, risksiz faiz ${\displaystyle r}$ ile yer değiştirir. Öbür taraftan, Ito önsavı sayesinde

${\displaystyle {\begin{aligned}d\ln(S_{t})&={\frac {1}{S_{t))}dS_{t}+{\frac {1}{2))\left({\frac {-1}{S_{t}^{2))}\right)dS_{t}dS_{t}\\&=\left({\frac {1}{S_{t))}\right)\left(rS_{t}dt+\sigma S_{t}dW_{t}^{\mathbb {Q} }\right)-\left({\frac {1}{2S_{t}^{2))}\right)\sigma ^{2}S_{t}^{2}dt\\&=\left(r-{\frac {1}{2))\sigma ^{2}\right)dt+\sigma dW_{t}^{\mathbb {Q} }\end{aligned))}$

hesaplanır. Her iki tarafta integraller uygun aralıklarda alındıktan sonra

${\displaystyle S_{T}=S_{t}e^{\left(r-{\frac {1}{2))\sigma ^{2}\right)(T-t)+\sigma \left(W_{T}-W_{t}\right)))$

bulunur. Geriye kalan, ${\displaystyle \tau =T-t}$ koyup bir alım opsiyonun ${\displaystyle t}$ zamanındaki fiyatının ${\displaystyle C_{t}(K,T)=\mathbb {E} ^{\mathbb {Q} }\left[e^{-r\tau }\max(S_{T}-K)\right]}$ olduğundan yola çıkarak, beklenen değeri tanımı gereği integrale çevirip hesaplamaktır. Yani,

${\displaystyle {\begin{aligned}C_{t}(K,T)&=\mathbb {E} _{t}^{\mathbb {Q} }\left[e^{-r\tau }(S_{T}-K)^{+}|{\mathcal {F))_{t}\right]\\&=e^{-r\tau }\int \max(S_{T}-K,0)\ d\mathbb {Q} \\&=e^{-r\tau }{\frac {1}{\sqrt {2\pi ))}\int _{-\infty }^{\infty }\max(S_{t}e^{\left(r-{\frac {1}{2))\sigma ^{2}\right)(T-t)+\sigma {\sqrt {T-t))x}-K,0)e^{-{\frac {x^{2)){2\pi ))}dx.\end{aligned))}$

Integral fonksiyonunu ${\displaystyle \max }$ fonksiyonundan kurtarmak için ${\displaystyle (S_{t}e^{\left(r-{\frac {1}{2))\sigma ^{2}\right)(T-t)+\sigma {\sqrt {T-t))x}-K)\geq 0\iff x\geq -d_{2))$ olduğu hesaplanır. Burada

${\displaystyle d_{2}:={\frac {1}{\sigma {\sqrt {T-t))))\left(\ln \left({\frac {S_{t)){K))\right)+\left(r-{\frac {1}{2))\sigma ^{2}\right)(T-t)\right)}$

alınmıştır. Buradan sonra kalkülüs teknikleri kullanılarak integraller ilk önce ${\displaystyle d_{1}=d_{2}+\sigma {\sqrt {T-t))}$ tanımlanarak

${\displaystyle C_{t}(K,T)=S_{t}\int _{-d_{1))^{\infty }e^{-{\frac {y^{2)){2))}dx-e^{-r\tau }K(1-N(-d_{2}))}$

haline getirilir. Burada ${\displaystyle N}$ standart normal dağılımın kümülatif dağılım fonksiyonudur. ${\displaystyle d_{1))$ ve ${\displaystyle d_{2))$ arasındaki ilişki ve ${\displaystyle N}$'nin özellikleri kullanılarak

${\displaystyle C_{t}(K,T)=S_{t}N(d_{1})-Ke^{-r\tau }N(d_{2})}$

bulunur. Black-Scholes formülünün orijinal hali ${\displaystyle t=0}$ alınarak elde edilir.

Black-Scholes formülü üzerinden bir Avrupa tipi opsiyonun risk hassasiyetleri analitik olarak hesaplanabilir.

		Açıklama	Alım opsiyonu	Satım opsiyonu
Delta	${\displaystyle {\frac {\partial V}{\partial S))}$	Dayanak varlığın spot fiyatına göre değişim	${\displaystyle N(d_{1})\,}$	${\displaystyle -N(-d_{1})=N(d_{1})-1\,}$
Gama	${\displaystyle {\frac {\partial ^{2}V}{\partial S^{2))))$	Deltaya göre değişim	${\displaystyle {\frac {N'(d_{1})}{S\sigma {\sqrt {T-t))))\,}$
Vega [4][5]	${\displaystyle {\frac {\partial V}{\partial \sigma ))}$	Volatiliteye göre değişim	${\displaystyle SN'(d_{1}){\sqrt {T-t))\,}$
Thita	${\displaystyle {\frac {\partial V}{\partial t))}$	Vade gününe kalan zamana göre değişim	${\displaystyle -{\frac {SN'(d_{1})\sigma }{2{\sqrt {T-t))))-rKe^{-r(T-t)}N(d_{2})\,}$	${\displaystyle -{\frac {SN'(d_{1})\sigma }{2{\sqrt {T-t))))+rKe^{-r(T-t)}N(-d_{2})\,}$
Ro	${\displaystyle {\frac {\partial V}{\partial r))}$	Faiz oranına göre değişim	${\displaystyle K(T-t)e^{-r(T-t)}N(d_{2})\,}$	${\displaystyle -K(T-t)e^{-r(T-t)}N(-d_{2})\,}$

• Black-Scholes denklemi
• Binom opsiyon fiyatlama modeli
• Black modeli
• Sıçramalı difüzyon modeli
• Monte Carlo opsiyon fiyatlama modeli
• Stokastik volatilite