We consider parallel single-server queues in heavy traffic with randomly split Hawkes arrival processes. The service times are assumed to be independent and identically distributed (i.i.d.) in each queue and are independent in different queues. In the critically loaded regime at each queue, it is shown that the diffusion-scaled queueing and workload processes converge to a multidimensional reflected Brownian motion in the non-negative orthant with orthonormal reflections. For the model with abandonment, we also show that the corresponding limit is a multidimensional reflected Ornstein–Uhlenbeck diffusion in the non-negative orthant.

We consider the random splitting and aggregating of Hawkes processes. We present the random splitting schemes using the direct approach for counting processes, as well as the immigration–birth branching representations of Hawkes processes. From the second scheme, it is shown that random split Hawkes processes are again Hawkes. We discuss functional central limit theorems (FCLTs) for the scaled split processes from the different schemes. On the other hand, aggregating multivariate Hawkes processes may not necessarily be Hawkes. We identify a necessary and sufficient condition for the aggregated process to be Hawkes. We prove an FCLT for a multivariate Hawkes process under a random splitting and then aggregating scheme (under certain conditions, transforming into a Hawkes process of a different dimension).

This study presents functional limit theorems for the Euler characteristic of Vietoris–Rips complexes. The points are drawn from a nonhomogeneous Poisson process on $\mathbb{R}^d$ , and the connectivity radius governing the formation of simplices is taken as a function of the time parameter t, which allows us to treat the Euler characteristic as a stochastic process. The setting in which this takes place is that of the critical regime, in which the simplicial complexes are highly connected and have nontrivial topology. We establish two ‘functional-level’ limit theorems, a strong law of large numbers and a central limit theorem, for the appropriately normalized Euler characteristic process.

In this paper we prove a quenched functional central limit theorem for a biased random walk on a supercritical Galton–Watson tree with leaves. This extends a result of Peres and Zeitouni (2008) where the case without leaves was considered. A conjecture of Ben Arous and Fribergh (2016) suggests an upper bound on the bias which we observe to be sharp.

Markov chain Monte Carlo (MCMC) methods provide an essential tool in statistics for sampling from complex probability distributions. While the standard approach to MCMC involves constructing discrete-time reversible Markov chains whose transition kernel is obtained via the Metropolis–Hastings algorithm, there has been recent interest in alternative schemes based on piecewise deterministic Markov processes (PDMPs). One such approach is based on the zig-zag process, introduced in Bierkens and Roberts (2016), which proved to provide a highly scalable sampling scheme for sampling in the big data regime; see Bierkens et al. (2016). In this paper we study the performance of the zig-zag sampler, focusing on the one-dimensional case. In particular, we identify conditions under which a central limit theorem holds and characterise the asymptotic variance. Moreover, we study the influence of the switching rate on the diffusivity of the zig-zag process by identifying a diffusion limit as the switching rate tends to ∞. Based on our results we compare the performance of the zig-zag sampler to existing Monte Carlo methods, both analytically and through simulations.

We consider the Δ(i)/G/1 queue, in which a total of n customers join a single-server queue for service. Customers join the queue independently after exponential times. We consider heavy-tailed service-time distributions with tails decaying as x-α, α ∈ (1, 2). We consider the asymptotic regime in which the population size grows to ∞ and establish that the scaled queue-length process converges to an α-stable process with a negative quadratic drift. We leverage this asymptotic result to characterize the head start that is needed to create a long period of uninterrupted activity (a busy period). The heavy-tailed service times should be contrasted with the case of light-tailed service times, for which a similar scaling limit arises (Bet et al. (2015)), but then with a Brownian motion instead of an α-stable process.

In this paper we investigate the functional central limit theorem (CLT) for stochastic processes associated to partial sums of additive functionals of reversible Markov chains with general spate space, under the normalization standard deviation of partial sums. For this case, we show that the functional CLT is equivalent to the fact that the variance of partial sums is regularly varying with exponent 1 and the partial sums satisfy the CLT. It is also equivalent to the conditional CLT.

Fractional Lévy processes generalize fractional Brownian motion in a natural way. We go a step further and extend the usual fractional Riemann-Liouville kernel to a regularly varying function. We call the resulting stochastic processes generalized fractional Lévy processes (GFLPs) and show that they may have short or long memory increments and that their sample paths may have jumps or not. Moreover, we define stochastic integrals with respect to a GFLP and investigate their second-order structure and sample path properties. A specific example is the Ornstein-Uhlenbeck process driven by a time-scaled GFLP. We prove a functional central limit theorem for such scaled processes with a fractional Ornstein-Uhlenbeck process as a limit process. This approximation applies to a wide class of stochastic volatility models, which include models where possibly neither the data nor the latent volatility process are semimartingales.

The Hawkes process is a self-exciting point process with clustering effect whose intensity depends on its entire past history. It has wide applications in neuroscience, finance, and many other fields. In this paper we obtain a functional central limit theorem for the nonlinear Hawkes process. Under the same assumptions, we also obtain a Strassen's invariance principle, i.e. a functional law of the iterated logarithm.

In this paper we study the functional central limit theorem (CLT) for stationary Markov chains with a self-adjoint operator and general state space. We investigate the case when the variance of the partial sum is not asymptotically linear in n, and establish that conditional convergence in distribution of partial sums implies the functional CLT. The main tools are maximal inequalities that are further exploited to derive conditions for tightness and convergence to the Brownian motion.

In this paper we study the asymptotic properties of the canonical plugin estimates for law-invariant coherent risk measures. Under rather mild conditions not relying on the explicit representation of the risk measure under consideration, we first prove a central limit theorem for independent and identically distributed data, and then extend it to the case of weakly dependent data. Finally, a number of illustrating examples is presented.

We establish heavy-traffic limits for nearly deterministic queues, such as the G/D/n many-server queue. Since waiting times before starting service in the G/D/n queue are equivalent to waiting times in an associated G n /D/1 model, where the G n interarrival times are the sum of n consecutive interarrival times in the original model, we focus on the G n /D/1 model and the generalization to G n /G n /1, where ‘cyclic thinning’ is applied to both the arrival and service processes. We establish different limits in two cases: (i) when (1 − ρ n )√n → β as n → ∞ and (ii) when (1 − ρ n )n → β as n → ∞, where ρ n is the traffic intensity in model n. The nearly deterministic feature leads to interesting nonstandard scaling.

We give a functional central limit theorem for spatial birth and death processes based on the representation of such processes as solutions of stochastic equations. For any bounded and integrable function in Euclidean space, we define a family of processes which is obtained by integrals of this function with respect to the centered and scaled spatial birth and death process with constant death rate. We prove that this family converges weakly to a Gaussian process as the scale parameter goes to infinity. We do not need the birth rates to have a finite range of interaction. Instead, we require that the birth rates have a range of interaction that decays polynomially. In order to show the convergence of the finite-dimensional distributions of the above processes, we extend Penrose's multivariate spatial central limit theorem. An example of the asymptotic normalities of the time-invariance estimators for the birth rates of spatial point processes is given.

For a linear random field (linear p-parameter stochastic process) generated by a dependent random field with zero mean and finite qth moments (q>2p), we give sufficient conditions that the linear random field converges weakly to a multiparameter standard Brownian motion if the corresponding dependent random field does so.

For a Markov renewal process where the time parameter is discrete, we present a novel method for calculating the asymptotic variance. Our approach is based on the key renewal theorem and is applicable even when the state space of the Markov chain is countably infinite.

This paper introduces a method of generating real harmonizable multifractional Lévy motions (RHMLMs). The simulation of these fields is closely related to that of infinitely divisible laws or Lévy processes. In the case where the control measure of the RHMLM is finite, generalized shot-noise series are used. An estimation of the error is also given. Otherwise, the RHMLM Xh is split into two independent RHMLMs, Xε,1 and Xε,2. More precisely, Xε,2 is an RHMLM whose control measure is finite. It can then be rewritten as a generalized shot-noise series. The asymptotic behaviour of Xε,1 as ε → 0+ is further elaborated. Sufficient conditions to approximate Xε,1 by a multifractional Brownian motion are given. The error rate in terms of Berry-Esseen bounds is then discussed. Finally, some examples of simulation are given.

Motivated by various applications in queueing systems, this work is devoted to continuous-time Markov chains with countable state spaces that involve both fast-time scale and slow-time scale with the aim of approximating the time-varying queueing systems by their quasistationary counterparts. Under smoothness conditions on the generators, asymptotic expansions of probability vectors and transition probability matrices are constructed. Uniform error bounds are obtained, and then sequences of occupation measures and their functionals are examined. Mean square error estimates of a sequence of occupation measures are obtained; a scaled sequence of functionals of occupation measures is shown to converge to a Gaussian process with zero mean. The representation of the variance of the limit process is also explicitly given. The results obtained are then applied to treat Mt/Mt/1 queues and Markov-modulated fluid buffer models.

Let X = (X(t):t ≥ 0) be a Lévy process and X∊ the compensated sum of jumps not exceeding ∊ in absolute value, σ2(∊) = var(X∊(1)). In simulation, X - X∊ is easily generated as the sum of a Brownian term and a compound Poisson one, and we investigate here when X∊/σ(∊) can be approximated by another Brownian term. A necessary and sufficient condition in terms of σ(∊) is given, and it is shown that when the condition fails, the behaviour of X∊/σ(∊) can be quite intricate. This condition is also related to the decay of terms in series expansions. We further discuss error rates in terms of Berry-Esseen bounds and Edgeworth approximations.

We establish new exponential inequalities for partial sums of random fields. Next, using classicalchaining arguments, we give sufficient conditions for partial sum processes indexed by large classes ofsets to converge to a set-indexed Brownian motion. For stationary fields of bounded random variables, thecondition is expressed in terms of a series of conditional expectations. For non-uniform ϕ-mixingrandom fields, we require both finite fourth moments and an algebraic decay of the mixing coefficients.

Processes that exhibit repeatability in their kth-order moments are frequently studied in signal analysis. Such repeatability can be conveniently expressed with the help of almost periodicity. In particular, almost periodically correlated (APC) processes play an important role in the analysis of repeatable signals. This paper presents a study of asymptotic distributions of the estimator of the spectral covariance function for APC processes. It is demonstrated that, for a large class of APC processes, the functional central limit theorem holds.