In r-neighbour bootstrap percolation on the vertex set of a graph G, a set A of initially infected vertices spreads by infecting, at each time step, all uninfected vertices with at least r previously infected neighbours. When the elements of A are chosen independently with some probability p, it is natural to study the critical probability pc(G, r) at which it becomes likely that all of V(G) will eventually become infected. Improving a result of Balogh, Bollobás and Morris, we give a bound on the second term in the expansion of the critical probability when G = [n]d and d ⩾ r ⩾ 2. We show that for all d ⩾ r ⩾ 2 there exists a constant cd,r > 0 such that if n is sufficiently large, then

$$p_c (\left[ n \right]^d ,{\rm{ }}r){\rm{\le }}\left( {\frac{{\lambda (d,r)}}{{\log _{(r - 1)} (n)}} - \frac{{c_{d,r} }}{{(\log _{(r - 1)} (n))^{3/2} }}} \right)^{d - r + 1} ,$$

where λ(d, r) is an exact constant and log(k) (n) denotes the k-times iterated natural logarithm of n.

Martingales constitute a basic tool in stochastic analysis; this paper considers their application to counting processes. We use this tool to revisit a renewal theorem and give extensions for various counting processes. We first consider a renewal process as a pilot example, deriving a new semimartingale representation that differs from the standard decomposition via the stochastic intensity function. We then revisit Blackwell’s renewal theorem, its refinements and extensions. Based on these observations, we extend the semimartingale representation to a general counting process, and give conditions under which asymptotic behaviour similar to Blackwell’s renewal theorem holds.

In this paper we treat a two-stage grouping procedure of building a k-out-of-n system from several clusters of components. We use a static framework in which the component reliabilities are fixed. Under such a framework, we address the impact of the selecting strategies, the sampling probabilities, and the component reliabilities on the constructed system’s reliability. An interesting finding is that the level of component reliabilities could be identified as a decisive factor in determining how the selecting strategies and the component reliabilities affect the system reliability. The new results generalize and extend those established earlier in the literature such as Di Crescenzo and Pellerey (2011), Hazra and Nanda (2014), Navarro, Pellerey, and Di Crescenzo (2015), and Hazra, Finkelstein, and Cha (2017). Several Monte Carlo simulation experiments are provided to illustrate the theoretical results.

A large and sparse random graph with independent exponentially distributed link weights can be used to model the propagation of messages or diseases in a network with an unknown connectivity structure. In this article we study an extended setting where, in addition, the nodes of the graph are equipped with nonnegative random weights which are used to model the effect of boundary delays across paths in the network. Our main results provide approximative formulas for typical first passage times, typical flooding times, and maximum flooding times in the extended setting, over a time scale logarithmic with respect to the network size.

Many wear processes used for modeling accumulative deterioration in a reliability context are nonhomogeneous Lévy processes and, hence, have independent increments, which may not be suitable in an application context. In this work we consider Lévy processes transformed by monotonous functions to overcome this restriction, and provide a new state-dependent wear model. These transformed Lévy processes are first observed to remain tractable Markov processes. Some distributional properties are derived. We investigate the impact of the current state on the future increment level and on the overall accumulated level from a stochastic monotonicity point of view. We also study positive dependence properties and stochastic monotonicity of increments.

This paper considers two variants of M/G/1 queues with impatient customers, which are denoted by M/G/1+Gw and M/G/1+Gs. In the M/G/1+Gw queue customers have deadlines for their waiting times, and they leave the system immediately if their services do not start before the expiration of their deadlines. On the other hand, in the M/G/1+Gs queue customers have deadlines for their sojourn times, where customers in service also immediately leave the system when their deadlines expire. In this paper we derive comparison results for performance measures of these models. In particular, we show that if the service time distribution is new better than used in expectation, then the loss probability in the M/G/1+Gs queue is greater than that in the M/G/1+Gw queue.

We couple a multi-type stochastic epidemic process with a directed random graph, where edges have random weights (traversal times). This random graph representation is used to characterise the fractions of individuals infected by the different types of vertices among all infected individuals in the large population limit. For this characterisation, we rely on the theory of multi-type real-time branching processes. We identify a special case of the two-type model in which the fraction of individuals of a certain type infected by individuals of the same type is maximised among all two-type epidemics approximated by branching processes with the same mean offspring matrix.

We analyse the jigsaw percolation process, which may be seen as a measure of whether two graphs on the same vertex set are ‘jointly connected’. Bollobás, Riordan, Slivken, and Smith (2017) proved that, when the two graphs are independent binomial random graphs, whether the jigsaw process percolates undergoes a phase transition when the product of the two probabilities is $\Theta({1}/{(n\ln n)})$. We show that this threshold is sharp, and that it lies at ${1}/{(4n\ln n)}$.

Based on a simple object, an i.i.d. sequence of positive integer-valued random variables {an}n∊ℤ, we introduce and study two random structures and their connections. First, a population dynamics, in which each individual is born at time n and dies at time n + an. This dynamics is that of a D/GI/∞ queue, with arrivals at integer times and service times given by {an}n∊ℤ. Second, the directed random graph Tf on ℤ generated by the random map f(n) = n + an. Assuming only that E [a0] < ∞ and P [a0 = 1] > 0, we show that, in steady state, the population dynamics is regenerative, with one individual alive at each regeneration epoch. We identify a unimodular structure in this dynamics. More precisely, Tf is a unimodular directed tree, in which f(n) is the parent of n. This tree has a unique bi-infinite path. Moreover, Tf splits the integers into two categories: ephemeral integers, with a finite number of descendants of all degrees, and successful integers, with an infinite number. Each regeneration epoch is a successful individual such that all integers less than it are its descendants of some order. Ephemeral, successful, and regeneration integers form stationary and mixing point processes on ℤ.

The nonstationary Erlang-A queue is a fundamental queueing model that is used to describe the dynamic behavior of large-scale multiserver service systems that may experience customer abandonments, such as call centers, hospitals, and urban mobility systems. In this paper we develop novel approximations to all of its transient and steady state moments, the moment generating function, and the cumulant generating function. We also provide precise bounds for the difference of our approximations and the true model. More importantly, we show that our approximations have explicit stochastic representations as shifted Poisson random variables. Moreover, we are also able to show that our approximations and bounds also hold for nonstationary Erlang-B and Erlang-C queueing models under certain stability conditions.

We study negative association for mixed sampled point processes and show that negative association holds for such processes if a random number of their points fulfils the ultra log-concave (ULC) property. We connect the negative association property of point processes with directionally convex dependence ordering, and show some consequences of this property for mixed sampled and determinantal point processes. Some applications illustrate the general theory.

We present the first algorithm that samples maxn≥0{Sn − nα}, where Sn is a mean zero random walk, and nα with $\alpha \in ({1 \over 2},1)$ defines a nonlinear boundary. We show that our algorithm has finite expected running time. We also apply this algorithm to construct the first exact simulation method for the steady-state departure process of a GI/GI/∞ queue where the service time distribution has infinite mean.

We consider a generalised Vervaat perpetuity of the form X = Y1W1 +Y2W1W2 + · · ·, where $W_i \sim {\cal U}^{1/t}$ and (Yi)i≥0 is an independent and identically distributed sequence of random variables independent from (Wi)i≥0. Based on a distributional decomposition technique, we propose a novel method for exactly simulating the generalised Vervaat perpetuity. The general framework relies on the exact simulation of the truncated gamma process, which we develop using a marked renewal representation for its paths. Furthermore, a special case arises when Yi = 1, and X has the generalised Dickman distribution, for which we present an exact simulation algorithm using the marked renewal approach. In particular, this new algorithm is much faster than existing algorithms illustrated in Chi (2012), Cloud and Huber (2017), Devroye and Fawzi (2010), and Fill and Huber (2010), as well as being applicable to the general payments case. Examples and numerical analysis are provided to demonstrate the accuracy and effectiveness of our method.

We consider the so-called frog model with random initial configurations. The dynamics of this model are described as follows. Some particles are randomly assigned to any site of the multidimensional cubic lattice. Initially, only particles at the origin are active and these independently perform simple random walks. The other particles are sleeping and do not move at first. When sleeping particles are hit by an active particle, they become active and start moving in a similar fashion. The aim of this paper is to derive large deviation and concentration bounds for the first passage time at which an active particle reaches a target site.

Inspired by a PDE–ODE system of aggregation developed in the biomathematical literature, we investigate an interacting particle system representing aggregation at the level of individuals. We prove that the empirical density of the individual converges to the solution of the PDE–ODE system.

In large storage systems, files are often coded across several servers to improve reliability and retrieval speed. We study load balancing under the batch sampling routeing scheme for a network of n servers storing a set of files using the maximum distance separable (MDS) code (cf. Li (2016)). Specifically, each file is stored in equally sized pieces across L servers such that any k pieces can reconstruct the original file. When a request for a file is received, the dispatcher routes the job into the k-shortest queues among the L for which the corresponding server contains a piece of the file being requested. We establish a law of large numbers and a central limit theorem as the system becomes large (i.e. n → ∞), for the setting where all interarrival and service times are exponentially distributed. For the central limit theorem, the limit process take values in ℓ2, the space of square summable sequences. Due to the large size of such systems, a direct analysis of the n-server system is frequently intractable. The law of large numbers and diffusion approximations established in this work provide practical tools with which to perform such analysis. The power-of-d routeing scheme, also known as the supermarket model, is a special case of the model considered here.

We study the percolation model on Boltzmann triangulations using a generating function approach. More precisely, we consider a Boltzmann model on the set of finite planar triangulations, together with a percolation configuration (either site-percolation or bond-percolation) on this triangulation. By enumerating triangulations with boundaries according to both the boundary length and the number of vertices/edges on the boundary, we are able to identify a phase transition for the geometry of the origin cluster. For instance, we show that the probability that a percolation interface has length $n$ decays exponentially with $n$ except at a particular value $p_{c}$ of the percolation parameter $p$ for which the decay is polynomial (of order $n^{-10/3}$). Moreover, the probability that the origin cluster has size $n$ decays exponentially if $p<p_{c}$ and polynomially if $p\geqslant p_{c}$.

The critical percolation value is $p_{c}=1/2$ for site percolation, and $p_{c}=(2\sqrt{3}-1)/11$ for bond percolation. These values coincide with critical percolation thresholds for infinite triangulations identified by Angel for site-percolation, and by Angel and Curien for bond-percolation, and we give an independent derivation of these percolation thresholds.

Lastly, we revisit the criticality conditions for random Boltzmann maps, and argue that at $p_{c}$, the percolation clusters conditioned to have size $n$ should converge toward the stable map of parameter $\frac{7}{6}$ introduced by Le Gall and Miermont. This enables us to derive heuristically some new critical exponents.

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 discuss percolation and random walks in a class of homogeneous ultrametric spaces together with similarities and differences in ultrametric and Euclidean spaces. We briefly outline the role of these models in the study of interacting systems. Several open problems are presented.

We examine a system of interacting random walks with leftward drift on ℤ, which begins with a single active particle at the origin and some distribution of inactive particles on the positive integers. Inactive particles become activated when landed on by other particles, and all particles beginning at the same point possess equal leftward drift. Once activated, the trajectories of distinct particles are independent. This system belongs to a broader class of problems involving interacting random walks on rooted graphs, referred to collectively as the frog model. Additional conditions that we impose on our model include that the number of frogs (i.e. particles) at positive integer points is a sequence of independent random variables which is increasing in terms of the standard stochastic order, and that the sequence of leftward drifts associated with frogs originating at these points is decreasing. Our results include sharp conditions with respect to the sequence of random variables and the sequence of drifts that determine whether the model is transient (meaning the probability infinitely many frogs return to the origin is 0) or nontransient. We consider several, more specific, versions of the model described, and a cleaner, more simplified set of sharp conditions will be established for each case.