The Black-Scholes option pricing formula gives the no-arbitrage price of a European call option as a function of current stock price S, strike K, time to expiration t, riskless discount factor RF, and σ (standard deviation of the log compound return of the stock):
C=S⋅N(z∗+σ)−K⋅RF−1⋅N(z∗)
where z∗=[ln(S/K)+lnRF−21σ2]/σ and N(⋅) is the standard normal cumulative distribution function (CDF). The Black-Scholes value is a function solely of observables — no preference parameters appear.
Key Ideas
Original derivation (Black-Scholes 1973): Continuous trading allows construction of a perfectly hedged portfolio (Δ shares long + short call) whose value is locally riskless; no-arbitrage forces it to earn rF. This yields the partial differential equation (PDE) 21v2x2w11+rxw1+w2−rw=0, solved analytically via a heat-equation substitution. Risk-neutral pricing — EQ[payoff]⋅RF−1 — is a consequence, not an assumption.
Capital Asset Pricing Model (CAPM) alternative derivation (Black-Scholes 1973): The option's instantaneous beta satisfies βw=(xw1/w)⋅β. Substituting into the CAPM expected-return equation, the preference parameters (α−r terms) cancel identically, recovering the same PDE without any hedging argument. Preference-independence is a consequence of the elasticity relationship, not a distributional assumption about the market risk premium.
Cox-Ross (1975) restatement: In a complete, frictionless market, the arbitrage-free price equals the price one would compute under risk neutrality. Continuous trading creates completeness; preference-independence follows from spanning.
Rubinstein (1976) discrete-time derivation: Under constant proportional risk aversion (CPRA) preferences, the pricing kernel is Yt=RMt−b. If St and Yt are jointly lognormal, the Black-Scholes formula holds exactly in discrete time with risk-averse investors — no continuous hedging required. The essential ingredient is the distributional assumption, not the trading technology.
How It Works
Delta Hedge and PDE (Black-Scholes 1973)
Assume stock follows geometric Brownian motion (GBM): dS=αSdt+vSdz. Consider a portfolio: long 1 share, short 1/w1 calls, where w1=∂w/∂x. The portfolio value is x−w/w1; its return has zero covariance with the market (the dz terms cancel exactly), so by no-arbitrage it must earn r. Writing out the return over dt via Itô's lemma and imposing =rdt gives:
21v2x2w11+rxw1+w2−rw=0(7)
With boundary condition w(x,t∗)=max(x−c,0) and substitution y=lnx+(r−21v2)(t∗−t), τ=v2(t∗−t), equation (7) reduces to the heat equation ∂F/∂τ=∂2F/∂y2. The standard convolution solution gives:
w(x,t)=xN(d1)−ce−r(t∗−t)N(d2)
where d1=[ln(x/c)+(r+v2/2)(t∗−t)]/(vt∗−t) and d2=d1−vt∗−t.
Corporate Liabilities as Options (Black-Scholes 1973)
Let x be the total market value of the firm's assets. With debt of face value B maturing at t∗:
Equity = w(x,t): a call on firm assets with strike B. Equity holders receive max(x−B,0) at maturity.
Bond value = x−w(x,t): the firm value minus the equity call.
Coupon bonds are compound options: each coupon payment is an option to continue (pay that coupon) and retain the right to pay future coupons.
This reframes all corporate securities as contingent claims on firm assets, independent of their priority structure, and directly precedes Merton's (1974) structural credit model.
From the Stochastic Discount Factor
Let Zt′=YtRF−1/E(Yt) be the normalized pricing kernel (risk-neutral probability kernel). The call price is:
C=E[(St−K)+Zt′]=E[(St−K)Zt′St≥K]
When (St,Yt) are jointly lognormal, the Stein–Rubinstein lemma evaluates the resulting bivariate integral:
Cov(x,g(y))=E[g′(y)]Cov(x,y)for bivariate normal (x,y)
Breaking the expectation into two integrals over the bivariate lognormal density and applying this identity yields the Black-Scholes formula analytically. See Stochastic Discount Factor for the full framework.
Put-Call Parity
The put payoff rewrites as max(0,K−S∗)=max(0,S∗−K)−S∗+K. Applying the valuation operator gives:
P=C−(S−D)+K⋅RF−1
where D is the present value of dividends through expiration. Holds by pure arbitrage regardless of the distributional assumption on S∗; puts are manufactured from calls in practice (Rubinstein 1987).
Binomial Method (Cox-Ross-Rubinstein 1979)
If the underlying can move only up or down by known percentages each discrete period, a portfolio of the index and cash exactly replicates the option. Working backward from expiration payoffs gives a unique option value at each node; for American options, replace the holding value with the exercise value whenever exercising is optimal. In the limit as the period length shrinks to zero, the binomial price converges to Black-Scholes.
Extensions
Stochastic interest rates (Merton 1973): If the bond price P follows any Itô process dP/P=μ(τ)dt+δ(τ)dq, the option price generalizes to H(S,P,τ;E)=EP(τ)⋅y[S/EP(τ),∫0τV2(s)ds] where V2=σ2+δ2−2ρσδ. CAPM is not needed; only no-sure-thing-profits among correlated securities and Itô's lemma are required.
American perpetual put (Merton 1973): Closed form G(S)=1+γE[(1+γ)S/γE]−γ, γ=2r/σ2; optimal exercise threshold C∗=γE/(1+γ).
Down-and-out barrier option (Merton 1973): For a knock-out barrier B[τ]=bEe−ητ, the price equals a standard call minus a discount term involving the complementary error function; first closed-form barrier option formula in the finance literature.
Known dividend yield: Replace S with SA−1 where A = product of (1+δt) over the option's life (Rubinstein 1976, Corollary to Theorem 5).
Stochastic volatility: When σ is random, Black-Scholes undervalues out-of-the-money (OTM) and in-the-money (ITM) options and generates a U-shaped implied volatility smile. The formula generalizes to E[BS(AV)] where AV = average variance over the option life; see Stochastic Volatility and Ball-Roma (1994).
Equilibrium foundation for stochastic volatility (Amin-Ng 1993): In a discrete-time CPRA framework with bivariate conditional lognormality of stock-consumption returns, Black-Scholes holds exactly when volatility is constant (Rubinstein 1976). When volatility is stochastic but predictable, the call price is E0[CBS(σ0T2)] — a preference-free formula with no risk-aversion parameters, providing an equilibrium basis for the Hull-White/Ball-Roma expected-Black-Scholes result. The interest rate is endogenous: rt=−lnρ+bμc,t+1−21b(1+b)hc,t+1.
Implied binomial trees (Rubinstein 1994): A nonparametric lattice calibrated to observed option prices across strikes, recovering the market's risk-neutral distribution without specifying a parametric model for σ.
Model error (Jacquier-Jarrow 2000): The gap between observed and model prices is treated as a random variable ηi with an estimated distribution, rather than a parameter to be fitted away. See Contingent Claim Model Error.
Nonparametric risk-neutral density (Aït-Sahalia-Lo 2000): The risk-neutral density (RND) is extracted from the implied-volatility surface via Nadaraya-Watson kernel regression, avoiding the lognormal distributional assumption of Black-Scholes. Constant-maturity (30-day) option prices are recovered by interpolating along the surface. Applied in Aloulou-Ellouze (2016) to study heterogeneous agent dynamics and the price formation process.
Why It Matters
Reference benchmark: Every derivative security is priced relative to Black-Scholes. Every empirical deviation — smile, skew, heavy tails — is measured against it.
Robustness: Rubinstein (1976) shows the formula does not depend on continuous trading; bivariate lognormality suffices. This explains its empirical durability even when continuous-trading assumptions clearly fail.
Spanning and completeness: The existence of a unique pricing kernel (complete markets) is the mathematical foundation; Black-Scholes is the closed-form solution in the lognormal case.
Preference-independence via two routes: Both the delta-hedge argument and the CAPM derivation confirm that the formula requires no knowledge of investors' risk preferences — the risk premium cancels algebraically. This makes it operationally useful even in markets where expected returns are unobservable.
Empirical baseline: Black-Scholes (1973) documented systematic over-pricing by option buyers relative to the formula, particularly for low-variance stocks. This evidence of market mis-pricing pre-dated the formula's adoption; subsequent empirical work (volatility smile, skew) refined the direction of deviations but confirmed the formula as the correct benchmark.
Open Questions
Empirically σ is stochastic and mean-reverting; extending to stochastic volatility (SV) while preserving tractability requires characteristic functions (Heston 1993) or the Moment-Generating Function (MGF)-of-average-variance approach (Ball-Roma 1994).
Nonzero correlation between stock and volatility breaks the joint lognormality argument and forces reliance on the full risk-neutral characteristic function.
The volatility smile (OTM implied vol > at-the-money (ATM)) cannot arise from the lognormal model; its sources (jumps, leverage, stochastic vol, fat tails) and the appropriate model extension remain actively debated.