We use a discrete-time analysis, giving necessary and sufficient conditions for the almost-sure convergence of ARCH(1) and GARCH(1,1) discrete-time models, to suggest an extension of the ARCH and GARCH concepts to continuous-time processes. Our ‘COGARCH’ (continuous-time GARCH) model, based on a single background driving Lévy process, is different from, though related to, other continuous-time stochastic volatility models that have been proposed. The model generalises the essential features of discrete-time GARCH processes, and is amenable to further analysis, possessing useful Markovian and stationarity properties.

A continuous-time random walk is a simple random walk subordinated to a renewal process used in physics to model anomalous diffusion. In this paper we show that, when the time between renewals has infinite mean, the scaling limit is an operator Lévy motion subordinated to the hitting time process of a classical stable subordinator. Density functions for the limit process solve a fractional Cauchy problem, the generalization of a fractional partial differential equation for Hamiltonian chaos. We also establish a functional limit theorem for random walks with jumps in the strict generalized domain of attraction of a full operator stable law, which is of some independent interest.

This paper studies recurrence properties of autoregressive (AR) processes with ‘super-heavy-tailed’ innovations. Specifically, we study the case where the innovations are distributed, roughly speaking, as log-Pareto random variables (i.e. the tail decay is essentially a logarithm raised to some power). We show that these processes exhibit interesting and somewhat surprising behaviour. In particular, we show that AR(1) processes, with the usual root assumption that is necessary for stability, can exhibit null-recurrent as well as transient dynamics when the innovations follow a log-Cauchy-type distribution. In this regime, the recurrence classification of the process depends, somewhat surprisingly, on the value of the constant pre-multiplier of this distribution. More generally, for log-Pareto innovations, we provide a positive-recurrence/null-recurrence/transience classification of the corresponding AR processes.

We present a necessary and sufficient condition for a stochastic exponential to be a true martingale. It is proved that the criteria for the true martingale property are related to whether a related process explodes. An alternative and interesting interpretation of this result is that the stochastic exponential is a true martingale if and only if under a ‘candidate measure’ the integrand process is square integrable over time. Applications of our theorem to problems arising in mathematical finance are also given.

We study a one-dimensional telegraph process (M t ) t≥0 describing the position of a particle moving at constant speed between Poisson times at which new velocities are chosen randomly. The exact distribution of M t and its first two moments are derived. We characterize the level hitting times of M t in terms of integro-differential equations which can be solved in special cases.

Let C 1, C 2,…,C m be independent subordinators with finite expectations and denote their sum by C. Consider the classical risk process X(t) = x + ct - C(t). The ruin probability is given by the well-known Pollaczek–Khinchin formula. If ruin occurs, however, it will be caused by a jump of one of the subordinators C i . Formulae for the probability that ruin is caused by C i are derived. These formulae can be extended to perturbed risk processes of the type X(t) = x + ct - C(t) + Z(t), where Z is a Lévy process with mean 0 and no positive jumps.

For a given k ≥ 1, subintervals of a given interval [0, X] arrive at random and are accepted (allocated) so long as they overlap fewer than k subintervals already accepted. Subintervals not accepted are cleared, while accepted subintervals remain allocated for random retention times before they are released and made available to subsequent arrivals. Thus, the system operates as a generalized many-server queue under a loss protocol. We study a discretized version of this model that appears in reference theories for a number of applications, including communication networks, surface adsorption-desorption processes, and reservation systems. Our primary interest is in steady-state estimates of the vacant space, i.e. the total length of available subintervals kX - ∑ℓ i , where the ℓ i are the lengths of the subintervals currently allocated. We obtain explicit results for k = 1 and for general k with all subinterval lengths equal to 2, the classical dimer case of chemical applications. Our focus is on the asymptotic regime of large retention times.

Large deviation results are obtained for the normed limit of a supercritical multitype branching process. Starting from a single individual of type i, let L[i] be the normed limit of the branching process, and let be the minimum possible population size at generation k. If is bounded in k (bounded minimum growth), then we show that P(L[i] ≤ x) = P(L[i] = 0) + x α F *[i](x) + o(x α ) as x → 0. If grows exponentially in k (exponential minimum growth), then we show that −log P(L[i] ≤ x) = x −β/(1−β) G*[i](x) + o (x −β/(1−β)) as x → 0. If the maximum family size is bounded, then −log P(L[i] > x) = x δ/(δ−1) H *[i](x) + o(x δ/(δ−1)) as x → ∞. Here α, β and δ are constants obtained from combinations of the minimum, maximum and mean growth rates, and F *, G * and H * are multiplicatively periodic functions.

Caching is an effective mechanism for reducing bandwidth usage and alleviating server load. However, the use of caching entails a compromise between content freshness and refresh cost. Constant refreshing allows a high degree of content freshness at a greater cost of system resource. Conversely, too little refreshing reduces content freshness but saves the cost of resource usage. To address the freshness-cost problem, we formulate the refresh scheduling problem with a generic cost model and use this cost model to determine an optimal refresh frequency that gives the best trade-off between refresh cost and content freshness. We prove the existence and uniqueness of an optimal refresh frequency under the assumptions that the arrival of updated content is Poisson and the age-related cost monotonically increases with decreasing freshness. In addition, we provide an analytic comparison of system performance under fixed refresh scheduling and random refresh scheduling, showing that two refresh schedulings with the same average refresh frequency are mathematically equivalent in terms of the long-run average cost.

Burn-in is a widely used engineering procedure useful for eliminating ‘weak’ items and consequently improving the quality of remaining items. The quality of items can be measured via various performance characteristics. In the present paper we develop new performance criteria for the burn-in method. Our criteria not only take into account the reliability of an item, they also incorporate covariates.

We consider a fluid queue controlled by a semi-Markov process and we apply the Markov-renewal approach developed earlier in the context of quasi-birth-and-death processes and of Markovian fluid queues. We analyze two subfamilies of semi-Markov processes. In the first family, we assume that the intervals during which the input rate is negative have an exponential distribution. In the second family, we take the complementary case and assume that the intervals during which the input rate is positive have an exponential distribution. We thoroughly characterize the structure of the stationary distribution in both cases.

We consider an infinite-capacity second-order fluid queue governed by a continuous-time Markov chain and with linear service rate. The variability of the traffic is modeled by a Brownian motion and a local variance function modulated by the Markov chain and proportional to the fluid level in the queue. The behavior of this second-order fluid-flow model is described by a linear stochastic differential equation, satisfied by the transient queue level. We study the transient level's convergence in distribution under weak assumptions and we obtain an expression for the stationary queue level. For the first-order case, we give a simple expression of all its moments as well as of its Laplace transform. For the second-order model we compute its first two moments.

This paper investigates the rate of convergence to the probability distribution of the embedded M/G/1 and GI/M/n queues. We introduce several types of ergodicity including l-ergodicity, geometric ergodicity, uniformly polynomial ergodicity and strong ergodicity. The usual method to prove ergodicity of a Markov chain is to check the existence of a Foster–Lyapunov function or a drift condition, while here we analyse the generating function of the first return probability directly and obtain practical criteria. Moreover, the method can be extended to M/G/1- and GI/M/1-type Markov chains.

For a K-stage cyclic queueing network with N customers and general service times, we provide bounds on the nth departure time from each stage. Furthermore, we analyze the asymptotic tail behavior of cycle times and waiting times given that at least one service-time distribution is subexponential.

We demonstrate that the residual lifetime distribution of a compound geometric distribution convoluted with another distribution, termed a compound geometric convolution, is again a compound geometric convolution. Conditions under which the compound geometric convolution is new worse than used (NWU) or new better than used (NBU) are then derived. The results are applied to ruin probabilities in the stationary renewal risk model where the convolution components are of particular interest, as well as to the equilibrium virtual waiting time distribution in the G/G/1 queue, an approximation to the equilibrium M/G/c waiting time distribution, ruin in the classical risk model perturbed by diffusion, and second-order reliability classifications.

An annihilating process is an interacting particle system in which the only interaction is that a particle may kill a neighbouring particle. Since there is no birth and no movement, once a particle has no neighbours its site remains occupied for ever. The survival probability is calculated for a random tree and for the square lattice. A connection is made between annihilating processes and the adsorption of molecules onto surfaces. A one-dimensional adsorption problem is solved in the case in which the two neighbours do not act independently.

We derive in this paper closed formulae for the joint probability generating function of the number of customers in the two FIFO queues of a generalized processor-sharing (GPS) system with two classes of customers arriving according to Poisson processes and requiring exponential service times. In contrast to previous studies published on the GPS system, we show that it is possible to establish explicit expressions for the generating functions of the number of customers in each queue without calling for the formulation of a Riemann–Hilbert problem. We specifically prove that the problem of determining the unknown functions due to the reflecting conditions on the boundaries of the positive quarter plane can be reduced to a Poisson equation. The explicit formulae are then used to derive some characteristics of the GPS system (in particular the tails of the probability distributions of the numbers of customers in each queue).

We introduce adaptive-simulation schemes for estimating performance measures for stochastic systems based on the method of control variates. We consider several possible methods for adaptively tuning the control-variate estimators, and describe their asymptotic properties. Under certain assumptions, including the existence of a ‘perfect control variate’, all of the estimators considered converge faster than the canonical rate of n −1/2, where n is the simulation run length. Perfect control variates for a variety of stochastic processes can be constructed from ‘approximating martingales’. We prove a central limit theorem for an adaptive estimator that converges at rate A similar estimator converges at rate n −1. An exponential rate of convergence is also possible under suitable conditions.

In this article, we analyse the error induced by the Euler scheme combined with a symmetry procedure near the boundary for the simulation of diffusion processes with an oblique reflection on a smooth boundary. This procedure is easy to implement and, in addition, accurate: indeed, we prove that it yields a weak rate of convergence of order 1 with respect to the time-discretization step.

We consider random dynamical systems with randomly chosen jumps on infinite-dimensional spaces. The choice of deterministic dynamical systems and jumps depends on a position. The system generalizes dynamical systems corresponding to learning systems, Poisson driven stochastic differential equations, iterated function system with infinite family of transformations and random evolutions. We will show that distributions which describe the dynamics of this system converge to an invariant distribution. We use recent results concerning asymptotic stability of Markov operators on infinite-dimensional spaces obtained by T. Szarek.