High order schemes for the tempered fractional diffusion equations

Lévy flight models whose jumps have infinite moments are mathematically used to describe the superdiffusion in complex systems. Exponentially tempering Lévy measure of Lévy flights leads to the tempered stable Lévy processes which combine both the α-stable and Gaussian trends; and the very large jumps are unlikely and all their moments exist. The probability density functions of the tempered stable Lévy processes solve the tempered fractional diffusion equation. This paper focuses on designing the high order difference schemes for the tempered fractional diffusion equation on bounded domain. The high order difference approximations, called the tempered and weighted and shifted Grünwald difference (tempered-WSGD) operators, in space are obtained by using the properties of the tempered fractional calculus and weighting and shifting their first order Grünwald type difference approximations. And the Crank-Nicolson discretization is used in the time direction. The stability and convergence of the presented numerical schemes are established; and the numerical experiments are performed to confirm the theoretical results and testify the effectiveness of the schemes.     

Multifractal analysis of foreign exchange data

In this paper we perform multifractal analyses of five daily Foreign Exchange (FX) rates. These techniques are currently used in turbulence to characterize scaling and intermittency. We show the multifractal nature of FX returns, and estimate the three parameters in the universal multifactal framework, which characterize all small and medium intensity fluctuations, at all scales. For large fluctuations, we address the question of hyperbolic (fat) tails of the distributions which are characterized by a fourth parameter, the tail index. We studied both the prices fluctuations and the returns, finding no systematic difference in the scaling exponents in the two cases. We discuss and compare our results with several recent studies, and show how the additive models are not compatible with data: Brownian, fractional Brownian, Lévy, Truncated Lévy and fractional Lévy models. We analyse in this framework the ARCH(1), GARCH(1,1) and HARCH (7) models, and show that their structure functions scaling exponents are undistinguishable from that of Brownian motion, which means that these models do not adequately describe the scaling properties of the statistics of the data. Our results indicate that there might exist a multiplicative ‘flux of financial information’, which conditions small‐scale statistics to large‐scale values, as an analogy with the energy flux in turbulence. Copyright © 1999 John Wiley & Sons, Ltd.     

Transition probabilities of Lévy‐type processes: Parametrix construction

We present an existence result for Lévy‐type processes which requires only weak regularity assumptions on the symbol with respect to the space variable x. Applications range from existence and uniqueness results for Lévy‐driven SDEs with Hölder continuous coefficients to existence results for stable‐like processes and Lévy‐type processes with symbols of variable order. Moreover, we obtain heat kernel estimates for a class of Lévy and Lévy‐type processes. The paper includes an extensive list of Lévy(‐type) processes satisfying the assumptions of our results.     

A comprehensive mathematical approach to exotic option pricing

This work illustrates how several new pricing expressions for exotic options can be derived within a Lévy framework by employing a unique pricing expression. To the purpose, a unifying formula is obtained by solving some nested Cauchy problem for pseudodifferential equations generalizing Black–Scholes PDE. The main result extends (Agliardi R. The quintessential option pricing formula under Lévy processes. Applied Mathematics Letters 2009; 22:1626‐1631) and is a powerful tool for generating new valuation expressions. Several examples of pricing formulas under the Lévy processes are provided to illustrate the flexibility of the method. Some of them are new in the financial literature. Finally, many existing pricing formulas of the traditional Gaussian model are easily obtained as a by‐product. Copyright © 2012 John Wiley & Sons, Ltd.     

Polynomial spectral collocation method for space fractional advection–diffusion equation

This article discusses the spectral collocation method for numerically solving nonlocal problems: one‐dimensional space fractional advection–diffusion equation; and two‐dimensional linear/nonlinear space fractional advection–diffusion equation. The differentiation matrixes of the left and right Riemann–Liouville and Caputo fractional derivatives are derived for any collocation points within any given bounded interval. Several numerical examples with different boundary conditions are computed to verify the efficiency of the numerical schemes and confirm the exponential convergence; the physical simulations for Lévy–Feller advection–diffusion equation and space fractional Fokker–Planck equation with initial δ‐peak and reflecting boundary conditions are performed; and the eigenvalue distributions of the iterative matrix for a variety of systems are displayed to illustrate the stabilities of the numerical schemes in more general cases. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 514–535, 2014     

Least squares estimator of Ornstein-Uhlenbeck processes driven by fractional Lévy processes with periodic mean

We deal with the least squares estimator for the drift parameters of an Ornstein-Uhlenbeck process with periodic mean function driven by fractional Lévy process. For this estimator, we obtain consistency and the asymptotic distribution. Compared with fractional Ornstein-Uhlenbeck and Ornstein-Uhlenbeck driven by Lévy process, they can be regarded both as a Lévy generalization of fractional Brownian motion and a fractional generaliza- tion of Lévy process.     

Estimates of densities for Lévy processes with lower intensity of large jumps

We obtain general lower estimates of transition densities of jump Lévy processes. We use them for processes with Lévy measures having bounded support, processes with exponentially decaying Lévy measures for large times and for processes with high intensity of small jumps for small times.     

Finite difference methods and their physical constraints for the fractional klein‐kramers equation

Incorporating subdiffusive mechanisms into the Klein‐Kramers formalism leads to the fractional Klein‐Kramers equation. Then, the equation can effectively describe subdiffusion in the presence of an external force field in the phase space. This article presents the finite difference methods for numerically solving the fractional Klein‐Kramers equation and does the detailed stability and error analyses. The stability condition, mvβ ≤ 16, shows the ratio between the kinetic energy of the particle and the temperature of the fluid can not be too large, which well agrees with the physical property of the subdiffusive particle, we call it “physical constraint.” The numerical examples are provided to verify the theoretical results on rate of convergence. Moreover, we simulate the fractional Klein‐Kramers dynamics and the simulation results further confirm the effectiveness of our numerical schemes. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1561–1583, 2010     

Spike‐Type Solutions to One Dimensional Gierer–Meinhardt Model with Lévy Flights

The Gierer–Meinhardt model with Lévy flights is shown to give rise to patterns of spikes with algebraically decaying tails. The spike shape is given by a solution to a fractional differential equation. Near an equilibrium formation the spikes drift according to the differential equations of the form known for Fickian diffusion, but with a new homoclinic. A nonlocal eigenvalue problem of a new type is formulated and studied. The system is less stable due to the Lévy flights, though the behavior of eigenvalues is changed mainly quantitatively.     

A risk model driven by Lévy processes

We present a general risk model where the aggregate claims, as well as the premium function, evolve by jumps. This is achieved by incorporating a Lévy process into the model. This seeks to account for the discrete nature of claims and asset prices. We give several explicit examples of Lévy processes that can be used to drive a risk model. This allows us to incorporate aggregate claims and premium fluctuations in the same process. We discuss important features of such processes and their relevance to risk modeling. We also extend classical results on ruin probabilities to this model. Copyright © 2003 John Wiley & Sons, Ltd.     

Closed‐form approximations for spread options in Lévy markets

We provide new closed‐form approximations for the pricing of spread options in three specific instances of exponential Lévy markets, ie, when log‐returns are modeled as Brownian motions (Black‐Scholes model), variance gamma processes (VG model), or normal inverse Gaussian processes (NIG model). For the specific case of exchange options (spread options with zero strike), we generalize the well‐known Margrabe formula (1978) that is valid in a Black‐Scholes model to the VG model under a homogeneity assumption.     

Spatial regularity of semigroups generated by Lévy type operators

We apply the probabilistic coupling approach to establish spatial regularity of semigroups associated with Lévy type operators, by assuming that the corresponding martingale problem is well posed. In particular, we can prove the Lipschitz continuity of the associated semigroups, when the coefficients are Hölder continuous but the corresponding Lévy kernel may be singular.     

Gradient estimates of Dirichlet heat kernels for unimodal Lévy processes

Under some mild assumptions on the Lévy measure and the symbol we obtain gradient estimates of Dirichlet heat kernels for pure‐jump isotropic unimodal Lévy processes in .     

Multiscaling comparative analysis of time series and geophysical phenomena

Different methods are used to determine the scaling exponents associated with a time series describing a complex dynamical process, such as those observed in geophysical systems. Many of these methods are based on the numerical evaluation of the variance of a diffusion process whose step increments are generated by the data. An alternative method focuses on the direct evaluation of the scaling coefficient of the Shannon entropy of the same diffusion distribution. The combined use of these methods can efficiently distinguish between fractal Gaussian and Lévy‐walk time series and help to discern between alternative underling complex dynamics. © 2005 Wiley Periodicals, Inc. Complexity 10: 51–56, 2005     

Lévy processes on a first order model

The classical notion of Lévy process is generalized to one that takes as its values probabilities on a first or‐der model equipped with a commutative semigroup. This is achieved by applying a convolution product on definable probabilities and the infinite divisibility with respect to it (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Uniform asymptotic estimates for ruin probabilities of renewal risk models with exponential Lévy process investment returns and dependent claims

This paper investigates the ruin probabilities of a renewal risk model with stochastic investment returns and dependent claim sizes. The investment is described as a portfolio of one risk‐free asset and one risky asset whose price process is an exponential Lévy process. The claim sizes are assumed to follow a one‐sided linear process with independent and identically distributed step sizes. When the step‐size distribution is heavy tailed, we establish some uniform asymptotic estimates for the ruin probabilities of this renewal risk model. Copyright © 2012 John Wiley & Sons, Ltd.     

Preconditioned iterative methods for fractional diffusion models in finance

Most recent qualitative models for financial assets assume that the dynamics of underlying equity prices follows a jump or Lévy process. It has been evident that some most intricate characteristics of such dynamics can be captured by CGMY and KoBoL procedures. The prices of financial derivatives with such models satisfy fractional partial differential equations or partial integro‐differential equations. This study focuses at aforementioned fractional equations and discretizes them via a monotone Crank–Nicolson procedure. A spatial extrapolation strategy is introduced to ensure an overall second‐order accuracy in approximations. Preconditioned conjugate gradient normal residual methods are incorporated for solving resulted linear systems. Numerical examples are given to illustrate the accuracy and efficiency of the novel computational approaches implemented. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1382–1395, 2015     

Fractional Lévy Processes as a Result of Compact Interval Integral Transformation

Fractional Brownian motion can be represented as an integral of a deterministic kernel w.r.t. an ordinary Brownian motion either on infinite or compact interval. In previous literature fractional Lévy processes are defined by integrating the infinite interval kernel w.r.t. a general Lévy process. In this article we define fractional Lévy processes using the com pact interval representation.

We prove that the fractional Lévy processes presented via different integral transformations have the same finite dimensional distributions if and only if they are fractional Brownian motions. Also, we present relations between different fractional Lévy processes and analyze the properties of such processes. A financial example is introduced as well.     

A comparison of numerical solutions of fractional diffusion models in finance

We compare the numerical solutions of three fractional partial differential equations that occur in finance. These fractional partial differential equations fall in the class of Lévy models. They are known as the FMLS (Finite Moment Log Stable), CGMY and KoBol models. Conditions for the convergence of each of these models is obtained.     

Approximation properties for solutions to non‐Lipschitz stochastic differential equations with Lévy noise

In this paper, we consider the non‐Lipschitz stochastic differential equations and stochastic functional differential equations with delays driven by Lévy noise, and the approximation theorems for the solutions to these two kinds of equations will be proposed respectively. Non‐Lipschitz condition is much weaker condition than the Lipschitz one. The simplified equations will be defined to make its solutions converge to that of the corresponding original equations both in the sense of mean square and probability, which constitute the approximation theorems. Copyright © 2014 John Wiley & Sons, Ltd.