Because of a string of coincidences involving my choice of courses for Part III and various lecturers’ choices about course content, I didn’t learn what a Levy process until a few weeks’ ago. Trying to get my head around the Levy-Khintchine formula took a little while, so the following is what I would have liked to have been able to find back then.
A Levy process is an adapted stochastic process started from 0 at time zero, and with stationary, independent increments. This is reminiscent, indeed a generalisation, of the definition of Brownian motion. In that case, we were able to give a concrete description of the distribution of . For a general Levy process, we have
So the distribution of is infinitely divisible, that is, can be expressed as the distribution of the sum n iid random variables for all n. Viewing this definition in terms of convolutions of distributions may be more helpful, especially as we will subsequently consider characteristic functions. If this is the first time you have seen this property, note that it is not a universal property. For example, it is not clear how to write a U[0,1] random variable as a convolution of two iid RVs. Note that exactly the same argument suffices to show that the distribution of
is infinitely divisible.
It will be most convenient to work with the characteristic functions
By stationarity of increments, we can show that this is equal to
This function is called the characteristic exponent. The argument resembles that used for Cauchy’s functional equations, by dealing first with the rationals using stationarity of increments, then lifting to the reals by the (right-)continuity of
As ever, uniquely determines the distribution of
, and so it also uniquely determines the distribution of Levy process. The only condition on
is that it be the characteristic function of an infinitely divisible distribution. This condition is given explicitly by the Levy-Khintchine formula.
Levy-Khintchine
is the characteristic function of an infinitely divisible distribution iff
for , Q a quadratic form on
, and
a so-called Levy measure satisfying
.
This looks a bit arbitrary, so first let’s explain what each of these terms ‘means’.
comes from a drift of
. Note that a deterministic linear function is a (not especially interesting) Levy process.
comes from a Brownian part
.
The rest corresponds to the jump part of the process. Note that a Poisson process is an example of a Levy process, hence why we might consider thinking about jumps in the first place. The reason why there is an indicator function floating around is that we have to think about two regimes separately, namely large and small jumps. Jumps of size bounded below cannot happen too often as otherwise the process might explode off to infinity in finite time with positive probability. On the other hand, infinitesimally small jumps can happen very often (say on a dense set) so long as everything is controlled to prevent an explosion on the macroscopic scale.
There is no canonical choice for where the divide between these regimes happens, but conventionally this is taken to be at . The restriction on the Levy measure near 0 ensures that the sum of the squares all jumps up some finite time converges absolutely.
gives the intensity of a standard compound Poisson process. The jumps are well-spaced, and so it is a relatively simple calculation to see that the characteristic function is
The intensity gives infinitely many hits in finite time, so if the expectation of this measure is not 0, we explode immediately. We compensate by drifting away from this at rate
To make this more rigorous, we should really consider then take a limit, but this at least explains where all the terms come from. Linearity allows us to interchange integrals and inner products, to get the term
If the process has bounded variation, then we must have Q=0, and also
that is, not too many jumps on an |x| scale. In this case, then this drift component is well-defined and linear , so can be incorporated with the drift term at the beginning of the Levy-Khintchine expression. If not, then there are some
for which it does not exist.
There are some other things to be said about Levy processes, including
- Stable Levy processes, where
, which induces the rescaling-invariance property:
. The distribution of each
is then also a stable distribution.
- Resolvents, where instead of working with the process itself, we work with the distribution of the process at a random exponential time.
Related articles
- Gaussian tail bounds and a word of caution about CLT (eventuallyalmosteverywhere.wordpress.com)
Pingback: Subordinators and the Arcsine rule | Eventually Almost Everywhere
Pingback: Random Mappings for Cycle Deletion | Eventually Almost Everywhere
So, first thanks for this note. Second, I would say that’s i’m not a specialist of stochastics calculus but i’m interested to this mathematical branch. So, I have just a small question and it’s as follows: can we define a similer calcul when psi(lambda)===psi(t,lambda).
Pingback: 100k Views | Eventually Almost Everywhere
Levy-Khintchine
\Psi(\lambda) is the characteristic function -> \Psi(\lambda) is the characteristic component.
Hi Dominic!
I have seen papers where the characteristic function (phi) for certain Lévy processes X_t have been defined as E(exp^{lambda*X_t}) = exp{phi(lambda)}. For what kind of Lévy processes does this apply?
Many thanks – Issy
so these levy processes are a 3-manifold. the semi-martingales form a proper orthogonal subgroup so you can shown that the quotient group is a vector space and the remainder is an ideal. note the zeroes of the riemann zeta function are the annihilator in this ring.
this is an excellent site. that gamma term we need to transform from the second derivative over the density to the arrival intensity integrated around the boundary and use the cauchy residue theorem. the arrival times are memoryless you need to show its positive definite symetric self adjoint with 0.5*sigma^2*t on the diag. note the real part of this is a semi-martingale and satisfies the normal flow equations with the -0.5*sigma^2*t the drift in the stochastic exponential to make it a martingale. by convolution of the gamma characteristic over the rationals you find the complex covariance matrix defines an operator the spectrum of which is the so-called zeroes of the riemann zeta function. the critical strip are those correlation terms which are in the spectrum but are not eigenvalues.