Bilateral gamma distributions and processes in financial mathematics

We present a class of Lévy processes for modelling financial market fluctuations: bilateral Gamma processes. Our starting point is to explore the properties of bilateral Gamma distributions, and then we turn to their associated Lévy processes. We treat exponential Lévy stock models with an underlying bilateral Gamma process as well as term structure models driven by bilateral Gamma processes, and apply our results to a set of real financial data (DAX 1996–1998).     

Tridiagonal implicit method to evaluate European and American options under infinite activity Lévy models

We propose an efficient implicit method to evaluate European and American options when the underlying asset follows an infinite activity Lévy model. Since the Lévy measure of the infinite activity model has the singularity at the origin, we approximate infinitely many small jumps by samples of a diffusion. The proposed methods to solve partial integro–differential equations for European options and linear complementarity problems for American options via an operator splitting method involve solving linear systems with tridiagonal matrices and so can significantly reduce the computations associated with the discrete integral operators. The numerical experiments verify that the proposed method has the second-order convergence rate under an infinite activity Lévy model.     

Characterization of dependence of multidimensional Lévy processes using Lévy copulas

This paper suggests Lévy copulas in order to characterize the dependence among components of multidimensional Lévy processes. This concept parallels the notion of a copula on the level of Lévy measures. As for random vectors, a version of Sklar's theorem states that the law of a general multivariate Lévy process is obtained by combining arbitrary univariate Lévy processes with an arbitrary Lévy copula. We construct parametric families of Lévy copulas and prove a limit theorem, which indicates how to obtain the Lévy copula of a multivariate Lévy process X from the ordinary copula of the random vector X_t for small t.     

BSDEs driven by time-changed Lévy noises and optimal control

We study backward stochastic differential equations (BSDEs) for time-changed Lévy noises when the time-change is independent of the Lévy process. We prove existence and uniqueness of the solution and we obtain an explicit formula for linear BSDEs and a comparison principle. BSDEs naturally appear in control problems. Here we prove a sufficient maximum principle for a general optimal control problem of a system driven by a time-changed Lévy noise. As an illustration we solve the mean–variance portfolio selection problem.     

Jump-adapted discretization schemes for Lévy-driven SDEs

We present new algorithms for weak approximation of stochastic differential equations driven by pure jump Lévy processes. The method uses adaptive non-uniform discretization based on the times of large jumps of the driving process. To approximate the solution between these times we replace the small jumps with a Brownian motion. Our technique avoids the simulation of the increments of the Lévy process, and in many cases achieves better convergence rates than the traditional Euler scheme with equal time steps. To illustrate the method, we discuss an application to option pricing in the Libor market model with jumps.     

The set-indexed Lévy process: Stationarity,Markov and sample paths properties

We present a satisfactory definition of the important class of Lévy processes indexed by a general collection of sets. We use a new definition for increment stationarity of set-indexed processes to obtain different characterizations of this class. As an example, the set-indexed compound Poisson process is introduced. The set-indexed Lévy process is characterized by infinitely divisible laws and a Lévy–Khintchine representation. Moreover, the following concepts are discussed: projections on flows, Markov properties, and pointwise continuity. Finally the study of sample paths leads to a Lévy–Itô decomposition. As a corollary, the semi-martingale property is proved.     

Continuous dependence estimates for viscosity solutions of integro-PDEs

We present a general framework for deriving continuous dependence estimates for, possibly polynomially growing, viscosity solutions of fully nonlinear degenerate parabolic integro-PDEs. We use this framework to provide explicit estimates for the continuous dependence on the coefficients and the “Lévy measure” in the Bellman/Isaacs integro-PDEs arising in stochastic control/differential games. Moreover, these explicit estimates are used to prove regularity results and rates of convergence for some singular perturbation problems. Finally, we illustrate our results on some integro-PDEs arising when attempting to price European/American options in an incomplete stock market driven by a geometric Lévy process. Many of the results obtained herein are new even in the convex case where stochastic control theory provides an alternative to our pure PDE methods.     

A stochastic linear–quadratic problem with Lévy processes and its application to finance

We study a Linear–Quadratic Regulation (LQR) problem with Lévy processes and establish the closeness property of the solution of the multi-dimensional Backward Stochastic Riccati Differential Equation (BSRDE) with Lévy processes. In particular, we consider multi-dimensional and one-dimensional BSRDEs with Teugel’s martingales which are more general processes driven by Lévy processes. We show the existence and uniqueness of solutions to the one-dimensional regular and singular BSRDEs with Lévy processes by means of the closeness property of the BSRDE and obtain the optimal control for the non-homogeneous case. An application of the backward stochastic differential equation approach to a financial (portfolio selection) problem with full and partial observation cases is provided.     

A stochastic heat equation with the distributions of Lévy processes as its invariant measures

We consider a linear heat equation on a half line with an additive noise chosen properly in such a manner that its invariant measures are a class of distributions of Lévy processes. Our assumption on the corresponding Lévy measure is, in general, mild except that we need its integrability to show that the distributions of Lévy processes are the only invariant measures of the stochastic heat equation.     

Analysis of generalized Lévy white noise functionals

Employing the Segal-Bargmann transform (S-transform for abbreviation) of regular Lévy white noise functionals, we define and study the generalized Lévy white noise functionals by means of their functional representations acting on test functionals. The main results generalize (Gaussian) white noise analysis initiated by T. Hida to non-Gaussian cases. Thanks to the closed form of the S-transform of Lévy white noise functionals obtained in our previous paper, we are able to define and study the renormalization of products of Lévy white noises, multiplication operator by Lévy white noises, and the differential operators with respect to a Lévy white noise and their adjoint operators. In the courses of our investigation we also obtain a formula for the products of multiple Lévy-Itô stochastic integrals. As applications, we discuss the existence of Hitsuda-Skorokhod integral for Lévy processes, Kubo-Takenaka formula for Lévy processes, and Itô formula for generalized Lévy white noise functionals.     

Adaptive nonparametric estimation for Lévy processes observed at low frequency

This article deals with adaptive nonparametric estimation for Lévy processes observed at low frequency. For general linear functionals of the Lévy measure, we construct kernel estimators, provide upper risk bounds and derive rates of convergence under regularity assumptions.     

White Noise of Poisson Random Measures

We develop a white noise theory for Poisson random measures associated with a pure jump Lévy process. The starting point of this theory is the chaos expansion of Itô. We use this to construct the white noise of a Poisson random measure, which takes values in a certain distribution space. Then we show, how a Skorohod/Itô integral for point processes can be represented by a Bochner integral in terms of white noise of the random measure and a Wick product. Further, based on these concepts we derive a generalized Clark–Haussmann–Ocone theorem with respect to a combination of Gaussian noise and pure jump Lévy noise. We apply this theorem to obtain an explicit formula for partial observation minimal variance portfolios in financial markets, driven by Lévy processes. As an example we compute the closest hedge to a binary option.     

Lévy-Driven Carma Processes

Properties and examples of continuous-time ARMA (CARMA) processes driven by Lévy processes are examined. By allowing Lévy processes to replace Brownian motion in the definition of a Gaussian CARMA process, we obtain a much richer class of possibly heavy-tailed continuous-time stationary processes with many potential applications in finance, where such heavy tails are frequently observed in practice. If the Lévy process has finite second moments, the correlation structure of the CARMA process is the same as that of a corresponding Gaussian CARMA process. In this paper we make use of the properties of general Lévy processes to investigate CARMA processes driven by Lévy processes {W(t)} without the restriction to finite second moments. We assume only that W (1) has finite r-th absolute moment for some strictly positive r. The processes so obtained include CARMA processes with marginal symmetric stable distributions.     

Analytical valuation of catastrophe equity options with negative exponential jumps

A catastrophe put option is valuable in the event that the underlying asset price is below the strike price; in addition, a specified catastrophic event must have happened and influenced the insured company. This paper analyzes the valuation of catastrophe put options under deterministic and stochastic interest rates when the underlying asset price is modeled through a Lévy process with finite activity. We provide explicit analytical formulas for evaluating values of catastrophe put options. The numerical examples illustrate how financial risks and catastrophic risks affect the prices of catastrophe put options.     

A Decomposition Theorem for Lévy Processes on Local Fields

The study of Lévy processes on local fields has been initiated by Albeverio et al. (1985)–(1998) and Evans (1989)–(1998). In this paper, a decomposition theorem for Lévy processes on local fields is given in terms of a structure result for measures on local fields and a Lévy–Khinchine representation. It is shown that a measure on a local field can be decomposed into three parts: a spherically symmetric measure, a totally non-spherically symmetric measure and a singular measure. We show that if the Radon–Nikodym derivative of the absolutely continuous part of a Lévy measure on a local field is locally constant, the Lévy process is the sum of a spherically symmetric random walk, a finite or countable set of totally non, spherically symmetric Lévy processes with single balls as support of their Lévy measure, end a singular Lévy process. These processes are independent. Explicit formulae for the transition function are obtained.     

Equilibruim approach of asset pricing under Lévy process

This work considers the equilibrium approach of asset pricing for Lévy process. It derives the equity premium and pricing kernel analytically for the stock price process, obtains an equilibrium option pricing formula, and explains some empirical evidence such as the negative variance risk premium, implied volatility smirk, and negative skewness risk premium by comparing the physical and risk-neutral distributions of the log return. Different from most of the current studies in equilibrium pricing under jump diffusion models, this work models the underlying asset price as the exponential of a Lévy process and thus allows nearly an arbitrage distribution of the jump component.     

Optimal control and dependence modeling of insurance portfolios with Lévy dynamics

In this paper we are interested in optimizing proportional reinsurance and investment policies in a multidimensional Lévy-driven insurance model. The criterion is that of maximizing exponential utility. Solving the classical Hamilton-Jacobi-Bellman equation yields that the optimal retention level keeps a constant amount of claims regardless of time and the company’s wealth level.A special feature of our construction is to allow for dependencies of the risk reserves in different business lines. Dependence is modeled via an Archimedean Lévy copula. We derive a sufficient and necessary condition for an Archimedean Lévy generator to create a multidimensional positive Lévy copula in arbitrary dimension.Based on these results we identify structure conditions for the generator and the Lévy measure of an Archimedean Lévy copula under which an insurance company reinsures a larger fraction of claims from one business line than from another.     

Multilevel Monte Carlo simulation for Lévy processes based on the Wiener–Hopf factorisation

In Kuznetsov et al. (2011) a new Monte Carlo simulation technique was introduced for a large family of Lévy processes that is based on the Wiener–Hopf decomposition. We pursue this idea further by combining their technique with the recently introduced multilevel Monte Carlo methodology. Moreover, we provide here for the first time a theoretical analysis of the new Monte Carlo simulation technique in Kuznetsov et al. (2011) and of its multilevel variant for computing expectations of functions depending on the historical trajectory of a Lévy process. We derive rates of convergence for both methods and show that they are uniform with respect to the “jump activity” (e.g. characterised by the Blumenthal–Getoor index). We also present a modified version of the algorithm in Kuznetsov et al. (2011) which combined with the multilevel methodology obtains the optimal rate of convergence for general Lévy processes and Lipschitz functionals. This final result is only a theoretical one at present, since it requires independent sampling from a triple of distributions which is currently only possible for a limited number of processes.     

Jump-diffusion processes in random environments

In this paper we investigate jump-diffusion processes in random environments which are given as the weak solutions of SDEs. We formulate conditions ensuring existence and uniqueness in law of solutions. We investigate the Markov property. To prove uniqueness we solve a general martingale problem for càdlàg processes. This result is of independent interest. Application of our results to generalized exponential Lévy model are present in the last section.     

Stochastic integration for Lévy processes with values in Banach spaces

A stochastic integral of Banach space valued deterministic functions with respect to Banach space valued Lévy processes is defined. There are no conditions on the Banach spaces or on the Lévy processes. The integral is defined analogously to the Pettis integral. The integrability of a function is characterized by means of a radonifying property of an integral operator associated with the integrand. The integral is used to prove a Lévy–Itô decomposition for Banach space valued Lévy processes and to study existence and uniqueness of solutions of stochastic Cauchy problems driven by Lévy processes.