Inspired by the works of Landriault et al. (2011), (2014), we study the Gerber–Shiu distribution at Parisian ruin with exponential implementation delays for a spectrally negative Lévy insurance risk process. To be more specific, we study the so-called Gerber–Shiu distribution for a ruin model where at each time the surplus process goes negative, an independent exponential clock is started. If the clock rings before the surplus becomes positive again then the insurance company is ruined. Our methodology uses excursion theory for spectrally negative Lévy processes and relies on the theory of so-called scale functions. In particular, we extend the recent results of Landriault et al. (2011), (2014).

In this paper we focus on a class of linear Hawkes processes with general immigrants. These are counting processes with shot-noise intensity, including self-excited and externally excited patterns. For such processes, we introduce the concept of the age pyramid which evolves according to immigration and births. The virtue of this approach that combines an intensity process definition and a branching representation is that the population age pyramid keeps track of all past events. This is used to compute new distribution properties for a class of Hawkes processes with general immigrants which generalize the popular exponential fertility function. The pathwise construction of the Hawkes process and its underlying population is also given.

We consider the distribution of the polygonal paths joining partial sums of classical Kloosterman sums $\text{Kl}_{p}(a)$, as $a$ varies over $\mathbf{F}_{p}^{\times }$ and as $p$ tends to infinity. Using independence of Kloosterman sheaves, we prove convergence in the sense of finite distributions to a specific random Fourier series. We also consider Birch sums, for which we can establish convergence in law in the space of continuous functions. We then derive some applications.

We establish a convexity property for the hitting probabilities of discrete random walks in ${\mathbb Z}^d$ (discrete harmonic measures). For d = 2 this implies a recent result on the convexity of the density of certain harmonic measures.

In this paper we consider an Ornstein–Uhlenbeck (OU) process (M(t))t≥0 whose parameters are determined by an external Markov process (X(t))t≥0 on a finite state space {1, . . ., d}; this process is usually referred to as Markov-modulated Ornstein–Uhlenbeck. We use stochastic integration theory to determine explicit expressions for the mean and variance of M(t). Then we establish a system of partial differential equations (PDEs) for the Laplace transform of M(t) and the state X(t) of the background process, jointly for time epochs t = t1, . . ., tK. Then we use this PDE to set up a recursion that yields all moments of M(t) and its stationary counterpart; we also find an expression for the covariance between M(t) and M(t + u). We then establish a functional central limit theorem for M(t) for the situation that certain parameters of the underlying OU processes are scaled, in combination with the modulating Markov process being accelerated; interestingly, specific scalings lead to drastically different limiting processes. We conclude the paper by considering the situation of a single Markov process modulating multiple OU processes.

We derive a new equation for the optimal investment boundary of a general irreversible investment problem under exponential Lévy uncertainty. The problem is set as an infinite time-horizon, two-dimensional degenerate singular stochastic control problem. In line with the results recently obtained in a diffusive setting, we show that the optimal boundary is intimately linked to the unique optional solution of an appropriate Bank–El Karoui representation problem. Such a relation and the Wiener–Hopf factorization allow us to derive an integral equation for the optimal investment boundary. In case the underlying Lévy process hits any point in R with positive probability we show that the integral equation for the investment boundary is uniquely satisfied by the unique solution of another equation which is easier to handle. As a remarkable by-product we prove the continuity of the optimal investment boundary. The paper is concluded with explicit results for profit functions of Cobb–Douglas type and CES type. In the former case the function is separable and in the latter case nonseparable.

The full-information duration problem with a random number N of objects is considered. These objects appear sequentially and their values Xk are observed, where Xk, independent of N, are independent and identically distributed random variables from a known continuous distribution. The objective of the problem is to find a stopping rule that maximizes the duration of holding a relative maximum (e.g. the kth object is a relative maximum if Xk = max{X1, X2, . . ., Xk}). We assume that N is a random variable with a known upper bound n, so two models, Model 1 and Model 2, can be considered according to whether the planning horizon is N or n. The structure of the optimal rule, which depends on the prior distribution assumed on N, is examined. The monotone rule is defined and a necessary and sufficient condition for the optimal rule to be monotone is given for both models. Special attention is paid to the class of priors such that N / n converges, as n → ∞, to a random variable Vm having density fVm(v) = m(1 - v)m-1, 0 ≤ v ≤ 1 for a positive integer m. An interesting feature is that the optimal duration (relative to n) for Model 2 is just (m + 1) times as large as that for Model 1 asymptotically.

We study the distribution Ex[exp(-q∫0t1(a,b)(Xs)ds); Xt ∈ dy], where -∞ ≤ a < b < ∞, and where q, t > 0 and x ∈ R for a spectrally negative Lévy process X. More precisely, we identify the Laplace transform with respect to t of this measure in terms of the scale functions of the underlying process. Our results are then used to price step options and the particular case of an exponential spectrally negative Lévy jump-diffusion model is discussed.

We provide exact computations for the drift of random walks in dependent random environments, including k-dependent and moving average environments. We show how the drift can be characterized and evaluated using Perron–Frobenius theory. Comparing random walks in various dependent environments, we demonstrate that their drifts can exhibit interesting behavior that depends significantly on the dependency structure of the random environment.

We propose a novel theoretical characterization of the optimal 'Gittins index' policy in multi-armed bandit problems with non-Gaussian, infinitely divisible reward distributions. We first construct a continuous-time, conditional Lévy process which probabilistically interpolates the sequence of discrete-time rewards. When the rewards are Gaussian, this approach enables an easy connection to the convenient time-change properties of a Brownian motion. Although no such device is available in general for the non-Gaussian case, we use optimal stopping theory to characterize the value of the optimal policy as the solution to a free-boundary partial integro-differential equation (PIDE). We provide the free-boundary PIDE in explicit form under the specific settings of exponential and Poisson rewards. We also prove continuity and monotonicity properties of the Gittins index in these two problems, and discuss how the PIDE can be solved numerically to find the optimal index value of a given belief state.

The problem of pricing arithmetic Asian options is nontrivial, and has attracted much interest over the last two decades. This paper provides a method for calculating bounds on option prices and approximations to option deltas in a market where the underlying asset follows a geometric Lévy process. The core idea is to find a highly correlated, yet more tractable proxy to the event that the option finishes in-the-money. The paper provides a means for calculating the joint characteristic function of the underlying asset and proxy processes, and relies on Fourier methods to compute prices and deltas. Numerical studies show that the lower bound provides accurate approximations to prices and deltas, while the upper bound provides good though less accurate results.

We discuss simulation of sensitivities or Greeks of multi-asset European style options under a special Lévy process model: that is, the subordinated Brownian motion model. The Malliavin calculus method combined with Monte Carlo and quasi-Monte Carlo methods is used in the simulations. Greeks are expressed in terms of the expectations of the option payoff functions multiplied by the weights involving Malliavin derivatives for multi-asset options. Numerical results show that the Malliavin calculus method is usually more efficient than the finite difference method for options with nonsmooth payoffs. The superiority of the former method over the latter is even more significant when both are combined with quasi-Monte Carlo methods.

We introduce a new model of competition on growing networks. This extends the preferential attachment model, with the key property that node choices evolve simultaneously with the network. When a new node joins the network, it chooses neighbours by preferential attachment, and selects its type based on the number of initial neighbours of each type. The model is analysed in detail, and in particular, we determine the possible proportions of the various types in the limit of large networks. An important qualitative feature we find is that, in contrast to many current theoretical models, often several competitors will coexist. This matches empirical observations in many real-world networks.

We present some correlated fractional counting processes on a finite-time interval. This will be done by considering a slight generalization of the processes in Borges et al. (2012). The main case concerns a class of space-time fractional Poisson processes and, when the correlation parameter is equal to 0, the univariate distributions coincide with those of the space-time fractional Poisson process in Orsingher and Polito (2012). On the one hand, when we consider the time fractional Poisson process, the multivariate finite-dimensional distributions are different from those presented for the renewal process in Politi et al. (2011). We also consider a case concerning a class of fractional negative binomial processes.

Superposition is a mapping on point configurations that sends the n-tuple into the n-point configuration , counted with multiplicity. It is an additive set operation such that the superposition of a k-point configuration in is a kn-point configuration in . A Poisson superposition process is the superposition in of a Poisson process in the space of finite-length -valued sequences. From properties of Poisson processes as well as some algebraic properties of formal power series, we obtain an explicit expression for the Janossy measure of Poisson superposition processes, and we study their law under domain restriction. Examples of well-known Poisson superposition processes include compound Poisson, negative binomial, and permanental (boson) processes.

In an optimal variance stopping problem the goal is to determine the stopping time at which the variance of a sequentially observed stochastic process is maximized. A solution method for such a problem has been recently provided by Pedersen (2011). Using the methodology developed by Pedersen and Peskir (2012), our aim is to show that the solution to the initial problem can be equivalently obtained by constraining the variance stopping problem to the expected size of the stopped process and then by maximizing the solution to the latter problem over all the admissible constraints. An application to a diffusion process used for modeling the dynamics of interest rates illustrates the proposed technique.

The Ginibre point process (GPP) is one of the main examples of determinantal point processes on the complex plane. It is a recurring distribution of random matrix theory as well as a useful model in applied mathematics. In this paper we briefly overview the usual methods for the simulation of the GPP. Then we introduce a modified version of the GPP which constitutes a determinantal point process more suited for certain applications, and we detail its simulation. This modified GPP has the property of having a fixed number of points and having its support on a compact subset of the plane. See Decreusefond et al. (2013) for an extended version of this paper.

Shot-noise processes are used in applied probability to model a variety of physical systems in, for example, teletraffic theory, insurance and risk theory, and in the engineering sciences. In this paper we prove a large deviation principle for the sample-paths of a general class of multidimensional state-dependent Poisson shot-noise processes. The result covers previously known large deviation results for one-dimensional state-independent shot-noise processes with light tails. We use the weak convergence approach to large deviations, which reduces the proof to establishing the appropriate convergence of certain controlled versions of the original processes together with relevant results on existence and uniqueness.

This paper concerns the asymptotic behavior of a random variable Wλ resulting from the summation of the functionals of a Gibbsian spatial point process over windows Qλ ↑ ℝd. We establish conditions ensuring that Wλ has volume order fluctuations, i.e. they coincide with the fluctuations of functionals of Poisson spatial point processes. We combine this result with Stein's method to deduce rates of a normal approximation for Wλ as λ → ∞. Our general results establish variance asymptotics and central limit theorems for statistics of random geometric and related Euclidean graphs on Gibbsian input. We also establish a similar limit theory for claim sizes of insurance models with Gibbsian input, the number of maximal points of a Gibbsian sample, and the size of spatial birth-growth models with Gibbsian input.

We give an account of the characteristics that result from reflecting a drifting local martingale (i.e. the sum of a local martingale and a multiple of its quadratic variation process) in 0 and b > 0. We present conditions which guarantee the existence of finite moments of what is required to keep the reflected process within its boundaries. Also, we derive an associated law of large numbers and a central limit theorem which apply when the input is continuous. Similar results for integrals of the paths of the reflected process are also presented. These results are in close agreement to what has previously been shown for Brownian motion.