Unit 7, lesson 53

Jump diffusion and incomplete hedging

How Poisson jumps enter a stock model under a pricing measure, why the drift between jumps is compensated by the expected jump, how jumps make the pricing generator nonlocal, and what jump residual a delta hedge leaves.

  • Video 53
  • Local volatility, stochastic volatility and jumps

A stock that can gap down by twenty percent overnight cannot be described by a diffusion, whose paths are continuous. Adding random jumps to a pricing model raises three questions. Which drift keeps the expected growth of the stock at the carry rate, which equation does an option price satisfy, and can a position in the stock still hedge the option? This lesson adds finite-activity Poisson jumps under a pricing measure, derives the compensated drift and the nonlocal pricing generator, and computes the exact residual that a delta hedge leaves when a jump occurs.

Key terms

Jump diffusion

Definition Under a chosen pricing measure Q\mathbb{Q}, let WQW^{\mathbb{Q}} be a Brownian motion, NN a Poisson process with constant intensity λ≥0\lambda \ge 0, and J1,J2,…J_1, J_2, \dots independent, identically distributed marks with Jk>0J_k > 0 and EQ[Jk]<∞\mathbb{E}^{\mathbb{Q}}[J_k] < \infty, where WQW^{\mathbb{Q}}, NN and the marks are independent. Let rr, qq and σ≥0\sigma \ge 0 be constants and write κJ=EQ[J−1]\kappa_J = \mathbb{E}^{\mathbb{Q}}[J - 1] for the mean proportional jump. The jump diffusion with proportional jumps is

dStSt−=(r−q−λκJ) dt+σ dWtQ+(JNt−1) dNt.\begin{aligned} \frac{dS_t}{S_{t-}} &= (r - q - \lambda\kappa_J)\, dt \\[0.4em] &\quad + \sigma \, dW^{\mathbb{Q}}_t \\[0.4em] &\quad + (J_{N_t} - 1)\, dN_t . \end{aligned}

Between jumps the price follows a geometric Brownian motion with drift r−q−λκJr - q - \lambda\kappa_J, and at the kkth jump time TkT_k it is multiplied by the kkth mark, STk=STk−JkS_{T_k} = S_{T_k-} J_k. The left limit St−S_{t-} appears because the size of a jump is proportional to the price just before it. The solution is

St=S0 e(r−q−λκJ)t×eσWtQ−σ2t/2×∏k=1NtJk,\begin{aligned} S_t &= S_0 \, e^{(r - q - \lambda\kappa_J)t} \\[0.4em] &\quad \times e^{\sigma W^{\mathbb{Q}}_t - \sigma^2 t/2} \\[0.4em] &\quad \times \prod_{k=1}^{N_t} J_k , \end{aligned}

with the empty product equal to 11 when Nt=0N_t = 0. It is positive because every mark is positive. A mark J≤0J \le 0 lies outside this multiplicative model.

In plain words The price moves continuously most of the time, and at rare random moments it is multiplied by a random factor. A factor below one is a drop and a factor above one is a rise.

Example The constructed values used throughout this lesson are λ=0.4\lambda = 0.4 per year, the fixed factor J=0.8J = 0.8, r=0.03r = 0.03, q=0q = 0 and S0=100S_0 = 100, with σ=0.2\sigma = 0.2 wherever a volatility is needed. Then κJ=−0.2\kappa_J = -0.2 and the drift between jumps is 0.03+0.08=0.110.03 + 0.08 = 0.11. Given nn jumps in the first year, the mean of S1S_1 is 100e0.11⋅0.8n100 e^{0.11} \cdot 0.8^n, about 111.63111.63 with no jump, 89.3089.30 with one jump and 71.4471.44 with two. Weighted with the Poisson probabilities of nn, these conditional means average to EQ[S1]=100e0.03≈103.05\mathbb{E}^{\mathbb{Q}}[S_1] = 100 e^{0.03} \approx 103.05.

Poisson jumps

Definition A Poisson process with intensity λ≥0\lambda \ge 0 is a counting process NN with N0=0N_0 = 0, right-continuous paths that increase only by jumps of size one, and independent increments with Nt−Ns∼Poisson⁡(λ(t−s))N_t - N_s \sim \operatorname{Poisson}(\lambda(t - s)) for s<ts < t, that is

Q(Nt−Ns=n)=e−λ(t−s)(λ(t−s))nn!\begin{aligned} &\mathbb{Q}(N_t - N_s = n) \\[0.4em] &\quad = e^{-\lambda(t - s)} \frac{(\lambda(t - s))^n}{n!} \end{aligned}

for n=0,1,2,…n = 0, 1, 2, \dots. The jump times are the times at which NN increases. Every bounded interval contains finitely many of them almost surely, which is what finite activity means. The probability generating function is

EQ[zNt]=eλt(z−1)\mathbb{E}^{\mathbb{Q}}\bigl[z^{N_t}\bigr] = e^{\lambda t (z - 1)}

for every real zz. Because the marks are independent of NN, conditioning on NtN_t gives

EQ[∏k=1NtJk]=EQ[(1+κJ)Nt]=eλtκJ.\begin{aligned} &\mathbb{E}^{\mathbb{Q}}\Bigl[\prod_{k=1}^{N_t} J_k\Bigr] \\[0.4em] &= \mathbb{E}^{\mathbb{Q}}\bigl[(1 + \kappa_J)^{N_t}\bigr] \\[0.4em] &= e^{\lambda t \kappa_J} . \end{aligned}

In plain words Jumps arrive at a constant average rate of λ\lambda per year, with no memory of when the last one occurred. Over a period of length tt the expected number of jumps is λt\lambda t.

Example With λ=0.4\lambda = 0.4 and t=1t = 1, the probability of no jump is e−0.4≈0.6703e^{-0.4} \approx 0.6703, of exactly one jump 0.4e−0.4≈0.26810.4 e^{-0.4} \approx 0.2681, and of two or more about 0.06160.0616. For the fixed factor J=0.8J = 0.8 the generating function gives EQ[0.8N1]=e0.4(0.8−1)=e−0.08≈0.9231\mathbb{E}^{\mathbb{Q}}[0.8^{N_1}] = e^{0.4(0.8 - 1)} = e^{-0.08} \approx 0.9231. The path with no jump is the most likely single outcome, but it is not the mean path, since the mean averages over every possible number of jumps.

Compensated jump process

Definition For a Poisson process with intensity λ\lambda, the compensated process

N~t=Nt−λt\widetilde{N}_t = N_t - \lambda t

is a martingale for any filtration to which NN is adapted and whose past is independent of the future increments of NN, because EQ[Nt−Ns∣Fs]=λ(t−s)\mathbb{E}^{\mathbb{Q}}[N_t - N_s \mid \mathcal{F}_s] = \lambda(t - s). Its variance is Var⁡(N~t)=λt\operatorname{Var}(\widetilde{N}_t) = \lambda t. For marks with EQ∣J−1∣<∞\mathbb{E}^{\mathbb{Q}}\lvert J - 1 \rvert < \infty, the compensated jump sum

Mt=∑k=1Nt(Jk−1)−λκJt\begin{aligned} M_t &= \sum_{k=1}^{N_t} (J_k - 1) - \lambda\kappa_J t \end{aligned}

is a martingale too, since the sum has independent increments and mean λκJt\lambda\kappa_J t. The compensator λκJt\lambda\kappa_J t is a deterministic drift that removes the mean of the jump sum. It does not replace the realized marks by their mean, and the jump times and the law of the marks stay as they were. In the language of random measures, the compensator of the jumps is λ dt F(dy)\lambda \, dt \, F(dy), where FF is the law of a mark.

In plain words A raw jump count only goes up, so when λ>0\lambda > 0 its increments have positive mean. Subtracting the expected growth λt\lambda t leaves a process with mean-zero increments that still jumps at the same random times by the same amounts.

Example With λ=0.4\lambda = 0.4, the variable N~1=N1−0.4\widetilde{N}_1 = N_1 - 0.4 equals −0.4-0.4 with probability about 0.67030.6703, equals 0.60.6 with probability about 0.26810.2681, and is 1.61.6 or more with probability about 0.06160.0616. Its mean is 00 and its variance 0.40.4. Marks equal to 0.70.7 or 0.90.9 with probability one half each have κJ=−0.2\kappa_J = -0.2, the same as the fixed factor 0.80.8, so the two jump sums have the same compensator −0.08 t-0.08\,t but different distributions. Replacing each mark by 0.80.8 would lose that difference.

Drift compensation

Definition Drift compensation is the choice of the drift between jumps that makes the discounted stock with reinvested dividends, e−(r−q)tSte^{-(r - q)t} S_t, a Q\mathbb{Q}-martingale. In the jump diffusion that drift is r−q−λκJr - q - \lambda\kappa_J. From the explicit solution, with τ=T−t\tau = T - t,

EQ[ST∣Ft]=St e(r−q−λκJ)τeλκJτ=St e(r−q)τ,\begin{aligned} &\mathbb{E}^{\mathbb{Q}}[S_T \mid \mathcal{F}_t] \\[0.4em] &= S_t \, e^{(r - q - \lambda\kappa_J)\tau} e^{\lambda\kappa_J\tau} \\[0.4em] &= S_t \, e^{(r - q)\tau}, \end{aligned}

because the Brownian factor eσ(WTQ−WtQ)−σ2τ/2e^{\sigma(W^{\mathbb{Q}}_T - W^{\mathbb{Q}}_t) - \sigma^2\tau/2} has mean one and the product of the marks over (t,T](t, T] has mean eλκJτe^{\lambda\kappa_J\tau}. Written with the compensated jump sum MM, the stock equation is

dStSt−=(r−q) dt+σ dWtQ+dMt.\begin{aligned} \frac{dS_t}{S_{t-}} &= (r - q)\, dt \\[0.4em] &\quad + \sigma \, dW^{\mathbb{Q}}_t + dM_t . \end{aligned}

The condition fixes only the product λκJ\lambda\kappa_J, not λ\lambda and the law of JJ separately. Suppose the physical model has the same form with drift μP\mu^{\mathbb{P}} between jumps, σ>0\sigma > 0, intensity λP>0\lambda^{\mathbb{P}} > 0 and the same law of the marks. Then every intensity λQ>0\lambda^{\mathbb{Q}} > 0 gives an equivalent martingale measure, with density

dQdP∣FT=e−ηWTP−η2T/2×(λQλP)NTe(λP−λQ)T,\begin{aligned} \frac{d\mathbb{Q}}{d\mathbb{P}}\Big\vert_{\mathcal{F}_T} &= e^{-\eta W^{\mathbb{P}}_T - \eta^2 T/2} \\[0.4em] &\quad \times \Bigl(\frac{\lambda^{\mathbb{Q}}}{\lambda^{\mathbb{P}}}\Bigr)^{N_T} e^{(\lambda^{\mathbb{P}} - \lambda^{\mathbb{Q}})T}, \end{aligned}

where WtQ=WtP+ηtW^{\mathbb{Q}}_t = W^{\mathbb{P}}_t + \eta t is the Q\mathbb{Q}-Brownian motion and η=(μP−r+q+λQκJ)/σ\eta = (\mu^{\mathbb{P}} - r + q + \lambda^{\mathbb{Q}}\kappa_J)/\sigma. This construction is a standard Girsanov result for jump diffusions, stated here without proof. The choice λQ=0\lambda^{\mathbb{Q}} = 0 is excluded, since a measure under which no jump occurs is not equivalent to one under which jumps have positive probability.

In plain words Jumps with a negative mean pull the expected price down, so the drift between jumps must be higher by the same amount for the total expected growth to equal the carry rate. The jump term is offset in the mean, not removed path by path.

Example With λ=0.4\lambda = 0.4, J=0.8J = 0.8, r=0.03r = 0.03 and q=0q = 0, the drift between jumps is 0.03−0.4(−0.2)=0.110.03 - 0.4(-0.2) = 0.11, well above the carry rate, and EQ[S1]=100e0.11e−0.08=100e0.03≈103.05\mathbb{E}^{\mathbb{Q}}[S_1] = 100 e^{0.11} e^{-0.08} = 100 e^{0.03} \approx 103.05. A simulation of two million paths with σ=0.2\sigma = 0.2 gives 103.03103.03 with a standard error of 0.0170.017. With the same factor 0.80.8 but intensity λQ=0.8\lambda^{\mathbb{Q}} = 0.8, the drift between jumps becomes 0.190.19 and the forward price is unchanged. The European put with strike 100100 and maturity 11 year, with σ=0.2\sigma = 0.2, is worth about 8.0068.006 at λQ=0.4\lambda^{\mathbb{Q}} = 0.4 and about 9.3649.364 at λQ=0.8\lambda^{\mathbb{Q}} = 0.8. Each value comes from an equivalent martingale measure, so neither creates an arbitrage with the stock and the money-market account, and the stock alone does not decide between them.

Nonlocal pricing generator

Definition For a twice continuously differentiable function ff on (0,∞)(0, \infty) for which the expectation below is finite, the generator of the jump diffusion under Q\mathbb{Q} is

Lf(s)=(r−q)sf′(s)+12σ2s2f′′(s)+λ EQ[f(sJ)−f(s)−s(J−1)f′(s)].\begin{aligned} \mathcal{L} f(s) &= (r - q) s f'(s) \\[0.4em] &\quad + \tfrac12 \sigma^2 s^2 f''(s) \\[0.4em] &\quad + \lambda \, \mathbb{E}^{\mathbb{Q}}\bigl[ f(sJ) - f(s) \\[0.4em] &\qquad - s(J - 1) f'(s) \bigr]. \end{aligned}

Expanding the last term gives the equivalent form (r−q−λκJ)sf′+12σ2s2f′′+λ EQ[f(sJ)−f(s)](r - q - \lambda\kappa_J) s f' + \tfrac12\sigma^2 s^2 f'' + \lambda\,\mathbb{E}^{\mathbb{Q}}[f(sJ) - f(s)] with the raw drift between jumps. The generator comes from Itô’s formula for jump processes, in which f(St)−f(S0)f(S_t) - f(S_0) is the sum of the usual diffusion terms and of the jump changes f(STk−Jk)−f(STk−)f(S_{T_k-} J_k) - f(S_{T_k-}), and the jump changes have the compensator λ∫0tEQ[f(Su−J)−f(Su−)] du\lambda \int_0^t \mathbb{E}^{\mathbb{Q}}[f(S_{u-} J) - f(S_{u-})] \, du, where the expectation is over JJ with Su−S_{u-} held fixed. The jump term evaluates ff at the point sJsJ, away from ss, which is what nonlocal means. For a European claim with payoff g(ST)g(S_T) the pricing equation is the partial integro-differential equation

∂tV+LV−rV=0,V(T,s)=g(s),\begin{gathered} \partial_t V + \mathcal{L} V - rV = 0, \\[0.4em] V(T, s) = g(s), \end{gathered}

with L\mathcal{L} acting on the variable ss. If it has a solution that is continuously differentiable in tt and twice in ss with suitable growth, and the payoff and the marks are integrable enough, then V(t,s)=EQ[e−r(T−t)g(ST)∣St=s]V(t, s) = \mathbb{E}^{\mathbb{Q}}[e^{-r(T - t)} g(S_T) \mid S_t = s], the Feynman-Kac representation for jump diffusions, stated here without proof. The price depends on the chosen Q\mathbb{Q} through λ\lambda and the law of JJ.

In plain words In a diffusion, the expected change of a value over the next instant depends on its slope and curvature at the current price. With jumps it also depends on the values at every price a jump can reach, weighted by the jump rate and the law of the jump size.

Example For f(s)=ln⁡sf(s) = \ln s the generator gives the expected growth rate of the log price,

Lln⁡s=r−q−12σ2+λ EQ[ln⁡J−(J−1)].\begin{aligned} \mathcal{L} \ln s &= r - q - \tfrac12\sigma^2 \\[0.4em] &\quad + \lambda \, \mathbb{E}^{\mathbb{Q}}\bigl[\ln J - (J - 1)\bigr]. \end{aligned}

With r=0.03r = 0.03, q=0q = 0, σ=0.2\sigma = 0.2, λ=0.4\lambda = 0.4 and J=0.8J = 0.8 this is 0.01+0.4(ln⁡0.8+0.2)≈0.01−0.00926=0.000740.01 + 0.4(\ln 0.8 + 0.2) \approx 0.01 - 0.00926 = 0.00074 per year, in agreement with EQ[ln⁡S1]−ln⁡S0=0.11−0.02+0.4ln⁡0.8\mathbb{E}^{\mathbb{Q}}[\ln S_1] - \ln S_0 = 0.11 - 0.02 + 0.4 \ln 0.8 from the explicit solution. Replacing the jump term by its local second-order Taylor approximation −12λ(J−1)2=−0.008-\tfrac12\lambda(J - 1)^2 = -0.008 would give 0.0020.002 instead, because a jump of twenty percent is not small.

Jump residual

Definition Let VV be an option value at a fixed time, differentiable at the current price ss, and hold one option together with −V′(s)-V'(s) shares, the delta hedge chosen just before a jump. If the price jumps from ss to sJsJ, the jump residual is the change in value of this hedged position,

RJ(s)=V(sJ)−V(s)−V′(s) s(J−1).\begin{aligned} R_J(s) &= V(sJ) - V(s) \\[0.4em] &\quad - V'(s)\, s (J - 1) . \end{aligned}

It is the exact change, not an approximation. If VV is twice continuously differentiable near ss, Taylor’s theorem gives RJ(s)≈12V′′(s)s2(J−1)2R_J(s) \approx \tfrac12 V''(s) s^2 (J - 1)^2 for small jumps, and for a convex VV the residual is nonnegative for every JJ. The residual vanishes for every mark in the support of the jump law only if all the jumped values V(sJ)V(sJ) lie on the tangent line of VV at ss. With σ>0\sigma > 0 the hedge ratio is already pinned to V′(s)V'(s) by the Brownian term, so unless every jumped value lies on that tangent line, the stock and the money-market account cannot hedge both risks, and the market made of these two assets is incomplete. Further traded options can span some of the jump risk, which changes the market. A degenerate model can differ. With σ=0\sigma = 0 and a single jump factor, the chord slope (V(sJ)−V(s))/(s(J−1))(V(sJ) - V(s))/(s(J - 1)) removes the jump risk exactly.

In plain words A delta hedge follows the option along its tangent line at the current price. A jump moves the price far along the curve in one step, and the gap between the curve and the tangent line at the new price is what the hedge misses.

Example Take the constructed value curve V(s)=s2/100V(s) = s^2/100 at s=100s = 100, where the delta is 22 shares. A jump with J=0.8J = 0.8 changes the option value from 100100 to 6464, a change of −36-36, while the short position of 22 shares gains 2⋅20=402 \cdot 20 = 40. The residual is

R0.8=−36+40=4,R_{0.8} = -36 + 40 = 4 ,

and for this curve RJ=100(J−1)2R_J = 100 (J - 1)^2 at s=100s = 100, so an upward jump with J=1.1J = 1.1 leaves 11, also positive. The Taylor formula 12V′′(s)s2(J−1)2\tfrac12 V''(s) s^2 (J - 1)^2 is exact here only because VV is quadratic.

Check yourself

  1. With λ=0.4\lambda = 0.4 per year, find the probability of no jump and of at least one jump during half a year, and the mean and the variance of N~0.5\widetilde{N}_{0.5}. For the fixed factor J=0.8J = 0.8, what is EQ[JN0.5]\mathbb{E}^{\mathbb{Q}}[J^{N_{0.5}}]?
Show the answer

The count N0.5N_{0.5} is Poisson with mean λt=0.2\lambda t = 0.2. The probability of no jump is e−0.2≈0.8187e^{-0.2} \approx 0.8187, so the probability of at least one is about 0.18130.1813. The compensated count N~0.5=N0.5−0.2\widetilde{N}_{0.5} = N_{0.5} - 0.2 has mean 00 and variance 0.20.2. The generating function gives

EQ[0.8N0.5]=e0.2(0.8−1)=e−0.04,\begin{aligned} \mathbb{E}^{\mathbb{Q}}\bigl[0.8^{N_{0.5}}\bigr] &= e^{0.2(0.8 - 1)} \\[0.4em] &= e^{-0.04}, \end{aligned}

which is about 0.96080.9608.

  1. A classmate writes the drift between jumps as r−q+λκJr - q + \lambda\kappa_J. With the constructed values r=0.03r = 0.03, q=0q = 0, λ=0.4\lambda = 0.4, J=0.8J = 0.8 and S0=100S_0 = 100, what is EQ[S1]\mathbb{E}^{\mathbb{Q}}[S_1] under that drift, and what goes wrong?
Show the answer

The classmate’s drift is 0.03−0.08=−0.050.03 - 0.08 = -0.05, and the jumps contribute the factor eλκJ=e−0.08e^{\lambda\kappa_J} = e^{-0.08} to the mean, so

EQ[S1]=100e−0.05e−0.08=100e−0.13≈87.81,\begin{aligned} \mathbb{E}^{\mathbb{Q}}[S_1] &= 100 e^{-0.05} e^{-0.08} \\[0.4em] &= 100 e^{-0.13} \approx 87.81, \end{aligned}

instead of 100e0.03≈103.05100 e^{0.03} \approx 103.05. The jump term is counted twice with the same sign instead of being offset. The discounted stock e−0.03tSte^{-0.03 t} S_t would then have mean 100e−0.16t100 e^{-0.16 t}, so it would not be a martingale and Q\mathbb{Q} would not be a pricing measure for the stock. The compensator must be subtracted, which gives the drift 0.110.11.

  1. Use the value curve V(s)=s2/100V(s) = s^2/100 at s=100s = 100. With σ=0\sigma = 0 and the single jump factor 0.80.8, which hedge ratio removes the jump risk? Why must the hedge use 22 shares when σ>0\sigma > 0, and what is then the residual of a jump? Finally, with σ=0\sigma = 0 and marks equal to 0.70.7 or 0.90.9 with probability one half each, is there one hedge ratio that removes both jumps, and what is the expected residual of the delta hedge per jump?
Show the answer

With one jump factor the chord slope removes the jump,

V(80)−V(100)80−100=64−100−20=1.8\begin{aligned} &\frac{V(80) - V(100)}{80 - 100} \\[0.4em] &\quad = \frac{64 - 100}{-20} = 1.8 \end{aligned}

shares, and with σ=0\sigma = 0 there is no other risk to match at that instant. With σ>0\sigma > 0 the Brownian shock moves the option by V′(100)=2V'(100) = 2 per unit of stock move to first order, so only 22 shares cancel it, and a jump with factor 0.80.8 then leaves the residual 64−100+40=464 - 100 + 40 = 4. With the two marks the chord slopes are (49−100)/(70−100)=1.7(49 - 100)/(70 - 100) = 1.7 and (81−100)/(90−100)=1.9(81 - 100)/(90 - 100) = 1.9, which differ, so no single hedge ratio removes both jumps even without a Brownian term. The delta hedge leaves 100(0.3)2=9100(0.3)^2 = 9 or 100(0.1)2=1100(0.1)^2 = 1, with mean 55. Replacing each mark by its mean 0.80.8 would give 44 and understate the expected residual, which is

100 E[(J−1)2]=100(Var⁡(J)+κJ2)=100(0.01+0.04)=5.\begin{aligned} &100\,\mathbb{E}\bigl[(J - 1)^2\bigr] \\[0.4em] &= 100\bigl(\operatorname{Var}(J) + \kappa_J^2\bigr) \\[0.4em] &= 100(0.01 + 0.04) = 5 . \end{aligned}

How this lesson connects