It is well known that between all processes with independent increments, essentially only the Brownian motion and the Poisson process possess the chaotic representation property (CRP). Thus, a natural question appears: What is an appropriate analog of the CRP in the case of a general Lévy process. At least three approaches are possible here. The first one, due to Itô, uses the CRP of the Brownian motion and the Poisson process, as well as the representation of a Lévy process through those processes. The second approach, due to Nualart and Schoutens, consists of representing any square-integrable random variable as a sum of multiple stochastic integrals constructed with respect to a family of orthogonalized centered power jumps processes. The third approach, never applied before to the Lévy processes, uses the idea of orthogonalization of polynomials with respect to a probability measure defined on the dual of a nuclear space. The main aims of this paper are to develop the three approaches in the case of a general (ℝ-valued) Lévy process on a Riemannian manifold and (what is more important) to understand a relationship between these approaches. We apply the obtained results to the gamma, Pascal, and Meixner processes, in which case the analysis related to the orthogonalized polynomials becomes essentially simpler and richer than in the general case.
Implementable coupling of Lévy process and Brownian motion
