Analytical Greeks for European Options

Medium·22 min read
Risk & GreeksBlack-ScholesGreeksP&L Decomposition

Setup

What Greeks Measure

A Greek is the partial derivative of an option's price with respect to one of its inputs. Greeks serve two purposes simultaneously:

  1. Hedging: each Greek tells you how much of a hedging instrument to hold to neutralise the corresponding risk.
  2. P&L attribution: the Taylor expansion of option P&L, decomposed by Greeks, explains where the desk's daily P&L came from.

Neither purpose is decorative. A trading desk that cannot decompose its P&L into Greeks contributions has no risk management — it has guessing.

Assumptions and Notation

All results in this module derive from the Black-Scholes formula. The assumptions are those of the Black-Scholes model: GBM dynamics, constant volatility σ\sigma, constant risk-free rate rr, no dividends.

Call price: C=SN(d1)KerτN(d2),C = S\,N(d_1) - Ke^{-r\tau}\,N(d_2), d1=ln(S/K)+(r+σ2/2)τστ,d2=d1στ,τ=Tt.d_1 = \frac{\ln(S/K) + (r + \sigma^2/2)\tau}{\sigma\sqrt{\tau}}, \qquad d_2 = d_1 - \sigma\sqrt{\tau}, \qquad \tau = T - t.

Put price: by put-call parity, P=KerτN(d2)SN(d1)P = Ke^{-r\tau}N(-d_2) - S\,N(-d_1).

Core identity (used throughout): Sn(d1)=Kerτn(d2),S\,n(d_1) = Ke^{-r\tau}\,n(d_2), where n(x)=12πex2/2n(x) = \frac{1}{\sqrt{2\pi}}e^{-x^2/2} is the standard normal pdf. This identity is proved by substituting the definitions of d1d_1 and d2d_2 and simplifying the exponentials.

Useful partial derivatives of d1d_1: d1S=1Sστ,d1σ=d2σ,d1τ=d22τ+rστ.\frac{\partial d_1}{\partial S} = \frac{1}{S\sigma\sqrt{\tau}}, \qquad \frac{\partial d_1}{\partial \sigma} = -\frac{d_2}{\sigma}, \qquad \frac{\partial d_1}{\partial \tau} = -\frac{d_2}{2\tau} + \frac{r}{\sigma\sqrt{\tau}}.

Note: d2/S=d1/S\partial d_2 / \partial S = \partial d_1 / \partial S (since d2=d1στd_2 = d_1 - \sigma\sqrt{\tau}, independent of SS).


Delta: Δ=C/S\Delta = \partial C / \partial S

Δcall=N(d1),Δput=N(d1)1=N(d1).\Delta_{\mathrm{call}} = N(d_1), \qquad \Delta_{\mathrm{put}} = N(d_1) - 1 = -N(-d_1).

Derivation: CS=N(d1)+Sn(d1)d1SKerτn(d2)d2S=N(d1)+Sn(d1)Kerτn(d2)=01Sστ=N(d1).\frac{\partial C}{\partial S} = N(d_1) + S\,n(d_1)\frac{\partial d_1}{\partial S} - Ke^{-r\tau}n(d_2)\frac{\partial d_2}{\partial S} = N(d_1) + \underbrace{S\,n(d_1) - Ke^{-r\tau}n(d_2)}_{=\,0} \cdot \frac{1}{S\sigma\sqrt{\tau}} = N(d_1).

The core identity eliminates the n(d1)n(d_1) and n(d2)n(d_2) terms.

Sign conventions:

  • Long call: Δ(0,1)\Delta \in (0, 1). At-the-money: Δ0.5\Delta \approx 0.5.
  • Long put: Δ(1,0)\Delta \in (-1, 0). At-the-money: Δ0.5\Delta \approx -0.5.
  • Δ1\Delta \to 1 as call goes deep ITM; Δ0\Delta \to 0 as call goes deep OTM.

Hedge interpretation: hold Δcall\Delta_{\mathrm{call}} shares per long call to be instantaneously delta-neutral.

Note: Δcall=N(d1)\Delta_{\mathrm{call}} = N(d_1) is not the risk-neutral probability Q(ST>K)=N(d2)\mathbb{Q}(S_T > K) = N(d_2). The difference N(d1)N(d2)>0N(d_1) - N(d_2) > 0 is the "probability gap" due to the lognormality of STS_T.


Gamma: Γ=2C/S2\Gamma = \partial^2 C / \partial S^2

Γ=n(d1)Sστ.\Gamma = \frac{n(d_1)}{S\sigma\sqrt{\tau}}.

Derivation: Γ=ΔS=N(d1)S=n(d1)d1S=n(d1)1Sστ.\Gamma = \frac{\partial \Delta}{\partial S} = \frac{\partial N(d_1)}{\partial S} = n(d_1)\frac{\partial d_1}{\partial S} = n(d_1) \cdot \frac{1}{S\sigma\sqrt{\tau}}.

Key properties:

  • Γcall=Γput\Gamma_{\mathrm{call}} = \Gamma_{\mathrm{put}}: identical for calls and puts with the same inputs. This follows from put-call parity: 2(CP)/S2=2(SKerτ)/S2=0\partial^2(C-P)/\partial S^2 = \partial^2(S - Ke^{-r\tau})/\partial S^2 = 0.
  • Γ>0\Gamma > 0 always: option value is convex in SS. Long options are long gamma.
  • Γ\Gamma is maximised at-the-money and decays away from the money.
  • Γ\Gamma \to \infty as τ0\tau \to 0 with S=KS = K: gamma spikes near expiry at-the-money (short-dated ATM options are extremely gamma-sensitive).

Trading interpretation: gamma income is the profit from rebalancing the delta hedge. A long gamma position earns 12Γ(dS)2\frac{1}{2}\Gamma (dS)^2 per unit time from spot moves — but pays theta to fund it.


Theta: Θ=C/t=C/τ\Theta = \partial C / \partial t = -\partial C / \partial \tau

Θcall=Sn(d1)σ2τrKerτN(d2).\Theta_{\mathrm{call}} = -\frac{S\,n(d_1)\,\sigma}{2\sqrt{\tau}} - rKe^{-r\tau}N(d_2).

Derivation (using τ=Tt\tau = T - t, so C/t=C/τ\partial C/\partial t = -\partial C/\partial\tau): Cτ=Sn(d1)d1τKerτ(rN(d2)+n(d2)d2τ).\frac{\partial C}{\partial \tau} = S\,n(d_1)\frac{\partial d_1}{\partial \tau} - Ke^{-r\tau}\left(-r\,N(d_2) + n(d_2)\frac{\partial d_2}{\partial \tau}\right). Using the core identity Sn(d1)=Kerτn(d2)S\,n(d_1) = Ke^{-r\tau}n(d_2), the d/τ\partial d/\partial\tau terms cancel, leaving: Θcall=Cτ=Sn(d1)σ2τrKerτN(d2).\Theta_{\mathrm{call}} = -\frac{\partial C}{\partial \tau} = -\frac{S\,n(d_1)\sigma}{2\sqrt{\tau}} - rKe^{-r\tau}N(d_2).

Sign conventions:

  • Θcall<0\Theta_{\mathrm{call}} < 0 always (for r>0r > 0): long calls lose time value.
  • Θput<0\Theta_{\mathrm{put}} < 0 for near-the-money puts; can be positive for deep ITM puts (reflecting that a very deep ITM put is almost equivalent to a bond position with positive carry).
  • Θ|\Theta| increases as τ0\tau \to 0: time decay accelerates near expiry.

The BS PDE as a P&L Identity

The Black-Scholes PDE can be written in Greek notation: Θ+12σ2S2Γ+rSΔrC=0.\Theta + \frac{1}{2}\sigma^2 S^2 \Gamma + rS\Delta - rC = 0.

Rearranged: Θ+12σ2S2Γ=r(CSΔ)=r(KerτN(d2))<0\Theta + \frac{1}{2}\sigma^2 S^2 \Gamma = r(C - S\Delta) = r \cdot (-Ke^{-r\tau}N(d_2)) < 0.

Interpretation: theta and gamma income sum to the risk-free financing cost of the position. A long gamma position (positive Γ\Gamma) must pay theta to hold it — this is the gamma-theta tradeoff.


Vega: ν=C/σ\nu = \partial C / \partial \sigma

ν=Sτn(d1).\nu = S\sqrt{\tau}\,n(d_1).

Derivation: Cσ=Sn(d1)d1σKerτn(d2)d2σ.\frac{\partial C}{\partial \sigma} = S\,n(d_1)\frac{\partial d_1}{\partial \sigma} - Ke^{-r\tau}n(d_2)\frac{\partial d_2}{\partial \sigma}. Using d1/σ=d2/σ\partial d_1/\partial\sigma = -d_2/\sigma and d2/σ=d1/στ=d2/στ\partial d_2/\partial\sigma = \partial d_1/\partial\sigma - \sqrt{\tau} = -d_2/\sigma - \sqrt{\tau}: Cσ=Sn(d1)(d2σ)Kerτn(d2)(d2στ).\frac{\partial C}{\partial \sigma} = S\,n(d_1)\left(-\frac{d_2}{\sigma}\right) - Ke^{-r\tau}n(d_2)\left(-\frac{d_2}{\sigma} - \sqrt{\tau}\right). Applying the core identity to the d2/σ-d_2/\sigma terms (they cancel), the surviving term is Kerτn(d2)τ=Sn(d1)τKe^{-r\tau}n(d_2)\sqrt{\tau} = S\,n(d_1)\sqrt{\tau}. Hence: ν=Sτn(d1).\nu = S\sqrt{\tau}\,n(d_1).

Properties:

  • νcall=νput\nu_{\mathrm{call}} = \nu_{\mathrm{put}}: same for both (from put-call parity).
  • ν>0\nu > 0: long options are long vega — they benefit from rising implied vol.
  • ν\nu is maximised at-the-money and decreases for deep ITM/OTM options.
  • ν\nu grows with τ\sqrt{\tau}: longer-dated options have more vega per unit spot.

Note on units: vega is usually quoted in price change per 1 vol point (i.e., per 0.01 change in σ\sigma in decimal, or equivalently per 1% change in annualised vol). Check units carefully — mixing percent and decimal conventions is a frequent source of desk errors.


Rho: ϱ=C/r\varrho = \partial C / \partial r

ϱcall=KτerτN(d2),ϱput=KτerτN(d2).\varrho_{\mathrm{call}} = K\tau e^{-r\tau}N(d_2), \qquad \varrho_{\mathrm{put}} = -K\tau e^{-r\tau}N(-d_2).

Derivation: Cr=Sn(d1)d1r+KτerτN(d2)Kerτn(d2)d2r.\frac{\partial C}{\partial r} = S\,n(d_1)\frac{\partial d_1}{\partial r} + K\tau e^{-r\tau}N(d_2) - Ke^{-r\tau}n(d_2)\frac{\partial d_2}{\partial r}. Since d1/r=d2/r=τ/σ\partial d_1/\partial r = \partial d_2/\partial r = \sqrt{\tau}/\sigma, the nn terms cancel by the core identity, leaving ϱcall=KτerτN(d2)\varrho_{\mathrm{call}} = K\tau e^{-r\tau}N(d_2).

Materiality: rho is small for short-dated equity options (where τ1\tau \ll 1) but material for long-dated options and for interest rate derivatives. For equity exotics with τ=5\tau = 5 years and r=4%r = 4\%, rho can dominate the Greeks P&L on an absolute basis.


Higher-Order Greeks

Beyond the five first-order Greeks, second-order cross-sensitivities are essential for managing options books, particularly for barrier options and structured products.

Vanna: 2C/Sσ\partial^2 C / \partial S \partial \sigma

Vanna=2CSσ=Δσ=νS=n(d1)d2σ.\mathrm{Vanna} = \frac{\partial^2 C}{\partial S \partial \sigma} = \frac{\partial \Delta}{\partial \sigma} = \frac{\partial \nu}{\partial S} = -\frac{n(d_1)\,d_2}{\sigma}.

Derivation: Vanna=(N(d1))/σ=n(d1)d1/σ=n(d1)(d2/σ)\mathrm{Vanna} = \partial(N(d_1))/\partial\sigma = n(d_1)\cdot\partial d_1/\partial\sigma = n(d_1)\cdot(-d_2/\sigma).

Interpretation: vanna measures how the delta changes as implied vol moves. If d2<0d_2 < 0 (OTM call), vanna is positive: as vol rises, the call's delta increases (the option becomes more likely to expire ITM). For a delta-hedged book, a simultaneous move in SS and σ\sigma creates a vanna P&L of vanna ΔSΔσ\cdot \Delta S \cdot \Delta\sigma.

Barrier options: for a knock-out barrier near the current spot, vanna can be very large — a small spot move near the barrier changes the option's delta dramatically, and simultaneously, implied vol changes amplify this.

Volga: 2C/σ2\partial^2 C / \partial \sigma^2

Volga=2Cσ2=νσ=νd1d2σ.\mathrm{Volga} = \frac{\partial^2 C}{\partial \sigma^2} = \frac{\partial \nu}{\partial \sigma} = \frac{\nu \cdot d_1 d_2}{\sigma}.

Derivation: Volga=(Sτn(d1))/σ=Sτ(d1n(d1))d1/σ=Sτn(d1)d1d2/σ=νd1d2/σ\mathrm{Volga} = \partial(S\sqrt{\tau}\,n(d_1))/\partial\sigma = S\sqrt{\tau}\,(-d_1\,n(d_1))\cdot\partial d_1/\partial\sigma = S\sqrt{\tau}\,n(d_1)\cdot d_1 d_2/\sigma = \nu \cdot d_1 d_2/\sigma.

Interpretation: volga is the convexity of the option price in volatility. Long volga positions benefit from large vol moves in either direction. OTM options (where d1,d2|d_1|, |d_2| are large and of the same sign) have higher volga than ATM options (where d1,d20d_1, d_2 \approx 0).

Charm: 2C/St=Δ/t\partial^2 C / \partial S \partial t = \partial \Delta / \partial t

Charm (delta decay) measures how the delta changes with time. Near expiry, ATM options have rapidly changing delta — a delta-hedged position must be rebalanced more frequently.


P&L Decomposition

A complete P&L decomposition for a single option position over a time interval [t,t+dt][t, t+dt], including both first-order and second-order contributions, uses Itô's lemma:

dCΔdS+12Γ(dS)2+Θdt+νdσ+VannadSdσ+12Volga(dσ)2.dC \approx \Delta\,dS + \frac{1}{2}\Gamma\,(dS)^2 + \Theta\,dt + \nu\,d\sigma + \mathrm{Vanna}\,dS\,d\sigma + \frac{1}{2}\mathrm{Volga}\,(d\sigma)^2.

For a delta-hedged portfolio (Π=CΔS\Pi = C - \Delta \cdot S, where the stock position is funded at rate rr):

d\Pi = \underbrace{\frac{1}{2}\Gamma\,(dS)^2}_{\text{gamma P&L}} + \underbrace{\Theta\,dt}_{\text{time decay}} + \underbrace{\nu\,d\sigma}_{\text{vega P&L}} + \underbrace{\mathrm{Vanna}\,dS\,d\sigma + \frac{1}{2}\mathrm{Volga}\,(d\sigma)^2}_{\text{cross-Greek P&L}}.

Under Black-Scholes assumptions: (dS)2=σ2S2dt(dS)^2 = \sigma^2 S^2\,dt, dσ=0d\sigma = 0. The portfolio earns zero P&L on average — the gamma income exactly offsets theta decay. In practice:

dΠactual=12S2Γ(σR2σ^2)dt+νdσ+VannadSdσ+,d\Pi_{\mathrm{actual}} = \frac{1}{2}S^2\Gamma\left(\sigma_R^2 - \hat{\sigma}^2\right)dt + \nu\,d\sigma + \mathrm{Vanna}\,dS\,d\sigma + \cdots,

where σ^\hat{\sigma} is implied vol and σR\sigma_R is realised vol. The first term is the gamma trading P&L (earn when realised vol exceeds implied), the second is the vega P&L (earn when implied vol moves in your favour), and the cross-Greek terms are the second-order vol P&L — significant for exotic options and barrier-heavy books.


Limitations

Constant vol assumption. Black-Scholes Greeks assume σ\sigma is constant. In practice, ΔBS\Delta_{\mathrm{BS}} with a flat-smile implied vol is not the correct hedge ratio — a sticky-strike or sticky-delta delta adjustment is needed, depending on which smile dynamics are assumed. The correct hedge ratio in a local vol model differs from N(d1)N(d_1).

Discrete hedging. Analytical Greeks are instantaneous rates. Real hedging occurs at discrete intervals, creating a replication error proportional to 12S2Γ((ΔS)2σ2S2Δt)\frac{1}{2}S^2\Gamma((\Delta S)^2 - \sigma^2 S^2\,\Delta t) per step.

Jump risk. Under Black-Scholes, Δ\Delta-hedging is exact (in the continuous limit). With jumps, a gap in SS creates an unhedgeable P&L of 12ΓJ2\frac{1}{2}\Gamma J^2 for a jump of size JJ. No hedge using only the underlying eliminates jump risk.

Model dependence. Greeks are model outputs. A position that is delta-neutral under Black-Scholes is not delta-neutral under Heston or SABR. The hedge ratio is model-dependent; using the wrong model produces systematic hedging errors.


Interview Angle

L1. State the five Greeks of a European call. Which is the same for calls and puts (by put-call parity)? What is the sign of each for a long call?

Γ and ν are identical for calls and puts (from 2(CP)/S2=0\partial^2(C-P)/\partial S^2 = 0 and (CP)/σ=0\partial(C-P)/\partial\sigma = 0, since CP=SKerτC - P = S - Ke^{-r\tau} contains neither σ\sigma nor S2S^2). For a long call: Δ > 0, Γ > 0, Θ < 0, ν > 0, ρ > 0.

L2. Derive Γ from first principles. Write the BS PDE in Greek notation and interpret the gamma-theta tradeoff. Show that a delta-hedged call earns P&L of 12S2Γ(σR2σ^2)dt\frac{1}{2}S^2\Gamma(\sigma_R^2 - \hat{\sigma}^2)dt when realised vol differs from implied.

Gamma derivation: Γ=N(d1)/S=n(d1)d1/S=n(d1)/(Sστ)\Gamma = \partial N(d_1)/\partial S = n(d_1)\cdot\partial d_1/\partial S = n(d_1)/(S\sigma\sqrt{\tau}) (using d1/S=1/(Sστ)\partial d_1/\partial S = 1/(S\sigma\sqrt{\tau})).

BS PDE: Θ+12σ^2S2Γ+rSΔrC=0\Theta + \frac{1}{2}\hat{\sigma}^2 S^2\Gamma + rS\Delta - rC = 0. The delta-hedged portfolio Π=CΔS\Pi = C - \Delta S earns dΠ=(Θ+12σR2S2Γ)dtd\Pi = (\Theta + \frac{1}{2}\sigma_R^2 S^2\Gamma)dt using actual vol. Subtracting the BS PDE identity evaluated at implied vol: dΠ=12S2Γ(σR2σ^2)dtd\Pi = \frac{1}{2}S^2\Gamma(\sigma_R^2 - \hat{\sigma}^2)dt.

L3. Derive vanna and volga analytically. Explain their role in the P&L of a barrier option book. Why does a knock-out barrier near the spot produce extreme vanna risk?

Vanna: 2C/Sσ=Δ/σ=n(d1)d1/σ=n(d1)d2/σ\partial^2C/\partial S\partial\sigma = \partial\Delta/\partial\sigma = n(d_1)\cdot\partial d_1/\partial\sigma = -n(d_1)d_2/\sigma. Also equals ν/S=(Sτn(d1))/S=τ(n(d1)Sd1n(d1)/(Sστ))=n(d1)τ(1d1/(στ))=n(d1)(στd1)/σ=n(d1)d2/σ\partial\nu/\partial S = \partial(S\sqrt{\tau}n(d_1))/\partial S = \sqrt{\tau}(n(d_1) - Sd_1n(d_1)/(S\sigma\sqrt{\tau})) = n(d_1)\sqrt{\tau}(1 - d_1/(σ\sqrt{\tau})) = n(d_1)(σ\sqrt{\tau} - d_1)/\sigma = -n(d_1)d_2/\sigma (consistent).

Volga: ν/σ=(Sτn(d1))/σ=Sτ(d1n(d1))(d2/σ)=νd1d2/σ\partial\nu/\partial\sigma = \partial(S\sqrt{\tau}n(d_1))/\partial\sigma = S\sqrt{\tau}(-d_1n(d_1))(-d_2/\sigma) = \nu d_1d_2/\sigma.

Knock-out barrier: as SBS \to B^- (spot approaching barrier from below), the knock-out call's value 0\to 0 discontinuously (the option is worth zero if knocked out). This creates a steep negative delta near the barrier (the delta hedge requires selling shares as spot rises toward the barrier). Delta changes rapidly with spot (Γ\Gamma is large near the barrier) and also with vol (vanna is large): a vol rise increases the probability of the spot reaching the barrier before expiry, further changing the delta. The combined spot-vol move creates large vanna P&L: vanna ΔSΔσ\cdot \Delta S \cdot \Delta\sigma. This is why barrier options require careful vanna and volga hedging — typically via OTM vanilla options whose vanna/volga profile partially offsets the barrier's exposure.

Read the theory? Verify your understanding.

Take the Quiz