PDF

Small-Time Asymptotics of Option Prices
and First Absolute Moments
βˆ—
Johannes Muhle-Karbe
Marcel Nutz
βˆ—
First version: June 11, 2010. This version: June 11, 2011.
Abstract
We study the leading term in the small-time asymptotics of atthe-money call option prices when the stock price process 𝑆 follows a
general martingale. This is equivalent to studying the rst centered
absolute moment of 𝑆 . We
√ show that if 𝑆 has a continuous part, the
leading term is of order 𝑇 in time 𝑇 and depends only on the initial
value of the volatility. Furthermore, the term is linear in 𝑇 if and only
if 𝑆 is of nite variation. The leading terms for pure-jump processes
with innite variation are between these two cases; we obtain their
exact form for stable-like small jumps. To derive these results, we use
a natural approximation of 𝑆 so that calculations are necessary only
for the class of Lévy processes.
Keywords option price, absolute moment, small-time asymptotics, approximation
by Lévy processes
AMS 2000 Subject Classications Primary 91B25, secondary 60G44.
JEL Classication G13.
Acknowledgements
We thank Martin Keller-Ressel, Sergey Nadtochiy and
Mark Podolskij for discussions, and an anonymous referee for detailed comments.
The rst author was partially supported by the National Centre
of Competence in Research Financial Valuation and Risk Management
(NCCR FINRISK), Project D1 (Mathematical Methods in Financial Risk
Management), of the Swiss National Science Foundation (SNF). The second
author acknowledges nancial support by Swiss National Science Foundation
Grant PDFM2-120424/1.
1
Introduction
We consider the problem of option pricing in mathematical nance where
the price of an option on a stock
βˆ—
𝑆
is calculated as the expectation under a
ETH Zurich, Department of Mathematics, Rämistrasse 101, 8092 Zurich, Switzerland,
email addresses
[email protected], [email protected].
1
risk-neutral measure. As usual we assume that the stock price is modeled
directly under this measure and set the interest rate to zero.
Therefore,
𝑆 on a ltered probability
(Ξ©, β„±, (ℱ𝑑 )𝑑β‰₯0 , 𝑃 ) satisfying the usual hypotheses; we assume for simplicity that β„±0 is trivial 𝑃 -a.s. Our main interest concerns the small-time
our basic model consists of a càdlàg martingale
space
asymptotics of
European call option prices,
𝐸[(𝑆𝑇 βˆ’ 𝐾)+ ]
where
π‘₯+ := max{π‘₯, 0}
and
𝐾 βˆˆβ„
for
at-the-money
case
𝐾 = 𝑆0 .
(1.1)
strike price of the option. More
in the leading order asymptotics
is the
precisely, we shall mostly be interested
for the
𝑇 ↓ 0,
Option price asymptotics are used in
nance to nd initial values for model calibration procedures and also for
model testing as explained below. From a probabilistic point of view one can
note that, up to a factor two, these are also the asymptotics of the centered
absolute rst moment and, at least for continuous
𝑆,
of the expected local
time at the origin.
In applications
𝑆
is often specied via a stochastic dierential equation
driven by a Brownian motion
with compensator
𝐹 (𝑑π‘₯)𝑑𝑑;
π‘Š and a Poisson random
𝑆 is of the form
measure
𝑁 (𝑑𝑑, 𝑑π‘₯)
i.e.,
𝑆 = 𝑆0 + 𝜎 βˆ™ π‘Š + πœ…(π‘₯) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯) 𝑑𝑑),
(1.2)
𝜎 = (πœŽβˆ«π‘‘ ) and πœ… = (πœ…π‘‘ (π‘₯)) are predictable integrands and we denote
𝑒
by 𝜎
π‘Šβˆ«π‘’ =
0 πœŽπ‘‘ π‘‘π‘Šπ‘‘ the Itô integral and by πœ…(π‘₯) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯) 𝑑𝑑)𝑒 the
∫
𝑒
integral
0 ℝ πœ…π‘‘ (π‘₯) (𝑁 (𝑑𝑑, 𝑑π‘₯) βˆ’ 𝐹 (𝑑π‘₯) 𝑑𝑑) of random measures. In view of
the applications, we shall adopt this representation for 𝑆 , but we recall in
where
βˆ™
Remark A.2 that this entails no essential loss of generality:
every càdlàg
martingale with absolutely continuous characteristics can be represented in
the form (1.2), and this includes all models of interest.
The
main results
for the at-the-money asymptotics (1.1) are obtained
under certain right-continuity conditions on the mapping
exclude the trivial case
√
𝑇
to
𝑇.
If
𝜎0 βˆ•= 0,
𝑆 ≑ 𝑆0 ,
𝑑 7β†’ (πœŽπ‘‘ , πœ…π‘‘ ).
If we
the possible convergence rates range from
i.e., in the presence of a Brownian component, the
√ √
(∣𝜎0 ∣/ 2πœ‹) 𝑇 irrespective of the jumps. On the
other hand, the leading term is 𝐢𝑇 if and only if 𝑆 is of nite variation, and
then 𝐢 is given explicitly in terms of πœ…0 and 𝐹 . For pure-jump processes
√
with innite variation the rate may be anywhere between
𝑇 and 𝑇 and it
need not be a power of 𝑇 . We consider a class of processes, containing most
relevant examples, whose small jumps resemble the ones of an 𝛼-stable Lévy
1/𝛼 while for
process, where 𝛼 ∈ [1, 2). For 𝛼 > 1 the leading term is 𝐢𝑇
𝛼 = 1 we nd 𝐢𝑇 ∣ log 𝑇 ∣; the constants are given explicitly.
leading term is given by
The
basic idea
to obtain these results is to calculate the option price
asymptotics for a simple model
𝑍
which approximates
2
𝑆
in a suitable sense.
More precisely, we obtain a natural approximation by freezing the coecients
𝜎
and
πœ…
in (1.2) at time
𝑑 = 0,
namely
𝑍 = 𝑆0 + 𝜎0 π‘Š + πœ…0 (π‘₯) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯) 𝑑𝑑).
Note that
the
𝜎0
and
πœ…0 (π‘₯)
Lévy process 𝑍
are deterministic since
β„±0
(1.3)
is trivial. We show that
has the same leading order asymptotics as
𝑆
under mild
regularity conditions. Therefore, an explicit treatment is necessary only in
the Lévy case, for which much ner arguments are possible. We prove that
we can pass with an error of order
𝑂(𝑇 )
from one pure-jump Lévy process
to another one when the small jumps have a similar behavior and use this
to reduce even further to very particular Lévy processes.
Literature.
Due to their importance for model calibration and testing,
small-time asymptotics of option prices have received considerable attention
in recent years; see [2, 5, 6, 7, 10, 11, 12, 13, 14, 15, 16, 17, 22, 24].
A
survey of recent literature is given in the introduction of Forde et al. [16].
We shall review only the works most closely related to our study; in particular, we focus on the at-the-money case. Here Carr and Wu [7] is an early
reference with results in the spirit of ours. The authors obtain by partially
heuristic arguments that the order of convergence for nite variation jumps is
√
𝑂( 𝑇 )
in the presence of a Brownian component and
𝑂(𝑇 )
otherwise. This
is in a general model including, e.g., exponential Lévy processes; however,
a boundedness assumption on the coecients of the log price excludes the
application to the Heston model, for example. For the pure-jump case with
innite variation the authors mention that there is a range of possibilities for
the order of convergence and they illustrate this by the so-called log stable
model. Given option price data, the results are used to study whether the
underlying process has jumps.
Durrleman [10] determines the rate
√
𝑂( 𝑇 )
and the corresponding coef-
cient in a similar model, again with bounded coecients and nite variation jumps. The result is stated in terms of the implied volatility, which
is an alternative parametrization for option prices (see also Corollary 3.3).
Forde [14] studies a class of continuous uncorrelated stochastic volatility
models and computes explicitly the rst two leading terms corresponding
to the orders
𝑇 1/2
and
𝑇 3/2
while Forde et al. [16] obtain the same two
coecients for the Heston model with correlation.
Forde and Figueroa-
López [13] determine the leading order term for the CGMY model.
More
generally, Tankov [27] obtains several results similar to ours in the setting of exponential Lévy models; detailed references are given below.
nally, for (arithmetic) Lévy processes
𝑆,
related to the power variation of order one:
𝐸[βˆ£π‘†π‘‡ ∣] =
𝐸[π‘›βˆ’1
βˆ‘
1β‰€π‘˜β‰€π‘› βˆ£π‘†π‘˜π›Ώ
βˆ’ 𝑆(π‘˜βˆ’1)𝛿 ∣]
Fi-
the at-the-money option price is
for
𝑇 = 𝛿 is
any 𝑛 ∈ β„• by
if
the mesh size,
the i.i.d. prop-
erty of the increments. Hence it is not surprising that we shall benet from
results of Jacod [19] on asymptotic properties of power variations.
3
The present paper is organized as follows. Section 2 contains the approximation result, which is stated for general martingales. Section 3 contains
√
the analysis for the order
𝑇
and Section 4 the higher leading orders corre-
sponding to the pure-jump case. An appendix contains some standard results
about Lévy processes that are used in the body of the text, often without
further mention. We refer to the monograph of Jacod and Shiryaev [20] for
any unexplained notion or notation from stochastic calculus.
2
Approximation of the Process 𝑆
In this section we compare two martingales
cients and study the distance
βˆ£π‘†π‘‡ βˆ’
𝑆
and
𝑆′
with dierent coef-
𝑆𝑇′ ∣ in mean as
application to be used in the sequel is the case where
proximation (1.3) of
𝑆.
𝑇 ↓ 0. The main
𝑆 β€² is the Lévy ap-
In that case, the assumption in the following result
becomes a Hölder-type condition in mean for the coecients
πœŽπ‘‘
and
πœ…π‘‘
of
𝑆.
Proposition 2.1. Let 𝑆 be a martingale of the form (1.2) and 𝑆 β€² a martingale of the analogous form 𝑆 β€² = 𝑆0β€² + πœŽβ€² βˆ™ π‘Š + πœ…β€² (π‘₯) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯)𝑑𝑑) with
𝑆0β€² = 𝑆0 . Let 𝛾 β‰₯ 0.
[∫
]
(i) If 𝐸[(πœŽπ‘‘ βˆ’ πœŽπ‘‘β€² )2 ] = 𝑂(𝑑𝛾 ) and 𝐸 ℝ βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…β€²π‘‘ (π‘₯)∣2 𝐹 (𝑑π‘₯) = 𝑂(𝑑𝛾 ),
then 𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑆𝑇′ ∣2 ] = 𝑂(𝑇 (1+𝛾) ) and 𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑆𝑇′ ∣] = 𝑂(𝑇 (1+𝛾)/2 ).
[∫
]
(ii) Let 𝛽 ∈ [1, 2]. If 𝐸 ℝ βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…β€²π‘‘ (π‘₯)βˆ£π›½ 𝐹 (𝑑π‘₯) = 𝑂(𝑑𝛾 ) and 𝜎 ≑ 0, then
𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑆𝑇′ βˆ£π›½ ] = 𝑂(𝑇 (1+𝛾) ) and 𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑆𝑇′ ∣] = 𝑂(𝑇 (1+𝛾)/𝛽 ).
The assertions remain valid if 𝑂(β‹…) is replaced by π‘œ(β‹…) throughout.
Proof.
We set
𝑋 := 𝑆 βˆ’ 𝑆 β€²
and denote by
𝑋 = 𝑋𝑐 + 𝑋𝑑
the decomposition
into the continuous and purely discontinuous martingale parts; note that
𝑋 𝑐 = (𝜎 βˆ’ 𝜎 β€² ) βˆ™ π‘Š
Fix
𝛽 ∈ [1, 2]
and
𝛾 β‰₯ 0.
and
𝑋 𝑑 = (πœ…(π‘₯) βˆ’ πœ…β€² (π‘₯)) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯)𝑑𝑑).
We use the Burkholder-Davis-Gundy inequality
(e.g., Protter [23, Theorem IV.48]) to obtain
( [
]
[
]
[
])
𝐸 βˆ£π‘‹π‘‡ βˆ£π›½ ≀ 2π›½βˆ’1 𝐸 βˆ£π‘‹π‘‡π‘ βˆ£π›½ + 𝐸 βˆ£π‘‹π‘‡π‘‘ βˆ£π›½
( [
)
[
𝛽/2 ]
𝛽/2 ]
≀ 𝐢𝛽 𝐸 βŸ¨π‘‹ 𝑐 , 𝑋 𝑐 βŸ©π‘‡ + 𝐸 [𝑋 𝑑 , 𝑋 𝑑 ]𝑇
𝐢𝛽 depending only on 𝛽 .
𝑋 𝑐 ≑ 0 and we only need to
(2.1)
for a universal constant
We rst treat the pure-
jump case (ii), then
estimate the second ex-
Recall that for a real sequence 𝑦 = (𝑦𝑛 ) the norms
βˆ‘
βˆ₯𝑦βˆ₯ℓ𝑝 = ( 𝑛 βˆ£π‘¦π‘› βˆ£π‘ )1/𝑝 satisfy βˆ₯𝑦βˆ₯ℓ𝑝 β‰₯ βˆ₯𝑦βˆ₯β„“π‘ž for 1 ≀ 𝑝 ≀ π‘ž < ∞. We apply
this for 𝑝 = 𝛽 and π‘ž = 2 to obtain that
(βˆ‘
)𝛽/2 βˆ‘
𝑑
𝑑 𝛽/2
2
[𝑋 , 𝑋 ]𝑇 =
βˆ£Ξ”π‘‹π‘‘ ∣
≀
βˆ£Ξ”π‘‹βˆ£π›½ = βˆ£πœ…(π‘₯) βˆ’ πœ…β€² (π‘₯)βˆ£π›½ βˆ— 𝑁𝑇 ,
pectation in (2.1).
𝑑≀𝑇
𝑑≀𝑇
4
hence using the denition of the compensator and Fubini's theorem we have
[
[
]
𝛽/2 ]
𝐸 [𝑋 𝑑 , 𝑋 𝑑 ]𝑇
≀ 𝐸 βˆ£πœ…(π‘₯) βˆ’ πœ…β€² (π‘₯)βˆ£π›½ βˆ— 𝑁𝑇
[
]
= 𝐸 βˆ£πœ…(π‘₯) βˆ’ πœ…β€² (π‘₯)βˆ£π›½ βˆ— (𝐹 (𝑑π‘₯)𝑑𝑑)𝑇
]
∫ 𝑇 [∫
β€²
𝛽
=
𝐸
βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…π‘‘ (π‘₯)∣ 𝐹 (𝑑π‘₯) 𝑑𝑑.
0
ℝ
𝑂(𝑑𝛾 ), hence the integral is of order
By assumption, the integrand is of order
𝑂(𝑇 (1+𝛾) ) and the rst assertion of (ii) follows by (2.1).
The second assertion
then follows by Jensen's inequality.
We now turn to the case (i). Of course, the previous estimates hold in
particular for
For
𝛽=2
𝛽 = 2,
so it remains to consider the continuous part in (2.1).
this is
[
]
𝐸 βŸ¨π‘‹ 𝑐 , 𝑋 𝑐 βŸ©π‘‡ = 𝐸
[∫
𝑇
]
(πœŽπ‘‘ βˆ’
πœŽπ‘‘β€² )2 𝑑𝑑
∫
𝑇
=
𝐸[(πœŽπ‘‘ βˆ’ πœŽπ‘‘β€² )2 ] 𝑑𝑑
0
0
and so the conclusion is obtained as before. Finally, we note that the proof
remains valid if
𝑂(β‹…)
is replaced by
π‘œ(β‹…)
throughout.
We illustrate the use of Proposition 2.1 by two applications to the approximation of stochastic dierential equations (SDEs). For the sake of clarity,
we do not strive for minimal conditions.
Corollary 2.2. Let 𝑓 : ℝ β†’ ℝ be continuously dierentiable with bounded
derivative and let 𝐿 be a square-integrable Lévy martingale. Then the SDE
𝑑𝑆𝑑 = 𝑓 (π‘†π‘‘βˆ’ ) 𝑑𝐿𝑑 ,
𝑆0 ∈ ℝ
has a unique solution 𝑆 and the Lévy process 𝑍𝑑 = 𝑆0 + 𝑓 (𝑆0 )𝐿𝑑 satises
𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑍𝑇 ∣] = 𝑂(𝑇 ) as 𝑇 ↓ 0.
Proof.
We recall from [23, Theorem V.67] that the SDE has a unique strong
solution
𝑆
𝑑 7β†’ 𝐸[𝑆𝑑2 ] is locally bounded. The Lévy process 𝐿
the form 𝐿 = π‘π‘Š + π‘₯ βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯)𝑑𝑑) and then
and that
a representation of
has
𝑆 = 𝑆0 + 𝑐𝑓 (π‘†βˆ’ ) βˆ™ π‘Š + (𝑓 (π‘†βˆ’ )π‘₯) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯)𝑑𝑑),
𝑍 = 𝑆0 + 𝑐𝑓 (𝑆0 ) βˆ™ π‘Š + (𝑓 (𝑆0 )π‘₯) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯)𝑑𝑑).
It suces to verify the conditions of Proposition 2.1(i) for
𝛾 = 1 and 𝑆 β€² = 𝑍 .
β€² 2
2
2
Since 𝑆𝑑 = π‘†π‘‘βˆ’ 𝑃 -a.s. for each 𝑑, we have 𝐸[(πœŽπ‘‘ βˆ’πœŽπ‘‘ ) ] = 𝑐 𝐸[βˆ£π‘“ (𝑆𝑑 )βˆ’π‘“ (𝑆0 )∣ ]
[∫
]
∫
β€²
2
2
2
and 𝐸
ℝ βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…π‘‘ (π‘₯)∣ 𝐹 (𝑑π‘₯) = 𝐸[βˆ£π‘“ (𝑆𝑑 ) βˆ’ 𝑓 (𝑆0 )∣ ] ∣π‘₯∣ 𝐹 (𝑑π‘₯). The
𝐿 is square-integrable (Lemma A.1(vi)), and so
𝐸[βˆ£π‘“ (𝑆𝑑 ) βˆ’ 𝑓 (𝑆0 )∣2 ] = 𝑂(𝑑). Now 𝑓 is Lipschitz-
last integral is nite since
it suces to show that
continuous by assumption, so it remains to prove that
𝐸[βˆ£π‘†π‘‘ βˆ’ 𝑆0 ∣2 ] = 𝑂(𝑑).
5
For this, in turn, it suces to verify the conditions of Proposition 2.1(i)
𝛾 = 0 and 𝑆 β€² ≑ 𝑆0 (and hence 𝜎 β€² ≑ 0 and πœ…β€² ≑ 0). Indeed, we have
2
2
2
that 𝐸[βˆ£πœŽπ‘‘ ∣ ] = 𝑐 𝐸[βˆ£π‘“ (𝑆𝑑 )∣ ] = 𝑂(1) as the linear growth of 𝑓 ensures that
2
2
𝐸[𝑓
bounded like 𝐸[𝑆𝑑 ], and similarly we have that
[ ∫(𝑆𝑑 ) ] is 2again locally
]
∫ 2
2
𝐸 ℝ βˆ£πœ…π‘‘ (π‘₯)∣ 𝐹 (𝑑π‘₯) = 𝐸[𝑓 (𝑆𝑑 ) ] ∣π‘₯∣ 𝐹 (𝑑π‘₯) = 𝑂(1) as 𝑑 ↓ 0.
for
Remark 2.3.
If the square-integrable Lévy martingale
Brownian component and if its Lévy measure
∫
𝐹
𝐿
does not have a
satises
∣π‘₯βˆ£π›½ 𝐹 (𝑑π‘₯) < ∞
ℝ
𝛽 ∈ [1, 2], then the
𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑍𝑇 ∣] = 𝑂(𝑇 2/𝛽 ).
for some
assertion in Corollary 2.2 can be strengthened
to
In particular, this applies for
𝛽=1
when
𝐿
is
of nite variation. The proof is as above, using part (ii) of Proposition 2.1
instead of part (i).
We give a second example where the coecient of the SDE is not Lipschitzcontinuous, as this sometimes occurs in stochastic volatility models.
Corollary 2.4. Assume that 𝑆 solves the SDE
√
𝑑𝑆𝑑 = 𝑆𝑑 π‘£π‘‘βˆ’ π‘‘π‘Šπ‘‘ ,
𝑆0 ∈ ℝ,
where 𝑣 β‰₯ 0 is a càdlàg adapted process. If 𝑑 7β†’ 𝐸[𝑣𝑑2+πœ€ ] and 𝑑 7β†’ 𝐸[𝑆𝑑4+πœ€ ]
are bounded in
a neighborhood of zero for some πœ€βˆš> 0, then the Lévy process
√
𝑍𝑑 = 𝑆0 + 𝑆0 𝑣0 π‘Šπ‘‘ satises 𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑍𝑇 ∣] = π‘œ( 𝑇 ).
In particular, this applies when 𝑣 is a square-root process, i.e., when 𝑆 is
the Heston model.
Proof.√By
Proposition 2.1(i) applied with
√
𝐸[(𝑆𝑑 𝑣𝑑 βˆ’ 𝑆0 𝑣0
)2 ]
= π‘œ(1);
π‘£βˆ’
not change if one replaces
𝛾 = 0,
it suces to verify that
notice that by continuity of
by
𝑣.
π‘Š
the SDE does
In view of
√
√
√
√
√
𝐸[(𝑆𝑑 𝑣𝑑 βˆ’ 𝑆0 𝑣0 )2 ] = 𝐸[𝑆𝑑2 𝑣𝑑 βˆ’ 𝑆02 𝑣0 ] + 2𝑆0 𝑣0 𝐸[𝑆0 𝑣0 βˆ’ 𝑆𝑑 𝑣𝑑 ]
√
√
2
2
it suces to check that 𝐸[𝑆𝑑 𝑣𝑑 ] β†’ 𝑆0 𝑣0 and 𝐸[𝑆𝑑 𝑣𝑑 ] β†’ 𝑆0 𝑣0 as 𝑑 ↓ 0.
√
Since 𝑆 𝑣 is right-continuous, this readily follows by the Cauchy-Schwarz
inequality and uniform integrability. That the assumptions are satised for
the Heston model follows from, e.g., Cox et al. [9, Section 3] and the proof
of Andersen and Piterbarg [3, Proposition 3.1].
3
Option Price of Order
√
𝑇
The main idea in this section is to calculate the option price for
𝑆
from (1.2)
via the approximation
𝑍 := 𝑆0 + 𝜎0 π‘Š + πœ…0 (π‘₯) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯)𝑑𝑑).
6
(3.1)
We rst have to ensure that this expression makes sense. Indeed, if
𝑆
is a
martingale, it follows that
∫
βˆ£πœ…π‘‘ (π‘₯)∣ ∧ βˆ£πœ…π‘‘ (π‘₯)∣2 𝐹 (𝑑π‘₯) < ∞ 𝑃 βŠ— 𝑑𝑑-a.e.,
(3.2)
ℝ
but of course this may fail on the nullset
standing assumption
that
𝑍
{𝑑 = 0}.
Hence we make the
is well dened and integrable; i.e., that
πœ…0 (π‘₯)
is
Borel-measurable and satises
∫
βˆ£πœ…0 (π‘₯)∣ ∧ βˆ£πœ…0 (π‘₯)∣2 𝐹 (𝑑π‘₯) < ∞.
(3.3)
ℝ
In any reasonable situation, one will be able to infer this condition from (3.2).
√
We can now prove our result for the at-the-money option price of order
𝑇.
In fact, we describe the slightly more general situation of
almost
the-money strikes by considering a deterministic strike function
at-
𝑇 7β†’ 𝐾𝑇
βˆšπΎπ‘‡ β†’ 𝑆0 as 𝑇 ↓ 0. The main observation is that the coecient
of order
𝑇 depends only on the initial value of 𝜎 and that the jumps are
irrelevant at this order. We denote by 𝒩 the Gaussian distribution.
such that
Theorem 3.1. Let 𝑆 be a martingale of the form
and lim 𝐸
2
lim 𝐸[(πœŽπ‘‘ βˆ’ 𝜎0 ) ] = 0
𝑑↓0
√
[∫
𝑑↓0
(1.2)
and assume that
]
βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…0 (π‘₯)∣ 𝐹 (𝑑π‘₯) = 0.
2
ℝ
√
If 𝐾𝑇 = 𝑆0 + πœƒ 𝑇 + π‘œ( 𝑇 ) for some πœƒ ∈ ℝ, then
𝐸[(𝑆𝑇 βˆ’ 𝐾𝑇 )+ ] = 𝐸[𝒩 (βˆ’πœƒ, 𝜎02 )+ ]
√
√
𝑇 + π‘œ( 𝑇 )
as 𝑇 ↓ 0.
(3.4)
In particular, for the at-time-money case 𝐾 ≑ 𝑆0 we have that
Remark
√
∣𝜎0 ∣ √
𝑇 + π‘œ( 𝑇 ) as 𝑇 ↓ 0.
𝐸[(𝑆𝑇 βˆ’ 𝑆0 )+ ] = √
2πœ‹
√
√
3.2. (a) The form 𝐾𝑇 = 𝑆0 + πœƒ 𝑇 + π‘œ( 𝑇 ) chosen
theorem is in fact the only relevant one.
√
𝐾𝑇 β†’ 𝑆0 is slower than ∼ 𝐢 𝑇 ,
+
of 𝐸[(𝑆𝑇 βˆ’ 𝐾𝑇 ) ] will simply be
in the
Indeed, if the convergence
then the leading order asymptotics
determined by
(𝑆0 βˆ’ 𝐾𝑇 )+ ,
and if
it is faster, we nd the same asymptotics as in the at-the-money case.
As usual, we write
𝑓 (𝑇 ) ∼ 𝑔(𝑇 ) if 𝑓 (𝑇 )/𝑔(𝑇 ) β†’ 1 as π‘‡βˆšβ†“ 0. One can
πœƒ satises πœƒ = (𝑆0 βˆ’ 𝐾𝑇 )/ 𝑇 + π‘œ(1) and
also note that the constant
can therefore be interpreted as a
degree of moneyness
for the option;
a similar notion was previously used by Medvedev and Scaillet [22].
√
(b) Even for the class of continuous martingales,
𝑇
is the highest order in
𝜎0 . Indeed,
1/2+πœ€
assume that this were also the case for the order 𝑇
and some πœ€ > 0.
which the option price depends only on the initial volatility
7
Set
𝑆
βˆ«π‘‘
(1 + 2πœ€) 𝑑2πœ€ and dene the martingale 𝑆𝑑 = 0 πœŽπ‘  π‘‘π‘Šπ‘  . Then
βˆ«π‘‡ 2
1+2πœ€ and hence
Gaussian with variance Var(𝑆𝑇 ) =
0 πœŽπ‘‘ 𝑑𝑑 = 𝑇
πœŽπ‘‘ =
is
√
1
𝐸[𝑆𝑇+ ] = √ 𝑇 1/2+πœ€ .
2πœ‹
This coecient is of course dierent from the one for
another martingale with
𝑆 β€² ≑ 0,
which is
𝜎0 = 0.
Let us formulate the theorem once more, in the language preferred by
practitioners. If
𝑆 > 0,
the
implied volatility πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) ∈ [0, ∞] of an at-the-
money call option with maturity
+
(
𝐸[(𝑆𝑇 βˆ’ 𝑆0 ) ] = 𝑆0 Ξ¦
where
Ξ¦
𝑇
is dened as the solution of
√
1
πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) 𝑇
2
)
(
)
√
1
βˆ’ 𝑆0 Ξ¦ βˆ’ πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) 𝑇 ,
2
denotes the standard normal distribution function.
(3.5)
This means
that the option price coincides with its counterpart in a Black-Scholes model
πœŽπ‘–π‘šπ‘π‘™ (𝑇 ). Note that πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) exits and is unique
π‘₯ 7β†’ Ξ¦(π‘₯) βˆ’ Ξ¦(βˆ’π‘₯) is strictly increasing and maps [0, ∞] to [0, 1] and
+
moreover 𝐸[(𝑆𝑇 βˆ’ 𝑆0 ) ] ∈ [0, 𝑆0 ]. The following generalizes the result of [10]
with volatility parameter
since
to unbounded coecients and innite variation jumps.
Corollary 3.3. Let 𝑆 > 0. Under the conditions of Theorem 3.1 we have
πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) β†’ ∣𝜎0 ∣/𝑆0 ; i.e., the implied volatility converges to the spot volatility
of the continuous martingale part of the log price.
Proof.
Ξ¦ (see, e.g.,
√ Abramowitz and SteΞ¦(π‘₯) βˆ’ Ξ¦(βˆ’π‘₯) ∼ 2/πœ‹ π‘₯ for small π‘₯. Since
By the asymptotic properties of
gun [1, Chapter 7]), we have that
√
πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) 𝑇 ↓ 0 by Theorem 3.1, we obtain that
)
(
)
(
√
√
πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) √
1
1
Ξ¦
πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) 𝑇 βˆ’ Ξ¦ βˆ’ πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) 𝑇 ∼ √
𝑇.
2
2
2πœ‹
∣𝜎0 ∣/𝑆0 √
+
√
As Theorem 3.1 yields 𝐸[(𝑆𝑇 βˆ’ 𝑆0 ) ]/𝑆0 ∼
𝑇 , the denition
2πœ‹
πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) shows that πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) β†’ ∣𝜎0 ∣/𝑆0 .
Proof of Theorem 3.1.
of
We may use the put-call parity to rewrite the call
price in terms of an absolute moment: since
𝑆
is a martingale,
𝐸[βˆ£π‘†π‘‡ βˆ’ 𝐾𝑇 ∣] = 𝐸[2(𝑆𝑇 βˆ’ 𝐾𝑇 )+ βˆ’ 𝑆𝑇 + 𝐾𝑇 ]
From this we also see that
= 2𝐸[(𝑆𝑇 βˆ’ 𝐾𝑇 )+ ] βˆ’ 𝑆0 + 𝐾𝑇 .
√
we may assume 𝐾𝑇 = 𝑆0 + πœƒ 𝑇
(√ )
π‘œ 𝑇 part does not matter). By a translation we
𝑆0 = 𝐾0 = 0. (We exploit here that we are working
martingales 𝑆 , not necessarily positive.)
8
(i.e., that the
may also assume that
with a general class of
Step 1: Continuous Lévy case.
Assume rst that 𝑆 is a continuous
𝑆𝑑 = 𝜎0 π‘Šπ‘‘ . Then
√
√
√
𝐸[(𝑆𝑇 βˆ’ πœƒ 𝑇 )+ ] = 𝑇 𝐸[(𝜎0 π‘Š1 βˆ’ πœƒ)+ ] = 𝑇 𝐸[𝒩 (βˆ’πœƒ, 𝜎02 )+ ].
Lévy
process, i.e., that
Step 2: General Lévy case.
denote by
𝑆 = 𝑆𝑐 + 𝑆𝑑
𝑆
Let
be a Lévy process (i.e.,
𝑆 = 𝑍)
and
its decomposition into the continuous and purely
discontinuous martingale parts. In view of
βˆ£π‘†π‘‡π‘ βˆ’ 𝐾𝑇 ∣ βˆ’ βˆ£π‘†π‘‡π‘‘ ∣ ≀ βˆ£π‘†π‘‡ βˆ’ 𝐾𝑇 ∣ ≀ βˆ£π‘†π‘‡π‘ βˆ’ 𝐾𝑇 ∣ + βˆ£π‘†π‘‡π‘‘ ∣
(√ )
𝑑
and Step 1, it suces to show that 𝐸[βˆ£π‘†π‘‡ ∣] is of order π‘œ
𝑇 . To relax the
𝑑
notation, let us assume that 𝑆 = 𝑆 . We can further decompose 𝑆 into
a martingale with bounded jumps, which is in particular square-integrable,
and a compound Poisson process
𝑋
which is integrable due to (3.3). One can
check by direct calculation or by an application of Theorem 4.1 below that
𝐸[βˆ£π‘‹π‘‡ ∣] = 𝑂(𝑇 ).
That is, we may even assume that
𝑆
is a square-integrable
pure-jump Lévy martingale. Then
√
βˆ£π‘†π‘‡ ∣/ 𝑇 β†’ 0
𝑇 ↓ 0;
√
𝐸[βˆ£π‘†π‘‡ ∣]/ 𝑇 β†’ 0,
in probability as
it suces to
√
this set is even bounded
{𝑆𝑇 / 𝑇 }𝑇 >0 . But
√ 2
[
]
2
in 𝐿 (𝑃 ) as 𝑆 is square-integrable; indeed, βˆ₯𝑆𝑇 / 𝑇 βˆ₯ 2
=
𝐸
[𝑆,
𝑆]
due
1
[
[ 2]
]
[
] 𝐿 (𝑃 )
to the relation 𝐸 𝑆𝑇 = 𝐸 [𝑆, 𝑆]𝑇 = 𝑇 𝐸 [𝑆, 𝑆]1 .
see, e.g., [19, Lemma 4.1]. To conclude that
show the uniform integrability of
Step 3: General case.
When
𝑆
is as in the theorem, we approximate
𝑆
by
𝑍 from (3.1). The assumptions of Proposition 2.1(i) are
β€² = 𝑍 and the order π‘œ(1), hence we obtain that 𝐸[βˆ£π‘† βˆ’ 𝑍 ∣]
𝑆
𝑇
𝑇
(√ )
π‘œ 𝑇 . In view of Step 2 applied to 𝑍 and
the Lévy process
satised for
is of order
𝐸[βˆ£π‘π‘‡ βˆ’ 𝐾𝑇 ∣] βˆ’ 𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑍𝑇 ∣] ≀ 𝐸[βˆ£π‘†π‘‡ βˆ’ 𝐾𝑇 ∣] ≀ 𝐸[βˆ£π‘π‘‡ βˆ’ 𝐾𝑇 ∣] + 𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑍𝑇 ∣],
this completes the proof.
4
Option Prices with Higher Leading Orders
In this section, we consider
pure-jump
martingales
𝑆 = 𝑆0 + πœ…(π‘₯) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯)𝑑𝑑),
i.e., we set
𝜎 ≑ 0
in (1.2).
Then the term of order
vanishes and the leading order is higher than
1/2.
(4.1)
√
𝑇
in Theorem 3.1
We shall again use the
approximation result from Section 2 to reduce to the Lévy case and consider
𝑍 = 𝑆0 + πœ…0 (π‘₯) βˆ— (𝑁 βˆ’ 𝐹 (𝑑π‘₯)𝑑𝑑).
(4.2)
However, this case is now more involved since the results depend on the
properties of the Lévy measure.
which ensures that
𝑍
We recall the standing assumption (3.3)
is well dened and integrable.
9
4.1 Finite Variation
We rst treat the case when
𝑆
is of nite variation, which leads to the
highest possible (nontrivial) convergence rate for the at-the-money option
price. Indeed, the following result shows that this class of price processes is
characterized by the rate
𝑂(𝑇 ).
Theorem 4.1. Let[ βˆ«π‘† be a pure-jump martingale
of the form (4.1) and as]
sume that lim𝑑↓0 𝐸 ℝ βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…0 (π‘₯)∣𝐹 (𝑑π‘₯) = 0. Then the following are
equivalent:
(i) 𝑆 is of nite variation on [0, 𝑇 ] for some 𝑇 > 0,
(ii) 𝐸[(𝑆𝑇 βˆ’ 𝑆0 )+ ] = 𝑂(𝑇 ) as 𝑇 ↓ 0.
In that case,
1
𝐸[(𝑆𝑇 βˆ’ 𝑆0 )+ ] = 𝐢 𝑇 + π‘œ(𝑇 )
2
∫
∫
𝐢 := ℝ βˆ£πœ…0 (π‘₯)∣ 𝐹 (𝑑π‘₯) + ℝ πœ…0 (π‘₯) 𝐹 (𝑑π‘₯).
where
Proof.
as 𝑇 ↓ 0,
As in the proof of Theorem 3.1 we may assume that
we have
𝐸[βˆ£π‘†π‘‡ ∣] = 2𝐸[𝑆𝑇+ ].
𝑆0 = 0
and then
Let us rst clarify the meaning of (i) under the
given conditions. The assumed convergence implies in particular that
[∫
𝑇
]
∫
βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…0 (π‘₯)∣ 𝐹 (𝑑π‘₯)𝑑𝑑 < ∞
𝐸
0
for some
𝑇 > 0.
(4.3)
ℝ
Due to this fact, the following are actually equivalent to (i):
(a)
𝑍
is of integrable variation,
(b)
𝑆
is of integrable variation on
[0, 𝑇 ]
for some
𝑇 > 0.
𝑆 = 𝑆 1 βˆ’ 𝑆 2 on [0, 𝑇 ] for two increasing processes
=
= 0. Since the jumps of the martingale 𝑆 are
1
2
integrable, the positive stopping time 𝜏 := inf{0 ≀ 𝑑 ≀ 𝑇 : 𝑆𝑑 + 𝑆𝑑 β‰₯ 1}
is such that 𝑆 is of integrable variation on [[0, 𝜏 ]]. By [20, II.1.33b] this im[∫𝜏 ∫
]
𝐹 (𝑑π‘₯)𝑑𝑑 < ∞. Now we can use that the product
plies that 𝐸
0 ℝ βˆ£πœ…βˆ«π‘‘ (π‘₯)∣
∫
∫
𝜏
𝜏 ℝ βˆ£πœ…0 (π‘₯)∣ 𝐹 (𝑑π‘₯) = 0 ℝ βˆ£πœ…0 (π‘₯)∣ 𝐹 (𝑑π‘₯)𝑑𝑑 is bounded by
∫ 𝜏∫
∫ 𝜏∫
βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…0 (π‘₯)∣ 𝐹 (𝑑π‘₯)𝑑𝑑 +
βˆ£πœ…π‘‘ (π‘₯)∣ 𝐹 (𝑑π‘₯)𝑑𝑑 < ∞ 𝑃 -a.s.
Indeed, assume (i). Then
𝑆 1 and
0
𝑆 2 with
𝑆01
𝑆02
0
ℝ
to conclude via Lemma A.1(vii) that
more, (a) together with (4.3) yields
𝐸
ℝ
𝑍 is of integrable variation. Further[βˆ«π‘‡ ∫
]
0 ℝ βˆ£πœ…π‘‘ (π‘₯)∣ 𝐹 (𝑑π‘₯)𝑑𝑑 < ∞, which is
(b) by [20, II.1.33b], and clearly (b) implies (i).
Step 1: Lévy case.
We rst assume that
the equivalent condition (a).
10
𝑆 = 𝑍.
Moreover, we replace (i) by
(ii) implies (a):
𝑓𝑛
on
ℝ
We consider an increasing sequence of continuous functions
satisfying
0 ≀ 𝑓𝑛 (π‘₯) ≀ ∣π‘₯∣ ∧ 𝑛,
for all
π‘₯βˆˆβ„
and
𝑛 β‰₯ 1.
𝑓𝑛 (π‘₯) = 0
For each
𝑛
for
∣π‘₯∣ < 1/𝑛,
lim 𝑓𝑛 (π‘₯) = ∣π‘₯∣
𝑛
we have
∫
1
1
lim inf 𝐸[βˆ£π‘π‘‡ ∣] β‰₯ lim inf 𝐸[𝑓𝑛 (𝑍𝑇 )] =
𝑇 ↓0 𝑇
𝑇 ↓0 𝑇
𝑓𝑛 (πœ…0 (π‘₯)) 𝐹 (𝑑π‘₯)
ℝ
where the equality follows from Sato [26, Corollary 8.9] since
𝑓𝑛
vanishes in a
neighborhood of the origin (this holds for any Lévy measure). By monotone
convergence as
π‘›β†’βˆž
we obtain
1
lim inf 𝐸[βˆ£π‘π‘‡ ∣] β‰₯
𝑇 ↓0 𝑇
∫
βˆ£πœ…0 (π‘₯)∣ 𝐹 (𝑑π‘₯).
The left hand side is nite by assumption, hence
(a) implies ∫(ii):
𝑑
If
𝑍
𝑍
is of integrable variation.
is of integrable variation, its total variation process
βˆ£π‘‘π‘π‘  ∣ is an integrable Lévy subordinator and by Lemma A.1(vii)
∫
∫
1
𝐸[Var(𝑍)𝑇 ] = 𝐸[Var(𝑍)1 ] = βˆ£πœ…0 (π‘₯)∣ 𝐹 (𝑑π‘₯) + πœ…0 (π‘₯) 𝐹 (𝑑π‘₯) = 𝐢.
𝑇
Var(𝑍)𝑑 :=
Since
0
βˆ£π‘π‘‡ ∣ ≀ Var(𝑍)𝑇 ,
we conclude that
lim sup
𝑇 ↓0
1
𝐸[βˆ£π‘π‘‡ ∣] ≀ 𝐢.
𝑇
For the converse inequality we consider the function
Since
𝑍
is of nite variation, we obtain for each
1
1
lim inf 𝐸[βˆ£π‘π‘‡ ∣] β‰₯ lim 𝐸[𝑔𝑛 (𝑍𝑇 )] =
𝑇 ↓0 𝑇
𝑇 ↓0 𝑇
∫
𝑛
𝑔𝑛 (π‘₯) = ∣π‘₯βˆ£βˆ§π‘› for 𝑛 β‰₯ 1.
that
∫
𝑔𝑛 (πœ…0 (π‘₯)) 𝐹 (𝑑π‘₯)+ πœ…0 (π‘₯) 𝐹 (𝑑π‘₯)
𝑍
are i.i.d.
(see also [19, Equation (5.8)]). Applying monotone convergence as
π‘›β†’βˆž
as a consequence of [19, Theorem 2.1(i)(c)] since the increments of
to the right hand side, we conclude that
lim inf
𝑇 ↓0
1
𝐸[βˆ£π‘π‘‡ ∣] β‰₯ 𝐢.
𝑇
Hence we have proved the claimed convergence rate for the case
Step 2: General case.
𝑆 = 𝑍.
Under the stated assumption, Proposition 2.1(ii) with
𝛽 = 1 and 𝛾 = 0 yields that 𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑍𝑇 ∣] = π‘œ(𝑇 ) as 𝑇 ↓ 0. Hence (ii) holds
𝑆 if and only if it holds for 𝑍 and so Step 1 yields the equivalence of (i)
and (ii). Moreover, the leading constant 𝐢 is the same as for 𝑍 since the
approximation error is of order π‘œ(𝑇 ).
for
11
Remark 4.2.
(a) For exponential Lévy processes, where
πœ…(π‘₯) = 𝑒π‘₯ βˆ’ 1,
the formula in (ii) was previously obtained by Tankov [27]. We thank
the referee for pointing out this reference.
𝐢 has the following feature. For a given absolute moment
βˆ£πœ…0 (π‘₯)∣ 𝐹 (𝑑π‘₯) of the Lévy measure
of 𝑍 , the value of 𝐢 may
∫
from πœ‡ to 2πœ‡ depending on
πœ…0 (π‘₯) 𝐹 (𝑑π‘₯). In particular, 𝐢 is
(b) The constant
πœ‡ :=
range
∫
minimal if the jumps are symmetric and maximal if all jumps have the
same sign.
(c) As in Corollary 3.3, it follows from Theorem 4.1 that the implied
volatility satises
πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) ∼
√
√
πœ‹/2 𝐢/𝑆0 𝑇
as
𝑇 ↓ 0.
(d) The following observation from Step 1 in the proof seems worth be-
𝑍 is a Lévy martingale of nite variation, the three
functions 𝑇 7β†’ 𝐸[βˆ£π‘π‘‡ ∣], 𝑇 7β†’ 𝐸[sup𝑑≀𝑇 βˆ£π‘π‘‘ ∣] and 𝑇 7β†’ 𝐸[Var(𝑍)𝑇 ] all
converge to zero as 𝑇 ↓ 0 with the same leading term 𝐢𝑇 .
ing recorded: if
4.2 Innite Variation
We now turn to pure-jump processes with innite variation. By the previous
results we know that the leading order for the at-the-money option price has
√
to be strictly between
𝑇
and
𝑇;
however, it need not be a power of
𝑇.
The
class of possible Lévy measures at time zero is very rich and it is unclear
how to compute the exact order in general.
A look at existing nancial
models suggests to impose additional structure which will pin down an order
of parametric form.
The base case is the
𝛼-stable
Lévy measure which is
given by
𝜈(𝑑π‘₯) =
for
𝜈
𝛼 ∈ (0, 2)
𝑔(π‘₯)
𝑑π‘₯,
∣π‘₯∣1+𝛼
𝑔(π‘₯) := π›½βˆ’ 1(βˆ’βˆž,0) (π‘₯) + 𝛽+ 1(0,∞) (π‘₯)
and two nonnegative constants
𝛽+
and
π›½βˆ’ .
More precisely,
𝐿𝑑 following a stable law with index of
𝛽 = (𝛽+ βˆ’ π›½βˆ’ )/(𝛽+ + π›½βˆ’ ), shift parameter
πœ‡ = 0 and scale parameter 𝑐 = (𝛽+ + π›½βˆ’ )1/𝛼 𝑑1/𝛼 . In the nondegenerate case
𝛽+ + π›½βˆ’ > 0 the process 𝐿 has innite variation if and only if 𝛼 ∈ [1, 2). For
𝛼 ∈ (1, 2) the rst moment 𝐸[βˆ£πΏπ‘‘ ∣] exists whereas for 𝛼 ∈ (0, 1] this is not
the case and in particular 𝐿 cannot be a martingale. On the other hand,
corresponds to a Lévy process
stability
𝛼,
skewness parameter
only the small jumps are relevant for the option price asymptotics of order
smaller than
𝑂(𝑇 ),
which leads us to the following denition.
Denition 4.3. Let 𝛼+ , π›Όβˆ’ ∈ (0, 2). A Lévy process is said to
(𝛼+ , π›Όβˆ’ )-stable-like small jumps if its Lévy measure 𝜈 is of the form
(
𝜈(𝑑π‘₯) =
)
𝑓 (π‘₯)
𝑓 (π‘₯)
1
(π‘₯) + 1+𝛼+ 1(0,∞) (π‘₯) 𝑑π‘₯
∣π‘₯∣1+π›Όβˆ’ (βˆ’βˆž,0)
∣π‘₯∣
12
have
for a Borel function
𝑓 β‰₯0
whose left and right limits at zero,
𝑓+ := lim 𝑓 (π‘₯)
π‘₯↓0
exist and satisfy
𝑓 (π‘₯) βˆ’ 𝑓+ = 𝑂(π‘₯)
and
as
π‘“βˆ’ := lim 𝑓 (π‘₯),
π‘₯↓0
π‘₯↑0
and
𝑓 (π‘₯) βˆ’ π‘“βˆ’ = 𝑂(π‘₯)
as
π‘₯ ↑ 0.
This class includes most of the processes used in nancial modeling; e.g.,
the tempered stable and in particular the CGMY processes, of which the
variance gamma process is a special case, or the normal inverse Gaussian
process. We refer to Cont and Tankov [8, Section 4.5] for more information
on these models. We observe that if the driving Lévy process of an SDE as
in Corollary 2.2 is chosen from this class, then the process
𝑍
dened in the
same corollary is again of the same type; merely the constants
𝑓+
and
π‘“βˆ’
change.
Since
𝛼+ ∨ π›Όβˆ’ ∈ (0, 1)
implies that the jumps are of nite variation
(cf. Theorem 4.1), we are interested here only in the case
In that case we shall see that the larger of the values
𝛼+
𝛼+ ∨ π›Όβˆ’ ∈ [1, 2).
π›Όβˆ’ determines
and
the leading order.
The statement of our main result requires the following constants. For
𝑓+ , π‘“βˆ’ β‰₯ 0
and
𝛼 ∈ (1, 2)
we set
𝐢(𝛼, 0, 0) := 0
and if
𝑓+ + π‘“βˆ’ > 0
then
]1
(
)[
(
)
( π›Όπœ‹ ) 2𝛼
1
2
𝑓+ βˆ’ π‘“βˆ’ 2
1
𝐢(𝛼, 𝑓+ , π‘“βˆ’ ) := (𝑓+ + π‘“βˆ’ ) 𝛼 Ξ“ 1 βˆ’
1+
tan2
πœ‹
𝛼
𝑓+ + π‘“βˆ’
2
(
(
))
(
)
1
𝑓+ βˆ’ π‘“βˆ’
π›Όπœ‹
× cos
arctan
tan
,
𝛼
𝑓+ + π‘“βˆ’
2
Ξ“ denotes
𝛼+ ∨ π›Όβˆ’ ∈ (1, 2)
where
the usual Gamma function. For
𝛼+ , π›Όβˆ’ ∈ (0, 2)
such that
we then dene
⎧

⎨𝐢(𝛼+ , 𝑓+ , π‘“βˆ’ ), 𝛼+ = π›Όβˆ’ ,
𝐢(𝛼+ , π›Όβˆ’ , 𝑓+ , π‘“βˆ’ ) := 𝐢(𝛼+ , 𝑓+ , 0),
𝛼+ > π›Όβˆ’ ,

⎩
𝐢(π›Όβˆ’ , 0, π‘“βˆ’ ),
𝛼+ < π›Όβˆ’ .
Theorem 4.4. Let 𝑆 be a pure-jump martingale of the form (4.1) such that
the Lévy process 𝑍 from (4.2) has (𝛼+ , π›Όβˆ’ )-stable-like small jumps and let
𝛼 := 𝛼+ ∨ π›Όβˆ’ . Assume that there exist 𝛽 ∈ [1, 2] and 𝛾 β‰₯ 0 such that
[∫
𝐸
]
βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…0 (π‘₯)βˆ£π›½ 𝐹 (𝑑π‘₯) = π‘œ(𝑑𝛾 )
ℝ
and
𝛽
≀ 𝛼.
1+𝛾
(i) If 𝛼 ∈ (1, 2), then
)
1 (
𝐸[(𝑆𝑇 βˆ’ 𝑆0 )+ ] = 𝐢 𝛼+ , π›Όβˆ’ , 𝑓+ , π‘“βˆ’ 𝑇 1/𝛼 + π‘œ(𝑇 1/𝛼 )
2
13
as 𝑇 ↓ 0.
(ii) If 𝛼+ = π›Όβˆ’ = 1 and 𝑓+ = π‘“βˆ’ , then
𝐸[(𝑆𝑇 βˆ’ 𝑆0 )+ ] =
Remark 4.5.
)
1(
𝑓+ + π‘“βˆ’ 𝑇 ∣ log 𝑇 ∣ + π‘œ(𝑇 ∣ log 𝑇 ∣)
2
as 𝑇 ↓ 0.
(a) For exponential Lévy processes, a result similar to part
(i) is obtained in [27]. There, the condition on the jumps in formulated
in a slightly dierent way and hence our theorem is not a strict generalization. More specically, in [27], the author deals with Lévy processes
whose characteristic function resembles the one of a stable process. On
the other hand, we consider processes whose jump measures resemble
the jump measure of a stable process around the origin, which seems
more convenient beyond Lévy models. For a stable process, one readily
veries that the two expressions for the small-time limit indeed coincide, once the relationship between the characteristic function and the
jump measure (cf., e.g., [25, Chapter I.1]) is taken into account.
The special case of a CGMY model is also treated in [13].
πœ…π‘‘ (π‘₯) is satised in particular
∫
𝐸[ ℝ βˆ£πœ…π‘‘ (π‘₯) βˆ’ πœ…0 (π‘₯)∣2 𝐹 (𝑑π‘₯)] = 𝑂(𝑑), as 𝛼 > 1. For example, this
(b) In part (i), the continuity assumption on
if
holds in the setting of the Lévy-driven SDE of Corollary 2.2 (see the
proof of that result).
√ 𝛼 ↑ 2, which corresponds to the Brownian case, we obtain
𝑇 as in Theorem 3.1. On the other hand, the limit 𝛼 ↓ 1 in
(c) In the limit
the order
part (i) does not yield the order obtained in part (ii) and the leading
constants explode since
limπ‘Žβ†“0 Ξ“(π‘Ž) = +∞.
𝑓 in Denition 4.3 is
of 𝑓 (π‘₯) βˆ’ 𝑓± = 𝑂(∣π‘₯∣), it is sucient to
some 𝜚 > 𝛼/2. The proof is identical.
(d) The theorem still holds true if the regularity of
weakened as follows: Instead
have
𝑓 (π‘₯) βˆ’ 𝑓± = 𝑂(∣π‘₯∣𝜚 )
for
(e) As in Corollary 3.3 we can deduce that the implied volatility satises
√
1/π›Όβˆ’1/2 in the setting of part
πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) ∼ πœ‹/2 𝐢(𝛼
√+ , π›Όβˆ’ , 𝑓+ , π‘“βˆ’ )/𝑆0 𝑇 √
(i) and πœŽπ‘–π‘šπ‘π‘™ (𝑇 ) ∼
πœ‹/2 (𝑓+ + π‘“βˆ’ )/𝑆0 𝑇 ∣ log 𝑇 ∣ in the setting of
part (ii) of Theorem 4.4.
(f ) The following consequence of a result due to Luschgy and Pagès (see
𝑆 is a (1, 1)-stable-like Lévy
𝐸[(𝑆𝑇 βˆ’ 𝑆0 )+ ] = 𝑂(𝑇 ∣ log 𝑇 ∣)
[21, Theorem 3]) complements part (ii): If
process, possibly with
𝑓+ βˆ•= π‘“βˆ’ ,
then
still holds. However, their method only yields an upper bound and not
∼ 𝐢𝑇 ∣ log 𝑇 ∣. As in the proof below
we can infer that the same bound holds if 𝑆 is not a Lévy process but
satises the assumptions of part (ii) excluding 𝑓+ = π‘“βˆ’ .
that the leading order is indeed
4.2.1 Proof of Theorem 4.4
The plan for the proof of Theorem 4.4 is as follows. We shall again reduce
from the general martingale to the Lévy case by Proposition 2.1.
14
Since
the rst absolute moment is known for stable processes, the main step for
part (i) will be to estimate the error made when replacing a stable-like Lévy
𝛼 = 𝛼+ ∨ π›Όβˆ’ > 1. For part (ii)
1-stable process fails to have a rst
process by a true stable process with index
the situation is slightly dierent as the
moment; in this case we shall instead use the normal inverse Gaussian as a
reference process.
In a rst step we show more generally that we can pass with an error
of order
𝑂(𝑇 )
from one Lévy process to another one when the small jumps
have a similar behavior.
Lemma 4.6. Let 𝐿 and 𝐿′ be pure-jump Lévy martingales with Lévy measures 𝜈 and 𝜈 β€² , respectively. Suppose that for some 𝛿 > 0, the Radon-Nikodym
derivative
( β€²
)
πœ“(π‘₯) =
𝑑 𝜈 (𝑑π‘₯)∣[βˆ’π›Ώ,𝛿]
(
)
𝑑 𝜈(𝑑π‘₯)∣[βˆ’π›Ώ,𝛿]
of the measures restricted to [βˆ’π›Ώ, 𝛿] exists and that
Then 𝐸[βˆ£πΏπ‘‡ ∣] = 𝐸[βˆ£πΏβ€²π‘‡ ∣] + 𝑂(𝑇 ) as 𝑇 ↓ 0.
Proof.
πœ‡
Let
𝐿 = π‘₯ βˆ— (πœ‡ βˆ’ 𝜈(𝑑π‘₯)𝑑𝑑)
βˆ«π›Ώ
βˆ’π›Ώ
βˆ£πœ“(π‘₯)βˆ’1∣2 𝜈(𝑑π‘₯) < ∞.
be the canonical representation of
is the random measure associated with the jumps of
into
𝐿=
𝐿≀𝛿
+
𝐿.
𝐿;
here
We decompose
𝐿
𝐿>𝛿 , where the Lévy process
𝐿≀𝛿 := π‘₯1{∣π‘₯βˆ£β‰€π›Ώ} βˆ— (πœ‡ βˆ’ 𝜈(𝑑π‘₯)𝑑𝑑)
is obtained by truncating the jumps at magnitude
𝛿.
Then
𝐿>𝛿 = 𝐿 βˆ’ 𝐿≀𝛿
is of nite variation and Theorem 4.1 yields
𝐸[∣𝐿>𝛿
𝑇 ∣] = 𝑂(𝑇 )
Hence we may assume that
𝛿
𝐿 = 𝐿≀𝛿 ,
as
𝑇 ↓ 0.
i.e., that the jumps are bounded by
𝜈 is concentrated on [βˆ’π›Ώ, 𝛿]. The
πœ“ ensures that π‘Œ := (πœ“(π‘₯)βˆ’1)βˆ—(πœ‡βˆ’πœˆ(𝑑π‘₯)𝑑𝑑) is a
martingale; moreover, πœ“ β‰₯ 0 implies that Ξ”π‘Œ β‰₯ βˆ’1.
in absolute value, or equivalently that
integrability assumption on
square-integrable Lévy
Hence the stochastic exponential
(
)
𝐷 := β„°(π‘Œ ) = β„° (πœ“(π‘₯) βˆ’ 1) βˆ— (πœ‡ βˆ’ 𝜈(𝑑π‘₯)𝑑𝑑)
is a nonnegative square-integrable martingale (cf. Lemma A.1(x)). We dene
the probability measure
𝑄β‰ͺ𝑃
on
β„±1
by
𝑑𝑄/𝑑𝑃 = 𝐷1 .
Then the Girsanov-
Jacod-Mémin theorem (cf. [20, III.3.24]) yields that under
(𝐿𝑑 )0≀𝑑≀1
β„Ž(π‘₯) = π‘₯,
the process
(𝑏𝑄 , 0, πœˆπ‘„ ) relative to the truncation
πœˆπ‘„ (𝑑π‘₯) = πœ“(π‘₯) 𝜈(𝑑π‘₯) = 1[βˆ’π›Ώ,𝛿] 𝜈 β€² (𝑑π‘₯) and
is Lévy with triplet
where
𝑄,
∫
𝛿
π‘₯(πœ“(π‘₯) βˆ’ 1) 𝜈(𝑑π‘₯);
𝑏𝑄 =
βˆ’π›Ώ
15
function
this integral is nite due to the assumption on
πœ“
and Hölder's inequality.
𝑂(𝑇 ), the 𝑄-distribution
𝑃 -distribution of 𝐿′≀𝛿 , which is obtained
β€²
′≀𝛿 .
As before we may assume that 𝐿 = 𝐿
After subtracting the linear drift which is of order
of
𝐿
from
therefore coincides with the
𝐿′
by truncating the jumps.
To summarize, we have
𝐸𝑄 [βˆ£πΏπ‘‡ ∣] = 𝐸[βˆ£πΏβ€²π‘‡ ∣] + 𝑂(𝑇 )
and hence it suces to show that
𝐸𝑄 [βˆ£πΏπ‘‡ ∣] βˆ’ 𝐸[βˆ£πΏπ‘‡ ∣] = 𝑂(𝑇 ).
Indeed,
𝐸𝑄 [βˆ£πΏπ‘‡ ∣] = 𝐸[𝐷𝑇 βˆ£πΏπ‘‡ ∣], Hölder's inequality and Lemma A.1(vi),(x) yield
𝐸𝑄 [βˆ£πΏπ‘‡ ∣] βˆ’ 𝐸[βˆ£πΏπ‘‡ ∣] ≀ 𝐸[βˆ£π·π‘‡ βˆ’ 1βˆ£βˆ£πΏπ‘‡ ∣]
{
}1/2
≀ 𝐸[(𝐷𝑇 βˆ’ 1)2 ] 𝐸[𝐿2𝑇 ]
{(
)
}1/2
= 𝐸[𝐷𝑇2 ] βˆ’ 1 𝐸[𝐿2𝑇 ]
{(
)
}1/2
= 𝑒𝑇 βŸ¨π‘Œ,π‘Œ ⟩1 βˆ’ 1 𝑇 ⟨𝐿, 𝐿⟩1
= 𝑂(𝑇 ),
where we have used that both
Proof of Theorem 4.4(i).
that
𝑆0 = 0
Step 1: 𝛼-stable Lévy case.
additional assumption that
Then
𝑍
and
𝐿 = 𝐿≀𝛿
are square-integrable.
As in the proof of Theorem 3.1 we may assume
and then we have
𝛼+ = π›Όβˆ’ =: 𝛼.
π‘Œ
𝐸[βˆ£π‘†π‘‡ ∣] = 2𝐸[𝑆𝑇+ ].
We rst note the result for
𝑍
and under the
𝑓 (π‘₯) = π‘“βˆ’ 1(βˆ’βˆž,0) (π‘₯) + 𝑓+ 1(0,∞) (π‘₯) and that
𝛼-stable Lévy motion with 𝛼 ∈ (1, 2)
is a centered
and in this case it is known that
𝐸[βˆ£π‘π‘‡ ∣] = 𝐢(𝛼, 𝑓+ , π‘“βˆ’ ) 𝑇 1/𝛼 ;
see Samorodnitsky and Taqqu [25, Property 1.2.17].
Step 2: Stable-like Lévy case with 𝛼+ = π›Όβˆ’ .
We again consider
𝑍
for the
𝛼+ = π›Όβˆ’ =: 𝛼 but now let 𝑓 be an arbitrary function satisfying
𝑓 (π‘₯) βˆ’ 𝑓+ = 𝑂(π‘₯) as π‘₯ ↑ 0 and 𝑓 (π‘₯) βˆ’ π‘“βˆ’ = 𝑂(π‘₯) as π‘₯ ↓ 0. For 𝛿 > 0 we
special case
dene the function
(
π‘“βˆ’
πœ“(π‘₯) := 1[βˆ’π›Ώ,0] (π‘₯)
+ 1{π‘“βˆ’ =0}
𝑓 (π‘₯)
where we use the convention
consider the case
𝑓 (π‘₯) β‰₯ 𝑓+ /2 > 0
𝑓+ > 0.
0/0 := 0.
)
(
)
𝑓+
+ 1(0,𝛿] (π‘₯)
+ 1{𝑓+ =0} ,
𝑓 (π‘₯)
We have
Then for
π‘₯ > 0
βˆ£πœ“(π‘₯) βˆ’ 1∣ = 𝑂(∣π‘₯∣).
small enough we have that
and thus
𝑓+ βˆ’ 𝑓 (π‘₯) ≀ 2 βˆ£π‘“+ βˆ’ 𝑓 (π‘₯)∣ = 𝑂(π‘₯)
βˆ£πœ“(π‘₯) βˆ’ 1∣ = 𝑓 (π‘₯) 𝑓+
The case
Indeed,
𝑓+ = 0 is trivial and π‘₯ < 0 is treated in the same
𝛿 small enough we can nd a constant 𝑀 > 0
by choosing
βˆ£πœ“(π‘₯) βˆ’ 1∣ ≀ 𝑀 ∣π‘₯∣,
16
π‘₯ ∈ [βˆ’π›Ώ, 𝛿].
as
π‘₯ ↓ 0.
way. As a result,
such that
πœ“
In particular,
satises the integrability assumption of Lemma 4.6. In the
sequel, we denote by
𝜈
the Lévy measure of
Case 2a: 𝑓+ > 0 and π‘“βˆ’ > 0.
πœˆβ€²
price only at the order
𝜈(β‹…) = 𝐹 (πœ…βˆ’1
0 (β‹…)).
In this case we have
𝛼-stable
is the Lévy measure of an
𝛿.
were truncated at magnitude
i.e.,
(
)
π‘“βˆ’
𝑓+
1[βˆ’π›Ώ,0) (π‘₯) 1+𝛼 + 1(0,𝛿] 1+𝛼 (π‘₯) 𝑑π‘₯.
∣π‘₯∣
∣π‘₯∣
𝜈 β€² (𝑑π‘₯) := πœ“(π‘₯) 𝜈(𝑑π‘₯) =
Note that
𝑍,
𝑂(𝑇 ).
Lévy motion whose jumps
As above, this truncation changes the option
Now Lemma 4.6 and Step 1 yield that
𝐸[βˆ£π‘π‘‡ ∣] = 𝐢(𝛼, 𝑓+ , π‘“βˆ’ ) 𝑇 1/𝛼 + 𝑂(𝑇 ).
Case 2b: 𝑓+ = 0 and π‘“βˆ’ > 0.
(4.4)
Note that in this case the formula for
the desired one in general, since we have
𝑓 (π‘₯)
instead of
𝑓+
πœˆβ€²
is not
in
(
)
π‘“βˆ’
𝑓 (π‘₯)
𝜈 (𝑑π‘₯) = πœ“(π‘₯) 𝜈(𝑑π‘₯) = 1[βˆ’π›Ώ,0) (π‘₯) 1+𝛼 + 1(0,𝛿] 1+𝛼 (π‘₯) 𝑑π‘₯.
∣π‘₯∣
∣π‘₯∣
β€²
However, the previous argument does apply if
and we shall reduce to this case.
that the positive jumps of
𝑍
𝑓 (π‘₯) = 𝑓+ = 0
𝑓 (π‘₯) = 𝑓+ = 0
π‘₯ > 0,
Indeed, Lemma 4.7 stated below shows
are of nite variation, hence of integrable vari-
ation by (3.3). By subtracting these jumps from
achieve that
for all
for all
π‘₯>0
𝑍
and compensating, we
and Theorem 4.1 shows that this
manipulation aects the option price only at the order
𝑂(𝑇 ).
Hence we may
conclude as in Case 2a to obtain (4.4).
The case where
𝑓+ > 0
π‘“βˆ’ = 0 is analogous. Finally, if 𝑓+ and π‘“βˆ’
𝑍 is of nite variation and Theorem 4.1
of order 𝑂(𝑇 ). Hence (4.4) again holds since
and
both vanish, Lemma 4.7 shows that
yields that the option price is
𝐢(𝛼, 0, 0) = 0.
Step 3: Stable-like Lévy case with 𝛼+ βˆ•= π›Όβˆ’ . We consider the case where
𝛼 := 𝛼+ > π›Όβˆ’ . The idea is to reduce to the case where 𝛼+ = π›Όβˆ’ but
π‘“βˆ’ = 0; i.e., we get rid of all the negative jumps and show that this induces
1/𝛼 ).
an error of order π‘œ(𝑇
Let πœ‡ be the random measure associated with the jumps of 𝑍 . By setting
𝑍 + := π‘₯+ βˆ— (πœ‡ βˆ’ 𝜈(𝑑π‘₯) 𝑑𝑑)
we decompose
and
𝑍 = 𝑍 + +𝑍 βˆ’ , where 𝑍 +
𝑍 βˆ’ := π‘₯βˆ’ βˆ— (πœ‡ βˆ’ 𝜈(𝑑π‘₯) 𝑑𝑑)
and
π‘βˆ’
are Lévy martingales having
only positive and negative jumps, respectively, and Lévy measures given by
𝜈 + (𝑑π‘₯) = 1(0,∞) (π‘₯) 𝜈(𝑑π‘₯) and 𝜈 βˆ’ (𝑑π‘₯) = 1(βˆ’βˆž,0) (π‘₯) 𝜈(𝑑π‘₯). Note the abuse of
+ and π‘₯βˆ’ refer to the positive and negative part of π‘₯ while 𝑍 ±
notation: π‘₯
± are new symbols.
and 𝜈
+ has (𝛼 , 𝛼 )-stable-like jumps, where the correWe observe that 𝑍
+
+
β€²
β€²
sponding function 𝑓 in Denition 4.3 satises 𝑓 (π‘₯) = 0 for π‘₯ < 0. In
17
particular the left limit is
erties for
π‘“βˆ’β€² = 0.
The martingale
π‘βˆ’
has analogous prop-
π›Όβˆ’ .
𝐸[βˆ£π‘π‘‡βˆ’ ∣] = 𝑂(𝑇 1/π›Όβˆ’ ) if π›Όβˆ’ ∈ (1, 2) and we
βˆ’
have also seen that 𝐸[βˆ£π‘π‘‡ ∣] = 𝑂(𝑇 ) if π›Όβˆ’ ∈ (0, 1) due to nite variation. For
βˆ’
the remaining case π›Όβˆ’ = 1 we have 𝐸[βˆ£π‘π‘‡ ∣] = 𝑂(𝑇 ∣ log 𝑇 ∣) by Remark 4.5(f ).
From Step 2 we know that
As
π›Όβˆ’ < 𝛼+
and
𝛼+ > 1
we therefore have
𝐸[βˆ£π‘π‘‡βˆ’ ∣] = π‘œ(𝑇 1/𝛼+ )
π›Όβˆ’ . On the other hand, we know from Step 2 applied
𝐸[βˆ£π‘π‘‡+ ∣] ∼ 𝐢𝑇 1/𝛼+ . Therefore the leading order coecient for 𝑍
+
same as for 𝑍 , i.e.,
(
)
𝐸[βˆ£π‘π‘‡ ∣] = 𝐢 𝛼+ , 𝑓+ , 0 𝑇 1/𝛼+ + π‘œ(𝑇 1/𝛼+ ).
in all three cases for
to
𝛼+
is the
that
The case
𝛼+ < π›Όβˆ’
is analogous.
Step 4: General case.
Proposition 2.1(ii) implies that
𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑍𝑇 ∣] = π‘œ(𝑇 (1+𝛾)/𝛽 )
and in particular
𝐸[βˆ£π‘†π‘‡ βˆ’ 𝑍𝑇 ∣] = π‘œ(𝑇 1/𝛼 ).
In view of the previous steps, this
completes the proof.
The following result was used in the preceding proof.
Lemma 4.7. Let 𝐿 be a Lévy process with (𝛼+ , π›Όβˆ’ )-stable-like small jumps
for some 𝛼+ , π›Όβˆ’ ∈ (0, 2). If the function 𝑓 from Denition 4.3 satises
βˆ‘ 𝑓+ = 0,+then the positive jumps of 𝐿 are of nite variation, that is,
𝑑≀𝑇 (Δ𝐿𝑑 ) < ∞ for any 𝑇 < ∞.
Proof.
Let
π‘₯ ↓ 0,
there exist
𝜈
be the Lévy measure of
𝛿 > 0
and
𝐿.
𝑀 > 0
Since
𝑓 (π‘₯) = 𝑓 (π‘₯) βˆ’ 𝑓+ = 𝑂(π‘₯) as
𝑓 (π‘₯) ≀ 𝑀 π‘₯ for π‘₯ ∈ (0, 𝛿].
such that
Therefore,
∫
𝛿
∫
π‘₯ 𝜈(𝑑π‘₯) =
0
0
𝛿
π‘₯𝑓 (π‘₯)
𝑑π‘₯ ≀ 𝑀
∣π‘₯∣1+𝛼+
𝛿
∫
0
π‘₯2
𝑑π‘₯ < ∞,
∣π‘₯∣1+𝛼+
showing that the small positive jumps are summable. Of course, the large
jumps are always summable.
We now come to the proof of the second part of the theorem.
main dierence to the above is that we cannot use the
a reference since it is not integrable.
1-stable
The
process as
Instead, we use the normal inverse
Gaussian process. It is the symmetry of its Lévy density around zero that
forces us to impose the condition
method.
𝑓+ = π‘“βˆ’
in the theorem to apply our
Indeed, we are not aware of a suitable process with suciently
asymmetric density for which the absolute moment asymptotics are known.
18
Proof of Theorem 4.4(ii). The proof has the same structure as for part (i).
Step 1: Normal inverse Gaussian case. First let 𝑍 be a symmetric normal
inverse Gaussian process with Lévy measure
𝜈 β€² (𝑑π‘₯) =
𝜌
𝐾1 (∣π‘₯∣),
πœ‹βˆ£π‘₯∣
𝜌 > 0 and πΎπœƒ denotes the modied Bessel function of the third kind of
πœƒ. Then we are in the setting of Theorem 4.4(ii) with 𝑓+ = π‘“βˆ’ = 𝜌/πœ‹
by the properties of 𝐾1 ; see, e.g., [1, Formula (9.6.9)]. The absolute moments
of 𝑍𝑑 were calculated explicitly and for all 𝑑 by Barndor-Nielsen and Stelzer
where
order
[4, Corollary 4]. By their formula and another property of Bessel functions
(see [1, Formula (9.6.8)]) we have
𝐸[βˆ£π‘π‘‡ ∣] =
2𝜌 πœŒπ‘‡
2𝜌
𝑒 𝑇 𝐾0 (πœŒπ‘‡ ) ∼
𝑇 ∣ log 𝑇 ∣ = (𝑓+ + π‘“βˆ’ ) 𝑇 ∣ log 𝑇 ∣.
πœ‹
πœ‹
Step 2: General Lévy case.
𝑍 be a (1, 1)-stable-like Lévy martingale
2
with Lévy measure 𝜈(𝑑π‘₯) = 𝑓 (π‘₯)/∣π‘₯∣ 𝑑π‘₯, where by assumption 𝑓 satises
𝑓0 := limπ‘₯β†’0 𝑓 (π‘₯) = 𝑓+ = π‘“βˆ’ . As the case 𝑓0 = 0 again follows from
Lemma 4.7, we may assume that 𝑓0 > 0. Then there exists a small 𝛿 > 0
such that 𝑓 is bounded away from zero on [βˆ’π›Ώ, 𝛿] and we can dene
Now let
πœ“(π‘₯) := 1[βˆ’π›Ώ,𝛿] (π‘₯)
𝑓0 𝐾1 (∣π‘₯∣)∣π‘₯∣
.
𝑓 (π‘₯)
𝜌 := πœ‹π‘“0 we have 𝜈 β€² (𝑑π‘₯) = πœ“(π‘₯) 𝜈(𝑑π‘₯) on [βˆ’π›Ώ, 𝛿], where 𝜈 β€² is as
As 𝑓0 𝐾1 (∣π‘₯∣)∣π‘₯∣ = 𝑓0 + 𝑂(∣π‘₯∣) by [1, Formula (9.6.11)], we have
By choosing
in Step 1.
πœ“(π‘₯) βˆ’ 1[βˆ’π›Ώ,𝛿] (π‘₯) = 1[βˆ’π›Ώ,𝛿] (π‘₯)
𝑓0 𝐾1 (∣π‘₯∣)∣π‘₯∣ βˆ’ 𝑓 (π‘₯)
= 𝑂(∣π‘₯∣)
𝑓 (π‘₯)
𝑓 (π‘₯) = 𝑓0 + 𝑂(∣π‘₯∣). Making 𝛿 > 0 smaller if
βˆ«π›Ώ
2
necessary, we conclude that
βˆ’π›Ώ (πœ“(π‘₯) βˆ’ 1) 𝜈(𝑑π‘₯) < ∞ and now the assertion
due to the assumption that
follows from Lemma 4.6 and Step 1. The last step to the general martingale
case is as in the proof of part (i).
A
Appendix
The following lemma collects some standard facts about Lévy processes that
are used throughout the text. A Lévy process
𝐿 is an adapted càdlàg process
𝐿0 = 0.
with independent and stationary increments and
Lemma A.1. Let 𝐿 be a Lévy process having triplet (𝑏, 𝑐, 𝜈) with respect to
the truncation function β„Ž(π‘₯) = π‘₯1∣π‘₯βˆ£β‰€1 .
∫
(i) ℝ 1 ∧ ∣π‘₯∣2 𝜈(𝑑π‘₯) < ∞.
19
(ii) There is a decomposition 𝐿𝑑 = πΏπ‘Žπ‘‘ + 𝐿𝑏𝑑 + 𝐿𝑐𝑑 + 𝐴𝑑 into independent Lévy
processes such that πΏπ‘Ž is a compound Poisson process,
𝐿𝑏 is a purely
√
discontinuous martingale with bounded jumps, 𝐿𝑐 = π‘π‘Š is a scaled
Brownian motion and 𝐴 ∈ ℝ.
∫
(iii) Let 𝑝 ∈ [1, ∞). Then ∣π‘₯∣>1 ∣π‘₯βˆ£π‘ 𝜈(𝑑π‘₯) < ∞ if and only if 𝐸[βˆ£πΏπ‘‘ βˆ£π‘ ] < ∞
for all 𝑑 β‰₯ 0. In particular, if 𝐿 has bounded jumps, or equivalently if
the support of 𝜈 is compact, then 𝐸[βˆ£πΏπ‘‘ βˆ£π‘ ] < ∞ for all 𝑝 ∈ [1, ∞).
∫
(iv) 𝐿 is a martingale if and only if the two conditions ∣π‘₯∣>1 ∣π‘₯∣ 𝜈(𝑑π‘₯) < ∞
∫
and 𝑏 + ∣π‘₯∣>1 π‘₯ 𝜈(𝑑π‘₯) = 0 hold.
(v) If 𝐿 is integrable, then 𝐸[𝐿𝑑 ] = 𝑑𝐸[𝐿1 ] and 𝐿𝑑 βˆ’ 𝑑𝐸[𝐿1 ] is a martingale.
(vi) If 𝐿 is a square-integrable martingale, then
(vii)
∫
𝐸[𝐿2𝑑 ] = 𝐸[[𝐿, 𝐿]𝑑 ] = ⟨𝐿, πΏβŸ©π‘‘ = π‘‘βŸ¨πΏ, 𝐿⟩1 = 𝑑𝑐 + 𝑑 ℝ ∣π‘₯∣2 𝜈(𝑑π‘₯) < ∞.
∫
If 𝑐 = 0, 𝐿 is of nite variation if and only if ∣π‘₯βˆ£β‰€1 ∣π‘₯∣ 𝜈(𝑑π‘₯) < ∞
∫
and of integrable variation if and ∫only if ℝ ∣π‘₯∣ 𝜈(𝑑π‘₯) < ∞. In that
𝑑
case the total variation
Var(𝐿)
∫
∫ 𝑑 := 0 βˆ£π‘‘πΏπ‘  ∣ is a Lévy process satisfying
𝐸[Var(𝐿)1 ] = ∣π‘₯∣ 𝜈(𝑑π‘₯) + ∣ π‘₯ 𝜈(𝑑π‘₯)∣.
(viii) For any 𝛿 > 0 there is a decomposition 𝐿 = 𝐿≀𝛿 + 𝐿>𝛿 into two
independent Lévy processes satisfying βˆ£Ξ”πΏβ‰€π›Ώ ∣ ≀ 𝛿 as well as βˆ£Ξ”πΏ>𝛿 ∣ > 𝛿
on the set {βˆ£Ξ”πΏ>𝛿 ∣} > 0. The corresponding Lévy measures are given
by 𝜈 ≀𝛿 (𝑑π‘₯) = 1[βˆ’π›Ώ,𝛿] (π‘₯) 𝜈(𝑑π‘₯) and 𝜈 >𝛿 (𝑑π‘₯) = 1β„βˆ–[βˆ’π›Ώ,𝛿] (π‘₯) 𝜈(𝑑π‘₯). If 𝐿>𝛿
has no Brownian component, it is a compound Poisson process with
drift.
(ix) If 𝐿 is a martingale, the stochastic exponential β„°(𝐿) is again a martingale.
(x) If 𝐿 is a square-integrable martingale, then so is β„°(𝐿) and moreover
𝐸[β„°(𝐿)2𝑑 ] = exp(π‘‘βŸ¨πΏ, 𝐿⟩1 ).
Proof.
(i)(viii) can be found in any advanced textbook about Lévy pro-
cesses; see, e.g., [26].
Statement (ix) is [8, Proposition 8.23].
One way to
deduce the formula in (x) is to use Yor's formula in
(
)
β„°(𝐿)2𝑑 = β„° 2𝐿 + [𝐿, 𝐿] βˆ’ ⟨𝐿, 𝐿⟩ + ⟨𝐿, 𝐿⟩ 𝑑
(
)
= β„° 2𝐿 + [𝐿, 𝐿] βˆ’ ⟨𝐿, 𝐿⟩ 𝑑 exp(π‘‘βŸ¨πΏ, 𝐿⟩1 ).
Noting that
2𝐿+[𝐿, 𝐿]βˆ’βŸ¨πΏ, 𝐿⟩ is a Lévy martingale, (ix) yields the result.
The following makes precise a remark from the introduction.
Remark A.2.
Let
𝑆
be any càdlàg martingale with absolutely continuous
predictable characteristics, then
𝑆
can be represented in the form (1.2).
Indeed, let
𝑑𝐡𝑑 = 𝑏𝑑 𝑑𝑑,
𝑑𝐢𝑑 = πœŽπ‘‘2 𝑑𝑑,
20
π‘‘πœˆπ‘‘ = 𝐾𝑑 (𝑑π‘₯) 𝑑𝑑
be the characteristics of
β„Ž(π‘₯) = π‘₯
𝑆
with respect to the trivial truncation function
(cf. [20, Chapter II] for background).
𝐡 = 0. Moreover, let 𝐹
𝐹 (ℝ) = ∞. Then there
exist a Brownian motion π‘Š and a Poisson random measure 𝑁 with compensator 𝐹 (𝑑π‘₯)𝑑𝑑 such that (1.2) holds. Moreover, πœ… and 𝐾 satisfy the relation
βˆ’1
𝐾𝑑 (𝐴) = 𝐹 (πœ…βˆ’1
denotes the preimage
𝑑 (𝐴)) for any Borel set 𝐴, where πœ…π‘‘
with respect to the spatial variable π‘₯. To be precise, the construction of π‘Š
and 𝑁 may necessitate an enlargement of the probability space, but this is
sible since
𝑆
The latter choice is pos-
is a martingale, which then implies
be any atomless
𝜎 -nite
measure on
ℝ
such that
harmless since we are interested only in distributional properties. We refer
to Jacod [18, Theorem 14.68(a)] for further details.
References
[1] M. Abramowitz and I. A. Stegun. Handbook of Mathematical Functions. U.S.
Government Printing Oce, Washington, D.C., 1964.
[2] E. Alòs, J. Léon, and J. Vives. On the short-time behavior of the implied
volatility for jump-diusion models with stochastic volatility. Finance Stoch.,
11(4):571589, 2007.
[3] L. B. G. Andersen and V. V. Piterbarg. Moment explosions in stochastic
volatility models. Finance Stoch., 11(1):2950, 2007.
[4] O. E. Barndor-Nielsen and R. Stelzer. Absolute moments of generalized
hyperbolic distributions and approximate scaling of normal inverse Gaussian
Lévy processes. Scand. J. Statist., 32(4):617637, 2005.
[5] H. Berestycki, J. Busca, and I. Florent. Asymptotics and calibration of local
volatility models. Quant. Finance, 2(1):6169, 2002.
[6] H. Berestycki, J. Busca, and I. Florent. Computing the implied volatility
in stochastic volatility models. Comm. Pure Appl. Math., 57(10):13521373,
2004.
[7] P. Carr and L. Wu. What type of process underlies options? A simple robust
test. J. Finance, 58(6):25812610, 2003.
[8] R. Cont and P. Tankov. Financial modelling with jump processes. Chapman
& Hall/CRC, Boca Raton, 2004.
[9] J. C. Cox, J. E. Ingersoll, and S. A. Ross. A theory of the term structure of
interest rates. Econometrica, 53(2):385407, 1985.
[10] V. Durrleman. Convergence of at-the-money implied volatilities to the spot
volatility. J. Appl. Probab., 45(2):542550, 2008.
[11] V. Durrleman. From implied to spot volatilities. Finance Stoch., 14(2):157
177, 2010.
[12] J. Feng, M. Forde, and J.-P. Fouque. Short-maturity asymptotics for a fast
mean-reverting Heston stochastic volatility model. SIAM J. Financial Math.,
1:126141, 2010.
[13] J. E. Figueroa-López and M. Forde. The small-maturity smile for exponential
Lévy models. Preprint, arXiv:1105.3180v1, 2011.
[14] M. Forde. Small-time asymptotics for a general local-stochastic volatility
model, using the heat kernel expansion. Preprint, 2009.
[15] M. Forde and A. Jacquier. Small-time asymptotics for implied volatility under
the Heston model. Int. J. Theor. Appl. Finance, 12(6):861876, 2009.
[16] M. Forde, A. Jacquier, and R. Lee. The small-time smile and term structure
of implied volatility under the Heston model. Preprint, 2011.
21
[17] J. Gatheral, E. P. Hsu, P. M. Laurence, C. Ouyang, and T.-H. Wang. Asymptotics of implied volatility in local volatility models. To appear in Math. Finance, 2011.
[18] J. Jacod. Calcul stochastique et problèmes de martingales, volume 714 of
Lecture Notes in Math. Springer, Berlin, 1979.
[19] J. Jacod. Asymptotic properties of power variations of Lévy processes. ESAIM
Probab. Stat., 11:173196, 2007.
[20] J. Jacod and A. N. Shiryaev. Limit Theorems for Stochastic Processes.
Springer, Berlin, 2nd edition, 2003.
[21] H. Luschgy and G. Pagès. Moment estimates for Lévy processes. Elect. Comm.
in Probab., 13:422434, 2008.
[22] A. Medvedev and O. Scaillet. Approximation and calibration of short-term
implied volatilities under jump-diusion stochastic volatility. Rev. Financial
Studies, 20(2):427459, 2007.
[23] P. Protter. Stochastic Integration and Dierential Equations. Springer, New
York, 2nd edition, version 2.1, 2005.
[24] M. Roper. Implied volatility explosions: European calls and implied volatilities
close to expiry in exponential Lévy models. Preprint, 2008.
[25] G. Samorodnitsky and M. Taqqu. Stable Non-Gaussian Random Processes.
Chapman & Hall, New York, 1994.
[26] K.-I. Sato. Lévy Processes and Innitely Divisible Distributions. Cambridge
University Press, Cambridge, 1999.
[27] P. Tankov. Pricing and hedging in exponential Lévy models: Review of recent
results. In Paris-Princeton Lectures on Mathematical Finance 2010, volume
2003 of Lecture Notes in Math., pages 319359, Springer, Berlin, 2010.
22