MX - Non Gaussian Merton Black Scholes Theory
MX - Non Gaussian Merton Black Scholes Theory
MX - Non Gaussian Merton Black Scholes Theory
NON-GAUSSIAN
MERTON-BLACK-SCHOLES
THEORY
Svetlana I. Boyarchenko
Sergei Z . Levendorskii
World Scientific
NON-GAUSSIAN
MERTON-BLACK-SCHOLES
THEORY
ADVANCED SERIES ON STATISTICAL SCIENCE &
APPLIED PROBABILITY
Published
Vol. 1: Random Walks of Infinitely Many Particles
by P. Revesz
Vol. 2: Ruin Probabilities
by S. Asmussen
Vol. 3: Essentials of Stochastic Finance: Facts, Models, Theory
by Albert N. Shiryaev
Vol. 4: Principles of Statistical Inference from a Neo-Fisherian Perspective
by L. Pace and A. Salvan
Vol. 5: Local Stereology
by Eva B. Vedel Jensen
Vol. 6: Elementary Stochastic Calculus — With Finance in View
by T. Mikosch
Vol. 7: Stochastic Methods in Hydrology: Rain, Landforms and Floods
eds. O. E. Barndorff-Nielsen et al.
Vol. 8: Statistical Experiments and Decisions: Asymptotic Theory
by A. N. Shiryaev and V. G. Spokoiny
Vol. 9: Non-Gaussian Merton-Black-Scholes Theory
by S. I. Boyarchenko and S. 2. Levendorskii
Advanced Series on
Applied Probability
NON-GAUSSIAN
MERTON-BLACK-SCHOLES
THEORY
Svetlana I. Boyarchenko
University of Texas at Austin, USA
Sergei Z . Levendorskii
Rostov State University of Economics, Russia
V f e World Scientific
wb New Jersey * London • Singapore • Hong Kong
Jersey-London-Sine
Published by
World Scientific Publishing Co. Pte. Ltd.
P O Box 128, Farrer Road, Singapore 912805
USA office: Suite IB, 1060 Main Street, River Edge, NJ 07661
UK office: 57 Shelton Street, Covent Garden, London WC2H 9HE
For photocopying of material in this volume, please pay a copying fee through the Copyright
Clearance Center, Inc., 222 Rosewood Drive, Danvers, MA 01923, USA. In this case permission to
photocopy is not required from the publisher.
ISBN 981-02-4944-6
vn
Vlll Preface
Preface vii
0.0.1 General notation x
Chapter 1 Introduction 1
1.1 The Gaussian Merton-Black-Scholes theory 1
1.1.1 Strong points 1
1.1.2 Drawbacks 2
1.1.2.1 Fat tails and anomalous skewness and kurtosis 3
1.1.2.2 The Black-Scholes market and formula . . . . 4
1.1.2.3 Volatility smile and volatility surface 5
1.1.3 Remedies for the MBS-theory 6
1.1.3.1 General remarks 6
1.1.3.2 Jump-diffusion models and stochastic volatil-
ity models 6
1.1.3.3 Levy processes 7
1.2 Regular Levy Processes of Exponential type 8
1.2.1 Characteristic function and exponent, Levy measure and
Levy-Khintchine formula 9
1.2.2 Definition of RLPE in ID 10
1.2.3 Infinitesimal generators of RLPE as PDO 11
1.3 Pricing of contingent claims 13
1.3.1 Discrete time models with a discrete space of states: No-
arbitrage and equivalent martingale measures 13
xi
XII Contents
Introduction
In this chapter, we list both strong and weak points of the Gaussian MBS-
theory, explain why it needs modification, describe informally the technique
used in the book and present briefly the main results of the book.
l
2 Introduction
1.1.2 Drawbacks
In the case of the MBS-theory, the motto: "Our vices are continuation of
our virtues" 1 should be rephrased in the reverse order: "Our virtues are
continuation of our vices". Both the simplicity and success of the MBS-
theory stem from the possibility of reduction to boundary value problems
for PDE or ODE, and this reduction becomes possible due to the choice of
the simplest class of stochastic processes - Gaussian processes - to model
the evolution of prices in financial and other markets. (The rival modern
martingale approach, which tries to remain in the realm of Probability
theory and avoid the usage of PDE and PDO whenever it can help it,
also capitalizes on the extraordinary nice properties of Gaussian processes).
In the basic Black-Scholes model, the price of a stock (or index, like the
Standard & Poor 500 Index) follows the Geometric Brownian motion: St =
expX t , where Xt is the Brownian motion. Then the probability density
PAt (%) of the increments Xt+At — Xt is given by
( \ 1 ( (x-^t)2\
PAt{x) exp (L1)
= 7^Tt {-^Ar)'
where p, and a2 are called the drift and volatility. As Eq. (1.1) shows, the
tails of the probability density decay faster than an exponenial function as
X -» ± 0 0 .
{{x-mf) _ ((x - m) 4 )
A3:- -3 ; «:- -4 .
gation, to buy the underlying asset for the specified price (the strike price),
K, at the specified expiry date, T. This contract is called the European call
option. Denote its price at time t < T, conditioned on the current price
(spot price), St, of the underlying asset by Fca,\\(St,t). The Black-Scholes
formula reads
and
0.194 ~
0.192 -
0.19 -
0.188 -
0.178 -
0.176 -
0.174^-
0.6 0.7 O.8 0.9 1 1.1 1.2 1.3 1.4
S/K
2
T h e last observation explains why the Gaussian MBS-theory performs much better for
options far from the expiry
8 Introduction
tails 3 have been used to model stock returns and price options: Variance
Gamma Processes (VGP), used by D. Madan with co-authors; Normal In-
verse Gaussian Processes (NIG), used by O.E. BarndorfF-Nielsen's group;
Hyperbolic Processes and Generalized Hyperbolic Processes (HP and GHP)
used by E. Eberlein's group; Truncated Levy Processes (TLP, constructed
by Koponen 4 , used by J.-P. Bouchaud and his group, and extended by Bo-
yarchenko and Levendorskii; and Normal Tempered Stable Levy Processes
(NTS Levy processes); we delegate the detailed discussion and references
to Chapter 3. As A.N. Shiryaev remarked, the name TLP was misleading,
and so we replace it with the name "KoBoL processes".
Processes of all the families listed above have been shown to fit better to
the dynamics of historic prices, and pricing formulas for European options,
based on these processes, also perform better than earlier models. Until
recently, almost no effective analytical formulas in more difficult situations
than the pricing of the European options and even simpler forwards and
futures have been known, though by now, there are many papers and a cou-
ple of books devoted to various aspects of the general theory of modelling,
pricing and hedging under Levy processes. The main goal of the book is to
partially fill in this gap, while remaining as close as possible to the initial
Merton-Black-Scholes framework.
When working with an empirical set of data, one specifies the type of the
process, that is, chooses a parameterized family of processes, and then
fits the parameters of the process to data. For a general theory like a
non-Gaussian analogue of the MBS-theory, too specific information is not
needed; in fact, unnecessary specification of the process makes the theory
harder to understand. There are two ways to describe the properties of the
process, which are sufficient but by no means necessary for our approach to
work, and the classes of Levy processes with exponentially decaying tails
used in empirical studies of Financial markets and described in Subsec-
tion 1.1.3 enjoy these properties. From the point of view of Probability
3
By semi-heavy or exponentially decaying tails we mean that the probability density
behaves, for x -> ±oo, as const • \x\p± e x p ( - c r ± | i | ) , for some p± € R and o± > 0.
4
Non-infinitely divisible truncated Levy distributions had been constructed earlier by
R.N. Mantegna
Regular Levy Processes of Exponential type 9
Theory, the most natural description comes in terms of the Levy density,
which can be visualized as the density of jumps of the process, but from the
point of view of PDO-theory and Analysis in general, it is more convenient
to work with the characteristic exponent of the Levy process.
Example 1.1 Let v € (0,2), and t/>(£) = c|£|". This is the characteristic
exponent of the stable Levy process of order v, and the Levy measure is
n(da?) = cV{-v)'l\x\-J-ldx.
5
For the list of basic definitions of Probability Theory, see Chapter 2
10 Introduction
"Almost" means that though processes of BM, NIG, HP, GHP, KoBoL and
NTS Levy families satisfy conditions of both definitions, VGP satisfies the
conditions of the first definition but not the second one, since the character-
istic exponent behaves like constln |£|, as £ —¥ oo. For pricing of contingent
claims of European type, the additional property Eq. (1.4) is not essential,
but it is needed to obtain effective explicit formulas for the factors in the
Wiener-Hopf factorization formula, which we need in the study of perpet-
ual American options and barrier options. This is the reason why we will
mainly use the second definition. The adjective "exponential" needs no
explanation, and "regular" indicates that from the analytical point of view,
RLPE is the most tractable subclass of Levy processes, if the Brownian
motion is not available (notice that BM is an RLPE). We will call v the or-
der of the process, A_ and A+ the steepness parameters, and c the intensity
parameter of the process. The A_ (resp., A+) characterizes the rate of the
exponential decay of the right (resp., left) tail of the probability densities,
and c plays the part similar to the variance of the Brownian motion.
( ^ ^ M L M , (i.6)
and Lf is in Co(R), the space of continuous functions vanishing at the
infinity. The map / t-t Lf is called the infinitesimal generator of the
process X. The infinitesimal generator admits an explicit representation in
terms of the generating triplet:
• From the point of view of the theory of PDO, Levy processes are
convenient since the symbols of the infinitesimal generators are in-
dependent of the state variable, x (in the PDO-language, these are
constant symbols), and therefore, many results can be obtained by
using the simplest tools of Complex Analysis.
• To study boundary value problems, explicit formulas for the factors
in the Wiener-Hopf factorization are needed, and Eq. (1.4) allows
one to derive relatively simple formulas. This explains why the
regularity condition in the definition of RLPE is important 7 .
• The derivatives of characteristic exponents of processes of model
classes of RLPE grow at the infinity slower than the characteristic
exponents themselves (cf. Eq. (1.5)); this property allows one to
generalize the class RLPE and construct a general class of Levy-like
Feller processes with good properties (see Chapter 14).
• Moreover, in empirical studies of Financial markets, the charac-
teristic exponents are characterized by a large parameter, like the
parameter a in Eq. (1.5), and so it is possible to develop approxi-
mate effective formulas. We will pursue this possibility in Chapters
7-8 and 10-11, where we study first-touch digital options, barrier
options, investment under uncertainty and endogenous default.
So, the characteristic exponents of RLPE's enjoy almost all desirable prop-
erties but one: the so-called transmission property, which ensures that a
solution of a regular boundary value problem, e.g., the Dirichlet problem,
6
Of course, all these definitions and constructions generalize to the n-dimensional case
7
Eq. (1.4) can be relaxed: see Chapter 15
Pricing of contingent claims 13
has the market value S • 6 = Yl^i Sj@j an^ payoff DT6. We say that a
portfolio 6 is an arbitrage portfolio if S • 6 < 0 and DT6 > 0, or 5 • 0 < 0 and
DT9 > 0. If an arbitrage portfolio exists, there exists an opportunity of
free lunches. In real financial markets, arbitrage opportunities may appear
but they are promptly eliminated due to the activity of arbitrageurs, who
make money by looking for those opportunities. Thus, the assumption of
no free lunches is sufficiently realistic, and as we will see, it can be used as
a cornerstone of the pricing theory.
Introduce an augmented payoff matrix
-S
TZ
DT
The no-arbitrage assumption implies that there does not exist a portfolio
6 such that 1Z6 > 0 and one of the components of 116 is positive. By the
separating hyperplane theorem, there exists a row vector A € R + + 1 such
that
XR. = 0. (1.12)
S = \DT. (1.13)
The A is called the vector of state prices. Notice that Eq. (1.13) deter-
mines prices but makes no use of the probabilities pj of the states of the
market tomorrow. In other words, for the no-arbitrage pricing, only the
information about possible future events but not their probabilities mat-
ter. Another interpretation is: if investors agree on the set of future events
(possible values of securities prices tomorrow) and there are no arbitrage
opportunities, they may disagree on probabilities of those events.
Suppose that our two-state model describes an investor who buys a
portfolio of securities at time 0, and liquidates (that is, sells) it at time
1, with no dividends paid in between. Then Djk admits a natural inter-
pretation as the price of the security Sj in state fc, at time 1. Denote
Sj(0) = Sj, and define a random variable Sj(l, •) on the probability space
fi = {wk | k — 1 , . . . ,m} by Sj(l,Uk) — Dkj- Assume further that one of
the securities, say, S„, is the riskless bond (usually denoted by B) yielding
Pricing of contingent claims 15
a riskless return r; thus, B(l, CJ) — (1 + r)B(0). In this case, the last of the
equations in the system Eq. (1.13) is
m
J3(0) = ^ A f c ( l + r ) B ( 0 )
k=i
(notice that S in Eq. (1.13) is a row vector, so the equations are written
in a row), and therefore, vector q = (1 + r)X satisfies the following two
conditions:
m
0<qk<l, allfc; ^ 9 f c = l. (1.14)
fc=i
S*(t) = (l + r)-tS(t).
The same sort of argument as in the two-period model above shows that
if there is no arbitrage opportunity, then there exists a new probability
measure Q on the same measure space (U,TT) such that for all 0 < s <
t <T,
S*{s)=E<*[S*(t)\Ts}, (1.16)
1.4 are binomial model and trinomial model, respectively. In the binomial
model, at each time step, S(t+1)/S(t) assumes the value u with probability
p G (0,1), and the value d with probability 1 — p; in the trinomial model,
S(t + l)/S{t) can assume three values, with probabilities pj £ (0,1), j =
1,2,3.
In the binomial model, the EMM is unique, and in the trinomial, there
are infinitely many EMM. Clearly, the latter is more flexible in the sense
that it is much easier to adjust its parameters to data.
Example 1.7 In the two securities-two states model, suppose that the
payoff on the risky security in the first state of the world tomorrow is greater
than the one in the second state, that is, du > du- Introduce a call option
on the risky security, with the strike price K, du < K < du. An option
owner will exercise the option and buy 1-security tomorrow if the first of
the possible states of the world materializes, and will not otherwise. Hence,
the payoff row F ( l ) is [du — K 0], and the price F(0) of the call option
today is
dn d\2 ^13
D dn - K d\2 - K 0
l +r l+r l+r
The reader can easily verify that rankD = 3 if and only if d\\ ^ dw Thus,
in the case d n = d\2 the option is redundant, and in the case d n > di2 its
introduction makes the market complete.
Had the real financial markets been complete, there would have been no
need in the creation of derivative securities. This observation implies that
models of incomplete financial markets are more realistic than models of
complete markets.
In multi-period models, the situation is similar, only the role of a portfo-
lio is played by a trading strategy (another name: dynamic portfolio). For
details, see references in the review of literature at the end of the chapter.
Pricing of contingent claims 19
8
That is, the measure inferred from the observations of returns
20 Introduction
r + ijjQ(-i)=0. (1.18)
We call Eq. (1.18) the EMM-condition. There are more subtle restrictions;
in particular, parameters c and v in Eq. (1.4) must be the same for the
historic measure and an EMM. Still, for any model class of RLPE, free
parameters remain, and one can introduce additional degrees of freedom
by considering mixtures of models processes. In Section 4, we produce
numerical results to show how one can change the price of an option and
the shape of the smile by playing with parameters of KoBoL.
Notice that if X is assumed to be Gaussian both under the historic
measure and an EMM, then the condition: c in Eq. (1.4) is fixed means
that the variance of the process does not change under the change of the
measure, and the EMM-requirement Eq. (1.18) fixes the drift (j, of an EMM:
r - / i - y = 0 ,
that is,
^Q(0 = y e 2 - i ( r - y ) € . (1-19)
The oldest variant of EMM is the Esscher transform, which have been
used in Actuarial Science for several decades, and in Financial Mathemat-
ics, from the beginning of the last decade: in terms of the characteristic
exponent, one looks for V Q m the form
For a chosen EMM Q and an European option with the expiry date T and
the terminal payoff g(Xr), the pricing formula Eq. (1.16) can be written as
Example 1.9 Consider a European call option with the strike price K
and the expiry date T. The terminal payoff is g(X{T)) = (e x ( T ) - K)+,
and the integral in Eq. (1.23) is well-defined for £ in the half-plane 3£ < — 1.
Take any a < - 1 , and calculate, for 9f = a:
r+oo
g(£) = / ( c (-«+i)* _ Ke-*')dx (1.24)
»€ - 1 %
_ Ke~^lnK
{% - m
~ (* + 0£ '
By substituting into Eq. (1-22), we obtain
The integral in Eq. (1.22), Eq. (1.25) being an example, can be calculated
by using the Fast Fourier Transform (FFT) or its modifications; there are
cases when FFT performs poorly, and other methods must be used. In
some cases, the probability densities of probability distributions of a process
can be calculated explicitly, and then the pricing formula can be written
22 Introduction
essentially in the same spirit as in the Gaussian case: only instead of the
tabulated normal distribution another distribution is used.
When fitting the model to real data, the observed prices of the option are
compared with the ones calculated from Eq. (1.22). There are many papers
devoted to the derivation of the formula for EMM from the formula for the
historic measure by using different heuristic and economic arguments (see
a discussion in Chapter 4). We believe that the parameters of EMM should
be inferred from data on both the stock and options, and not derived by
some formal reasoning.
The perfect hedging in the Levy market is impossible, and an investor
can only try to minimize the risk, after some measure of the risk is chosen.
and expanding into the power series in Eq. (1.29) rather than in Eq. (1.26).
If we add to Eq. (1.28) the terminal condition
or
2 / 2\
dxf
dtf + \dlf +[r-Y) ~rf = °- C 1 - 31 )
This is the Black-Scholes equation. In terms of the variable S = ex,
Eq. (1.31) assumes the standard form
2
ft/+yS2a|/ + r S 0 s / - r / = O. (1.32)
Apply Eq. (1.38) with r = T(B)AT, and compare with Eq. (1.37). In view
of Eq. (1.34) and Eq. (1.35), we conclude that Eq. (1.33) must hold (this
requires justification: see Chapter 2). Thus, the naive way of writing the
boundary value problem for the price of an option is correct.
In many cases, for RLPE in particular, the justification is simplified
significantly by using the potential theory of Levy processes.
with the density g°(X(t),t). The typical examples are bonds of corpora-
tions: a typical bond pays some interest during its life-time; when the bond
matures, it pays the principal, and in the (random) event of the default,
the bond expires worthless or some amount is paid. In this case, the price
of the security is the sum of the stochastic integral which expressed the
present value of the stream of revenues, and the expectation of the payoff
ge(X(r),T) on the expiration date:
+ £[e-rV(*(r),r) \X{0)=x].
The problem Eq. (1.44)-Eq. (1.45) is called the Wiener-Hopf equation, and
it can be solved by the Wiener-Hopf factorization method. The following
scheme can be realized under fairly weak regularity assumptions on the
data and the symbol a(£) := r + V' Q (0; m the case of an RLPE, the latter
is sufficiently regular.
Step 1. Factorize a(£), that is, represent it in the form
a(0=a+(0M0. ^ R , (1-46)
Analytical methods used in the book 29
where a+ (resp., a_) admits the analytic continuation into the half-plane
9£ > 0 (resp., 9£ < 0), and does not vanish there. If a± and l / a ±
are polynomially bounded in the corresponding half-plane, the factors are
uniquely denned up to scalar multiples, due to the Liouville theorem, so in
some cases, they can be guessed.
The factorization problem can be solved for fairly wide classes of sym-
bols, by using analytical methods, and these constructions admit general-
izations for the multidimensional case, when one considers the Wiener-Hopf
equation in the half-space x„ > 0, and the analytic continuation is made
with respect to the dual variable £ n .
When ip® is the characteristic exponent of a one-dimensional Levy pro-
cess, the Wiener-Hopf factorization has an important probabilistic interpre-
tation in terms of the supremum and infimum processes Mt = s u p 0 < s < t Xs
and Nt =mf0<s<tXs:
r(r + ^ Q ( 0 ) _ 1 = ^ ( 0 ^ ( 0 , (1-47)
where
Certainly, Eq. (1.48)-Eq. (1.49) are not explicit, and some additional efforts
are to be made to produce analytical formulas.
Step 2. The solution to the problem Eq. (1.44)-Eq. (1.45) in a natural
function class exists (for instance, this is the case if the data are smooth and
of compact support, and the solution is sought in the class of continuous
bounded functions); the solution is unique and given by
u = 0r-(D)l(h>+oo)0+(D)r-1/G, (1.50)
(2) The third step requires the knowledge of some technique of the
reader, and in the end the final answer has to be computed numer-
ically by using some integration procedure.
(3) Fortunately, for wide classes of processes and parameters values
observed in real financial markets, effective approximate formulas
are available, which we demonstrate in the main part of the book.
(4) The very form of Eq. (1.50) allows us to derive an algebraic equation
for the optimal exercise price of the American perpetual options for
fairly general class of payoffs, and for the investment threshold in
Real Option theory; this equation is new even in the Gaussian case.
(5) The same form allows us to formulate the correct form of the Mar-
shallian Law and naturally separate the two factors which influ-
ence the capital accumulation. These results are important for
Economics, and are new in the Gaussian case as well.
(6) Finally, if the reader is willing to work with non-Gaussian Levy
processes, there is no hope to use differential equations as in the
Black-Scholes theory anyway: Steps 1-3 are the simplest general
scheme available.
In the case of American options with the finite time horizon, the explicit
analytical methods are not available, but after the discretization of time,
the Wiener-Hopf method can be applied as well.
idea of locally risk minimizing hedging. We show that for RLPE of order
less than 2, the hedging ratio is Holder continuous until expiry, even at
the strike, whereas the Gaussian delta-hedge is discontinuous at the strike,
at the expiry. This makes non-Gaussian RLPE-hedging much more stable
than the Gaussian one. We produce numerical examples and compare the
hedging ratio for different processes.
In Chapter 5, we consider perpetual American options. We formulate
the optimal stopping problem, and the corresponding free boundary prob-
lem. We explicitly solve the problem for a fairly general payoff g, and show
that in the case of the put options and similar more general options, the
optimal exercise price H = eh is determined from the equation
(^(D)-19)(h) = 0,
where <f)~ is determined from Eq. (1.49) with q = r + A, where A > 0 is the
dividend rate. In particular, for the put, the equation reduces to
/•+oo
[ / e->t+Nidt
(<P+(D)-1g)(h)=0.
By using these explicit formulas, we show that in some cases, the smooth
pasting principle fails, and discuss its generalizations.
In Chapter 6, we consider the American options with the finite time
horizon. The explicit formulas being non-available even in the Gaussian
model, we consider RLPE-analogues of several approximate methods used
in the Gaussian case. It is shown that the behaviour of the RLPE- price
of the American put near expiry drastically differs from the one in the
Gaussian case.
In Chapters 7 and 8, we derive explicit pricing formulas for touch-and-
out options and barrier options, respectively. We consider the asymptotics
An overview of the results covered in the book 33
of the price near the barrier, and explain why this is the place, where the
Gaussian model differs most from non-Gaussian ones.
In Chapter 9, we consider simplest of the multi-asset contracts, and
multi-asset hedging. We show that in the case of highly correlated assets
in the portfolio, the difference between the Gaussian (seemingly) riskless
portfolio and the RLPE-locally risk-minimizing one can be quite substan-
tial. Notice that the main cause for the Long Term Capital Management
disaster was the mispricing of risks, and the hedging model we suggest is
not much more difficult to implement than the Gaussian hedging.
1.6.6 Extensions
In Chapter 13, we consider a discrete-time model, corresponding to perpet-
ual American options, that is, Bermudan options. From the analytic point
of view, there is not much difference with the corresponding continuous
time models, but the model becomes (infinitely) more flexible.
In Chapter 14, we consider the simplest extension to the case of Feller
processes: NIG-like Feller processes constructed by Ole E. Barndorff-Nielsen
and Levendorskii (2001). This is the first step for developing non-Gaussian
analogues of interest rate models.
36 Introduction
1.7 Commentary
Levy processes
In the first three sections of this chapter, we basically follow Bertoin (1996)
and Sato (1999) ([B] and [S]). The account is by no means complete, and the
chapter is not meant to substitute for any standard text on Levy processes.
We only aim to provide the list of basic definitions and results, which are
necessary should the reader want to go into the details of our constructions.
In Section 2.4, we consider the stochastic integral of a Levy process
stopped at the exterior of an open set, and deduce a boundary value problem
for the generalized Black-Scholes equation, which the integral solves. This
reduction is the cornerstone of the modified MBS-theory, and the main tool
in Chapters 4-11 and 14. The reduction is based on the Dynkin formula
and the representation of the q-order harmonic measure of a set relative to
a point in terms of the q-potential measure. It admits a generalization to
the case of a strong Markov process with absolutely continuous q-potential
measure.
39
40 Levy processes
»[UAj} = J2rtAn).
The triple (fi, .F, /i) is called a roeasune space. A measure space (il, T, fi)
and the cr-algebra J- are said to be p- complete if .F contains any B C fi
of outer measure 0: inf/x[{A G F\ A D B}] — 0. If .7* is not complete,
the completion of F" is the smallest u-algebra containing all subsets of fi of
outer measure 0.
The measure fj, is called finite if fi(Q) < 00, and er-finite if there is a
sequence An of elements of T such that fi(An) < 00, and L)An = fl.
Let (fi, J7, P ) and (Q',.F',P') be measure spaces. We say that / : Q, —>
ST is a {IF/!F')-measurable function or measurable function for short if for
a n y y l ' G - F , . T 1 ^ ' ) <E-T7-
E[g(X)) = [ g(x)Px(dx)
P[Ajln-.-nAjn\ = P[Ah]...-p[Ajn}.
MO = A(-0-
The function
1
The inverse Fourier transform T is denned by
The dual measure to (i, denoted ft,, is defined by jl(B) = fi(—B), where
-B = {x\ -x£B}.
The next theorem lists several properties of the characteristic functions,
or equivalently, the Fourier transform of measures. For the proofs and
additional properties, see Section 2 in [S].
Px{t) = PxAti)---Pxm(U),
for f = ( f i , . . . , f m ) e R n , w/iere n = m H h nm.
f#) / / // ftas a ./irate absolute moment of order m € N , i/iai zs,
can be found as
limP[|X.-Xt|>e]=0.
s—>t
A Levy process on R™ is called an n-dimensional Levy process. If the first
condition is not required, X is called a Levy process in law.
Denote by /z m the m-fold convolution of a probability measure with
itself:
fim = n * • • • * fi (m times).
/i* = F~l exp[£0] is well-defined and infinitely divisible, and there is a Levy
process in law, X, such that Pxt = A**-
(ii) If {-Xt}t>o is a Levy process in law on R n , then for any t>0, Pxt
is infinitely divisible and Pxt = /•**, where /x = P x i •
(Hi) If X and Y are Levy processes in law on R n such that Pxt = P ^ ,
then X and Y are identical in law.
(iv) Let X be a Levy process in law on R n . Then it has a modification
which is a Levy process.
Theorem 2.5 (i) Let X be a Levy process on R". Then its characteristic
exponent admits the representation
The triple {A, II, 7) is called the generating triplet of X. The A and II
are called the Gaussian covariance matrix and Levy measure of X. When
II = 0, X is Gaussian, and if A = 0, X is called purely non-Gaussian.
Notice that, essentially, the term —i(x,£)l£>(x) in Eq. (2.15) is needed
in order to ensure the convergence of the integral, and hence other functions
can be (and are) used instead of c(x) := 1 D ( X ) , for instance, c(x) = 1/(1 +
|x| 2 ); the A and II are independent of the choice of c. If II satisfies the
condition
which is stronger than Eq. (2.16), then Eq. (2.15) can be simplified
where
7o=7- / xlD(x)Yl(dx).
For the proof and formulas for the generating triplet of V> in terms of the
drift and Levy measure of Z and the generating triplet of Y, see [S], The-
orem 30.1.
Au = UAU', Uv = IK/" 1 ,
Theorem 2.8 (i) Let q > 0. There exists a unique pair of infinitely
divisible distributions p+ and p~ having drift 0 supported on (—oo, 0] and
[0, +oo), respectively, such that their Fourier transforms <j>+ and <j>~ satisfy
(ii) The functions <f>+ and (j>~ admit the following representations
r+oo r+oo
<t>t(£.)=q e-qtE[e^Mt}dt =q e-gtE[ei(<-Xt-Nt)}dt, (2.24)
Jo Jo
f+oo r+oo
<t>q{€)=q\ e-^E[e^N*\dt = qj e-^E\e^Xt-Mt)\dt, (2 25)
Jo Jo
and
r+oo r+oo
<£+(£) = exp / t-xe~qtdt I (eixii - l)ii*{dx) (2.26)
Jo Jo
r+OO rO
Example 2.3 Let X be a Brownian motion with the drift 7 and variance
a2. Then the characteristic exponent is
m -r - i-rt.
It is clear that for q > 0, the equation q +1/>(£) = 0 has two roots —i/?_
and -i/3+ in the upper and lower half-planes, respectively, and therefore,
52 LSvy processes
X
q(q + V>(£)) admits the factorization Eq. (2.23) with
Since 4>^(Q has no zeroes in the lower half-plane 3?£ < 0, Eq. (2.23) implies
that — i/3 is the root of q + ip(£). After this root is found, which is a simple
numerical procedure, we can calculate 0+ from Eq. (2.29), and after that
4>~ from Eq. (2.23).
The reader can consider spectrally positive Levy processes by using the
symmetry argument.
The simple Markov property can be reinforced by allowing for certain ran-
dom times, that is, random variables with values in [0, +oo] instead of the
deterministic time t.
54 Levy processes
A random time T is called a stopping time if for every t > 0, the event
{T < t} belongs to the <r-algebra Tt- Since the filtration is right-continuous,
we have that T is a stopping time if and only if {T < t) g Tt for alH > 0.
For an arbitrary random time T, we denote by TT the P-completion of the
tr-algebra generated by the process X killed at time T, X o k r , 1 and the
variable 1{T<+<X>}XT-
It can be checked that when T is a stopping time, !FT consists of events
A such that {T < t} n A G Tt for all t > 0. (Intuitively, a random variable
T is a stopping time if one can decide whether T < t or not by making
observations of the evolution of X up to the moment t; A 6 TT if A is an
event expressible by observations up to time T).
Note that if B is closed, the last equality may fail for some x on the bound-
ary of B.
J
kt denotes the killing operator, ktw(s) = w{s) if s < t, and d otherwise, where d denotes
the cemetery point of the extended probability space fi U d (for details, see [B], Section
0)
Levy processes as Markov processes 55
Uqf(x) = f f(y)U«(x,dy)
Uq - UT + (q - r)UqUr = 0, Vq,r>0,
and the resolvent equation implies that the image of Co := Co(R-") under Uq
does not depend on q. We denote it by V. The Feller property implies that
for every function / € Co, limq-^+oo qUqf = / in the uniform topology.
Thus, V is dense in Co- By using the last observation and the resolvent
equation once again, one deduces that the map Uq : Co —> V is a bijection.
We define the infinitesimal generator L : T> —* Co as the inverse to Uq :
Co -+ V. Thus,
and
(q - L)Uqf = f, V / € V. (2.32)
56 Levy processes
F(Ptf)(t) = E
VR"
= E f emxt-v).e>f(y)dy
.JRn
= E J{Xt,S)
JR"
and therefore,
(2.34)
Uy = {q + ^{D))-lf. (2.35)
(9 + S t y ) - 7 e L i . (2.36)
Lemma 2.2 If f G L\ and Eq. (2.36) holds, then Eq. (2.35) is valid.
(U"f)(x) = / e-qt(Ptf)(x)dt
Jo
/•+oo p
= / e-«*(27r)- n / e- < < I ' € >-iK«/(0d£<ft.
JO JR»
Due to Eq. (2.36), the last integral computed in the reverse order dtd£
converges absolutely, and hence we can apply the Fubini theorem and obtain
Eq. (2.35). •
and then
+ f (f(x + y)-f(x)-lD(y){y,f'(x)))U(dy).
2.3.5 Duality
The process X — —X is called the dual process of X; it is a Levy process
for the same triple (CI, T, P ) . The sign~is used to denote the ^-resolvent Uq
and the infinitesimal generator L of X, etc. For x e R n , P 1 denotes the
law of x + X under P , that is, the law X under P _ x . It is clear that ifi is
the complex adjoint to tp: $(£) = VK-£)>£ € R n , and hence, L = —ip(-D).
The name "dual" stems from the following proposition.
/ Ptf(x)g(x)dx = / f(x)Ptg(x)dx,
JR" JRn
and for every q > 0,
/ Uqf(x)g(x)dx = / f(x)Uqg(x)dx.
JR" JR"
/ Lf(x)g{x)dx = I f{x)Lg{x)dx.
JR" ./R"
Proof. The first statements is obtained by a change of variables, see [B],
Proposition II. 1.1 and [S], Proposition 41.7, and the last two follow. •
(i) For some q > 0, the measure V is absolutely continuous w.r.t. the
Lebesgue measure;
(ii) For any q > 0, the measure Vq is absolutely continuous w.r.t. the
Lebesgue measure;
(Hi) For some q > 0, and some x € R n , the measure Ug(x,dx) is
absolutely continuous w.r.t. the Lebesgue measure;
(iv) For every q > 0, and every x 6 R n , the measure Uq{x,dx) is
absolutely continuous w.r.t. the Lebesgue measure;
Levy processes as Markov processes 59
(v) The resolvent operators have the strong Feller property, that is, for
every q > 0 and f G Loo(R n ), the function Uqf is continuous.
When the assertions of Theorem 2.11 hold, we say that the resolvent kernel
is absolutely continuous, or that the (ACP)-condition holds. The following
Lemma gives convenient sufficient conditions which are satisfied by Regular
Levy Processes of Exponential type.
Lemma 2.4 Let an n-dimensional Levy process X satisfy (ACT)-condi-
tion, that is, the transition measures Pxt are absolutely continuous for all
t > 0. Then (ACP)-condition holds for X and for (n + 1)-dimensional
process {(Xt,t)}.
where vq is the density of Vq, and vq(x) = vq(—x) is the density of Vq.
2.3.7 Operators P%
n
Let B C R be an F a -set (for the general definition, see [S], p.279; for our
purposes, it suffices to notice that Borel sets are Fo-sets), and let T'B be
the hitting time of B by X. For g € Loo(B), define
Plg{x):=Ex[e-qT'Bg{XTlB)]. (2.43)
The map g i-> Pgg(x) defines the measure which is called the g-order
harmonic measure of B relative to x. By gathering Proposition 42.13,
Theorem 42.5 and Definition 42.6 from [S], we obtain
Lemma 2.5 Let B be an Fa-set. There exists a a-finite measure dfis
supported on B such that
PqBl(x)= f vq(y-x)d»B(y).
60 Levy processes
Proof. We derive Eq. (2.45) from Eq. (2.44) by using an additional as-
sumption which holds for RLPE's: the characteristic exponent ift of X
admits the analytic continuation into a tube domain R n 4- iU, where U is
an open subset of R n , containing 0.
For 7 G R n , define a function u 7 by u7(x) = e^1,x^ . As the first step of
the proof, we show that Eq. (2.45) holds with g = u 7 provided q+ip^i^) >
0. Change the probability measure: P 7 = e ^ ^ H ' ^ - n O p . Under the new
measure, X is the Levy process with the characteristic exponent ip7(£) —
VK£ ~~ *7) — iJ{~i'y) since
therefore
,-{q+H-i-l))T'B 9+V>(-t7) 1(2:).
e-b'*) ptuyix) = E% =:Pi 7,B
or equivalently,
and finally,
1
( P | u 7 , (q - u ' ^ - V H 7 ) K K ) = 0, V «,' e C0°°(SC).
= -V>(#)
= i,
and hence -u~ 1 (L 7 — T/>(—i7))u7 = L. Eq. (2.45) with g — u7 has been
proved.
For the next step, consider a compact B and an arbitrary bounded mea-
surable g. We can approximate g in the L°°-norm by continuous functions,
and each continuous function by polynomials (the Stone-Weierstrass theo-
rem). Fix e > 0 such that the ball V(e) := {7 | ||7|| < e} C U. Since for
any multi-index a,
iimon7-^(e^-ir-n^,
3=1 3=1
gN(x)= Yl cNi~/e{'y'x),
-7€iv
and
Ex
'e-iT'Bg(X(TB))\ <EX \e-^*™g(X(rB{R)))\+\\g\\o0E' g 1TB\B(R)
For a Levy process, T'(B \ B{R)) -+ +00 as R -> +00, a.s., hence
-iT'g
Ex 0 as 72 —> +00,
/ PqB(R)9(x){q-L)w{x)dx = Q, VweC™(Bc).
This section is the focal point, where Probability and Analysis meet to
produce in the end analytical formulas for prices of contingent claims and
value functions of firms in Real Options theory. By definition, each of those
prices and value functions is a sum of the stochastic integral which represent
the expected discounted stream of payments (or revenues) during the life-
time of the contingent claim or the firm, and the expected terminal payoff in
the end. By using Dynkin's formula and the potential theory, we deduce the
pseudo-differential equation for the price in the region, where the contingent
claim remains alive, and the terminal payoff translates into the boundary
condition. Unlike in the Gaussian case, the boundary condition is non-local:
the unknown function must be specified not only at the boundary but the
other side of the boundary as well. On solving the boundary problem in an
appropriate function class and obtaining a formula for the unique solution
in this class, we obtain the price we have been looking for. In the last step,
we use standard technical tools from elements of the theory of generalized
functions, the theory of PDO and Complex Analysis.
Boundary value problems for the Black-Scholes-type equation 63
U,{l-L)w)L2={gc,w)L„ (2.49)
-E
-B^e-^'BE*™
r
Jo
e-^B[/V(Xri3)l+P|50(x)
" 9 * ncc
9 (Xt)dt + PqB9°{x)
P%h{x)
64 Levy processes
We will see that the prices of contingent claims with finite time horizon,
e.g., barrier options, can be expressed in this form.
Theorem 2.13 Let X be an n-dimensional Levy process satisfying (ACP)-
condition. Then the stochastic integral Eq. (2.50) is a bounded solution to
the following boundary value problem
{ gr(x,t),
0,
t < T a n d (x,t) eD
t > T and x £ R n .
We can write Eq. (2.50) as Eq. (2.46), with (x,t) and X in place of a; and
X:
f(x,t) = Ex>t
Jo
f e-qsgc(Xs)ds + Ex>t[e-qTBg0(XT')} (2.55)
Remark 2.2 More general payoffs, and processes which do not satisfy
(ACP)-condition, can also be considered.
Remark 2.3 Since the proof of Theorem 2.12 is based on the results from
the theory of Markov processes having absolutely continuous potential mea-
sures, the main results of this section admit the straightforward generaliza-
tions for Markov processes.
2.5 Commentary
For the systematic treatment of the general theory of Levy processes, see
[B], [S], and for a shorter exposition, Rogers and Williams (1994); for the
general theory of Markov processes, see Blumenthal and Getoor (1968),
Rogers and Williams (1994), Sharpe (1989) and the bibliography there.
Chapter 3
3.1.1 An overview
D.B. Madan with corroborators were the first to use Levy processes differ-
ent from the stable ones, namely, Variance Gamma Processes (VGP)—see
Madan and Seneta (1990), Madan, Carr and Chung (1998)—in financial
modelling. A host of processes has been introduced, starting from the
family of Generalized Hyperbolic distributions constructed by Barndorff-
Nielsen (1977): Hyperbolic Processes (HP) and Normal Inverse Gaussian
(NIG) processes were constructed in Eberlein and Keller (1995) and Barndorff-
Nielsen (1998), respectively, and their generalizations: Generalized Hyper-
bolic Processes (GHP), and Normal Tempered Stable Levy Processes (NTS
Levy) - in Eberlein and Prause (2001), Prause (1998), Raible (2000), and
Barndorff-Nielsen and Levendorskii (2001), Bamdorff-Nielsen and Shep-
hard (2001c), respectively. By now the list of papers, where processes of
this cluster have been used for modelling behaviour of asset prices and
term structure of interest rates, pricing of derivative securities, and risk
management, is fairly long: Barndorff-Nielsen and Jiang (1998), Eberlein,
Keller and Prause (1998), Prause(1998), Eberlein and Raible (1999), Eber-
67
68 Regular Levy Processes of Exponential type in ID
lein et al. (2000a, 2000b), Eberlein (2001a, 2000b), Eberlein and Prause
(2001), Barndorff-Nielsen and Prause (2001), Barndorff-Nielsen and Shep-
hard (2001c). In Barndorff-Nielsen and Shephard (2001c), Normal Modified
Stable laws were constructed, thereby generalizing the class of Generalized
Hyperbolic distributions. For extensions to Feller processes and stochastic
volatility models, see Chapter 14 and Commentary to Chapter 1, respec-
tively.
(Non-infinitely divisible) Truncated Levy distributions were constructed and
used for financial modelling purposes by Mantegna and Stanley (1994,
1997); see also the monograph Mantegna and Stanley (2000). Novikov
(1994) and Koponen (1995) introduced infinitely divisible analogues of
truncated Levy distributions. The extension of Koponen's family was con-
structed in Boyarchenko and Levendorskii (1999a, 1999b, 2000); we will call
it KoBoL family. Later this family was used by Carr et al. (2001) under
the name CGMY-model. It was shown that processes of the KoBoL family
provide good fit for many equities.
Below, we provide technical details of the derivation of necessary formulas
mainly for KoBoL processes, since for these processes, all proofs can be
made in full by using very simple analytical tools, and Variance Gamma
Processes can be regarded as members of KoBoL family. The reader can
find technical details concerning processes of the other families in op. cit.
Clearly,
+0O
min{l,x 2 }n ± (i/,A;dx) < +00,
/ -00
Lemma 3.1 For u < 2 and A > 0, there exists C{v,X) such that
where
a) if v < 2, v ^ 0 , 1 , then
0(i/,A;O=r(-i/)[A"-(A-iO"]; (3-3)
b) if v = 0, then
0(O,A;O = l n ( A - i £ ) - l n A ; (3.4)
c) if v — 1, then
</>(l,A;£) = A l n A - ( A - * £ ) l n ( A - i O . (3.5)
(3.6)
(ii) if v — 0, then
(Hi) if v — 1, then
iMO = - i / ^ + c+[(-A_)ln(-A_)-(-A_-*01n(-A_-i0]
+c_ [A+ In A+ - (A+ + ifl ln(A + + if J], (3.8)
if v = 0, then
if v — 1, then
mi : = mi.-r/T =
#'(°); (3-12)
Tn2 : = m 2 , r / r = V"(0); (3.13)
(3)
m3 : = lim m 3 r / r = -n/> (0); (3.14)
T—*0 '
The proof will be given at the end of the subsection. Of course, Lemma 3.3 is
valid not only for KoBoL processes, but we will use it mainly for them. For
a KoBoL process of order v ^ 1, by using the property T(l +a;) = xF(x) of
the Gamma function, we deduce from Eq. (3.9)-Eq. (3.10) and Eq. (3.12)-
Eq. (3.15)
' X V—2.0
v-1.5
x
• v-1.0 .
." —- v-O.5
x
•
' ". x X
_"
• •
X
* _*
". * .*
'. X
-'
s -«-
X x
-* •_
X
\ •- " x
x ."
.•
\ . x X .* /
^\^.x
jp'
-300 -200 —100 10O 200 300
Fig. 3.1 The characteristic exponent and log-density of a KoBoL process: dependence
on the order of the process. Parameters: n = 0,X+ = —A_ = 24,rri2 = 0.18, r = 1 day,
v £ {0.5; 1.0; 1.5; 2}
1 1 1 1 i
Fig. 3.2 Probability density of a KoBoL process: dependence on time and the order of
the process. Parameters: n — 0,A+ = —A_ = 24,nx% = 0.3, v G {0.5; 1.0; 1.5;2}. Left
panel: j — \ day, right panel: T = 5 days.
Fig. 3.4 Skewness and kurtosis of the probability density of a KoBoL process: depen-
dence on the steepness parameters. Parameters: fi — 0, v = 1.4, mi = 0.18, T = 1
day.
+oo .
/
(\e~Xx - (A - i O c - ^ " ^ * ) x~vdx
r(1
" V - (A - »0"] + i£C(i/, A),
where
Since -i/r(-i>) = T(l - v), Eq. (3.2)-Eq. (3.3) have been proved.
Finally, if v E (1,2), we integrate by part twice, and arrive at the same
formulas Eq. (3.2)-Eq. (3.3), with
dx
__1 l-e~A
~ /x + A A
76 Regular Levy Processes of Exponential type in ID
= lnV +
^~A—
c) If v = 1, we differentiate ^ + (1,A,£) with respect to A and get
—i(>+(0,\,£). Clearly, liniA-n-oo V^Clj ^>£) — 0, therefore we can find
V>+(1,A,£) by integration:
+oo
V»+(0,/i,OdM
/A
/"+oo
fo
° /_ e~^ -_ 1 \
= / Iln(/x-i0-ln^-^ J d/z
= (0* - *0 1Q(M - *£) - / i + i ^ - Z i l n / i + zi + ^ l n / i ) ! ^ 0 0
—z£ / d/x
J\ V
= AlnA-(A-i£)ln(A-i£)-i£ l+lnA+/ dfx
(here we have used lim^_,. +00 [(^ — i£)\n((i — i£) — filnfi — ifc(l+ln/i)] = 0),
and Eq. (3.2) and Eq. (3.5) follow.
Proof of Lemma 3.3. Since
+oo
e-faS-nKO^ (3 20)
/
-oo
we have
+oo
xpT(x)dx
/ -oo
p + OO /" + 00
-t-oo /-t-oo
/ xe- i x «- r , ^>d£d:r
/ -oo J —oo
/+oo /.-foo
/ e-^^de-^dx
-oo J—oo
-f-OO /- + 0 0
/
/ e-"^'(0e-r*(€)dfda:.
-oo J—oo
Model Classes 77
-t-OO /-+O0
/ i2e-fa«-T*«)deda;-m?iT>
/ -00 J — OO
and integrating by part twice and using Eq. (3.21) and Eq. (3.12), we obtain
Eq. (3.13). Eq. (3.14) and Eq. (3.15) are derived similarly.
The distribution with density Eq. (3.22) {v G (0,2), 8 > 0,7 > 0) is referred
to as a tempered stable law, and we denote it by TS(v,8,j). TS(v,8,y) is
infinitely divisible, and its Levy density and characteristic exponent are
given by
V2V/2 X
u(x) - A ~ -l-u/2 -\i2x
78 Regular Levy Processes of Exponential type in ID
Let Y(B) be the Brownian motion with drift B, and let X be the subor-
dination of Y(B) by TS(u, 5,7). Then the characteristic exponent of X,
denote it (f>, is of the form
.2
and by using Eq. (3.23), we obtain
a
Pt(X) = -eMAS^^f2 - M + 8X]K^.6{jf-^6)\ (3.26)
$(0) = c[ln(A + 2 0 ) - l n A ] ,
where Ky denotes the modified Bessel function of the third kind with index
A. An integral representation of Ky is given by
1 f°°
K^{z) = ^ J V1'1 exp[-0.52(y + y^jdy. (3.29)
The a,/3,n and 5 play the same parts as in the case of normal inverse
gaussian distributions, which obtain as a subclass of the class of generalized
80 Regular Levy Processes of Exponential type in ID
d{x) =
o Jfx rr—^ exp[-oy<52 + ( x _ M ) 2 + p{x _ rfi (3.30)
2aoKi{dy/a2 - /32)
_ ,w I »' - ^ "i V2 * » ( V « 2 - ( / » + «) 2 )
(3 3I)
«o —UrcmtFJ w-w '
for |/3-Hf| < a. A generalized hyperbolic distribution is infinitely divisible,
and the density of the Levy measure is: if A > 0, then
0x f /-co e x p ( - \ / 2 y + a2|x|) \
stabilize to positive constants as y -> +0, and have the following asymp-
totics, as y -> +00:
+0(!rl) = (1+0to 2) ( 5)
K^ ) ^ ~" ' "
Hence, by making the change of variable y = \x\~2z2/2 in Eq. (3.32) and
Eq. (3.33) and using Eq. (3.35), we obtain Eq. (3.34). D
From Eq. (3.34), we conclude that all Generalized Hyperbolic Processes are
of infinite variation; in particular, NIG and HP are.
and
For K0B0L, NTS Levy processes and VGP, Eq. (3.36)-Eq. (3.37) easily
follow from the formulas for the characteristic exponents; for GHP, the
verification of Eq. (3.36)-Eq. (3.37) requires more effort. It follows from
Eq. (3.36)-Eq. (3.37) that if two model processes are characterized by dif-
ferent A_(l) < A_ (2), then the right tails of probability densities of the
second process are fatter, independently of the classes the processes belong
to, and of the orders of the processes.
If A_(l) = A_(2), then the behaviour may differ: for instance, Eq. (3.36)
holds with u = A_ for K0B0L processes of order v > 0 and NTS Levy
processes of any order v e (0,2), whereas Eq. (3.37) holds with the same
UJ = A_ for VGP and HP. Thus, the right tails of probability densities of the
latter group of processes are fatter than those of the former group, but this
observation has hardly any practical significance: it is very difficult to infer
A_ from any empirical data set, and hence an error in A_ will dominate an
additional subtle difference in the tail behaviour.
82 Regular Levy Processes of Exponential type in ID
It can be shown that of two KoBoL processes with the same A_, the
one of the lower order has the fatter right tail, and the same is true of NTS
Levy processes.
In the same manner, A+ is characterized by the behaviour of the left
tail.
/
e- A + x n(da;) + / e-x-xU{dx) < oo. (3.38)
-oo Jl
Lemma 3.5 Let X be a Levy process of exponential type [A_, A+]. Then
(i) the characteristic exponent ip is holomorphic in the strip 9£ €
(A_, A+), and continuous up to the boundary of the strip;
(ii) there exist C and v > 0 such that for all £ in the strip 9£ G [A_, A+],
mO\<C(l + W; (3.39)
(Hi) for any q > 0, there exist S > 0 and er_ < 0 < cr+ such that for any
[w_,o;+] C ((T-,<T+) and all £ in the strip 9£ € [CJ_,CJ + ],
has at most one purely imaginary root in the lower (resp., upper)
half-plane, call it —i(i+ (resp., —ift-);
Two definitions of Regular Levy Processes of Exponential type 83
g + lK»(A±=FO))<0, (3.43)
and if it exists, it is a simple root.
Proof. (i) is immediate from Eq. (3.38), and (ii) can be easily deduced
from the Levy-Khintchine formula, by considering separately the integrals
over |x| < |£|- x and \x\ > | £ | - J .
(iii)-(iv) Set M\{a) — f*™ e~axfi1(dx), where ^(dx) is the probability
distribution of Xi. By differentiating twice, we conclude that Mi is convex,
and clearly, Mi(0) = 1 < eq. Hence, there exist w_ < 0 < u+ and 5 > 0
such that for all o e [w_,w+], M\{a) < eq~s.
Now, for any £ € R , and these a,
+oo I /-+oo
/
-oo I J— oo
therefore Eq. (3.40) holds with cr_ = infw_, a+ — supo;+, and Eq. (3.41)
implies Eq. (3.40).
(v)-(vi) Notice that by the proof of (iii), a <-* q + ip(ia) is concave and
equals to q > 0 at 0. •
For the sake of brevity, from now on we consider processes with Levy
measures (almost) symmetric in a neighbourhood of the origin.
Definition 3.3 Let A_ < 0 < A+ and v € [0,2). A purely non-Gaussian
Levy process is called a Regular Levy Process of Exponential type [A_, A+]
and order v if its Levy measure satisfies Eq. (3.38) and, in a neighbourhood
of zero, admits a representation II(dx) = f(x)dx, where / satisfies the
following condition:
there exist v' < i/, c > 0, and C > 0 such that
Notice that a regular Levy process of exponential type has bounded varia-
tion if and only if v < 1, since this is equivalent to f_ (\x\ Al)Il(dx) < +oo.
All model classes in Section 3.1 are RLPE in the sense of Definition 3.3:
• KoBoL processes with c+ = C-, of order u > 0, v ^ 1, are RLPE
of order v and exponential type [A_, A + ];
84 Regular Levy Processes of Exponential type in ID
lMO = - * ^ + # 0 . (3-45)
where <fi is holomorphic in the strip 5£ £ (A_,A + ), is continuous
up to the boundary of the strip, and admits a representation
One can easily generalize both definitions by using c± > 0 in Eq. (3.44) on
the half-axis ±x > 0, and in Eq. (3.46), as *ft£ -4 ±oo.
From now on, we exclude RLPE of order 0, VGP in particular, from our
consideration, since the logarithmic behaviour of the characteristic expo-
nent as £ —¥ +oo leads to technical complications at many places; in fact,
Properties of the characteristic exponents and probability densities of RLPE 85
only for pricing of options of the European type, VGP are as convenient as
the other RLPE's. Each time one has to consider barriers or early exercise
boundaries, one has to use the Wiener-Hopf factorization, which can be
done much simpler when the characteristic exponent grows polynomially
at the infinity. At the same time, we include Brownian motions: according
to Definition 3.4, they are RLPE of order 2 and any exponential type. No-
tice that a mixture of independent RLPE is an RLPE, the order being the
largest one of the orders of the components, and the exponential type being
the intersection of the exponential types of the components. In particular,
any RLPE with a non-trivial Gaussian component is of order 2.
/ -oo
86 Regular Levy Processes of Exponential type in ID
-t-oo—iu
-co—iu
/
/+oo
i
e -««-'-*K+ -)df. (3.50)
-oo
The last equality obtains when we shift the line of integration. To justify
the shift, we use the Cauchy theorem and Eq. (3.48). Differentiate w.r.t. x
under the integral sign in Eq. (3.50) s = 0 , 1 , . . . times; from Eq. (3.48), we
conclude that the resulting integrand admits the estimate via
Cs(l + | £ | ) s e x p [ - < m
Lemma 3.7 Let X be one of the model processes, of order v > 0 and
exponential type [A_,A + ].
Then
(i) the <f> in Eq. (3-45) admits the analytic continuation into the com-
plex plane with two cuts: (—zco,iA_] and [i\+,+ioo), and outside
any neighbourhood of iA_ and i\+ satisfies the following estimate:
(%%) all the roots in the plane with the cuts are purely imaginary.
I? + T K O I < C ( I + KI)
is of a greater order than the lower bound Eq. (3.49). This leads to ad-
ditional problems with solutions near the boundary, in particular, to the
failure of the smooth pasting principle for optimal stopping problems.
the formula for the infinitesimal operator of the NIG Levy process. Here
is a list of observations which naturally lead to NIG Levy:
2
(i) the generator of the Brownian motion is ^-A, but the tails of a
Gaussian distribution decay faster than an exponential function;
(ii) stable non-Gaussian Levy processes have generators of the form
L = -S\ — A|"/ 2 , where S > 0 and v 6 (0,2), and their symbols are
non-smooth at the origin: — <S|£|", which leads to polynomial decay
of the tails of the density functions;
(iii) the tails of the density functions observed in financial markets, in
turbulence and in many other fields of study usually have exponen-
tial decay; this means that the symbols of generators must be not
only smooth but holomorphic in a strip of the form 9£ G (A_, A + ),
where A_ < 0 < A+ (in the multi-dimensional case, in a tube do-
main 9£ € U, where U € R n is an open set containing 0).
L =-8[{a2 - ^ 2
- a% (3.53)
where a > 0, 6 > 0, a > |/?|, and D — —id is the standard notation
in the theory of PDO; it has the symbol £ (The last terms in Eq. (3.53)
and Eq. (3.54) are needed to ensure that L • 1 = 0, which is necessary for
a process without killing). Eq. (3.54) gives the infinitesimal generator of
a (multi-dimensional) NTS Levy process, and in the simplest case of the
square root, of a (multi-dimensional) NIG.
(3.55)
Then
a) (ptiO admits the analytic continuation into a half-plane Sr£ > w_ and
can be calculated as follows:
1
exp (27U)- j _ dr) ; (3.58)
oo+iuj- lit - V)
90 Regular Levy Processes of Exponential type in ID
b) <j>q (£) admits the analytic continuation into a half-plane 9£ < w+ and
can be calculated as follows:
c) </\j~(£) 1 (resp., <j>q (£) 1) admits the analytic continuation into a wider
half-plane 3£ > A_ (resp., Off < X+) by
_ 1
<(0 = 9 _ 1 (9 + ^ ( 0 ) ^ 7 ( 0 . ^ e (A_,w_]; (3.61)
1
W = 9 _ 1 (g + ^ ( 0 X ( 0 , ^ 6 K,A+). (3.62)
Proof, a) Consider the expression under the exponent sign in Eq. (2.26):
r+oo p—qt r+oo
HO • = -J- (eixi ~ l)^{dx)dt
r
Jo Jo
r+oo „-qt r+oo r+oo
= / / (eixt - 1)(2TT)- 1 / e-^-^^drjdxdt.
JO * JO J-oo
On the strength of Eq. (3.55), we may apply the Cauchy theorem and shift
the line of integration:
r+oo „—qt r+oo f+oo+iui-
f(0 = —— (eix* - 1)(2TT)- 1 / e-i">-t*l'l)dndxdt.
Jo <• Jo J-oo+iui-
Now the inner double integral converges absolutely, hence we can apply the
Fubini theorem and integrate w.r.t. x first:
r+oo „—qt r+oo+iuj-
= —-(27TZ)-1 / e-'^ttrt - Z)-1 - ri-^dndt.
Jo * J-oo+iui-
Integrate by part:
r+oo -qt r+oo+iu- _ r
\nri—^txl)'(n)e-t^dndt
z
Jo J-oo+iu- V
/•+OO
r+oo r+ ™
/-t CO+iuJ-
, _ ,
= (27T2)-1 / / v-£.,.,,
ln-!—^xlj'{T])e^''+^^dr1dt.
Jo J-oo+iui- V
The Wiener-Hopf factorization 91
From Eq. (3.56), the integral above calculated in the reverse order dtdr) con-
verges absolutely. Hence we can apply the Fubini theorem once again and
obtain Eq. (3.57); integrating in Eq. (3.57) by part, we arrive at Eq. (3.58).
b) The dual process X is of exponential type [—A+, — A_], its character-
istic exponent is i/>, and [—CJ+, - w _ ] plays the part of [a;_,a;+] in Lemma
3.5. Write down the Wiener-Hopf factorization for X and apply the com-
plex conjugation; then the "+"-factor for X becomes the "-"-factor for X,
and Eq. (3.57) for X becomes Eq. (3.59) for X.
c) follows from Eq. (2.23) and Lemma 3.5, (i). •
Remark 5.1 If X is an RLPE in the sense of Definition 3.4, then Eq. (3.56)
and Eq. (3.55) hold; hence, Theorem 3.2 holds as well.
Lemma 3.8 Let CJ_ and w+ be as in Theorem 3.2.
Then there exists C > 0 such that in the half-plane ± 3 £ > ± w T , </>*
admits estimates
satisfies
lim 5 ( 0 = 1- (3.65)
In all cases, Eq. (3.65) follows from Eq. (3.45)-Eq. (3.46). In the first three
cases, Eq. (3.65) is immediate, and in the last case, the simplest way is to
check that lnf?(£) -t 0 as £ -> ±oo:
hm lnJB(0
KS = ± — K++ q = - ~ - / c _ + l n
? A±oo ' - 2 ^ 2 ' (C2 + M 2)l/2
+ lim ( - K + - K _ + 1) In |f|
f-+±oo
, . 7ri . u
= ± (vK ++ - « - ) —- T larctan— = 0
' 2 c
by our choice of K + and rc_.
The last factor in Eq. (3.64) assumes values in a half-plane Rz > 0 by
Eq. (3.40), and the same is true of the product of the first three factors,
since the first one is positive, A_(f) and A + (f) assume values in the half-
plane but in different quadrants, and 0 < K± < 1. Hence, for all ( £ R ,
—-K < argB(f) < 7r, and therefore, b — In B is well-defined on R. Fix
w_ < 0 < w+ such that <r_ < u>-, ui+ < a+, where cr± are from Eq. (3.40),
and notice that all the arguments above are valid on any line 3f = a 6
[w_,w + ]. Next, for r > <*;_, T\ € [w_,r) and real f, set
• r + OO+in b/ \
It follows from Eq. (3.45), Eq. (3.46), Eq. (3.64) and Eq. (3.65) that there
exist C,Ci,p > 0 such that for any 77 in a strip ^sr] G [w_,ui+],
Hence, the integrals in Eq. (3.66)-Eq. (3.67) converge, and 6+(£) (resp.,
&_(£)) is well-defined in a half-plane §£ > w_ (resp., Q£ < LJ+).
We set
c^ f/iere exist C, c > 0 and pi > 0 swc/i that in a half-plane Sr£ > w_,
1
MO* - A+(0 ± K + 1 < C ( i + i e i ) ± , , + - p l ; (3.73)
1
MO* - A - ( 0 ± K - I < C ( l + \Z\)±K—Pl; (3.75)
^(O-1=a±(O/a±(0). (3.76)
b) By the residue theorem, we have for T\ € (w_, 9£) and TI £ (9£, w+)
2yr
\J-oo+iTi J-oo+iT2 ) S — »7
+0O+iT2
i /• 6(T/) ,
27T,/_ 0 0 . H r 2 ^ — 77
= 6(0-Mfl.
Hence, exp 6+ (£) exp 6_ (£) = £?(£), and Eq. (3.71) is immediate on a nar-
row strip a;_ < 9f£ < u>+; on a wider strip 9£ € [A_,A+], it holds by
construction: see Eq. (3.69) and Eq. (3.70).
c) By using Eq. (3.68), we obtain
0+
#(0 : exp $+(£)•
(3.79)
Similarly, from Eq. (3.59), we deduce, for £ in the lower half-plane,
-P-
4>7(0 = exp*,(0, (3.80)
-P- + tf
where
r+oo i\>'{iz - 0) il>'{iz + 0)
* r ( 0 = (2 In dz.
Jx+ [q + ip(iz-0) q + ip(iz + 0)
(3.81)
If —A_ (resp., A+) is large, |$2"(f)| (resp., | $ J ( 0 I ) is small uniformly in
£ in the upper (resp., lower) half-plane, which can be easily seen from the
explicit formulas for the characteristic exponents and Eq. (3.78) (resp.,
Eq. (3.81)). Hence, we may calculate the integrals in Eq. (3.78) and
Eq. (3.81) with a large relative error and still obtain </>£(£) from Eq. (3.79)
and 0~(£) from Eq. (3.80) with good accuracy. This observation can be
used to develop effective numerical procedures. In fact, even the simple
approximations
97
98 Pricing and hedging of contingent claims of European type
^|*=exp[0X(t)-d(M)], (4-1)
where d(6, t) is a constant. Thus, the density process of the Esscher trans-
form enjoys the following special property: it depends on the current value
X(u), t) but not on the whole prehistory of the process, that is, on X(u, s), s 6
[0, r]; for a general density process, this property may fail. The discounted
price process of the stock S*(t) — e~rtS(t) must be a martingale under Q,
therefore for t > 0, we have
By applying the same consideration to the riskless bond with the dynamics
B(t) = B{0)ert, we obtain
- # P H 0 ) - d ( 0 , t ) = O. (4.3)
p
Prom Eq. (4.3), we find d(8, t) = —ttp (—id), and by substituting into
Eq. (4.2) and dividing by t, we obtain the equation for 8:
Let a solution to Eq. (4.4) exist. Then the Esscher transform exists, and
Eq. (4.1) can be written as
^\rt=exp[eX(t) + WP(-i6)}.
^Q(0 = ^ P ( e - ^ ) - ^ P ( - ^ ) - (4.5)
E x a m p l e 4.1 Let W(t) be the standard Brownian Motion, and let X(t) =
fit + aW(t) be a Brownian Motion with a drift.
Then </>p(£) = ^ £ 2 - */•»£> and Eq. (4.4) turns into
r = - ( - t ( l + 8)a)2/2 + n{l + 9) + (-i6a)2/2 - fi$. (4.6)
Eq. (4.6) has a unique solution 6 — -(/x + CT2/2 - r)/a2, and Eq. (4.5)
gives the well-known formula for the characteristic exponent of a Brownian
motion under the equivalent martingale measure:
tfQ(fl = y £ a - t f ( r - Y ) . (4-7)
Eq. (4.7) can also be obtained from Girsanov's theorem.
For an arbitrary Levy process, Eq. (4.4) may have no real solutions,
and hence, the Esscher transform may not exist. The following lemma
shows that if it exists, it is unique, and provides a necessary and sufficient
condition for the existence of the Esscher transform.
Then
Proof. Parts (ii) and (iii) follow from (i), and to prove (i), it suffices to
show that f'(9) > 0, Vffe (-A+, - A _ - 1). Let pi(x)dx be the law of Xi.
Set
-t-oo
x 6 e OI pi(x)da;,
/ -OO
F(a) = In 1(0, a),
We have
+oo p+oo
F"(6) =
/ oo x2e0x (x)dx
Pl / e8xPl(x)dx
-a
=
+oo
-oo
xe0xpi(x)dx
J—oo
x
r /.+oo
/
U -oo
(I(2,a)I(0,a)-I(l,a)2)/I(0.a)2,
e0xpi(x)dx
Notice that condition Eq. (4.8) is easy to verify once an explicit formula for
ipp is chosen. For processes in the empirical studies of financial markets,
Eq. (4.8) is usually satisfied.
Equivalent Martingale Measures in a L€vy market 101
/
(y/yjx) - l) 2 II(da;) + / (ex - l)y(x)U{dx) < +00, (4.10)
-00 J{x>l]
and
+00
and
r + ipQ(-i)=0. (4.15)
102 Pricing and hedging of contingent claims of European type
Proof, a) and the first part of b) are Proposition 1 in Eberlein and Jacod
(1997); the converse part of b) is from Proposition 2.19 in Raible (2000).
Notice that the integral in Eq. (4.12) converges due to the Cauchy-Schwarz
inequality:
The first factor in the RHS is finite due to Eq. (4.10), and the second factor
is bounded since both II(dx) and y{x)Xi{dx) are Levy measures,
c) We use Eq. (4.12) and Eq. (4.13)
+oo
(l-e»«+^l[.1,1](i)Mz)n(.fa)
/ -oo
/+oo
(1 - eixi + iz£l[-i,i](a;))n(dx)
-OO
+oo
(1 - e*t)(y(x) - l)Tl(dx)
/ -OO
+oo
(1 - e^)(y(x) - l)U(dx),
/ -OO
Then
^ Q ( O = - « X + 0Q(O, (4-19)
where fi is the same as \i in Eq. (3-45) for V»p, and <jfi enjoys the
property Eq. (3.46):
(Hi) there exists v-^<v such that as £ —> oo in the strip 9£ G [fj,-,fi+],
0Q(O=cK|"+O(|Sr), (4.20).
Proof. From Eq. (4.16)-Eq. (4.18) and the conditions on the rate of the
decay of the Levy density II(dx), we see that Eq. (4.10) holds, and (i)
follows from Theorem 4.1. In view of Eq. (3.45)-Eq. (3.46) and Eq. (4.14),
it suffices to show that the RHS in Eq. (4.14) is holomorphic in the strip
9£ € [/*_,/*+], and satisfies an estimate
+oo
(I - eix^(y(x) - l)U(dx) o(ier2). (4.2i)
!/
as £ -> oo in the strip. The analyticity in the strip has been discusses
already, and Eq. (4.21) is proved by using Eq. (4.16) in the same way as
the conditions of the second definition of RLPE were deduced from the
conditions of the first one. (See Chapter 3). D
(The second term makes the left tail of the Levy density fatter). Then
Eq. (4.10) holds, and direct computations (essentially the same as in the
calculation of a K0B0L process in Section 3.1) show that
If
where r = T—t is the time to the expiry; by the change of variable y i-» y—x,
+oo
Pr(y-x)g(y)dy. (4.25)
/
-oo
Let g be measurable, and let there exist a € (A_, A+) such that the function
ga defined by ga{x) := e<rxg(x) belongs to the class Z-i(R). Then we can
define
+oo
e-^g(y)dy
/ -oo
on the line $J£ = a. On the strength of Eq. (3.48), we can apply the
Cauchy theorem and shift the line of the integration in the inner integral
in Eq. (4.27):
+oo p+oo+itr
/ -oo J— oo-\-i(T
Now the integral converges absolutely, and we can change the order of
integration and obtain
-f-00+2<7
e ^ ~ r ( r + ^ ) ) m d t ( 4 2 8 )
/
-oo+i<T
x
Set f(x,t) = F(e ,t). Then by using the notation of PDO, we can write
Eq. (4.28) as
For the call option with the strike price K, the terminal payoff is g(Xx) =
(eXr - K)+, and hence ga e Li(R) for any a < - 1 . Since A_ < - 1 , we
can choose a £ (A_, - 1 ) , and apply Eq. (4.28). The direct calculation gives
r+oo
e - » V - K)+dx
•oo
f+OO
= / ( e ^(i-'«) _ Ke-*x)dx
UnK
e(l-i£)\nK e(l-i£)\nK
Ke-iSlnK
(4.30)
and by substituting into Eq. (4.28), we obtain the pricing formula for the
European call option:
FcMSut) - - ^ L ^ ^ - 5 *, (4-31)
+Fput{St,t).
By taking into account that V>(0) = 0 and r + rp(—i) — 0, since we price the
options under an EMM, we obtain the so-called put-call parity relation
where fi € R, v £ (0,2],i/ ^ l,c > 0, and A+ > 0 > A_ (in the case
v = 1, the formula for the characteristic exponent is different: see Eq. (3.8)).
When v = 2, we obtain a Brownian motion. Parameters are chosen so that
the Esscher transform exists. In each series, we vary the strike-price ratio
S/K, and one of the parameters, the remaining ones being fixed; the spot
price of the stock is 100. To facilitate the comparison with the Black-Scholes
formula, in each series we fix mi = iip'(Q) and m 2 = i>"{0). On the left
panel, we plot the difference between the Black-Scholes price and the KoBoL
price of the put, and on the right panel, we plot the implied volatility.
Were the Black-Scholes model correct, the implied volatility would have
been equal to m 2 , and the implied volatility surface flat. The reader is
advised to compare the pictures below with similar pictures in Eberlein
et al. (1998) and Eberlein and Prause (2001) for option pricing under
Hyperbolic Processes, Normal Inverse Gaussian processes and Generalized
Hyperbolic processes.
108 Pricing and hedging of contingent claims of European type
N_-__ ^ v v 1 , 2 v. ,. c , 1.1
T(daya) ° OS °-9 T (days) ° o.e °9
S/K S/K
Fig. 4.1 Black-Scholes price vs. KoBoL price: dependence on the moneyness and time
to expiry. Parameters: r = 0.05, mi = 0.12, m 2 = 0.18,1/ = 1.4, A_ = -24, A+ = 20.
Fig. 4.2 Black-Scholes price vs. KoBoL price: dependence on the moneyness and
the left steepness parameter, A+. Parameters: r = 0.05,mi = 0.12,7712 = 0.18,v —
1.4, \- = -24, T= 1 (day).
In Fig. 4.1, it is easily seen that the volatility smile becomes more
pronounced as r -> 0. Fig. 4.2 demonstrates that in the region where
the steepness parameters are large, the dependence of the price on these
parameters is weak. This is the good news since the steepness parameters
Pricing of European options and the generalized Black-Scholes formula 109
Fig. 4.3 Black-Scholes price vs. KoBoL price: dependence on the moneyness and the
order of the process, v. Parameters: r = 0.05, m i = 0.12,7712 = 0.18, A_ = — 24, A+ =
20, T = 1 (day).
are especially difficult to infer from the data with a reasonable degree of
accuracy. In the region where the steepness parameters are small, the
dependence of the price on these parameters can be quite sizable. Since A+
and — A_ do not differ much, the smile is almost symmetric. By playing
with steepness parameters, it is not difficult to obtain smiles of different
shapes.
Fig. 4.3 shows that the order of the process, v, is the most important
parameter. As v —> 2 and/or min{A+,— A_} —• +00 so that mi and 7712
remain fixed, the KoBoL price converges to the Black-Scholes price; the
difference is especially large for small v and/or A+ and - A _ .
In the last series of numerical examples, we consider the process which
is characterized by the same parameters as in Fig. 4.1 (in particular, u —
1.4, A+ = 20, A_ = —24), but instead of the Esscher transform, we use the
family of EMM constructed in Example 4.2, with fi+ = 0.5,7 = li a n d
a = 9,6 + 0.05, ...,0 + 0.3. For each value of a, we find d so that the EMM-
condition r + ^(—i) = 0 holds. Thus, a = 6 gives the Esscher transform.
The larger the a, the more weight the EMM places on extreme events.
By comparing Fig. 4.1-4.3 with Fig. 4.4, we see that it is easier to
obtain asymmetric smiles with EMM different from the Esscher transform.
110 Pricing and hedging of contingent claims of European type
Fig. 4.4 Black-Scholes price vs. KoBoL price: EMM from Example 4.2. Parameters:
r - 0.05,mi = 0.12, m 2 = 0.18,1/ = 1.4, A_ = - 2 4 , A+ = 2 0 , ^ + = 0.5,7 = 1,T = 1
(day).
not only the order 1 and the parameter S(= c in the general definition of
RLPE) are the same for P and Q, but the drift parameter \i as well (see
Raible (2000), Proposition 2.20). Since the class GHP is a 5-parametric
family, the EMM-condition leaves 2 free parameters. If one wishes to work
with a subclass of NIG, only one free parameter remains.
Finally, we would like to notice that by using the asymptotic formulas
for the prices of the up-and-out and down-and-out options near the barrier,
which we obtain in Chapter 7 (or similar formulas for barrier options),
and the observed prices of these options near the barrier, we can infer
whether there is a diffusion component or not. Further, if there is none and
hence, the process is of order v 6 (0,2), then one can obtain the following
information about the EMM chosen by the market without using the data
on the stock itself: first, whether v belongs to (1,2), (0,1) or {1}, and then
(i) if v e (1,2), then to find v;
(ii) if v — 1, then to obtain a simple relation between the parameters c
and (j, in Eq. (3.45)-Eq. (3.46);
(iii) if v € (0,1), then decide whether the drift /x = 0 o r not, and
• if fi = 0, then to find v\
• if n ^ 0, then to find the sign of fi.
Thus, the data on the touch-and-out options and barrier options near the
barrier can be used to decide whether the chosen class of processes is suit-
able for pricing of options on a given stock or not.
For more details, see the next section.
= (r + 1>(Dx))f(x,t).
We conclude that the price / of a contingent claim of the European type
satisfies the following equation
We consider several examples; the reader can easily calculate the prices of
other options.
_ exp[—i£,lnK]
Similarly, for the digital put option with the strike price K, the terminal
payoff is 1, if ST < K, and 0 otherwise, and we obtain
+00+t<72
eiSMSr/K)^-!^ (438)
/
-00+l<72
for any er2 € (0, A+).
4.4.3 Combinations
There are several types of options whose payoffs are linear combinations
of various calls, puts, and constants. For instance, a straddle is the option
with the payoff of the form
e^-lR±(-M-)€Li(R).
For any XQ, decompose the terminal payoff g in the sum of the following
two functions
•oo+iu+
4.5 Hedging
mmE[(H-Vr'6)2}
found in the references at the end of the chapter. For hedging under pro-
cesses of Koponen's family, see Bouchaud and Potters (2000) and Matacz
(2000).
We have chosen the locally-risk minimizing hedging in order to demon-
strate in the simplest situation how the PDO-technique allows one to obtain
analytical results.
where
and
L = dt + Lp =dt-ipP(Dx). (4.44)
= (dt-i>P(Dx))f(x,t)exAt + o(At),
and
= (dt-ipp(Dx))f(x,t)At + o(At).
Since
and
+ 4>p{-i)]F(ex,t) + o{At).
Now we can substitute the result and Eq. (4.43) into Eq. (4.42), and obtain
B (4 48)
«) _^(_M)+2*P(-0 • -
<KO = - ^ + <72£2/2,
and therefore,
Hence,
and
9(t,x) = (dxf)-e-x=dsF(S(t),t).
Ciu + Kir1-',
where e > 0, and C\ can be chosen the same for each t <T and S/K on any
compact. Since the integrand depends on t continuously, we conclude that
9(t,x) is continuous (in fact, Holder continuous) on (—oo,+oo) x (—oo,T].
This makes the hedging under purely discontinuous RLPE models more
stable than the Gaussian hedging, near expiry; and at the strike, much
more stable indeed.
Fig. 4.5 Black-Scholes 6 vs. KoBoL 8: dependence on the moneyness and time to expiry.
Parameters: r = 0.05, m i = 0.12,m 2 = 0.18,1/= 1.4, A_ = - 2 4 , X+ = 20, K = 100.
4.6 Commentary
121
122 Perpetual American options
K(x)=sup^[e-«rfl(XT)]. (5.1)
M
5.1.2 Free boundary value problem for the price of the per-
petual American option
Suppose that the optimal stopping time is the hitting time of the exterior
of an open set C C R:
On the strength of Theorem 3.1 and Lemma 2.4, an RLPE satisfies (ACP)-
condition, hence, in the case of g € Loo(B), where B = R \ C , we can apply
Theorem 2.12 and conclude that V, is a bounded solution to the following
boundary value problem:
formula for u = u(h,x), and by analysing the properties of the latter, find
the unique candidate for the optimal h. The verification of the sufficient
optimality conditions in the class MQ is relatively simple: we will be able to
directly compare V — V(h, x) for different h, and check that the candidate
for the optimal boundary is indeed optimal.
Lemma 5.1 Let (C,V) be a solution to Eq. (5.3)-Eq. (5.6) such that V
is continuous, let r* be the hitting time of C, and let
From Eq. (5.10), Eq. (5.6) and Eq. (5.3), we conclude that for any stopping
time r,
qT
V(x) > Ex e' V(XT)] , (5.11)
124 Perpetual American options
a.e. (Eq. (5.6) and Eq. (5.3) say nothing about W at dC but this is
unnecessary: dC has the Lebesgue measure 0, and therefore for a.e. x,
U«{x,dC) = 0 - see [S], Proposition 41.9). Prom Eq. (5.10) and Eq. (5.3),
for a chosen stopping time r»,
By using Eq. (5.5) and Eq. (5.4), we deduce from Eq. (5.11) and Eq. (5.12)
Theorem 5.1 The optimal exercise region for the perpetual American
put is S < H* or equivalently x < h*, where
H*=eh'=Kcj>;(-i), (5.13)
v
*(x) = ~ « r / etc.* — d
^ (5-14)
27r J_00+iu+ £(£ +1)
for any u>+ € (0, a+).
Example 5.1 Let X be a Brownian motion, let the volatility be a2, and
let there be no dividends, that is, q = r. Let /?_ = —2r/cr2 be the negative
Perpetual American put: the optimal exercise price and the rational put price 125
y/32 + ( r - y ) / 3 - r = 0
?lv"
2
+K(r-—)V'-rV
2'
= 0.
Then from Eq. (2.28), </>-(£) = ( - £ _ ) / ( - / ? _ + i£), and Eq. (5.13) and
Eq. (5.14) assume the form
and
K
V M - /+0O+iu,+ exp[i(x-/i,)Cl(-/3-),,
respectively. The condition w+ G (0, o+) implies that the line Q£ = u+ lies
below the pole £ = — i/3- of the integrand, therefore by pushing the line
of the integration up to the infinity, we cross this only pole in the upper
half-plane. By applying the residue theorem, we derive from Eq. (5.16)
Kexp\p-(x-h.)}
V.{x) = j - ^ - ,
or
K(x) =
G3c)' (-^-r^
yr^j:) • ( 5 - 17 )
Eq. (5.15) and Eq. (5.17) are due to Merton (1973), who derived them by
using the ODE technique and smooth pasting principle.
Remark 5.2 On comparing Eq. (5.13) and Eq. (2.26), one is tempted to
write Eq. (5.13) in the form
(Uf4(q-L)g)(x) = Ot (5.19)
We will also see that the general form for Eq. (5.14) is
(1) show how to guess the unique candidate for the optimal stopping
time in the class Mo of hitting times -r(a) of segments (—oo, a],
(2) show how to prove, relatively simply, for any RLPE, that this is
indeed the optimal stopping time in the class Mo, and then
(3) produce more involved proof of optimality in the class M, for any
mixture of independent model RLPE,
(4) discuss the smooth pasting condition, and
(5) obtain simple approximate formulas.
After Theorem 5.1 was obtained, Mordecki (2000) proved it (in the form
Eq. (5.18)) for any Levy process by using the different technique. Notice
that in Section 5.5 we manage to prove the optimality of Eq. (5.19) in
the class Mo for wider class of payoffs, and the more practically oriented
reader may be satisfied with a straightforward proof of the optimality in
the intuitively natural class Mo.
V(h,x):=Ex[e-«rWg(XrW)]. (5.21)
Lemma 5.2 Let there exist /i, with the following properties:
V(h,x)<g(x); (5.22)
V{K,x)>V(h,x). (5.24)
Proof. Clearly, the rational price of the option must satisfy Eq. (5.23),
hence Eq. (5.22) excludes h < h*. Due to Eq. (5.23), h* is an admissible
choice, and Eq. (5.23)-Eq. (5.24) ensure that a choice h > h* is no better
than h*. •
To apply Lemma 5.2, we need an explicit formula for V(h, x). We derive it
by applying the Wiener-Hopf method to the problem Eq. (5.3)-Eq. (5.4),
which now assumes the form
Set
w(x) = ((/>-(D)-lg)(x).
Theorem 5.2 For any h G R, a solution to the problem Eq. (5.3)-
Eq. (5.4) in the class of measurable functions bounded on [h, +00) exists,
and
V = <j>-{D)l(_ooM<l>-{D)-1g; (5.27)
(for the sign " + " , it is immediate from the definition of <£~, and the formula
for the sign "-" follows), and <^"(0) = 1. Second, write Eq. (5.27) for x > h,
as
and then as
r+oo
V(h, x) = g{x) - p~{x- y)w(y)dy. (5.30)
Jh
h = ln(K<p-(-i)). (5.31)
If w(h) < 0, then w{x) < w(h) < 0 for each x > h, and therefore from
Eq. (5.30), we conclude that for each x > h,
V(h,x) >g{x)-w(h).
v(h,$) = fh e-ixt(K-ex<f>;(-i))dx
J—OO
V(h,x) = ct>-(D)v(h,x).
Perpetual American put: the optimal exercise price and the rational put price 129
V(h,x)<V(h,x). (5.32)
But for y > h, w(y) < 0, and therefore from Eq. (5.29), we conclude that
V(h,x)-V{h,x)= f p-{x-y)w(y)dy<0.
Jh
W: - {q-L)V
= qiq-'iq + iPm^iDWh,-)
X
= q^{D)- v{h,.).
Under the choice Eq. (5.31), v(h,£) decays as | £ | - 2 , as ^ —> oo, therefore
in the case K+ < 1, the Fourier transform of W admits a bound via C(l +
K l ) - 1 - 6 , where e > 0. Hence, W € C 0 (R), and Eq. (5.7)-Eq. (5.8) hold.
If K+ = 1, then from Eq. (3.73)-Eq. (3.76), we deduce the representation
^(0 _ 1 =9- 1 (9 + # ( 0 W 0 .
130 Perpetual American options
We formulate sufficient conditions for Eq. (5.34) to hold; they are satisfied
by any model RLPE, and any mixture of independent model RLPE. Of the
process, we require
(i) the function (f> in Eq. (3.46) admits the analytic continuation into
the lower half-plane with the cut (—ioo, iA_], and admits the bound
Eq. (3.51) in this half-plane, outside a neighbourhood of iA_; and
if v = 2, there exist c and u\ < 2 such that <j>{tz) — c£2 satisfy
Eq. (3.51) with v\ instead of v\
(ii) in a neighbourhood of zA_, <j> may have a weak singularity:
\m\<C\Z-i\-\-a, (5.35)
for some a < 1;
(iii) for any z € (—oo,A_), the limit
W(+0) = W(0)
-)-ooi-ia di
/ (0(-t^)(i-i0
•oo+iff *Pq
The integrand is holomorphic in the upper half-plane 3?£ > 0 and admits
an estimate via C(l + \£\)~2+K+, for 9£ > a > 0. Hence, we can push the
line of integration up, and in the limit a —• +oo obtain zero. This finishes
the proof in the case K+ < 1.
If K+ = 1, we can represent (q + V>(£))<?V(£) = # J ( £ ) _ 1 m the form
(q + tf(0 W 0 = 9 a + ( 0 ) - x ( - i 0 + X (€),
J —oo+io -
•
5.2.4 Failure of the smooth pasting principle for some
RLPE's and its substitute
Theorem 5.3 Let (j>Q satisfy
9
+oo+icr
107(0(1-*0_1|de<+oo, (5.38)
/. oo+ia
Perpetual American put: the optimal exercise price and the rational put price 133
for some a € (0, <r+). Then the price of the perpetual American put satisfies
the smooth pasting principle.
** J-oo+io
I ix(o(i-*o_1|de<+oo.
—oo+icr
Then the smooth pasting principle fails.
which is discontinuous. •
Notice that for an RLPE, Eq. (5.38) fails if and only if [i > 0 and v € (0,1);
this is the case of a process of bounded variation, with positive drift.
As we already saw, the natural candidate for the optimal exercise price is
determined from the equation
w(x):=<l>-(D)-1g{x)=0,
and this candidate can be singled out formally in one of the following forms:
The second principle works for purely non-Gaussian RLPE, i.e., for RLPE
of order v < 2.
In all cases, one may say that the optimal choice of h makes V(h, •) "more
regular" at h than generically.
Hv = eh*>=K/3-/(l3--l),
and
respectively, which look exactly like Eq. (5.15)-Eq. (5.16) in the Brownian
motion case. However, if one fits an RLPE by a Brownian motion, one
obtains the different /3_, and the different optimal exercise price and ra-
tional put price. Generically, the larger the A + , the smaller the differences
between the exact solution (ft#,K) and the approximation (/i a p,K P ), and
between (/i*, Vr) and the Gaussian result.
a(D):=q + i>(D)=q-L,
substitute V(x) = e'^V^x) into Eq. (5.25), after that multiply Eq. (5.25)
and Eq. (5.26) by e 7 *, and use the equality
We obtain
\g\a\x)\<Cse-^, (5.42)
for each s = 0 , 1 , 2 , . . . (for the next steps of the proof, Eq. (5.42) for s < 2
suffices). Construct G 7 which coincides with g 7 on R_ and admits a bound
Eq. (5.42) on R, for s < 2, and set u 7 = Vy - G 7 , F1 = -a(D + ry)G 7 .
Then u 7 solves the problem
Eq. (5.42) implies that G1 e H2(R), and from Eq. (3.45)-Eq. (3.46) we
conclude that F 7 € i 7 2 _ " ( R ) , where P = i/, if i/ > 1 or /i = 0, and P =
maxfi/, 1} otherwise. Recall that we are looking for V which is measurable
and bounded on (0, +oo). Hence, we are looking for « 7 which is measurable
and admits a bound via Celx. We want to reduce the problem to the case
of an unknown function of the class L2CR+). Since cr_ < 0, we can choose
7' 6 (<T_ — 7, —7). Set
substitute u^(x) — e - 7 x u 7 i 7 '(x) into Eq. (5.43) and Eq. (5.44), and after
that multiply Eq. (5.43) by e 7 x. By using Eq. (5.39), we obtain
Hence, the expressions on both sides of Eq. (5.45) belong to i J - < / ( R ) , and
the support of their difference, call it F-, is a subset of (—00,0]. This
o
means that F- eH _ 1 / ( R - ) , and we can write the Wiener-Hopf equation
Eq. (5.45) in the form
a(D +1(7 + 7 ' ) K , 7 ' = -F7,7' + F-. (5.47)
136 Perpetual American options
Multiply by q~l, and then apply <f>^{D + 2(7 + 7')). Since in the strip
3£e(A_,A + )(D(<7_,<r + )D(<7_,0))
?-1a(0 = ^(0-V7(0- 1
and by our choice, 7 + 7' G (cr_, 0), we obtain
^ - g - V + P + i d + T'))^-
By construction, G 7 € i / 2 ( R ) , and
u 7 , y e L 2 (R+) = t f ° ( R + ) . F_ € f f - p ( R _ ) .
Prom Theorem 3.3, we know that for any a e (CT_,CT + ),
Ul = -4>-{D + n)9+cl>-(D + i 7 ) _ 1 G 7 .
= 0-(£> + t 7 )0_0-(L> + » 7 ) - 1 a T .
Since G 7 coincides with g~, on R_, Eq. (3.74) and Eq. (3.76) ensure that
supp <j>-{D + J 7 ) - 1 ( G 7 - 9-y) C [0, +oo). Thus,
and in the formula for Vy, we may replace G 7 with g 7 . By using Eq. (5.39),
we finally arrive at
V = 4,~{D)9-<l>-{D)-1g, (5.51)
which is Eq. (5.27). This proves the first two statements in (ii)-(iii), in
case a). To study the continuity of V, consider w := 4>~{D)~xg. Clearly,
w € H2~K~ (R), and since 2 — K_ > 1/2, we can apply Theorem 15.15, and
obtain
V = w(0)V1+V2, (5.53)
where
^^-(DXl-iD)-1*,
and
V2 = 4>"(D)(1 - iD)-l9.{\ - iD)4>-{D)-lg.
Notice that for any e > 0, 8 £ H'^-^R), and 0_(1 ~ iD)^ {D)~l g &
H°{R). Hence, V2 6 H\R). If K_ > 0, we obtain Vx e H 1 / 2 + P ( R ) J for
any p G (0, K_), and therefore, the same is true of V. But for s > 1/2,
i P ( R ) C Co(R), and therefore, V is continuous. This finishes the proof of
(ii) for the case K_ 6 (0,0.5).
138 Perpetual American options
/ -oo
Theorem 5.5 Let the dividend rate be positive, and let X be a regular
Levy process of exponential type [A_, A+], where A_ < —1 < 0 < A + . Then
the optimal stopping time is the hitting time of the interval [h*(K), +oo),
where
Remark 5.3 a) As in the case of the put, our method proves the optimality
in the class .Mo, for any RLPE, and in the class M, for any mixture of
independent model RLPE.
b) After Theorem 5.5 was obtained, Mordecki (2000) proved that it holds
for any Levy process satisfying E[eXt] < +oo, in the form
c) Our proof shows that the natural general form for Eq. (5.56) is
E[({q-L)g)(NT)) = 0, (5.59)
where g(x) = ex — K, and N is the infimum process, and for Eq. (5.57),
the general form is
In the next Section, we will obtain Eq. (5.59)-Eq. (5.60) for more general
payoffs.
d) Prom the point of view of the analytical technique, which we use, the
principal difference between the case q > r and the case q = r is that in the
former case, the condition Eq. (3.40) holds with u>_ < — 1 < 0 < u+, and
in the latter case, due to the EMM-condition r + i/)Q(—i) = 0, it holds with
u;_ > — 1 only. x A less important issue is the unboundedness of the payoff,
which forces one to approximate the payoff by a sequence of bounded payoffs
{gn} in order that Theorem 2.12 be applicable (cf. Remark 2.1). Modulo
this approximation, the proofs and statements for the put can easily be
transformed into their analogs for the call by changing the direction on
the real axis and the reflection of the complex plane w.r.t. the origin. In
particular, the essentially new elements appear only in the proof of the
analog of Eq. (5.27).
V(h,x)<g(x); (5.61)
1
To be more specific, in the former case, the real part of the operator q — L = q + t/>Q (D)
is positive-definite on the L2-space of function which grow as ex as x —• +00, whereas
in the latter case, e x is an eigenfunction of q — L.
Perpetual American call 141
V(h*,x)>g(x); (5.62)
V{h*,x)>V(h,x). (5.63)
Since gn is bounded, we can apply Theorem 2.12 and conclude that V(gn; ft, •)
is a bounded solution to the following problem
Further, by replacing all the signs in the statement and the proof of Theo-
rem 5.2, we find that the unique bounded measurable solution (and in the
case K + = 1, the unique continuous bounded solution) is given by
V(gn; h, •) = ^ ( l O l ^ + ^ + t D ) - V - (5.67)
By using Eq. (3.72) and Eq. (3.76), it is easy to obtain
'i-(h,+oo)(t>q{D)~1(gn ~ l(/i,+oo)5n) = 0.
The result is
V(gn; h, •) = e-*x4>t{D + t 7 ) l ( h i + o o ) # ( 0 + i 7 ) _ W 1
; •)• ( 5 -™)
V(gn;h,-)->W(g;h,-), (5.71)
where
W(g; h, •) = <t>+{D)\(h}+oo)4>${D)-lg.
Compare Eq. (5.64) with Eq. (5.71). Since the limit in the sense of gen-
eralized functions is unique, we have W(g; h, •) = V(g;h,-) in the sense of
the generalized functions. By explicitly computing W(g; h, •):
+00+17
^•(-*)_1e(1-^ Ke-*h
e f a C #(0 d£
/ -oo+it 1-tf -^
We set w(x) = <j>+{D)~lg{x), and prove that the optimal h* in Lemma 5.5
is the solution to the equation w(x) = 0, that is, <f>%(—i)~1eh' — K = 0.
This and the rest of the proof can be done in essentially the same manner
as in the case of puts.
Put-like and call-like options: the case of more general payoffs 143
where w'_ < LJ'+ < a+. Then all the statements of Theorem 5.2 hold.
Proof. Eq. (5.73)-Eq. (5.74) are the very conditions we used in the proof
of Theorem 5.2, with only one exception: if tv'+ > 0, then in the proof, we
must take 7 S (LJ'+,O-+). •
Remark 5.4 Notice that Eq. (5.74) is not a real restriction: if g fails to
satisfy this condition, we can replace g with any Ig, which coincides with g
on (—00, h] and satisfies Eq. (5.74). In particular, for any such Ig,
V = rq{D)h-ooM^{D)-Hg, (5.75)
which is understood as
+O0 + tUJ+
eix%{0mdi, (5-76)
/
-OO-HCJ-L
1
where v = l(_ oo/l ](/) q (D) lg, and w + S (u'+,a+) is arbitrary.
1
Setw(x)=(4>g(D)- g){x).
T h e o r e m 5.7 Let Eq. (5.73)-Eq. (5.74) hold with LJ'_ < u'+ < a+ and
some h, and let there exist h\ < Ji2 such that the following conditions are
satisfied:
Then for any h G [hi, hi], T(/I) is an optimal stopping time in the class
Mo, o,nd the rational option price is given by
e fa V7(0«(Ode, (5-80)
where -o+ < 71 < 72 < • • •• Then Eq. (5.74) and Eq. (5.73) hold, and in
concrete cases, conditions Eq. (5.77)-Eq. (5.79) are easy to verify since
In particular, if they are satisfied then hi = /12; call it ht. Notice that
the necessary condition for Eq. (5.77) and Eq. (5.79) is ci > 0 and c; < 0,
respectively.
Example 5.2 Suppose, the option owner has the right to sell a share of
the stock for K + ay/S, where S is the spot price. Then
g(x) = K + aex'2 - ex,
w(x) = tf + a < £ - H / 2 ) - V / 2 - 4 > - H ) - V ,
and the optimal exercise boundary is h* = y/Y, where Y is the unique
positive root of the equation
K + a 0 - ( - i / 2 ) " 1 r - ^(-iy'Y2 = 0.
When /i* is found, one can easily calculate v = l(_0o,/i»]u; a n d *K£)> f° r a n v
£ in the strip 3£ e (0,o+):
and after substituting into Eq. (5.80), obtain the rational option price:
-t-0O-|-JU>+
e i(z-h.)^-(0
/ -OO-f
-OO-f 1W+
ILJ+
K ae05h- eh'
X I-TT
-i£ + < ^ H / 2 ) ( 0 . 5 - i £ ) ^-(-i)(l-^)J
where u>+ G (0,0+) is arbitrary. By using the approximation Eq. (3.82) for
<^>~(£), we obtain the approximation to the rational option price:
f+oo+iu+ _a
e<(x-M€.
•oo+iw+
„„0.5/u ph. 1
+
-t£ ^(-i/2)(0.5-iO ^(-0(1 ~^)J
The integrand has the only pole £ = — i/3_ hanging above the line of inte-
gration, and this pole is simple. By pushing the line of integration up and
using the residue theorem, we obtain the following approximation:
where
-/?_ae°- 5/l - ~(3-eh'
D =K +
&-(-x/2)(0.5-/3_) ^-(-t)(l-/J_)
The next theorem shows that for model RLPE, the stopping time obtained
in Example 5.2 is optimal not only in the class Mo but in the class M as
well.
i i
x c ex
9( ) =Y^ t vllk } ~^2ck exphk*]'
x
( 5 - 83 )
fe=i fc=i
146 Perpetual American options
where cjj: are positive, 7 ^ > — cr+, k = l,...,l, are not necessarily
different, as well as 7^" > — o~+, k = 1 , . . . , /, and satisfy
1 h 1 1 1
ct^{-iltr ^ =^{-ilkT ^ (5-85)
Then r(/i*) is the optimal stopping time in the class M, and the rational
option price is given by Eq. (5.80) with any w+ € (maxj(—7,), cr+), h — /i»,
and
and Eq. (5.85) implies that h, = h*(gj) is the same for all j . But when /i*
is the same for all payoffs gj, Eq. (5.6) is evidently "additive" w.r.t.
g in the sense that if it holds for all V*(gj,x) defined by gj, then it holds
for V*(g,x). Thus, it remains to consider the payoff Eq. (5.87), and we
can repeat the proof of Eq. (5.6) for the put with the straightforward
changes. D
Let cr_ < 0 be from Lemma 3.5, and consider the problem
Theorem 5.9 Let g and its derivatives g(3' = D3g, s = 0,1,2, be measur-
able, and satisfy the estimates Eq. (5.73)-Eq. (5.74) with o~- < w'_ < UJ'+ .
Then for any h € R, o solution to the problem Eq. (5.88)-Eq. (5.89) in the
class of functions bounded on (—oo,/i] exists, and
V = tfiD^+^+iDr'g; (5.90)
Remark 5.5 Notice that Eq. (5.73) is not a real restriction: if g fails to
satisfy this condition, we can replace g with any Ig which coincides with g
on [h, +oo) and satisfies Eq. (5.73). In particular, for any such Ig,
which is understood as
+oo+iw_
e^+TOSK, (5.92)
/
•oo+iui_
where v — l[h^+00)<j)+ (D)~xlg, and w_ € (cr-,uj'_) is arbitrary.
Set
w(x) = (<^+(Z))- 1 g )(x).
Theorem 5.10 Let Eq. (5.7S)-Eq. (5.74) hold with a~ <u'_ < u'+ and
some h, and let there exist hi < hi such that the following conditions are
satisfied:
Then for any h € [^1,^2], T(JI) is an optimal stopping time in the class
Aio, and the rational option price is given by
e 1 ( l - f c ) c #(0*(Ode. (5-96)
Notice that Eq. (5.94) can be written in the form Eq. (5.59).
Consider payoffs of the form Eq. (5.81), where 71 < 72 < • • • < —<r_. Then
Eq. (5.74) and Eq. (5.73) hold, and in concrete cases, conditions Eq. (5.93)-
Eq. (5.95) are easy to verify since
l
w{x)=YJCj<t>t(-ilirl^jX- (5-97)
In particular, if they are satisfied then h\ — J12; call it h*. Notice that
the necessary condition for Eq. (5.93) and Eq. (5.95) is ci < 0 and c; > 0,
respectively.
where cj~ are positive, 7^" < — <J_, k = l,...,l, are not necessarily
different, as well as 7^ < — c _ , k — 1 , . . . , /, and satisfy
Then r(/i*) is the optimal stopping time in the class M, and the rational
option price is given by Eq. (5.96) with any <J+ G (<r_,minj(—jj)), h = h*,
and
7
HO = -i24<t>i(-nf)-^t-w + ' "7£ (5.101)
5.5 Commentary
In the finite horizon case, there are no explicit analytical pricing formulas
even in the Brownian motion case. We formulate the free boundary value
problem for the price of an American option, and consider several analogues
of approximate methods used in the Gaussian option pricing theory. We
also derive an approximate formula for the early exercise boundary near the
maturity and show that its behaviour in the RLPE-case drastically differs
from the behaviour in the Gaussian case.
where supremum is taken over the set of stopping times satisfying T(U/) <
T.Vwefl.
Let L = L® be the infinitesimal generator of X under Q. Assume that
the optimal stopping time is of the form r'B AT, where r'B is the hitting time
151
152 American options: finite time horizon
Then B is the optimal early exercise region, and V* is the rational option
price.
Proof. Repeat the proof of Lemma 5.1 with the process X instead of X,
taking into account that a sample path of X cannot cross the line t = T by
a jump. •
For the detailed study of the free boundary problem for the American
put in the Gaussian case, see monographs Musiela and Rutkowski (1997),
Karatzas and Shreve (1998), Shiryaev (1999) and the bibliography in op.
cit.
r + y>QH) = 0 (6.5)
for the characteristic exponent ip** of the process X, under Q. If the stock
pays dividends at the constant rate A > 0, then {e~(r~x)tSt} is a martingale
Approximations of the American put price 153
r - A + V>Q(-i)=0. (6.6)
r + ^Q(_j)>0; (6.7)
and from the technical point of view, this is the difference between Eq. (6.5)
and Eq. (6.7) which leads to Merton's (1973) result on the no-early-exercise
of the American call on a non-dividend-paying stock. Merton's no-arbitrage
argument is universal; in particular, it is valid not only in the Gaussian case,
but also for any Levy process.
Similarly, if the interest rate is zero: r ~ 0, then it is non-optimal to
exercise the American put before maturity.
Let Vj(x) = V(ex,tj) be the value of the option at time t = tj. For
j = n, when tn = T, vn is the terminal payoff: vn(x) = max{.K" — e x , 0 } .
At time tj, j = n — 1,n — 2 , . . . , 1, the option owner chooses the optimal
exercise boundary hj as the solution to the equation
where Tj = tj+i — tj. In Eq. (6.8), the LHS is the payoff at the current
level of the log-price, x, and the RHS is the discounted expected value
of keeping the option alive. If an analytical formula for the probability
density pt of the process Xt under Q is known, as is the case with NIG, or
tabulated, which can be done for any model RLPE, then Eq. (6.8) can be
solved numerically, and after the optimal exercise boundary hj is found, Vj
is calculated as
In the Brownian motion case, it is known that Eq. (6.8) has a unique
solution, and the sequence {VJ} converges to the option exact price when
A :— maxTj —* 0; one should expect that the same holds in the RLPE-case
as well. To estimate the limit, one may use Richardson's approximation
scheme (see, e.g., Musiela and Rutkowski (1997)).
for j = n — 1, n — 2 , . . . , 0,
or
a A (0 = (l + rA)^+(0-V7(0 _1 .
where <j>^ is the Fourier transform of a probability distribution P^ (dx) sup-
ported on R T (see Chapter 5), and solve the problem Eq. (6.9)-Eq. (6.10)
by using the Wiener-Hopf factorization method.
Denote by u,i the solution to the problem Eq. (6.9)-Eq. (6.10) with
Vj+\ = 0. Theorem 5.2 gives
v
ji =<t>q{D)l(-00,hi}<l>q{D)~1g,
and reasoning the same way as in the proof of Theorem 5.2, we can show
that the solution to the problem Eq. (6.9)-Eq. (6.10) with g = 0, call it Vj2,
is given by
vj2 = (1 + rA)-1<j>-(D)l[hjt+oo)<l>+(D)vj+1.
Vj = ^•(D)l(_OOifcil07(D)-1«;
+{l + T^)-1d>-{D)llhii+eo)^+(D)vj+1
+(l+rA)-1<t>-(D)l[hjt+oo)cl>+(D)vj+1
= (K - ex) + (1 +rA)-1<t>q-(D)l[h.<+oo)wj,
156 American options: finite time horizon
where
Wj = 4>+(D)vj+1-(l + rA)<f,-(D)-\K-e')
= <l>+(D)vj+1 _ ( l + r A ) ( ^ - ^ ( - i ) - V )
= 4>+(D)[vi+1 + e'-(l + rA)K].
where
(Eq. (6.12) will be proved at the end of the proof of Lemma). By using
Eq. (6.12), we derive
< 4>+(D)K-K<K-K = 0,
and
Vj{x) = Ex
f
Jo
e-^A^Vj+^X^dt + EX ,-QTi
(K ex^)
where r, is the hitting time of (—oo, hj}. Since vn(x) — TOSX.{K — e x ,0} G
[0, K), we find that vn_i is non-negative, and it admits an upper bound
via
that is, vn-i satisfies Eq. (6.12). By iterating, we derive Eq. (6.12) for
j=n-2,...,0. •
(and for some x, the inequality is strict), which finishes the proof of Theo-
rem. D
To sum up: the algorithm for the discrete time approximation to the opti-
mal exercise boundary and the rational put price is as follows:
1. Set hn<* = h\K, vn(x) = max{K — e T ,0};
2. For j = n — 1, n — 2 , . . . , define
Wj = <j>+(D)[vj+1 + ex - (1 + rA)K],
Wj (x) = 0,
and set
Vj(x) =K-ex + (l+ rA)~V7(I?)l[fcil.,+«x,)«'j(a:)-
Vj(x) = ma,xV(vj+i,hj-,x),
hj
where V(v, h; x) is the value of a barrier option with the terminal payoff v,
the down-and-out barrier h, and the rebate g (we may assume that only
American put near expiry 159
barriers h s.t. g(x) > 0 for all x < h are considered); the maturity date, r, is
a random variable which is exponentially distributed with scale parameter
A = 1/A. Hence,
r+oo
Vj(x) = maxA / e~xtU(vj+i,hj-,x,t)dt,
h
i Jo
where U(VJ+I, hf,x, t) is the value (at time 0) of a barrier option with the
maturity date t, the terminal payoff t>j+i, down-and-out barrier hj, and
rebate g. Prom Theorem 2.13, we conclude that U(vj+\,hj;x,t) is equal
to W(x,0) (the arguments Vj+\,hj,t are suppressed), where W(x,s) is a
bounded solution to to the following boundary value problem
Divide Eq. (6.13) by A, and take into account that A - 1 = A; we see that the
problem Eq. (6.13)-Eq. (6.14) turns into the problem Eq. (6.9)-Eq. (6.10),
and as in the preceding subsection, the free boundary hj must to be chosen
so as to ensure that the solution is optimal. The result is the procedure
described at the end of Subsection 6.2.2
' T h e Laplace-Carson transform differs from the standard Laplace transform by an ad-
ditional A-factor.
160 American options: finite time horizon
as t —> T—0 (see, e.g., Barles et al. (1995) and Lamberton (1995)). We have
been unable to study the behaviour of the optimal early exercise boundary
near the expiry in the proper continuous time setting; instead, we study the
asymptotics of the boundary point /i„_i = h(A) at the last moment before
expiry, as A —» 0, in the discretization of the continuous time model. We
believe that the asymptotics of h(A) can be used as a good proxy for the
asymptotics of the boundary near expiry.
We study NTS Levy processes, NIG in particular; similarly, KoBoL and
Hyperbolic processes can be treated. The main result is that in many cases,
in particular, for processes of order v G [1,2), the limit of h(A) is negative,
and not zero, as in the Gaussian case. Notice the opposite effect for the
boundary far from expiry: the numerical results for the optimal exercise
price h» of the perpetual American put show that typically, the Gaussian
model gives smaller /i*, and we may conclude that in the non-Gaussian
case, the optimal exercise boundary is more flat than in the Gaussian case.
This observation shows that there is some compensation for the apparent
difficulties with the design of appropriate numerical procedures for pricing
of the American put in the non-Gaussian case.
For processes of order between 0 and 1, the behaviour of h(A) of various
types is possible, depending on all the parameters of the process; we will
not list the corresponding formulas here.
4>+(D)[(l-e°)++e*-l-rA}=0.
Rewrite it as
< # ( Z ? ) [ ( e * - l ) + - r A ] = 0,
0 + ( 0 ) . l = # ( O ) - l = l,
American put near expiry 161
^(0 =
(9 + ^(O)0,-(O'
we continue the transformation of the equation for a: = h(A):
rA. (6.17)
oo+iW_ (9 + V » ( 0 ) ^ ( 0 ( - ^ ) ( l - i O
In the NTS Levy case, A± = ±a + 3, and the formula for the characteristic
exponent of an NTS Levy process is
and
We conclude that
2
[^(iz + 0) - V(w - 0)]* = -a2),//2(e-i'r^2-ei7rl//2)z
<5((/3-z)
= 2Jsin^((^-z)2-a2r/2, (6.20)
and
(q + iP(iz + 0))(q + i>(iz-0)) = (9 + / i 2 - t f ( a a - / 3 2 ) " / 2 ) 2 (6.21)
2 2 2 v
+<5 ((/3 - z ) - a )
162 American options: finite time horizon
+25cos^((/3-z)2-a2y/2
2 2
.(q + ^ z - S i a -^ ).
a) i / € (1,2);
b) u = 1;
c) v G (0,1) and \x = 0;
d) v G (0,1) and ^ < 0;
e) v G (0,1) and /* > 0.
(6 24)
/ z(z + l) ^ - J s i n f "
Before proving Theorem 6.3, we make two remarks. First, the solution to
Eq. (6.24) exists, and it is unique since the LHS decreases monotonically
American put near expiry 163
Proof. First we consider the case u e (1,2). Suppose that the following
two assumptions hold:
a) h = /i(A) ->• 0 as q = A - 1 + r -> +00;
b) — h(A) > Aq~xlv, where A > 0 is independent of q.
Then on (2h~1,h~1), we have e - 1 2 < e - 2 , q+ip(iz±0) x q, and <f>~ (iz) x 1
due to Eq. (6.22) 2 Therefore, the LHS in Eq. (6.18) is bounded away from
zero by
,i/h
c {-h)2-vdz=C]_{-h)l-v,
J2/h
2
For two functions / and g we write / x g if / = O(g) and vice versa.
164 American options: finite time horizon
Now, let v = 1. Then the LHS in Eq. (6.18) admits a bound from below
via
fq _ fA dz
I exzz~1dz = / —z r - —• +oo as x —• - 0 ,
7i 7-1 e z'
and the proof finishes as in the case v G (1,2). •
Chapter 7
First-touch digitals
7.1 A n overview
165
166 First-touch digitals
the case when the dynamics of the stock price is modelled as the geometric
Brownian motion, see, e.g., Ingersoll (2000); general formulas for the RLPE-
case were obtained in Boyarchenko and Levendorskii (2001a, 2002b), where
more general rebates than here were considered, namely, first-touch power
options, which pay not $ 1 but S&; and first-touch contracts which pay a
non-zero amount iff the first barrier has been crossed but the second one
(situated farther) has not. Notice that the contracts of this sort do not
make sense in the Gaussian theory since a continuous trajectory stops at
the barrier.
In addition to exact albeit complicated formulas, we obtain several sim-
pler approximate formulas. The first formula is valid near the boundary
H, and the next two are valid for options far from the terminal date. From
the asymptotic formula for the price of the first-touch digital as the price of
the underlying approaches the barrier we conclude that if there is no Gaus-
sian component and the process is purely discontinuous then the difference
between the option prices in the Gaussian model and in the RLPE model
is especially large near the boundary. Moreover, different RLPE models
usually produce essentially different types of behaviour, which can be used
to decide, which type of processes fits the data better.
The riskless rate r > 0 is constant, and an EMM Q is chosen so that under
Q, {Xt} = {InSt} is an RLPE of order v € (0,2] and exponential type
[A_, A+], where A_ < — 1 < 0 < A + ; V denotes the characteristic exponent
of X under the risk-neutral measure chosen by the market.
Set x = ln(S/H), u(x,t) = Vd(H,T;S,t). Then for t < T and x € R,
To solve the problem Eq. (7.5)-Eq. (7.7), make the Fourier transform w.r.t.
T; since the terminal condition Eq. (7.7) is homogeneous, we obtain the
following family of the problems on R, parameterised by A with 3A < 0:
(In Chapters 3 and 5, we used Eq. (7.10) with c0 = 0; clearly, we can choose
-a- > 0 and o~+ > 0 so that Eq. (7.10) holds with c 0 > 0). It follows that if
0o > 0 is sufficiently small, then Eq. (7.10) holds uniformly in A in the half-
168 First-touch digitals
v(-,X) = ( M A ^ X i A r 1 ^ - . (7.11)
(Cf. Eq. (5.27)). If K_ < 1, the solution is unique in the class of bounded
functions, and if K_ = 1, it is unique in the class of bounded continuous
functions. Notice that the condition K_ = 1 holds (or fails) for all A simul-
taneously. The equality K_ = 1 implies that either the process is of order
2 and hence has a diffusion component, or it is a process of order u 6 (0,1)
and fi < 0, that is, a process of bounded variation, with negative drift. In
both cases, by using Eq. (7.1), we obtain that u is continuous. Explicitly,
Eq. (7.11) is
+oo+ia; +
e»^_(A,0(A0- 1 «,
/ -oo-f-iu;-)-
/ e^A+a:«)^_(A,0(AO-1^dA, (7.12)
-oo+ia J— oo+iu>+
C U + K D - ^ I + IAI)-1-',
for some C,e > 0. We conclude that v is continuous on (0, +oo) 2 , therefore
Eq. (7.12) gives the price we are looking for. By returning to the initial
variables, we obtain
Theorem 7.2 For S > H and t < T,
+oo+i<r /•+oo+iw+
/ (7-13)
/
-oo-\-i<7 J — oo+iu;+
.e*((T-*)A+ln(S/ff)C)^_(Ai^(A^rl^dAj
By making the inversion of the x-axis w.r.t. the origin, which leads to the
change X H. X, S i-» -S, H >-> -H, £ H-> - £ , and (f>-{\,0 >-» <MA, - 0 .
The Wiener-Hopf factorization with a parameter 169
and then making the change of variable £ i-» — £ in the integral, we obtain
the following theorem.
The reader may notice that Eq. (7.13) and Eq. (7.14) are not computa-
tionally effective, though relatively short. Below, we derive fairly effective
formulas; some of them are long, however.
The reader interested in the main formulas but not in the technical details
of the proofs can skip the next section, where we study the analytic con-
tinuation of the factors in the Wiener-Hopf factorization formula w.r.t. an
additional parameter A, and obtain necessary estimates.
E(e,0) = t e + E(0),
where
with e, 9 > 0, and the best situation is when a choice 6 > n/2 is possible.
Moreover, we need to have certain uniform estimates for the factors not
only w.r.t. £ in appropriate half-planes, but w.r.t. A g £(e,0) as well.
In such situations, the standard terminology in the theory of PDO is "the
factorization with a parameter17.
170 First-touch digitals
For the convenience of the reader, we rewrite Eq. (3.71) and Eq. (3.57)-
Eq. (3.60):
r + i\
= <MA,£)MA,0; (7.15)
r + i\ + il>(£)
the formulas for $ + (A,£), Q£> a-, are
+oo+ip-
exp -F
^nij-oo+ip. »A + r + V(»?)
In^Zidr,
t]
(7.17)
where p- € (cr_, 9f), and the formulas for <^>_(A,f), 3£ < a+, are
+oo+ip+
= exp In v-t dr\ ,(7-19)
-oo+ip+ i\ + r + ip(r)) r)
where p+ € (9^, cr+). Recall that in Chapter 3, the formulas have been
proved for i\ + r > 0. Clearly, each of the expressions in Eq. (7.15) and
Eq. (7.16)-Eq. (7.19) admits the analytic continuation w.r.t. A into some
domain, but to study the estimates for the factors in the Wiener-Hopf
factorization formula, the second representation of the factors is more ap-
propriate.
Eq. (3.71) assumes the form
o(A,0 = d a + ( A , 0 a - ( A , 0 , (7.20)
In the second case, a good estimate for the factors a± can be obtained on
A € £(e,0) with 9 < 7r/2, which is not sufficient for certain purposes. The
choice 9 > n/2 is still possible but only when a±(A, £) are constructed as
functions which are not holomorphic w.r.t. A (see Chapter 8).
Let c,n,u be as in Eq. (3.45)-Eq. (3.46). If v € (1,2] or v G (0,1)
and fi = 0, set d = c, K± — u/2; if v = 1, set d = (p? + c 2 ) 1 / 2 , n± =
0.5 ± 7T_1 arctan(/j/c). Thus, d and K± are the same as in Chapters 3 and
5. For s G R, set
and fix 9 G (0,7ri//2). If u > 1, we can (and will) choose 9 > ir/2.
If p- < 0, and e > 0, are sufficiently small, then A+(A, £) s is well-defined
in the region A € £(e,0), 9£ > p _ , and satisfies there an estimate
B(A (7 24)
'^-dA+(A,0«+A_(A,0«-- '
Due to Eq. (7.10), b(A,£) := ln.B(A,f) is well-defined for £ in a sufficiently
narrow strip 3£ € [p-,p+], where p_ < 0 < p+, and A € S(e,9), provided
e > 0 is sufficiently small. For £ in the half-plane 3 ^ > p_ and A e S(e, 9),
set
+
MA , f) i r "-^,, (7.25)
2?ri J-00+ip_ J] - £
(
*-™-$$- ™>
7.3.2 Analytic continuation of the factors w.r.t. A
If v > 1, fix 9 G (7r/2,7ri//2); if v G (0,1], fix 9 G (0,TTI//2), and assume
that Eq. (7.21) holds.
Lemma 7.1 Let X be an RLPE. Then there exist C,C\ > 0, p- < 0 < p+
and e > 0 such that for all £ in the strip 3£ € [p_,p+] and all A G £(e,0)
t/ie following three estimates hold:
and
Proof. It follows from Eq. (3.45)-Eq. (3.46) and Eq. (7.21), that as A ->
oo, and £ -¥ oo in the strip 9£ G [A_, A+],
\i\ + c | £ | f = |3?A|2 + | - 3A + c | £ | f
From now on, we assume that the conditions of Lemma 7.1 hold, and use
e, p~, p+ constructed in this lemma. We also fix arbitrary w_ G {p~, 0) and
uj+e{0,p+).
L e m m a 7.2 a) There exist C,C\,c\ > 0 such that for all £ in the half-
plane Sf > w_ and aZi A € Y,(e,8),
|fc+(A,OI<C, (7.34)
and
ftj There exists C,C\,Ci > 0 suc/i i/iai for all f in i/ie half-plane Sf < w+
and all A G £(e,0),
l&-(A,OI<C (7-36)
and
Proof. We prove Eq. (7.34)-Eq. (7.35); Eq. (7.36)-Eq. (7.37) are proved
similarly. It suffices to prove Eq. (7.34); Eq. (7.35) is an evident corollary.
By using Eq. (3.45)-Eq. (3.47) and Eq. (7.33), we easily obtain the
following estimates
(IAI + 1) 1 /", and for each pair (A, £), introduce intervals Jj C R:
Ji = {fl\\ri-S\<K},
h = {ri\\n-Z\>K, \v\<K},
J3 = fo||»7-fl>M, \V\>K},
J4 - {fl\K<\r,-i\<\rt\, \r,\>K}.
b(X,V) b(X,0
(7.40)
£ + ia — T] £ + ia — r\
where
l^(A^,r7,a)|<C3(l + | A | ) ^ - ^ . (7.41)
Since
f
J-K
drj
-AT ™ - 7J
In-
-K-ia
K-ia
<2?r,
b{X,j])dr]
\[ *{X\v)dTI <2nlnC + C3 [ (1 + \X\)^'-^"df, = C4. (7.42)
b(\,T))dT)
\[P±^<C5[ (l + |A|)-1/"c67 = C6. (7.43)
Further, we infer from Eq. (7.38) that b admits an estimate of the same
form as B — 1, and using this estimate on J3, we obtain
b(X,r])dT]
/ £ + ia - 7]
<C- f ~
1+W+ W dr). (7.44)
JJ3 J\n\>K \v\0-
By changing the variable 77 = Kr)', we see that the RHS in Eq. (7.44) is
bounded uniformly in A £ E(e, 0), £ £ R, a > 0. Since v' £ [0, v), a function
is decreasing on [0, +00), and therefore, we deduce from Eq. (7.38) an esti-
mate, for 7] G J4,
and the change of variable 77 = £ + Jf T/ shows that the RHS in Eq. (7.46)
is bounded uniformly in A G £(e, 0).
By gathering bounds Eq. (7.42)-Eq. (7.44) and Eq. (7.46), we obtain
Eq. (7.34). •
Lemma 7.3 a) For any N, there exist C, p > 0 such that if 9£ > w_,
A G £(e,0), and |A| < iVln(|f| + 1), then
6j For any N, there exist C, p > 0 suc/i tftai i/ 3£ < LJ+, A G £(e,#), and
|A| < JVln(|f| + l ) , i/ien
and
Proof, a) and b) are proved similarly, and Eq. (7.48) is immediate from
Eq. (7.47). To prove the latter, take e G (0,1), and consider the integrand in
Eq. (7.25) in the regions |ri| < (f) e , and \n\ > (£)£, where {£) = (l + l ^ 2 ) 1 / 2 .
For |r?| < (£) e , we have
K- r/T1 < c ^ r 1
and for any a > 0,
i\ + r + il>{n)
In < CAv)°
dA+{\,r])K+A_(\,n)K-
176 First-touch digitals
Hence, for a positive c-i and 9£ > c-i + p-, we obtain an estimate
r+°°+*p- b(\,T))
w
L -oo+ip- £ — *?
M<(Oe^drl <c(0~ 1+€(1+<7) -
If \v\ > (0% e use Eq. (3.45)-Eq. (3.46) to obtain that
1+ U l W
dA+(A,n)«+A_(A,r?)«- ''
where pi = v — v\ > 0, therefore
b(\,n) = 0((r,)-").
Now, to finish the proof of Eq. (7.47), it remains to notice that
+oo dn
< c(0~pl/2
/ -oo WHS-V) n
Lemma 7.4 a) There exist C, p > 0 such that for all £ in the half-plane
3f >w_ and all X 6 £(e,0),
\dxb+(X,0\< 0(1 + 1X1^ + ^)-", (7.51)
and
|9 A (a + (A,a ±1 )l < C(l + IAI1/" + | ^ | ) ± K + - p ; (7.52)
fc) There exists C, p > 0 suc/i t/iat /or aH £ in ihe half-plane S£ < w+ and
aHAeE(e,0),
From Eq. (7.55), we see that to finish the proof, it suffices to apply the
estimate
+°° dr]
I 2
(\)Kt WHZ-Tl)
<C(A>-pi/2(£)~pl/4
D
Theorem 7.4 Let v > 1. Then there exists 8 > 0 such that as x -> +0,
a) The asymptotics Eq. (7.57) holds not only for processes of order u > 1
but for processes of order 1, and for processes of order v G (0,1), if /x — 0,
as well; only the factors in the Wiener-Hopf factorization formula have
to be defined differently, as in Chapter 8. The proof remains the same.
b) The asymptotics Eq. (7.57) is uniform on each segment [TO, +OO) of the
positive half-axis in the sense that for each To > 0, the constant in the
O-term can be chosen the same for all T >TQ.
c) If there is a Gaussian component, then K_ = 1, and from Eq. (7.57)
we conclude that the solution is of the class C1 up to the boundary. If
there is no Gaussian component and Eq. (7.21) holds, then «;_ 6 (0,1),
and the price of the option behaves not as in the Gaussian case.
d) Still, the relative error of the Gaussian approximation is not large here
since the leading term is 1 in the both models. Similar effect is more
pronounced for barrier options without a rebate and junk bonds, since
in these cases, the ratio of the RLPE-price to the Gaussian price may
be very large indeed.
W(X,T) := (1-V(X,T))/XK-
is approximately constant near the barrier, and Fig. 7.3 demonstrates this
effect.
Price near the barrier 179
Fig. 7.1 The NIG-price of the down-and-out option, at T = 3,5,10 and 15 days to
expiry. Parameters: r — 0.05,5 = 7.2, a = 40, /3 — 8.
v = 1 + /,
+00+MT p + CO+iul-
/ e i < r * + *«0_(A,O(AO- x d£dA, (7-58)
/
-oo+itr J — oo+iu>-
0.16
• x=3
x T=5
- - T=10
x=15
xxx*xx**x><
Fig. 7.2 The Gaussian model vs. NIG model: the price of the down-and-out option, at
T = 3,5,10 and 15 days to expiry. Parameters: r = 0.05, <S = 7.2, a — 40,0 = 8.
+ QO+IU1-
f(x,T) = (2n)-2 jj_ ri(rx+xe) a-(A,0) d£d\. (7.59)
'CJ-oo+iw- a~(A,f)A£
(Here we have used Eq. (7.30)). Notice that there exist C, c0 > 0 such that
r+oo+iu)-
f(x,r) = (27T)'2 / + 0 O + " " f e**+*> a
~n'l!l*W. (7.60)
-oo+iw-
Price near the barrier 181
x x x x x x x x x x x
KXXXX X X X X X x x x x x x x
• t=3
X T=5
- - t=10
T=15
10
In (S/H)
Fig. 7.3 Graph of W(x,r) near the barrier. Parameters: r — 0.05,5 = 7.2, a = 40,
8, at T = 3,5,10 and 15 days to expiry.
In addition, there exist N,C such that if |A| > JVln(|£| + 1), then
JXr O-(A,0)
<Cexp[-r|A|co/2]($)- K — 2 , V^<w_,
JXr O-(A,0)
{l + iOK-+1X'
From Eq. (7.50), on the set |A| < iVln(|£| + 1), 9£ < w_ we have the
estimate
for some p > 0, C2 > 0. By taking the last three estimates together, we
182 First-touch digitals
a
I ^
Jc
n)vtxdX
o-(A,f)*f A
= ^l + i^~K"1 I eiXT^^-dX
Jc A
+ «&(0, (7-61)
with some p > 0. It follows from Eq. (7.62), that for any e > 0, w belongs
to the Sobolev space HK-+1/2+p~e(R) and vanishes on (—oo, 0]. But for
any s > 1/2 and e 6 (0, s - 1/2), HS(R) is embedded into the Holder space
C*-i/2-e(H,) j therefore w € C K - + " _ € ( R ) , V e > 0. Since w vanishes on
(—oo,0], we conclude that
Further,
+oo-|-iu»_
eixi{l + i O - " - - 1 ^ = T(K_ + l j ^ i ^ - e - * , (7.64)
/ -oo-j-iu'-
and by combining Eq. (7.61), Eq. (7.63) and Eq. (7.64), we arrive at
The second term in Eq. (7.66) exponentially vanishes as r —> +oo, and we
obtain from Eq. (7.57):
Asymptotics as r —y +00 183
Theorem 7.5 There exist 5 > 0 and o+ > 0 suc/i £/ia£ as x -» +0 and
r -> +00,
U(X,T) =
* " rfc^+V"" + ^ " " ^ + °( e ~' + r )- (7-67)
The result of this section is intuitively clear: when the terminal date is far
away, the price is essentially independent of time, the generalized Black-
Scholes equation reduces to the stationary one, and therefore, the price
is approximately equal to the solution of the stationary boundary value
problem on the half-axis. The formal proof is as follows.
In Eq. (7.58), transform the line of the integration 9 A — a into the
contour C, next change the order of the integration as in the beginning of
the proof of Theorem 7.4, and then transform C into C\. Then, similarly
to Theorem 7.5, we obtain
Theorem 7.6 There exists a+ > 0 such as r —* +00, uniformly in
x € [0, + 0 0 ) ,
+oo+iu>+
e ^ < / » _ ( 0 , 0 ( - ^ ) - 1 ^ + O( e - f f + r ), (7.68)
/ -oo+iw+
where w+ > 0 is sufficiently small.
Note that in the proofs of Theorem 7.5 and Theorem 7.6, an additional
assumption on v is used but they can be proved for any 1/, in the same way
as the pricing formula for the perpetual American options was derived in
Chapter 5.
This page is intentionally left blank
Chapter 8
Barrier options
(1) the down-and-out call option with the terminal payoff (ST — K)+\
we denote by V&0tCa\\(K, H, T; S, t) the price of this contract at the
185
186 Barrier options
An "in" option becomes the European option when the specified barrier
is crossed (before or on the terminal date); otherwise it expires worth-
less. For instance, the up-and-in call option becomes the European call
option when the barrier is crossed from below; we denote its price by
Vui;ca.u(K, H,T;S,t). Similarly, one defines the up-and-in put option, and
down-and-in puts and calls. In all states of nature, a European option pays
the same amount as the portfolio of the down-and-out option and up-and-
in option with the same barrier (of course, this is pertinent to a portfolio
of the up-and-out option and down-and-in option with the same barrier
as well), therefore the standard no-arbitrage consideration shows that the
pricing problem for an "in" option reduces to the pricing problem of the
corresponding "out" option and a European option, for instance,
where Vca\\(K, T; S, t) is the price of the European call with the same strike
K and expiry date T. The equalities similar to Eq. (8.1) hold for other
pairs of barrier options, therefore it suffices to write the explicit formulas
for "out" options only.
the general theory of Chapter 15. We illustrate the usage of the general
theory by explicitly solving the pricing problem for the European down-
and-out call V<i0;ca,ii(K,H,T;S,t) in the case K < H. Similarly, other
barrier options can be considered. For the pricing formulas of other barrier
options, more general payoffs and rebates including, see Boyarchenko and
Levendorskii (2002b), where simpler and more direct approach is used. No-
tice, however, that the general theory of Chapter 15 admit a straightforward
generalization for the case of option pricing under Levy-like Feller processes,
when only approximate analytical formulas can be obtained (see the end
of Chapter 16), whereas the approach in Boyarchenko and Levendorskii
(2002b) relies on the explicit analytical representation of the solution, and
therefore, it is more difficult to generalize the latter approach.
v(x,t) = E^[e^T^g(XT)lT0>T
where To is the hitting time of (—oo, 0] by the process X. According to the
general scheme described in Section 2.4, we approximate g by a sequence
of non-negative functions gn with compact support, which converges to g
pointwise: <7«(x) T g{x) f° r an Y x >Q. The standard construction of such
a sequence is as follows: take a non-increasing function \ S C°°(R) such
that x(x) = 1 for x < 1, x(x) = 0 for x > 2, and 0 < x(x) < 1 for x S [1,2],
and set gn{x) = g{x)x{x/n). Set
vn{x,t)]v{x,t). (8.2)
o*'
Here the initial data gn eH (R-+), for any s' < 1/2, since gn is smooth
up to the boundary x = 0, and has the compact support, and the unknown
wn is bounded and measurable. Take 7 G (A_,—1), and define wnn(x) —
eyxwn(x). Similarly, define 5,1,7,57, and w7. The choice of 7 is explained
as follows: <77(x) must decay as x —• +00 so that <77 S L2(R+). Thus, if we
consider a payoff with the different rate of growth at the infinity, another
choice of 7 may be needed; and for up-and-out put option, say, when the
problem on R_ is to be considered, and the payoff stabilizes to a constant
as x —> —00, one should choose 7 € (0, A + ).
Substitute wn(x) — e~yxwnn and gn(x) — eT",xgn~1 into Eq. (8.6)-
yx
Eq. (8.8), next multiply by e , and use the equality
the result is
Recall that we understand the problem Eq. (8.9)-Eq. (8.11) in the sense
of generalized functions, and we look for a solution which is bounded and
measurable on the strip [0, +oo) x [0,T]. The standard argument shows
that such a solution is unique, and we will explicitly construct it by using
the representation theorem for analytical semigroups.
Here / denotes the identity operator, and C is large enough so that on the
line 3£ = 7,
with some CQ > 0; the existence of such C and Co follows from Eq. (3.45)-
Eq. (3.46).
The inverse (Ay 4- CI)-1 is constructed as follows. Under condition
Eq. (8.14), the symbol a(£ + ij) + C admits a factorization
and we define
(A 7 + CI)~lv = o7,C;-(I')-1e+a7,c;+(£')-1^,
1
1) on the strength of Theorem 15.13, 0+a7tC;+{D) lv is independent
of the choice of Iv e Z<2(R), and the map
£ 2 ( R + ) 9 » H 6+a7^+(D)-Hv 6 L 2 (R+)
is bounded;
2) on the strength of Theorem 15.12, the map
L 2 (R+) 9 / H+ ^ . . ( Z ? ) - 1 / € L 2 (R+)
is bounded.
Hence, (A 7 + CI)-1 is a bounded operator in L 2 ( R + ) , and a direct
calculation shows that it is the inverse to
Ay + C J : 2 > ( A y ) - * L 2 ( R + ) .
(see Theorem 15.22). On the formal level, the system Eq. (8.12)-Eq. (8.13)
is easy to solve:
but in order to justify the answer, and to ensure the calculation of the
limit Eq. (8.2), some additional efforts are to be made. In order that
the representation theorem for analytical semigroups (Theorem 15.18) be
applicable, we have to impose an additional restriction on the order v and
fi in Eq. (3.45):
if i/ S (0,1), t h e n / i = 0. (8.16)
where £e,c is the boundary of Ee,c- Moreover, for any r > 0, the RHS
defines an operator, which is bounded in L2(R+). From our construction
of the sequence {gn,-y}n>i, it follows that it converges in Z,2(R+) to <?7;
hence, by passing to the limit in Eq. (8.18), and multiplying by e _ 7 X , we
obtain
When we construct the resolvent, we will be able to show that the RHS in
Eq. (8.19) is a continuous function; hence, the limit on the LHS is the price
of the down-and-out call option we have been looking for:
e- 7X a±(A, D + i 7 ) ~ V x
= a±(A, D)'1,
192 Barrier options
d+a+^Dy'g = 8+a+{\,D)-l{ex-K)
= ^+(a+(A,-i)-1e--JPfa+(A,0)-1),
and since the Fourier transform of the function on the RHS is equal to
o+(A)-i)-1(^-l)-1-/fo+(A>0)-1(tO-1,
we arrive at the pricing formula
We have proved
Theorem 8.1 The rational price of the down-and-out call option with
the barrier H = 1 and the strike K < H is given by
1 /• r+oo+ia
T t A+i
Vdo;cM(K,l,T;S,t) = ^- / e( - ) « l n S a _ ( A , 0 - 1 (8-24)
2m
JC»,C J-oo+ia
•[a+(A, - i ) ( * £ - I ) - 1 ~ ^ a + ( A , 0 ) - 1 ( i O " " 1 ] ^ A .
_1
Proof. From Eq. (3.45)-Eq. (3.46), we conclude that there exist R, c\ > 0
such that for all f in the strip $£ G [A_, A+], outside the ball of radius R,
centered at 0,
If C is sufficiently large then on U(C, c\), we have |£| > R, hence Eq. (8.26)
holds, and therefore, Eq. (8.25) holds as well, with some Co > 0. If C\ > 0
is sufficiently large, and |A| > Cj(l + |£|) y , then from Eq. (3.45)-Eq. (3.46),
we conclude that Eq. (8.25) holds, with possibly another CQ > 0. Finally,
let ci and C\ be chosen as above. Then on a set V(ci, C\, C), defined by
c l (i+Kir<iAi<ci(i+Kir, IAI>C,
we have, as C —> +oo:
and
A + a ( 0 ~ A - i / i £ + c|£|, if v = 1. (8.28)
These two conditions allow us to derive the pricing formula for any con-
tinuous payoff which grows as S13 as S —• +oo, for any /? 6 (—A+, —A_) J
Notice that due to the EMM-condition r + ip{-i) = 0, (ii) fails if C — 0
and p- = — 1. This observation highlights the importance of the flexibility
of the contour Ee,c- In particular, with the choice C = 0, it is impossible
to use the pricing formula for the down-and-out call.
With our choice of K±, function
and
w
|6_(A,0l<Ci, VAe£e,c ^ + - (8.35)
By the residue theorem, for A_ < a\ < a < ui < A+,
&(A, 77)
- /
a +oo+io-i /•+oo+ia'2\
£ + ia - 77
= b(\,£ + ia).
Hence, B± = expb± satisfy B = B+B- on E ^ c x {^ | S^ £ (w_,w + )},
and if we set
A + a ( 0 = a+(A,Oa-(A,0- (8-36)
Lemma 8.4
aj For any A e Se,c, a+(A, £) is holomorphic in the half-plane $>£ > A_,
and satisfies an estimate
c(l + IAI1/" + |e|) K - < |a_(A,0| < C{\ + IAI1/" + |£|)«-, (8.38)
where C,c > 0 are independent of A e £fl,c a n ^ £ ^ ^ e half-plane
S£ < w+ (7m£ may depend on w+ < A_J;
c) /or a// A € Se,c flnd £ in the strip A_ < S£ < A + , Eq. (8.36) holds;
d) factors in Eq. (8.36) are uniquely defined by properties a) and b), up to
scalar multiples, depending on A.
Proof. Clearly, A±(A, £) K± satisfy a) and b), and since b± are holomor-
phic and bounded on the same sets due to Eq. (8.34)-Eq. (8.35), a) and b)
are proved; Eq. (8.36) has already been proved.
To prove d), fix A, and suppose, A + a(£) = a+(A,£) a -(A,£) is an-
other factorization with the same properties. Then a'+(A, £ ) / a + (A, 0 (resp.,
a'_(A,£)/a_(A,£)) is holomorphic in the upper half-plane 3£ > A_ (resp.,
the lower half-plane S£ < A + ). Both functions are bounded and non-zero,
and coincide on R. Hence, the analytic continuation of any of them is a
bounded holomorphic function on C. By the Liouville theorem, it must be
constant. •
A; ) _ = (l + | A | r - a _ ( A , D + i 7 r 1 ,
and
Prom Eq. (8.37) and Eq. (8.38), we conclude that both operators are bounded
in £2(R+) (see Theorem 15.12 and Theorem 15.13), uniformly in A, and
hence, Eq. (8.17) has been proved.
To prove the continuity, we notice that gy G/f(R+), for any s < 0.5 since
g is the restriction to [0, +00) of a smooth function, which exponentially
decays at the infinity together with its derivative. Take 5 < 0.5 such that
s + K- > 0.5. By using Eq. (8.37) and Eq. (8.38) once again, we can show
Q 3 Q S-TK —
Consider the price of a barrier option without a rebate near the barrier.
It can be shown that if the time to the expiry is fixed then its behavior is
qualitatively the same as of W{X,T) := 1 — V(x,r), where V(X,T) is the
price of the corresponding first-touch digital. See Fig.7.1-7.3, which clearly
demonstrate how significant the difference from the Gaussian case can be.
8.4 Commentary
et al. (1995), Musiela and Rutkowski (1997), Hull (2000) and the bibliog-
raphy there, but to the best of our knowledge only Gaussian processes have
been allowed.
For the double barrier options with both upper and lower barriers (also
in the Gaussian case) see Geman and Yor (1996), Pelsser (2000), and the
bibliography there.
The results of this chapter, in a more general set-up, and by using a bit
different technique, are obtained in Boyarchenko and Levendorskii (2002b).
Chapter 9
Multi-asset contracts
The variety of derivative products on several assets is huge - see the ref-
erences at the end of the chapter, where fairly comprehensive overviews
for the Gaussian case can be found; in many cases, it is possible to obtain
RLPE-analogues of the results of the Gaussian theory.
We construct multi-dimensional analogues of RLPE, and calculate prices
of (power) forwards and several types of European options in the RLPE
market of several assets. We calculate weights of a locally risk-minimizing
portfolio of several underlying assets and a contingent claim of the Euro-
pean type, and the prices of the Asian options of the simplest type.
In the Gaussian case, the evolution of the tuple of several assets can be
modelled as
n
1=1
where bj, j = 1 , . . . , n, are constants, E := \<rij\ is a symmetric matrix, and
W(t) = (W\(t),..., Wn(t)) is an R n -valued standard Brownian motion pro-
cess. In other words, each component of a multi-dimensional Brownian mo-
tion is represented as a mixture of independent Brownian one-dimensional
motions. For a generic non-Gaussian Levy process, the similar reduction
to finite mixtures of one-dimensional processes is impossible.
The definition of multi-dimensional KoBoL uses the decomposition into
199
200 Multi-asset contracts
By repeating the proof of Lemma 3.1, which treats the ID-case, we conclude
that there exists 7 € R n such that the generating triplet (0,7, H(dx))
defines the characteristic exponent
Definition 9.1 Let v € (0,2), v ^ 1. The Levy process with the charac-
teristic exponent Eq. (9.1) is called a KoBoL process of order v.
Multi-dimensional Regular Levy Processes of Exponential type 201
The class of characteristic exponents Eq. (9.1) is too large for parameters
fitting purposes since the function A and the measure II' belong to infinite-
dimensional spaces. The following subclass of the KoBoL family depends
on the finite number of parameters, and still seems to be sufficiently rich
for practical applications. Another possibility is to use a measure U'(d(f>)
supported on a finite subset of the unit sphere.
We call a Levy process with the characteristic exponent Eq. (9.2) a KoBoL
process of order v, with the intensity parameter c, steepness parameter A,
asymmetry parameter /?, pseudo-covariance matrix S, and pseudo-drift /i.
dhp_
(0) = cT(2 -u)f (A + (/?, «A)) ,/ - 2 (E«A) j (E^) fc ^. (9.4)
If a KoBoL with the characteristic exponent Eq. (9.1) fits well to empirical
data, but one tries to fit the Gaussian process to the same data, the drift
and covariance matrix of the Gaussian process will be given by the RHS in
Eq. (9.3) and Eq. (9.4), respectively. This remark can be used when one
compares results in KoBoL and Gaussian cases.
Let Y(/3°) be the n-dimensional Brownian motion with the drift /3° 6 R™
and (non-degenerate) variance-covariance matrix A, and let X be the sub-
ordination of Y(/3°) by TS(u, 6,7). Then the characteristic exponent of X,
denote it (j>, is of the form
/3 = A-1lP, a 2 = 7 2 + (A/3,/3),
we obtain
Here (•,•) denotes the standard scalar product both in R n and in C n re-
garded as R 2 n with the basis {ei, • • •, en; ie\, • • •, ien}. With this identifica-
tion, A acts in R 2 n as a block-diagonal matrix, each of the diagonal blocks
being equal to the n x n matrix A.
By adding the drift, we obtain the characteristic exponent of a Normal
Tempered Stable (NTS) Levy process
following conditions:
(a) there exists a bounded open set U £ R n containing the origin such that
ip admits the analytic continuation into the tube domain Hn+iU C C n ,
and the continuous extension up to the boundary of the tube domain;
(b)
We consider forward contracts, European puts and calls, and exchange and
basket options. The market consists of n > 1 stocks with the price processes
Sj(t) = expX,(£), j = 1 , . . . ,n, and a riskless bond yielding the constant
rate of return r > 0. We assume that X = (Xi,..., Xn) is an RLPE under
the historic measure P .
204 Multi-asset contracts
dQ = e{K'x^-ctdP,
The EMM-condition means, in effect, that for the price of the riskless bond
and the price for any of the underlying stocks, we must have
By applying Eq. (9.12) and Eq. (9.13) to the riskless bond, whose price is
deterministic: f(Xt,t) = ert, and using the definition of the characteristic
exponent, we obtain the following equation
P
e-t(V (-i*)-c) _ 1
whence
c = iPP{-in); (9.14)
and by applying Eq. (9.12) and Eq. (9.13) to any of the underlying stocks,
that is, with f(x, t) = eXj, j = 1 , . . . , n, and taking Eq. (9.14) into account,
we derive
P
exj-t(V- ( - i e i - i K ) - i / ' P ( - i K ) + 7 - ) _ gXj j = I ... n.
and if a solution exists, we can define the Esscher transform of the historic
measure P by
European-style contracts 205
From Eq. (9.16), we derive the following formula for the characteristic ex-
ponent rfrQ of the process X under Q:
V>Q(£)=VP(£-«)-V>PHK)- (9-17)
Notice that at the initiation, the price of the forward contract must be 0,
and therefore, the delivery price is defined at time 0 by
= (2TT)-" / / ei<I'«>-T<r+*9<°'«"»<V=o5n(£n)de
+oo+icr
+oo-
ei*"t"-TW*«>*"»gn{£n)dSni
/ -oo+
or
r+oo+ia
-t-oo-t-ia
ei£nlnS.-r(r+^(0,Jn))5n($n)dSni (Q 2 2 )
/ -oo+to-
European-style contracts 207
In other words, the pricing formula is the same as in the case of one stock
S(t) = expX(t), where X is an ID Levy process with the characteristic
exponent ip^(0, •). We leave to the reader the verification of the fact that
1/^(0, •) is the characteristic exponent of a Levy process in ID.
The other types of European options on one asset or on weighted geo-
metric averages of several assets are priced in the same fashion.
Certainly, "far" and "near" do not have the exact meaning; presumably, "at
the distance more than 0.25" and "at the distance less than 0.25" can be
used, and if the results do not agree with the desired error at the boundary,
where "far" becomes "near", one can change the definitions of "far" and
"near".
In the case a), we set f(X(0),t) — 0.
In the case b), we expand G(S) in the Taylor series around S(0):
G(S)~£ C a (S(0))(S-S(0)) Q ,
a
208 Multi-asset contracts
G(S)~ £ ca(S(0))(S-S(0)r,
\a\<m
and after multiplying out, obtain an approximation for g(x) = G(S) of the
form
Now, by substituting the approximation Eq. (9.23) into Eq. (9.19), and
using Eq. (9.20), we obtain
Notice that the construction in the case b) presupposes that ZQ, where
max{<7(x), 0} loses smoothness, is sufficiently far so that the influence of Z0
can be neglected.
In the case c), we cannot neglect it, and the procedure becomes more
complicated. Let Xoo £ ZQ be the closest point to X(0). In a neighbour-
hood U of Xoo, which contains X(0) but is not too large (probably, for
typical parameters' values, the neighbourhood of radius 0.4-0.5 will do),
we introduce the new coordinate system as follows. Let x — x(x'l) be the
parameterisation of ZQ in the neighbourhood of Xoo, s.t. Xoo = ^(0), and
set x'2 = \n(g(x)+Ki), where K\ > 0 is sufficiently large so that g(x) + Ki
is positive and bounded away from zero in the chosen neighbourhood U.
Then in U, g{x) — ex* — K\, and the payoff becomes the payoff of a Eu-
ropean call: max{ex2 — Ki,0}. By making the change of variables in the
Cauchy problem for the generalized Black-Scholes equation
we obtain
where ^(x', Dxi) is the PDO ^ ( A r ) written in the new coordinate sys-
tem. The calculus of PDO allows one to write down the asymptotic expan-
Locally risk-minimizing hedging with a portfolio of several assets 209
sion of the symbol, and after that, obtain the approximate solution to the
problem Eq. (9.25)-Eq. (9.26). For details, see Section 16.5.
(F(S(t),t); e1(t)S1(t),...,0n(t)Sn(t))
The problem and construction of the hedging portfolio can easily be mod-
ified for the case when F is a basket containing not only derivatives but
some of the basic securities Sjt, I = 1 , . . . , m; in this case, the corresponding
dj, (t) are assumed to be zero. The modification for the case of the restricted
borrowing is also possible.
The locally risk-minimizing formulas for the weights of the hedging port-
folio, constructed in the paper, admit a natural modification for the case
of Levy-like Feller processes constructed by Barndorff-Nielsen and Leven-
dorskii (see Chapter 14). This modification can be applied to portfolios
containing bonds and stocks.
As in the one-asset-case in Chapter 4, we show that the weights of the
hedging portfolio for a European option are continuous up to the expiry
date even at the strike price, where they are discontinuous in the Gaussian
case.
210 Multi-asset contracts
and
(ii) F, Sj, FSj, and SjSk, j,k = l,...,n, are priced under P .
Let X be an RLPE both under the historic measure P and the EMM Q
chosen by the market, and let Up and U^ be the corresponding open sets
in Definition 9.2. Let {ej}J=1 be the standard basis in R n . We will use the
identity
and
for all j , k = 1,..., n, then conditions (i) and (ii) are satisfied. Notice that
Eq. (9.29) holds for European calls and puts.
Fix t, and denote Xj = \nSj(t), f(x,t) = F(eXl,... ,ex";t).
The price of the portfolio at time t + At is
n
W(t + At) = F{S{t + A),t + At) + ^BjflSjit + At) + erAtw0(t),
3= 1
where wo{t) is the residual in Eq. (9.27), and we assume that the investor
chooses 6(t) to minimize the variance of W (t + At) under the historic mea-
sure P, conditioned on the information available at the date t. We calculate
Ef[(W(t + At) - Ef[W(t + At)})2)
= Ef[(F(S(t + At),t + At) - Ef[F(S(t + At),t + At)})2]
n
+2 J2 Oj(t)E?[(F(S(t + At),t + At) - E?[F(S(t + At),t + At)})
3=1
x (£,-(* + At) - Ef[Sj(t + At)})}
Minimizing w.r.t. 6(t), we obtain the following first order conditions: for
j = l,...,n,
71
The explicit formula for the terms on the LHS is immediate from the defi-
nition of the characteristic exponent (and from Eq. (9.20)):
EP[Sj(t+At)}Ep[Sk(t+At)} = exp(-At^p(-iej)+iJp(-iek)))Sj(t)Sk(t).
(9.33)
By applying the Taylor formula to Eq. (9.32) and Eq. (9.33) around At = 0,
we see that the LHS admits the representation
n
^ ^ ( t ) 5 , ( i ) 5 f c ( t ) [ - V P H f e +ek)) + 1>p(-iej) + ^p(-iek)} At+o(At).
fc=i
(9.34)
Notice an evident formula for the infinitesimal generator L of the (n + 1)-
dimensional process {Xt} — {{Xt, t)}, when X is regarded as a Levy process
under P :
L = dt + Lp = dt-rpp(Dx). (9.35)
By using Eq. (9.35), we obtain for the terms in the RHS of Eq. (9.31)
and
= (dt-^p{Dx))f(x,t)At + o(At).
By inserting Eq. (9.34) and Eq. (9.39) into Eq. (9.31), dividing by At and
Sj(t) and passing to the limit At —> +0, we obtain, for j = 1 , . . . , n:
n
J^ek(t)Sk(t) [-4>p(-i(ej +ek)) + iPP(~iej) + VPHefc)] (9-40)
fc=i
COV:=COV(X):=[A(ei,ek)]lk=l,
If
COV ^ 0, (9.42)
The direct calculation shows that COV(X) = [djk] is the covariance matrix,
and
n
A(iD,ej) = y^gjfcdj.
k=i
Thus, in the Gaussian case, Eq. (9.41) reduces to the standard delta-hedging
result, which gives the riskless portfolio. In the non-Gaussian case, the
locally risk-minimizing formula Eq. (9.43) does not determine a riskless
portfolio, the latter being non-existent.
If Eq. (9.42) holds, and an explicit formula for F(S, t) = f(x, t) is avail-
able, Eq. (9.43) can be used to calculate w = w(t). Notice that Eq. (9.43)
gives locally risk-minimizing weights for any contingent claim, not neces-
sarily of the European type.
/ Q
... / e *<*.€>-'-('-+* (fl)$(£)(£1.. . # „ ,
-oo+i<7i J — oo+i<rn
(9.44)
where g(£) is the Fourier transform of g, and the choice of real dj, j =
1,2,..., n, is determined by the type of the contingent claim (for instance,
for the put on the last stock, one takes <rn > 0 such that —anen e C/Q, and
<7i = • • • = <7n_! = 0). Therefore in Eq. (9.43),
... (9.45)
-OO-HiTi J —00+i(Tn
exp [t(x, 0 - r(r + V>Q(0)] Aft, e ^ S t f K i • • • # « •
If the terminal payoff depends on the price of the last stock only: g(Xx) =
g n ( X r ( n ) ) , for instance, if F is the price of a European put on the last
stock, then g(£) = <?„(£„) ® %=o, where £' = ( £ i , . . . ,£n-i), and hence
Eq. (9.45) reduces to a much simpler expression
+ 0 0 + t<7„
(9.46)
/ •oo+icrn
• exp [ixnn - r(r + ipQ(0, T)en))] A(ir)en, ej)gn(jf)dr}.
Locally risk-minimizing hedging with a portfolio of several assets 215
i/-e<l, (9.47)
where e > 0 is from Eq. (9.9), let Eq. (9.30) holds, and let F(S(t),t) be the
price of the European put or call.
Then the weights of the locally risk-minimizing hedging portfolio given
by Eq. (9.43) and Eq. (9.46) are continuous up to the expiry date, for all
strikes.
Proof. For European puts and calls, the Fourier transform of the termi-
nal payoff decays as |J7| -2 , as 7] —> 00, and hence, under condition Eq. (9.47),
the integrand in Eq. (9.46) is bounded (in modulus) by an integrable func-
tion C(l + |T?|) - 2 + l / - e . Since the integrand is a continuous function in
(X,T) 6 R x [0,+00), we conclude that the integral is continuous on the
same set, too. ' •
As the proof shows, the result is valid for straddles, strangles and essentially
any other European contingent claim used in real Financial markets, with
the payoff, which is continuous, piecewise smooth, and satisfies (i) and (ii)
of the previous section (but it fails for digitals, of course).
Since the Gaussian delta-hedge produces discontinuous weights at the
strike and expiry, our result Eq. (9.43)-Eq. (9.46) in the non-Gaussian
RLPE-case differs drastically from the Gaussian hedging, near expiry and
strike. The difference can be substantial in other regions of (x, £)-space as
well, especially in the market, where some of the assets are highly correlated
stocks.
and taking into account that A(Q,ej) = 0, we derive from Eq. (9.43)
An Asian option is the name used to denote options whose terminal payoff
depend on the average asset values during some period of the option life
time. There are Asian options both of the European style and the American
style; the averaging can be in continuous time, e.g., with the payoff of the
form
where
Consider an option with the terminal payoff depending on both X(T) and
J(T): g(X(T),J(T)).
In addition to the observable vector of log-spot prices at time £, X(t),
we have the observable J(£), and clearly, the price of the option at time t
depends on the triple (X(t),J(t),t); denote this price f(X(t),J(t);t). At
a sampling date tm,
F ( J + x) = ^(sir'd'jF^x3 (9.53)
s>0
= exp[xdj}F(J),
we can move all the exp[(5s, x)ir)]- factors to the left, and obtain
where
of the Asian option with the payoff function g is given by Eq. (9.56)-
Eq. (9.58). If g is independent of x, then the simpler formulas Eq. (9.59)-
Eq. (9.60) hold.
This page is intentionally left blank
Chapter 10
221
222 Investment under uncertainty and capital accumulation
process, for example, the geometric Brownian motion. We model the price,
Pt, as Pt = exp(-Xt), where X = {Xt}t>o is an RLPE. We demonstrate that
RLPE are not less tractable in models of investment under uncertainty than
Gaussian processes, the former being more realistic than the latter.
At each moment of time, the firm's manager has to decide whether to
increase or not the stock of capital given the current stock and the market
price. If the firm suffered an adverse demand shock in the past, its capital
stock may be too big for the current price level, but since investment is
irreversible, the firm has to stay put. Should the price increase, it may
become optimal to invest, and in this case, the manager needs to determine
the optimal amount of new capital. We assume that it is optimal to invest
when the price crosses a certain barrier, which is called the investment
threshold. As soon as the investment threshold is determined, the optimal
amount of investment follows automatically, as we will see later.
There are at least two different prescriptions how to choose the trigger
price of investment. One of them is known for almost a century, and is
called the net present value (NPV) rule, or the Marshallian law. This rule
(in marginal terms) says that investment is optimal as long as the present
value of the expected marginal revenue does not exceed the marginal cost of
investment. The second rule was spelled out by Dixit and Pindyck (1996),
and it is based on the option pricing approach to the problems of investment
under uncertainty. Dixit and Pindyck (1996) pointed out that the NPV
rule was incorrect since it did not take into consideration the option-like
nature of investment decisions. To correct the Marshallian prescription,
the marginal cost of investment has to be multiplied by a certain factor.
This correction factor comes from a solution to a so-called fundamental
quadratic equation and has no economic meaning.
We use the Real Options approach as well and solve the problem essen-
tially by the method of Chapter 5, because the choice of the investment
threshold is similar to the choice of the optimal exercise price of a per-
petual American call option. Our approach not only allows deriving an
analytical solution for the case of output price following an RLPE, but also
yields several new insights. The most operational one from the standpoint
of economic applications is a modified version of the Marshallian law.
We consider only such processes that the price can move both up-
wards and downwards with positive probability. If the price process is
non-decreasing almost surely, the irreversibility does not matter, and the
Marshallian law is correct. On the other hand, if the price process is non-
The investment threshold 223
increasing almost surely in t, then the firm will never increase the capital
stock further, because the state of the market can only deteriorate. When
the price can move in both directions with non-zero probability, starting
with the original price process, it is possible to introduce a non-trivial in-
fimum price process. N = {Nt = info<s<t Pt}t>o- The solution to the
problem under consideration allows us to resurrect the Marshallian law to
a certain extent. We state that it is optimal to invest as long as the present
value of expected marginal revenue, computed for the infimum price pro-
cess instead of the original price process, does not exceed the marginal
cost of investment. Notice that this investment rule reflects the "bad news
principle", which was first formulated by Bernanke (1983). The critical
price which triggers new investment depends on downward moves, because
the ability to avoid the consequences of "bad news" leads us to postpone
investment decisions. The main advantage of our investment rule is that
it obviates the need to introduce the correction factor like in Dixit and
Pindyck (1996). In other words, the firm's manager may remain in the
Marshallian world as long as she keeps in mind infimum processes instead
of real ones.
Here we specify the objective of the firm and sketch the solution technique.
Details can be found in Boyarchenko (2001).
The firm's objective is to maximize the total expected discounted profit:
r+oo
r-roo
IL(K,x) = max Ex / e-rt(eXtG(Kt)-rKt)dt (10.1)
{Kt} Jo
where the meaning of the stochastic integral above is explained in Bo-
yarchenko (2001). Here we treat the spot log-price, x, and the level of
accumulated capital, K, as state variables, and Kt as a control variable.
Due to irreversibility of investment, in each state (K,x), Kt > K.
In order the stochastic integral in Eq. (10.1) were bounded, we impose
the following conditions:
Following the tradition in the literature, we are going to view the space
224 Investment under uncertainty and capital accumulation
The next important problem, which can be solved in the framework of our
model is the accumulation of capital in the long run. As it was stressed by
Hubbard (1994), the benchmark models of investment under uncertainty
considered in Dixit and Pindyck (1996), do not suggest specific predictions
about the level of investment. Since the investment rule itself is not observ-
able, one has to use the data on investment and capital stock to evaluate
investment models; this shows the importance of theoretical results on the
long-run capital accumulation. Abel and Eberly (1999) examine the be-
haviour of the capital stock in the long run and calculate explicitly the
impacts of irreversibility and uncertainty on the expected long-run capital
stock. They assume that an exogenous demand shock follows a geometric
Brownian motion. Two types of effect are revealed: the user cost effect,
which tends to reduce the capital stock, and the hangover effect, which
226 Investment under uncertainty and capital accumulation
arises because the irreversibility prevents a firm from selling capital when
its marginal revenue product is low. Neither of these effects dominates
globally, so the effect of increased uncertainty cannot be determined unam-
biguously.
We address a similar problem in the case when the output price follows
an RLPE. To this end, we specify the production function as a Cobb-
Douglas one, i.e., G{K) — dKe, 6 6 (0,1); the marginal cost is normalized
to unity. We consider a new born firm with the price process P starting at
PQ = 1. As before, our reasoning is valid when prices move both upwards
and downwards. We introduce a (non-trivial) supremum price process as
M = {MT = sup 0 < s < t P t }t>o and derive the following formula for the
expected value -£[1^] of the capital, Kt, at time t:
where W_ and W+(£) are determined by the infimum process and supre-
mum one, respectively:
/•+oo
W_ = E / e-rtNtdt | No = 1 ,
Jo
W+(t) = E{MlK1~6) | Mo = 1].
E[KtR] = {edW^/^-^W^t),
where
r r+°o
/•-t-OO
W* = E / e-rtPtdt | Po = 1 ,
Jo
W*(t) = E[Pt1/{1~9) | Po = 1] = ea\
by considering the ratio of the expected levels of accumulated capital
K(t) = § ^ = UCHE(t),
where
uc = (W-/w*)1/(1-v
is the user-cost effect, and
mt) = E[M!Z~lIM° = 1]
- the hangover effect of the irreversibility, and study how both these effects
and the ratio n(t) itself depend on parameters of the process. Clearly,
UC < 1 (unless the process is non-decreasing), and HE(t) > 1 (unless
the process is non-increasing), so that both effects work in the opposite
directions, and the joint effect on the capital accumulation is ambiguous:
n(t) can be larger or smaller than 1, depending on parameters of the process.
Dixit and Pindyck (1996) show that the Marshallian law prescribes very
low trigger price of investment. A natural question which arises in our
model is how the investment threshold changes if one considers a regular
Levy process instead of a Gaussian process with the same first and second
228 Investment under uncertainty and capital accumulation
not differ too much from a Gaussian one, and the decrease of A + (respec-
tively, —A_) can be interpreted as the taking into account the more and
more of negative (respectively, positive) large jumps. We believe that the
result of this computational experiment is instructive. We allow A_ to vary
from —30 to —22, and A+ from —A_ to — A_—8, and we see that though the
variance does not change much, the threshold changes four-fold. The other
observation is that the threshold is much more sensitive to the increase in
skewness than in variance and kurtosis. Consider two cases: 1) A_ = —30
remains fixed and A + decreases from 30 to 22, and 2) both A+ and —A_
decrease from 30 to 22. In the first case, the variance and kurtosis increase
less than in the second case, and the threshold and modulus of skewness in-
creased more in the first case. So even this insignificant increase of variance:
by less than 7 % (and more significant increase of skewness and kurtosis),
which can easily be disregarded in practice since they come mostly from the
tails of the probability density, and hence, can be left over as too extreme
events, can have a dramatic impact on the investment threshold.
The results of the numerical example highlight the importance of ob-
taining the correct estimates for the parameters of an RLPE governing the
underlying stochastic variable. In principle, it is possible to infer these
parameters from the moments (up to the fourth order) of the correspond-
ing probability distribution. However, this will provide a good fit for the
central part of the distribution only, but not for the tails of the distribu-
tion. To obtain parameter estimates fitting the whole distribution, more
sophisticated methods should be used.
Now we consider the errors which may arise when one uses Gaussian
models in non-Gaussian situations, both in non-Gaussian and Gaussian
models. Consider the same firm and process as in the example for the
investment threshold above, with r = 0.06, 8 = 0.33 and d = r/6 in the
Cobb-Douglas function, c = 0.5, v = 1.1. We change A+ in a wider range
from 2 to 30; A_ changes in the range from -30 to -10. So, here we allow for
much fatter left tail. Following Abel and Eberly (1999), we determine the
last parameter, fi, from the requirement that the expected rate of growth
of the price remains constant. For each set of parameters, which define
a non-Gaussian KoBoL process Xt — In Pt, we take the Gaussian process
with the same variance as that of Xt.
The numerical results show (see Boyarchenko (2001) for details) that the
investment threshold in the non-Gaussian case decreases up to 15 percent,
230 Investment under uncertainty and capital accumulation
Let —A_ and A_ be large, which means that the rate of the exponential
decay of the tails of the probability density is large. Then for typical param-
eters' values for model classes of RLPE (Normal Inverse Gaussian processes,
Hyperbolic processes and KoBoL family), the equation r + if>(£) = 0 has
exactly two solutions; both of them are purely imaginary: —i[3±, where
/?_ < 0 < /3+. In this respect, model classes of RLPE are similar to
Brownian motion, but their characteristic exponents are not holomorphic
in the whole complex plane. Further, the approximate formulas Eq. (3.82)
are valid, which look like the exact formulas in the Brownian motion case.
It can be shown that for typical parameters' values, the relative error of
Eq. (3.82) is several percent or less. For details, see Chapter 5, where the
next correction term is given. Hence, these formulas can be used for com-
parative statics purposes, and in this approximate form, RLPE-modelling
is not more difficult than in the Brownian motion case - and more realistic
at the same time.
Chapter 11
11.1 A n overview
Two approaches to credit risk modelling can mainly be divided into two
categories: intensity based models, where the primitive is a stochastic model
for the default event (see the monograph Bielecki and Rutkowski (2001),
and in the non-Gaussian setting, Eberlein et al. (2000a, 2000b), Eberlein
and Ozkan (2001), and the bibliography there), and the structural approach.
In this Chapter, we apply the structural approach to credit risk modelling
for the case of non-Gaussian processes. The central feature of this approach
231
232 Endogenous default and pricing of the corporate debt
is an attempt to model explicitly the evolution of the assets of the firm. The
structural approach was pioneered by Black and Scholes (1973) and Merton
(1974) and extended by Black and Cox (1976), Longstaff and Schwartz
(1995). In a series of papers, Leland (1994a, 1994b) and Leland and Toft
(1996) focus on the choice of optimal capital structure and bankruptcy level
by a firm which keeps a constant profile of debt.
All the aforementioned structural models of credit risk have relied on
diffusion process to model the evolution of the firm's assets. While the
diffusion approach is convenient, it cannot capture the basic features of
credit risk observed in practice. In particular, Jones, Mason and Rosenfeld
(1984) find that the credit spreads on corporate bonds are too high to be
matched by diffusion approach. Also, under a diffusion process, firms never
default by surprise because a sudden drop in the firm's value is impossible.
Therefore, if a firm is not currently in financial distress, the probability
of its default on very short-term debt is zero. This fact implies that the
credit spreads tend to zero as the maturity of debt tends to zero, which
contradicts empirical evidence (see, for example, Fons (1994) and Sarig
and Warga (1989)).
The first thing that must be said about realistic models of the assets
of the firm is that jumps of essentially any size have to be incorporated.
Moreover, both negative and positive jumps have to be admitted to cap-
ture the evolution of assets in the information economy. One of the ways to
divert from the diffusion approach is to model the underlying asset process
as a jump-diffusion process. Zhou (1997) uses a jump-diffusion model in
the framework of the structural approach to valuing default-risky securi-
ties. Under a jump-diffusion approach, a firm can default instantaneously
because of a sudden drop in its value. The model allows for matching the
size of credit spreads on corporate bond to empirical data. However, the
bankruptcy level is given exogenously in Zhou's model, and the author does
not study the optimal capital structure. It is also necessary to mention that
in order that jump-diffusion models could be fitted well to empirical data,
many jumps have to be incorporated, which makes the models not very
efficient from both the analytical and computational points of view.
Here we study the optimal capital structure and endogenous default by
modelling the evolution of the firm's assets as a Levy process. One of the
advantages of models based on Levy processes is that such models admit
any number of jumps (infinite number including). Realistic models can be
obtained within a family of regular Levy processes of exponential type. Our
Endogenous default 233
goal is to show that RLPE are analytically tractable for problems of credit
risk modelling and endogenous default. Roughly speaking, we are going to
demonstrate that RLPE are the "second-best" if the Brownian motion fails
to produce realistic results.
We solve the endogenous default problem by using essentially the same
approach as in Chapters 5 and 10. We take any candidate for the level
of bankruptcy and calculate the corresponding candidate for the value of
equity by applying the Wiener-Hopf method. After this, we show that there
is only one candidate for the level of bankruptcy for which the found equity
satisfies certain natural conditions. The optimality can be proved as in
Chapter 5.
The plan of the Chapter is as follows. In Section 11.2, we specify and
solve a simple endogenous default problem. In Sections 11.3 and 11.4, we
study the firm near the bankruptcy, and a solvent firm, respectively, and
suggest two types of approximate formulas for all endogenous variables. The
first set of approximate formulas can be used very close to the default level,
and the second set of approximate formulas is valid for Levy processes with
large truncation parameters, i.e., for the processes with not very fat tails,
and the firm not very close to the bankruptcy level. These approximations
do not neglect the effects of non-Gaussianity on one hand, and on the other
hand, they allow to perform comparative statics analysis. In Section 11.5,
using the approximate solutions, we derive the optimal leverage for the firm.
We produce several series of numerical examples which show that for the
firm not very far from the default level, the difference between the Gaussian
and non-Gaussian models can be sizable, and very near the bankruptcy level
- very large indeed. In Section 11.6 we discuss the results obtained in the
Chapter and address the comparative statics issue. Technical results are
proved in Section 11.7.
Vt=expXt, (11.1)
erature (see, for example, Leland (1994a)) we assume that a riskless asset
exists and pays a constant interest rate r > 0.
The firm issues debt, which attracts coupon payment Cdt when the firm
is solvent. In the event of default at time t, a fraction a G (0,1) of the
firm's value is lost to bankruptcy costs, therefore the debt holders receive
(1 - a)Vt = (1 - a) expX t . Let VB = eh be the bankruptcy level, i.e., the
firm defaults should the value of its assets fall below VB • Let r(h) be the
hitting time of (—oo,h] by X. For D(Xt), the value of the debt, we have
rr{h)
D(x) = Ex 'Cdt + Ex \e-rT(h\\ - a)e x '< h >] . (11.2)
Eq. (11.5) is the analogue of (8i) in Leland (1994a); the condition (8ii) in
that paper is valid here as well. In fact, it suffices to require that BC is
continuous and bounded; such a solution is unique.
Now we consider tax benefits, TB(x), associated with debt financing.
These benefits resemble a security that pays a constant coupon TC1 as
long as the firm is solvent, and pays nothing at the bankruptcy. Hence,
1
T here should not be confused with the stopping time r(/i) above; the latter will not
appear any more
Endogenous default 235
The value of equity is the total value of the firm less the value of debt:
it can be represented as
E(x) = ex + / ( x ) , (11.11)
where f(x) = TB(x) - BC(x) - D(x). Notice that problems Eq. (11.3)-
Eq. (11.4), Eq. (11.5)-Eq. (11.6) and Eq. (11.7)-Eq. (11.8) are boundary
problems for a linear operator on the same half-axis, each having a unique
solution. This means that / is a unique bounded solution to the boundary
problem of the same sort, only with the data being the sum of the data of
the above problems, i.e.,
Now it is easy to see that f(x), D(x), TB(x),BC(x) are solutions to the
following boundary problem:
^ - ( D J l f . o o ^ f P ) - 1 ^ = A^e' - e f c ^ ( - 0 _ 1 * i ( x - h)],
where
*1(x)=<t>;(D)9+(x)ex, (11.14)
*o = (l>7(D)0+- (11.15)
We have proved
C(r-l)
*o(x -h) + eh(f>-{-i)-1yl{x - h)
simplifying, we obtain
where
w(y) = r-xC{T - 1) + e ^ - ( - i ) - 1 .
The value of the firm's equity having being obtained, the next step should
be to determine the level of bankruptcy for the case, when the firm can
choose this level, i.e., in the absence of a covenant (for example, positive
net-worth requirement). Leland (1994a) computes the endogenous level
=
of bankruptcy by using a smooth pasting condition: dE/dV\v=vB 0-
However, as we already saw in Chapter 5, the smooth pasting fails for
some non-Gaussian processes. Therefore, we are going to find the optimal
bankruptcy level as in Chapter 5, without resorting to the smooth pasting
principle.
Set
Notice that due to Eq. (11.20), d{h) > 0 iff w(x) > 0 V x > h. If d(h) < 0,
we deduce from Eq. (11.19)-Eq. (11.20) that E{x) is negative at some points
x > h, which is impossible. Hence, h> /i„, where /i* solves d(h) — 0, or
VB=eh' = C{l T)
~ <t>;{-i), (11-21)
r
and it remains to show that with any choice h > /i„, the value of the equity
Since w(y) > 0 on [/i,,+oo), we conclude that E{h*;x) > E(h;x), and the
optimality 2 has been verified.
Thus, we have proved the following theorem.
VB = C(l - T)E
f
Jo
a~r^ont nt
e dt (11.24)
where {nt} is the infimum process. It follows from Eq. (11.24), that the
asset value, VB, at which bankruptcy occurs is proportional to the coupon,
C, independent of bankruptcy costs, a, and decreases as the corporate tax
rate, r, increases. Moreover, we can see that the bankruptcy level does not
depend on the current asset value, V, but it depends on the parameters of
the infimum process for the firm's assets.
In the case of the Brownian motion, the bankruptcy level Eq. (11.21) is
exactly the same as determined by the smooth pasting condition in Leland
(1994a); it can be shown exactly as in Chapter 5, where we have shown that
our general optimal exercise price for the perpetual American put coincides
with Merton's result if X is the Brownian motion.
Now we determine the values of the firm's equity, debt, tax benefits,
bankruptcy costs and the value of the firm itself, given the value of the
assets and the level of bankruptcy, the former being a linear function of
the coupon C. By substituting Eq. (11.21) into Eq. (11.22), we obtain the
value of equity under the optimal choice of the default level
By using the explicit formulas for the factors <j)f, first, h» can be calcu-
lated from Eq. (11.21), then $ 0 ( z - M and $ i ( z - ft,) from Eq. (11.15)-
Eq. (11.14). Finally, the values in Eq. (11.25) and Eq. (11.27)-Eq. (11.33)
can be computed. However, since the formulas for the factors involve inte-
grals, it is difficult to use them for comparative statics purposes. Therefore,
we are going to derive approximate formulas for all the relevant values.
Consider the firm near the bankruptcy, i.e., the case of V/VB near 1, or
equivalently, x :— \II(V/VB) > 0 small. Clearly, this is the situation when
the presence of jumps must be felt in the first place, and hence the Gaussian
models can be expected to produce serious errors, and we show that the
mistake in the valuation of the firm equity near the bankruptcy level may
240 Endogenous default and pricing of the corporate debt
E = C (1-T),,(0W VA~»
rr(K_ + 2) \ VB)
Equity of a firm near bankruptcy level and the yield spread for junk bonds 241
Since the error of the approximation Eq. (11.40) is of the same order as the
value of equity in Eq. (11.39) or even larger, the usage of these formulas
for the calculation of the leverage makes no sense: the only result will be:
the leverage is approximately one, which is evident. (With the next terms
in Eq. (11.34)- Eq. (11.38), more information can be obtained).
For the return on the debt and the yield spread more informative for-
mulas are available:
C C r
R= = (U 41)
D~ (l-a)Vb (l-a)(l-T)#r(-»)' '
c r[i-(i-tt)(i-T)^(-o]
S r (1L42)
~D ~ (i_a)(i_r)0-H) •
They give the limit values of the return on the debt and the yield spread of
the firm near the bankruptcy level. In Fig. 11.1 we plot the default level VB
and the limit values of the yield spread in the KoBoL model, and the ratios
of these quantities in the non-Gaussian model to the ones in the Gaussian
model with the same volatility, for the interest rate r = 0.06, corporate tax
level r = 0.35, and bankruptcy cost a = 0.5; the characteristic exponent
is given by Eq. (11.35) with c_ = c + = c, u = 1.2, and the steepness
parameters A_ and A + varying from -30 to -8 and from 30 to 8, respectively.
The volatility vol=0.35 is fixed, which defines c as a function of A_ and A+,
242 Endogenous default and pricing of the corporate debt
from
r + ip(-i) = 0 (11.43)
We see that for the limit value of the yield spread, the difference between
the non-Gaussian case and the Gaussian one is much smaller, only several
percent.
If the firm is not near the bankruptcy level, i.e., h = ln(V/Ve) is not very
small, the approximate formulas above are not good or even not applicable
at all. Fortunately, there exists another set of approximate formulas, which
give good approximations, if h is not very small - and if it is large, these
approximations are very good indeed. Namely, as we noticed in Lemma 3.7
and used in Chapters 5 and 7, the models classes of RLPE: TLP, NIG, NTS
Levy and HP, have a very useful property, namely, a(£) = r + VK£) admits
the analytic continuation into the complex plane with cuts (—ioo,iA_] and
[i\+, +oo), has no zeroes outside the imaginary axis, and at most 1 zero in
(0, iA_|_) and at most one zero in (iA_,0).
Moreover, for parameter values typically observed in Financial Markets,
both zeroes exist, and the steepness parameters are large. Denote the zero
in the upper half-plane by i{3 (in Chapters 3, 5 and 7, it was denoted by
—i/?_). In the last Section, we will show that under these conditions, the
following approximate equalities hold:
where
(11.47)
ii>'(if3)<t>+(ipy
The case of a solvent firm 243
1.04v.-
0.31s
Fig. 11.1 Left panels: the default level Vg and the limit value of the yield spread. Right
panels: KoBoL model vs. Gaussian model. Parameters: r = 0.06, T = 0.35, a — 0.5;
vol=0.35, v = 1.2
Below, we will see that in order to get the correct dependence (in terms of
curvature) of the value of the firm and the value of equity on asset value,
K(p) should be negative, which is the case here. To see this, notice first
that
r+oo
/•-t-oo
tfW) E / e-rte~'3rntdt\mt =0
Jo
where {mt = s u p 0 < s < t Xt} is the supremum process for Xt (see Eq. (2.24)),
and so it is positive. Hence, it suffices to prove that iip'{i[)) is negative.
244 Endogenous default and pricing of the corporate debt
0 (e* - e-px).
*i(s)
0 +1
It is easy to verify that the following approximations for Eq. (11.27)-
Eq. (11.33), which we derive by substituting Eq. (11.45)-Eq. (11.46), are
similar to exact answers in Leland (1994a) (compare Eq. (7), (9) and (11)-
(13) there (it may be not so easily seen since the formulas in op. cit. are
written in different variables):
C ( l - r)
E = £7jz7) - 1 ~ tfW1 - C1 + P)-1)^^'-^
D
= V-
C(l-r)
(1 - a)ex + Cr~l{\ +
1 + K(fi)
m +1) (0
p~xK{p)e-0{-x-h^)
(11.48)
C(l-r)(l-a)(c,
+ (/? + ly^K^e-^-V)
C(l-T)
1 (l-r)(l-a) -0'
£ 1 + ^(0)
^ 0+ 1 (£) (11.49)
i> = F +
Cr (1 - T)a0 V_ -P (11.50)
0(0 + 1) 1 + 0 +
1 +
VB
KtJ$) (V
TB =
r 1+ 0
(0 (11.51)
BC =
C(l-r) Fr + K(fi) ( V \~ff
V
C(l-r) /3 + 1 UBJ
N _/3
0+ H0-
a C ( l - r ) AT(/3) / V
(11.52)
The case of a solvent firm 245
coupon
Fig. 11.2 Left panels: the value of the debt and the yield spread. Right panels: KoBoL
model vs. Gaussian model. Parameters: r = 0.06, T = 0.35, a = 0.5; vol=0.35, v = 1.2,
A_ = -30.
Fig. 11.3 Left panels: equity and debt value. Right panels: KoBoL model vs. Gaussian
model. Parameters: V/VB = 1.2, C = 10, r = 0.06, a = 0.5, r = 0.35, vol=0.35, u = 1.2.
to the bankruptcy level (see Fig.ll.3-11.4). We fix the ratio of the firm
assets to the bankruptcy level V/VB = 1.2, which makes the approximate
formulas Eq. (11.48)-Eq. (11.51) applicable, with a small error, if A+ and
—A_ are not small; and this is the case we consider here: both A+ and
—A_ vary from 30 to 8. We see that the non-Gaussian equity value and the
Gaussian one can differ by 23 percent. In this series of examples, the value
of the debt in non-Gaussian model can be larger by 6 percent than in the
Gaussian model, the same holds for the firm's value, and the yield spread;
the leverage differs insignificantly.
Endogenous debt level and endogenous leverage 247
v* = V +— - ^ - , (11.54)
r p+ 1
D* = — ^ r , (11-55)
r p +1
where
1
m = 1+
p + 1 + (1 - r ) a / 3 / r '
r(/3 + 1)
R* = C*/D* =
mp
which is again the same expression as in the diffusion case, but for a different
P-
Numerical examples show that as far as —A_ and A+ are not small,
only the endogenous coupon and the default level can differ by 3 percent;
the values of other endogenous variables change less when we replace the
Gaussian model with the non-Gaussian one, of the same variance.
248 Endogenous default and pricing of the corporate debt
11.6 Conclusion
*i(0 = ^ ( 0 ( ^ - j ) _ 1 -
250 Endogenous default and pricing of the corporate debt
(ic-jr^a+^r'+oder2),
therefore it follows from Eq. (11.34) that for j = 0,1,
= r(/c)(i + iO"K,
and from the definition of \&j, j = 0,1:
-oo+i<r
for any a < -j, and from Eq. (11.56)-Eq. (11.57) it follows that
We know that \Pj(x) = 0 for any x < 0, hence from Eq. (11.58), Fj also
vanishes on (-oo,0]. Since Fj G CK-+P(R), for any p G (0, v - i>i), we
conclude that a s i - 4 +0,
Eq. (11.58)-Eq. (11.57) prove Eq. (11.36)-Eq. (11.37), and Eq. (11.38) is
proved similarly, the equality
-oo+i<r
-00+t<7
V>(£) admits the analytic continuation into the complex plane with cuts
(—ioo-NA_], [i\+, -H'oo), and it has only one zero i/3 in the upper half-plane.
Therefore we can use the Wiener-Hopf factorization formula Eq. (2.23) to
define the meromorphic continuation of (j>~ (£) into the complex plane with
the cut [i\+, +ioo) and the only pole at k = i/3. Notice that as £ -*• i/3,
0r(O =
r
+
(r + WP) + i>'(il3)(Z ~ iP))&W) °(1)
T
= K{P){Z-iPrl+0{l). (11.60)
On the strength of the residue theorem, and from Eq. (11.60), we obtain
and similarly,
-1 J eix^;
(27TJ)- 1 / e ^ - ^ + ia)-1^
252 Endogenous default and pricing of the corporate debt
r+oo
= (27r)-1/ e-"*{z){z + a)-ldz,
where
For model classes of RLPE, $(z) grows at the infinity not faster than Cz"
- see, e.g., Eq. (11.35) - and hence, we conclude that the integrals in
Eq. (11.61)-Eq. (11.62) admit a bound via
+oo
r+oo r+oo
/-t-oo
zx l l V-\J0-ZX
l
C / e- z"- dz = -Cx~ j z"- de
+oo
= Cx-1A+"-1e-A+I+Cx-1(j/-l) f e-"z"-2dz.
JxA+
If A+ is large and x > 0 is not too small, then we conclude that the integrals
in Eq. (11.61)-Eq. (11.62) are small. For instance, if x = 0.2 (that is, we
are 20 percent from the default level), v is close to 1 and A + = 40, which
is quite often the case in the real financial markets, then the upper bound
for the integrals is less than 0.01.
Auxiliary results 253
Fig. 11.4 Left panels: firm's value, yield spread, and leverage. Right panels: KoBoL
model vs. Gaussian model. Parameters: V/VB = 1.2, C = 10, r = 0.06, a = 0.5,
T = 0.35, vol=0.35, v = 1.2.
This page is intentionally left blank
Chapter 12
Fast pricing of European options
We consider the Fast Fourier Transform (FFT), and a new method of op-
tion pricing based on the transformation of the line of integration into the
integral over an appropriate cut in the complex plane (integration-along-cut
method - IAC). The method is applicable for option pricing under VGP,
NTS Levy and KoBoL of order less than 1, and NIG.
12.1 Introduction
255
256 Fast pricing of European options
after reduction to the cut, the integrand grows at infinity, and hence, I AC
is not applicable. We have not checked whether IAC is applicable for all
GHP. Since the formula for the characteristic exponent of a GHP, which is
not a NIG, is rather complicated, and the relative advantage of IAC is small
for processes of order 1, we surmise that it is not rational to apply IAC to
the other GHP's (though this question deserves a more careful study).
The implementation of the method requires the detailed calculations in-
volving the characteristic exponent. We illustrate the usage of the method
by calculating the price of the European put option under the KoBoL pro-
cesses of order (0,1).
Consider the model market of a riskless bond and a stock. Let r > 0 be
a (constant) riskless rate, and let St = exp Xt be the price process of the
stock. We assume that X is a Levy process under an EMM Q chosen by
the market; ip is the characteristic exponent of X under Q. We assume
that the steepness parameters satisfy A _ < — 1 < 0 < A + .
Consider the European put with the strike price K and the terminal
date T. The payoff at the expiry equals max.{K — ST, 0}, and without loss
of generality, we may assume that K = 1. Then
+oo
e-fa€(i_e*)+dz
/ -OO
Set
dv = A+" + (-A_)",
K = A+ + (-A_),
X+ = A+/K,
X- = -A-/K,
£/ = K(X + T ^ ) ,
V = —K"TCU,
P = 1/K,
and denote:
R(x;a) = Rie-K^x+Tfi\
and
expjizTii - V[(X- -iz + a)" + (x+ +iz- <*)"]}
1p(Z:a) = ; ; r~; : r
V ;
(-iz + <r)(-iz + a + p)
then
+oo
eixzip(z; a)dz. (12.8)
/
-oo
For a chosen step of numerical integration r), we use the trapezoidal rule
for the integral Eq. (12.8) and obtain an approximation for f(x,t):
N
f{x, t) » R(x; a) £ e*"' ip(Zj; a)v, (12.9)
FFT and IAC 261
xk = -b + X{k-l), k = l,2,..,N,
N
f(xk,t) » R(xk;<T) £iei*'<--b+x<k-1H(zi;<r)TI,
y
(12-10)
j=i
N
f(xk,t) » Rix^a)^^-1™-1^-^^;^- (12-11)
j=i
A*? = f . (12.12)
(12.13)
where 5n is the Kronecker symbol.
The sum in Eq. (12.13) is of the form Eq. (12.6), hence it can be com-
puted by FFT. Notice that it is rather difficult to estimate the error of this
computational scheme.
12.3.2.1 Case a)
If U > 0, we push the contour of integration to the cut in the upper half-
plane, the choice of a cut being defined by a natural requirement: iU£ is
negative at the cut. As a result, we obtain
r+<
/•-t-tc
°° exp[iUt] ,-V>(£+0) _ -VM-O)
f(x,t) i?, d£,
- * ( - « £ + p)
where <£(£) = ( x - - i O " + (x+ + *0"-
By changing variables £ —> iz, z S R, we obtain:
/•+0
+°° exp[-£/z] e-V^>(iz+0) _ e-V<t>(iz-0)\ dz_
/(x,t) = «!
*(* + P)
For z > x+, we have
Set
then
v
= - 2 t e x p ( - V ( X - + z) - V(z - x+)"a„) sin(V(z - X+YK),
and
+°° exp[-E/z] s i n ( y ^ ( z - X+Y)
f(x,t) 2i? x
x+ L z z
( + P)
exp ( - V [ ( x - + z)" + au{z - X+T]) dz.
(12.15)
+
" °° exp[-£/z] sm(V&„z")
2Rie-u*+ / exp ( - V [ ( l + z)v + auzv\) dz.
Jo (z + x+)(z + X+ + P)
FFT and I AC 263
(12.16)
where u > 0, v > 0, pi > 0, pi > pi, v e (0,1). Thus, in the case U > 0,
the computation of f(x,t) reduces to the computation of the integrals in
Eq. (12.16):
12.3.2.2 Case b)
In the case U < 0, we push the contour of integration in Eq. (12.2) to
the cut in the lower half-plane. In the process of the transformation, the
contour crosses two simple poles £ — 0 and £ = — i. Using the residue
theorem, we obtain
Hence
f(x,t) = e x p ( - r r ) — exp(x)
r+ v
r °° exp\Uz] sm(VK(z - X-) )
Jy-
• exp (-V[(X+ + z)v + av(z- x _)"]) dz.
f J. V j . vexp(-V[(l + z r + a ^ ] ) d z ,
Jo (z + X-)(-2 + X - - P )
and use the notation of Eq. (12.16):
We conclude that any numerical procedure for the integral Eq. (12.16),
together with Eq. (12.19), can be used to compute f(x, t) in the case U < 0.
Example. In the case v = 0.5, we have av = 0, bv — 1, and the integral
becomes especially simple:
+oo
l(0.5,u,v,Pl,P2)= f * ' " s m ^ ) exp {-vy/T+7) dz. (12.20)
7 (^ + Pl)(^ + P2)
0
12.4 C o m p a r i s o n of F F T a n d I A C
(a)
0.9 - (b)
0.8
0.6 - / • • "
0.5 - / /
0.4 -
^ ^ ^ ^
0.3 - ^ — - ^ ^ ^ - - • • - • • • " " • "
0.2 - '"
0.1 -"
'
In discrete time, one can make very weak assumptions on the underlying
stochastic process, and yet obtain analytical solutions; at the same time,
the discrete time approach makes it possible to solve the pricing problem
for non-standard American options known as Bermudan options. A Bermu-
dan option can be exercised on specified dates only, during the life of the
option. An example of a Bermudan option is a bond option whose exer-
cise is restricted to coupon payment dates. Consider a Bermudan option
which can be exercised at the prescribed dates tj = jAt only, and the cor-
responding perpetual American option in discrete time: the answer for the
latter option gives the optimal exercise price and the price of the Bermudan
option on the prescribed dates, the price of the Bermudan option at time
t € (tj,tj+i) can be calculated as the price of the European option with
the expiry date i,+i, and the price of the Bermudan option at time tj+i as
the terminal payoff.
The state space in our model is R, and not a discrete set as in multino-
mial models (for these models, see, for example, Duffle (1996) and Shiryaev
267
268 Discrete time models
where A_ < — 1 < 0 < A + , and there exist positive constants e 0 ,Cb, and p
270 Discrete time models
such that
Notice that Eq. (13.4) implies that P(dx) is absolutely continuous: P(dx) =
p(x)dx, and Eq. (13.4) is satisfied if p is continuous, differentiable a. e., and
+oo
e£°lxy(x)|dx < + o o .
/ •oo
Eq. (13.4) can be relaxed but we believe that this condition is weak enough
for practical purposes.
where M denotes the set of all stopping times r = T(U;), satisfying r(w) <
o o , {J € Q.
As in Chapter 5, first we will find a solution to Eq. (13.5) in the class
Mo of hitting times r{a) of segments (—oo,a]. To this end, we consider
V(h; x), given by Eq. (5.21), which is the price of the put, when h is chosen
as the exercise boundary. The sufficient conditions for the optimal exercise
boundary are given by Lemma 5.2. Since this lemma is stated in terms of
the value function V(h,x), we need an explicit formula for it. To obtain
such a formula, one can proceed as follows.
Fix h, a prospective candidate for the optimal exercise boundary, then
the option pricing problem can be written as
Thus, the original problem has been reduced to the Wiener-Hopf equa-
tion. As in the previous chapters, we can obtain the solution to Eq. (13.10)-
Eq. (13.11) using the Wiener-Hopf method. After the solution is obtained,
it is easy to guess the necessary condition for the exercise boundary. In
fact, under our standing assumption on the process, the optimal exercise
boundary is determined by the requirement: V(h; •) is continuous at h. We
call this requirement the continuous pasting condition as opposed to the
smooth pasting condition, which never holds in discrete time.
It can be shown (for details see Boyarchenko and Levendorskii (1998))
that the necessary condition is similar to Eq. (5.31), so that an explicit
formula for h obtains. After that, one has only to check that the sufficient
conditions hold for this h. This is done exactly as in Subsection 5.2.2. Fi-
nally, by using the following sufficient condition from Darling et al. (1972),
one can prove the optimality in the same fashion as it was done in op. cit.
for calls.
272 Discrete time models
and
then
and
°° ~k yO
(13.18)
Set
C«*T1 t (13.19)
(13.20)
0+(O = E[e*Xr]
/•+00
= Y V l - z)zn / eixidP(Xn < x)
= / e-<*tp+(dx),
J — OO
where
This proves a) for the sign "+"; the proof for the sign "-" is similar, and
b) follows from a). •
6) <A+(£) is holomorphic in the half-plane S£ > <r_, and can 6e defined, for
any w_ > <r_ and S£ > u;_, 6j/
1
ex /-+°°+^- e - V ( - ^ ) .!?-£.
<M£) = P (13.23)
oo+iu,_ l - e - r p ( - n ) n
+0o+ia,_ gin[i_e-rp(-T7)]
= 6XP dn (13.24)
p/- «(£ - v)
«+(0) (13.25)
MO'
where
•1 +oo+tu;_
/-t-OO- lnfl-e-^-n)]
a+(£)=exp dr] (13.26)
2fft . / _ „ , +
c) $-{£,) is holomorphic in the half-plane Off < o~+, and can be defined, for
any w+ < cr+ and S£ < u>+, by
ex (13.27)
<M0 = P 27T*
m 7y__0 0 + j w + l_e-pp(-T7) n
+oo+iw+ glnfl-e-^-n)]
exp d77 (13.28)
2iriJ_ oo+iw+
MO) (13.29)
MO'
where
a_(f) = exp
1
2m f/•+°°+i"+ ln[l - e " p ( - " ) ]
-oo+iuq. "-£
dn (13.30)
°° zk f°° f+°°
= E l / (e"C - 1 ) ( 2 f f ) " 1 / _ e-ixnP(-v)hdvdx
Assume that Eq. (13.4) holds with p > 1. Then the integrand admits a
bound via the integrable function Ceu-X(l + |»/|) - ''. Hence, we may apply
the Fubini theorem and integrate w.r.t. x first:
(In the case p G (0,1], we have to resort to the theory of oscillatory integrals
in order to justify the change of order of the integration). As o;_ —> 0,
hence, for any z G (0,1), there exists o>_ G ( c - , 0 ) such that zC(wJ) <
1. This implies that ^2 J converges absolutely, and for such an u;_, we
can apply the Fubini theorem and change the order of the integration and
summation:
Inll-zpi-vMrj-t)-1-^1^
•oo+iw-
This proves Eq. (13.24)-Eq. (13.25), and Eq. (13.23) is deduced from
Eq. (13.24) by integration by part.
276 Discrete time models
c) is proved as b). •
a(0=o+(Oa_(0 (13.31)
b)For any w € (<x_,0), there exist C,pi > 0 such that for all £ in the
half-plane 9f£ > u>,
M O ^ - l l ^ C a + iel)-" 1 ; (13.32)
c) For any LJ G (0,cr+), there exist C,pi > 0 such that for all £ in the
half-plane 3£ < w,
Ci(l + \T)\rp\r)-tr\
where C\ and p > 0 are constants. By considering separately cases \TJ\ <
If 1/2) \T}\ > 2|£| and |f |/2 < \rj\ < 2|f |, and using the triangle inequality, it
is easy to prove that there exists C 2 such that for all rj on the line Sty = w_
and all £ in the half-plane 3?£ > w,
Ci(i+M)-"|77-f|-1<c2(i+|£|)-"/2[(i+|^-f|)-1-''/2+(i+H)-1-''/2].
The RHS is integrable over the line QTJ = w_, and the result is C(l +
|f \)~p/2, with some constant C; Eq. (13.32) follows.
c) Eq. (13.33) is derived in the same fashion. •
The discrete time model is more flexible than the continuous time model
since an empirical probability density can easily be approximated by an
appropriate model one so that the factors in the Wiener-Hopf factorization
The Wiener-Hopf factorization 277
p_(-tj) = f eixr,p_(x)dx
J —OO
oo
r0 m+
ixr,+\+x(_xydx
•°° 3= 1
+
m
= 5 ^ c + r ( j + i)(A+ + tf 7 )-'- 1 ,
3=1
1_ e-rp(-r))
l n ^ i = ln(l-^/r,) = 0(|r,r1)
as T] —> oo, we conclude that on any line S77 = w, which does not pass
through a root of P Q , the integrand in Eq. (13.27) admits an estimate via
C(l + |T?|) - 2 , and the same holds in the half-plane S77 > u>, if there are
no roots of PQ in this half-plane. Therefore, we can apply the residue
theorem and shift the line of the integration: w+ —* +00. Each zero of PQ
contributes to the resulting formula, and in the result, we obtain that 0_
is a rational function. Similarly, 4>+ is, and since <£_ (resp., <f>+) is bounded
and does not vanish in the half-plane 3£ < 0 (resp., S£ > 0), we conclude
that
* - < « = n r ^ n ^ . (13.38)
where {—z/3j} (resp., {—ijj}) are all the zeroes of P (resp., Q) in the upper
half-plane. Similarly, <£+(£) is defined by zeroes in the lower half-plane.
In this Section, we present the analytical formulas for the optimal exercise
price and the rational price of the perpetual American put option in discrete
time.
27ry_0o+ia,, (-if)(l-i£)
(i) g is continuous on R;
(ii) g is bounded and positive in some neighbourhood of - c o ;
(iii) there exists /i* such that
for all x < h*, (4>-(D)~1g)(x) > 0, and
for all x > h,, ((A-(J5) _1 /)(x) < 0.
Hence the result can be generalized for any payoff functions satisfying the
properties above. These properties are most easily verified for the payoffs
given by Eq. (5.81). It is also possible to impose additional conditions on
the parameters jj in Eq. (5.81) to obtain the optimal solution in the class
M.\ the conditions and proof (based on the reduction to the case of the put)
are essentially the same as in Chapter 5, where the continuous time model
is studied.
Below we give the simplest example, when a discrete time process is not
supposed to correspond to the observation of a Levy process (in continuous
time).
p ± (*) = ^ l R ± ( x ) e A ^ ,
1 H-oo \
and
a(r)) = 1 - e r p(-7?)
280 Discrete time models
= 1-
2e r A_ — it]+ A+ + irj
2(1 - e r )A_A + + (1 - 2er)(A+ + \J)vt\ + 2er„2
r
rj
r
~ 2e (A + + iii){-\- - irf)
Denote by — i/3_ (resp., — i/3+) the root of the numerator in the upper
half-plane (resp., lower half-plane); then
i ( 2 e r - l ) ( A + + A_)±iy/]9
-ifc = 4e r
D := (2e -l) (A+ + A_)2-16er(er-l)A+A_,
r 2
and we have
(77 + »/3_)(ij + i{3+)
ofo) = 2er(r]-i\ )(T]-i\-)'
+
* /« (A + + i Q ( - / 3 - )
=
*-«> A+(-/3_+i0-
Insert into the formula for the optimal exercise price:
(A+ + l)(-/3_)
Hm = eh' = K
A + (-/3_ + l) '
and then into the formula for the price of the perpetual American put: for
x > /&*,
K ,+oo+iW+ e <(«-J».)€(A + +ifl(-/?_)
«.
'•w-5/. ,+iui+ A+HL+ifl(-iO(l-*)
Under condition Eq. (13.22), io+ hes below — i/3_, hence iu; + does as well,
and when we push the hne of the integration up, it crosses the pole — i/3_.
By the residue theorem, we obtain
K{\+ + p_)e^^-h')
K(x) =
A+(l-/3_)
Chapter 14
14.1 Introduction
281
282 Feller processes of normal inverse Gaussian type
1
There is one useful property which generators of NIG and other RLPE processes fail to
have: the transmission property, Boutet de Monvel (1971), or a weaker smoothness-in-
a-half-space property, Eskin (1973), which simplifies the treatment of boundary value
problems. This property ensures that a solution to a "sufficiently good" boundary value
problem, e.g., the Dirichlet problem, with smooth data is smooth up to the boundary.
Introduction 283
may either decrease the scale parameter 8 (which is in direct analogy with
a standard device in Gaussian modelling) or increase a. For instance, one
may expect that if for some c > 0, p > 1,
then the trajectories of the process never reach 0 from above. Processes of
this type are suitable for interest rate modelling purposes. To treat such
processes, the calculus of PDO with degenerate symbols is needed; for such
a calculus, see Levendorskii (1993) and Levendorskii and Paneah (1994).
In this book, we consider the simpler case of bounded a, 0, 5, fj., assum-
ing that 6 and a±0 are bounded away from 0. In addition, we assume that
these parameters are smooth and their derivatives are small. The simplest
example which we have in mind, and which is capable of reproducing the
mean-reverting effect, is the case of constant 5 and a, and 0 given by
where e > 0, x > 0, and \0o ± xl < a- The derivatives of 0 admit estimates
\al;](x,0\<Cap(Om-M, (14.4)
Theorem 14.1 Let Eq. (14-4) be valid on R n x R", and let there exist
c> 0 such that on R " x R n
Then there exists A > 0 such that a(x, D) + A : <S(R") —> <S(Rn) is invert-
ible.
2
T h e more standard notation would be S J " 0 ( R n x ( R n + iU)); we have omitted t h e lower
indices to simplify the notation.
Constructions of NIG-like Feller process via pseudodifferential operators 285
Notice that it is a process of order 1 and exponential type [A_, A+] for
any A_ < A+ satisfying
Here (•, •) is a bilinear form in C n , which extends the standard scalar prod-
uct in R n .
Let a be defined by
K(W)=
/ (l-e^'Xtis). (14.11)
Jo
To show that A is a PDO and calculate its symbol, we notice that Eq. (14.9)
can be written as
/•OO
Then
a) n(-L)+d is a PDO with the symbol of the class 5 1 ( R n x (Rn+iU)),
and the symbol admits an asymptotic expansion Eq. (14-13) in the sense
that for any N > 0,
N-l
b
rN = b-J2 i-i e S1-N(Rn x (Rn + iU)), (14.15)
j=0
where
For the proof and inductive formulae for bj, j < 0, see Chapter 16.
Theorem 14.4 means, in particular, that NIG-like Feller processes con-
structed by subordination from diffusions with the uniformly elliptic — L,
where L is the generator of the diffusion, are iVIG-like processes in the
sense of the "naive" definition. Moreover, if we use the symbol a(x, £) =
(d 2 — L(x, ^)) 1 / ' 2 — d to construct a NIG-like process, and if the coeffi-
cients of L are slowly varying and/or d is large, we obtain a process whose
generator differs insignificantly from -K(-L), the generator of the process
obtained via subordination.
(dt + L-r))f(t,x)=0,
Lemma 14.1 For any t > 0, exp[—ta(x, D)} is an integral operator with
the kernel of the class C°°.
The proof will be given in Chapter 16. Lemma 14.1 implies that we can
apply the reduction of Chapter 2 to NIG-like Feller processes, and deduce
the generalised Black-Scholes equation
f(t,x) = e-Tr(QTg)(x),
Psg(x) = esr<l>(T-s,x),
under a diffusion process, and then integrating e~rTPsg(x) w.r.t. the mea-
sure XT(ds) to obtain
f(t,x)= f e^-T>4>{T-s,x)\T{ds).
7(0,00)
Recall that one should choose the parameters of the model so that the
EMM-requirement
be met.
formula for QT. That can be derived from the representation theorem for
analytic semigroups Theorem 15.18 similarly to Chapter 8, where the tran-
sition semigroup was computed for the case of an operator of a boundary
problem for a generator of NIG and some other Levy processes, in order
to compute the rational prices of barrier options.
In order that the RHS in Eq. (14.20) and Eq. (14.22) be finite, some
regularity conditions on g are needed, and these conditions are related to
the domain U in the definition of the generator of a NIG-like Feller process;
for instance, in ID when U is an interval (A_, A + ), it suffices to assume that
g is piecewise continuous and satisfies estimates
with A_ < w_ < w+ < A + . In the case of European puts and calls, it
amounts to the restriction on A_,A + : A_ < — 1 < 0 < A + . Due to the
put-call parity, it suffices to calculate the price of the European put, for
which g(x) = (K — e x ) + , where K is the strike price. Thus, below, we
discuss the pricing of the European put.
We consider a NIG-like Feller process in ID with the generator — a(x, D),
where a is defined by Eq. (14.1), and n,5,a,(3 satisfy Eq. (14.6) and
Eq. (14.7). The EMM-requirement reduces to
r + a(x,-i)=0. (14.24)
If 5,a,f3 are already chosen, one can use Eq. (14.24) and Eq. (14.1) to find
the drift fi.
The construction of the transition semigroup starts with
Proof. Fix large C and consider first the region Hc,e — {^ I \M > C, arg A G
[-6,0]}. Since
for C fixed, we can easily find 6 G (7r/2,7r) and C\ such that for all A 6 Sc,e
and all (ar.f) G R x (R +1{0}),
lACA + r + a f o O r 1 ! ^ ! . (14.26)
By using Eq. (14.26), and arguing as in the proof of Theorem 14.1 (cf. also
Theorem 16.12), we obtain that Eq. (14.25) holds for A G £c,e, w i t h some
Co independent of these A.
To prove that for C fixed, there exists 6 G (TT/2,TV) such that Eq. (14.25)
holds for A satisfying |A| < C and argA G [—6,9} it suffices to notice that
since — a(x, D) is a generator of the Feller process, r + Ua(x, D) is positive-
definite.
Remark 14.1. Since the boundedness theorem is valid for much more
general spaces, the estimate Eq. (14.25) and the formulas and estimates
below are also valid in these spaces.
The estimate Eq. (14.25) means that the representation theorem for the
semigroup with the generator —r — a(x, D) is applicable, and the following
formula is valid:
Here £# is the contour A — A(CT), —OO < a < +oo, where arg X(a) = —6, a <
0, argA(o-) = 0,er > 0.
By using an asymptotic expansion of the symbol R\{x,£) of the resol-
vent (A + r + a ( x , D ) ) - 1 :
(for the proof and the full asymptotic expansion, see Chapter 16), we can
compute the RHS in Eq. (14.27) in the form of a series. Since
(27ri)" 1 / eAT(A + r + a ( a ; , 0 ) _ ^ A - e x p [ - r ( r + a ( x , 0 ) ] ,
JCe
the leading term in the asymptotic expansion of the RHS in Eq. (14.27) is
-foo+itr
exp[ia£-T(r + a(x,0)]5(0#.
/ •OO+tlT
Applications for financial mathematics 293
To sum up: a formula for the leading term for the price of the European
put with the strike price normalized to 1 is
In Chapter 16, we give formulas for the next terms of the asymptotic expan-
sion Eq. (14.29), and show that the first omitted term is already fairly small:
if the symbol a(x, £) depends on parameters (e, d) and satisfies estimates
for any a, (3, where Cap is independent of (e, d), then the first omitted term
in Eq. (14.29) admits an estimate via Cr2e, where C is independent of
(T, e, d), and the next term is less still: of order r 2 e 2 / d .
Notice that the leading term looks exactly as in the case of the option
pricing under Levy processes, only when we calculate f(t,x), we use the
characteristic exponent depending on x. Thus, the leading term can be
calculated as easily as in the case of Levy processes.
Remark 14.2. Consider the NIG-hke Feller process defined by the symbol
Eq. (14.1) with constant 8, parameters fi,a,P being bounded and having
bounded derivatives, and satisfying
15.1 Introduction
295
296 Pseudo-differential operators with constant symbols
explicit solutions are not available, and this is a situation, where a general
theory can establish the existence and uniqueness of the solution in an ap-
propriate function class, prove the regularity of the solution, and obtain
effective approximate formulas for the solution. In Chapter 8, we demon-
strated a realization of this programme for pricing of barrier options in the
simple case of RLPE; this Chapter contains additional potentially useful
constructions and treats more general classes of symbols and spaces. If one
repeats the same constructions in the case of Levy-like Feller processes and
boundary value problems for contingent claims (and uses necessary results
of Chapter 16 on PDO with variable symbols, like in Chapter 14), one
obtains approximate pricing formulas.
We give a collection of main definitions and results of the theory of PDO
with constant symbols, without the connection to Probability, since this
connection was illustrated in many places in the main text. The exposi-
tion is essentially self-contained; of the reader, only the knowledge of basic
Functional Analysis and Complex Analysis is required. All necessary re-
sults from the theory of generalized functions and function spaces are listed;
in many cases, proofs of the results of the theory of generalized functions
are given. We omit some long proofs if the technical tricks in the proofs are
not used later, and the proof can easily be found in many monographs. We
basically follow the exposition in the monograph Eskin (1973), with some
modifications.
SU
\\f\\c»(K;B) = £ P ll/ (Q) (*)lle.
|a|<m
By using the Leibnitz rule and the triangle inequality, it is easy to show that
k*v £ <S(Rn), and that a fixed k € <S(Rn) defines a continuous operator
K : 5 ( R n ) 3 v H4 k * v € <S(Rn).
In the literature, one can find many versions of the definition of the Fourier
transform. For instance, Eskin (1973) uses Eq. (15.2) with ix£ in the expo-
nent (hence, Eq. (15.3) with — ix£), and the factor {2ir)~n can be shared be-
tween T and T~1 differently; but in all cases, if dx and d£ are the Lebesgue
measures in R ^ and R n £, respectively, the product of the factors in front
of the integral signs in Eq. (15.2) and Eq. (15.3) must be (27r) _n .
The following properties of the Fourier transform are well-known.
in particular,
HfillL(R-) = (2*)n|Mli2(R")-
The Fourier transform Eq. (15.2) is well-defined on a wider space L i ( R n ) ,
and by using the Parceval identity Eq. (15.4), T can be extended from a
dense subset <S(Rn) C Z<2(Rn) to a topological isomorphism of L 2 (R")- In
the next subsection, a natural generalization of Eq. (15.4) is used to define
the Fourier transform in the space of generalized functions. This construc-
tion has a disadvantage of being indirect; in some important cases, one can
obtain analogues of the direct formulas Eq. (15.2)-Eq. (15.3); the corre-
sponding procedures are considered in Subsection 15.3.3-Subsection 15.3.4.
L e m m a 15.1
(i) If two functions, f and f\, define the same regular functional and are
continuous on an open set Q, then f(x) = / i ( x ) , V i 6 f i .
(ii) If f and / „ , n = 1,2,... are measurable, fn(x) -4 f{x) a.e., f satisfies
Eq. (15.5), and |/„(x)| < f(x), Vn,x, then / „ - > • / in <S'(R n ).
\a^(x)\<Ca(x)m°, (15.6)
(a/,u) = ( / , a u ) , Vue5(R").
302 Pseudo-differential operators with constant symbols
A(x,D)= £ aa(x)Da,
\a\<m
f*g = F-1(fg),
Space S'(R n ) of generalized functions on R" 303
and then
Hf*g) = fg-
|fi(OI<C(0', (15-9)
where C, s are independent of f. This is exactly the situation, which we
encounter in the next Section, where the solutions to pseudo-differential
equations will be found by means of the Fourier transform and its inverse.
Regularization 1. The generalized function u, the inverse Fourier trans-
form of it, can be defined by
(u.i;) = (2ir)-n{u,v)
= (2TT)-"/ fi(0«m
= (2rr)- n f d£ f dxe-ix^u(0vW),
JR" JR"
where the integral is understood as the iterated one. Take a positive integer
N, notice that
(0-2N(l-&x)Ne-ixS=e-ixZ,
304 Pseudo-differential operators with constant symbols
The formula Eq. (15.10) allows one to explicitly calculate the pairing of
the generalized function u, given by Eq. (15.3), and any test function v €
<S(Rn). This defines u. It is easy to see that Eq. (15.10) with 2N > s + n
defines a continuous linear functional over <S(Rn).
Regularization 2. Take a sequence of uniformly bounded measurable
functions tpn, converging pointwise to 1. Then consider integrals
/„ = (27r)-" f f e-**il>n(du{CMx)dZdx.
JRn JRn
(/i®/2,w) = (fi,ui)(f2,u2),
and extended by linearity to 5(R", _ 1 ) ®<S(Ra;n).
Pseudo-differential operators with constant symbols on R n 305
KOI<<7<0"V (15.11)
where C,m are independent of £. Then we write a € S' m (R n ). The-
multiplication-by-a operator is continuous in <S'(R n ), therefore the compo-
sition Jr~1aJr is a continuous operator in 5 ' ( R n ) . It is denoted by a(D),
and A = a(D) := !F~1aJr is called a pseudodifferential operator (PDO)
with the symbol a, of order m. The action of A = a(D) on <S(Rn) can be
306 Pseudo-differential operators with constant symbols
defined by
or equivalently,
Proof. Since the Fourier transform of elxy is equal to (27r)n<57, and a<57 =
a(7)(57, we have
T-xaTeix~< = T-Xa(j)(2v)n5y
= a(7)JP-1(27r)n<57
= a(i)eix"<.
b) Let the Fourier transform of k € S'(Rn) satisfy Eq. (15.11). Then the
convolution operator K with the generalized kernel k is a PDO with the
symbol k.
Lemma 15.3 For m,s 6 R, A S (D) : S(Rn) -> <S(Rn), and A S (D) :
<S'(Rn) -» S'(R n ) ore topological automorphisms, and A S (D) : i f m ( R n ) ->
i f m _ s ( R n ) is an isometry, with inverses defined by \~S(D).
Proof. Evident. •
In the following theorem, the main properties of the scale of the Sobolev
spaces are collected; the proofs can be found in, e.g., Triebel (1978).
|a|<m
Since all the terms are non-negative, we conclude that ||u|| s < oo iff
||Z) a u||| ,Rn* < oo, \a\ < m, i.e., the first statement of Lemma holds.
The second statement is immediate from Eq. (15.19). •
where C is independent o/£ € R n . Then the operator Eq. (15.20), and the
operator
A:S'(Rn)^S'(Rn) (15.22)
are invertible, the inverse A~* being a PDO with the symbol a ( £ ) - 1 .
Proof. By construction, AA~l and A-1 A act as the identity, and from
Eq. (15.21) and Theorem 15.4, A " 1 : H"-m(Rn) -> H°{Rn) is bounded.
Hence, it is the bounded inverse of the operator Eq. (15.20). Since both
A, A'1 are continuous in S ' ( R n ) , the operator Eq. (15.22) is invertible as
well. D
A symbol a € S m ( R n ) satisfying Eq. (15.21) and the operator A = a{D)
will be called elliptic. Usually one calls a and A = a(D) elliptic if Eq. (15.21)
is satisfied outside some compact.
Au = f,
u = a(£>)- 1 /,
or
The proof is essentially the same as the proof of the following theorem in
Eskin (1973), Theorem 4.2.
Theorem 15.7 Let s > 1/2 and n > 1. Then the restriction-to-the
hyperplane -xn = 0- operator
Theorem 15.8 (Eskin (1973), Theorem 4.3) Let s > r + n/2, where
r > 0 is an integer. Then Hs(Rn) C Co(R n ), densely and continuously.
IMIc= E H^Hc-
\a\<m
These spaces are called Holder spaces. The Holder spaces constitute a part
of the scale of the Holder-Zygmund spaces {C s (R n )} s 6 R (for the definition,
see, e.g., Triebel (1978)), namely, for non-integer positive s, C s (R n ) =
C s ( R n ) 1 . The scale of the Holder-Zygmund spaces enjoys the standard
property C s (R n ) C C ( R n ) , if s > r, and the Sobolev embedding theorem is
valid for the Holder-Zygmund spaces (see, e.g., Triebel (1978)):
Theorem 15.9 Let s > n / 2 . Then for any s' < s - n / 2 , H9(Rn) C
Cs ( R n ) , densely and continuously.
^ o t e that for integer s > 0, neither C s ( R n ) nor C^(R n ) coincide with C"(R n )
Pseudo-differential operators with constant symbols on R n 311
Theorem 15.10 Let s,m e R, and let a € C°°(R n ) satisfy, for any
multi-index a,
where Ca is independent of £.
Then A = a{D) : C s (R n ) -> C s _ m ( R n ) is bounded, with the norm
admitting an estimate
|a|<JV«eR"
Proof. By Theorem 15.5, A is invertible in S' (R"), with the inverse given
by A'1 = a ( D ) - 1 . Prom the Leibnitz rule, and Eq. (15.23) and Eq. (15.21),
it follows that a ( £ ) - 1 satisfies Eq. (15.23) with — m instead of m. By
Theorem 15.10, A : C S (R") -+ C s - m ( R n ) and A'1 : C 8 - m ( R n ) -> C s (R n )
are bounded, hence, they are mutual bounded inverses. •
o *
15.5.1 Spaces i f ( R n ± ) and PDO in a half-space
Lemma 15.5 (Eskin (1973), Lemmas 4-2 and 4-3). H(R™±) is equal to
the closure ofC$°(Rn±) in H"(Rn).
Theorem 15.12 Let a £ S m ( R n ) admit the analytic continuation w.r.t.
£n into the lower half-plane r = 9 £ n < 0 (resp. r > 0) and satisfy an
estimate
± W
I ((OT^n) m , ifn>l.
Clearly, A™ is holomorphic w.r.t. £„ + i r in the half-space ± r > 0, and
satisfies Eq. (15.25), hence
o * o *—m
Am(D) :ff(R" T ) - > # (R%)
is bounded. The same holds with —m instead of m, and A™A±m = 1;
therefore, it is an isomorphism.
by
||/lh.=inf||//||.,
where infimum is taken over all extensions If G H s ( R n ) . Set p± = PR.»±.
It maps # s ( R n ) onto H s ( R n ± ) .
Theorem 15.13 Let a G 5 m ( R n ) admit the analytic continuation w.r.t.
£n into the upper half-plane r = 9£„ > 0 (resp. r < 0) and satisfy an
estimate Eq. (15.25) for all £ G R n , and r > 0 (resp., r < 0).
Then the operator
Hs(Kn+) 3 u 4 p+a(D)lu € if s _ m ( R n + ) (15.26)
(resp., the operator
Hs(Rn_) 3 u H-> p_a(D)lu € tfs-m(Rn-))
is well-defined and bounded.
Proof. Consider the case "+". Let I'u 6 Hs(Rn) be another extension
o «
of u. Then supp (lu-l'u) C R " _ , hence lu-l'u G H ( R " - ) , and a(D)(lu-
o s-m
I'u) 6ff ( R n - ) , by Theorem 15.12. Hence, p+a(D)(lu - I'u) = 0, and
the operator Eq. (15.26) is well-defined. To show the boundedness, take lu
with the norm ||Zu||s < 2 | | U | | R " + ; S ; then from Theorem 15.4,
||p + a(D)iu|| R n + ; 8 _ m < ||o(D)Ju||,_ m < C\\lu\\B < 2C||u|| R « + ; .. Q
II''«II'.<C|MIR..
314 Pseudo-differential operators with constant symbols
(see, e.g., Subsection 2.9.3 in Triebel (1978) or Lemma 4.5 in Eskin (1973)).
Hence,
iMiR«±,.<r«n;<ciMiR. ±.s-
The opposite inequality is valid for any extension operator: by Lemma 15.4,
the norms || • || s and || • \\'s in i ? s ( R n ) are equivalent, therefore
T{9+{xn)e-TX")(£n) = 1 / ( T + i&,),
By using the Lebesgue theorem, we see that the limit of the LHS in Eq. (15.28)
is the Fourier transform of 6+f, and therefore,
F{0+f) = n + / , (15.29)
where
I I + / ( 0 := (2W)- 1 f ^ /^' V n )
drjn.
Similarly,
where
+o
° ke,vn)
n-/(e = - M - / J&iO-n,
•dr]n.
By using Eq. (15.29) and Eq. (15.30), and a series of rather technical esti-
mates, Eskin (1973) deduces a series of important results.
Theorem 15.14 (Eskin (1973), Theorem 5.1) For \s\ < 1/2, 0± admits
a unique continuous extension 8± : HS(R") —¥H(R-n±)-
Lemma 15.7 (Eskin (1973), Lemma 5.4) For \s\ < 1/2, any function
o *
/ G Hs(Rn) admits a unique representation / = / + + / _ , where f± EH
(R"±), and f± = 6±f.
Define the extension-by-zero operator e± : <S(Rn±) —>• I/2(R n ) by
e±f{x)
= { 0, ±xn < 0.
Prom Theorem 15.14 and Lemma 15.7, we deduce
Lemma 15.8 Let \s\ < 1/2. Then e± admits a unique bounded extension
s n
e± : H (R ±) -> Hs(Rn), and Hs(Rn±) can be identified with # ( R n ± ) .
Theorem 15.15 (Eskin (1973), Lemma 5.5) Let f <E Hm+s(Rn), m > 0
integer, \s\ < 1/2. Then
m
9+f = £ AI*(D) [ ( p ' A ^ 1 (£>)/) ® So] + Al m (D)0 + A™(£>)/, (15.31)
fc=i
and
m
A
9-f = Y, +k(D) [ ( p ' A ^ 1 (D)f) ® So] + A+m(D)6_A™{D)f, (15.32)
jt=i
f = Y,fk®Dk~ls°> (15-33)
316 Pseudo-differential operators with constant symbols
where fk G ff*+*-1/2(Rn-1).
f(x,+0)=go(x). (15.37)
f{x,t)=e«>tr°{x,t)
g(x,t)=e°°tg''°(x,t)
where aCTo(£) := a(£) + a0 satisfies Eq. (15.34) and Eq. (15.38). Having
in mind the possibility of such a change of the unknown and the data, we
assume in the sequel that a satisfies Eq. (15.34) and Eq. (15.38).
/(£,0)=3o(£), (15.40)
(here / denoted the "partial" Fourier transform w.r.t x). By solving the
Cauchy problem Eq. (15.39)-Eq. (15.40):
\H\2mMs,r) = 1 1 (M + (0m)2S(02rm,V)\2^dr].
JR" JR
o(«-r) of
Spaces Hm,i (R n x R+) a r e defined by the analogy with H(R-n x R + ) , and
similarly to Theorem 15.12, we deduce from Eq. (15.44), that for any s,r,
0(-l/2-e,S')
9o ® <5o € H m , i (Rn x R+). (15.46)
o U/2-e,0)
belongs to Hm\ ( R n x R+)- This statement can be made more exact,
but even this form suffices to prove that / is well-defined as a general-
ized function. The initial condition can also be understood in the sense of
generalized functions: take any v £ 5 ( R n ) , multiply Eq. (15.41) by 6(£),
and integrate w.r.t. f over R n . By using Eq. (15.34)-Eq. (15.38) and the
Lebesgue theorem, it is not difficult to show that
= (9o,v).
Under additional conditions on the data, it is possible to show that /(£, x) ->
go(x) as t —¥ +0, a. e.
To prove that Eq. (15.37) holds in even stronger sense, and make more
exact statements on the regularity of / , the following method is more con-
venient.
Asu = a(Dx)u.
We can interpret the problem Eq. (15.36)-Eq. (15.37) as the Cauchy prob-
lem for the ordinary differential equation on R + with the operator-valued
320 Pseudo-differential operators with constant symbols
coefficient. For the sake of brevity, we consider and solve it in the case
9 = 0:
f'(t)+A*f(t) = 0, t>0, (15.47)
/ ( + 0 ) = go- (15.48)
Theorem 15.17 For any sell andg0 G Hs(Rn), the problem Eq. (15.47)-
Eq. (15.48) has a solution f € C°([0,+oo);Hs(R+)).
The solution is unique, and for any to > 0,r € R , j € Z+, and 7 <
Co, where CQ > 0 is the constant in Eq. (15.38), it satisfies the following
estimate:
sup||e77(j)(t)||r<+oo. (15.49)
t>t0
• Tt = exp(-L4), t > 0;
• Tt+r = TtTTfort,r>0;
• for any t > 0 and go € B, lim T _^ TTgo = Ttgo in the topology of B;
To deduce Theorem 15.17 from Theorem 15.18, we need two auxiliary re-
sults.
Lemma 15.9 There exist Co > 0 and 9 G (n/2,n) such that if (A,£) G
E# x R , then
Proof. If c > 0 is sufficiently small and (0™ < c|A|, then from Eq. (15.34)
and the triangle inequality we deduce
|A + a ( 0 | > ( | A | + ( 0 m ) / 2 ,
and Eq. (15.52) holds for these (f,A). Fix c > 0; then from Eq. (15.38),
there exists 9 e (w/2,ir) such that if |A| < c - 1 ( £ ) m and A € £<?, then
SRA > — 5fta(£)/2, and hence for some ci,C2 > 0,
lA + a ^ l ^ c x ^ r ^ c a d A l + ^ n .
From Theorem 15.4 and Lemma 15.9, we conclude that there exist C\ > 0
and 9 6 (TT/2, it) such that if A e Tig, then Eq. (15.50) hold. Hence, Theorem
15.18 is applicable, and the solution to the problem Eq. (15.47)-Eq. (15.48)
is given by Eq. (15.51) with A = As.
By making the Fourier transform w.r.t. x, then changing the variable
A = in, and deforming the contour {n € C | in G Cg} into the real line
$sn = 0, we obtain that Eq. (15.51) gives the solution Eq. (15.43) to the
problem Eq. (15.36)-Eq. (15.37) with g = 0. In particular, for any go G
Hs(Rn), f[t) = exp[-tAs]g0 is of the class C 0 ( [ 0 , + o o ) ; # s ( R " ) ) , and /
satisfies the initial condition Eq. (15.48), where the limit is understood in
the sense of the topology of i f s ( R n ) .
It remains to prove the estimate Eq. (15.49). First, notice that if
we make the change of the unknown f(t) = e _ 7 t / 7 ( t ) in Eq. (15.47)-
322 Pseudo-differential operators with constant symbols
/;(*) + ( A s - 7 ) / 7 ( t ) = 0, (15.53)
A(0) = 5o- (15.54)
Due to Eq. (15.38) and the condition 7 < Co, the problem Eq. (15.53)—
Eq. (15.54) enjoys the same properties as the problem Eq. (15.47)-Eq. (15.48),
and estimates Eq. (15.49) for / are equivalent to estimates Eq. (15.49) with
7 = 0 for / T . Hence, it suffices to prove Eq. (15.49) with 7 = 0.
By differentiating under the integral sign in Eq. (15.51) and using Eq. (15.52),
we obtain
+00
(l + l a D ' V ^ d a l l s o l U .
/ -00
For any to > 0, the integral converges uniformly in t > to, since 8 >
7r/2, and we deduce that Eq. (15.49) holds with r = s. By applying to
Eq. (15.47) a(D)-1, which maps Hs(Rn) to Hs+m(Rn) due to Eq. (15.42)
and Theorem 15.4, we find that
f(t) = -a{D)-lf'(t)
satisfies Eq. (15.49) with r = s+m. By induction, we obtain that Eq. (15.49)
holds with any r, and Theorem 15.17 is proved.
Lemma 15.10 Let a satisfy the estimates Eq. (15.23) and Eq. (15.38).
Then there exist 8 € (ir/2,n) and Cj,j > 0, such that for all A G £#, and
fell"
£ \D^((\+aior'my <cj(i+ 1X1)^, j = o,i,.... (15.55)
Parabolic equations 323
Proof. For j = 0, this is Lemma 15.9, and Eq. (15.55) for j > 1 can be
easily proved by induction, the Leibnitz rule, Eq. (15.50) and Eq. (15.23)
being used. D
Lemma 15.11 Let s,m G R, and let a satisfy the estimates Eq. (15.23)
and Eq. (15.38). Then there exist 6 G (7r/2,7r) and C such that for all
AGE*,
and by replacing in the proof of Theorem 15.17 the scale of the Sobolev
spaces with that of the Holder-Zygmund spaces, we obtain the Holder-
Zygmund analogue of Theorem 15.17.
Theorem 15.19 Let a satisfy the estimates Eq. (15.23) and Eq. (15.38).
Then for anysGK and g0 G C s ( R n ) , the problem Eq. (15.47)-Eq. (1548)
has a solution f G C°([0,+oo);C s (R n )).
The solution is unique. It is given by the RHS of Eq. (15.51), and
for any to > 0, r, j G Z + , and 7 < CQ, where CQ > 0 is the constant in
Eq. (15.38), the solution satisfies the following estimate:
Notice that due to the remark after Theorem 15.10, the results of this
subsection are valid for any other scale of Besov spaces.
Theorem 15.20 Let a satisfy the estimates Eq. (15.23) and Eq. (15.38).
Then for any s € R, and functions g 6 C°([0,+oo);C s (R n )) and go e
C*(R n ), the following statements hold:
a) the problem Eq. (15.58)-Eq. (15.59) has a unique solution of the class
C°([0,+oo);C s (R n )), which is given by
f(t) = fi + h, (15.61)
where
f2 = exp[-tAs](g0 - A);
d) for any i 0 > 0, positive non-integer r, j 6 Z+, and 7 < CQ, where CQ > 0
is the constant in Eq. (15.38), fo satisfies the estimate Eq. (15.57).
Proof. a)-b) Direct verification; the details are left to the reader as an
exercise.
c) Set f2(t) = f{t) - A . The problem Eq. (15.58)-Eq. (15.59) becomes
the problem Eq. (15.47)-Eq. (15.48) for f2 with g0 - fx £ C s (R n ) on the
RHS of Eq. (15.48). After solving this problem and applying Eq. (15.57),
we obtain c). D
p+a{D)u = g. (15.66)
This is the Wiener-Hopf equation. In order that the problem Eq. (15.66)
be well-posed, the symbol a must satisfy some regularity conditions.
(i) a e 5 m ( R ) ;
(ii) there exists C > 0 such that
(iii) there exist non-zero complex constants d± and p\ £ (0,1] such that
as £ —> ±oo,
(iv) there exist C > 0 and /92 € (0,1] such that for all £ € R,
j=0
be a polynomial of degree m, with no real zeroes. Clearly, conditions
Eq. (15.67)-Eq. (15.69) are satisfied, and if {fj}i<j<ro is the set of zeroes
of a, then we may set
M O = am I I (£ - &)>
a-(o = n «-&•)•
In the general case, the construction is rather lengthy, so the reader must
be patient.
Choose d € C \ {0} and K0 € C so that the function
B 0 (O := d - J ( l + i O " ° ( l - i O _ " ° a ( 0 ( l + C 2 ) _ m / 2
satisfies
lim B 0 ( 0 = 1- (15-71)
£-+±oo
By taking the logarithm of B o ( 0 and using Eq. (15.68) and Eq. (15.71),
we obtain that d and KQ solve the system
ln(d±/d) ± i7TKo = 0.
(We have chosen the branch of the logarithm so that In a € R for a positive
a.) We see that d must be equal to
and therefore
K0 = iln(d+/d_)/27r. (15.73)
Due to Eq. (15.67) and Eq. (15.69), Bo is a continuous function, which does
not vanish on R, hence from Eq. (15.71), we conclude that the winding
number of the curve {Bo(0 I £ £ ( _ 0 ° , +oo)} around zero
+oo
dargB0(0 (15.74)
/
-oo
The Wiener-Hopf equation on a half-line I 327
is an integer. Set
hence
6-(£-ia) = (27TZ) _1
/" °° 7 - ^ dry. (15.82)
J-bo i - w -i)
Due to Eq. (15.80), the integrals in Eq. (15.81)-Eq. (15.82) converge. Set
where t+ admits the following estimate: for any e > 0, there exists Ce > 0
SMC/I i/iat /or a// £ in i/ie half-plane 9 £ > 0,
o_(C) ±1 = ( l + i O ± K - + * - ( 0 , (15-91)
where tt admits the following estimate: for any e > 0, there exists C( > 0
SMC/I that for all £ in the half-plane 3£ < 0,
Proof. First, we prove (o+) - {d+); (o_) - (d_) are proved similarly.
Since the function d(—1)'(1 — i£)K+ satisfies (a+) — (d+), it follows from
Eq. (15.83), that it suffices to prove the corresponding properties of b+,
namely
b+ (£) is holomorphic in the upper half-plane 9£ > 0, admits the continuous
extension into the closed half-plane and satisfies estimates Eq. (15.87) and
Eq. (15.88) with K+ = 0.
By differentiating under the integral sign in Eq. (15.81) and using Eq. (15.80),
we obtain, for £ in the half-plane 9£ > T > 0 and j = 1,2,...:
-t-oo
/ (l + M)-' 1 K-'?r 1 - i *J- (15-94)
-oo
From the triangle inequality, for rj e R and the same £,
(i+iiir'^i.ir'i^-i (15.95)
therefore for any UJ 6 [0, p^, the RHS in Eq. (15.94) admits an estimate via
+oo
le-ur1-^"^. (15.96)
/
-oo
We may take any ui < 1, u € [0, pi}\ with this choice, for j > 1, the integral
in Eq. (15.96) converges, and we obtain that b+(£) satisfies Eq. (15.87) with
K+ = 0 .
Now we prove Eq. (15.88) for b+ with K+ = 0, in the half-plane 3f > 0.
Fix such £, and introduce the following subsets of R:
MO = { » j e R | | e - f j | > i , lfl<M/2},
MO = {TjeR||£-»/|>l, \£\>2\r,\},
MO = {77 G R | K - rj| > 1, N / 2 < |^| < 2|TJ|},
Then for j - 1,2,3,4 and 3£ > 0, define b+tj(0 by Eq. (15.81) with the
integration over Jj(£). Clearly,
4
330 Pseudo-differential operators with constant symbols
= taln^1)"^' j
<r->+0 _(a-2+1)1/2 +
= in, (15.105)
we conclude that b+ i4 (£) satisfies the estimate Eq. (15.99) with p = minjpi,
P2} instead of pi.
We have shown that b+j ( 0 , .7 = 1,2,3,4, satisfy the estimate Eq. (15.99)
with p instead of p\, hence &+(£) also does, and to finish the proof of
(a + ) - (d + ), it remains to show that b+ admits the continuous extension
up to the boundary 3£ = 0. We have proved the uniform boundedness of
b+j(£), as Qf£ —>• +0, 6 + ^ ( 0 being represented as a sum of two terras. Of
those, the first one has the limit
-{2-Ki)-1b(C)i^ = 6(0/2
as a —¥ +0; by the inspection of the proof, it is seen that the integrands in
the integrals defining the second term and 6+,j(f) for j = 1,2,3 have limits
as a —¥ +0 and £ remains fixed. Prom the Lebesgue theorem, all these
integrals also have limits as a ->• +0. This finishes the proof of (a + ) - (d+).
We see that Eq. (15.70) follows from Eq. (15.76) and the following equal-
ity:
6±(C ± u r ) - • 6(C)/2 ± 6 ( 0 ,
where b is the same for both signs ± , and Eq. (15.106) is proved.
It remains to prove part (f). Suppose, there exists another pair o'±
satisfying (o±) - (b±) and Eq. (15.70). Then on the real axis
<(O/MO = M0/a'-(0,
and the LHS (resp., the RHS) admits a bounded holomorphic extension
into the upper (resp., lower) half-plane. Thus, we have a bounded holo-
332 Pseudo-differential operators with constant symbols
p+a+(D)a+{D)-1lg - p+a+(D)e-a+(D)-1lg = g.
The Wiener-Hopf equation on a half-line I 333
0 s—RK_
The second term on the LHS is equal to zero since 0_a + (.D) -1 Zg £H
O 5 — 771
g M- a-(D)-H+a+(D)-Hg
where \s' — JJ/s_| < 1/2, find a unique solution in H (R+) DH(R+), and
after that check whether the solution defined by Eq. (15.109) belongs to
H(R.+)- If it does, we have a unique solution of the Wiener-Hopf equation
Eq. (15.66), of the class ^ ( R + ) , and if it does not, Eq. (15.66) has no
solutions of this class. In the latter case, one may introduce into Eq. (15.66)
potential operators in order to obtain a solvable problem (see, e.g., Eskin
(1973), Section 8), but we will not need this construction here.
Case III. s — 3t/c_ = —q+s', where q is a positive integer and \s'\ < 1/2.
Due to Theorem 15.4, A.-(D)-qa+{D)-1lg e # S ' ( R ) , hence we can
apply Lemma 15.7:
By inserting Eq. (15.113) into Eq. (15.108) and rearranging terms, we ob-
tain
a_(D)u - A _ ( £ » ) ^ + A _ ( £ » ) - ? a + ( ^ ) - 1 ^ (15.114)
, 1
= A_(D)^_A_(D)-'a+(£ )" ^ + 9i,-
a_(D)u - A_(D)«6»+A_(D)-«a+(I»)- 1 Zg = ^ T c ^ - 1
^.
u = a^(D)-1A-(D)"e+A^(D)-qa+(D)-llg + a^(D)'1^2cjDj-160.
j=i
(15.115)
One of standard ways to obtain a well-posed problem is to add to Eq. (15.66)
appropriate boundary conditions (see, e.g., Eskin (1973), Section 8); for
our purposes, a more appropriate way to choose a unique solution is a
specification of the behaviour of the solutions to the Wiener-Hopf equation
Eq. (15.66), near the boundary. Possible types of the behaviour are studied
in the next Subsection.
K_ G ( - 1 / 2 , 3 / 2 ) ; pj+K-> 1. (15.116)
u = ca^(D)-1S0+c1a-(D)-1AZ1(D)S0 (15.121)
1 1 1
+a-{D)- AZ (D)9+A-(D)a+(D)- lg,
where
O »K_+p-l/2-£
vi = ctZ(D)S0eH (R+), V e > 0 ;
0»K_+l/2-£
v2 = cia-{D)-1AZ1{D)S0eH (R+), V e > 0;
v3 = a-(D)-1AZ1(D)6+A-(D)a+(D)-1lgeH K +12
~ (R+).
By using Theorem 15.9 and Eq. (15.116), we conclude that for any e > 0,
vi € C H K - + p _ 1 _ e ( R ) , and v2,v3 e C R K - _ £ ( R ) . Since they are supported
on [0, +oo) and p < 1, we conclude that Vj(x) = 0(xUK-+p~1"t), for
any e > 0. By applying Eq. (15.124) with v = K_ to the first term in
Eq. (15.122), we obtain Eq. (15.117).
This proves a), and to prove b), we set c = 0 in Eq. (15.121), and apply
Eq. (15.91) to the first remaining term on the RHS:
where
0K-+p+l/2-e
vn = c1(lg)tZ(D)AZ1(D)80 £H (R+), V e > 0,
The Wiener-Hopf equation on a half-line I 337
and v3 is the same as in the proof of a). The same argument as above shows
that
and by applying Eq. (15.124) with v = K _ + 1 to the first term in Eq. (15.125),
we obtain Eq. (15.118).
II. If 5R/c_ = 1/2, take s < 0 such that s - /c_ £ ( - 3 / 2 , - 1 / 2 ) ; then the
same argument as in the part I shows that the solution u £H(R+) admits
either the representation Eq. (15.117) or Eq. (15.118). Since x+ ' is not
square integrable in the neighbourhood of zero, Eq. (15.117) is impossible.
It remains to consider the case 3JK_ £ ( - 1 / 2 , 1 / 2 ) , when the solution is
unique and given by Eq. (15.109). By applying Eq. (15.31) to Eq. (15.109)
and arguing as in the part I, we arrive at Eq. (15.118). •
In the following corollary of Theorem 15.24, we list cases when the bounded
and continuous solution exists and it is unique.
d) if 1 < K_ < 3/2, there exists a unique u with bounded u'; it is given
by Eq. (15.118) and belongs to &(&).
Proof. Note that the behaviour of u' near the origin is the same as the
behaviour of (1 + iD)u = u + u', apply (1 + iD) to Eq. (15.121) and argue
exactly as in the proofs of Theorem 15.24 and Theorem 15.25. •
15.8 P a r a b o l i c e q u a t i o n s o n [0,T] x R +
u{x,t) = eCT0'u<r°(x,f),
A^(A,0 = (l + | A | 1 / m ± i O s , (15.135)
since 5(0,£) = B(£) satisfies Eq. (15.79). Therefore, the winding number
around zero of the curve {B(X, £) | - oo < £ < +00} is independent of
A e S 9 ; on the strength of Eq. (15.78),
+00
dargfl(A,0=0. (15.139)
/
-00
By using Eq. (15.137) and Eq. (15.139), we conclude that b = \nB is well-
defined on Hg x R.
340 Pseudo-differential operators with constant symbols
B(A,0 = l + O ( ( 0 - p i ) , (15.141)
( c ) / o r a / i ( A , 0 G £« x R,
Proof. Parts (a±), (c) and (d) are immediate from parts (a±),(e) and
(/) of Theorem 15.21, and the proof of {b±) is a modification of the proof
of the same parts of Theorem 15.21. As in the proof of Theorem 15.21, it
suffices to show that 6± are uniformly bounded: there exists C such that
for all A G £# and all £ in an open half-plane ± 3 ^ > 0,
\b±(\,0\<C (15.149)
We prove Eq. (15.149) for the sign "+"; the proof for the sign "-" is essen-
tially the same. Let A e S« and ^ > 0. Set K(X) = (1 + \X\)1/m, and
introduce the following subsets of R:
Then for j = 1,2,3,4 and Sf > 0, define b+J{X, f) by Eq. (15.142) with the
integration over Jj(X, Q. Prom Eq. (15.34), and Eq. (15.69) and Eq. (15.38),
we conclude that B(X, f) is uniformly bounded and bounded away from zero
on £9 x R:
3J/c±e[0,m]. (15.154)
p+a+(D)0+a+(D)-1lg = p+lg = g
(for more details, see the proof of Theorem 15.22), hence the kernel of
{A")'1 is trivial, and we can introduce A" as the left inverse to ( A 8 ) - 1 .
Thus, A3 is the operator in H"~m(R+) with the domain
and
and
0 s—m 0 s—m
(l + I A D ^ - ^ a - t A , ! ) ) - 1 ^ (R+)-^fl- (R+) (15.160)
With the estimate Eq. (15.157) at hand, we can use Theorem 15.18 with
B = Hs_m(R+), define Tt = exp[-L4 s ] by Eq. (15.51) and obtain analogues
of Theorem 15.19 and Theorem 15.20. We state the inhomogeneous version.
/ G C1((0,+oo);tfs-m(R+)) nC°([0,+oo);^1/2+SK--€(R+));
d) for any to > 0, j £ Z + , and 7 < CQ, where CQ > 0 is the constant in
Eq. (15.38), fi satisfies the estimate
In this Section, real A± and u± are assumed to satisfy A_ < 0 < A+ and
w_ < 0 < UJ+.
Definition 15.3 «S(R; [A_, A+]) denotes the space of C°°-functions such
that for any 7 G [A_, A + ], elxu G <S(R). The topology in <S(R; [A_, A+]) is
PDO in the Sobolev spaces with exponential weights, in ID 345
is a topological isomorphism;
f) Co°(Ii) is embedded into S(R; [A_, A+]), densely and continuously.
9) If X- < A_ < A+ < A'+, then <S(R; [A'_, A'+]) C S(R;[A_,A+]),
densely and continuously; in particular, S(R; [X'_, X'+]) C «S(R), densely
and continuously;
346 Pseudo-differential operators with constant symbols
h) If\'_ < A_ < A+ < X'+, thenS'(R; [-X+, -A_]) C S'(R; [-X'+, -X'_]),
densely and continuously.
Proof, a), b) and d) are immediate from the definition, and c) and e)
follow from b) and d), respectively, by duality.
f) follows from b) and Theorem 15.3, part c).
g) Observe that
a) u(£) admits the analytic continuation into the strip Q£ G [A_,A+], and
the continuous extension up to the boundary;
b) for each r e [A_,A + ], u{- + IT) 6 <S(R), and u 6 C°°([A_, A + ];5(R)).
Definition 15.5 For u G <S'(R; [—A+, — A_]) and 7 G [—A+, —A_], define
u(- + n) eS'(R) by
IMI.,7 = IKII».
where || • || s is the norm in the Sobolev space HS(R).
is finite.
I*i7 a r +oo+i7
+00+17
-00+17
-00+17
(l + k | 2 ) s K 0 | 2 ^
\N 1X '/ *2
/
(15.163)
PDO in the Sobolev spaces with exponential weights, in ID 347
-t-00
e _il "u 7 (a;)da;
/ -00
/+00
e- ia "'e' 1 ' x u(x)dx
-00
/+00
e - "*u(a;)dx,
-00
which is Eq. (15.2), the only difference being that in Section 15.2, £ is
assumed real.
Example 15.8 Let u{x) = 9+(x)e0x. Then the Fourier transform u(£)
is well-defined for 3 £ < -K/3:
f+co
u(£) = / e-^eP'dx = - ( 0 - tO"1,
JO
and it follows that u G H S ' T (R), if and only if s < 1/2 and 7 < -3fy?.
For s 6 R, and A > 0, define \'X{D) as the PDO with the symbol (1 + A2 +
£2y/2
Lemma 15.14 Let m,s G R, p G R, 7 G [A_,A+] and A > max{|A_ —
p|> •*+-/>}•
Tften
is a topological isomorphism.
By Lemma 15.3, Asx(D + iX±) is continuous in <5>(R), which implies that the
RHS in Eq. (15.165) admits an estimate via a linear combination of terms
of the form Mx±-k',N';o- D
Corollary 15.2 Let conditions of Lemma 15.14 hold, and let
be defined by duality: foru € «S'(R; [—A++p, —A_+p]) andv € <5(R; [A_, A+]),
(e"xAUD)u,v) = (u,A{(D)e"xv).
Proof. By definition,
-t-oo
/
(n)2su(n + »7)ii(fj + n)dn; (15.168)
-OO
S
c) C0°°(R) C <S(R;[A_,A+]) C H (R) C <S'(R; [-A+, -A_]) densely, and
these embeddings are continuous;
d) i f s ' 7 ( R ) C i f s ' 7 ( R ) , densely and continuously;
e) i ? - s ' ~ 7 ( R ) can be naturally identified with the dual space (if s ' 7 (R))*,
by
+00+I7
w(t)v(Z)<%- (15.169)
/
-00+17
Proof, a) Evident.
b) We can define (u,v)s,-y = (u 7 ,w 7 ) s , and then
+00
/
(Ti)2s^(T))€y (n)dTj
-00
+00
/
{T})2SU(TJ + 17)1/(77 + i~f)dn.
-00
c) By using topological isomorphisms defined by the multiplication-by-
e7X operator, we reduce to the case s = 0. After that we apply part c) of
Theorem 15.3 and parts f)-h) of Lemma 15.13.
d) Lemma 15.15 with s = 0 and part d) of Lemma 15.13 reduce to the
part d) of Theorem 15.3.
e) Set A = |7|, and introduce an equivalent norm in HS(R) by
H i ; = ||AX(D)«||0.
The corresponding norm in i ? s ' 7 (R) is
||u||;, 7 = \\Asx(D)u,\\0,
(u,v)' = {Als(D)uy,Vl)0.
350 Pseudo-differential operators with constant symbols
Since # * ' 7 ( R ) is the Hilbert space, any <j> G (if s ' 7 (R))* can be denned by
some u G # S ' 7 ( R ) :
(w,v) = (iy_7,i;7)
+oo
/
-oo
+00-M7
/ w(S)H0dt.
-00+17 •
The following lemma and theorem are evident corollaries of Lemma 15.4
and Theorem 15.6, respectively.
Lemma 15.16 Let s > 0 be an integer, and 7 G R. Then an equivalent
norm in ff s ' 7 (R) can be defined by
p':C0°°(R)9uH«(0)6C
PDO in the Sobolev spaces with exponential weights, in ID 351
by
a(5)«(0 = a(Ofi(0, ^ = 7;
if a(- + ry)u(- + 17) G I a ( R ) , we can use the explicit formula
+oo+i7
eix^a(0u(0dx. (15.171)
/
-00+17
Theorem 15.31 Let m,s G R and 7 G [a;_,w+], and let a G S m ( R +
i[u-,uj+]). Then the operator
|a(0-1|<C<fl-m, 9£ = 7, (15-173)
d±=lima(0|ei"m, (15.174)
Corollary 15.3 Let a satisfy conditions of Lemma 15.17, and let Z(a)
be the set defined in Lemma 15.17.
Then for any 7 G (CJ_, oj+)\Z(a), the operator Eq. (15.172) is invertible.
PDO in the Sobolev spaces with exponential weights, in ID 353
Theorem 15.34 Let s > 1/2, and 7 G R. Then for any s1 < s - 1/2,
i f s ' 7 ( R ) C C ' 7 ( R ) , densely and continuously.
Theorem 15.35 Let m G R, and let a G C°°(R + i[w_, w+]) satisfy, for
any multi-index a,
\a{a)(0\<Ca(Om~M, (15.175)
where Ca is independent of £ € Ti. + i[u!-,u)+].
Then for any s G R and 7 G [u-,ui+], A = a(D) : C ' T ( R ) -»
s m 7
C ~ ' (R) is bounded, with the norm admitting an estimate
SU
\\a(D)\\<C^2 P «>-m+|a,|o(o)(OI, (15-176)
|a|<JV«6R+i7
Theorem 15.36 Letm G R, and leta satisfy Eq. (15.173) and Eq. (15.176).
Then for any s,7 G R, a(D) : C S ' 7 (R) -¥ C s - m ' 7 ( R ) is invertible.
is a topological isomorphism.
Proof. Clearly, A™± satisfies the condition of Theorem 15.38, for any m,
hence the operator Eq. (15.179) and its counterpart with —m instead of
m are continuous. Since A™±(£>)A^™(D) = I and the same holds for the
product in the reverse order, Eq. (15.179) is the topological isomorphism^
Theorem 15.39 Let s' < s and 7' < 7. Then
o».7 o*'.7' o S>V o»'>7
H (R+) CH (R+) and H ( R - ) CH (R-),
densely and continuously.
Proof. With 7 = 7 ' , this is an evident corollary of Eq. (15.29) d), there-
fore it suffices to consider the case s = s'. Take A > I7I, and use Lemma
15.19 to reduce to the case s = 0, i.e., to the statements:
0
+f = £ fi,-K!'-(D)S0 + A*™ (£>)*+A£_(D)/, (15.183)
k=i
where fjt+~p'Akx^(D)f€C.
Theorem 15.44 (cf. Theorem 15.16) Let 7 G R and s = —m + s', where
m > 0 is an integer and \s'\ < 1/2.
0 ».7 o s.7
T/ien any f £H (R+)n H (R~) admits a decomposition
m
/ = £/tI>*-%, (15.185)
jfe=i
358 Pseudo-differential operators with constant symbols
where fk € C .
This is the Wiener-Hopf equation. In order that the problem Eq. (15.190)
be well-posed, the symbol a must satisfy some regularity conditions.
Definition 15.11 We say that a 6 5™ g (R + i[A_,A + ]) and call it a
regular symbol of order m, if the following conditions hold:
(i) a e S m ( R + i[A_,A+]);
(ii) there exist non-zero complex constants d± and p\ e (0,1] such that
\a'(0\<C(0m'p2- (15.192)
Note that a 7i _ (resp., a 7 ,_) may have zeroes in the strip 9£ € (7>A+]
(resp., 9 { e L , 7 ) ) .
360 Pseudo-differential operators with constant symbols
(c_) for any r < A+, e > 0 and j = l,2,..., there exists CT,e,j such that
for all f in the half-plane Q£ < T,
where
O »iT
Then u €H (R-+), the solution to Eq. (15.190), admits a representation
« = £ cj(lg)A^AD)a1AD)-16o+A^_a^(D)-1e+Al_(D)a^+{D)-1lg,
2=1
(15.208)
where A = max{-A_,A+}, and
1
u= a 7 ; _(D)- A A ,_(Z))^ + A A ,_(D)-«a 7 ) + (Z>)- 1 Z f f +a 7 ,_(Z3)- 1 ^ c , ^ - 1
^.
3=1
(15.210)
growing data]
Parabolic equations o n R x R.+ with exponentially growing data 363
Assume that m > 0, a € S™g(R + i[A_, A+]), and Eq. (15.191) holds with
Kd± > 0, and consider the problem
f(x,t) = e°°t-^f(ro„(x,t)
g(x,t) = e'lt-^g^yfat),
lx
go(x) = e~ gon(x).
In terms of /<r0l7 and <7<70)7, the problem Eq. (15.212)-Eq. (15.37) becomes
account that aVo^ e S^nreg(R), we can use the results of Section 15.8 to
solve the problem Eq. (15.216)-Eq. (15.218), and then obtain the solution
to the problem Eq. (15.212)-Eq. (15.214). The details are left to the reader.
Chapter 16
Elements of calculus of
pseudodifferential operators
In this Chapter, we Ust basic definitions and facts for simple classes ofPDO
in R n , which are necessary for the proofs of the results on NIG-like Feller
processes in Chapter 14, and calculation of approximate prices of European
options in Chapter 14 and basket options and exchange options in Chapter
9. The results of this Chapter can be used to calculate prices of barrier
options under Levy-hke Feller processes.
In contrast to the case ofPDO with constant symbols, where the Fourier
transform establishes the equivalence of the action of PDO and the mul-
tiplication operator, and so rather weak regularity conditions on symbols
suffice, here appropriate estimates on the derivatives of symbols must be
imposed. For simplicity, one usually assumes that symbols are infinitely
smooth, and imposes conditions on all derivatives, though weaker condi-
tions suffice to obtain the main results. In Section 16.1, we explain the
main ingredients and structure of the general theory of PDO for the sim-
plest class. We omit the most technical proofs like the composition theorem
and the theorem on the boundedness in appropriate scales ofBanach spaces
but we include several proofs, which demonstrate how the theory of PDO
works in relatively simple situations.
In Section 16.2, PDO depending on a parameter are considered; this
variant of the theory of PDO is needed in order to obtain approximate so-
lutions to pseudo-differential equations, when the exact analytical formulas
are non-avaiiable, and in Section 16.3, we consider PDO with symbols holo-
morphic in a tube domain; this variant of the theory of PDO is necessary to
treat problems with exponentially growing data. In Section 16.4, the proofs
of the technical results of Chapter 14 are presented. The transformation of
365
366 Elements of calculus of pseudodifferential operators
We start with this by now quite classical case since results for PDO with
symbols of the class 5 m ( R n x ( R n + iU)) can be deduced quite easily
from the results for S^0(Rn x R n ) . Most of the results below, with the
exception of the boundedness theorem for the Holder spaces, can be found
Grubb (1996) and Hormander (1985) (and local versions in Shubin (1978));
however, the form and proofs of some results, which we give below, do not
necessarily coincide with the proofs in these monographs.
if for any N
N-l
rN:=a~Y, «j £ S^N{Rn x R"). (16.3)
o=o
Basics of the theory of PDO with symbols of the class S™0(R.n x R n ) 367
aj(x,t£)=tm-jaj(x,t),
(i) there exists a £ 5j^ 0 (R n x R n ) such that Eq. (16.2) holds, and
(ii) if a1 £ Sj™0(Rn x R") is another symbol satisfying Eq. (16.2), then
a1-a<E5-°°(RnxRn).
\a\<m
D)u{x)= J2 aa(x)(Dau)(x).
|a|<m
368 Elements of calculus of pseudodifferential operators
ar(x,D)u= ^2 (Da(aau))(x).
\a\<m
/ -OO J —CO
where the integral is understood as the iterated one. PDO with the right
symbols naturally arise when one considers the adjoint operators. The
following theorem ensures that the adjoint operator belongs to the same
class of PDO, and gives the formula for its symbol. Its proof and the proof
of the Composition Theorem are rather technical, and we refer the reader
to Shubin (1978), Hormander (1985) and Grubb (1996).
Theorem 16.3 Let a e ^ ( R " x R " ) .
Then the (complex) adjoint a(x,D)* defined by
(a(x,D)u,v) = (u,a(x,D)*v)
a*(x,S)~ X>!rXa£i(:r,0.
|o|>0
Theorem 16.3 and Theorem 16.2 allow one to define the action of a{x, D)
on <S'(Rn) by duality:
Clearly, if <f> € <S(R") then the new definition of a(x, D)<f> coincides with the
old one. Further, a(x, D) is a continuous map from <S'(Rn) into itself. This
map can be defined as a unique continuous extension of the map defined by
Eq. (16.4) on 5 ( R n ) , and the standard agreement in the theory of PDO is
to denote this map by the same symbol a(x, D); the restriction of a(x, D)
on any subspace L C <S'(Rn) is also denoted by the same symbol.
1
In this context, Eq. (16.4) is the definition of the PDO a(x,D) with the left symbol a.
Basics of the theory of PDO with symbols of the class S J " 0 ( R " x R n ) 369
The next basic fact and the most important tool in the theory is the
Composition Theorem.
T h e o r e m 16.4 Let a G S^ 0 (R n x R n ) and b G S^(Rn x Rn).
Then C = a(x, D)b(x, D) is a PDO with the symbol c € S™ 0 +m '(R n x R " ) ,
which admits an asymptotic expansion
c(x,0~ £ (a!)"1^(x,06(Q)(x,0, ( 16 - 5 )
|a|>0
in the sense that for any integer N > 0,
l a ^ r ^ C ^ r ™ . (16-7)
where C is independent of (x,£).
We call a an elliptic symbol, and A = a(x, D) an elliptic operator.
R e m a r k 16.2. The standard definition of an elliptic symbol requires
Eq. (16.7) outside a compact set in £-space. We use Eq. (16.7), since it
suffices for our purposes and admits a trivial generalization for the case of
operators depending on a parameter.
The following definition formalizes an intuitive notion of an approximate
inverse.
Definition 16.4 We say that B — b(x, D) is a parametrix of A — a(x, D)
if and only if
b ~ 6_ m + 6_ m _! + • • •, (16.9)
370 Elements of calculus of pseudodifferential operators
where
i4A(x,D) = / + t ( i ) i ? ) , (16.12)
«~$3*., (16-14)
S>1
j+\a{=s
where
&m/a(*,0 = a ( * , 0 1 / 2 : = exp[-Ino(a;,Pl-
under assumption Eq. (16.16), the symbol bm/2 is well-defined, and by using
Eq. (16.7), one can easily check that it belongs to S™Q2(Rn x R n ) . Fix any
sequence of 6 m / 2 _ s G S™ 0 /2_S (R" x R n ) , s > 1, define b G S™ 0 /2 (R n x R»)
by Eq. (16.18), by using Theorem 16.1, and calculate the symbol of b(x, D)2
by using the composition theorem. We obtain that it admits an asymptotic
expansion
I M I o = J2 \\Dau\\c...
\a\<m
Theorem 16.8 Let a G I S^0(Rn xRn), and let for some s,A = a(x,D) :
Hs{Rn) -» H3-m(Rn) be invertible.
Then the inverse is a PDO of the class S^™(R n x R n ) , and its symbol
admits the same asymptotic expansion Eq. (16.9) as the parametrix.
Proof. <S(Rn) C Hr(Rn), for any r, hence both A and A'1 maps <S(Rn)
into itself. Therefore, they are mutual inverses as operators in 5 ( R n ) . Since
5 ( R n ) C Hr{Rn) densely, and A " 1 : Hr-m{Rn) -> Hr{Rn) is bounded,
we conclude that A' is the inverse to A : H (R ) -> F r _ m ( R n ) .
1 r n
•
Remark 16.3. Corollary 16.1 is valid not for the scale H"(Rn) only but
for many other scales of spaces as well.
Theorem 16.9 Let the conditions of Theorem 16.6 hold, and let $la(x, D)
be positive-definite, and%a(x,D)($la(x,D))~1 be bounded.
Then there exists a PDO B = b(x, D) such that B2 = A, and its symbol
b e 5^o ( R n x R " ) admits the same asymptotic expansion Eq. (16.18) as
an approximate square root.
call it d, which describes the smallness of the derivatives w.r.t. £. The case
when only derivatives w.r.t. x (respectively £) are small obtains by fixing
d (respectively e).
So, let £ C (0,1] and V C [1, +oo) be the sets on which e and d live; we
assume that the statement e/d —• 0 makes sense on the set £ x T>, and we
regard a pair (e, d) as a parameter.
The starting point is a natural generalization of the class S™ 0 (R n x R n ) .
S-°°((£ x V);Rn x R n ) := n m € R S f t ( ( £ x D ) ; R n x R n ) .
instead of £ x P .
N-l
rN:=a-Y, <*j G e ^ S J ^ f x XJ; R n x R " ) . (16.25)
j=o
|a(e,d,s,0-1|<C,(d2 + |e|2)-m/2.
where C is independent of (e, d, x, £) with e/d < eo-
We call A — a(e,d; x, D) an elliptic operator with a parameter.
6_ m = l / o € S^J*((f x 2J) £ l ;R n x R n ) ,
and
and
\\(d/e)t2(x,D)\\c--tc-<C. (16.29)
When e/d < ei := min{e 0 ,1/C}, we conclude from Eq. (16.27)-Eq. (16.29)
that AB° and B°A are invertible in C s - m ( R n ) and C s ( R n ) , respectively.
By the boundedness theorem, B° : C3~m(Rn) -+ C*(R n ) is bounded, and
hence, B°{1+ ti(x,D))~1 and (1 + t2{x,D))-1B° are bounded right and
left inverses to A : C"(R ) — C s - m ( R n ) .
n
Remark 16.4. Theorem 16.12 is valid for other spaces as well; for instance,
in the scale of Sobolev spaces H3(Rn), it is valid without the restrictions
on s and m.
andb admits the asymptotic expansion Eq. (16.18) in the sense of Definition
16.6, with the terms given by Eq. (16.19)-Eq. (16.20).
Proof. Exactly the same argument as in the proof of Theorem 16.6 allows
us to construct b G S^0((£ x V)eo\Rn x R"), which satisfies Eq. (16.17)
with t G S~°°{{£ x V)e0;Rn x R n ) ; Eq. (16.18) is understood in the sense
of Definition 16.6, and all the terms are given by Eq. (16.19)-Eq. (16.20).
By using a straightforward modification of the results on fractional powers
in Seeley (1967), Shubin (1978) and Grubb (1996), we obtain that if e0 > 0
is sufficiently small and (e,d) G (£ x V)eo, then B = a{x,D)1/2 is well-
defined, its symbol b belongs to S™£2((£ x V)£o;Rn x R n ) , and b - a1'2 G
Operators with symbols holomorphic in a tube domain 377
b = b + k, (16.31)
and show that then Eq. (16.32) holds with j +1 instead of j ; this will finish
the proof of the Theorem.
By calculating b(x, D)2 with the help of Theorem 16.11 and taking
Eq. (16.30)-Eq. (16.32) and Eq. (16.17) into account, we find that
Let U C R n be an open set such that its closure contains the origin:
U 3 {0}, and consider the classes 5 m ( R n ) x ( R n +iU)) = Sm0(Rn x (R n +
iU)) introduced in Definition 14.1. If we work with the same spaces as in
Subsection 16.1, we can repeat all the statements there word by word by
simply adding +iU in the notation of each class of symbols. This means that
we make no use of the fact that symbols admit holomorphic continuation
into the tube domain, and if we consider European options with bounded
payoffs, e.g., European puts, we may do it and lose no essential information.
When the payoff is exponentially growing at infinity, as is the case with
European calls, appropriate spaces are spaces with exponential weights,
and one can consider the action of PDO in such spaces if and only if their
symbols admit analytic continuation into an appropriate tube domain. For
PDO of the class S^0{Rn x (R n + iU)), the part of 5 ( R " ) is played by
5 ( R n ; U) = {u | e ^ u e <S(Rn), V y € U},
e<x'v>a(x,D)u = (e<x'y)a(x,D)e-<XtV)")e<x'v)u
= a(x,D + iy)e<x'y)u.
where C can be chosen in the form of the RHS in Eq. (16.34). The sum
in the brackets defines the norm equivalent to the norm Eq. (16.33), which
finishes the proof of the Theorem. •
We live to the reader as an exercise the proof for spaces defined by Eq. (16.33)
with the weight (eu-x + e ^ 1 ) - 1 instead of e w - x + e"+ x .
It follows from the boundedness theorem and Eq. (16.36) that if A is suffi-
ciently large, the RHS in Eq. (16.35) is an invertible operator in I,2(R n ).
Hence, for these A, a(x, D) + A : Hs(Rn) —> L2(R n ) has a bounded left in-
verse (/ 4- t(A; x, D)) _1 6(A; x, D). Similarly, it has a bounded right inverse,
hence it is invertible. By Corollary 6.10, a(x, D) + A is invertible in <S(Rn).
Theorem 14.1 has thus been proved.
Proof of Theorem 14.4. a) We have a(x,£) := d2 - L(x,£) e 5j > 0 (R n x
( R n +iU)) for any U, and if U satisfies Eq. (14.14), we see that Eq.' (16.16)
holds for all (x,£) G R n x (R n + iU). Hence, the analogues of Theorem
16.9 and Theorem 16.6 for PDO with symbols holomorphic in a tube do-
main are valid, which gives not only Eq. (14.15)-Eq. (14.17) but the other
380 Elements of calculus of pseudodifferential operators
terms of the asymptotic expansion of the symbol as well (see Eq. (16.18)-
Eq. (16.20)).
b)-d) are special cases of the analogues of Theorems 6.11 and 6.7 for
PDO depending on a parameter.
Theorem 14.4 has thus been proved.
Proof of Lemma 14-1- By using Eq. (14.26) and repeating the proof of
Theorem 16.8, it is possible to show that the symbol of (\+r + ip(x, .D)) - 1 ,
call it R(\;x,£), satisfies estimates
decays as £ —> oo faster than any power of (£), together with all derivatives.
But the kernel kt(x, y) of a PDO exp[—t(r+ip(x, D)] is expressed via symbol
as follows:
+oo
/
-oo
and by differentiating under the integral sign, we conclude that kt € C°°(R n x
Rn).
Calculation of the omitted terms of the asymptotic expansion Eq. (14-28)
We use Theorem 16.8 and Theorem 16.5, and Eq. (16.9)-Eq. (16.13) with
\ + r + a(x, £) instead of a(x, £). For instance, the RHS of Eq. (14.28) with
two terms (the only case which may have some relevance to practice) is
By substituting Eq. (16.38) into Eq. (14.27) we can obtain the next terms
in Eq. (14.29). In particular, the first omitted term is equal to
transform the contour C$ by pushing it to the left. There is only one pole
of order three, and when it is crossed, by the residue theorem the term
(16.40)
appears. After that we can push the contour to infinity. Since the integrand
in Eq. (16.39) admits a bound via
CeexP(-r|£|)|£|-2.
exp(-r|£|)|£r2^ < d ,
-oo-H<r
K* : S(Rn) B u H u o / c e <S(Rn),
Then (see, e.g., Theorem 4.2 in Shubin (1978)) the following theorem is
valid.
Theorem 16.16 Ai is a PDO with the symbol <zi € 5 ^ 0 ( R n x R " ) ,
which admits an asymptotic expansion
Notice that Theorem 16.16 is valid for the other classes of symbols, which
we consider.
By using Theorem 16.16 and approximate calculations, similar to the
ones made in the preceding section, we can obtain approximate formulas
for the solution to the problem Eq. (9.25)-Eq. (9.26), and hence, approxi-
mations to prices of basket options and exchange options.
ii) then, assuming that a is a symbol from a good class (that is, its deriva-
tives are small in an appropriate sense) and using the explicit formulas
for the factors, obtained in Chapter 3, we show that the factors also be-
long to classes of symbols with sufficiently good properties, uniformly
in A € £e,c, for some 0 € {TT/2, IT) and C > 0;
iii) next, define the initial approximation to the resolvent by
R°X(X + A) = I-TX,
where T\ is an operator with the small norm, uniformly in A G Ee,c,
for some 8 e (7r/2, n) and C > 0;
iv) finally, we can calculate the resolvent with any precision, by summing
the following series:
(\ + A)-1 = (I + Tx+T$ + ---)R°x. (16.42)
By using Eq. (16.42) and Theorem 15.18, we can find an approximate price
for a barrier option under a Levy-like Feller process.
Of course, the concrete realization of the program requires long calcu-
lations and additional constructions.
This page is intentionally left blank
Bibliography
Abel, A.B. and Eberly, J.C. (1999) "The effects of irreversibility and uncertainty
on capital accumulation", Journal of Monetary Economics 44, 339-377.
Barles, G., Burdeau, J., Romano, M. and Samsoen, N. (1995) "Critical stock
price near expiration", Mathem. Finance 5, 77-95.
Barndorff-Nielsen, O.E. (1977) "Exponentially decreasing distributions for the
logarithm of particle size", Proc. Roy. Soc. London. Ser. A 3 5 3 , 401-419.
Barndorff-Nielsen, O.E. (1998) "Processes of normal inverse Gaussian type", Fi-
nance and Stochastics 2, 41-68.
Barndorff-Nielsen, O.E. and Jiang, W. (1998) "An initial analysis of some Ger-
man stock price series", Working Paper Series, 15, Aarhus: CAF Univ. of
Aarhus/Aarhus School of Business.
Barndorff-Nielsen, O.E. and Levendorskii, S.Z. (2001) "Feller processes of Normal
Inverse Gaussian type", Quantitative Finance 1, 318-331.
Barndorff-Nielsen, O.E., Nicolato, E. and Shephard, N. (2001) "Some recent de-
velopments in stochastic volatility modelling", Working paper
Barndorff-Nielsen, O.E. and Prause, K. (2001) "Apparent scaling", Finance and
Stochastics 5, 103-113.
Barndorff-Nielsen, O.E. and Shephard, N. (2001a) "Modelling by Levy processes
for financial econometrics", in Levy processes: Theory and Applications,
O.E. Barndorff-Nielsen, T. Mikosch and S. Resnik (Eds.), Birkhauser, 283-
318.
Barndorff-Nielsen, O.E. and Shephard, N. (2001b) "Non-Gaussian Ornstein-
Uhlenbeck-based models and some of their uses in financial economics (with
discussion)", J. Royal Stat. Soc. 6 3 . 167-241.
Barndorff-Nielsen, O.E. and Shephard, N. (2001c) "Normal modified stable pro-
cesses", Working paper
Barndorff-Nielsen O.E. and Shephard N. (2002) Levy based dynamics models for
financial economics, forthcoming.
Bernanke, B.S. (1983) "Irreversibility, uncertainty, and cyclical investment",
385
386 Bibliography
Eberlein, E. and Jacod, J. (1997) "On the range of options prices", Finance and
Stochastics 1, 131-140.
Eberlein, E., Kallsen, J. and Kirsten, J. (2001) "Risk management based on
stochastic volatility", FDM Preprint 72, University of Freiburg.
Eberlein, E. and Keller, U. (1995) "Hyperbolic distributions in Finance",
Bernoulli 1, 281-299.
Eberlein, E., Keller, U. and Prause, K. (1998)
"New insights into smile, mispricing and value at risk: The hyperbolic
model", Journ. of Business 7 1 , 371-406.
Eberlein, E. and Ozkan, F. (2001) "The defaultible Levy term structure: ratings
and restructuring", FDM Preprint 72, University of Freiburg.
Eberlein, E. and Prause, K. (2001) "The generalized hyperbolic model: finan-
cial derivatives and risk measures", in Mathematical Finance- Bachelier
Congress 2000, H. Geman, D. Madan, S. Pliska, T. Vorst (Eds.), Springer-
Verlag
Eberlein, E. and Raible, S. (1999) "Term structure models driven by general Levy
processes", Math. Finance 9, 31-53.
Eberlein, E. and Raible, S. (2001) "Some analytic facts on the generalized hyper-
bolic model", in Proceedings of the 3rd European Meetings of Mathemati-
cians, Birkhauser
Elliot, D.F. and Rao, K.R. (1982) Fast transforms: Algorithms, analyses, appli-
cations, Academic Press, New York
Eskin, G.I. (1973) Boundary problems for elliptic pseudo-differential equations,
Nauka, Moscow (Transl. of Mathematical Monographs, 52, Providence,
Rhode Island: Amer. Math. Soc, 1980).
Ethier, St. N. and Kurtz, Th. G. (1986) Markov processes, Wiley Series in Prob-
ability and Mathematical Statistics, John Wiley & Sons, New York
Fama, E.F. (1965) "The behavior of stock market prices", Journ. of Business 38,
34-105.
Fedoryuk, M.V. (1988) The saddle point method, Nauka, Moscow (In Russian)
Follmer, H. and Leukert, P. (1999) "Quantile hedging", Finance and Stochastics
3, 251-273.
Follmer, H. and Schweizer, M. (1991) "Hedging of contingent claims under in-
complete information", in Applied Stochastic Analysis, M.H.A. Davies and
R.J .Elliott (Eds.), New York, Gordon and Bleach, 389-414.
Follmer, H. and Sondermann, D. (1986) "Hedging of non-redundant contin-
gent claims", in Contributions to Mathematical Economics, A. Mas-Collel,
W. Hildebrand (Eds.), North-Holland, Amsterdam, 205-223.
Fons, J.S. (1994) "Using default rates to model the term structure of credit risk",
Financial Analysts Journal, September/October, 25-32.
Geman, H. and Yor, M. (1996) "Pricing and hedging double-barrier options: A
probabilistic approach", Mathem. Finance 6:4, 365-378.
Gerber, H. and Shiu, E. (1994) "Option pricing by Esscher transforms (with
discussion)", Trans. Soc. Actuaries 46, 99-191.
Bibliography 389
Geske, R. and Johnson, H.E. (1984) "The American put option valued analyti-
cally", Journ. of Financed 1511-1524.
Goll, T. and Kallsen, J. (2000) "Optimal portfolios for logarithmic utility",
Stochastic Processes and their Applications 89, 31-48.
Goll, Th. and Riischendorf, L. (2001) "Minimax and minimal distance martingale
measures and their relationship to portfolio optimization", Finance and
Stochastics 5:4, 557-581.
Gourieroux, C., Laurent, J.P. and Pham, P. "Mean-variance hedging and
numeraire", Mathem. Finance
Grubb, G. (1996) Functional calculus of pseudo-differential boundary problems,
Birkhauser, Boston Basel Berlin
Heston, S.L. (1993) "A closed-form solution for options with stochastic volatility,
with applications to bond and currency options", Rev. Finan. Stud. 6, 327-
343.
Hormander, L. (1985) The analysis of linear partial differential operators III,
Springer-Verlag, Berlin
Hubbard, R.G. (1994) "Investment under uncertainty: keeping one's options
open", Journal of Economic Literature 32, 1816-1831.
Hull, J.C. (2000) Options, futures & other derivatives. 4th ed., Prentice Hall, New
Jersey.
Hyakawa, K. and Kumano-go, H. (1971) "Complex powers of a system of pseudo-
differential operators", Proc. Japan Acad. 47, 359-364.
IngersoU, J.E., Jr. (2000) "Digital contracts: Simple tools for pricing complex
derivatives", Journal of Business 73:1, 67-88.
Jacob, N. (1996) Pseudo-differential operators and Markov processes, Academie
Verlag, Berlin
Jones, E.P., Mason, S.P. and Rosenfeld, E. (1984) "Contingent claims analysis
of corporate capital structures: An empirical investigation", Journal of
Finance 39, 611-627.
Kallsen, J. (1998) "Semimartingale modelling in Finance. Dissertation", Mathe-
matische Fakultat, Universitat Freiburg im Breisgau.
Kallsen, J. (2000) "Optimal portfolios for exponential Levy processes", Mathe-
matical Methods of Operations Research 51:3, 357-374.
Karatzas, I. and Shreve, S.E. (1998) Methods of mathematical Finance, Springer-
Verlag, Berlin Heidelberg New York.
Keller, U. (1997) "Realistic modelling of financial derivatives. Dissertation",
Mathematische Fakultat, Universitat Freiburg im Breisgau.
Koponen, I. (1995) "Analytic approach to the problem of convergence of truncated
Levy flights towards the Gaussian stochastic process", Physics Review E 52,
1197-1199.
Kramkov, D. and Schachermayer, W. (1999) "The asymptotic elasticity of util-
ity functions and optimal investment in incomplete markets", Ann. Appl.
Probab. 9, 904-950.
Kwok, Y.-K. (1998) Mathematical models of financial derivatives, Springer-
390 Bibliography
Sarig, O. and A. Warga, A. (1989) "Some empirical estimates of the risk structure
of interest rates", Journal of Finance 44, 1351-1360.
Sato, K. (1999) Levy processes and infinitely divisible distributions, Cambridge
University Press, Cambridge.
Schweizer, M. (1996) "Approximation pricing and the variance-optimal martin-
gale measure", Ann. Probab. 24:1, 206-236.
Schweizer, M. (1999) "A minimality property of the minimal martingale mea-
sure", Stat. Probab. Lett. 42, 27-31.
Seeley, R. (1967) "Complex powers of an elliptic operator", Amer.Math.Soc. Proc.
Symp. Pure Math. 10 288-307.
Sharpe, M.J. (1989) General theory of Markov processes, Academic Press, New
York.
Shiryaev, A.N. (1999) Essentials of stochastic Finance. Facts, models, theory,
World Scientific, Singapore New Jersey London Hong Kong.
Shubin, M.A. (1978) Spectral theory and pseudo-differential operators, Nauka,
Moskow (Transl. 1980)
Spitzer, F. (1964) Principles of random walk, Van Nostrand, Princeton.
Taylor, M. (1981) Pseudo-differential operators, Princeton University Press,
Princeton NJ
Triebel, H. (1978) Interpolation theory, function spaces, differential operators.
VEB Deutscher Verlag der Wissenschaften, Berlin.
Wiener, N. and Hopf, E. (1931) "Uber eine classe singularer Integralgleichungen",
Sitzber. Deutsch. Akad. Wiss. Berlin, Kl. Math. Phys. Tech., 696-706.
Wilmott P., Dewynne J.N. and Howison S. (1993) Option pricing: Mathematical
models and computation, Oxford Financial Press, Oxford.
Wilmott, P., Howison, S. and Dewynne, J. (1995) The mathematics of finan-
cial derivatives. A student introduction, Cambridge University Press, Cam-
bridge.
Yamazaki, M. (1986) "A quasi-homogeneous version of paradifferential operators,
I. Boundedness on spaces of Besov type", J.Fac.Sci.Univ. Tokyo Sect. I A
Math. 33, 131-174.
Yang, S.R. and Brorsen, B.W. (1992): "Nonlinear Dynamics of Daily Cash
Prices," American Journal of Agricultural Economics, 74:3, 706-715.
Yosida, K. (1964) Functional analysis, Springer-Verlag, Berlin, Heidelberg, New-
York
Zhou, C. (1997) "A jump-diffusion approach to modelling credit risk and valuing
defaultable securities", Board of Governors of the Federal Reserve System,
Finance and Economics Discussion Papers Series: 1997/27.
Zhu, J. (2000) Modular pricing of options (Lecture Notes in Economics and Math-
ematical Systems. Vol. 493), Springer-Verlag, Berlin Heidelberg New York.
Index
393
394 Index
variance, 3
volatility, 2, 4
implied, 5
historic, 5
smile, 5
surface, 5
Taking account of a diverse audience, the book has been written In such
a way that it is simple at the beginning and more technical in further
chapters, so that it is accessible to graduate students in relevant areas
and mathematicians without prior knowledge of finance or economics.
I I I M H HI 4*44 I
www.worldscientific.com
4955 he J!