This paper considers the family of invariant measures of Markovian mean-field interacting particle systems on a countably infinite state space and studies its large deviation asymptotics. The Freidlin–Wentzell quasipotential is the usual candidate rate function for the sequence of invariant measures indexed by the number of particles. The paper provides two counterexamples where the quasipotential is not the rate function. The quasipotential arises from finite-horizon considerations. However, there are certain barriers that cannot be surmounted easily in any finite time horizon, but these barriers can be crossed in the stationary regime. Consequently, the quasipotential is infinite at some points where the rate function is finite. After highlighting this phenomenon, the paper studies some sufficient conditions on a class of interacting particle systems under which one can continue to assert that the Freidlin–Wentzell quasipotential is indeed the rate function.

We show that for arrival processes, the ‘harmonic new better than used in expectation’ (HNBUE) (or ‘harmonic new worse than used in expectation’, HNWUE) property is a sufficient condition for inequalities between the time and customer averages of the system if the state of the system between arrival epochs is stochastically decreasing and convex and the lack of anticipation assumption is satisfied. HNB(W)UE is a wider class than NB(W)UE, being the largest of all available classes of distributions with positive (negative) aging properties. Thus, this result represents an important step beyond existing result on inequalities between time and customer averages, which states that for arrival processes, the NB(W)UE property is a sufficient condition for inequalities.

This paper studies the input-queued switch operating under the MaxWeight algorithm when the arrivals are according to a Markovian process. We exactly characterize the heavy-traffic scaled mean sum queue length in the heavy-traffic limit, and show that it is within a factor of less than 2 from a universal lower bound. Moreover, we obtain lower and upper bounds that are applicable in all traffic regimes and become tight in the heavy-traffic regime.

We obtain these results by generalizing the drift method recently developed for the case of independent and identically distributed arrivals to the case of Markovian arrivals. We illustrate this generalization by first obtaining the heavy-traffic mean queue length and its distribution in a single-server queue under Markovian arrivals and then applying it to the case of an input-queued switch.

The key idea is to exploit the geometric mixing of finite-state Markov chains, and to work with a time horizon that is chosen so that the error due to mixing depends on the heavy-traffic parameter.

We calculate the mean throughput, number of collisions, successes, and idle slots for random tree algorithms with successive interference cancellation. Except for the case of the throughput for the binary tree, all the results are new. We furthermore disprove the claim that only the binary tree maximizes throughput. Our method works with many observables and can be used as a blueprint for further analysis.

We study heterogeneously interacting diffusive particle systems with mean-field-type interaction characterized by an underlying graphon and their finite particle approximations. Under suitable conditions, we obtain exponential concentration estimates over a finite time horizon for both 1- and 2-Wasserstein distances between the empirical measures of the finite particle systems and the averaged law of the graphon system.

This article considers the individual equilibrium behavior and socially optimal strategy in a fluid queue with two types of parallel customers and incomplete fault. Assume that the working state and the incomplete fault state appear alternately in the buffer. Different from the linear revenue and expenditure structure, an exponential utility function can be constructed to obtain the equilibrium balking thresholds in the fully observable case. Besides, the steady-state probability distribution and the corresponding expected social benefit are derived based on the renewal process and the standard theory of linear ordinary differential equations. Furthermore, a reasonable entrance fee strategy is discussed under the condition that the fluid accepts the globally optimal strategies. Finally, the effects of the diverse system parameters on the entrance fee and the expected social benefit are explicitly illustrated by numerical comparisons.

We study 2-stage game-theoretic problem oriented 3-stage service policy computing, convolutional neural network (CNN) based algorithm design, and simulation for a blockchained buffering system with federated learning. More precisely, based on the game-theoretic problem consisting of both “win-lose” and “win-win” 2-stage competitions, we derive a 3-stage dynamical service policy via a saddle point to a zero-sum game problem and a Nash equilibrium point to a non-zero-sum game problem. This policy is concerning users-selection, dynamic pricing, and online rate resource allocation via stable digital currency for the system. The main focus is on the design and analysis of the joint 3-stage service policy for given queue/environment state dependent pricing and utility functions. The asymptotic optimality and fairness of this dynamic service policy is justified by diffusion modeling with approximation theory. A general CNN based policy computing algorithm flow chart along the line of the so-called big model framework is presented. Simulation case studies are conducted for the system with three users, where only two of the three users can be selected into the service by a zero-sum dual cost game competition policy at a time point. Then, the selected two users get into service and share the system rate service resource through a non-zero-sum dual cost game competition policy. Applications of our policy in the future blockchain based Internet (e.g., metaverse and web3.0) and supply chain finance are also briefly illustrated.

We consider bond percolation on high-dimensional product graphs $G=\square _{i=1}^tG^{(i)}$, where $\square$ denotes the Cartesian product. We call the $G^{(i)}$ the base graphs and the product graph $G$ the host graph. Very recently, Lichev (J. Graph Theory, 99(4):651–670, 2022) showed that, under a mild requirement on the isoperimetric properties of the base graphs, the component structure of the percolated graph $G_p$ undergoes a phase transition when $p$ is around $\frac{1}{d}$, where $d$ is the average degree of the host graph.

In the supercritical regime, we strengthen Lichev’s result by showing that the giant component is in fact unique, with all other components of order $o(|G|)$, and determining the sharp asymptotic order of the giant. Furthermore, we answer two questions posed by Lichev (J. Graph Theory, 99(4):651–670, 2022): firstly, we provide a construction showing that the requirement of bounded degree is necessary for the likely emergence of a linear order component; secondly, we show that the isoperimetric requirement on the base graphs can be, in fact, super-exponentially small in the dimension. Finally, in the subcritical regime, we give an example showing that in the case of irregular high-dimensional product graphs, there can be a polynomially large component with high probability, very much unlike the quantitative behaviour seen in the Erdős-Rényi random graph and in the percolated hypercube, and in fact in any regular high-dimensional product graphs, as shown by the authors in a companion paper (Percolation on high-dimensional product graphs. arXiv:2209.03722, 2022).

We study three classes of shock models governed by an inverse gamma mixed Poisson process (IGMP), namely a mixed Poisson process with an inverse gamma mixing distribution. In particular, we analyze (1) the extreme shock model, (2) the δ-shock model, and the (3) cumulative shock model. For the latter, we assume a constant and an exponentially distributed random threshold and consider different choices for the distribution of the amount of damage caused by a single shock. For all the treated cases, we obtain the survival function, together with the expected value and the variance of the failure time. Some properties of the inverse gamma mixed Poisson process are also disclosed.

We continue with the systematic study of the speed of extinction of continuous-state branching processes in Lévy environments under more general branching mechanisms. Here, we deal with the weakly subcritical regime under the assumption that the branching mechanism is regularly varying. We extend recent results of Li and Xu (2018) and Palau et al. (2016), where it is assumed that the branching mechanism is stable, and complement the recent articles of Bansaye et al. (2021) and Cardona-Tobón and Pardo (2021), where the critical and the strongly and intermediate subcritical cases were treated, respectively. Our methodology combines a path analysis of the branching process together with its Lévy environment, fluctuation theory for Lévy processes, and the asymptotic behaviour of exponential functionals of Lévy processes. Our approach is inspired by the last two previously cited papers, and by Afanasyev et al. (2012), where the analogue was obtained.

The purpose of this study is to present a subgeometric convergence formula for the stationary distribution of the finite-level M/G/1-type Markov chain when taking its infinite-level limit, where the upper boundary level goes to infinity. This study is carried out using the fundamental deviation matrix, which is a block-decomposition-friendly solution to the Poisson equation of the deviation matrix. The fundamental deviation matrix provides a difference formula for the respective stationary distributions of the finite-level chain and the corresponding infinite-level chain. The difference formula plays a crucial role in the derivation of the main result of this paper, and the main result is used, for example, to derive an asymptotic formula for the loss probability in the MAP/GI/1/N queue.

We determine almost sure limits of rescaled intrinsic volumes of the construction steps of fractal percolation in ${\mathbb R}^d$ for any dimension $d\geq 1$. We observe a factorization of these limit variables which allows one, in particular, to determine their expectations and covariance structure. We also show the convergence of the rescaled expectations and variances of the intrinsic volumes of the construction steps to the expectations and variances of the limit variables, and we give rates for this convergence in some cases. These results significantly extend our previous work, which addressed only limits of expectations of intrinsic volumes.

A spread-out lattice animal is a finite connected set of edges in $\{\{x,y\}\subset \mathbb{Z}^d\;:\;0\lt \|x-y\|\le L\}$. A lattice tree is a lattice animal with no loops. The best estimate on the critical point $p_{\textrm{c}}$ so far was achieved by Penrose (J. Stat. Phys. 77, 3–15, 1994) : $p_{\textrm{c}}=1/e+O(L^{-2d/7}\log L)$ for both models for all $d\ge 1$. In this paper, we show that $p_{\textrm{c}}=1/e+CL^{-d}+O(L^{-d-1})$ for all $d\gt 8$, where the model-dependent constant $C$ has the random-walk representation

\begin{align*} C_{\textrm{LT}}=\sum _{n=2}^\infty \frac{n+1}{2e}U^{*n}(o),&& C_{\textrm{LA}}=C_{\textrm{LT}}-\frac 1{2e^2}\sum _{n=3}^\infty U^{*n}(o), \end{align*}

$U^{*n}$

$n$

$d$

$\{x\in{\mathbb R}^d\;: \|x\|\le 1\}$

$p$

The switch process alternates independently between 1 and $-1$, with the first switch to 1 occurring at the origin. The expected value function of this process is defined uniquely by the distribution of switching times. The relation between the two is implicitly described through the Laplace transform, which is difficult to use for determining if a given function is the expected value function of some switch process. We derive an explicit relation under the assumption of monotonicity of the expected value function. It is shown that geometric divisible switching time distributions correspond to a non-negative decreasing expected value function. Moreover, an explicit relation between the expected value of a switch process and the autocovariance function of the switch process stationary counterpart is obtained, leading to a new interpretation of the classical Pólya criterion for positive-definiteness.

We study the asymptotic limit of random pure dimer coverings on rail yard graphs when the mesh sizes of the graphs go to 0. Each pure dimer covering corresponds to a sequence of interlacing partitions starting with an empty partition and ending in an empty partition. Under the assumption that the probability of each dimer covering is proportional to the product of weights of present edges, we obtain the limit shape (law of large numbers) of the rescaled height functions and the convergence of the unrescaled height fluctuations to a diffeomorphic image of the Gaussian free field (Central Limit Theorem), answering a question in [7]. Applications include the limit shape and height fluctuations for pure steep tilings [9] and pyramid partitions [20; 36; 39; 38]. The technique to obtain these results is to analyze a class of Macdonald processes which involve dual partitions as well.

We study the $R_\beta$-positivity and the existence of zero-temperature limits for a sequence of infinite-volume Gibbs measures $(\mu_{\beta}(\!\cdot\!))_{\beta \geq 0}$ at inverse temperature $\beta$ associated to a family of nearest-neighbor matrices $(Q_{\beta})_{\beta \geq 0}$ reflected at the origin. We use a probabilistic approach based on the continued fraction theory previously introduced in Ferrari and Martínez (1993) and sharpened in Littin and Martínez (2010). Some necessary and sufficient conditions are provided to ensure (i) the existence of a unique infinite-volume Gibbs measure for large but finite values of $\beta$, and (ii) the existence of weak limits as $\beta \to \infty$. Some application examples are revised to put in context the main results of this work.

In 2008, Tóth and Vető defined the self-repelling random walk with directed edges as a non-Markovian random walk on $\unicode{x2124}$: in this model, the probability that the walk moves from a point of $\unicode{x2124}$ to a given neighbor depends on the number of previous crossings of the directed edge from the initial point to the target, called the local time of the edge. Tóth and Vető found that this model exhibited very peculiar behavior, as the process formed by the local times of all the edges, evaluated at a stopping time of a certain type and suitably renormalized, converges to a deterministic process, instead of a random one as in similar models. In this work, we study the fluctuations of the local times process around its deterministic limit, about which nothing was previously known. We prove that these fluctuations converge in the Skorokhod $M_1$ topology, as well as in the uniform topology away from the discontinuities of the limit, but not in the most classical Skorokhod topology. We also prove the convergence of the fluctuations of the aforementioned stopping times.

A comparison theorem for state-dependent regime-switching diffusion processes is established, which enables us to pathwise-control the evolution of the state-dependent switching component simply by Markov chains. Moreover, a sharp estimate on the stability of Markovian regime-switching processes under the perturbation of transition rate matrices is provided. Our approach is based on elaborate constructions of switching processes in the spirit of Skorokhod’s representation theorem varying according to the problem being dealt with. In particular, this method can cope with switching processes in an infinite state space and not necessarily of birth–death type. As an application, some known results on the ergodicity and stability of state-dependent regime-switching processes can be improved.

We consider a simple random walk on $\mathbb{Z}^d$ started at the origin and stopped on its first exit time from $({-}L,L)^d \cap \mathbb{Z}^d$. Write L in the form $L = m N$ with $m = m(N)$ and N an integer going to infinity in such a way that $L^2 \sim A N^d$ for some real constant $A \gt 0$. Our main result is that for $d \ge 3$, the projection of the stopped trajectory to the N-torus locally converges, away from the origin, to an interlacement process at level $A d \sigma_1$, where $\sigma_1$ is the exit time of a Brownian motion from the unit cube $({-}1,1)^d$ that is independent of the interlacement process. The above problem is a variation on results of Windisch (2008) and Sznitman (2009).