Power law condition for stability of Poisson hail

We consider the Poisson hail model introduced by Baccelli and Foss. We give a power law condition for the tails (spatial and temporal) of the distribution of jobs to ensure stability as the rate parameter $\lambda $ tends to zero. We then show that in a weak sense it is optimal.


Introduction
The purpose of this article is to loosen conditions for stability in the "Poisson hail" interacting queueing model introduced by [1]. In the discrete setting for this model, there are countably many jobs (identified by countably many points in space-time). Job i requires service τ i from a subset B i ⊂ Z d . As in the preceding paper, we associate to a job i a (semi-arbitrary) server x = x(i) ∈ Z d who is in some sense central in the group B i . We suppose that for each site w ∈ Z d the jobs i with x(i) = w arrive according to a Poisson process N w of rate λ. The jobs i arriving at w will have their subsets B i and service times τ i distributed as i.i.d. vectors, so the arrivals at site w may be considered as a marked Poisson process w . In other words, w is a Poisson process on R × R + × 2 Z d . Points in its support are typically denoted by (t, τ, B), the t's forming the aforementioned rate-λ Poisson process. The pair (τ, B) is referred to as the mark of the point t. We also assume that the w are obtained as follows: Let, for each w, w be an independent copy of 0 . Then let w contain all points of the form (t, τ, B + w) where (t, τ, B) is a point of w . Thus, the arrival process (including marks) is translation invariant. Physically, we can think of the system as a model of hailstones of cylindrical shape B × [0, τ ] ⊂ Z d+1 , where τ is the height of the stone and B its base. When a hailstone appears at some point of time at which all sites w ∈ B are free, it starts melting at rate 1. If there is at least one w ∈ B occupied by a previously arrived stone, then the current stone will not start melting before all sites in w become free; at the first moment of time this happens, the hailstone starts melting at rate 1. (Only the ground, Z d is hot and heat is not transmitted upwards!) At each time t, we let W (t, x) be the total work required for x to be free of hailstones provided no stones arrive after t. In queueing terms, W (t, x) is a workload. In hailstone terms, W (t, x) is the sum of the heights of all hailstones which contain x in their base and have not been melted yet. Since the superposition of N w , w ∈ Z d , has infinite rate, it follows that within any time interval of positive length there are infinitely many stones arriving. Thus, W (t, ·) will change infinitely many times in any right neighborhood of t. However, typically, for fixed x ∈ Z d , and any ε > 0, W (t, x) depends only on W (t − ε, y), for y ranging in a finite (but random) number of sites. This is due to the fact that the only have to look at those w with points (s, τ, B) such that t − ε ≤ s ≤ t and x ∈ B.
Fix x ∈ Z d and suppose there is w ∈ Z d such that (t, τ, B) is a point of w . Then By convention, we shall assume that t → W (t, x) is right-continuous: W (t, x) = W (t+, x). On the other hand, if there is no w such that (t, τ, B) is a point of w with x ∈ B, then W (s, t), s ≥ t, decreases linearly for a interval of positive length until either it reaches zero or there is job arriving at some s > t at some site w whose base contains x. We have thus completely specified the dynamics of the system. The system considered here differs from that of [3] in that the latter (i) considers only finitely many sites (Z d is replaced by a finite set) but (ii) works for stationary and ergodic arrival processes.
The system is said to be stable if (starting from full vacancy at time 0) the distribution of W (t, x) is tight as t varies for fixed x. The central question to be addressed is when is the system stable (for λ sufficiently small). More precisely, for which laws on (τ, B) for jobs arriving at the origin is it the case that there exists λ 0 ∈ (0, ∞) so that the system is stable for all arrival rates λ < λ 0 . To avoid trivialities, we assume that B is a finite set, a.s. The founding article [1] showed that the system was indeed stable provided that there is c ∈ (0, ∞) so that where diam B is the diameter of set B, i.e., the maximum of |x − y| ∞ over all x, y ∈ B, and where |x| ∞ := max 1≤i≤d |x i |. The proof in [1] is based on a comparison with an auxiliary branching process with weights requiring the condition stated in the last display.
Our purpose in this paper is to slacken this condition to the existence of the (d +1+ ε)th moment for τ + diam B. We then (easily) show that this condition is (in a certain weak sense) almost optimal. The key idea is to use ideas on laws of large numbers for lattice animals. This was first proved in [2], though for this paper we take as reference the article by James Martin [4]. In analogy to [3], one could also ask whether stability is possible for more general arrival processes. This question, however, is outside the scope of our paper as our method explicitly uses the Poissonian assumptions.
Our principal result is Then there exists λ 0 > 0 so that for job arrival rate λ < λ 0 the system is stable.
That this result is (in a weak sense) the best possible is shown by Theorem 2 For any d +1 > ε > 0, we can find a (spatially homogeneous) job arrival process so that and the system is unstable.

Remark 1
The condition of Theorem 1 is equivalent to the following: There exists ε > 0 and C > 0 such that, for any x ≥ 0, Similarly, the condition of Theorem 2 is equivalent to the same thing with −ε in place of ε.
Given stability, it is easy to see that starting from complete vacancy (that is, no workload at any site), the system converges in distribution to an explicitly describable equilibrium. It is natural to ask whether the system possesses other, not necessarily spatially homogeneous, equilibria. While not definitively answering this we show, Theorem 3 Under the conditions of Theorem 1, there exists λ 0 > 0 so that for arrival rate 0 < λ < λ 0 , the only equilibrium for the system that is spatially translation invariant is the limit measure obtained by starting from zero workload.
We now assemble some observations and techniques from earlier papers, [1,3]. Start the system at time −n from full vacancy and consider how the workload W n (t, x) at time t ≥ −n and site x ∈ Z d is obtained.
Definition 1 Let n (x, t) be the set of locally constant cadlag (= piecewise constant, continuous on the right with left limits at every point) paths γ : Associate to such a γ ∈ n (x, t) the score where the sum is over jobs (τ i , B i ) which are arrive at time s i with γ (s i ) ∈ B i . Based on the way that the workload evolves [see discussion around Eq. (1)] we obtain that W n (t, x) = sup γ ∈ n (x,t) V (γ ). See Fig. 1.
There are three monotonicity properties that the system possesses and which we take into account when analyzing its stability. We start from full vacancy at time −n and consider W n (t, x) for some t ≥ −n. Then W n (t, x) will increase if we (i) delay all arrivals between −n and t, or (ii) increase the heights of the stones, or (iii) enlarge their bases.
Thanks to the monotonicity, it was deduced in [1] that it is enough, for the results sought, to consider the case where the sets B for the team of servers required for a job i with x(i) = 0 is a cube centered at the origin and we write (for a job arriving at server We will work with time doubly infinite, notwithstanding the fact that we consider the process on [−n, ∞).
The first step is the discretization of the Poisson processes. We consider for m ∈ Z and x ∈ Z d the random variables Fig. 1 Graphical representation of the part of the process responsible for the computation of the profile W (t, ·) when the system starts with W (−n, ·) ≡ 0. Consider the evaluation of W (t, 0) at site x = 0. Horizontal intervals represent hailstone (job) arrivals with heights τ i . Only those arrivals which can potentially influence W (t, 0) are shown. Consider a path γ as indicated, from (u, 2) to (t, 0). Its score is We will deal with the discretized system where at server w at time n a job requiring service time T w,n from each server in the cube [w − R w,n , w + R w,n ] d . By the monotonicity (see [1] for detail), this discretization is effective in that if we can show stability for the discretized system of jobs, then we will have shown stability for the original system: The workload at time m for this system will dominate that arising from the nondiscretized model. It is also as well to note that we have not given up too much here. In principle, if we have multiple w jobs arriving during interval (m −1, m], then we could, again in principle, lose if one job required a long service but only from w while a second job required a very short service from a large cube of servers centered at w. However, this will be rare for small λ, where our analysis is most relevant. As noted, we wish to exploit the celebrated results (see [2]) on greedy lattice animal systems. Recall that a lattice animal of Z r is simply a connected subset (when Z r is considered as a graph with the standard edge set). We are given a collection of i.i.d. positive random variables {X (x)} x∈Z r . We suppose the existence of ε > 0 so that EX (0) r +1+ε < ∞ or, equivalently, the existence 1 of ε > 0 and C < ∞ so that for all t > 1, (The 1 in the power is unnecessary but it is in this case that we will use our results.) We will then parametrize our system by taking i.i.d. random variables X λ (x) to be equal to X (x) with probability λ and otherwise 0. For ζ ⊂ Z r , its X λ value (or score) is simply The size |ζ | of the lattice animal is its cardinality. Note that X λ (ζ ) → 0, as λ → 0, in probability, for any lattice animal of finite size. We fix positive integer k and c 1 > 0 and consider the event We wish to prove the following upper bound on the probability of A k , a result which may be of independent interest. We note that P(A k ) depends both on λ and c 1 .

Proposition 5
Given any c 1 > 0, there exists a λ 0 > 0 and a function C : [0, λ 0 ) → (0, ∞) so that C(λ) → 0 as λ → 0 and so that, for λ < λ 0 and for all positive integers k, Remark 2 (i) We can use the above to bound the probability that there is a lattice animal of size u ≥ 2 k containing the origin whose value is ≥ c 1 u, when λ is small, by considering ≥k A c 1 /2 whose probability, by the above, is less than The above formalism will certainly apply to our situation with random variables R x,m + T x,m at each site x ∈ Z r . Indeed, if X (x) denotes the random variable at site x for rate λ = 1 conditioned on there being at least one arrival, then it is easy to see that with rate λ < 1, R x,m + T x,m is stochastically less than X 1−e −λ (x).
Some notation used in the proof and elsewhere. If We also let |x| 1 := r i=1 |x i |. We use 1 A for the indicator of A.
Proof of Proposition 5 We split the value X λ (ζ ), see Eq. (4), into three parts: where The constants q and v appearing in the splitting are chosen as Define next four events: where m is a positive integer satisfying where Note that if a lattice animal ζ of size |ζ | = 2 k contains 0, then ζ ⊂ B k . We obtain an upper bound for P(A k ) via

Bound for P(A k,d ):
Since A k occurs, there is a lattice animal ζ of size 2 k containing the origin and having value X λ (ζ ) ≥ c 1 (5), But this, together with the fact that A k,c does not occur, implies that there is for some constant K m .

Bound for P(A k,c ):
The event A k,c is the event that the sum of at most (2 × 2 k + 1) r Bernoulli random variables, each taking value 1 with probability at most p k = λC2 −kv(r +1+ε) , exceeds m. To bound this probability, we observe that if S n ( p) is the sum of n i.i.d. Bernoulli( p) random variables, then P(S n ( p) ≥ m) is upper bounded by the probability that there is a set A ⊂ {1, . . . , n} of size m such that all Bernoulli random variables are equal to 1 on A, so for some constant K m and thanks to the choice (6) for m.
Lemma 6 (Lemma 1 in [2], Lemma 2.1 in [4]) For any lattice animal ζ of size n containing the origin and any 1 ≤ ≤ n we can find a sequence 0 = u 0 , u 1 , . . . , u h of points in Z r , h being the integer part of 2n/ , and If ζ is a lattice animal containing 0 we can find a sequence π = (π 0 , . . . , π 2n ) such that successive elements are either identical or neighbors in Z r (π is a path) and such that {π 0 , . . . , π 2n } = ζ . (To do this, consider a spanning tree of ζ and form π by traversing the tree "from the bottom.") Then From this, it is immediate that for given there are at most 9 r 2n/ such 2 ball coverings. We use this result with n = 2 k . We consider of "scale" 2 i with 2 i ≤ 2 qk /k 2 . Let i 0 be the maximal such value. For given i, we choose the value = l(i) to equal the integer part of λ −1/2r 2 i/q . With this value the probability that a L ∞ ball of radius 2 contains a site having an X λ value is small for λ small but not (in principle) negligible. From this, it is easily seen that given a sequence u 0 , u 1 , . . . , u h satisfying the above (and therefore given an covering), the probability that the number of sites within the covering having value at least 2 i is at least 2(2 k / )c 1 is bounded above by 20 −2r ·2 k / for λ small. Thus, we see that outside an event of probability 9 hr 20 −2·2 k r/ , this bound will hold for all coverings. Summing over i such that 2 i ≤ 2 kq /k 2 we have that outside probability for each such i and for each corresponding (i) covering, the number of sites in the covering whose X λ value at least 2 i is at most 2(2 k / )c 1 . Thus (outside of probability 2( 1 2 ) c(λk 2/q ) for some universal c), we have, for any lattice animal ζ of size 2 k , which is bounded by Constant(ε)2 k λ (1+ε)/2r . The conclusion follows for large k. Thus, we have shown the proposition.

Corollary 7 Define
For c 1 < 1 fixed, there exists a constant λ 1 = λ 1 (c 1 ) and a function H defined on [0, λ 1 ) tending to zero as λ tends to zero, so that for all 0 < λ < λ 1 and all positive integers R, Proof We treat the case u ≥ R only as that for u ≤ R is essentially the same. Fix By Proposition 5 and Remark 2(i), P(B c 2 u (x)) ≤ C(λ) /2 k(1+ε) where k is the largest integer with 2 k ≤ u. Therefore, , we can create a new lattice animal ξ containing both y and ξ by adding at most (c 1 − c 2 )R/c 2 points to ξ . Since we assumed R ≤ u ≤ |ξ |, we have |ξ | ≤ |ξ | ≤ (c 1 /c 2 )|ξ |. By positivity of the random variables Thus, the event x∈[−R,R] r B c 1 u (x) is a subset of the event that random variable N defined at the start of the proof is at least ( (c 1 −c 2 )R rc 2 ) r . Our result now follows from Markov's inequality.
In this section, we construct clusters for our Poisson hail corresponding to integer intervals (m − 1, m]. The clusters themselves will follow a clustering procedure of [1] and will depend only on the random variables {R x,m } x . Our departure will consist in the temporal (or workload) variable we associate to each cluster. Our clusters will have the property that if C ⊂ Z d is a cluster and γ : (m − 1, m] → Z d is a path satisfying property (ii) of Definition 1, then for some ε > 0 and some finite constant C = C(ε) for any ε conforming to the hypotheses of Theorem 1.
For fixed "time" m and y ∈ Z d let D y,m := B(y, R y,m ) be the L ∞ ball centered at y and having radius R y,m . The cluster C(x, m) containing x is defined as the union of such D y,m over y having the property that there exists integer K and sites y = z 0 , . . . , Let D(x, m) be the diameter of the cluster C(x, m). It is clear that these clusters have the property (8) above. What is not a priori clear is that even with very small rate λ the clusters will be a.s. finite. However, the preceding section enables us to prove Then there exists a function K (λ) tending to zero as λ tends to zero so that, for λ sufficiently small and all positive integers z, Proof Consider the GLA system with random variables {Z (x)} x∈Z d for If the diameter D(0, m) of the cluster C(0, m) containing the origin exceeds z, then there must exist L and 0 = x 0 , x 1 , . . . , x L so that, for all 1 ≤ i ≤ L, and |x L | ≥ z. We choose ζ to be the lattice animal (Recall that |x| 1 denotes the L 1 norm.) Then ζ is a lattice animal in Z d containing the origin for which y∈ζ Z (y) The result follows from Proposition 5 applied to c 1 < 1 4d .

Arguing as in Corollary 7, we obtain
Corollary 9 There is a function C(λ) tending to zero as λ → 0 so that for λ small, for all L, and for R ≤ L/2, while for λ small and R ≥ L/2 We now consider the "time" T (x, m) associated with the cluster C(x, m). This definition is a little less direct than that for D(x, n): Given x ∈ Z d and integer m (and so given cluster C(x, m)), T (x, m) is equal to the maximum value of over sequences x 0 , x 1 , . . . , x L ∈ C(x, m) and m ≥ t 0 ≥ t 1 ≥ · · · ≥ t L ≥ m − 1 so that, for all i a job arrives at x i at time t i = t x,m j (i) having work time τ x i ,m j (i) and We remark that, under the latter two conditions, if x 0 ∈ C(x, m), then necessarily the "subsequent" x i are also in this cluster. We note also that this definition (which requires more information than the discretized data) ensures that, for any site in C(x, m), the waiting time accrued during time interval (m − 1, m] is less than or equal to T (x, m).
Lemma 10 There exists function K (λ) which tends to zero as λ tends to zero so that for λ sufficiently small for all z ≥ 1, Proof By the previous lemma we may suppose that D(0, n) ≤ z/100. Again if T (0, m) takes a value exceeding z, then there must exist L and a sequence x 0 , x 1 , . . . , x L ∈ C(0, m) and times m ≥ t 0 ≥ t 1 ≥ · · · ≥ t L ≥ m − 1 so that, for all i ≤ L − 1, there is a job arrival at x i at time t i and It is important to note that we do not assume that the x i are distinct. Indeed, it is for this reason that we use that bound involving and equally T y,m ≥ k τ . Thus as before we obtain, as in Lemma 8, but with Z (x) = R x,m + T x,m , that for a lattice animal ζ = i P(x i−1 , x i ) that the GLA(Z ) score (i.e., x∈ζ (R x,m + T x,m )) will exceed (z + |ζ |)/4d.

Again we have
Corollary 11 There exists function C(λ) which tends to zero as λ tends to zero so that for so that for all L and for R ≤ L/2, while for λ small and R ≥ L/2

Workload Bounds and Stability
We now apply the foregoing to analyze the workload stability for small values of λ.
It is enough to show tightness of the workload W n (0, 0) at time 0 when the system starts empty at time −n. Recall that W n (0, 0) is obtained as the maximum of scores V (γ ) where γ ranges in the set of paths n (0, 0). See Definition 1.
Due to the monotonicity properties of the system, W n (0, 0) is readily seen to be bounded above by the quantity W n,D (0, 0) which corresponds to the discretized system and is given by where the supremum is taken over discrete time indexed paths γ : [−r, 0] → Z d for some 0 ≤ r ≤ n satisfying The score V D (γ ) of γ is given by We now consider a cube H of length R in Z d+1 = Z d ×Z where the first d coordinates are considered as "spatial" and the last one temporal. Accordingly, we write H as H ×I where I is a temporal interval of length R and H is a cube of length R in Z d . We Proposition 12 There exists constant K λ so that, for u ≤ R, V (H, u) is stochastically less than Poisson of parameter K λ R d+1 (u+1) d+1+ε . For u ≥ R it is bounded by a Poisson of parameter K λ R (u+1) 1+ε . Furthermore, as λ tends to zero, K λ tends to zero. Given the above for C a cluster corresponding to temporal interval (m − 1, m], we define a value X m (C) to signify the value D(x, m) + T (x, m) for a (and so any) x in C.
In analyzing V D (γ ) we will consider a (nonstandard) lattice animal system on Z d+1 . Instead of having i.i.d. random variables indexed by points (x, m) ∈ Z d+1 , we will consider a lattice animal model based on the random variables X m (C) which, while independent for distinct collections of index m, are not independent. Given a lattice animal ⊂ Z d+1 , we write m to denote ∩ Z d × {m} ( of course in general m will not be a lattice animal). Obviously given , all but finitely many m will be empty. The value V C ( ) associated with such a lattice animal will be To analyze sup V C ( ) over all lattice animals containing the origin of Z d+1 and of cardinality N , we proceed as in Sect. 2. For a positive integer k, we let L k := 2 k(1+ε/2(d+1)) .
(We are primarily interested in k with L k ≤ N / log 2 (N )). We know from Sect. 2 that there are less than K . . , C k 2N /L k +1 }, so that each considered is contained in the union of the C k i for one of these collections. Given such a collection we have, by Proposition 12, that for any j, V (C k j , 2 k ) is stochastically less than a Poisson random variable of parameter K λ /2 ε/2 . Furthermore (again by Lemma 2 of [1]), we have that, having identified all clusters intersecting i< j C k i , then conditional number of "extra" clusters intersecting C k j is stochastically less than this Poisson random variable. Thus, just as in Sect. 2, we obtain the following: If N k ( ) is the number of clusters of value more than 2 k that intersect , then, for all N /L k ≤ log 2 (N ) and N large, Thus for every lattice animal of size N containing the origin, outside of probability bounded by Const × e − log 2 (N ) , the contribution to V C (A) from clusters C having for some finite C(ε), where (we stress) K λ tends to zero as λ tends to zero.
Given this bound we easily deal with the clusters having value greater than N 0 using Proposition 12 and the arguments of Sect. 2 and obtain Proposition 13 For the above system, for each δ > 0, there exists a sufficiently small λ > 0, such that the probability that there exists a lattice animal of size at least n, containing the origin of Z d+1 and with V C ( ) > δ| | is less than C δ /n ε/2 for all positive integers n.
Proof We associate to each path γ (−m), γ (−m + 1), . . . , γ (0) the score 0 j=−m+1 Note that, by definition of T (x, n), for each j, the sum of durations for jobs which arrive at time s ∈ ( j − 1, j] and so that the job requires service from both server γ (s) and γ (s−) must be less than T (γ ( j), j) = T (γ ( j − 1), j). Thus, in particular for a continuous time curve γ ∈ m (0, 0), We can associate the path γ with the "lattice animal" γ in Z d+1 consisting of points (γ (i), i) for i = −m, −m + 1, . . . 0 together with for each −m < i ≤ 0 the points (y, i) which lie on a path P i from (γ (i), i) to (γ (i + 1), i) which lies within Z d × {i} and has length |γ (i) − γ (i + 1)| 1 . Thus, this lattice animal ζ has size We note that, while obviously | γ | ≥ m + 1, there are no nonrandom upper bounds for the cardinality. Thus, the inequality above can be rewritten as We now invoke Proposition 13 with δ = 1 20d to deduce that (with λ sufficiently small) the (bad) event B = ∃ lattice animal , 0 ∈ , | | ≥ m, V C ( ) ≥ δ| | has probability bounded by C/m ε/2 for some universal C. Our analysis of now splits into two cases. In both cases, we suppose that bad event B does not occur.
Firstly suppose that | γ | ≤ 4dm + m. In this case 20d . This implies that On the other hand, suppose that | γ | > 4dm + m. Now we have Proof of Theorem 1 From Proposition 14 we have that the workload W n,D (0, 0) of the discretized system is tight as n varies. On the other hand, by monotonicity, the limit as n → ∞ of W n (0, 0) exists a.s. Tightness ensures that this limit is finite. If we now start the system in full vacancy at time 0 and consider the workload W (x, n) for some n > 0, we have that W (x, n) is in distribution equal to W n (0, 0). Therefore, W (x, n) converges in distribution as n → ∞.

Remark 3
We have actually shown something stronger than tightness: Namely, that, starting with an initially empty system, the workload profile at time t converges in distribution, as t → ∞, to some distribution which we will denote by μ. Standard arguments show that μ is an invariant measure: If we start with W (0, ·) distributed according to μ, then W (t, ·) also has distribution μ. Since W (t, ·) is translation invariant in space, this is the case for the limit μ. We have thus proved the existence of an invariant probability measure which is also spatially invariant.

Necessity and Proof of Theorem 2
Let 0 < ε < d + 1. We consider the case where the stone heights (job service times) τ satisfy, for t positive integer, and the stone basis B is the cube We consider the number of job arrivals in space time cube [0, t) d+1 of duration at least 2t for integer t, that is, the number of arrivals (B, τ ) so that (1) τ ≥ 2t, (2) B is a cube of side length 2τ + 1 centered at a site in [0, t) d and (3) the job arrives at a time in [0, t).
This random variable has expectation t ε λ/2 d+1−ε which for fixed λ tends to infinity as t becomes large. In particular for t large this expectation strictly exceeds 1 2 . We fix such a t now.
Obviously this applies to any translation of the cube and the random variables associated to disjoint space time cubes are independent.
We consider the path γ n in nt (0, 0) which is identically the origin for all s ∈ [−tn, 0]. We note that under our assumptions on (B, τ ) any job arriving at a site in [0, t) d having τ ≥ 2t requires service from the origin. Hence the path γ n has value at least where X j is the number of jobs arriving during interval (− jt, −( j −1)t] and satisfying (1) and (2) above. By the law of large numbers, V (γ n ) tends to infinity a.s., as n tends to infinity. This is enough to establish instability of the workload in this case no matter what the value of λ > 0 might be.

Remark 4
In fact, this argument can easily be generalized to show that if for each arriving job, τ = R and if, for some ε > 0, E(τ d+1−ε ) = ∞, then the system does not have stability.

Uniqueness
We now briefly address the question of unicity of invariant measures for the workloads when the power law condition holds and when λ is sufficiently small. We know that if condition (2) is satisfied and parameter λis sufficiently small, then the distribution μ of workloads, obtained by starting the system at time −n with the workloads identically zero and letting n → ∞, is invariant. The question that naturally arises is whether other equilibria for the workload, under Poisson arrival of jobs, are possible.
We consider systems that are stationary under spatial translations and show the following.
Theorem 15 Under the condition (2) above, there exists λ 0 so that if the arrival rate λ is less than λ 0 , and ν is an invariant probability for the system on the space of workloads that is preserved by spatial translation, then ν = μ.
In this section, the assumption that all jobs require service from cubes of servers is not "without loss of generality" so we remark that we only use the weak "irreducibility" condition that for every neighbor e of the origin there exist sequences 0 = x 0 , x 1 , . . . , x r = e 123 and bases B 1 , B 2 , . . . , B r so that for all i, x i−1 , x i ∈ B i and jobs B i occur with strictly positive probability.
To show the claimed uniqueness it suffices to show that for such a measure ν and any bounded cylinder function h, we have Assuming that ν is invariant this is equivalent to for any n (and so, in particular, for n large). Given this and our construction of the measure μ, it will be enough to show that for ε > 0 and h as above, both fixed, for n large. This will be our objective in the following.
As ν is temporally invariant, we have, by the ergodic theorem [5] that, for every M, a.s., where I T is the σ -field of events that are invariant under temporal shifts. Thus, for an ε > 0 fixed, we can find an M so large that, where I S is the sigma field of spatially shift invariant events. Thus, we have that, for k 0 sufficiently large, with probability at least 1 − 2ε, We now note that at time 0, say, the existence of a large workload V at a site 0, say, implies that with reasonable probability the workload will be of order V for a time of order V in the time interval Proof From our "irreducibility" assumptions on the distribution of jobs, it is clear that there exist for each neighbor e of the origin 0 a sequence of jobs with bases B 1 , B 2 , . . . , B R so that, for each i, B i ∩ B i+1 = ∅, 0 ∈ B 1 and e ∈ B R and the rate at which job with base B i arrives is strictly positive. Taking R 1 to be the maximum over the Rs as the neighbor e varies and c to be the minimum over the rates B i as e and i vary, we obtain that, for any x, there exists a "path" B 1 , B 2 , . . . , B R so that R ≤ d R 1 |x| ∞ , for each i, B i ∩ B i+1 = ∅, 0 ∈ B 1 and x ∈ B R and the rate at which job B i arrives is at least c. Thus, for every x ∈ [−c 1 V, c 1 V ] d , the probability that W (x, V /2) ≤ V /2 is bounded by the probability that a parameter V c/2 Poisson process is less than dc 1 V R 1 . The result now follows easily from Poisson tail probabilities.
We note that if V > 4M (assuming as we may that ε < 1/3), then W (x, t) > V /4, for t ∈ [V /2, 3V /4], implies that event B x M does not occur. This implies that Proposition 17 If for some x with |x| ≤ K V we have W (x, 0) ≥ V > 4M, then with probability at least c 1 , for ε fixed small enough.
This yields the simple corollary.
Corollary 18 For M and ε as above, let A(V, K ) be the event that W (x, 0) > t for some t > V and some |x| ≤ K t. Then, under measure ν, the probability that A(V, K ) occurs is less than ε/c 1 provided c 1 /(2K ) d > ε.
From this result, our claim is straightforward.