arXiv:2007.11390v2 [math.PR] 2 Nov 2020 eae tcatcmdl.Teeoe nteecss tis it cases, these in is Therefore, CTMC models. stochastic acc behaved the and conclusions the compromising without ignored be xml,i nepeso o the for expression an if example, [ matrices [ structures [ (BPDs) [ networks type as to referred ttoaydsrbtos us-ttoaydistributio quasi-stationary distributions, stationary the called (usually extinction itiuini rval h ia et esr tzr.I zero. at measure popula delta the sure for Dirac almost example the trivially extinct for is goes cases, distribution eventually other population In the run. ration, long the in system xeti e ae.I h neligsohsi rcs ha process stochastic underlying the If cases. few in except ol sdt oe ope ellrbhvos eeexpres gene behaviors, cellular complex [ model to used monly h ttoaydsrbto ae rdc-om hsi f is This product-form. [ networks a queueing takes distribution stationary the 21 nmn ae ti esnbet xetwl-eae biolog well-behaved expect to reasonable is it cases many In Date l general, in known be not might expression explicit an While 2010 h ttoaydsrbto poie teit)i genera is exists) it (provided distribution stationary The tcatcbooia oesbsdo otnostm Mar continuous-time on based models biological Stochastic e od n phrases. and words Key .Nieaeihrn xrni n nrni rprisof properties intrinsic and extrinsic inherent are Noise ]. us-ttoaydistribution quasi-stationary r oebr2020. November 3rd : ahmtc ujc Classification. Subject Mathematics Abstract. erdcin oeo hc r it-et processes. birth-death po are stochastic which and of processes, none branching reproduction, of s class single-cell extended general biochemical an a to chains), results Markov continuous-time our as apply asymptotics We tail approach. establish mimartingale to approach re The Similar distributions. distributions. thr heavy-tailed (iii) (light-ta by and distributions types butions, Conley-Maxwell-Poisson three (i) into classified ra ers: transition are power-law distributions asymptotic tionary with establi chains is Markov measures stationary time for identity new A chain Markov integers. continuous-time of distributions stationary H SMTTCTISO II ITIUIN OF DISTRIBUTIONS LIMIT OF TAILS ASYMPTOTIC THE 41 45 ergodic 24 ,a ela o eeaiain fsc rcse ihmore with processes such of generalizations for as well as ], , 49 , ii distributions limit 25 HAGX,MD HITA ASN N ASE WIUF CARSTEN AND HANSEN, CHRISTIAN MADS XU, CHUANG 30 .QD ihepii xrsin peri vnrrrcases rarer even in appear expressions explicit with QSDs ]. hti,teeeit a exists there is, that , ,cmlxblne ecinntok [ networks reaction complex-balanced ], ,i admevrnet [ environments random in ], hspprivsiae alaypoiso ttoaydis stationary of asymptotics tail investigates paper This 54 OTNOSTM AKVCHAINS MARKOV TIME CONTINUOUS , 32 otnostm akvcan,iett o ttoaymea stationary for identity chains, Markov time Continuous ,Jcsnntok n akt-hnyMnzPlco (B Baskett-Chandy-Muntz-Palacios and networks Jackson ], Q -process ntepeetpaper. present the in QD,ta s h ogtm eairo h rcs before process the of behavior long-time the is, that (QSD), tail 1. itiuini nw,te h xsec n relative and existence the then known, is distribution rmr 02;scnay6J8 07,90E20. 60J74, 60J28, secondary 60J27; Primary [ ) Introduction unique 41 .Jity ttoaydsrbtosadQD are QSDs and distributions stationary Jointly, ]. s tcatcrato network. reaction stochastic ns, 1 20 ,o ihtiignlboktasto rate transition block tridiagonal with or ], ttoaydsrbto hc ecie the describes which distribution stationary ut r eie o quasi-stationary for derived are sults ld,(i xoeta-alddistri- exponential-tailed (ii) iled), hd npriua,frcontinuous- for particular, In shed. hs ae,i ae es ostudy to sense makes it cases, these n nasbe ftenon-negative the of subset a on s 3 ohsi eeepeso model, expression gene tochastic sdffrn rmtecasclse- classical the from different is a s l icl osaei xlctform, explicit in state to difficult lly , y n hsteegdcstationary ergodic the thus and ly, eesl optbeparamet- computable easily ee uhbooia ytm n cannot and systems biological such uainpoesswt bursty with processes pulation reapetecs o reversible for case the example or 38 aua oasm h modelling the assume to natural e,ti smttc o sta- for asymptotics tail tes, ealdbalanced detailed ecinntok (modeled networks reaction , inadteeouino DNA of evolution the and sion rc ftemodels. the of uracy o his(TC)aecom- are (CTMCs) chains kov inpoesswtotimmig- without processes tion 11 s uc nmn ae.For cases. many in suffice ess clsses hsas well- also thus systems, ical ,adbrhdahprocesses birth-death and ], rbtosadquasi- and tributions ue,ti distribution, tail sures, eea graphical general [ 53 tutr,then structure, ]. CMP)- 2 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF sizes of moments might be assessed from the decay rate of the tail distribution. Additionally, relative recurrence times might be assessed for stationary distributions of CTMCs. With this in mind, we aim to establish results for the tail behavior of a stationary distri- bution or a QSD, provided either one exists. In particular, we concentrate on CTMCs on the non-negative integers N0 with asymptotic power-law transition rate functions (in a sense to be made specific) [2, 7]. Our approach is based on a simple observation, namely a new and generic identity (and thus an equivalent definition) for limit distributions as well as stationary measures (Theorems 3.2 and 3.4). The identity is a consequence of the master equation. For BDPs it coincides with the product-form expression for the stationary distribution, namely, π(n)= cnπ(0), where cn is a constant depending on the birth and death rates. Furthermore, the identity allows us to study the tail behavior of limit distributions, provided they exist, and to characterize their forms (Theorems 4.1 and 4.4). Specifically, in Section 5, for CTMCs with transition rate functions that are asymptotically power law, we show that there are three regimes: Either the decay follows (i) a Conley-Maxwell-Poisson distribution (light-tailed), (ii) a geometric distribution, or (iii) a heavy-tailed distribution. We apply our main results to biochemical reaction networks, a general single-cell stochastic gene expression model, an extended class of branching processes, and stochastic population processes with bursty reproduction, none of which are BDPs. Tail asymptotics of stationary distributions have been intensively investigated in queueing theory [9, 5, 33], as well as for discrete time Markov chains (DTMCs) [42, 6, 19]. For this reason, we provide a brief comparison between our results and existing results. The Lyapunov function approach (also called the martingale approach) is widely taken as a standard method to prove stability (ergodicity) of Markov processes [43]. Additionally, the approach has been used to obtain the tail asymptotics of stationary distributions (assuming ergodicity) [42, 6, 10]. For DTMCs with bounded jumps and zero asymptotic drift, heavy tail asymptotics has been obtained [42, 6, 19]. The motivation arises from Lamperti’s problem [34] when the drift is of order x−1 for large x. By constructing suitable Lyapunov functions [42, 19, 56] (see also [6]), asymptotic tail estimates have been established when drift is critical (order x−1) and when slowly decaying (order x−r, 0 2. Preliminaries 2.1. Sets and functions. Denote the set of real numbers, positive real numbers, integers, positive integers, and non-negative integers by R, R+ Z, N, and N0, respectively. For m,n N, let Rm×n denote the set of m by n matrices over R. Further, for any set B, let #B∈ denote its cardinality and ½B the corresponding indicator function. For b R, A R, let bA = ba: a A , and A + b = a + b: a A . Given A R, let min∈ A and⊆ max A denote{ the minimum∈ } and maximum{ of the set,∈ respectively.} By⊆ convention, max A = , min A = + if A = ∅, and max A = + if A is unbounded from above, and min A =−∞ if A is unbounded∞ from below. ∞ −∞ Let f and g be non-negative functions on an unbounded set A R+. We denote f(x) . g(x) if there exists C,N > 0 such that ⊆ f(x) Cg(x), for all x A, x N, ≤ ∈ ≥ that is, f(x)=O(g(x)) since f is non-negative. Here O refers to the standard big O notation. The function f is said to be asymptotic power-law (APL) if there exists r1 R such that f(x) log f(x) ∈ limx→∞ xr1 = a exists and is finite. Hence r1 = limx→∞ log x . An APL function f is called hierarchical (APLH) on A with (r1, r2, r3) if there further exists r2, r3 with r2 +1 r > r > r r 2, and a> 0, b R, such that for all large x A, ≥ 1 2 3 ≥ 1 − ∈ ∈ f(x)= axr1 + bxr2 + O(xr3 ). The requirement r2 +1 r1 and r3 r1 2 comes from the analysis in Sections 6-7, where asymptotic Taylor expansion≥ of functions≥ − involves the powers of first few leading terms. Here r1, r2, and r3 are called the first, second and third power of f, respectively. All rational 4 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF functions, polynomials, and real analytic APL functions are APLH. Not all APL functions are APLH, e.g., f(x) = (1 + (log(x + 1))−1)x on N. For an APLH function f, its first power is uniquely determined, while the other two powers ∗ ∗ 2 are not. Let r2 and r3 be the infimums over all (r2, r3) R such that f is APLH on A with ∈ ∗ ∗ (r1, r2, r3). For convention, we always choose the minimal powers (r1, r2 , r3 ) whenever f is an ∗ ∗ 2 APLH with (r1, r2, r3 ). As an example, f(x)= x +3x +4 is APLH on N0 with (2, r2, r3) for any 2 > r2 > r3 1 (in which case b =0)or1= r2 > r3 0 (b = 3). In this case, f is APLH ≥ ∗ ∗ ≥ 1/3 on N0 with minimal powers (r1, r2 , r3 ) = (2, 1, 0). In contrast, take f(x) = x + x log x. Then f is APLH on N with (1, r2, r3) for any r2 > r3 > 1/3 (b = 0). In this case, r1 = 1 ∗ ∗ and r2 = r3 =1/3, but f is not APLH on N with (1, 1/3, 1/3). For any real analytic APLH function f on N , f is APLH on N with (r , r∗, r∗), where r = lim log f(x) , r∗ = r 1 0 0 1 2 3 1 x→∞ log x 2 1 − and r∗ = r 2. 3 1 − 2.2. Measures. Any positive measure µ on a set A N0 can be extended naturally to a positive measure on N with no mass outside A, µ(N ⊆A) = 0. For any positive measure µ 0 0 \ on N0, let ∞ T : N [0, 1], x µ(y), µ 0 → 7→ y=x X be the tail distribution (or simply the tail) of µ. Let be the set of probability distributions on A. For a,b > 0, define P 1+ = µ : T (x) . exp( ax log x(1 + o(1))) , Pa { ∈P µ − } 1− = µ : T (x) & exp( ax log x(1 + o(1))) , Pa { ∈P µ − } 2+ = µ : T (x) . exp( bxa(1 + o(1))) , Pa,b { ∈P µ − } 2− = µ : T (x) & exp( bxa(1 + o(1))) , Pa,b { ∈P µ − } 3+ = µ : T (x) . x−a , Pa { ∈P µ } 3− = µ : T (x) & x−a , Pa { ∈P µ } where o refers to the standard little o notation. Furthermore, define 2+ = 2+, 2− = 2−, Pa ∪b>0Pa,b Pa ∪b>0Pa,b 2+ = 2+, 2− = 2−, P>1 ∪a>1Pa P<1 ∪00Pa P ∪a>0Pa i+ i− The sets a , i = 1, 2, 3, are decreasing in a, while a , i = 1, 2, 3, are increasing in a. P 2+ 2−P Similarly, a,b is decreasing in both a and b, while a,b is increasing in both a and b. P 2+ 2− P The probability distributions in 1 1 decay as fast as exponential distributions and P ∩P 1+ 2+ are therefore exponential-tailed. Similarly, those in >1 are heavy-tailed probability 3− 2− P ∪P distributions, and those in <1 are light-tailed [31]. P ∪P 2 The Conley-Maxwell-Poisson (CMP) distribution on N0 with parameter (a,b) R+ has probability mass function given by [31]: ∈ −1 ∞ ax aj CMP (x)= , x N . (a,b) (x!)b (j!)b ∈ 0 j=0 X 1+ In particular, CMPa,1 is a Poisson distribution. For every probability distribution µ 1−, there exists (a ,b ), (a ,b ) R2 , such that ∈P ∩ P 1 1 2 2 ∈ + CMP CMP (a1,b1) . Tµ(x) . (a2,b2). THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 5 The Zeta distribution on N0 with parameter a> 1 has probability mass function given by [31]: 1 Zeta (x)= x−a, a ζ(s) ∞ −a where ζ(a)= i=1 i is the Riemann zeta function of a. For every probability distribution 3+ 3− µ , there exists a1, a2 > 1 such that ∈P ∩P P Zeta Zeta a1 (x) . Tµ(x) . a2 (x). 2.3. Markov chains. Let (Yt : t 0) (or Yt for short) be a minimal CTMC with state space N and transition rate matrix≥ Q = (q ) , in particular each entry is finite. Recall Y ⊆ 0 x,y x,y∈Y that a set A is closed if qx,y = 0 for all x A and y A [46]. Let ∂ ( be a finite closed absorbing⊆ Y set, and define ∂c = ∂. Furthermore,∈ ∈Y\ define Y Y\ Ω= y x: q > 0, for some x, y , { − x,y ∈ Y} and the transition rate functions by, λ (x)= q , x , ω Ω. ω x,x+ω ∈ Y ∈ Let Ω± = ω Ω: sgn(ω)= 1 be the sets of forward and backward jump vectors, respect- ively. { ∈ ± } For any probability distribution µ on ∂c, define P ( )= P ( )dµ(x), µ · x · ZY where P denotes the probability measure of Y with initial condition Y = x . x t 0 ∈ Y A (probability) measure π on is a stationary measure (distribution) of Yt if it is a non- negative equilibrium of the so-calledY master equation [26]: (2.1) 0= λ (x ω)π(x ω) λ (x)π(x), x . ω − − − ω ∈ Y ωX∈Ω ωX∈Ω (Here and elsewhere, functions defined on are put to zero when evaluated at x Z.) Any stationary distribution π of Y satisfiesY for all t 0, 6∈ Y ⊆ t ≥ P (Y A)= π(A), A 2Y , π t ∈ ∈ where 2Y is the power set of [14]. Y Let τ∂ = inf t > 0: Yt ∂ be the entrance time of Yt into the absorbing set ∂. We say Y admits certain{ absorption∈ }if τ < almost surely (a.s.) for all Y ∂c. Moreover, the t ∂ ∞ 0 ∈ process associated with Yt conditioned to be never absorbed is called a Q-process [13]. A c probability measure ν on ∂ is a quasi-stationary distribution (QSD) of Yt if for all t 0, c ≥ P (Y A τ >t)= ν(A), A 2∂ . ν t ∈ | ∂ ∈ 3. Identities for limit distributions Let ω∗ = gcd(Ω) be the (unique) positive greatest common divisor of Ω. Define the scaled largest positive and negative jump, respectively: −1 −1 ω+ := max ωω∗ , ω− := min ωω∗ . ω∈Ω+ ω∈Ω− Furthermore, define for j ω ,...,ω +1 , ∈{ − + } ω Ω− : jω∗ >ω , if j ω−,..., 0 , (3.1) Aj = { ∈ } ∈{ } ω Ω : jω ω , if j 1,...,ω +1 . ({ ∈ + ∗ ≤ } ∈{ + } Hence, ∅ = A A A A =Ω for ω Proposition 3.1 ([17]). Assume ∂ = ∅. Let ν be a QSD of Y on ∂c. Then for x N ∂, 6 t ∈ 0 \ θ ν(x)+ λ (x ω)ν(x ω) λ (x)ν(x)=0, ν ω − − − ω ωX∈Ω ωX∈Ω where θν = ν(y)λω(y) c ωX∈Ω− y∈∂ X∩(∂−ω) is finite. c Proof. The identity for x ∂ is stated in [17]. For x N0 , the identity trivially holds (with both sides being zero),∈ which can be argued similarly∈ to\ Ythe proof of Theorem 3.2. The following generic identities provide an equivalent definition of stationary distributions but in a more handy form that turns to be useful in estimating tails of stationary distributions. Theorem 3.2. The following statements are equivalent. (1) π is a stationary measure (stationary distribution) of Yt on . (2) π is a positive measure (probability distribution) on satisfying,Y for x N , Y ∈ 0 −1 0 ωω∗ λω(x jω∗)π(x jω∗)= λω(x jω∗)π(x jω∗) < . − − − − − ∞ ω∈Ω− 1 ω∈Ω j=1 X j=ωωX∗ +1 X+ X (3) π is a positive measure (probability distribution) on satisfying, for x N , Y ∈ 0 0 ω+ π (x jω ) λ (x jω )= π (x jω ) λ (x jω ) < . − ∗ ω − ∗ − ∗ ω − ∗ ∞ j=ω−+1 j=1 X ωX∈Aj X ωX∈Aj Proof. We only prove for the equivalent representations for stationary measures. Then the equivalent representations for stationary distributions follow. We use LHS (RHS) as shorthand for left (right) hand side of an equation. First, we claim that π is a stationary measure if and only if (2.1) holds for all x N0. It suffices to show (2.1) holds for x N , provided π is stationary measure for Y on ∈. Since ∈ 0 \Y t Y x N0 and π(N0 ) = 0, the LHS of (2.1) equals zero. Since is closed, x ω implies∈ \λ Y(x ω) = 0\ (otherwise, Y x = x ω + ω ). Hence π(x Y ω)λ (x ω)− = 0∈ for Y ω − − ∈ Y − ω − all ω Ω, which means the RHS of (2.1) is also zero. This shows (2.1) holds for x N0 , provided∈ π is a stationary measure on . ∈ \ Y Moreover, we can rewrite (2.1) as follows:Y (3.2) λ (x ω)π(x ω) λ (x)π(x) = λ (x)π(x) λ (x ω)π(x ω) . ω − − − ω ω − ω − − ω∈Ω− ω∈Ω X X+ Indeed, since Q is a transition rate matrix, row sums of off diagonal entries are finite. Hence also λ (x)π(x) < , x N , ω ∞ ∈ 0 ωX∈Ω with the adopted convention that functions are zero for x Z . From (2.1) it follows that ∈ \ Y λ (x ω)π(x ω) < , x N . ω − − ∞ ∈ 0 ωX∈Ω This allows us to divide the sum of non-negative real numbers into two finite sums: λ (x ω)π(x ω)= λ (x ω)π(x ω)+ λ (x ω)π(x ω) < , ω − − ω − − ω − − ∞ ωX∈Ω ωX∈Ω− ωX∈Ω+ λ (x)π(x)= λ (x)π(x)+ λ (x)π(x) < . ω ω ω ∞ ωX∈Ω ωX∈Ω− ωX∈Ω+ THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 7 Hence (3.2) is valid since all sums are finite. Assume w.o.l.g. that ω∗ = 1 and 0 N0. Otherwise, make the transformation −1 ∈Y ⊆ Xt = (Yt min )ω∗ . Then the process (Xt : t 0) is a CTMC on , and 0 N0 with − Y ≥ X∗ −1 ∈X ⊆ ω∗. (The same construction is applied elsewhere.) Moreover, Ω = (ω ) Ω is the set of jump vectors for the process Xt with gcd(Ω) = 1. (1) (2). We prove it by induction. e ⇒ Base case. Let x = 0. The LHSe of the equation in (2) is zero, since λω(j) = 0, for 0 j < ω, ω Ω . Similarly, the RHS of the equation in (2) also vanishes, since π( j)=0 ≤ − ∈ − − for 1 j ω+. Induction≤ ≤ step. Assume 0 ω λ (x j)π(x j)= λ (x j)π(x j) < ω − − ω − − ∞ − j=ω+1 j=1 ωX∈Ω X ωX∈Ω+ X holds for some x N . Since π is a stationary measure, it follows from (2.1) that ∈ 0 (λ (x ω)π(x ω) λ (x)π(x)) = (λ (x)π(x) λ (x ω)π(x ω)) < . ω − − − ω ω − ω − − ∞ ωX∈Ω− ωX∈Ω+ Hence 0 λ ((x + 1) j)π((x + 1) j) ω − − − j=ω+1 ωX∈Ω X 0 = λ (x (j 1))π(x (j 1)) ω − − − − − j=ω+1 ωX∈Ω X 0 = λ (x j)π(x j)+ (λ (x ω)π(x ω) λ (x)π(x)) ω − − ω − − − ω − j=ω+1 − ωX∈Ω X ωX∈Ω ω = λ (x j)π(x j)+ (λ (x)π(x) λ (x ω)π(x ω)) ω − − ω − ω − − j=1 ωX∈Ω+ X ωX∈Ω+ ω−1 = λ (x j)π(x j) ω − − j=0 ωX∈Ω+ X ω = λ (x +1 j)π(x +1 j) < , ω − − ∞ j=1 ωX∈Ω+ X and the proof is completed. (1) (2). It is a direct result of Fubini’s theorem. For x N : ⇔ ∈ 0 0 0 π(x j) λ (x j)= λ (x j)π(x j), − ω − ω − − j=ω−+1 − j=ω+1 X ωX∈Aj ωX∈Ω X ω+ ω π(x j) λ (x j)= λ (x j)π(x j). − ω − ω − − j=1 j=1 X ωX∈Aj ωX∈Ω+ X (3) (1). Since both sides of the equation are non-negative and finite for all x N0, then by⇒ induction, subtracting the equation in (3) for x from the equation for x + 1 yields∈ the desired conclusion. Corollary 3.3. Assume Ω+ = ∅ or Ω− = ∅. Let π be a stationary measure of Yt on , then supp π ( supp λ )= ∅. Y ∩ ∪ω∈Ω ω 8 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF Proof. Assume w.o.l.g. that Ω+ = ∅ and ω∗ = 1. (The other case Ω− = ∅ is similar.) From Theorem 3.2(3), the RHS is always zero, and hence π (x jω ) λ (x jω )=0, x N . − ∗ ω − ∗ ∈ 0 ωX∈Aj Let j = 0. Then π(x) λ (x)=0, x N , ω ∈ 0 ωX∈Ω− which implies that supp π x N0 : λω(x)=0 , ⊆ω ∩∈Ω{ ∈ } and the conclusion holds. This corollary shows that if a CTMC jumps uni-directionally (e.g., a pure birth or a pure death process), then all stationary measures, if such exists, are concentrated on absorbing states [55]. A special form of Theorem 3.2 under more assumptions has been stated in the context of stochastic reaction networks [28, Prop.5.4.9]. The following identities for QSDs also present equivalent definitions of the latter. Theorem 3.4. Assume ∂ = ∅. Then the following statements are equivalent. 6 c (1) ν is a QSD of Yt on ∂ . (2) ν is a probability measure on ∂c, and for x N ∂, ∈ 0 \ −1 0 ωω∗ λ (x jω )ν(x jω )= θ T (x)+ λ (x jω )ν(x jω ) < . ω − ∗ − ∗ ν ν ω − ∗ − ∗ ∞ ω∈Ω− −1 ω∈Ω j=1 X j=ωωX∗ +1 X+ X (3) ν is a probability measure on ∂c, and for x N ∂, ∈ 0 \ 0 ω+ ν (x jω ) λ (x jω )= θ T (x)+ ν (x jω ) λ (x jω ) < . − ∗ ω − ∗ ν ν − ∗ ω − ∗ ∞ j=ω−+1 j=1 X ωX∈Aj X ωX∈Aj Proof. The proof is similar to that of Theorem 3.2 and thus omitted. c Corollary 3.5. Let ν be a QSD of Yt on ∂ . If Ω− = ∅, then supp ν ( ω∈Ω supp λω)= ∅. In particular, if ∂c = supp λ then there does not exist a QSD of∩ Y∪ on ∂c. ∪ω∈Ω ω t Proof. Ω− = ∅ implies θν = 0. The rest of the proof is similar to that of Corollary 3.3. The difference between Corollaries 3.3 and 3.5 lies in the fact that a QSD ν may exist with positive probability on ω∈Ω supp λω, provided Ω+ = ∅, while Ω− = ∅, as illustrated by the following example. ∪ 6 Example 3.6. Consider a pure death process Yt on N0 with linear death rates dj = dj for j N . Let 0 1, d Γ(1 a/d) − where Γ( ) is the Gamma function. Then ν is a QSD of Y on N (∂ = 0 ), and supp ν = N. · t { } Formulae for stationary distributions and QSDs of BDPs follow directly from Theorems 3.2 and 3.4. Corollary 3.7 ([4, 12]). (i) Let Yt be a BDP on N0 with birth and death rates bj and dj , respectively, such that b > 0 and d > 0 for all j N. If π is a stationary distribution for j−1 j ∈ Yt, then j−1 b π(j)= π(0) i , j N. d +1 ∈ i=0 i Y THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 9 (ii) Let Yt be a BDP on N0 with birth and death rates bj and dj , respectively, such that b = 0, and b > 0 and d > 0 for all j N. Then a probability distribution ν on N is a QSD 0 j j ∈ trapped into 0 for Yt if and only if j−1 (3.3) d ν(j)= b ν(j 1)+ d ν(1) 1 ν(i) , j 2. j j−1 − 1 − ≥ i=1 ! X Proof. Here Ω = 1, 1 , ω∗ = ω+ = 1, ω− = 1, and = N0. Moreover, λ−1(j) = dj and λ (j)= b for j {−N. } − Y 1 j ∈ (i) ∂ = ∅. Since π is a stationary distribution on N0, it follows from Theorem 3.2 that π(j)λ (j)= π(j 1)λ (j 1), j N. −1 − 1 − ∈ Hence the conclusion is obtained by induction. (ii) ∂ = 0 and ∂c = N. It follows from Theorem 3.4 that a probability measure ν is a QSD on N if{ and} only if θ = λ (1)ν(1), ν(j)λ (j)= θ T (j)+ ν(j 1)λ (j 1), j N 1 , ν −1 −1 ν ν − 1 − ∈ \{ } that is, (3.3) holds. Regarding the tail distributions, we have the following identities. Corollary 3.8. Assume Ω is finite and ∂ = ∅. Let π be a stationary distribution of Yt on . Then, for x N , Y ∈ 0 −1 T (x) λ (x)+ λ (x ω ) + T (x jω ) π ω ω − ∗ π − ∗ j=ω− ωX∈A0 ωX∈A1 X λ (x jω ) λ (x (j + 1))ω ) · ω − ∗ − ω − ∗ ωX∈Aj ω∈XAj+1 ω+ = T (x jω ) λ (x jω ) λ (x (j + 1)ω ) π − ∗ ω − ∗ − ω − ∗ j=1 X ωX∈Aj ω∈XAj+1 where Aj is defined in (3.1). Proof. Assume w.l.o.g. ω =1 and 0 . The LHS of the equation in Theorem 3.2(3) is ∗ ∈ Y 0 LHS = (T (x j) T (x j + 1)) λ (x j) π − − π − ω − j=ω−+1 X ωX∈Aj 0 −1 = T (x j) λ (x j) T (x j) λ (x j 1) π − ω − − π − ω − − j=ω−+1 j=ω− X ωX∈Aj X ω∈XAj+1 −1 = T (x j) λ (x j) λ (x j 1) + T (x) λ (x), π − ω − − ω − − π ω j=ω− X ωX∈Aj ω∈XAj+1 ωX∈A0 while the RHS equals ω+ RHS = (T (x j) T (x j + 1)) λ (x j) π − − π − ω − j=1 X ωX∈Aj ω+ = T (x j) λ (x j) λ (x j 1) T (x) λ (x 1), π − ω − − ω − − − π ω − j=1 X ωX∈Aj ω∈XAj+1 ωX∈A1 which together yield the desired identity. 10 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF c Corollary 3.9. Assume Ω is finite, ∂ = ∅, and let ν be a QSD of Yt on ∂ . Then for all x N ∂, 6 ∈ 0 \ −1 T (x) λ (x)+ λ (x ω ) + T (x jω ) ν ω ω − ∗ ν − ∗ j=ω− ωX∈A0 ωX∈A1 X λ (x jω ) λ (x (j + 1))ω ) · ω − ∗ − ω − ∗ ωX∈Aj ω∈XAj+1 ω+ = θ T (x)+ T (x jω ) λ (x jω ) λ (x (j + 1)ω ) , ν ν ν − ∗ ω − ∗ − ω − ∗ j=1 X ωX∈Aj ω∈XAj+1 where Aj is defined in (3.1). Proof. Similar to that of Corollary 3.8. 4. Asymptotic tails of limit distributions To establish the asymptotic tails of limit distributions, we assume the following. (A1) #Ω < . ∞ 1 2 3 (A2) is unbounded, and for ω Ω, λω is an APLH function on with (Rω, Rω, Rω) and whichY strictly positive for all large∈x . Y ∈ Y (A3) ∂c is irreducible. Assumption (A1) guarantees that the chain have bounded jumps. Assumption (A2) is common in applications. Moreover, it implies that both ∂c and are unbounded (see Pro- position A.1). In particular, (A2) is satisfied provided the followingY assumption holds: (A2)’ For ω Ω, λ is a strictly positive polynomial for all large x . ∈ ω ∈ Y Assumption (A3) is assumed for the ease of exposition and to avoid non-essential technic- alities. Moreover, (A3) means that either Yt is irreducible or the conditional process of Yt before entering ∂ is irreducible. This assumption is satisfied for many known one-dimensional infinite CTMCs modeling biological processes (e.g., for population processes). In addition, (A3) implies that Ω+ = ∅ and Ω− = ∅ (otherwise there are no non-singleton communicating classes). 6 6 The following parameters are well-defined and finite. Let 1 2 3 1 2 3 E− = Rω, Rω, Rω , E+ = Rω, Rω, Rω , ω∈∪Ω−{ } ω∈∪Ω+{ } and define R = max E E , R = max E , R = max E , − ∪ + − − + + σ = min R R1 , R R1 , σ = min R R2 , R R2 , 1 { − − − + − +} 2 { − − − + − +} where R1 = max r E : r < R , R1 = max r E : r < R , − { ∈ − −} + { ∈ + +} R2 = max r E : r < R1 , R2 = max r E : r < R1 . − { ∈ − −} + { ∈ + +} (These values are only used to define σ1, σ2.) 1 2 1 Hence R− Rω R− 1, for some ω Ω− with Rω = R−. Similarly, R+ R+ 1, 2 ≥ ≥2 − ∈ ≥ − R− R− 2, and R+ R+ 2. This implies that 0 < σ1 1 and σ1 < σ2 2. If all ≥ − ≥ − ≤2 1 3 ≤ 1 transition rate functions are real analytic, then by convention, Rω = Rω 1, Rω = Rω 2, for all ω Ω, and hence σ = 1 and σ = 2. Furthermore, let − − ∈ 1 2 λ (x) ω λ (x)ω ω∈Ω λω(x)ω ω∈Ω− ω ω∈Ω+ ω α = lim , α− = lim | | , α+ = lim . x→∞ xR x→∞ xR− x→∞ xR+ P P P THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 11 R 2 − λω(x) λω(x)ω αx 1 λω(x)ω β = lim ω∈Ω , γ = lim ω∈Ω − , ϑ = lim ω∈Ω , x→∞ xR− x→∞ xR−σ1 2 x→∞ xR P P P −1 γ(α+ω∗) if σ1 < 1, −1 ∆= − −1 , δ = ∆(ω+ ω− 1) . ( γ + Rϑ)(α+ω∗) if σ1 =1. − − − We emphasize that α, α−, α+,β and ϑ do not depend on the choice of second and third powers of the transition rate functions, whereas σ1, σ2, γ, ∆ and δ do depend on the powers. In the next three statements we further characterize the tail of the stationary distributions. Theorem 4.1. Assume (A1)-(A2) and ∂ = ∅. Let π be a stationary distribution of Yt on with unbounded support. Then α 0, and in particular when α =0, Y ≤ – if ∆=0 then σ1 < 1, – if σ1 =1 then ∆ > 1. Moreover, if 1+ 1− (i) R > R , then π − − , − + (R−−R )(ω ω∗) 1 1 ∈P + + ∩P(R−−R+)ω∗ 2+ 2− (ii) R− = R+ and α< 0, then π 1 1 . ∈P ∩P2+ 2− (iii) α =0, ∆ > 0, and σ1 < 1, then π 1−σ1 1−σ1 , 3+∈P ∩P (iv) α = 0 and σ1 = 1, then π ∆−1. In particular, if in addition, (iv)’ δ > 1, then 3+ 3− ∈ P π ∆−1 δ−1. ∈P ∩P 3− (v) α =0, ∆=0, and σ2 < 1, then π . ∈P1−σ2 (vi) α =0, ∆=0, and σ 1, then π 3−. 2 ≥ ∈P As a consequence, a trichotomy regarding the tails of the stationary distributions are derived. Corollary 4.2. Assume (A1)-(A2) and ∂ = ∅. Any unboundedly supported stationary distribution of Y on is t Y – light-tailed and its tail decays like a CMP distribution if Theorem 4.1(i) holds, – exponential-tailed if Theorem 4.1(ii) holds, – heavy-tailed if one of Theorem 4.1(iii), (iv)’, and (vi) holds. In particular the tail decays like a Zeta distribution if (iv)’ holds. ′ Corollary 4.3. Assume (A1), (A2) , (A3), ∂ = ∅, R 3, and (R 1)ϑ α+ 0. Any unboundedly supported stationary distribution of Y on ≥is ergodic. − − ≤ t Y ′ Proof. By (A2) , R N0 and σ1 =1. By[56, Theorem 3.1], Yt is explosive if either (1) R 2 and α > 0, or (2) R∈ 3, α = 0 and γ ϑ > 0. By Theorem 4.1, α 0. If a stationary≥ ≥ − ≤ distribution exists and if Yt is non-explosive, then Yt is is positive recurrent and the stationary distribution is unique and ergodic [46]. When α = 0 and R 3, it follows from Theorem 4.1 ≥ that ∆ > 1, i.e., γ ϑ < (R 1)ϑ α+ 0 as assumed. Hence Yt is always non-explosive and thus the stationary− distribution− − is ergodic.≤ c Theorem 4.4. Assume (A1)-(A3) and ∂ = ∅. Let ν be a QSD of Yt on ∂ . Then α 0 R, and if 6 ≤ ≤ – R = 0, then α− θν . If, in addition, R− = R+, then α θν , and if R− > R+, then β θ . ≥ ≤ − ≥ ν – R− = R+ > 0 and α =0, then R > σ1. Moreover, if R− > R+, and 1− (i) R =0 and β>θν , then ν −1 , ∈P(R−−R+)ω∗ 2− (ii) R =0, β = θν and R− R+ 1, then ν 1 , − ≤ ∈P2− (iii) R =0, β = θν and R− R+ 1, then ν 1 , 1− − ≤ ∈P 1+ (iv) R> 0, then ν − . If, in addition R> 1, then ν − . 1 (R−−R )(ω ω∗) 1 ∈P(R−−R+)ω∗ ∈P + + 12 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF If R+ = R−, and 2− 2+ (v) R> 0 and α< 0, then ν 1 . If, in addition, R> 1, then ν 1 , (vi) R> 0, α =0, and ∈P ∈P 2− – R = σ1 < 1, then ν 1−R, ∈P3− – R σ1 =1, then ν , ≥ ∈P 2− – min 1, R > σ1, then ν 1< , (vii) R =0, and{ } ∈P 2− – α + θν =0 and σ1 < 1, then ν 1−σ1 , ∈P3− – α + θν =0 and σ1 =1, then ν , – α + θ < 0, then ν 2−. ∈P ν ∈P1 Furthermore, if 3+ (viii) R =1, then ν − , θ α 1 ∈P ν + 2+ (ix) 0 λ (x)= ½N(x), λ (x)= ab , j N, x N , −1 j j−1 ∈ ∈ 0 where a> 0, and bj 0 for all j N0. Here the backward jump 1 represents the degradation of mRNA with unity≥ degradation∈ rate, and the forward jumps− j N account for bursty ∈ THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 13 production of mRNA with transcription rate a and burst size distribution (bj )j∈N0 . When j bj = (1 δ)δ for 0 <δ< 1, the stationary distribution is the negative binomial distribution [50]: − Γ(x + a) π(x)= δx(1 δ)a, x N . Γ(x + 1)Γ(a) − ∈ 0 When a = 1, π is also geometric. While Theorem 4.1, if it did apply, would seem to suggest Tπ decays like a CMP distribution, since 1 = R− > R+ = 0. The technical reason behinds this, is that the proof applies Corollary 3.8 which requires (A1). It does not seem possible to directly extend the result in Theorem 4.1 to CTMCs with unbounded jumps. Example 4.7. Consider a BDP with birth and death rates: b λ (x)= S(b, j)xj , λ (x)= a, x N , −1 1 ∈ 0 j=1 X where a> 0, b N, and S(i, j) is the Stirling numbers of the second kind [1]. Here R+ =0 < R = b, and α ∈= S(b,b)= 1 < 0. This BDP has an ergodic stationary distribution on N − − − 0 [56], and the unique stationary distribution is π = CMPa,b. By Theorem 4.4(v), the tail of a QSD decays no faster than exponential distributions when α< 0 and 0 R = R 1, which is also confirmed by the examples below. ≤ − + ≤ Example 4.8. Consider the linear BDP on N0 with birth and death rates: λ (x)= bx, x N , and λ (1) = d, λ (x)= d 2−1 + b (x + 1), x N 1 , 1 ∈ 0 −1 −1 · ∈ \{ } where b and d are positive constants [47]. For this process, a QSD νis 1 ν(x)= , x N. x(x + 1) ∈ Hence T decays as fast as the Zeta distribution with parameter 2. Here α = d/2 < 0, and ν − R = R− = R+ = 1. Example 4.9. Consider the linear BDP on N0 with bj = bj and dj = d1j with 0 Example 4.10. Consider a BDP on N0: x 1 λ (x)= , x N ; λ (x)= x 1+2 , x N. 1 x +2 ∈ 0 −1 − x ∈ Here R = R− = 1 > R+ = 0, α = 1, σ1 = 1, σ2 = 2. Using the same Lyapunov function constructed in the proof of [56, Theorem− 4.4], it can be shown that there exists a uniquely ergodic QSD. By Theorem 4.4(iv), the tail of the QSD decays no more slowly than a CMP distribution. Indeed, the QSD is given by 1 1 ν(x)= ,T (x)= , x N. (x 1)!(x + 1) ν x! ∈ − The tail of the QSD decays like a Poisson distribution. 14 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF Example 4.11. Consider a BDP on N0: λ (x)= x2, λ (x)= x2 + x, x N . 1 −1 ∈ 0 Here R = R− = R+ = 2 > 1, and α = 0. Corollary 4.5 states that any QSD (if it exists) is heavy-tailed and its tail decays no faster than a Zeta-like distribution. Indeed, a QSD of the process is given by 1 1 ν(x)= ,T (x)= , x N. x(x + 1) ν x ∈ Example 4.12. Consider a quadratic BDP on N0: λ (x)= x(x + 3)/2, λ (x)= x(x + 1), x N . 1 −1 ∈ 0 Then R− = R+ = R =2 > 1, α = 1/2. Hence there exists a uniquely ergodic QSD [56]. By Corollary 4.5, this QSD decays exponentially.− Indeed, the QSD is given by ν(x)=2−x,T (x)=2−x+1, x N. ν ∈ 5. Applications In this section, we apply the results on asymptotic tails to diverse models in biology. We emphasize that for all models/applications, the transition rate functions are real analytic APLH on a subset of N0, and thus σ1 = 1 and σ2 = 2, which we will not further mention explicitly. 5.1. Biochemical reaction networks. In this section, we apply the results of Section 5 to some examples of stochastic reaction networks (SRNs) with mass-action kinetics. These are use to describe interactions of constituent molecular species with many applications in systems biology, biochemistry, genetics and beyond [23, 48]. A SRN with mass-action kinetics is a CTMC on Nd (d 1) encoded by a labeled directed graph [2]. We concentrate on SRNs 0 ≥ on N0 with one species (S). In this case the graph is composed of reactions (edges) of the form κ nS mS, n,m N0 (n molecules of species S is converted into m molecules of the same species),−−→ encoding a∈ jump from x to x+m n with propensity λ(x)= κx(x 1) . . . (x n+1), κ> 0. Note that multiple reactions might− result in the same jump vector.− − In general little is known about the stationary distributions of a reaction network, let alone the QSDs, provided either such exist [3, 27, 29, 30]. Special cases include complex balanced networks (in arbitrary dimension) which have Poisson product-form distributions [3, 11], reaction networks that are also birth-death processes, and reaction networks with irreducible components, each with a finite number of states. Example 5.1. To show how general the results are we consider two SRNs, none of which are birth-death processes. (i) Consider a reaction network with a strongly connected graph [56]: κ3 κ κ S 1 2S 2 3S For this reaction network, Ω = 1, 2 , and { − } λ (x)= κ x + κ x(x 1), λ (x)= κ x(x 1)(x 2). 1 1 2 − −2 3 − − Hence, α = κ3 < 0. It is known that there exists s a unique exponentially ergodic stationary distribution−π on N [56]. (The state 0 is neutral [55].) By Theorem 4.1, π 1+ 1−. Hence ∈P1 ∩P1 π is light-tailed and Tπ decays as fast as a Poisson distribution. since ω+ = ω∗ = 1, R+ = 2, 2 R− = 3. However, the stationary distribution is generally not Poisson. If κ2 = κ1κ3, then the reaction network is complex-balanced, hence the stationary distribution is Poisson [3]. If the parameter identity is not fulfilled then the distribution cannot be Poisson in this case [11]. THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 15 (ii) Consider a similar reaction network including direct degradation of S [56]: κ3 κ κ κ ∅ 4 S 1 2S 2 3S The threshold parameters are the same as in (i), and it follows from [56] that the reaction network has a uniformly exponentially ergodic QSD ν. By Theorem 4.4, Tν decays like a CMP distribution. Example 5.2. The following bursty Schl¨ogl model was proposed in [22]: κ0 κ κ ∅ ↽⇀ S, 3S 3 2S 2 (2 + j)S, −−κ−−−−1 −−→ −−→ where j N. When j = 1, it reduces to the classical Schl¨ogl model. ∈ The associated process has a unique ergodic stationary distribution π on N0 [56]. Bi- furcation with respect to patterns of the ergodic stationary distribution is discussed in [22], based on a diffusion approximation in terms of the Fokker-Planck equation. Using the res- ults established in this paper the tail distribution can be characterized rigorously. In fact 1+ 1− π − . Hence π is light-tailed and T decays like a CMP distribution. ∈Pj 1 ∩P1 π Proof. We have Ω = 1, 1 j . The ergodicity follows from [56, subsection 4.3] as a special {− }∪{ } case. It is straightforward to verify that (A1)-(A3) are all satisfied. Moreover, R+ = 2, R− = 3. Then the conclusion follows from Theorem 4.1. Example 5.3. Consider the following one-species S-system modelling a gene regulatory net- work [16]: (κ , ξ ) (κ , ξ ) ∅ 1 1 S 2 2 3S −−−−−→ ←−−−−− with the following generalized mass action kinetics (GMAK): Γ(x + ξ ) Γ(x + ξ ) λ (x)= κ 1 , λ (x)= κ 2 , x N , 1 1 Γ(x) −2 2 Γ(x) ∈ 0 where κ1,κ2 > 0 are the reaction rate constants, and ξ2 > ξ1 > 0 are the indices of GMAK. By Stirling’s formula log Γ(x) = (x 1/2)log x x + log √2π + O(x−1), − − hence λ is APLH with (ξ , ξ 1, ξ 2) and λ is APLH with (ξ , ξ 1, ξ 2). Then 1 1 1 − 1 − −2 2 2 − 2 − R− = ξ2 > R+ = ξ1, ω∗ = 1, ω− = 2 and ω+ = 1. Using the same Lyapunov function constructed in the proof of [56, Theorem− 4.4], it can be shown that there exists a uniquely 1+ 1− ergodic stationary distribution π on N0 with support N. By Theorem 4.1, π . ∈Pξ3−ξ1 ∩Pξ3−ξ1 5.2. An extended class of branching processes. Consider an extended class of branching processes on N0 [15] with transition rate matrix Q = (qx,y)x,y∈N0 : r(x)µ(y x + 1), if y x 1 0 and y = x, r(x)(1− µ(1)), if y ≥= x − 1,≥ 6 − − ≥ qx,y = q0,y, if y > x =0, q0, if y = x =0, −0, otherwise, where µ is a probability measure on N , q = q , and r(x) is a positive finite function 0 0 y∈N 0,y on N0. Assume P (H1) µ(0) > 0, µ(0) + µ(1) < 1. (H2) N q0,yy< , M = N kµ(k) < . y∈ ∞ k∈ 0 ∞ (H3) rP(x) is a polynomial ofP degree R 1 for large x. ≥ 16 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF The tail asymptotics of infinite stationary measures in the null recurrent case is investigated in [37] under (H1)-(H2) for general r. Here we assume r is polynomial (H3). The following is a consequence of the results of Section 5. Theorem 5.4. Assume (H1)-(H3), Y =0, and that µ has finite support. 0 6 (i) Assume q0 > 0. Then there exists an ergodic stationary distribution π on N0 if (i-1) M < 1 or (i-2) M =1 and R> 1. Moreover, Tπ decays like a geometric distribution if (i-1) holds while like a Zeta distribution if (i-2) holds. (ii) Assume q0 =0. Then there exists an ergodic QSD ν on N if (i-1) M < 1 and R> 1 or (i-2) M =1 and R> 2. Moreover, Tν decays like a geometric distribution if (ii-1) holds while no faster than a Zeta distribution if (ii-2) holds. Proof. For all k Ω, let ∈ r(x)µ(k + 1), if x N, λk(x)= ∈ (q0k, if x =0. By (H1), µ(k) > 0 for some k N. Hence regardless of q , by positivity of r, (A1)- ∈ 0 (A3) are satisfied with Ω− = 1 and Ω+ = j N: j +1 supp µ or q0j > 0 . Let R R−1 R−2 {− } { ∈ ∈ } r(x) = ax + bx + O(x ) with a> 0. It is straightforward to verify that R+ = R− = R, α = a(M 1). The ergodicity follows from [56], and the tail asymptotics follow from Theorems 4.1 −and 4.4. 5.3. Stochastic population processes under bursty reproduction. Two stochastic pop- ulation models with bursty reproduction are investigated in [8]. The first model is a Verhulst logistic population process with bursty reproduction. The process Yt is a CTMC on N0 with transition rate matrix Q = (qx,y)x,y∈N0 satisfying: cµ(j)x, if y = x + j, j N, c 2 ∈ qx,y = x + x, if y = x 1 N , K − ∈ 0 0, otherwise, where c > 0 is the reproduction rate, K N is the typical population size in the long-lived metastable state prior to extinction [8], and∈ µ is the burst size distribution. Approximations of the mean time to extinction and QSD are discussed in [8] against various different burst size distributions of finite mean (e.g., Dirac measure, Poisson distribution, geometric distribution, negative-binomial distribution). The existence of an ergodic QSD for this population model is established in [56]. Nevertheless, the tails of QSD is not addressed therein. Theorem 5.5. Assume µ has a finite support. Let ν be the unique ergodic QSD on N trapped to zero for the Verhulst logistic model Yt. Then Tν decays like a CMP distribution. c 2 Proof. We have Ω = 1 supp µ, λ−1(x)= K x + x, λk(x)= cµ(k)x, for k supp µ and x N. Since µ has a{− finite} ∪support, (A1)-(A3) are satisfied. Moreover, since supp∈ µ = ∅, we ∈ 6 have R− = 2 and R+ = 1. Again, the ergodicity result follows from [56]. The tail asymptotics follow directly from Theorem 4.1. THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 17 6. Proof of Theorem 4.1 Let λ (x) ∈ ω lim ω Aj , if j = ω +1,..., 0, x→∞ xR− − P λω (x) αj = ω∈Aj limx→∞ R+ , if j =1,...,ω+, x 0, P otherwise, R− λω (x)−αj x ω∈Aj lim − , if j = ω +1,..., 0, x→∞ xR− σ1 − R P λω (x)−αj x + γj = ω∈Aj limx→∞ R+−σ1 , if j =1,...,ω+, x 0, P otherwise. 0 ω+ A3 Note that β = α0. From Lemma B.1, α− = ω∗ j=ω−+1 αj , α+ = ω∗ j=1 αj . By( ), R− −σ1 −σ2 x (αj + γj x + O(xP )), if j = ω− +1P,..., 0, (6.1) λω(x)= R+ −σ1 −σ2 (x (αj + γj x + O(x )), if j =1,...,ω+. ωX∈Aj Since λ (x) = sgn(j) λ (x) λ (x) , j = ω ,..., 1, 1,...,ω , jω∗ ω − ω − − + ωX∈Aj ω∈XAj+1 we have R− −σ1 −σ2 x ((αj+1 αj ) + (γj+1 γj )x + O(x )), if j = ω−,..., 1, λjω∗ (x)= − − − xR+ ((α α ) + (γ γ )x−σ1 + O(x−σ2 )), if j =1,...,ω . ( j − j+1 j − j+1 + Since (A1)-(A2) imply x : λ (x)=0 is finite, then from Corollary 3.3, it ∩ω∈Ω { ∈ Y ω } follows that both Ω− = ∅ and Ω+ = ∅ since supp π is unbounded. Hence α− α0 > 0, α α > 0, and 6 <ω <ω <6 . ≥ + ≥ 1 −∞ − + ∞ For the ease of exposition and w.o.l.g., we assume throughout the proof that ω∗ = 1 (recall the argument in the proof of Theorem 3.2). Hence N0 + b N0 for some b N0 by Lemma A.1. ⊆Y ⊆ ∈ Most inequalities below are based on the identities in Theorem 3.2 and Corollary 3.8. Therefore, we use LHS (RHS) with a label in the subscript as shorthand for the left (right) hand side of an equation with the given label. The claims that ∆ = 0 implies σ1 < 1 and that σ1 = 1 implies ∆ > 1 are proved in (iii)-(vi) Step I below. We first show α 0. Suppose by means of contradiction that α > 0. Then either (1) ≤ R+ > R− or (2) R+ = R− and α+ > α− holds. Define the auxiliary function f (x)= λ (x j), j = ω +1,...,ω . j ω − − + ωX∈Aj From (6.1) it follows that R− −σ1 −1 − min{σ2,σ1+1} x (αj + γj x αj jRx + O(x )), if j = ω− +1,..., 0, fj(x)= − xR+ (α + γ x−σ1 α jRx−1 + O(x− min{σ2,σ1+1})), if j =1,...,ω . ( j j − j + Let −R− x fj (x) αj , if j = ω− +1,..., 0, βj (x)= − x−R+ f (x) α , if j =1,...,ω . ( j − j + Then there exist N , N N with N >N ,N such that for all x N , 3 4 ∈ 3 1 4 ≥ 3 (6.2) β (x) N x−σ1 , j = ω +1,...,ω , | j |≤ 4 − + 18 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF From Theorem 3.2(iii), we have 0 ω+ (6.3) xR−−R+ (α + β (x)) π (x j)= (α + β (x)) π (x j) , j j − j j − j=ω−+1 j=1 X X Since R R and T (x) 1 for all x N , summing up in (6.3) from x to infinity yields − ≤ + π ≤ ∈ 0 ∞ 0 ∞ ω+ (6.4) yR−−R+ (α + β (y)) π (y j)= (α + β (y)) π (y j) . j j − j j − y=x j=ω−+1 y=x j=1 X X X X R−−R+ In the light of the monotonicity of Tπ(x) and x , it follows from (6.2) that there exists C = C(N ) > 0 and N N with N N such that for all x N , 4 5 ∈ 5 ≥ 3 ≥ 5 ∞ 0 LHS yR−−R+ α + N y−σ1 π (y j) (6.4) ≤ j 4 − y=x j=ω−+1 X X 0 ∞ xR−−R+ α + N x−σ1 π (y j) ≤ j 4 − j=ω−+1 y=x X X 0 = xR−−R+ α + N x−σ1 T (x j) j 4 π − j=ω−+1 X 0 xR−−R+ T (x) α + N x−σ1 ≤ π j 4 j=ω−+1 X xR−−R+ T (x) α + Cx−σ1 . ≤ π − Similarly, with a possibly larger C and N , for all x N , 5 ≥ 5 ∞ ω+ RHS α N y−σ1 π (y j) (6.4) ≥ j − 4 − y=x j=1 X X ω+ ∞ α N x−σ1 π (y j) ≥ j − 4 − j=1 y=x X X ω+ = α N x−σ1 T (x j) j − 4 π − j=1 X ω+ T (x 1) α N x−σ1 ≥ π − j − 4 j=1 X α Cx−σ1 T (x 1), ≥ + − π − which further implies that for all x large enough, −σ1 Tπ(x) α+ Cx 1 xR+−R− − > 1, ≥ T (x 1) ≥ α + Cx−σ1 π − − since either (1) R+ > R− or (2) R+ = R− and α+ > α− holds. This contradiction shows that α 0. ≤ Next, we provide asymptotics of Tπ(x) case by case. (i) R− > R+. Recall Stirling’s formula for the Gamma function [1]: log Γ(x)= x log x x + O(log x), − THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 19 where log is the natural logarithm. Based on this Stirling’s formula, it suffices to prove that there exists C > 0 such that x+Cx1−(R−−R+)+O(log x) e R+−R− α1 (6.5) Tπ(x) & Γ(x) , α0 e −1 1−σ1 ω+ x+Cx +O(log x) −1 R+−R− α+ R+−R− (6.6) Tπ(x) . Γ(xω+ ) ω+ . α0 e Next, we prove (6.5) and (6.6) one by one. We first show (6.5). Recall that (A2) ensures that there exists N N such that λω is a strictly positive non-decreasing polynomial on N + N for all ω Ω.∈ Moreover, A = ω 0 ∈ j { ∈ Ω− : ω < j if j 0, and Aj = ω Ω+ : ω j if j > 0. It follows from Corollary 3.8 that for all x }N + N≤ ω , { ∈ ≥ } ∈ 0 − − T (x) λ (x)+ λ (x 1) π ω ω − ωX∈A0 ωX∈A1 −1 + T (x j) λ (x j) λ (x (j + 1)) π − · ω − − ω − j=ω− X ωX∈Aj ω∈XAj+1 ω+ (6.7) = T (x j) λ (x j) λ (x (j + 1)) . π − ω − − ω − j=1 X ωX∈Aj ω∈XAj+1 Furthermore, we have the following estimates for both sides of the above equality: −1 LHS = T (x) λ (x)+ λ (x 1) + T (x j) (6.7) π ω ω − π − j=ω− ωX∈A0 ωX∈A1 X λ (x (j +1))+ λ (x j) λ (x (j + 1)) · − j − ω − − ω − ωX∈Aj −1 T (x) λ (x)+ λ (x 1) + λ (x j) λ (x (j + 1)) ≤ π ω ω − ω − − ω − ω∈A ω∈A j=ω− ω∈A X0 X1 X Xj R− R+ R−−1 = Tπ(x) α0x + α1x +O x R− −σ˜ = Tπ(x)x α0 +O x , whereσ ˜ = min 1, R R > 0. { − − +} By the monotonicity of λω, ω+ RHS T (x j) λ (x j) λ (x j) (6.7) ≥ π − ω − − ω − j=1 X ωX∈Aj ω∈XAj+1 ω+ T (x 1) λ (x j) ≥ π − j − j=1 X = T (x 1) λ (x) λ (x) λ (x ω) π − ω − ω − ω − ω∈A ω∈A X1 X1 T (x 1) α xR+ +O xR+−1 . ≥ π − 1 20 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF Then there exist N1 >N2 > 0 such that for all x N1, (6.8) ≥ −1 Tπ(x) α1 + O(x ) α1 α1 xR+−R− = xR+−R− + O(x−σ˜ ) xR+−R− 1 N x−σ˜ . T (x 1) ≥ α + O(x−σ˜ ) α ≥ α − 2 π − 0 0 0 Hence if 0 < R− R+ < 1, thenσ ˜ = R− R+ < 1, and there exists C = C(N1,N2) > 0 such that for all x N− , − ≥ 1 e e x−N1 α1 N2 T (x) T (N 1) (x j)R+−R− 1 π ≥ π 1 − α − − (x j)σ˜ j=0 0 Y − x x α1 R+−R− N2 = Tπ(N1 1) j 1 σ˜ − α0 − j jY=N1 jY=N1 x α x+1−N1 Γ (x + 1)R+−R− & 1 exp 2N j−σ˜ R −R− 2 α0 Γ(N1) + − jX=N1 x−Cx1−σ˜ α1 & Γ(x)R+−R− xR+−R− , α0 e since 1 x−1 exp( 2x−1) for large x, and we employ the fact that T (N 1) > 0 since − ≥ − π 1 − supp π = is unbounded. Hence (6.5) holds. Similarly, if R− R+ 1, then δ = 1, and analogousY arguments can be applied. − ≥ Next we show (6.6). Rewrite (6.3), e 0 ω+ (6.9) (α + β (x)) π (x j)= xR+−R− (α + β (x)) π (x j) , j j − j j − j=ω−+1 j=1 X X Summing up in (6.9) from x to infinity yields ∞ 0 ∞ ω+ (6.10) (α + β (y)) π (y j)= yR+−R− (α + β (y)) π (y j) . j j − j j − y=x j=ω−+1 y=x j=1 X X X X It follows from (6.2) that there exists C = C(N ) > 0 and N N such that for all x N , 4 5 ∈ ≥ 5 ∞ 0 LHS α N y−σ1 π (y j) (6.10) ≥ j − 4 − y=x j=ω−+1 X X 0 ∞ α N x−σ1 π (y j) ≥ j − 4 − j=ω−+1 y=x X X 0 = α N x−σ1 T (x j) j − 4 π − j=ω−+1 X α Cx−σ1 T (x), ≥ 0 − π THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 21 ∞ ω+ RHS yR+−R− α + N y−σ1 π (y j) (6.10) ≤ j 4 − y=x j=1 X X ω+ ∞ xR+−R− α + N x−σ1 π (y j) ≤ j 4 − j=1 y=x X X ω+ = xR+−R− α + N x−σ1 T (x j) j 4 π − j=1 X ω+ xR+−R− T (x ω ) α + N x−σ1 ≤ π − + j 4 j=1 X xR+−R− T (x ω ) α + Cx−σ1 , ≤ π − + + which together further imply that −σ1 Tπ(x) α+ + Cx α+ xR+−R− = xR+−R− + O(x−σ1 ) . T (x ω ) ≤ α Cx−σ1 α π − + 0 − 0 The remaining arguments are analogous to the arguments for (6.5). (ii) R− = R+ and α− > α+. Analogous to (i), we will show that there exist real constants δ+,δ− and C > 0 such that for all δ>δ+ and δ <δ−, x+Cx1−σ1 +O(log x) e α+ (6.11) Tπ(x) & , α− e −1 1−σ1 (ω+−ω−−1) x+Cx +O(log x) α+ (6.12) Tπ(x) . . α− e We first prove (6.11). Since R = R− = R+, R −σ1 fj(x)= x (αj + βj (x)), βj (x)=O(x ), j = ω− +1,...,ω+. Moreover, α = α α < 0 implies that + − − 0 ω+ αj < αj . j=ω−+1 j=1 X X From (6.14) it follows that 0 0 ω+ ω+ π (x j) α + π (x j) β (x)= π (x j) α + π (x j) β (x). − j − j − j − j j=ω−+1 j=ω−+1 j=1 j=1 X X X X Summing up the above equality from x to yields ∞ 0 ∞ 0 ∞ π (y j) α + π (y j) β (y) − j − j j=ω−+1 y=x j=ω−+1 y=x X X X X ω+ ∞ ω+ ∞ = π (y j) α + π (y j) β (y). − j − j j=1 y=x j=1 y=x X X X X Since each double sum in the above equality is convergent, we have 0 ∞ ω+ ∞ 0= α (π(y) π (y j))+ α (π (y j) π(y)) j − − j − − j=ω−+1 y=x j=1 y=x X X X X 22 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF 0 ω+ α α T (x) − j − j π j=ω−+1 j=1 X X 0 ∞ + (π(y)β (y + j) π (y j) β (y)) j − − j j=ω−+1 y=x X X ω+ ∞ + (π (y j) β (y) π(y)β (y + j)) − j − j j=1 y=x X X ω+ ∞ 0 ∞ + π(y)β (y + j) π(y)β (y + j). j − j j=1 y=x j=ω−+1 y=x X X X X This further yields the following equality 0 ∞ ω+ ∞ (α α )T (x)+ π(y)β (y + j) π(y)β (y + j) − − + π j − j j=ω−+1 y=x j=1 y=x X X X X −1 −j−1 ω+ j (6.13) = π (x + ℓ) f (x + j + ℓ)+ π (x ℓ) f (x + j ℓ), j − j − j=ω−+1 j=1 X Xℓ=0 X Xℓ=1 −R e e where fj (x) = αj + βj(x) = x fj (x) 0, j = ω− +1,...,ω+. From (6.2), it follows that there exist C,N > 0 such that ≥ e0 ∞ ω+ ∞ ∞ π(y)β (y + j) π(y)β (y + j) C π(y)y−σ1 Cx−σ1 T (x), j − j ≤ ≤ π j=ω−+1 y=x j=1 y=x y=x X X X X X for x N. Hence ≥ LHS (α α )+ Cx−σ1 T (x). (6.13) ≤ − − + π Using Fubini’s theorem, we have: 0 j−1 ω+ ω+ (6.14) RHS = π(x j) f (x + ℓ j)+ π(x j) f (x + ℓ j). (6.13) − ℓ − − ℓ − j=ω−+2 − j=1 X ℓ=Xω +1 X Xℓ=j Hence further choosing larger N and C, we havee for all x N, e ≥ ω+ ω+ RHS π (x 1) f (x + j 1) = π (x 1) (α + β (x + j 1)) (6.13) ≥ − j − − j j − j=1 j=1 X X (6.15) (α Cx−σ1 )eπ (x 1) , ≥ + − − This implies from π(x 1) = T (x 1) T (x) that − π − − π T (x) α Cx−σ1 π + − . T (x 1) ≥ α π − − Using similar arguments as in the proof of (i), one obtains (6.11). Next we show (6.12). We establish the reverse estimates for both sides of (6.13). Similarly, there exists some C,N > 0 such that for all x N, ≥ LHS (α α Cx−σ1 )T (x) (α α Cx−σ1 )T (x ω 1), (6.13) ≥ − − + − π ≥ − − + − π − − − RHS (α + Cx−σ1 ) (T (x ω ) T (x ω 1)) . (6.13) ≤ + π − + − π − − − This implies that for a possibly larger C, N, for all x N, ≥ T (x ω 1) α + Cx−σ1 π − − − + . T (x ω ) ≤ α π − + − THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 23 The remaining arguments are similar to those in the proof of (i). (iii)-(vi) α = 0. Hence α+ = α−. Recall γ, if σ < 1, ∆= α−1 − 1 + · γ + Rϑ, if σ =1. (− 1 Let δ = ∆(ω ω 1)−1. For j = ω +1,...,ω , + − − − − + γj , if σ1 < 1, rj = γ jRα , if σ =1, ( j − j 1 ϑ (x)= β (x) r x−σ1 . j j − j Hence we have −σ2 ϑj (x)=O(x ), j = ω− +1,...,ω+, where min 1, σ , if σ < 1, σ = 2 1 2 σ , { } if σ =1. 2 1 Let η = σ σ and ε = min σ , η . Hence 0 < ε η 1. If σ < 1, then η 1 σ . If 2 − 1 { 1 } ≤ ≤ 1 ≤ − 1 σ1 = 1, then ε = η. To show (iii)-(vi), it suffices to prove that there exists C > 0 such that (6.16) ∆ 1−σ1 1−σ1−ε exp x +O x + log x , if σ1 < 1, σ1 + ε =1, ∆ > 0 − 1−σ1 6 ∆ 1−σ1 exp x + O(log x) , if σ1 < 1, σ1 + ε =1, ∆ > 0 − 1−σ1 −(∆−1) x , if σ1 =1, ∆ > 0, Tπ(x) & C 1−σ2 σ exp x +O x + log x , if σ2 < 1, σ1 + σ2 =1, ∆=0 − 1−σ2 6 C 1−σ2 exp x + O(log x) , if σ2 < 1, σ1 + σ2 =1, ∆=0 − 1−σ2 −C x , if σ2 1, ∆=0, ≥ where σ = max 1 σ1 σ2, 0 , and if ∆ > 0, then (6.17) { − − } δ 1−σ1 max{1−σ1−ε,0} exp x +O x + log x , if σ1 < 1, σ1 + ε =1, − 1−σ1 6 T (x) . δ 1−σ1 π exp x + O(log x) , if σ1 < 1, σ1 + ε =1, − 1−σ1 − max{δ−1,0} x , if σ1 =1. To show (6.16) and (6.17) for a probability distribution µ on N0, define its weighted tail distribution on N0 as ∞ W : N [0, 1], x y−σ1 µ(y). µ → 7→ y=x X In the following, we will show there exist constants C > σ1 such that (6.18) −∆ 1−σ1 1−σ1 −ε exp x +O x , if σ1 < 1, σ1 + ε =1, ∆ > 0, 1−σ1 6 −∆ 1−σ1 exp x + O(log x) , if σ1 < 1, σ1 + ε =1, ∆ > 0, 1−σ1 −∆ x , if σ1 =1, ∆ > 0, Wπ(x) & −C 1−σ2 max{1−σ1−σ2,0} exp x +O x , if σ2 < 1, σ1 + σ2 =1, ∆=0, 1−σ2 6 −C 1−σ2 exp x + O(log x) , if σ2 < 1, σ1 + σ2 =1, ∆=0, 1−σ2 −C x , if σ2 1, ∆=0, ≥ 24 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF with ∆ > 1, when σ1 = 1. Moreover, if ∆ > 0, then −δ 1−σ1 max{1−σ1−ε,0} exp x + O(x ) , if σ1 < 1, σ1 + ε =1, 1−σ1 6 (6.19) W (x) . −δ 1−σ1 π exp x + O(log x) , if σ1 < 1, σ1 + ε =1, 1−σ1 −δ x , if σ1 =1. Then we will prove (6.16) and (6.17) based on (6.18) and (6.19). Step I. Prove (6.18) and (6.19). Here we will also show ∆ 0, and in particular ∆ > 1 when ≥ σ1 = 1. We first show (6.18). Since α+ = α−, 0 ∞ −σ1 LHS(6.13) = π(y) rj (y + j) + ϑj (y + j) j=ω−+1 y=x X X ω+ ∞ π(y) r (y + j)−σ1 + ϑ (y + j) − j j j=1 y=x X X 0 ω+ = r r W (x) j − j π j=ω−+1 j=1 X X 0 ∞ + π(y) ϑ (y + j)+ r ((y + j)−σ1 y−σ1 ) j j − j=ω−+1 y=x X X ω+ ∞ π(y) ϑ (y + j)+ r ((y + j)−σ1 y−σ1 ) . − j j − j=1 y=x X X By Lemma B.1, 0 ω+ r r = α ∆. j − j + j=ω−+1 j=1 X X Moreover, 0 ∞ π(y) ϑ (y + j)+ r ((y + j)−σ1 y−σ1 ) j j − j=ω−+1 y=x X X ω+ ∞ π(y) ϑ (y + j)+ r ((y + j)−σ1 y−σ1 ) − j j − j=1 y=x X X ∞ . π(y)y−σ1−η x−ηW (x). ≤ π y=x X Since RHS(6.13) 0 for all large x, we have ∆ 0. From (6.15) it≥ follows that there exist C, N≥ N such that for all x N, ∈ ≥ LHS α ∆+ Cx−η W (x), (6.13) ≤ + π while RHS α (1 Cx−σ1 )π(x 1) (6.13) ≥ + − − = α (1 Cx−σ1 )(x 1)σ1 (W (x 1) W (x)) + − − π − − π = α (xσ1 C σ xσ1−1 + O(xσ1−2))(W (x 1) W (x)). + − − 1 π − − π Further choosing larger N and C, by the monotonicity of W , for all x N, π ≥ (xσ1 C σ xσ1−1 + O(xσ1−2))W (x 1) − − 1 π − xσ1 C + ∆ σ xσ1−1 + Cx−η + O(xσ1−2) W (x). ≤ − − 1 π THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 25 If σ < 1, then η 1 σ , and hence η + σ 2 σ . If σ = 1, then ε = η. Recall 1 ≤ − 1 1 − ≤ − 1 1 σ2 = η + σ1. Then we have W (x) xσ1 C σ xσ1−1 + O(xσ1−2) π − − 1 W (x 1) ≥ xσ1 C +∆+ Cx−η σ xσ1 −1 + O(xσ1−2) π − − − 1 1 ∆x−σ1 (1 + O(x−ε)), if ∆ > 0, = − 1 Cx−σ2 (1 + O(xmax{−σ1,σ2−2})), if ∆=0. − First assume ∆ > 0. Since ε 1, by Euler-Maclaurin’s formula, ≤ x Wπ(x) log log(1 ∆j−σ1 + O(j−σ1−ε)) Wπ(N 1) ≥ − − jX=N −∆ x1−σ1 + O(xmax{0,1−σ1−ε}), if σ < 1, σ + ε =1, 1−σ1 1 1 −∆ 1−σ1 6 = x + O(log x), if σ1 < 1, σ1 + ε =1, 1−σ1 ∆log x + O(1), if σ =1, − 1 which implies that −∆ 1−σ1 1−σ1−ε exp x + O(x ) , if σ1 < 1, σ1 + ε =1, 1−σ1 6 W (x) & −∆ 1−σ1 π exp x + O(log x) , if σ1 < 1, σ1 + ε =1, 1−σ1 −∆ x , if σ1 =1, i.e., (6.18) holds. Moreover, since xσW (x) T (x) 0 as x , we have 1 π ≤ π → → ∞ 0, if σ < 1, ∆ > 1 1, if σ =1. 1 Now assume ∆ = 0, then −C 1−σ2 max{1−σ1−σ2,0} exp x + O(x ) , if σ2 =1, σ1 + σ2 =1, 1−σ2 6 6 W (x) & −C 1−σ2 π exp x + O(log x) , if σ2 =1, σ1 + σ2 =1, 1−σ2 −C 6 x , if σ2 =1, where we use the fact that σ1 + σ2 = 1 implies 0 < σ1, σ2 < 1 and σ1 + σ2 = 1. Moreover, σ also due to x1 Wπ(x) Tπ(x) 0 as x , we have σ2 1, which implies that σ1 < 1. In addition, C > σ when≤ σ = 1,→ i.e., σ <→1 ∞ σ . Hence for≤ some C > σ , 1 2 1 ≤ 2 1 −C 1−σ2 max{1−σ1−σ2,0} exp x + O(x ) , if σ2 < 1, σ1 + σ2 =1, 1−σ2 6 W (x) & −C 1−σ2 π exp x + O(log x) , if σ2 < 1, σ1 + σ2 =1, 1−σ2 −C x , if σ2 1. ≥ Next we show (6.19) by establishing the reverse estimates for both sides of (6.13). From (6.14) it follows that there exist positive constants N and C (i =1, 2) such that x N, i ≥ LHS α ∆ Cx−η W (x) α ∆ Cx−η W (x (ω + 1)), (6.13) ≥ + − π ≥ + − π − − whereas ω+ RHS α 1+ C x−σ1 π (x j) (6.13) ≤ + 1 − j=ω−+2 X ω+ α (xσ1 + C ) π (x j) (x j)−σ1 ≤ + 2 − − j=ω−+2 X = α (xσ1 + C ) (W (x ω ) W (x (ω + 1))), + 2 π − + − π − − 26 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF Hence, when ∆ > 0, then for all x N, ≥ σ1 Wπ(x (ω− + 1)) x + C2 − =1 ∆x−σ1 + O(x−σ1−ε). W (x ω ) ≤ xσ1 + C + ∆ Cx−η − π − + 2 − Analogous to the above analysis, one might show −δ 1−σ1 max{1−σ1−ε,0} exp x + O(x ) , if σ1 < 1, σ1 + ε =1, 1−σ1 6 W (x) . −δ 1−σ1 π exp x + O(log x) , if σ1 < 1, σ1 + ε =1, 1−σ1 −δ x , if σ1 =1, −1 in particular, one can further choose δ 1 (not necessarily δ = ∆(ω+ ω− 1) ) when σ = 1, also due to xσ W (x) T (x) ≥0 as x . − − 1 1 π ≤ π → → ∞ Moreover, one can always show W (x) x−σ1 T (x) x−σ1 , hence (6.19) also holds when π ≤ π ≤ σ1 = 1. Step II. Prove (6.16) and (6.17) based on (6.18) and (6.19). −σ1 Since Wπ(x) x Tπ(x), (6.16) follows directly from (6.18). Next, we prove≤ (6.17) based on (6.19). Recall that π(x) xσW (x). ≤ 1 π Assume ∆ > 0. We only prove the case σ1 < 1 and σ1 + ε = 1. The other two cases can be proved using analogous arguments. Then there exist N N and C > σ such that ∈ 1 1 −δ 1−σ1 exp y + C1 log y is decreasing on [N, + ), and for all x N, 1−σ1 ∞ ≥ ∞ ∞ T (x)= π(y) yσW (y) π ≤ 1 π y=x y=x X X ∞ −δ 1−σ1 . exp y + C1 log y 1−σ1 y=x X∞ −δ 1−σ1 . exp y + C1 log y dy 1−σ1 x−1 Z ∞ . yC1+σ1 d exp −δ y1−σ1 1−σ1 x−1 − Z = (x 1)C1+σ1 exp −δ (x 1)1−σ1 − 1−σ1 − ∞ −δ 1−σ1 + (C1 + σ1) exp y + (C1 + σ1 1)log y dy 1−σ1 x−1 − Z (x 1)C1+σ1 exp −δ (x 1)1−σ1 ≤ − 1−σ1 − ∞ σ1−1 −δ 1−σ1 + (C1 + σ1)(x 1) exp y + C1 log y dy, 1−σ1 − x−1 Z which further implies that for all x N, ≥ (x 1)C1+σ1 exp −δ (x 1)1−σ1 − 1−σ1 − Tπ(x) . 1+O (x 1)σ1−1 − . exp −δ (x 1)1−σ1 + O(log x) 1−σ1 − = exp −δ x1−σ1 + O(log x) . 1−σ1 This shows (6.17) in that case. THE ASYMPTOTIC TAILS OF LIMIT DISTRIBUTIONS OF CTMCS 27 7. Proof of Theorem 4.4 Again (A3) implies supp ν = ∂c [17], which is unbounded by (A2). Comparing the identities for stationary distributions and QSDs, the unique difference comes from an extra term on the RHS of the identity of QSDs with coefficient θν > 0. This makes the identity in Theorem 3.2(3) for stationary distributions to be an inequality with its LHS greater than its RHS for QSDs. Hence all arguments in the proof of Theorem 4.4 establishing α 0 as well as the lower estimates for T (the tail of the stationary distribution) carry over ≤ π to Tν. Next, we show R 0. The proof is in a similar spirit to that for α 0. Since α 0, ≥ ≤ ≤ R− = R. Again, assume w.o.l.g. that ω∗ = 1 such that ∂ contains all large positive integers by Proposition A.1. From Theorem 3.4(3), similar to (6.3), we have for all large x, 0 xR(α + Cx−σ1 )T (x) xR (α + β (x)) ν (x j) − ν ≥ j j − j=ω−+1 X ω+ = θ T (x)+ xR+ (α + β (x)) ν (x j) ν ν j j − j=1 X θ T (x), ≥ ν ν which yields xR(α + Cx−σ1 ) θ 0, − − ν ≥ This shows R 0, since θ > 0. Moreover, if R = 0, then α θ . The claim that ≥ ν − ≥ ν R− = R+ = 0 implies α θν is proved below in (vii). Similar to (6.7) and the≤− inequality (6.8) based on it, one can also obtain α θ if R =0 0 ≥ ν and R− > R+. Moreover, there exists C > α1 > 0 such that for all large x, −1 R+−R− α1−Cx Tν (x) x −R− −σ˜ , if R− > R+, α0−θν x +Cx (7.1) −1 T (x 1) α1−Cx ν ≥ ( −R R−1 , if R− = R+, − α0+α1−θν x +Cx where we recallσ ˜ = min 1, R− R+ . Similar to (6.10), we establish{ − } (7.2) ∞ 0 ∞ ∞ ω+ (α + β (y)) ν (y j)= θ y−R− T (y)+ yR+−R− (α + β (y)) ν (y j) . j j − ν ν j j − y=x j=ω−+1 y=x y=x j=1 X X X X X ∞ −R− Since LHS(7.2) is finite, we have y=x y Tν (y) is also finite. Furthermore, by similar analysis as in the proof of Theorem 4.1 that there exists C > 0 such that for all large x N, P ∈ ∞ θ y−R− T (y) (α + Cx−σ1 )T (x) 0, if x . ν ν ≤ − ν → → ∞ y=x X Step I. Establish lower estimates for Tν based on the above inequality, using similar asymp- totic analysis demonstrated repeatedly in the proof of Theorem 4.4. (1) R− = R > R+. R = 0 (cases (i)-(iii)). Then α θ . If α >θ , then there exists C > 0 such that • 0 ≥ ν 0 ν α0−θν 1−(R−−R+) Tν (x) & exp (R− R+) log Γ(x) log x Cx + O(log x) , − − − α1 − e 1− i.e., ν . Hence case (i) is proved. If α0 = θν , then ∈PR−−R+ e Tν (x) xmin{0,1+R+−R−}( α1 x−1), T (x 1) ≥ C − ν − 28 CHUANG XU, MADS CHRISTIAN HANSEN, AND CARSTEN WIUF which yields that C C Tν (x) & exp min 0, R+ R− +1 log Γ(x) log x log x , { − } − α1 − α1 2− 1− i.e., ν if 0 > R+ R− 1, and ν if R+ R− < 1. Hence the cases ∈P1 − ≥− ∈PR−−R+−1 − − (ii) and (iii) are also proved. R> 0 (case (iv)). Based on (7.1), there exists C > 0 such that • α0 1−min{σ,R−} log Tν (x) & exp (R− R+) log Γ(x) log x Cx + O(log x) , − − − e α1 − 1− i.e., ν . Hence the former part of (iv) is proved. The seconde part of (iv) is proved ∈PR−−R+ e below in Step II. (2) R− = R+ = R. Then (7.2) is ∞ 0 ∞ ∞ ω+ (α + β (y)) ν (y j)= θ y−R− T (y)+ (α + β (y)) ν (y j) . j j − ν ν j j − y=x j=ω−+1 y=x y=x j=1 X X X X X from which it implies that there exists C > 0 and N N such that for all large x N, ∈ ∈ T (x) α Cx−σ1 (7.3) ν + − , T (x 1) ≥ α θ x−R− + Cx−σ1 ν − − − ν Based on which we establish the following lower estimates for Tν(x). (v) R> 0 and α< 0. We can show α+ exp (log α− )x + O(log x) , if min R, σ1 =1, Tν (x) & { } α+ 1−min{R,σ1} exp (log )x + O(x ) , if min R, σ1 =1, ( α− { } 6 i.e., ν 2−. The latter part is proved in Step II below. ∈P1 (vi) R> 0 and α = 0. We prove the conclusions case by case. 0
Arxiv:2007.11390V2 [Math.PR] 2 Nov 2020 Eae Tcatcmdl.Teeoe Nteecss Tis It Cases, These in Is Therefore, CTMC Models
![Arxiv:2007.11390V2 [Math.PR] 2 Nov 2020 Eae Tcatcmdl.Teeoe Nteecss Tis It Cases, These in Is Therefore, CTMC Models](http://data.docslib.org/img/c8d0bb65755e01f09a6cd0170ac94855-1.webp)