This paper considers the first passage times to constant boundaries and the two-sided exit problem for Lévy process with a characteristic exponent in which at least one of the two jumps having rational Laplace transforms. The joint distribution of the first passage times and undershoot/overshoot are obtained. The processes recover many models that have appeared in the literature such as the compound Poisson risk models, the perturbed compound Poisson risk models, and their dual ones. As applications, we obtain the solutions for popular path-dependent options such as lookback and barrier options in terms of Laplace transforms.

We continue with the systematic study of the speed of extinction of continuous-state branching processes in Lévy environments under more general branching mechanisms. Here, we deal with the weakly subcritical regime under the assumption that the branching mechanism is regularly varying. We extend recent results of Li and Xu (2018) and Palau et al. (2016), where it is assumed that the branching mechanism is stable, and complement the recent articles of Bansaye et al. (2021) and Cardona-Tobón and Pardo (2021), where the critical and the strongly and intermediate subcritical cases were treated, respectively. Our methodology combines a path analysis of the branching process together with its Lévy environment, fluctuation theory for Lévy processes, and the asymptotic behaviour of exponential functionals of Lévy processes. Our approach is inspired by the last two previously cited papers, and by Afanasyev et al. (2012), where the analogue was obtained.

We develop a novel Monte Carlo algorithm for the vector consisting of the supremum, the time at which the supremum is attained, and the position at a given (constant) time of an exponentially tempered Lévy process. The algorithm, based on the increments of the process without tempering, converges geometrically fast (as a function of the computational cost) for discontinuous and locally Lipschitz functions of the vector. We prove that the corresponding multilevel Monte Carlo estimator has optimal computational complexity (i.e. of order $\varepsilon^{-2}$ if the mean squared error is at most $\varepsilon^2$) and provide its central limit theorem (CLT). Using the CLT we construct confidence intervals for barrier option prices and various risk measures based on drawdown under the tempered stable (CGMY) model calibrated/estimated on real-world data. We provide non-asymptotic and asymptotic comparisons of our algorithm with existing approximations, leading to rule-of-thumb principles guiding users to the best method for a given set of parameters. We illustrate the performance of the algorithm with numerical examples.

In this paper, we analyze a two-queue random time-limited Markov-modulated polling model. In the first part of the paper, we investigate the fluid version: fluid arrives at the two queues as two independent flows with deterministic rate. There is a single server that serves both queues at constant speeds. The server spends an exponentially distributed amount of time in each queue. After the completion of such a visit time to one queue, the server instantly switches to the other queue, i.e., there is no switch-over time.

For this model, we first derive the Laplace–Stieltjes transform (LST) of the stationary marginal fluid content/workload at each queue. Subsequently, we derive a functional equation for the LST of the two-dimensional workload distribution that leads to a Riemann–Hilbert boundary value problem (BVP). After taking a heavy-traffic limit, and restricting ourselves to the symmetric case, the BVP simplifies and can be solved explicitly.

In the second part of the paper, allowing for more general (Lévy) input processes and server switching policies, we investigate the transient process limit of the joint workload in heavy traffic. Again solving a BVP, we determine the stationary distribution of the limiting process. We show that, in the symmetric case, this distribution coincides with our earlier solution of the BVP, implying that in this case the two limits (stationarity and heavy traffic) commute.

We consider a Lévy process Y(t) that is not continuously observed, but rather inspected at Poisson( $\omega$ ) moments only, over an exponentially distributed time $T_\beta$ with parameter $\beta$ . The focus lies on the analysis of the distribution of the running maximum at such inspection moments up to $T_\beta$ , denoted by $Y_{\beta,\omega}$ . Our main result is a decomposition: we derive a remarkable distributional equality that contains $Y_{\beta,\omega}$ as well as the running maximum process $\bar Y(t)$ at the exponentially distributed times $T_\beta$ and $T_{\beta+\omega}$ . Concretely, $\overline{Y}(T_\beta)$ can be written as the sum of two independent random variables that are distributed as $Y_{\beta,\omega}$ and $\overline{Y}(T_{\beta+\omega})$ . The distribution of $Y_{\beta,\omega}$ can be identified more explicitly in the two special cases of a spectrally positive and a spectrally negative Lévy process. As an illustrative example of the potential of our results, we show how to determine the asymptotic behavior of the bankruptcy probability in the Cramér–Lundberg insurance risk model.

Growth-fragmentation processes describe the evolution of systems in which cells grow slowly and fragment suddenly. Despite originating as a way to describe biological phenomena, they have recently been found to describe the lengths of certain curves in statistical physics models. In this note, we describe a new growth-fragmentation process connected to random planar maps with faces of large degree, having as a key ingredient the ricocheted stable process recently discovered by Budd. The process has applications to the excursions of planar Brownian motion and Liouville quantum gravity.

We consider fragmentation processes with values in the space of marked partitions of $\mathbb{N}$, i.e. partitions where each block is decorated with a nonnegative real number. Assuming that the marks on distinct blocks evolve as independent positive self-similar Markov processes and determine the speed at which their blocks fragment, we get a natural generalization of the self-similar fragmentations of Bertoin (Ann. Inst. H. Poincaré Prob. Statist.38, 2002). Our main result is the characterization of these generalized fragmentation processes: a Lévy–Khinchin representation is obtained, using techniques from positive self-similar Markov processes and from classical fragmentation processes. We then give sufficient conditions for their absorption in finite time to a frozen state, and for the genealogical tree of the process to have finite total length.

Consider the strong subordination of a multivariate Lévy process with a multivariate subordinator. If the subordinate is a stack of independent Lévy processes and the components of the subordinator are indistinguishable within each stack, then strong subordination produces a Lévy process; otherwise it may not. Weak subordination was introduced to extend strong subordination, always producing a Lévy process even when strong subordination does not. Here we prove that strong and weak subordination are equal in law under the aforementioned condition. In addition, we prove that if strong subordination is a Lévy process then it is necessarily equal in law to weak subordination in two cases: firstly when the subordinator is deterministic, and secondly when it is pure-jump with finite activity.

We define a new family of multivariate stochastic processes over a finite time horizon that we call generalised Liouville processes (GLPs). GLPs are Markov processes constructed by splitting Lévy random bridges into non-overlapping subprocesses via time changes. We show that the terminal values and the increments of GLPs have generalised multivariate Liouville distributions, justifying their name. We provide various other properties of GLPs and some examples.

We investigate the long-time behavior of the Ornstein–Uhlenbeck process driven by Lévy noise with regime switching. We provide explicit criteria on the transience and recurrence of this process. Contrasted with the Ornstein–Uhlenbeck process driven simply by Brownian motion, whose stationary distribution must be light-tailed, both the jumps caused by the Lévy noise and the regime switching described by a Markov chain can derive the heavy-tailed property of the stationary distribution. The different role played by the Lévy measure and the regime-switching process is clearly characterized.

In a recent article, Döring et al. (2018) conditioned a stable process with two-sided jumps to avoid an interval. As usual, the strategy was to find an invariant function for the process killed on entering the interval and to show that the corresponding h-transformed process is indeed the process conditioned to avoid an interval in a meaningful way. In the present article we consider the case of a completely asymmetric stable process. It turns out that the invariant function found by Döring et al. does not exist or is not invariant, but nonetheless we will characterise the conditioned process as a Markov process.

We characterize various forms of positive dependence, such as association, positive supermodular association and dependence, and positive orthant dependence, for jump-Feller processes. Such jump processes can be studied through their state-space dependent Lévy measures. It is through these Lévy measures that we will provide our characterization. Finally, we present applications of these results to stochastically monotone Feller processes, including Lévy processes, the Ornstein–Uhlenbeck process, pseudo-Poisson processes, and subordinated Feller processes.

The fractional nonhomogeneous Poisson process was introduced by a time change of the nonhomogeneous Poisson process with the inverse α-stable subordinator. We propose a similar definition for the (nonhomogeneous) fractional compound Poisson process. We give both finite-dimensional and functional limit theorems for the fractional nonhomogeneous Poisson process and the fractional compound Poisson process. The results are derived by using martingale methods, regular variation properties and Anscombe’s theorem. Eventually, some of the limit results are verified in a Monte Carlo simulation.

In the literature on optimal stopping, the problem of maximizing the expected discounted reward over all stopping times has been explicitly solved for some special reward functions (including (x+)ν, (ex − K)+, (K − e− x)+, x ∈ ℝ, ν ∈ (0, ∞), and K > 0) under general random walks in discrete time and Lévy processes in continuous time (subject to mild integrability conditions). All such reward functions are continuous, increasing, and logconcave while the corresponding optimal stopping times are of threshold type (i.e. the solutions are one-sided). In this paper we show that all optimal stopping problems with increasing, logconcave, and right-continuous reward functions admit one-sided solutions for general random walks and Lévy processes, thereby generalizing the aforementioned results. We also investigate in detail the principle of smooth fit for Lévy processes when the reward function is increasing and logconcave.

In this paper we present a set of results relating to the occupation time α(t) of a process X(·). The first set of results concerns exact characterizations of α(t), e.g. in terms of its transform up to an exponentially distributed epoch. In addition, we establish a central limit theorem (entailing that a centered and normalized version of α(t)∕t converges to a zero-mean normal random variable as t→∞) and the tail asymptotics of ℙ(α(t)∕t≥q). We apply our findings to spectrally positive Lévy processes reflected at the infimum and establish various new occupation time results for the corresponding model.

We introduce and study the class of branching-stable point measures, which can be seen as an analog of stable random variables when the branching mechanism for point measures replaces the usual addition. In contrast with the classical theory of stable (Lévy) processes, there exists a rich family of branching-stable point measures with a negative scaling exponent, which can be described as certain Crump‒Mode‒Jagers branching processes. We investigate the asymptotic behavior of their cumulative distribution functions, that is, the number of atoms in (-∞, x] as x→∞, and further depict the genealogical lineage of typical atoms. For both results, we rely crucially on the work of Biggins (1977), (1992).

In this paper we consider the factor analysis for Lévy-driven multivariate price models with stochastic volatility. Our main aim is to provide conditions on the volatility process under which we can possibly reduce the dimension of the driving Lévy motion. We find that these conditions depend on a particular form of the multivariate Lévy process. In some settings we concentrate on nondegenerate symmetric α-stable Lévy motions.

Motivated by a common mathematical finance topic, we discuss the reciprocal of the exit time from a cone of planar Brownian motion which also corresponds to the exponential functional of Brownian motion in the framework of planar Brownian motion. We prove a conjecture of Vakeroudis and Yor (2012) concerning infinite divisibility properties of this random variable and present a novel simple proof of the result of DeBlassie (1987), (1988) concerning the asymptotic behavior of the distribution of the Bessel clock appearing in the skew-product representation of planar Brownian motion, as t→∞. We use the results of the windings approach in order to obtain results for quantities associated to the pricing of Asian options.

Let X be a jump-diffusion process and X* its running supremum. In this paper we first show that for any t > 0, the law of the pair (X*t, Xt) has a density with respect to the Lebesgue measure. This allows us to show that for any t > 0, the law of the pair formed by the random variable Xt and the running supremum X*t of X at time t can be characterized as a weak solution of a partial differential equation concerning the distribution of the pair (X*t, Xt). Then we obtain an expression of the marginal density of X*t for all t > 0.