Forward Equations for Option Prices in Semimartingale Models
View metadata, citation and similar papers at core.ac.uk brought to you by CORE
provided by Research Papers in Economics
Forward equations for option prices in semimartingale models
Amel Bentata and Rama Cont
Laboratoire de Probabilit´eset Mod`elesAl´eatoires, CNRS - Universit´ede Paris VI. and IEOR Department, Columbia University, New York.
First version: Nov 2009. Revised: March 2011.
Abstract We derive a forward partial integro-differential equation for prices of call options in a model where the dynamics of the underlying asset under the pricing measure is described by a -possibly discontinuous- semimartin- gale. This result generalizes Dupire’s forward equation to a large class of non-Markovian models with jumps. MSC Classification Numbers: 60H30 91G20 35S10 91G80
Contents
1 Forward PIDEs for call options 3
hal-00445641, version 3 - 21 Apr 2011 1.1 General formulation of the forward equation ...... 3 1.2 Derivation of the forward equation ...... 6 1.3 Uniqueness of solutions of the forward PIDE ...... 9
2 Examples 15 2.1 Itˆoprocesses ...... 15 2.2 Markovian jump-diffusion models ...... 16 2.3 Pure jump processes ...... 17 2.4 Time changed L´evyprocesses ...... 18 2.5 Index options in a multivariate jump-diffusion model ...... 20 2.6 Forward equations for CDO pricing ...... 26
1 Since the seminal work of Black, Scholes and Merton [7, 30] partial differen- tial equations (PDE) have been used as a way of characterizing and efficiently computing option prices. In the Black-Scholes-Merton model and various exten- sions of this model which retain the Markov property of the risk factors, option prices can be characterized in terms of solutions to a backward PDE, whose variables are time (to maturity) and the value of the underlying asset. The use of backward PDEs for option pricing has been extended to cover options with path-dependent and early exercise features, as well as to multifactor models (see e.g. [1]). When the underlying asset exhibit jumps, option prices can be computed by solving an analogous partial integro-differential equation (PIDE) [2, 14]. A second important step was taken by Dupire [15, 16, 18] who showed that when the underlying asset is assumed to follow a diffusion process
dSt = St(t, St)dWt
prices of call options (at a given date t0) solve a forward PDE in the strike and maturity variables: ∂C ∂C K2(T,K)2 ∂2C t0 (T,K) = −r(T )K t0 (T,K) + t0 (T,K) ∂T ∂K 2 ∂K2
on [t0, ∞[×]0, ∞[ with the initial condition: ∀K > 0 Ct0 (t0,K) = (St0 − K)+. This forward equation allows to price call options with various strikes and ma- turities on the same underlying asset, by solving a single partial differential equation. Dupire’s forward equation also provides useful insights into the in- verse problem of calibrating diffusion models to observed call and put option prices [6]. Given the theoretical and computational usefulness of the forward equation, there have been various attempts to extend Dupire’s forward equation to other types of options and processes, most notably to Markov processes with jumps [2, 10, 12, 26, 9]. Most of these constructions use the Markov property of the hal-00445641, version 3 - 21 Apr 2011 underlying process in a crucial way (see however [27]). As noted by Dupire [17], the forward PDE holds in a more general context than the backward PDE: even if the (risk-neutral) dynamics of the underlying asset is not necessarily Markovian, but described by a continuous Brownian martingale dSt = SttdWt then call options still verify a forward PDE where the diffusion coefficient is given by the local (or effective) volatility function (t, S) given by q 2 (t, S) = E[t ∣St = S]
This method is linked to the “Markovian projection” problem: the construction of a Markov process which mimicks the marginal distributions of a martingale [5, 23, 29]. Such “mimicking processes” provide a method to extend the Dupire equation to non-Markovian settings.
2 We show in this work that the forward equation for call prices holds in a more general setting, where the dynamics of the underlying asset is described by a - possibly discontinuous - semimartingale. Our parametrization of the price dynamics is general, allows for stochastic volatility and does not assume jumps to be independent or driven by a L´evyprocess, although it includes these cases. Also, our derivation does not require ellipticity or non-degeneracy of the diffusion coefficient. The result is thus applicable to various stochastic volatility models with jumps, pure jump models and point process models used in equity and credit risk modeling. Our result extends the forward equation from the original diffusion setting of Dupire [16] to various examples of non-Markovian and/or discontinuous pro- cesses and implies previous derivations of forward equations [2, 10, 9, 12, 16, 17, 26, 28] as special cases. Section 2 gives examples of forward PIDEs obtained in various settings: time-changed L´evyprocesses, local L´evymodels and point processes used in portfolio default risk modeling. In the case where the under- lying risk factor follows, an Itˆoprocess or a Markovian jump-diffusion driven by a L´evyprocess, we retrieve previously known forms of the forward equation. In this case, our approach gives a rigorous derivation of these results under precise assumptions in a unified framework. In some cases, such as index options (Sec. 2.5) or CDO expected tranche notionals (Sec. 2.6), our method leads to a new, more general form of the forward equation valid for a larger class of models than previously studied [3, 12, 34]. The forward equation for call options is a PIDE in one (spatial) dimension, regardless of the number of factor driving the underlying asset. It may thus be used as a method for reducing the dimension of the problem. The case of index options (Section 2.5) in a multivariate jump-diffusion model illustrates how the forward equation projects a high dimensional pricing problem into a one-dimensional state equation.
1 Forward PIDEs for call options hal-00445641, version 3 - 21 Apr 2011 1.1 General formulation of the forward equation Consider a (strictly positive) price process S whose dynamics under the pricing measure ℙ is given by a stochastic volatility model with jumps: Z T Z T Z T Z +∞ y ˜ ST = S0 + r(t)St− dt + St− tdWt + St− (e − 1)M(dt dy) (1) 0 0 0 −∞
where r(t) > 0 represents a (deterministic) bounded discount rate, t the (ran- dom) volatility process and M is an integer-valued random measure with com- pensator (dt dy; !) = m(t, dy, !) dt, representing jumps in the log-price, and M˜ = M − is the compensated random measure associated to M (see [13] for further background). Both the volatility t and m(t, dy), which represents the intensity of jumps of size y at time t, are allowed to be stochastic. In particular,
3 we do not assume the jumps to be driven by a L´evyprocess or a process with independent increments. We assume the following conditions:
Assumption 1 (Full support). For every t, supp(St) = [0, ∞[. Assumption 2 (Integrability condition). " !# 1 Z T Z T Z ∀T > 0, exp 2 dt + dt (ey − 1)2m(t, dy) < ∞ (H) E 2 t 0 0 ℝ
The value Ct0 (T,K) at time t0 of a call option with expiry T > t0 and strike K > 0 is given by
− R T r(t) dt t0 ℙ Ct0 (T,K) = e E [max(ST − K, 0)∣ℱt0 ] (2) As argued in Section 1.2, under Assumption (H), the expectation in (2) is finite. Our main result is the following:
Theorem 1 (Forward PIDE for call options). Let t be the exponential double tail of the compensator m(t, dy) ( R z dx ex R x m(t, du) z < 0 (z) = −∞ −∞ (3) t R +∞ x R ∞ z dx e x m(t, du) z > 0
and define, for t ∈ [t0, ∞[, z > 0,
( p 2 (t, z) = [ ∣S − = z]; E t t (4) t,y(z) = E [ t (z) ∣St− = y]
Under assumption (H), the call option price (T,K) 7→ Ct0 (T,K), as a function of maturity and strike, is a solution (in the sense of distributions) of the partial integro-differential equation: hal-00445641, version 3 - 21 Apr 2011 2 2 2 Z +∞ 2 ∂Ct0 ∂Ct0 K (T,K) ∂ Ct0 ∂ Ct0 K = −r(T )K + 2 + y 2 (T, dy) T,y ln ∂T ∂K 2 ∂K 0 ∂K y (5)
on [t0, ∞[×]0, ∞[ with the initial condition: ∀K > 0 Ct0 (t0,K) = (St0 − K)+.
Remark 1. Recall that f :[t0, ∞[×]0, ∞[7→ ℝ is a solution of (5) in the sense ∞ of distributions if for any test function ' ∈ C0 (]0, ∞[, ℝ) and for any T ≥ t0, Z ∞ ∂f ∂f K2(T,K)2 ∂2f Z +∞ ∂2f K dK'(K) − − r(T )K + 2 + y 2 (T, dy) T,y ln = 0 0 ∂T ∂K 2 ∂K 0 ∂K y ∞ where C0 (]0, ∞[, ℝ) is the set of infinitely differentiable functions with compact support in ]0, ∞[. This notion of generalized solution allows to separate the discussion of existence of solutions from the discussion of their regularity (which may be delicate, see [14]).
4 Remark 2. The discounted asset price
− R T r(t)dt SˆT = e 0 ST , is the stochastic exponential of the martingale U defined by Z T Z T Z y UT = t dWt + (e − 1)M˜ (dt dy). 0 0 Under assumption (H), we have 1 ∀T > 0, exp ⟨U, U⟩d + ⟨U, U⟩c < ∞ E 2 T T ˆ and [32, Theorem 9] implies that (ST ) is a ℙ-martingale. The form of the integral term in (5) may seem different from the integral term appearing in backward PIDEs [14, 25]. The following lemma expresses T,y(z) in a more familiar formin terms of call payoffs: Lemma 1. Let n(t, dz, y, !) dt be a random measure on [0,T ]×ℝ×ℝ+ verifying Z ∞ ∀t ∈ [0,T ], (ez ∧ ∣z∣2)n(t, dz, y, !) < ∞ a.s. −∞
Then the exponential double tail t,y(z) of n, defined as ( R z dx ex R x n(t, du, y) z < 0 (z) = −∞ −∞ (6) t,y R +∞ x R ∞ z dx e x n(t, du, y) z > 0 verifies Z K [(yez − K)+ − ez(y − K)+ − K(ez − 1)1 ]n(t, dz, y) = y ln {y>K} t,y y ℝ hal-00445641, version 3 - 21 Apr 2011 Proof. Let K, T > 0. Then: Z z + z + z [(ye − K) − e (y − K) − K(e − 1)1{y>K}]n(t, dz, y) Zℝ z z z = [(ye − K)1 K − e (y − K)1{y>K} − K(e − 1)1{y>K}]n(t, dz, y) {z>ln ( y )} Zℝ z z = [(ye − K)1 K + (K − ye )1{y>K}]n(t, dz, y). {z>ln ( y )} ℝ ∙ If K ≥ y, then Z z z 1{K≥y}[(ye − K)1 K + (K − ye )1{y>K}]n(t, dz, y) {z>ln ( y )} ℝ Z +∞ z ln ( K ) = y(e − e y ) n(t, dz, y). K ln ( y )
5 ∙ If K < y, then Z z z 1{K
Using integration by parts, t,y can be equivalently expressed as ( R z (ez − eu) n(t, du, y) z < 0 (z) = −∞ t,y R ∞ u z z (e − e ) n(t, du, y) z > 0 Hence: Z K [(yez − K)+ − ez(y − K)+ − K(ez − 1)1 ]n(t, dz, y) = y ln . {y>K} t,y y ℝ
1.2 Derivation of the forward equation In this section we present a proof of Theorem 1 using the Tanaka-Meyer formula for semimartingales [24, Theorem 9.43] under assumption (H). Proof. We first note that, by replacing by the conditional measure given ℙ ℙ∣ℱt0
ℱt0 , we may replace the conditional expectation in (2) by an expectation with respect to the marginal distribution pS (dy) of S under . Thus, without T T ℙ∣ℱt0 loss of generality, we set t0 = 0 in the sequel and consider the case where ℱ0 is hal-00445641, version 3 - 21 Apr 2011 the -algebra generated by all ℙ-null sets and we denote C0(T,K) ≡ C(T,K) for simplicity. (2) can be expressed as Z − R T r(t) dt + S C(T,K) = e 0 (y − K) pT (dy). (7) + ℝ By differentiating with respect to K, we obtain: Z ∞ ∂C − R T r(t) dt S − R T r(t) dt (T,K) = −e 0 p (dy) = −e 0 1 , ∂K T E {ST >K} K (8) 2 ∂ C − R T r(t) dt S (T, dy) = e 0 p (dy). ∂K2 T
K K Let Lt = Lt (S) be the semimartingale local time of S at K under ℙ (see [24, Chapter 9] or [33, Ch. IV] for definitions). For ℎ > 0, applying the Tanaka-
6 + Meyer formula to (St − K) between T and T + ℎ, we have
Z T +ℎ 1 (S − K)+ = (S − K)+ + 1 dS + (LK − LK ) T +ℎ T {St−>K} t 2 T +ℎ T T (9) X + + + (St − K) − (St− − K) − 1{St−>K}ΔSt. T As noted in Remark 2, the integrability condition (H) implies that the dis- ˆ − R t r(s) ds counted price St = e 0 St = ℰ(U)t is a martingale under ℙ. So (1) can be expressed as dSt = r(t)St−dt + dSˆt and Z T +ℎ Z T +ℎ Z T +ℎ ˆ 1{St−>K}dSt = 1{St−>K}dSt + r(t)St−1{St−>K}dt T T T where the first term is a martingale. Taking expectations, we get: R T +ℎ r(t) dt R T r(t) dt e 0 C(T + ℎ, K) − e 0 C(T,K) "Z T +ℎ # 1 K K = E r(t)St 1{St−>K}dt + (LT +ℎ − LT ) T 2 ⎡ ⎤ X + + + E ⎣ (St − K) − (St− − K) − 1{St−>K}ΔSt⎦ . T + Noting that St−1{St−>K} = (St− − K) + K1{St−>K}, we obtain " T +ℎ # T +ℎ Z Z R t ∂C 0 r(s) ds E r(t)St−1{St−>K}dt = r(t)e C(t, K) − K (t, K) dt T T ∂K using Fubini’s theorem and (8). As for the jump term, ⎡ ⎤ hal-00445641, version 3 - 21 Apr 2011 X + + E ⎣ (St − K) − (St− − K) − 1{St−>K}ΔSt⎦ T "Z T +ℎ Z # x + + x = E dt m(t, dx)(St−e − K) − (St− − K) − 1{St−>K}St−(e − 1) T Z T +ℎ Z h x + + = E dt m(t, dx) (St−e − K) − (St− − K) T + x x i −(St− − K) (e − 1) − K1{St−>K}(e − 1) Z T +ℎ Z h x + = dtE m(t, dx) (St−e − K) T x + x i −e (St− − K) − K1{St−>K}(e − 1) 7 Applying Lemma 1 to the random measure m we obtain: Z x + x + x K m(t, dx) (St−e −K) −e (St−−K) −K1{St−>K}(e −1) = St− t,St− ln St− holds, leading to: ⎡ ⎤ X + + E ⎣ (St − K) − (St− − K) − 1{St−>K}ΔSt⎦ T Z T +ℎ K = dt E St− t,St− ln T St− Z T +ℎ K = dt E St−E t,St− ln ∣St− T St− Z T +ℎ K = dt E St− t,St− ln (10) T St− ∞ Let ' ∈ C0 (]0, ∞[) be an infinitely differentiable function with compact sup- port. The occupation time formula (see [24, Theorem 9.46]) yields: Z +∞ Z T +ℎ Z T +ℎ K K c 2 2 dK '(K)(LT +ℎ − LT ) = '(St−)d[S]t = dt '(St−)St−t 0 T T ˆ + (ST ) is a martingale, hence : E[ST ] < ∞. Since (ST +ℎ − K) < ST +ℎ, + P + + (ST −K) < ST , ∣ T Z ∞ "Z T +ℎ # hal-00445641, version 3 - 21 Apr 2011 K K 2 2 E dK'(K)(LT +ℎ − LT ) = E '(St−)St−t dt 0 T Z T +ℎ 2 2 = dt E '(St−)St−t T Z T +ℎ 2 2 = dt E E '(St−)St−t ∣St− T "Z T +ℎ # 2 2 = E dt '(St−)St−(t, St−) T Z ∞ Z T +ℎ 2 2 S = '(K)K (t, K) pt (dK) dt 0 T T +ℎ ∞ 2 Z R t Z ∂ C 0 r(s) ds 2 2 = dt e '(K)K (t, K) 2 (t, dK) T 0 ∂K 8 where the last line is obtained by using (8). Gathering together all the terms, we obtain: Z ∞ h R T +ℎ r(t) dt R T r(t) dt i dK '(K) e 0 C(T + ℎ, K) − e 0 C(T,K) 0 Z T +ℎ Z ∞ R t r(s) ds ∂C = dt r(t) e 0 dK '(K)[C(t, K) − K (t, K)] T 0 ∂K T +ℎ ∞ 2 Z R t Z '(K) ∂ C 0 r(s) ds 2 2 + dt e K (t, K) 2 (t, dK) T 0 2 ∂K Z T +ℎ Z ∞ K + dt dK '(K) E St− t,St− ln (11) T 0 St− Dividing by ℎ and taking the limit ℎ → 0 yields: Z ∞ R T r(t) dt ∂C dK'(K)e 0 (T,K) + r(T )C(T,K) 0 ∂T Z ∞ R T r(t) dt ∂C = dK'(K)e 0 r(T )C(T,K) − r(T )K (T,K) 0 ∂K ∞ 2 Z '(K) R T ∂ C 2 2 0 r(t) dt + K (t, K) e 2 (T, dK) 0 2 ∂K Z ∞ K + dK '(K)E St− T,ST − ln 0 ST − Z ∞ R T r(t) dt ∂C = dK'(K)e 0 r(T )C(T,K) − r(T )K (T,K) 0 ∂K ∞ 2 Z '(K) R T ∂ C 2 2 0 r(t) dt + K (t, K) e 2 (T, dK) 0 2 ∂K Z ∞ Z +∞ 2 hal-00445641, version 3 - 21 Apr 2011 R T ∂ C K 0 r(t) dt + dK '(K)e y 2 (T, dy) T,y ln (12) 0 0 ∂K y ∞ Since this equality holds for any ' ∈ C0 (]0, ∞[, ℝ), C(., .) is a solution of ∂C K2(t, K)2 ∂2C ∂C (T,K) = (T,K) − r(T )K (T,K) ∂T 2 ∂K2 ∂K Z +∞ ∂2C K + y 2 (T, dy) T,y ln 0 ∂K y in the sense of distributions on [0,T ]×]0, ∞[. 1.3 Uniqueness of solutions of the forward PIDE Theorem 1 shows that the call price (T,K) 7→ Ct0 (T,K) solves the forward PIDE (5). Uniqueness of the solution of such PIDEs has has been shown using 9 analytical methods [4, 21] under various types of conditions on the coefficients . We give below a direct proof of uniqueness for (5) using a probabilistic method, under explicit conditions which cover most examples of models used in finance. Define, for u ∈ ℝ, t ∈ [0,T [, z > 0 the measure n(t, du, z) by −u ∂ n(t, [u, ∞[, z) = −e [ t,z(u)] u > 0 ∂u (13) ∂ n(t, ] − ∞, u], z) = e−u [ (u)] u < 0 ∂u t,z Throughout this section, we make the following assumption: and Assumption 3. ∀T > 0, ∀B ∈ ℬ(ℝ) − {0}, (t, z) → (t, z), (t, z) → n(t, B, z) are continuous in z ∈ ℝ+, uniformly in t ∈ [0,T ] and Z + 2 ′ ∃KT > 0, ∀(t, z) ∈ [0,T ] × ℝ , ∣(t, z)∣ + (1 ∧ ∣z∣ ) n(t, du, z) ≤ KT (H ) ℝ Note that (H′) implies our previous assumption (H). Theorem 2. Under Assumption 3, if either (i) ∀R > 0 ∀t ∈ [0,T [, inf (t, z) > 0 {0≤z≤R} or (ii) (t, z) ≡ 0 and ∃ ∈]0, 2[, ∃C > 0, ∀R > 0, ∀(t, z) ∈ [0,T [×[0,R], Z C du ∀f ∈ C0( − {0}, ), n(t, du, z) − f(u) ≥ 0 0 ℝ ℝ+ ∣u∣1+ Z ′ C du ′ ∃KT,R > 0, ∣u∣ n(t, du, z) − 1+ dt ≤ KT,R {∣u∣≤1} ∣u∣ Z T and (iii) lim sup n (t, {∣u∣ ≥ R}, z) dt = 0 R→∞ 0 z∈ + hal-00445641, version 3 - 21 Apr 2011 ℝ then the call option price (T,K) 7→ Ct0 (T,K), as a function of maturity and strike, is the unique solution (in the sense of distributions) of the partial integro- differential equation (5) on [t0, ∞[×]0, ∞[ with the initial condition: ∀K > 0 Ct0 (t0,K) = (St0 − K)+. The proof uses the uniqueness of the solution of the forward Kolmogorov equation associated to a certain integro-differential operator. We start with the following result, which has some independent interest: ∞ Proposition 1. Define for t ∈ [0,T ] and f ∈ C0 (ℝ), the integro-differential operator Lt given by x2(t, x)2 L f(x) = r(t)xf ′(x) + f ′′(x) t 2 Z (14) + [f(t, xey) − f(t, x) − x(ey − 1).f ′(x)] n(t, dy, x) ℝ 10 Under Assumption 3, if either conditions (i) or (ii) and (iii) of Theorem 2 + hold, then for each x0 in ℝ , there exists a unique family (pt(x0, dy), t ≥ 0) of bounded measures such that Z dp Z ∀g ∈ C∞(]0, ∞[, ), g(y) (x , dy) = p (x , dy)L g(y) p (x ,.) = 0 ℝ dt 0 t 0 t 0 0 x0 (15) where x0 is the point mass at x0. Furthermore, pt(x0,.) is a probability measure on [0, ∞[. Proof. Denote by (Xt)t∈[0,T ] the canonical process on D([0,T ], ℝ+). Under assumptions (i) (or (ii)) and (iii), Lt verifies Assumptions 1–4 in [31] and by [31, Theorem 1], the martingale problem for (Lt)t∈[0,T ] is well-posed: for any x0 ∈ ℝ, s ∈ [0,T [, there exists a unique probability measure ℚs,x0 on D([0,T ], ℝ+) ∞ + such that ℚs,x0 (Xs = x0) = 1 and for f ∈ C0 (ℝ ): Z t f(Xt) − f(x0) − Luf(Xu) du s is a ℚs,x0 -martingale. Furthermore, (Xt) is a Markov process under ℚx0 , and (Pt)t∈[0,T ] defined by 0 + ℚx0 ∀f ∈ Cb (ℝ ) Ptf(x0) = E [f(Xt)] (16) is a (non-homogeneous) positive strongly continuous contraction semigroup on 0 + Cb (ℝ ) [19, Chapter 1]. If pt(x0, dy) denotes the law of (Xt) starting from x0 under ℚ, the martingale property shows that pt(x0, dy) satisfies the equation (15) that we simply rewrites after integration with respect to time t: Z Z t Z pt(x0, dy)g(y) = g(x0) + ps(x0, dy)Lsg(y) ds (17) 0 hal-00445641, version 3 - 21 Apr 2011 This solution of (15) is in particular positive with mass 1. ∞ + 1 To show uniqueness, let f ∈ C0 (ℝ ) and ∈ C ([0,T ]) and consider the non- time dependent operator A mapping functions of the form (t, x) ∈ [0,T ] × ℝ → ∞ + 1 f(x) (t), which will be denoted C0 (ℝ ) ⊗ C ([0,T ]), into : ′ A(f )(t, x) = (t)Ltf(x) + f(x) (t) (18) Using [19, Theorem 7.1 and Theorem 10.1,Chapter 4]), uniqueness holds for ∞ + the martingale problem associated to the operator L on C0 (ℝ ) if and only if ∞ + uniqueness holds for the martingale problem associated to the A on C0 (ℝ ) ⊗ 1 + C ([0,T ]). For any x0 in ℝ , if (X, ℚx0 ) is a solution of the martingale problem L, then the law of t = (t, Xt) is a solution of the martingale problem for A: for ∞ + any f ∈ C0 (ℝ ) and ∈ C([0,T ]): Z Z t Z pt(x0, dy)f(y) (t) = f(x0) (0) + ps(x0, dy)A(f )(s, y) ds (19) 0 11 Assume there exists a measure qt(dy) such that q0(dy) = x0 (dy) solution of (17), then after integration by parts: Z Z t Z qt(dy)f(y) (t) = f(x0) (0) + qs(dy)A(f )(s, y) ds (20) 0 holds. ∞ + 1 Define, for t in [0,T ], g ∈ C0 (ℝ ) ⊗ C ([0,T ]) Z Ptg(x0) = pt(x0, dy)g(t, y) Z Qtg = qt(dy)g(t, y) Given (19) and (20), for all > 0: Z t Z Z t Pt(f )(x0) − P(f )(x0) = pu(x0, dy)A(f )(u, y) du = Pu(A(f ))(x0) du Z t Z Z t Qt(f ) − Q(f ) = qu(dy)A(f )(u, y) du = Qu(A(f )) du (21) Since the functions t → (t, .), and t → n(t, B, .) for any B ∈ ℬ(ℝ) − {0} are ∞ + bounded in t on [0,T ], it implies that for any fixed f ∈ C0 (ℝ ) and any fixed ∈ 1 C ([0,T ]), t → QtA(f ) and t → PtA(f )(x0) are bounded on [0,T ] and shows ∞ + 1 that Qt. ans Pt.(x0) are weakly right-continuous in t on C0 (ℝ ) ⊗ C ([0,T ]), i.e, for T ≥ t′ ≥ t: lim Pt′ (f )(x0) = Pt(f )(x0) lim Qt′ (f ) = Qt(f ) t′→t t′→t Fix > 0, we have hal-00445641, version 3 - 21 Apr 2011 Z ∞ Z ∞ Z t −t −t e Pt(f )(x0) dt = f(x0) (0) + e Ps(A(f ))(x0) ds dt 0 0 0 Z ∞ Z ∞ −t −t = f(x0) (0) + e e dt Ps(A(f ))(x0) ds 0 s Z ∞ −s = f(x0) (0) + e Ps(A(f ))(x0) ds 0 Consequently, Z ∞ Z ∞ −t −t e Pt( − A)(f )(x0) dt = f(x0) (0) = e Qt( − A)(f ) dt (22) 0 0 0 + Since Pt is a positive strongly continuous contraction semigroup on Cb (ℝ ) for ∞ + the operator L on C0 (ℝ ), one can easily show that it holds for the operator A ∞ + 1 on the domain C0 (ℝ ) ⊗ C ([0,T ]). Applying the Hille-Yosida theory (see [19, 12 0 + Proposition 2.1 and Theorem 2.6]), for all > 0, ℛ( − A) = Cb (ℝ × [0,T ]), ∞ + 1 where ℛ( − A) denotes the image of C0 (ℝ ) ⊗ C ([0,T ]) by the mapping 0 + g → ( − A)g. Hence, since (22) holds then for all ℎ in Cb (ℝ × [0,T ]) : Z ∞ Z ∞ −t −t e Ptℎ (x0) dt = e Qtℎ dt (23) 0 0 0 + Since Cb (ℝ ×[0,T ]) is separating (see [19, Proposition 4.4, Chapter 3]), Pt.(x0) and Qt. are weakly right-continuous and (22) holds for any > 0, the flows 0 + 0 + qt(dy) and pt(x0, dy) are the same on Cb (ℝ × [0,T ]) and obviously on Cb (ℝ ). This ends the proof. We can now study the uniqueness of the forward PIDE (5) and prove The- orem 2 Proof. of Theorem 2. If one decomposes Lt into a differential and an integral component: Lt = At + Bt y2(t, y)2 A f(y) = r(t)yf ′(y) + f ′′(y) t 2 Z z z ′ Btf(y) = [f(ye ) − f(y) − y(e − 1)f (y)]n(t, dz, y) ℝ ∂ + then using the fact that y ∂y (y − x) = x1{y>x} + (y − x)+ = y 1{y>x} and ∂2 + ∂y2 (y − x) = x(y) where x is a unit mass at x, we obtain y2(t, y)2 A (y − x)+ = r(t)y 1 + (y) t {y>x} 2 x and Z hal-00445641, version 3 - 21 Apr 2011 + z + + z + BT (y − x) = [(ye − x) − (y − x) − (e − 1) x1{y>x} + (y − x) ]n(t, dz, y) Zℝ z + z + z = [(ye − x) − e (y − x) − x(e − 1)1{y>x}]n(t, dz, y) ℝ Then, using Lemma 1 for the random measure n(t, dz, y) and t,y its exponential double tail: x B (y − x)+ = y ln t t,y y Hence, the following identity holds: y2(t, y)2 x L (y − x)+ = r(t) x1 + (y − x) + (y) + y ln t {y>x} + 2 x t,y y (24) 13 Let f :[t0, ∞[×]0, ∞[7→ ℝ be a solution in the sense of distributions of (5) with + the initial condition : f(0, x) = (S0 − x) . Integration by parts yields Z ∞ 2 ∂ f + 2 (t, dy)Lt(y − x) 0 ∂x Z ∞ ∂2f y2(t, y)2 x = 2 (t, dy) r(t)(x1{y>x} + (y − x)+) + x(y) + y t,y ln 0 ∂x 2 y Z ∞ 2 Z ∞ 2 ∂ f ∂ f + = −r(t)x 2 (t, dy)1{y>x} + r(t) 2 (t, dy)(y − x) 0 ∂x 0 ∂x x2(t, x)2 ∂2f Z ∞ ∂2f x + 2 + 2 (t, dy)y t,y ln 2 ∂x 0 ∂x y ∂f x2(t, x)2 ∂2f Z ∞ ∂2f x = −r(t)x + r(t)f(t, x) + 2 + 2 (t, dy) y t,y ln ∂x 2 ∂x 0 ∂x y Hence given (5), the following identity holds: Z ∞ 2 ∂f ∂ f + (t, x) = −r(t)f(t, x) + 2 (t, dy)Lt(y − x) (25) ∂t 0 ∂x or equivalently after integration with respect to time t: ∞ 2 R t Z R t ∂ f 0 r(s) ds 0 r(s) ds + e f(t, x) − f(0, x) = e 2 (t, dy)Lt(y − x) (26) 0 ∂x After integration by parts, one shows that: Z ∞ 2 ∂ f + f(t, x) = 2 (t, dy)(y − x) (27) 0 ∂x Hence (25) rewrites: hal-00445641, version 3 - 21 Apr 2011 ∞ 2 t ∞ 2 Z R t ∂ f Z Z R s ∂ f 0 r(s) ds + + 0 r(u) du + e 2 (t, dy)(y−x) −(S0−x) = e 2 (s, dy)Ls(y−x) ds 0 ∂x 0 0 ∂x (28) R t ∂2f 0 r(s) ds Define qt(dy) ≡ e ∂x2 (t, dy), we have q0(dy) = S0 (dy) = p0(S0, dy). ∞ Take g in C0 (]0, ∞[, ℝ), after integration by parts, one shows that: Z ∞ g(y) = g′′(z)(y − z)+ dz (29) 0 14 R Replacing in g(y)qt(dy) and using (28) ℝ ∞ ∞ 2 Z Z R t ∂ f 0 r(s) ds g(y)qt(dy) ≡ g(y) e 2 (t, dy) 0 0 ∂x ∞ ∞ 2 Z Z R t ∂ f ′′ 0 r(s) ds + = g (z) e 2 (t, dy)(y − z) dz 0 0 ∂x ∞ ∞ t ∞ 2 Z Z Z Z R s ∂ f ′′ + ′′ 0 r(u) du + = g (z)(S0 − z) dz + g (z) e 2 (s, dy)Ls(y − z) dz 0 0 0 0 ∂x t ∞ 2 ∞ Z Z R s ∂ f Z 0 r(u) du ′′ + = g(S0) + e 2 (s, dy)Ls[ g (z)(y − z) dz] 0 0 ∂x 0 Z t Z ∞ = g(S0) + qs(dy)Lsg(y) ds 0 0 R t ∂2f 0 r(s) ds This is the equation (17). Proposition 1 yields to e ∂x2 (t, dy) is unique R t ∂2f 0 r(s) ds that is e ∂x2 (t, dy) ≡ pt(S0, dy) (with the notations in Proposition 1)and leads to the uniqueness of the solution of the forward PIDE (5). 2 Examples We now give various examples of pricing models for which Theorem 1 allows to retrieve or generalize previously known forms of forward pricing equations. 2.1 Itˆoprocesses When (St) is an Itˆoprocess i.e. when the jump part is absent, the forward equation (5) reduces to the Dupire equation [16]. In this case our result reduces to the following: Proposition 2 (Dupire PDE). Consider the price process (St) whose dynamics hal-00445641, version 3 - 21 Apr 2011 under the pricing measure ℙ is given by: Z T Z T ST = S0 + r(t)Stdt + SttdWt (30) 0 0 Define q 2 (t, z) = E [t ∣St = z] If " Z T !# 1 2 E exp t dt < ∞ a.s. (A1a) 2 0 the call option price (2) is a solution (in the sense of distributions) of the partial differential equation: ∂C ∂C K2(T,K)2 ∂2C t0 = −r(T )K t0 + t0 (31) ∂T ∂K 2 ∂K2 15 on [t0, ∞[×]0, ∞[ with the initial condition: ∀K > 0 Ct0 (t0,K) = (St0 − K)+. Notice in particular that this result does not require a non-degeneracy con- dition on the diffusion term. Proof. It is sufficient to take ≡ 0 in (1) then equivalently in (5). We leave the end of the proof to the reader. 2.2 Markovian jump-diffusion models Another important particular case in the literature is the case of a Markov jump-diffusion driven by a Poisson random measure. Andersen and Andreasen [2] derived a forward PIDE in the situation where the jumps are driven by a compound Poisson process with time-homogeneous Gaussian jumps. We will now show here that Theorem 1 implies the PIDE derived in [2], given here in a more general context allowing for a time- and state-dependent L´evymeasure, as well as infinite number of jumps per unit time (“infinite jump activity”). Proposition 3 (Forward PIDE for jump diffusion model). Consider the price process S whose dynamics under the pricing measure ℙ is given by: Z T Z T Z T Z +∞ y St = S0 + r(t)St−dt + St−(t, St−)dBt + St−(e − 1)N˜(dtdy) 0 0 0 −∞ (32) where Bt is a Brownian motion and N a Poisson random measure on [0,T ] × ℝ with compensator (dz) dt, N˜ the associated compensated random measure. Assume: ⎧ ′ ⎨ (., .) is bounded (A1a) R 2y ′ ⎩ {∣y∣>1} e (dy) < ∞ (A2a) Then the call option price R T hal-00445641, version 3 - 21 Apr 2011 − r(t) dt t0 ℙ Ct0 (T,K) = e E [max(ST − K, 0)∣ℱt0 ] is a solution (in the sense of distributions) of the partial integro-differential equation: ∂C ∂C K2(T,K)2 ∂2C t0 = −r(T )K t0 + t0 ∂T ∂K 2 ∂K2 (33) Z ∂C + (dz) ez C (T, Ke−z) − C (T,K) − K(e−z − 1) t0 t0 t0 ∂K ℝ on [t0, ∞[×]0, ∞[ with the initial condition: ∀K > 0 Ct0 (t0,K) = (St0 − K)+. Proof. As in the proof of Theorem 1, by replacing ℙ by the conditional measure given ℱ , we may replace the conditional expectation in (2) by an expecta- ℙℱt0 t0 tion with respect to the marginal distribution pS (dy) of S under . Thus, T T ℙ∣ℱt0 without loss of generality, we put t0 = 0 in the sequel, consider the case where ℱ0 16 is the -algebra generated by all ℙ-null sets and we denote C0(T,K) ≡ C(T,K) for simplicity. By differentiating (2) in the sense of distributions with respect to K, we obtain: ∞ 2 ∂C R T Z ∂ C R T − 0 r(t) dt S − 0 r(t) dt S (T,K) = −e pT (dy), 2 (T, dy) = e pT (dy). (34) ∂K K ∂K In this particular case, m(t, dz) dt ≡ (dz) dt and t is simply defined by: ( R z dx ex R x (du) z < 0 (z) ≡ (z) = −∞ −∞ t R +∞ x R ∞ z dx e x (du) z > 0 Then (4) yields t,y(z) ≡ (z) = (z). Let now focus on the term Z +∞ ∂2C K y 2 (T, dy) ln 0 ∂K y in (5). Applying Lemma 1: Z +∞ ∂2C K y 2 (T, dy) ln 0 ∂K y Z ∞ Z − R T r(t) dt S z + z + z = e 0 pT (dy) [(ye − K) − e (y − K) − K(e − 1)1{y>K}](dz) 0 ℝ Z Z ∞ z − R T r(t) dt S −z + + −z = e e 0 pT (dy)[(y − Ke ) − (y − K) − K(1 − e )1{y>K}] (dz) ℝ 0 Z ∂C = ez C(T, Ke−z) − C(T,K) − K(e−z − 1) (dz) (35) ∂K ℝ This ends the proof. 2.3 Pure jump processes hal-00445641, version 3 - 21 Apr 2011 We now consider price processes with no Brownian component. Assumption (H) then reduces to " Z T Z !# y 2 ∀T > 0, E exp dt (e − 1) m(t, dy) < ∞ (A2a) 0 and the forward equation for call option becomes ∂C ∂C Z +∞ ∂2C K + r(T )K = y 2 (T, dy) T,y ln (36) ∂T ∂K 0 ∂K y It is convenient to use the change of variable: v = ln y, k = ln K. Define, c(k, T ) = C(ek,T ). Then one can write this PIDE as: Z +∞ 2 ∂c ∂c 2(v−k) ∂ c ∂c + r(T ) = e 2 − (T, dv) T,v(k − v) (37) ∂T ∂k −∞ ∂k ∂k 17 where T,v is defined by: v T,v(z) = E [ T (z)∣ST − = e ] with: ( R z dx ex R x m(T, du) z < 0 (z) = −∞ −∞ T R +∞ x R ∞ z dx e x m(T, du) z > 0 In the case, considered in [9], where the L´evydensity mY has a deterministic separable form: mY (t, dz, y) dt = (y, t) k(z) dz dt (38) Equation (37) allows us to recover1 equation (14) in [9]: Z +∞ 2 ∂c ∂c 2(v−k) v ∂ c ∂c + r(T ) = (k − v)e (e ,T ) 2 − (T, dv) ∂T ∂k −∞ ∂k ∂k where is defined as the exponential double tail of k(u) du, i.e: ( R z dx ex R x k(u) du z < 0 (z) = −∞ −∞ R +∞ x R ∞ z dx e x k(u) du z > 0 The right hand side can be written as a convolution of distributions: ∂c ∂c ∂2c ∂c + r(T ) = [a (.) − ] ∗ g where (39) ∂T ∂k T ∂k2 ∂k −2u u g(u) = e (u) aT (u) = (e ,T ) (40) Therefore, it implies that from the knowledge of c(., .) and a choice for (.) we can recover aT hence (., .). As noted by Carr et al. [9], this equation is analogous to the Dupire formula for diffusions: it enables to “invert” the struc- ture of the jumps–represented by – from the cross-section of option prices. hal-00445641, version 3 - 21 Apr 2011 Note that, like the Dupire formula, this inversion involves a double deconvolu- tion/differentiation of c which illustrates the ill-posedness of the inverse problem. 2.4 Time changed L´evyprocesses Time changed L´evy processes were proposed in [8] in the context of option pricing. Consider the price process S whose dynamics under the pricing measure ℙ is given by: t R t Z 0 r(u) du St ≡ e Xt Xt = exp (LΘt )Θt = sds (41) 0 2 where Lt is a L´evyprocess with characteristic triplet (b, , ), N its jump measure and (t) is a locally bounded positive semimartingale. We assume L 1Note however that the equation given in [9] does not seem to be correct: it involves the double tail of k(z) dz instead of the exponential double tail. 18 and are ℱt-adapted. − R t r(u) du Xt ≡ e 0 St is a martingale under the pricing measure ℙ if exp (Lt) is a martingale which requires the following condition on the characteristic triplet of (Lt): 1 Z b + 2 + (ez − 1 − z 1 )(dy) = 0 (42) 2 {∣z∣≤1} ℝ Define the value Ct0 (T,K) at time t0 of the call option with expiry T > t0 and strike K > 0 of the stock price (St): R T − 0 r(t) dt ℙ Ct0 (T,K) = e E [max(ST − K, 0)∣ℱt0 ] (43) Proposition 4. Define (t, x) = E[t∣Xt− = x] and the exponential double tail of (du) ( R z dx ex R x (du) z < 0 (z) = −∞ −∞ (44) R +∞ x R ∞ z dx e x (du) z > 0 Assume = 1 2 + R (ey − 1)2(dy) < ∞ holds and 2 ℝ E [exp ( ΘT )] < ∞ (45) Then the call option price Ct0 :(T,K) 7→ Ct0 (T,K) at date t0, as a function of maturity and strike, is a solution (in the sense of distributions) of the partial integro-differential equation: ∂C ∂C K2 (T,K)2 ∂2C = −r (T,K)K + ∂T ∂K 2 ∂K2 (46) Z +∞ ∂2C K hal-00445641, version 3 - 21 Apr 2011 + y 2 (T, dy) (T, y) ln 0 ∂K y on [t, ∞[×]0, ∞[ with the initial condition: ∀K > 0 Ct0 (t0,K) = (St0 − K)+. Proof. Using [5, Theorem 4], (LΘt ) writes t t Z p Z LΘt = L0 + sdBs + bsds 0 0 Z t Z Z t Z + s zN˜(ds dz) + s zN(ds dz) 0 ∣z∣≤1 0 {∣z∣>1} where N is an integer-valued random measure with compensator (dz) dt, N˜ its 19 compensated random measure. Applying the Itˆoformula yields Z t 1 Z t X X = X + X dL + X 2 ds + X − X − X ΔL t 0 s− Ts 2 s− s s s− s− Ts 0 0 s≤t Z t Z t 1 2 p = X0 + Xs− bs + s ds + Xs− s dBs 0 2 0 Z t Z Z t Z + Xs−s zN˜(ds dz) + Xs−s zN(ds dz) 0 {∣z∣≤1} 0 {∣z∣>1} Z t Z z + Xs−s(e − 1 − z)N(ds dz) 0 ℝ R z Under our assumptions, (e − 1 − z 1{∣z∣≤1})(dz) < ∞, hence: Z t 1 Z Z t √ X = X + X b + 2 + (ez − 1 − z 1 ) (dz) ds + X dB t 0 s− s 2 s {∣z∣≤1} s s− s 0 ℝ 0 Z t Z z + Xs−s(e − 1)N˜(ds dz) 0 ℝ Z t Z t Z p z = X0 + Xs− s dBs + Xs−s(e − 1)N˜(ds dz) 0 0 ℝ and (St) may be expressed as: Z t Z t Z t Z p z St = S0 + Ss−r(s) ds + Ss− s dBs + Ss−s(e − 1)N˜(ds dz) 0 0 0 ℝ Assumption (45) implies that assumption (H) of Theorem 1 and (St) is now in the suitable form (1) to apply Theorem 1, which yields the result. 2.5 Index options in a multivariate jump-diffusion model hal-00445641, version 3 - 21 Apr 2011 Consider a multivariate model with d assets: Z T Z T Z T Z i i i i i i yi ˜ ST = S0 + r(t)St− dt + St− tdWt + St− (e − 1)N(dt dy) d 0 0 0 ℝ where i is an adapted process taking values in ℝ representing the volatility of asset i, W is a d-dimensional Wiener process, N is a Poisson random measure on [0,T ] × ℝd with compensator (dy) dt, N˜ denotes its compensated random measure. i i j The Wiener processes W are correlated: for all 1 ≤ (i, j) ≤ d, ⟨W ,W ⟩t = i,jt, with ij > 0 and ii = 1. An index is defined as a weighted sum of the asset prices: d X i It = wiSt i=1 20 The value Ct0 (T,K) at time t0 of an index call option with expiry T > t0 and strike K > 0 is given by − R T r(t) dt t0 ℙ Ct0 (T,K) = e E [max(IT − K, 0)∣ℱt0 ] (47) The following result is a generalization the forward PIDE studied by Avellaneda et al. [3] for the diffusion case: Theorem 3. Forward PIDE for index options. Assume ⎧ h T i ∀T > 0 exp 1 R ∥ ∥2 dt < ∞ (A ) E 2 0 t 1b ⎨ R d (1 ∧ ∥y∥) (dy) < ∞ a.s. (A2b) (48) ℝ ⎩ R 2∥y∥ {∥y∥>1} e (dy) < ∞ a.s. (A3b) Define ⎧ R z x R dx e 1 P i y (dy) z < 0 −∞ d 1≤i≤d−1 wiS e i ℝ ln t− ≤x ⎨ It− t(z) = R ∞ x R (49) dx e 1 P i y (dy) z > 0 z d 1≤i≤d−1 wiS e i ℝ ln t− ≥x ⎩ It− and v u ⎡⎛ ⎞ ⎤ u d 1 X i j i j (t, z) = u w w S S ∣I − = z ; (50) z ⎷E ⎣⎝ i j ij t t t− t−⎠ t ⎦ i,j=1 t,y(z) = E [t (z) ∣It− = y] (51) The index call price (T,K) 7→ Ct0 (T,K), as a function of maturity and strike, is a solution (in the sense of distributions) of the partial integro-differential equation: hal-00445641, version 3 - 21 Apr 2011 ∂C ∂C (T,K)2 ∂2C Z +∞ ∂2C K = −r(T )K + 2 + y 2 (T, dy) T,y ln ∂T ∂K 2 ∂K 0 ∂K y (52) on [t0, ∞[×]0, ∞[ with the initial condition: ∀K > 0 Ct0 (t0,K) = (It0 − K)+. Proof. (Bt)t≥0 defined by Pd w Si idW i dB = i=1 i t− t t t 1/2 Pd i j i j i,j=1 wiwjij tt St−St− is a continuous local martingale with quadratic variation t: by L´evy’stheorem, 21 B is a Brownian motion. Hence I may be decomposed as 1 ⎛ ⎞ 2 d Z T Z T d X i X i j i j IT = wiS0 + r(t)It− dt + ⎝ wiwjij tt St−St−⎠ dBt 0 0 i=1 i,j=1 (53) Z T Z d X i yi ˜ + wiSt−(e − 1)N(dt dy) d 0 ℝ i=1 The essential part of the proof consists in rewriting (It) in the suitable form (1) to apply Theorem 1. Applying the Itˆoformula to ln (IT ) yields: ln (IT ) − ln (I0) d Z T h 1 X = r(t) − w w ij Si Sj 2I2 i j ij t t t− t− 0 t− i,j=1 P i yi P i yi Z wiS e wiS e i − 1≤i≤d t− − 1 − ln 1≤i≤d t− (dy) dt It− It− 1 ⎛ d ⎞ 2 Z T 1 X + w w ij Si Sj dB I ⎝ i j ij t t t− t−⎠ t 0 t− i,j=1 T P i yi ! Z Z wiS e + ln 1≤i≤d t− N˜(dt dy) 0 It− The last equality is obtained since P w Si eyi P w Si eyi R 1≤i≤d i t− − 1 − ln 1≤i≤d i t− (dy) < ∞: using the convexity It− It− w Si property of the logarithm (one recalls that P i t− = 1), and the H¨older 1≤i≤d It− inequality: hal-00445641, version 3 - 21 Apr 2011 P i yi ! i 1≤i≤d wiSt−e X wiSt− X ln ≤ yi ≤ ∣yi∣ ≤ ∥y∥, It− It− 1≤i≤d 1≤i≤d P w Si eyi P w Si eyi hence the functions y → ln 1≤i≤d i t− and y → 1≤i≤d i t− are in- It− It− tegrable with respect to (dy) by assumptions (A2b) and (A3b). We furthermore observe that ! Z P w Si eyi 1≤i≤d i t− 1 ∧ ln (dy) < ∞ a.s. It− (54) P i yi T 1≤i≤d wiSt−e Z Z 2 ln e It− (dy) dt < ∞ a.s. 0 {∥y∥>1} 22 R yi Similarly, since (A2b) and (A3b), for all 1 ≤ i ≤ d, e − 1 − 1{∣yi∣≤1}yi (dy) < i ∞ and ln (ST ) rewrites: Z T 1 Z i i i 2 yi ln (ST ) = ln (S0) + r(t) − (t) − e − 1 − 1{∣yi∣≤1}yi (dy) dt 0 2 Z T Z T Z i i ˜ + t dWt + yi N(dt dy) 0 0 1 d−1 Define the d-dimensional martingale Wt = (Wt , ⋅ ⋅ ⋅ ,Wt ,Bt). For all 1 ≤ i j i, j ≤ d − 1, ⟨W ,W ⟩t = i,jt; Pd j j wjijS ⟨W i,B⟩ = j=1 t− t t t 1/2 Pd i j i j i,j=1 wiwjij tt St−St− Define ⎛ Pd j j ⎞ j=1 wj 1j St−t 1 ⋅ ⋅ ⋅ 1,d−1 1/2 Pd w w ij Si Sj ⎜ ( i,j=1 i j ij t t t− t−) ⎟ ⎜ ⎟ ⎜ . .. . . ⎟ ⎜ . . . . ⎟ Θt = ⎜ Pd j j ⎟ ⎜ j=1 wj d−1,j St−t ⎟ d−1,1 ⋅ ⋅ ⋅ 1 ⎜ Pd i j i j 1/2 ⎟ ⎜ ( wiwj ij S S ) ⎟ ⎜ i,j=1 t t t− t− ⎟ Pd w Sj j Pd w Sj j ⎝ j=1 j 1,j t− t ⋅ ⋅ ⋅ j=1 j d−1,j t− t 1 ⎠ Pd i j i j 1/2 Pd i j i j 1/2 ( i,j=1 wiwj ij tt St−St−) ( i,j=1 wiwj ij tt St−St−) There exists a standard Brownian motion (Zt) such that Wt = AZt where A is a d × d matrix verifying Θ = tA A.