Random bridges have gained significant attention in recent years due to their potential applications in various areas, particularly in information-based asset pricing models. This paper aims to explore the potential influence of the pinning point’s distribution on the memorylessness and stochastic dynamics of the bridge process. We introduce Lévy bridges with random length and random pinning points, and analyze their Markov property. Our study demonstrates that the Markov property of Lévy bridges depends on the nature of the distribution of their pinning points. The law of any random variables can be decomposed into singular continuous, discrete, and absolutely continuous parts with respect to the Lebesgue measure (Lebesgue’s decomposition theorem). We show that the Markov property holds when the pinning points’ law does not have an absolutely continuous part. Conversely, the Lévy bridge fails to exhibit Markovian behavior when the pinning point has an absolutely continuous part.

We consider the problem of obtaining effective representations for the solutions of linear, vector-valued stochastic differential equations (SDEs) driven by non-Gaussian pure-jump Lévy processes, and we show how such representations lead to efficient simulation methods. The processes considered constitute a broad class of models that find application across the physical and biological sciences, mathematics, finance, and engineering. Motivated by important relevant problems in statistical inference, we derive new, generalised shot-noise simulation methods whenever a normal variance-mean (NVM) mixture representation exists for the driving Lévy process, including the generalised hyperbolic, normal-gamma, and normal tempered stable cases. Simple, explicit conditions are identified for the convergence of the residual of a truncated shot-noise representation to a Brownian motion in the case of the pure Lévy process, and to a Brownian-driven SDE in the case of the Lévy-driven SDE. These results provide Gaussian approximations to the small jumps of the process under the NVM representation. The resulting representations are of particular importance in state inference and parameter estimation for Lévy-driven SDE models, since the resulting conditionally Gaussian structures can be readily incorporated into latent variable inference methods such as Markov chain Monte Carlo, expectation-maximisation, and sequential Monte Carlo.

Given a spectrally negative Lévy process, we predict, in an $L_1$ sense, the last passage time of the process below zero before an independent exponential time. This optimal prediction problem generalises [2], where the infinite-horizon problem is solved. Using a similar argument as that in [24], we show that this optimal prediction problem is equivalent to solving an optimal prediction problem in a finite-horizon setting. Surprisingly (unlike the infinite-horizon problem), an optimal stopping time is based on a curve that is killed at the moment the mean of the exponential time is reached. That is, an optimal stopping time is the first time the process crosses above a non-negative, continuous, and non-increasing curve depending on time. This curve and the value function are characterised as a solution of a system of nonlinear integral equations which can be understood as a generalisation of the free boundary equations (see e.g. [21, Chapter IV.14.1]) in the presence of jumps. As an example, we numerically calculate this curve in the Brownian motion case and for a compound Poisson process with exponential-sized jumps perturbed by a Brownian motion.

We propose non-asymptotic controls of the cumulative distribution function $\mathbb{P}(|X_{t}|\ge \varepsilon)$ , for any $t>0$ , $\varepsilon>0$ and any Lévy process X such that its Lévy density is bounded from above by the density of an $\alpha$ -stable-type Lévy process in a neighborhood of the origin.

We consider two-dimensional Lévy processes reflected to stay in the positive quadrant. Our focus is on the non-standard regime when the mean of the free process is negative but the reflection vectors point away from the origin, so that the reflected process escapes to infinity along one of the axes. Under rather general conditions, it is shown that such behaviour is certain and each component can dominate the other with positive probability for any given starting position. Additionally, we establish the corresponding invariance principle providing justification for the use of the reflected Brownian motion as an approximate model. Focusing on the probability that the first component dominates, we derive a kernel equation for the respective Laplace transform in the starting position. This is done for the compound Poisson model with negative exponential jumps and, by means of approximation, for the Brownian model. Both equations are solved via boundary value problem analysis, which also yields the domination probability when starting at the origin. Finally, certain asymptotic analysis and numerical results are presented.

Let $X_t^\sharp$ be a multivariate process of the form $X_t =Y_t - Z_t$ , $X_0=x$ , killed at some terminal time T, where $Y_t$ is a Markov process having only jumps of length smaller than $\delta$ , and $Z_t$ is a compound Poisson process with jumps of length bigger than $\delta$ , for some fixed $\delta>0$ . Under the assumptions that the summands in $Z_t$ are subexponential, we investigate the asymptotic behaviour of the potential function $u(x)= \mathbb{E}^x \int_0^\infty \ell\big(X_s^\sharp\big)ds$ . The case of heavy-tailed entries in $Z_t$ corresponds to the case of ‘big claims’ in insurance models and is of practical interest. The main approach is based on the fact that u(x) satisfies a certain renewal equation.

The signature of a path can be described as its full non-commutative exponential. Following T. Lyons, we regard its expectation, the expected signature, as a path space analogue of the classical moment generating function. The logarithm thereof, taken in the tensor algebra, defines the signature cumulant. We establish a universal functional relation in a general semimartingale context. Our work exhibits the importance of Magnus expansions in the algorithmic problem of computing expected signature cumulants and further offers a far-reaching generalization of recent results on characteristic exponents dubbed diamond and cumulant expansions with motivations ranging from financial mathematics to statistical physics. From an affine semimartingale perspective, the functional relation may be interpreted as a type of generalized Riccati equation.

In this paper, we construct operator fractional Lévy motion (ofLm), a broad class of infinitely divisible stochastic processes that are covariance operator self-similar and have wide-sense stationary increments. The ofLm class generalizes the univariate fractional Lévy motion as well as the multivariate operator fractional Brownian motion (ofBm). OfLm can be divided into two types, namely, moving average (maofLm) and real harmonizable (rhofLm), both of which share the covariance structure of ofBm under assumptions. We show that maofLm and rhofLm admit stochastic integral representations in the time and Fourier domains, and establish their distinct small- and large-scale limiting behavior. We also characterize time-reversibility for ofLm through parametric conditions related to its Lévy measure. In particular, we show that, under non-Gaussianity, the parametric conditions for time-reversibility are generally more restrictive than those for the Gaussian case (ofBm).

In this paper, we study the pricing of vulnerable Asian options with liquidity risk. We employ general Lévy processes to capture the changes in the liquidity discount factors and the information processes of all assets. In the proposed pricing model, we obtain the closed-form pricing formula of vulnerable Asian options using the Fourier transform methods. Finally, the derived pricing formula is used to illustrate the effects of asymmetric jump risk, and the effects are relatively stable on (vulnerable) Asian options with different moneynesses.

In this paper we consider the one-dimensional, biased, randomly trapped random walk with infinite-variance trapping times. We prove sufficient conditions for the suitably scaled walk to converge to a transformation of a stable Lévy process. As our main motivation, we apply subsequential versions of our results to biased walks on subcritical Galton–Watson trees conditioned to survive. This confirms the correct order of the fluctuations of the walk around its speed for values of the bias that yield a non-Gaussian regime.

This paper formulates and solves a new problem of global practical inverse optimal exponential path-tracking control of mobile robots driven by Lévy processes with unknown characteristics. The control design is based on a new inverse optimal control design for nonlinear systems driven by Lévy processes and ensures global practical exponential stability almost surely and in the pth moment for the path-tracking errors. Moreover, it minimizes cost function that penalizes tracking errors and control torques without having to solve a Hamilton–Jacobi–Bellman or Hamilton–Jaccobi–Isaacs equation.

Motivated by a recent paper (Budd (2018)), where a new family of positive self-similar Markov processes associated to stable processes appears, we introduce a new family of Lévy processes, called the double hypergeometric class, whose Wiener–Hopf factorisation is explicit, and as a result many functionals can be determined in closed form.

We consider an insurance entity endowed with an initial capital and a surplus process modelled as a Brownian motion with drift. It is assumed that the company seeks to maximise the cumulated value of expected discounted dividends, which are declared or paid in a foreign currency. The currency fluctuation is modelled as a Lévy process. We consider both cases: restricted and unrestricted dividend payments. It turns out that the value function and the optimal strategy can be calculated explicitly.

We consider the approximate Euler scheme for Lévy-drivenstochastic differential equations.We study the rate of convergence in law of the paths.We show that when approximating the small jumps by Gaussianvariables, the convergence is much faster than when simplyneglecting them.For example, when the Lévy measure of the driving processbehaves like |z|−1−αdz near 0, for some α ∈ (1,2), we obtain an error of order 1/√n with a computational cost of order nα. For a similar error when neglecting the small jumps, see[S. Rubenthaler, Numerical simulation of the solution of a stochastic differential equation driven by a Lévy process. Stochastic Process. Appl. 103 (2003) 311–349], the computational cost is of ordernα/(2−α), which is huge when α is close to 2. In the same spirit, we study the problem of the approximation of a Lévy-driven S.D.E.by a Brownian S.D.E. when the Lévy process has no large jumps. Our results rely on some results of [E. Rio, Upper bounds for minimal distances in the central limit theorem. Ann. Inst. Henri Poincaré Probab. Stat. 45 (2009) 802–817] about the central limit theorem, in the spirit of the famous paper by Komlós-Major-Tsunády [J. Komlós, P. Major and G. Tusnády, An approximation of partial sums of independentrvs and the sample df I. Z. Wahrsch. verw. Gebiete 32 (1975) 111–131].

Let (Xt, t ≥ 0) be a Lévy process started at 0, with Lévymeasure ν. We consider the first passage time T x of(Xt, t ≥ 0) to level x > 0, and Kx := XTx - x theovershoot and Lx := x- XTx- the undershoot. We first provethat the Laplace transform of the random triple (Tx,Kx,Lx )satisfies some kind of integral equation. Second, assuming thatν admits exponential moments, we show that $(\widetilde{T_x},K_x,L_x)$ converges in distribution asx → ∞, where $\widetilde{T_x}$ denotes a suitablerenormalization of T x .


In this paper we study asymptotic behavior of convexrearrangements of Lévy processes. In particular weobtain Glivenko-Cantelli-type strong limit theorems for the convexifications when the corresponding Lévy measureis regularly varying at + with exponent α ∈ (1,2).

We continue the investigation started in a previous paper, onweak convergence to infinitely divisible distributions with finitevariance. In the present paper, we study this problem for someweakly dependent random variables, including in particularassociated sequences. We obtain minimal conditions expressed interms of individual random variables. As in the i.i.d. case, wedescribe the convergence to the Gaussian and the purelynon-Gaussian parts of the infinitely divisible limit. We alsodiscuss the rate of Poisson convergence and emphasize the specialcase of Bernoulli random variables. The proofs aremainly based on Lindeberg's method.

Arbitrage-free prices u of European contracts on risky assets whoselog-returns are modelled by Lévy processes satisfya parabolic partial integro-differential equation (PIDE) $\partial_t u + {\mathcal{A}}[u] = 0$ .This PIDE is localized tobounded domains and the error due to this localization isestimated. The localized PIDE is discretized by theθ-scheme in time and a wavelet Galerkin method with N degrees of freedom in log-price space. The dense matrix for ${\mathcal{A}}$ can be replaced by a sparse matrix in the wavelet basis, and the linear systemsin each implicit time step are solved approximativelywith GMRES in linear complexity. The total work of the algorithm for M time steps is bounded byO(MN(log(N))2) operations and O(Nlog(N)) memory.The deterministic algorithm gives optimal convergence rates (up to logarithmic terms) for the computed solutionin the same complexity as finite difference approximationsof the standard Black–Scholes equation.Computational examples for various Lévy price processes are presented.

In this paper, the class of periodic Ornstein-Uhlenbeck processes is defined. It is shown that periodic Ornstein-Uhlenbeck processes are stationary Markov random fields and the class of stationary distributions is characterized. In particular, any self-decomposable distribution is the stationary distribution of some periodic Ornstein-Uhlenbeck process. As examples, gamma periodic Ornstein-Uhlenbeck processes and Gaussian periodic Ornstein-Uhlenbeck processes are considered.

A class of continuous-time models is developed for modelling data with heavy tails and long-range dependence. These models are based on the Green function solutions of fractional differential equations driven by Lévy noise. Some exact results on the second- and higher-order characteristics of the equations are obtained. Applications to stochastic volatility of asset prices and macroeconomics are provided.