For spectrally negative Lévy processes, we prove several fluctuation results involving a general draw-down time, which is a downward exit time from a dynamic level that depends on the running maximum of the process. In particular, we find expressions of the Laplace transforms for the two-sided exit problems involving the draw-down time. We also find the Laplace transforms for the hitting time and creeping time over the running-maximum related draw-down level, respectively, and obtain an expression for a draw-down associated potential measure. The results are expressed in terms of scale functions for the spectrally negative Lévy processes.

The eigentime identity for one-dimensional diffusion processes on the halfline with an entrance boundary at ∞ is obtained by using the trace of the deviation kernel. For the case of an exit boundary at ∞, a similar eigentime identity is presented with the aid of the Green function. Explicit equivalent statements are also listed in terms of the strong ergodicity or the uniform decay for diffusion processes.

A recent result of Ding, Lee and Peres (2012) expressed the cover time of the random walk on a graph in terms of generic chaining for the commute distance. Their argument is based on Dynkin's isomorphism theorem. The purpose of this article is to present an alternative approach to this problem, based only on elementary hitting time estimates and chaining arguments.

In this note we compute the Laplace transform of hitting times, to fixed levels, of integrated geometric Brownian motion. The transform is expressed in terms of the gamma and confluent hypergeometric functions. Using a simple Itô transformation and standard results on hitting times of diffusion processes, the transform is characterized as the solution to a linear second-order ordinary differential equation which, modulo a change of variables, is equivalent to Kummer's equation.

We study the two-sided exit problem of a time-homogeneous diffusion process with tax payments of loss-carry-forward type and obtain explicit formulae for the Laplace transforms associated with the two-sided exit problem. The expected present value of tax payments until default, the two-sided exit probabilities, and, hence, the nondefault probability with the default threshold equal to the lower bound are solved as immediate corollaries. A sufficient and necessary condition for the tax identity in ruin theory is discovered.

The inverse stable subordinator provides a probability model for time-fractional differential equations, and leads to explicit solution formulae. This paper reviews properties of the inverse stable subordinator, and applications to a variety of problems in mathematics and physics. Several different governing equations for the inverse stable subordinator have been proposed in the literature. This paper also shows how these equations can be reconciled.

Two classical stochastic processes are considered, the Ehrenfest process, introduced in 1907 in the kinetic theory of gases to describe the heat exchange between two bodies, and the Engset process, one of the early (1918) stochastic models of communication networks. In this paper we investigate the asymptotic behavior of the distributions of hitting times of these two processes when the number of particles/sources goes to infinity. Results concerning the hitting times of boundaries in particular are obtained. We rely on martingale methods; a key ingredient is an important family of simple nonnegative martingales, an analogue, for the Ehrenfest process, of the exponential martingales used in the study of random walks or of Brownian motion.

This paper provides a simple proof for the fact that the hitting time to an infrequently visited subset for a one-dependent regenerative process converges weakly to an exponential distribution. Special cases are positive recurrent Harris chains and Harris processes. The paper further extends this class of limit theorems to ‘rewards’ that are cumulated to the hitting time of such a rare set.

We consider the level hitting times τy = inf{t ≥ 0 | Xt = y} and the running maximum process Mt = sup{Xs | 0 ≤ s ≤ t} of a growth-collapse process (Xt)t≥0, defined as a [0, ∞)-valued Markov process that grows linearly between random ‘collapse’ times at which downward jumps with state-dependent distributions occur. We show how the moments and the Laplace transform of τy can be determined in terms of the extended generator of Xt and give a power series expansion of the reciprocal of Ee−sτy. We prove asymptotic results for τy and Mt: for example, if m(y) = Eτy is of rapid variation then Mt / m-1(t) →w 1 as t → ∞, where m-1 is the inverse function of m, while if m(y) is of regular variation with index a ∈ (0, ∞) and Xt is ergodic, then Mt / m-1(t) converges weakly to a Fréchet distribution with exponent a. In several special cases we provide explicit formulae.

In a stationary ergodic process, clustering is defined as the tendency of events to appear in series of increasedfrequency separated by longer breaks. Such behavior, contradicting the theoretical “unbiased behavior” with exponential distribution of the gaps between appearances, is commonly observed in experimental processes and often difficult to explain. In the last section we relate one such empirical example of clustering, in the area of marine technology. In the theoretical part of the paper we prove, using ergodic theory and the notion of category, that clustering (even very strong) is in fact typical for “rare events” defined as long cylinder sets in processes generated by a finite partition of an arbitrary (infinite aperiodic) ergodic measure preserving transformation.

Geometric convergence to 0 of the probability that the goal has not been encountered by the nth generation is established for a class of genetic algorithms. These algorithms employ a quickly decreasing mutation rate and a crossover which restarts the algorithm in a controlled way depending on the current population and restricts execution of this crossover to occasions when progress of the algorithm is too slow. It is shown that without the crossover studied here, which amounts to a tempered restart of the algorithm, the asserted geometric convergence need not hold.

It is well established that resource variability generated by spatial patchiness and turbulence is an importantinfluence on the growth and recruitment of planktonic fish larvae. Empirical data show fractal-like prey distributions, and simulationsindicate that scale-invariant foraging strategies may be optimal. Here we show how larval growth and recruitment in a turbulent environment canbe formulated as a hitting time problem for a jump-diffusion process. We present two theoretical results. Firstly, if jumps are of a fixed sizeand occur as a Poisson process (embedded within a drift-diffusion), recruitment is effectively described by a diffusion process alone.Secondly, in the absence of diffusion, and for “patchy” jumps (of negative binomial size with Pareto inter-arrivals), the encounter processbecomes superdiffusive. To synthesise these results we conduct a strategic simulation study where “patchy” jumps are embedded in adrift-diffusion process. We conclude that increasingly Lévy-like predator foraging strategies can have a significantly positive effect onrecruitment at the population level.

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 .


We consider a birth, death and catastrophe process where the transition rates are allowed to depend on the population size. We obtain an explicit expression for the expected time to extinction, which is valid in all cases where extinction occurs with probability 1.

Consider N players, respectively owning x1, x2, …, xN monetary units, who play a sequence of games, winning from and losing to each other integer amounts according to fixed rules. The sequence stops as soon as (at least) one player is ruined. We are interested in the ruin process of these N players, i.e. in the probability that a given player is ruined first, and also in the expected ruin time. This problem is called the N-player ruin problem. In this paper, the problem is set up as a multivariate absorbing Markov chain with an absorbing state corresponding to the ruin of each player. This is then discussed in the context of phase-type distributions where each phase is represented by a vector of size N and the distribution has as many absorbing points as there are ruin events. We use this modified phase-type distribution to obtain an explicit solution to the N-player problem. We define a partition of the set of transient states into different levels, and on it give an extension of the folding algorithm (see Ye and Li (1994)). This provides an efficient computational procedure for calculating some of the key measures.

The problem of computing the moment generating function of the first passage time T to a > 0 or −b < 0 for a one-dimensional Wiener process {X(t), t ≥ 0} is generalized by assuming that the infinitesimal parameters of the process may depend on the sign of X(t). The probability that the process is absorbed at a is also computed explicitly, as is the expected value of T.

We consider growth-collapse processes (GCPs) that grow linearly between random partial collapse times, at which they jump down according to some distribution depending on their current level. The jump occurrences are governed by a state-dependent rate function r(x). We deal with the stationary distribution of such a GCP, (Xt)t≥0, and the distributions of the hitting times Ta = inf{t ≥ 0 : Xt = a}, a > 0. After presenting the general theory of these GCPs, several important special cases are studied. We also take a brief look at the Markov-modulated case. In particular, we present a method of computing the distribution of min[Ta, σ] in this case (where σ is the time of the first jump), and apply it to determine the long-run average cost of running a certain Markov-modulated disaster-ridden system.

A cut-off phenomenon is shown to occur in a sample of n independent, identically distributed Ornstein-Uhlenbeck processes and its average. Their distributions stay far from equilibrium before a certain O(log(n)) time, and converge exponentially fast after. Precise estimates show that the total variation distance drops from almost 1 to almost 0 over an interval of time of length O(1) around log(n)/(2α), where α is the viscosity coefficient of the sampled process. The distribution of the hitting time of 0 by the average of the sample is computed. As n tends to infinity, the hitting time becomes concentrated around the cut-off instant, and its tails match the estimates given for the total variation distance.

The eigentime identity is proved for continuous-time reversible Markov chains with Markov generator L. When the essential spectrum is empty, let {0 = λ0 < λ1 ≤ λ2 ≤ ···} be the whole spectrum of L in L2. Then ∑ n≥1 λ n -1 < ∞ implies the existence of the spectral gap α of L in L∞. Explicit formulae are presented in the case of birth–death processes and from these formulae it is proved that ∑ n≥1 λ n -1 < ∞ if and only if α > 0.

The birth, death and catastrophe process is an extension of the birth–death process that incorporates the possibility of reductions in population of arbitrary size. We will consider a general form of this model in which the transition rates are allowed to depend on the current population size in an arbitrary manner. The linear case, where the transition rates are proportional to current population size, has been studied extensively. In particular, extinction probabilities, the expected time to extinction, and the distribution of the population size conditional on nonextinction (the quasi-stationary distribution) have all been evaluated explicitly. However, whilst these characteristics are of interest in the modelling and management of populations, processes with linear rate coefficients represent only a very limited class of models. We address this limitation by allowing for a wider range of catastrophic events. Despite this generalisation, explicit expressions can still be found for the expected extinction times.