Consumption-investment problem with transaction costs for Lévy-driven price processes

We consider an optimal control problem for a linear stochastic integro-differential equation with conic constraints on the phase variable and with the control of singular–regular type. Our setting includes consumption-investment problems for models of financial markets in the presence of proportional transaction costs, where the prices of the assets are given by a geometric Lévy process, and the investor is allowed to take short positions. We prove that the Bellman function of the problem is a viscosity solution of an HJB equation. A uniqueness theorem for the solution of the latter is established. Special attention is paid to the dynamic programming principle.     

A Feynman–Kac-type formula for Lévy processes with discontinuous killing rates

The challenge to fruitfully merge state-of-the-art techniques from mathematical finance and numerical analysis has inspired researchers to develop fast deterministic option pricing methods. As a result, highly efficient algorithms to compute option prices in Lévy models by solving partial integro-differential equations have been developed. In order to provide a solid mathematical foundation for these methods, we derive a Feynman–Kac representation of variational solutions to partial integro-differential equations that characterize conditional expectations of functionals of killed time-inhomogeneous Lévy processes. We allow a wide range of underlying stochastic processes, comprising processes with Brownian part as well as a broad class of pure jump processes such as generalized hyperbolic, multivariate normal inverse Gaussian, tempered stable, and \(\alpha\)-semistable Lévy processes. By virtue of our mild regularity assumptions as to the killing rate and the initial condition of the partial integro-differential equation, our results provide a rigorous basis for numerous applications in financial mathematics and in probability theory. We implement a Galerkin scheme to solve the corresponding pricing equation numerically and illustrate the effect of a killing rate.     

On Kolmogorov equations for anisotropic multivariate Lévy processes

For d-dimensional exponential Lévy models, variational formulations of the Kolmogorov equations arising in asset pricing are derived. Well-posedness of these equations is verified. Particular attention is paid to pure jump, d-variate Lévy processes built from parametric, copula dependence models in their jump structure. The domains of the associated Dirichlet forms are shown to be certain anisotropic Sobolev spaces. Singularity-free representations of the Dirichlet forms are given which remain bounded for piecewise polynomial, continuous functions of finite element type. We prove that the variational problem can be localized to a bounded domain with explicit localization error bounds. Furthermore, we collect several analytical tools for further numerical analysis.     

Optimal dividend strategy with transaction costs for an upward jump model

In this paper, we consider the optimal dividend problem with transaction costs when the incomes of a company can be described by an upward jump model. Both fixed and proportional costs are considered in the problem. The value function is defined as the expected total discounted dividends up to the time of ruin. Although the same problem has already been studied in the pure diffusion model and the spectrally negative Lévy process, the optimal dividend problem in an upward jump model has two different aspects in determining the optimal dividends barrier and in the property of the value function. First, the value function is twice continuous differentiable in the diffusion case, but it is not in the jump model. Second, under the spectrally negative Lévy process, downward jumps will not cause any payment actions; however, it might trigger dividend payments when there are upward jumps. In deriving the optimal barriers, we show that the value function is bounded by a linear function. Using this property, we establish the verification theorem for the value function. By solving the quasi-variational inequalities associated with this problem, we obtain the closed-form solution to the value function and hence the optimal dividend strategy when the income sizes follow a common exponential distribution. In the presence of a fixed transaction cost, it is shown that the optimal strategy is a two-barrier policy, and the optimal barriers are only dependent on the fixed cost and not the proportional cost. A numerical example is used to illustrate how the fixed cost plays a significant role in the optimal dividend strategy and also the value function. Moreover, an increased fixed cost results in larger but less frequent dividend payments.     

Wavelet Galerkin pricing of American options on Lévy driven assets

The price of an American-style contract on assets driven by a class of Markov processes containing, in particular, Lévy processes of pure jump type with infinite jump activity is expressed as the solution of a parabolic variational integro-differential inequality (PIDI). A Galerkin discretization in logarithmic price using a wavelet basis is presented. Log-linear complexity in each time-step is achieved by wavelet compression of the moment matrix of the price process’ jump measure and by wavelet preconditioning of the large matrix LCPs at each time-step. Efficiency is demonstrated by numerical experiments for pricing American put contracts on various jump-diffusion and pure jump models. Failure of the smooth pasting principle is observed for American put contracts for certain finite variation pure jump price processes.     

Financial inverse problem and reconstruction of?infinitely divisible distributions with Gaussian component

We consider an inverse problem of partial integro-differential equations of market prices of call options with many maturities and strike prices for geometric Lévy processes. We show the well-posedness (reconstruction, uniqueness, and stability) of the inverse problem among the class of infinitely divisible distributions with analyticity.     

Application of Homotopy Analysis Method to Option Pricing Under Lévy Processes

Option pricing under the Lévy process has been considered an important research direction in the field of financial engineering, where a closed-form expression for the standard European option is available due to the existence of analytically tractable characteristic function according to the Lévy–Khinchin representation. However, this approach cannot be applied to exotic derivatives (such as barrier options) directly, although a large volume of exotic derivatives are actively traded in the current options market. An alternative approach is to solve the corresponding partial integro-differential equation (PIDE) numerically, which is, in fact, time-consuming and is not computationally tractable in general. In this paper, we apply the so-called homotopy analysis method (HAM) to solve the corresponding PIDE in a semi analytic form, being obtained from the following three steps: (1) Apply the Fourier transform to convert the PIDE to an ordinal differential equitation (ODE), and construct a differential system of ODEs. (2) Solve the system of ODEs, where each differential equation is shown to have an analytical solution. (3) Express the option price using the sum of infinite series, where each term may be expressed analytically and derived by applying Steps (1) and (2) recursively. To illustrate our technique more precisely, we take the variance gamma model as an example and provide the semi-analytic form. Numerical examples demonstrate a fast convergence of our proposed method to the prices of European and down-and-out call options with a few number of terms. Note that this method is easy to implement and can be applied to other types of options under general Lévy processes.     

Forward rate models with linear volatilities

The existence of solutions to the Heath?CJarrow?CMorton equation of the bond market with linear volatility and general Lévy random factor is studied. Conditions for the existence and non-existence of solutions in the class of bounded fields are presented. For the existence of solutions, the Lévy process should necessarily be without a Gaussian part and without negative jumps. If this is the case, then necessary and sufficient conditions for the existence are formulated either in terms of the behavior of the Lévy measure of the noise near the origin or the behavior of the Laplace exponent of the noise at infinity.