![Arxiv:1707.07350V4 [Math.PR] 11 Oct 2018 Pet.Semr Eeecso Akrud Oiain N a 16]](http://data.docslib.org/img/bcb32a6f788bb7772e083b809d260dab-1.webp)
arXiv:1707.07350v4 [math.PR] 11 Oct 2018 pet.Semr eeecso akrud oiain n a 16]. and motivations 15, 3] background, [4, Section on papers [34, survey references Liggett more also literature, See see mathematics [43], opments. In Spitzer by [35]. investigated h biology first actually in was model synthesis The protein transitions. modeling of phase in features th directio particular representative right in from including the captures removed towards models, definition, and particles simple end the its left driving despite the field to external mode an injected investigated and one-dimen particles a widely has over process also repulsion most The hardcore which the with literature. independently physics of jumping the one cles in is systems non-equilibrium dimension open one in boundaries 1.1. h tt space state the (1.1) in se–isnpoes rwinecrin Brownian excursion; process. Brownian process; Askey–Wilson tion; IPEECUINPOESSWT PNBOUNDARIES OPEN WITH PROCESSES EXCLUSION SIMPLE 2010 e od n phrases. and words Key h SPwt pnbudre sa reuil nt-tt akvp Markov finite-state irreducible an is boundaries open with ASEP The II LCUTOSFRDNIYO ASYMMETRIC OF DENSITY FOR FLUCTUATIONS LIMIT Background. ahmtc ujc Classification. Subject Mathematics Abstract. inr itiuini em fjitmmnso akvpr Markov Askey–Wi the a of of measures moments orthogonality re joint from j with of constructed the is terms process of in exclusion representation distribution simple a tionary asymmetric on depend the rely latter, proofs of the Our function of Browni time-reversal boundary). a the a processes, of p (or independent current meander two Brownian maximal of Br a the sum are of the fluctuations is boundary limit fluctuation the the at process interestingly, phases, exclusion density res Most simple low/high earlier processes asymmetric an the totally independent extends In for two This [19] of al. excursion. maximal sum et Brownian the a the In and infinity. is fl motion to the fluctuation goes of limit sites limits of scaling the number the the of as picture density complete a establish we rbto as nw stesed tt) sasohsi p stochastic a as state), t steady to the [0 respect as with known process (also exclusion tribution simple asymmetric the in , ] ntrepae ftemdladterbudre ihnt within boundaries their and model the of phases three In 1]. > α { 0 , eivsiaeteflcutoso uuaiedniyo par of density cumulative of fluctuations the investigate We 1 0 } h smercsml xlso rcs AE)wt open with (ASEP) process exclusion simple asymmetric The > β , n OZMEZBY N IA WANG YIZAO AND BRYC W LODZIMIERZ 1. ihparameters with nrdcinadmi results main and Introduction smercsml xlso rcs;saiglmt phase limit; scaling process; exclusion simple asymmetric 0 γ , ≥ 0 δ , 00;60K35. 60F05; 1 ≥ 0 , n 0 and ene;Lpaetasom tangent transform; Laplace meander; ntesm phase. same the in oesidxdby indexed rocess ≤ esainr dis- stationary he snpolynomials. lson cutoso the of uctuations pc otesta- the to spect winmotion. ownian n nteside the on ing itgenerating oint nmto and motion an ae h limit the hase, < q l yDerrida by ult cse,which ocesses, urn phase, current Brownian a , efnregion, fan he oecomplicated more piain nthe in pplications 1 . o al devel- early for .TeASEP, The n. oesparti- models inllattice, sional ticles ih end, right e sisorigin its as h model the oeson rocess transi- sfor ls 2 WLODZIMIERZ BRYC AND YIZAO WANG
Informally, the process models the evolution of the particles located at sites 1,...,n that can jump to the right with rate 1 and to the left with rate q, if the target site is unoccupied. Furthermore, particles arrive at site 1 (respectively, n), if empty, at rate α (respectively, δ), and leave site n (respectively, 1), if occupied, at rate β (respectively, γ). The transitions are summarized in Figure 1. For q < 1, particles move in an asymmetric way, with higher rate to the right than to the left; in the special case q = 0, particles move only to the right and the model is known as the totally asymmetric simple exclusion process (TASEP). We let πn denote the stationary distribution of the ASEP as a Markov process on 0, 1 n, which is also called the steady state in the physics literature. We let { } τ1,...,τn denote the occupations of each corresponding location: τj = 1 if the j-th location is occupied by a particle, and τj = 0 otherwise. All statistics of the ASEP are then expressed in terms of τ1,...,τn.
α β q 1 q 1
··· ··· 1 2 3 k 1+1 k +1k +1 n 1 +n1 ··· − ··· − γ δ Figure 1. Transition rates of the asymmetric simple exclusion process with open boundaries, with parameters α, β, γ, δ, q. Black disks represent occupied sites.
Throughout we assume (1.1) and work with the following parameterization of the ASEP, which dates back at least to the 90s in the physics literature (e.g. [40]): set
1 2 κ± = 1 q x + y (1 q x + y) +4xy , x,y 2x − − ± − − and denote p + + (1.2) A = κβ,δ, B = κβ,δ− , C = κα,γ, and D = κα,γ− . By definition, A, C 0 and it is easy to check that 1 < B,D 0, compare [10, 44]. The phase transition≥ of the ASEP is known to− be characterized≤ by A and C only. For example, it has been known since Derrida et al. [17], Sandow [40] that the ASEP has the following three phases: (1) maximal current phase A< 1,C < 1, (2) low density phase C > 1,C>A, (3) high density phase A> 1,A>C. Derrida et al. [21, 22] distinguish also the two regions (1) fan region AC < 1, (2) shock region AC > 1. Figure 2 illustrates the three phases and the two regions. In this paper, we restrict ourselves to the fan region and its boundary, as our approach does not work for the shock region. See [21, 22] for more discussions of the properties of ASEP in the shock region. LIMITFLUCTUATIONSOFOPENASEP 3
+ C = κα,γ shock region
B LD HD LD B + Bme 1 boundary of fan/shock region b
B + Bex B
MC HD
(0, 0) 1 + A = κβ,δ B + Bme
Figure 2. A phase diagram for the limit fluctuations of den- sity of the ASEP in the steady state. LD, HD, MC stand for low density, high density and maximal current (phases), respec- tively. The fan region consists of the shaded area. B, Bex and Bme stand for Brownian motion, excursion and meander, respectively, and Bme stands for the reversed Brownian meander defined by Bme d Bme Bme x x [0,1] = 1 x 1 x [0,1]. Processes in the sums are as- { } ∈ { − − } ∈ sumedb to be independent. The multiplicative constants (possibly dependingb on A, C) in front of processes are omitted.
Another set of commonly used parameters is the pair (ρa,ρb) with
1 A ρ = and ρ = . a 1+ C b 1+ A
For example, the fan region and the shock region are often characterized equiva- lently by ρa > ρb and ρa < ρb, respectively (e.g. [21] and [22, (1.1)]). From the point of view of modeling a non-equilibrium system with open boundaries, the two parameters represent the densities of the two reservoirs connected to the left and the right of the system. For convenience, we shall use A and C exclusively in the sequel. We are interested in the cumulative density function
nx 1 ⌊ ⌋ [0, 1] x τ , ∋ 7→ n j j=1 X 4 WLODZIMIERZ BRYC AND YIZAO WANG which we consider as a random process under πn. The following limits in probability for the cumulative density function are well known: 1 x A< 1,C < 1 (maximal current phase) 2 nx 1 ⌊ ⌋ 1 (1.3) lim τj = x C> 1,C>A (low density phase) n n 1+ C →∞ j=1 X A x A> 1,A>C (high density phase), 1+ A see for example [21, 42, 44]. Phase diagram affects the behavior of many other statistics, including current [20, 40], correlation functions of the density [25, 45], and the large deviation functionals of the density [21] or the current [14]. See [16] and more references therein. The fluctuations of the cumulative density function with appropriate normaliza- tion are easy to describe for the boundary of the fan region (AC = 1). Since it is known that in this case, τ1,...,τn are i.i.d. Bernoulli random variables with mean A/(1 + A)=1/(1 + C) (see [24] and Remark 2.4), the scaling limit is the Brownian motion, an immediate consequence of Donsker’s theorem [3]. Theorem 1.1 (Boundary of fan region). When AC =1,
nx ⌊ ⌋ √ 1 A A B τj x x [0,1] √n − 1+ A ⇒ 1+ A { } ∈ j=1 x [0,1] X ∈ as n in the space D([0, 1]). → ∞ One might expects naturally that the Brownian-motion behavior at the boundary of the fan region persists if one looks at the phase that is close to the boundary. However, this intuition is not entirely correct, as the limit is non-Gaussian for A < 1,C < 1 arbitrarily close to the point (A, C) = (1, 1) at the boundary, a remarkable result due to Derrida et al. [19], who showed this in the special case q =0,γ = δ = 0.
1.2. Main results. In this paper, we provide a complete picture of the limit fluc- nx tuations of the cumulative density function, that is, of the process j⌊=1⌋ τj x [0,1] { } ∈ with appropriate normalization as n , in the fan region. First, as conjectured in [19, Section 3], we show that the limit→ ∞ fluctuation for a full rangeP of parameters (1.1) in the maximal current phase is the same as for the case of q =0,γ = δ =0 studied in [19]. Second, and most interestingly, we identify two different limit fluc- tuations at the boundary of the maximal current phase. Third, in the low/high density phases in the fan region, we show that the scaling limit of fluctuations is a Brownian motion. Our results are stated in terms of Brownian motion, Brownian excursion and Brownian meander, denoted by B, Bex and Bme respectively throughout this paper. One may think of Brownian excursion as the Brownian bridge conditioned to stay strictly positive until time t = 1, and Brownian meander as the Brownian motion conditioned to stay strictly positive over time interval (0, 1]. See for example [23, 32, 37–39, 47] for more background and applications. LIMITFLUCTUATIONSOFOPENASEP 5
We first state our results on the maximal current phase and its boundary. In- troduce nx ⌊ ⌋ 1 h (x)= τ , x [0, 1], n j − 2 ∈ j=1 X and view hn(x) x [0,1] as a stochastic process with law induced by πn. The follow- ing theorem{ extends} ∈ the already mentioned result of Derrida et al. [19] to a larger range of parameters (1.1) confirming the conjecture in [19, Section 3]. We let ‘f=.d.d.’ denote convergence of finite-dimensional distributions. Recall that definition (1.2)⇒ gives A 0, C 0. ≥ ≥ Theorem 1.2 (Maximal current phase). If A< 1,C < 1 then 1 f.d.d. 1 B Bex hn(x) x [0,1] = x + x x [0,1] , √n { } ∈ ⇒ 2√2 { } ∈ as n , where the Brownian motion B and the Brownian excursion Bex are independent→ ∞ stochastic processes. The boundary of the maximal current phase, see Figure 2, splits into three regions with different limit fluctuations: the corner point where A = 1, C = 1 with asymptotically Brownian fluctuations described in Theorem 1.1, and two line- segments corresponding to A < 1, C = 1 and A = 1,C < 1 with the following fluctuations. Theorem 1.3 (Boundary of maximal current phase). We have, 1 B Bme x + x x [0,1] A =1,C < 1 2√2 { } ∈ 1 f.d.d. hn(x) x [0,1] = √n { } ∈ ⇒ 1 B Bme Bme x + 1 x 1 x [0,1] A< 1, C =1 2√2 − − ∈ B Bme as n , where the Brownian motion and the Brownian meander are independent→ ∞ stochastic processes. For the low/high density phases, we use centering as indicated in (1.3). For x [0, 1], introduce ∈ nx ⌊ ⌋ 1 hL(x)= τ n j − 1+ C j=1 X nx ⌊ ⌋ A (1.4) hH(x)= τ , n j − 1+ A j=1 X and view both as stochastic processes with laws induced by πn. Theorem 1.4 (Low/high density phases of fan region). Suppose AC < 1. In the low density phase, C > 1, we have √ 1 L f.d.d. C B hn(x) x [0,1] = x x [0,1] as n . √n ∈ ⇒ 1+ C { } ∈ → ∞ In the high density phase, A> 1, we have √ 1 H f.d.d. A B hn (x) x [0,1] = x x [0,1] as n . √n ∈ ⇒ 1+ A { } ∈ → ∞ 6 WLODZIMIERZ BRYC AND YIZAO WANG
The paper is organized as follows. In Section 1.3 we describe informally the basic ideas behind the proof. Section 2 provides technical background on Askey–Wilson processes and generating functions of ASEP. In Section 3 we prove Theorem 1.4. Section 4 presents proofs of Theorems 1.2 and 1.3. In Appendix A we discuss the Laplace transform criterion for weak convergence that we use.
1.3. Overview of the proof. Our starting point is the identity
n E n j=1(1 + tj +2√tj Ytj ) τj (1.5) tj = n n , for all 0
(n) 1/2 of Y0 converges to an infinite measure ν(du) which is proportional to either u du 1/2 (in the case A< 1,C < 1) or u− du (in the case A =1,C < 1). bUp to a normalizing constant, the Laplace transform of the finite-dimensional E (n) (n) distributions of hn(x) x [0,1] takes the form of Gn(Ys1 ,..., Ysd , Y0) for some { } ∈ d+1 s1 > >sd > 0, where the sequence of functions Gn : R R+, converges to function··· G defined in (4.23). Convergence in (1.6) is ab key step→ tob showb that these Laplace transforms converge to the limit given by the functional
E Z Z Z (G( s1 ,..., sd ,u) 0 = u)ν(du) R | Z + of the tangent process Z, see (4.31). The second ingredient of our proof consists of some recently developed duality formulas [8] that express the Laplace transforms of Brownian excursion and mean- der in terms of the tangent process Zt t 0 ((4.2) and (4.3)). We recognize that the integral above has two factors: the{ Laplace} ≥ transform of the Brownian motion, and a functional of Zs s 0 which we identify as the Laplace transform of Brownian excursion (A< {1,C} <≥ 1), see (4.32), or of Brownian meander (A = 1,C < 1), see (4.33). A delicate issue actually arises here as due to the use of Markov process we establish convergence of the Laplace transforms only in an open region away from the origin in Rd. We clarify how this leads to the desired weak convergence in Appendix A. Technical difficulties arise in the above approach when transitions probabilities of Yt t 0 are of mixed type near t = 1. We avoid this issue by applying the so- called{ } particle-hole≥ duality, which is a well known symmetry feature of the ASEP. In particular, the case A< 1, C = 1 and the low density phase will be derived from the case A =1,C < 1 and the high density phase, respectively, by this duality. Remark 1.5. In principle, our approach might work for the weakly asymmetric exclusion process, as in [18], where the authors consider the case q 1 at a rate that may depend on n and show that the fluctuations are Gaussian.↑ This would require to determine→ ∞ first the relevant tangent process as q 1. We also mention that there is a huge literature on the asymptotic↑ behavior of ASEP as a temporal-spatial process, by letting the ASEP to evolve from a non- stationary distribution, and possibly with q 1 at the rate that may depend on n . See for example [12, 26, 27, 30] and↑ references therein. Such results are beyond→ ∞ the scope of our methods.
2. Askey–Wilson process and ASEP 2.1. Askey–Wilson process. Askey–Wilson processes are a family of Markov processes based on Askey–Wilson measures, which we recall first. The Askey– Wilson measures are the probability measures that make the Askey–Wilson poly- nomials orthogonal. We do not use these polynomials here, and instead we write directly the orthogonality measure as given in [1], see also [33, Section 3.1] where a typo to weight of higher atoms is corrected. The formulas below incorporate this correction and probabilistic normalization, and come from [9]. The Askey–Wilson probability measure ν(dy; a,b,c,d,q) depends on five param- eters a,b,c,d,q. It is assumed that q ( 1, 1). For the parameters a,b,c,d, it is assumed that they are all real, or two∈ of− the parameters are real and the other 8 WLODZIMIERZ BRYC AND YIZAO WANG two form a complex conjugate pair, or the parameters form two complex conjugate pairs, and in addition (2.1) ac,ad,bc,bd,qac,qad,qbc,qbd,abcd,qabcd [1, ). 6∈ ∞ The Askey–Wilson measure is invariant with respect to permutations of a,b,c,d. More precisely, the measure is of mixed type
ν(dy; a,b,c,d,q)= f(y; a,b,c,d,q)dy + p(z)δz(dy), z F (a,b,c,d,q) ∈ X with the absolutely continuous part supported on [ 1, 1] and with the discrete part supported on a finite or empty set F . For certain choices− of parameters, the measure can be only discrete or only absolutely continuous. The absolutely continuous part is (2.2) 2 (q,ab,ac,ad,bc,bd,cd; q) (e2iθy ; q) f(y; a,b,c,d,q)= ∞ ∞ , 2π(abcd; q) 1 y2 (aeiθy ,beiθy ,ceiθy , deiθy ; q) ∞ ∞ − where y = cos θ (with the convention that f( y; a,b,c,d,q) = 0 when y > 1). y p Here and below, for complex α, n N and q < 1 we use the q-Pochhammer| | symbol ∈ ∪{∞} | | n 1 k − (2.3) (α; q) = (1 αqj ), (a , ,a ; q) = (a ; q) . n − 1 ··· k n j n j=0 j=1 Y Y The set F = F (a,b,c,d,q) of atoms of ν(dy; a,b,c,d,q) is non-empty if there is a parameter α a,b,c,d with α > 1. In this case, necessarily α is real and generates atoms:∈ for { example,} if a| >| 1 then it generates the atoms | | 1 1 (2.4) y = aqj + for j =0, 1,... such that aqj 1, j 2 aqj | |≥ and the corresponding masses are (2.5) 2 (a− ,bc,bd,cd; q) p(y ; a,b,c,d,q)= ∞ , 0 (b/a,c/a,d/a,abcd; q) ∞ 2 2 2j j (a ,ab,ac,ad; q)j (1 a q ) q p(yj ; a,b,c,d,q)= p(y0; a,b,c,d,q) − 2 , j 1. (q,qa/b,qa/c,qa/d; q)j(1 a ) abcd ≥ − The formula of p(y ; a,b,c,d,q) given here only applies for a,b,c,d = 0, and takes a j 6 different form otherwise. We shall however only need p(y0; a,b,c,d,q) in this paper. The Askey–Wilson process is a time-inhomogeneous Markov process introduced in Bryc and Wesolowski [9], based on Askey–Wilson measures. It is then explained in Bryc and Wesolowski [10] how each ASEP with parameters α,β > 0,γ,δ 0, q [0, 1) is associated to an Askey–Wilson process Y , the parameters of which≥ are denoted∈ by A,B,C,D,q, with A,B,C,D given in (1.2). As we already noted, (1.2) implies A, C 0 and 1 < B,D 0. So for the Askey–Wilson process to exist, the restriction≥ (2.1) becomes− AC <≤1, which we assume throughout in the sequel. Then, the Askey–Wilson process with parameters (A,B,C,D,q) is introduced as the Markov process with marginal distribution P(Y dy)= ν dy; A√t,B√t,C/√t,D/√t, q , 0 p 2 p p piθy for 0 n E n j=1(1 + tj +2√tj Ytj ) (2.8) tτj = , j h 2nE(1 + Y )n i *j=1 + Q 1 Y n where Yt t 0 is the Askey–Wilson process with parameters (A,B,C,D,q). { } ≥ Now, to establish our limit theorems, it suffices to analyze the asymptotics of the two expectations that appear in the numerator and in the denominator on the right-hand side of (2.8). For this purpose, we shall see that asymptotically, only the law of Yt t [1 ε,1] matters for arbitrarily small ε > 0. We first proceed in Section 3 with{ } the∈ − proof of the low/high density phases, in which case the law of Yt near upper boundary of its support is easy to analyze. Remark 2.2. The connection between Askey–Wilson polynomials and the ASEP has been known for a long time, see for example [41, 44, 45]. In [44, 45], using Askey–Wilson polynomials and complex integrals, the asymptotics of most com- monly investigated statistics are computed, including current, density, partition 10 WLODZIMIERZ BRYC AND YIZAO WANG function and the multiple-point correlation function, for results in both fan and shock regions (except the case A = C > 1 where the steady state does not have constant density). The identification of the Askey–Wilson Markov process in The- orem 2.1 turned out to be convenient for our proofs, at the expense of restriction of parameters of ASEP to the fan region AC < 1. Notice that in general, Askey– Wilson polynomials do not necessarily admit a positive orthogonality measure, and conditions on the coefficients for its existence are subtle (see [9]). Remark 2.3. The version of the matrix ansatz method that we use is more analytic so our method differs from the usual applications of the matrix ansatz that seem to have more combinatorial flavor. For example, a formula for the joint distribution of the increments of hn is given in [24, Eq. (3.7)] and used to derive the large devi- ation principle via a combinatorial argument [24, Eq. (3.16), (3.17)], essentially by expressing the probability of interest as a sum of probabilities indexed by different paths and then counting the number of paths that asymptotically have the same order of probabilities. This argument is of a completely different nature of ours. The combinatorial nature of the matrix ansatz method has also been exploited in applications to problems on combinatorial enumeration [11]. Remark 2.4. At the boundary of the fan region, AC = 1, one can read from [9, Eq. (2.14) and (2.15)] that [A + B AB(C + D)]t + C + D CD(A + B) 1 1 Yt = − − = A√t + , 2√t(1 ABCD) 2 A√t − is deterministic. Now, from (2.8) we can read out that τj j=1,...,n are independent, τj { } and t = t A/(1 + A)+1/(1 + A). So these are Bernoulli random variables j n j with A τ = . h j in 1+ A Theorem 1.1 now is a consequence of the well known Donsker’s theorem [3]. 3. Proofs for low/high density phases In this section, we investigate the case A > 1,AC < 1 and C > 1,AC < 1. In the representation (2.8), the law of the associated Askey–Wilson process with parameters (A,B,C,D; q) may have atoms. It turns out that we will only need the point mass on the largest atom. We shall only use this representation for the high density phase (A > 1,AC < 1). For the low density phase, the result shall follow by the particle-hole duality. Fix A > 1 and C < 1/A. Recall that B,D ( 1, 0] and atoms are only generated by parameters that have absolute value∈ larger− than 1, so possibly by A√t,C/√t and D/√t. When t (max 1/A2,D2 , 1], all the atoms are generated ∈ { } by A√t by (2.4) with a = A√t, and in this case we let yj (t) denote the (j + 1)-th largest atom of the law of Yt. In particular, we have 1 1 1 2 y0(t)= A√t + > 1 for t max ,D , 1 . 2 A√t ∈ A2 We shall need the mass of Y1 on y0(1), which is denoted by, recalling (2.5), (1/A2,BC,BD,CD; q) p = p(y (1); A,B,C,D,q)= ∞ . 0 0 (B/A, C/A, D/A, ABCD; q) ∞ In the sequel, we write an bn if limn an/bn = 1. ∼ →∞ LIMITFLUCTUATIONSOFOPENASEP 11 3.1. Proof of Theorem 1.4 for the high density phase A > 1,AC < 1. We prove the convergence of corresponding Laplace transform. We first recall the Laplace transform of the finite-dimensional distribution of the Brownian motion. For x = 0 < x < < x x = 1, c ,...,c > 0, s = c + + c , 0 1 ··· d ≤ d+1 1 d k k ··· d k =1,...,d, and sd+1 = 0, we have d d E B E B (3.1) exp ck xk = exp (sk sk+1) xk − ! − − ! Xk=1 Xk=1 d d+1 E B B 1 2 = exp sk( xk xk 1 ) = exp sk(xk xk 1) . − − − ! 2 − − ! Xk=1 Xk=1 For the ASEP in the high density phase, consider the centered cumulative density c Rd function (1.4) and its Laplace transform with argument = (c1,...,cd) + defined by ∈ d H c H ϕx,n( )= exp ckhn (xk) . * − !+ Xk=1 n Note that in Theorem 1.4, the limit Brownian motion is scaled by √A/(1 + A). Therefore, by Theorem A.1 to prove the high density phase of Theorem 1.4 it suffices to prove c d+1 H A 2 lim ϕx,n = exp sk(xk xk 1) . n 2 − →∞ √n 2(1 + A) − ! kX=1 To do so, we first write d nxk ⌊ ⌋ H c A ϕx,n( )= exp τj (ck + + cd) * − − 1+ A ··· + k=1 j= nxk 1 +1 X ⌊ X− ⌋ n d d+1 nk A sk τj = exp sk(nk nk 1) (e− ) , 1+ A − − k=1 !*k=1 j=nk 1+1 + X Y Y− n where n = nx , k =1,...,d + 1. By (2.8), k ⌊ k⌋ d+1 nk nk 1 sk sk/2 − − 1 1+ e− +2e− Y sk H c E e− (3.2) ϕx,n( )= n s A/(1+A) 2 Zn " e− k # kY=1 n with Zn = E(1 + Y1) . Lemma 3.1. If A> 1 and AC < 1, then (1 + A)2n Z p . n ∼ 2nAn 0 This result has been known in the literature. See Remark 4.6. We provide a proof here for completeness. Proof. Let y1∗(1) = max(y1(1), 1) denote the upper bound of the support of the law of Y on R y (1) . Note that y∗(1) < y (1) (recall that y (1) and y (1) are the 1 \{ 0 } 1 0 0 1 12 WLODZIMIERZ BRYC AND YIZAO WANG locations of the point masses of the first and second atom of Y1, see (2.4), and that the first atom lies above 1 and above the second atom). It then follows that n n Zn = (1 + y) ν(dy; A,B,C,D,q)+ (1 + y) ν(dy; A,B,C,D,q). y0(1) [ 1,y (1)] Z{ } Z − 1∗ The first term equals (1 + A)2n p (1 + y (1))n = p , 0 0 0 2nAn n and the second term is bounded from above by (1 + y1∗(1)) , which converges to 0 n when divided by (1 + y0(1)) . In view of the asymptotics of Zn, we introduce s s/2 1+ e− +2e− y A ψ(s,y)= sA/(1+A) 2 , e− (1 + A) and have c d+1 nk nk 1 H Mn sk − − ϕx with M = E ψ , Y s /√n . ,n n e− k √n ∼ p0 " √n # kY=1 The desired result now follows from the following. Lemma 3.2. With the notation above, d A 2 lim Mn = p0 exp sk(xk xk 1) . n 2 − →∞ 2(1 + A) − ! Xk=1 Proof. We have that (3.3) 2 s 1 sA/(1+A) A s/(1+A) As 2 ψ(s,y (e− )) = e + e− =1+ + o(s ) 0 1+ A 1+ A 2(1 + A)2 s /√n as s 0. Introduce sk,n = sk/√n and tk,n = e− k . Note that due to our choice of s ↓, we have t d+1 nk nk 1 E − − Mn,1 = ψ sk,n, Ytk,n 1 . Yt1,n =y0(t1,n) "k=1 { }# Y We shall show that Mn Mn,1 as n . Indeed, we have that Ytk,n y0(tk,n) and hence ψ(s , Y ) ∼ ψ(s ,y (t→ ∞)) almost surely. First, observe≤ that k,n tk,n ≤ k,n 0 k,n P P (Yt1,n = y0(t1,n)) = (Ytk,n = y0(tk,n), k =1,...,d + 1). That is, once the process Ys reaches the highest point y0(s) at some time s, nec- 2 essarily s> 1/A , Ys stays on the deterministic trajectory (y0(t))t s. This follows by computing P(Y = y (t) Y = y (s)) for 1/A2 So by (2.6) and (2.5), P(Y = y (t) Y = y (s)) = ν y (t) ; A√t,B√t, s/tA√s, s/t/(A√s) t 0 | s 0 { 0 } p p = p y0(t); A√t,B√t,As/√t, 1/(A√t), q (1/(A2t),ABs,B/A,s/t; q) = ∞ =1. (B/A, s/t, 1/(A2t),ABs; q) ∞ P Introduce also p0,n = (Yt1,n = y0(t1,n)). Recalling (2.5), we have p0,n = p y0(t1,n); A t1,n,B t1,n,C/ t1,n,D/ t1,n, q 2 (1/(A t1,n),BC,BD,CD/t1,n; q) = p p p ∞ p (B/A, C/(At1,n),D/(At1,n),ABCD; q) ∞ (1/A2,BC,BD,CD; q) ∞ = p → (B/A, C/A, D/A, ABCD; q) 0 ∞ as n . Therefore, by (3.3), → ∞ d+1 d nk nk 1 A 2 p − − p Mn,1 = 0,n ψ (sk,n,y0(tk,n)) 0 exp 2 sk(xk xk 1) → 2(1 + A) − − ! kY=1 kX=1 as n . On the other hand, introducing y∗(s) = max y (s), 1 , on the event → ∞ 1 { 1 } Y = y (s) , we have Y y∗(s) almost surely. Then, { s 6 0 } s ≤ 1 d+1 n1 nk nk 1 M (1 p )ψ(s ,y∗(t )) ψ (s ,y (t )) − − . n,2 ≤ − 0,n 1,n 1 1,n k,n 0 k,n kY=2 Now, one sees immediately that, since by continuity limn y1∗(t1,n)/y0(t1,n) = →∞ y∗(1)/y (1) (0, 1), 1 0 ∈ n1 M 1 p ψ(s ,y∗(t )) n,2 − 0,n 1,n 1 1,n M ≤ p ψ(s ,y (t )) n 0,n 1,n 0 1,n 2s1 s1 n1 1 p 1+ e− ,n +2e− ,n y∗(t ) = − 0,n 1 1,n 0 p 1+ e 2s1,n +2e s1,n y (t ) → 0,n − − 0 1,n as n . Therefore M M , and the desired result follows. → ∞ n ∼ n,1 3.2. Proof of Theorem 1.4 for the low density phase C > 1,AC < 1. The result for the low density phase is an immediate consequence of the result for the high density phase, by the particle-hole duality which we now explain. We have seen the definition of an ASEP with parameters (α, β, γ, δ, q). Instead of thinking of particles jumping around, we view the particles as background and allow the holes to jump around (viewing Figure 1 as white particles jumping around among black unoccupied sites). In this way, equivalently a hole jumps to the unoccupied left and right sites with rates 1 and q, respectively, and disappears at site 1 with rate α and at site n with rate δ, and enters site n if unoccupied with rate β and site 1 if unoccupied with rate γ. This is the ASEP with parameters (β, α, δ, γ, q), if we relabel the sites 1,...,n by n, . . . , 1 . Fix q [0, 1). Let{ πA,B,C,D} denote{ the} stationary distribution of the ASEP ∈ n with parameters (α, β, γ, δ, q). Let τ1,...,τn be as before, and set εj =1 τn j+1. − − 14 WLODZIMIERZ BRYC AND YIZAO WANG Introduce nx ⌊ ⌋ C hL(x)= ε . n j − 1+ C j=1 X b L A,B,C,D The above argument shows that hn(x) x [0,1] with respect to πn has the H { } ∈ C,D,A,B same law as hn (x) x [0,1] (defined in (1.4)) with respect to πn . Therefore, the high density{ phase} ∈ of Theoremb 1.4 tells that √ 1 L f.d.d. C B hn(x) = x x [0,1] . √n x [0,1] ⇒ 1+ C { } ∈ n o ∈ L nx A,B,C,D We are interested in h (bx) = ⌊ ⌋(τ 1/(1 + C)) with respect to π . n j=1 j − n Observe that C P 1 εj + τn j+1 =0, − 1+ C − − 1+ C so 0 nx = nx ⌊ ⌋ (3.4) hL(x)+ hL(1) hL(1 x) = . n n − n − C ε n(1 x) +1 nx = nx 1+ C ⌊ − ⌋ b b − 6 ⌊ ⌋ Since the error term is uniformly bounded, the finite-dimensional distributions of 1/2 L n− hn(x) x [0,1] have the same limit as the finite-dimensional distributions of ∈ 1/2{ L } L n− hn(1 x) hn(1) x [0,1], and we arrive at { − − } ∈ 1 √C √C b L b f.d.d. B B d B hn(x) x [0,1] = 1 x 1 x [0,1] = x x [0,1] . √n ∈ ⇒ 1+ C { − − } ∈ 1+ C { } ∈ This proves the low density phase of Theorem 1.4. 4. Proofs for Maximal current phase and its boundary The proofs for the two cases are very similar and are hence unified. We need some preparation for the proof. In Section 4.1 we review an important auxiliary Markov process Z, and in particular how this Markov process shows up in the Laplace representations of Brownian excursion and meander. Another important role of this Markov process is that it is the tangent process of the Askey–Wilson process at the boundary. This result, playing a central role in the proof, will be established first in Proposition 4.1 in Section 4.2. The case A 1,C < 1 is then proved in Theorem 4.4 in Section 4.3. The case A < 1, C = 1≤ is proved by the particle-hole duality in Section 4.4. 4.1. An auxiliary Markov process. An auxiliary Markov process, denoted by Z in the rest of the paper, will play a crucial role in the proof for the maximal current phase and its boundary. This is a positive self-similar Markov process with values in [0, ) and transition probability density function ∞ 2(t s)√y (4.1) gs,t(x, y)= − 1 x 0,y 0 , s < t. π [(t s)4 + 2(t s)2(x + y) + (x y)2] { ≥ ≥ } − − − This process is self-similar in the sense that, letting Px denote the law of Z starting at Z0 = x, Z P d 2 Z P λt t 0 , x = λ t t 0 , x/λ2 for all λ,x > 0. { } ≥ { } ≥ LIMITFLUCTUATIONSOFOPENASEP 15 This process has not been much investigated in the literature, except for a series of recent papers [7, 8, 46]. In [7], when investigating the path properties of so- called q-Gaussian processes, we proved that the process Z arises as their tangent process at the boundary. (We also proved in Section 3 therein that the transformed 2 process Z via Zt = Zt/2 + t /4 has already shown up in the literature: in a general framework connecting non-commutative stochastic process and classical Markov process developede e by Biane [2]. In this framework, Z as a classical Markov process corresponds to the free 1/2-stable process, the knowledge of which we do not need here.) e Recently, we also found out in [8] that the process Z plays an intriguing role in representations of the Laplace transforms of finite-dimensional distributions of Brownian excursion and Brownian meander. Write E ( ) = E( Z = u). For all u · · | 0 d N, s1 >s2 > >sd >sd+1 =0, and 0= x0 < x1 < < xd xd+1 = 1, we have∈ shown in [8]··· that ··· ≤ d E Bex (4.2) exp (sk sk+1) xk − − ! Xk=1 d+1 1 1/2 E 1 Z = u du u exp sk (xk xk 1) , R − √2π + −2 − ! Z kX=1 and d E Bme (4.3) exp (sk sk+1) xk − − ! Xk=1 d+1 1 1 E 1 Z = du u exp sk (xk xk 1) . R 1/2 − √2π + u −2 − ! Z kX=1 Note that the formulae here are obtained via the changes of variables sk = sd+1 k and xk =1 td+1 k, k =1,...,d +1, − − − and x = 0, where s and t correspond to variables s ,t used in [8]. 0 b k k k k The above identities can be obtained by direct computation using the joint den- sity functions of Brownianb excursion and meander. These explicit densities of the two processes will not be used in this paper. Standard references about Brownian excursions and meanders include [37–39]. 4.2. Tangent process of Askey–Wilson process with A 1 and C < 1. The ≤ proof is essentially based on the laws of the Askey–Wilson process Yt t [1 ε,1] for some ε > 0 small enough. With A 1 and C < 1, for this range{ } of∈ t−the marginal and transition probability laws≤ are absolutely continuous with respect to the Lebesgue measure, with compact support on [ 1, 1]. For the purpose of computing asymptotics, we express the formula by regrouping− factors into those that tend to a non-zero constant as t 1 and x 1, and those that tend to zero. (Some factors go to zero only when A↑ = 1, and↑ we include them in the second group.) In particular, the Askey–Wilson process Y has the marginal probability 16 WLODZIMIERZ BRYC AND YIZAO WANG density function (4.4) πt(x)= f x; A√t,B√t,C/√t,D/√t, q (q; q) (ABt,AC,AD,BC,BD,CD/t ; q) = ∞ ∞ 2π(ABCD; q) (B√teiθx ,Ceiθx /√t,Deiθx /√t; q) 2 ∞| ∞| (e2iθx ; q) 2 | ∞| , × √1 x2 (A√teiθx ; q) 2 − | ∞| with x = cos θ , for x [ 1, 1], and transition probability density function x ∈ − iθx iθx (4.5) ps,t(x, y)= f y; A√t,B√t, s/te , s/te− (q; q) (ABt; q) p(B√seiθxp; q) 2 = ∞ ∞| ∞| 2π(ABs; q) (B√teiθy ; q) 2 ∞| ∞| (A√seiθx ,e2iθy ; q) 2(s/t; q) | ∞| ∞ , 2 iθ i(θ +θ ) i( θ +θ ) 2 × 1 y (A√te y , s/te x y , s/te − x y ; q) − | ∞| where y = cos θy, for x,p y [ 1, 1]. (Recallp (2.2) and qp-Pochhammer notation in (2.3).) ∈ − It turns out to be helpful to consider Yt t [1 ε,1] in the reversed time direction { } ∈ − when Y1 is close to 1. For this purpose, we introduce Markov process (n) (4.6) Y =2n (1 Y 2s/√n ) ,s 0. s − e− ≥ (The parameterization for Y (n) is chosen so that we have the corresponding con- b s vergence to the tangent process (4.7) below.) The following limit theorem of Y (n) (n) is at the core of the proof ofb Theorems 1.2 and 1.3. Let πs denote the probability (n) (n) (n) b density function of Ys and ps,t the transition probability density of Y from time s to t. b b b Proposition 4.1. Under theb notation above, for A 1 and C < 1, 0 s u v 2s/√n 2t/√n x =1 ,y =1 ,s = e− ,t = e− , n − 2n n − 2n n n 2c 1 √u A< 1,C < 1 √n · (4.10) π (x ) sn n ∼ √u 2c √n A =1,C < 1, 2 · s2 + u and there exists some constant c such that u c A< 1,C < 1 n r (4.11) π (x ) sn n ≤ n c A =1,C < 1, u r for all n large enough. Moreover, 2n g (v,u) A< 1,C < 1 · s,t (4.12) p (y , x ) tn,sn n n ∼ t2 + v 2n g (v,u) A =1,C < 1. · s,t s2 + u Proof. We first establish 2 iθx iθy 1 2 2 2 (4.13) tn/sne± n e n ; q (t s) + (√u √v) (q; q) . ∼ n − ± ∞ ∞ p 18 WLODZIMIERZ BRYC AND YIZAO WANG For this, write 2 iθx iθy (4.14) tn/sne± n e n ; q ∞ 2 2 p i( θ +θ ) i( θ +θ ) = 1 t /s e ± xn yn t /s e ± xn yn q; q . − n n n n ∞ Since t 1,s 1,θ p0,θ 0, the second p factor above is asymptotically n n xn yn equivalent→ to (q; q→)2 . For→ the first factor,→ ∞ (4.15) 2 i( θx +θy ) 2 2 1 t /s e ± n n =1+ t /s 2 t /s x y 1 x 1 y − n n n n − n n n n ∓ − n − n p 1 p p p (t s)2 + (√u √v)2 . ∼ n − ± This proves (4.13). As special cases, we have 2 iθx 1 2 2 (4.16) √sne n ; q s + u (q; q) , ∞ ∼ n ∞ 2(t s) (tn/sn; q) − (q; q) , ∞ ∼ √n ∞ 2 u eiθxn ; q (q; q)2 , ∞ ∼ n ∞ 2 4u e2iθxn ; q (q; q)2 , ∞ ∼ n ∞ as n . → ∞ We now examine πsn (xn) using (4.4). The pointwise asymptotics (4.10) are straightforward to obtain. We prove the upper bounds. Recall that πsn has support on [ 1, 1]. Therefore from now on we assume x [ 1, 1], or equivalently u − n ∈ − ∈ [0, 4n]. We first focus on the second fraction of πsn (xn) in the expression (4.4), 2 (e2iθxn ; q) ∞ πsn (xn)= 2 . 2 iθx 1 x n (A√sne n ; q) − ∞ Let c denote a constante independentp of n, but may change from line to line. Suppose 2 iθx 2 A < 1 first. Observe that (A√sne n ; q) (A; q) . Furthermore, we have ∞ ≥ ∞ 2iθxn 2 2 2iθxn 2 (e ; q) = 4(1 xn) (e q; q) (see e.g. (4.14) and (4.15)). Therefore, | ∞| − | ∞| 2 (e2iθxn ; q) 2 u πs (xn) c ∞ c 1 x c . n ≤ 2 ≤ − n ≤ n 1 xn r − p Now suppose A =e 1. Then p 2 (e2iθxn ; q) 1 x2 ∞ n πsn (xn) c 2 c − 2 . ≤ 2 iθxn ≤ iθxn 1 xn (√sne ; q) 1 p√sne − ∞ − e iθ 2 2 By (4.15), 1 √s e xpn = 1+ s 2√s x = (1 √s ) +2√ s (1 x ) > − n n − n n − n n − n 2√s1(1 xn). So, we see that πs (xn) is bounded from above by, − n 1 x2 n ec − n c . 1 x ≤ u p − n r LIMITFLUCTUATIONSOFOPENASEP 19 So we see that πsn (xn) can be controlled by the bounds in (4.11). For the first fraction of πsn (xn) in (4.4), it suffices to control e (CD/sn; q) ∞ iθ iθ . (Ce xn /√sn,De xn /√sn; q) ∞ For the numerator, since CD ( 1, 0], we have (CD/sn; q) (CD/s1; q) . ∈ − ∞ ≤ ∞ For the denominator, for n large enough so that C/√sn < (1 + C)/2 < 1, we 2 iθx 2 iθx 2 have (Ce n /√sn; q) ((1 + C)/2; q) . Similarly, (De n /√sn; q) > ∞ ≥ ∞ | ∞| ((1 D)/2; q) for n large enough so that D/√sn < (1 D)/2 < 1. So the first − ∞ − − fraction can be bounded by some constant c. This completes the proof of (4.11). Now to show (4.12), observe first that by (4.13) applied to each factor, 2 i(θx +θy ) i( θx +θy ) tn/sne n n ; q tn/sne − n n ; q ∞ ∞ p1 p (q; q)4 (t s)2 + (√u + √v)2 (t s)2 + (√u √v)2 ∼ n2 ∞ − − − 1 = (q; q)4 (t s)4 + 2(t s)2(u + v) + (u v)2 . n2 ∞ − − − So (4.5) now yields 2 iθy (q; q) (A√tne n ; q) ∞ ∞ ptn,sn (yn, xn) 2 ∼ 2π iθxn (A√sne ; q) ∞ 2 2iθx (e n ; q) (tn/sn; q) ∞ ∞ 2 × 2 i(θy +θx ) i( θy +θx ) 1 xn ( tn/sn e n n ; q) ( tn/sne − n n ; q) − ∞ ∞ 2 iθy (q;pq) (A√ ptne n ; q) p ∞ ∞ 2 ∼ 2π iθxn (A√sne ; q) ∞ 4u 2 2(t s) (q; q) − (q; q) n ∞ √n ∞ × u 1 (q; q)4 (t s)4 + 2(t s)2(u + v) + (u v)2 n n2 ∞ − − − r 2 iθy (A√tne n ; q) ∞ = 2 2n gs,t(v,u). iθxn · · (A√sne ; q) ∞ By (4.16), we have 1 A< 1 2 iθy (A√tne n ; q) lim ∞ = n iθ 2 2 →∞ (A√sne xn ; q) t + v ∞ A =1. s2 + u This completes the proof. Proof of Proposition 4.1. Recall the relation between Y and Y in (4.6). We have (n) that πs (u) = πsn (xn)/(2n), and so the second part of the proposition follows from (4.10) and (4.11). Fix 0 b 20 WLODZIMIERZ BRYC AND YIZAO WANG (n) 1 πtn (yn) 1 ps,t (u, v)= psn,tn (xn,yn) = ptn,sn (yn, xn) , 2n πsn (xn) 2n where psn,tn denotes the transition probability of Y in the reversed time direction. In the case A 0, which completes the proof for the first part of the proposition. 4.3. Proof of Theorems 1.2 and 1.3 for the case A 1,C < 1. For d N, ≤ ∈ c1,...,cd > 0, x0 = 0 < x1 < < xd xd+1 = 1, introduce the Laplace transform ··· ≤ d ϕx,n(c)= exp ckhn(xk) * − !+ Xk=1 n d nxk ⌊ ⌋ 1 = exp τj (ck + + cd) . * − − 2 ··· + k=1 j= nxk 1 +1 X ⌊ X− ⌋ n This time it will be more convenient to write 1 (4.17) s = (c + + c ), k =1,...,d k 2 k ··· d (notice the extra 1/2 compared to sk in previous sections), sd+1 = 0, and nk = nx , k =1,...,d + 1. We have, by Theorem 2.1 again, ⌊ k⌋ d d+1 nk 2sk τj ϕx,n(c) = exp sk(nk nk 1) (e− ) − − k=1 !*k=1 j=nk 1+1 + X Y Y− n d+1 2sk sk nk nk 1 d E 2s − k=1(1 + e− +2e− Ye− k ) − = exp s (n n ) k k k 1 nE n − − ! hQ 2 (1 + Y1) i Xk=1 d+1 nk nk 1 E 2s − − k=1 (cosh(sk)+ Ye− k ) = , hQ Zn i n where we write Zn = E(1 + Y1) . In order to establish the convergence of finite- dimensional distributions of hn, we compute the limit of d+1 nk nk 1 E cosh(s /√n)+ Y 2s /√n − − c k=1 k e− k (4.18) ϕx = , ,n √n h Z i Q n as n , and identify the limit with the Laplace transform of the corresponding process.→ ∞ Now in view of Theorem A.1, Theorems 1.2 and 1.3 are a consequence of the following. LIMITFLUCTUATIONSOFOPENASEP 21 Theorem 4.4. For all d N, 0 < x < < x 1 and c = (c ,...,c ) (0, )d, ∈ 1 ··· d ≤ 1 d ∈ ∞ c (4.19) lim ϕx,n n √n →∞ d d ck ck ex E exp Bx E exp B A< 1,C < 1, − √ k − √ xk k=1 2 2 ! k=1 2 2 ! = X X d d ck ck me E exp Bx E exp B A =1,C < 1. − 2√2 k − 2√2 xk k=1 ! k=1 ! X X n We first determine the asymptotics of E(1 + Y1) in (4.18). Lemma 4.5. We have 2n 4√π c A< 1,C < 1 n3/2 1 n · · Zn = E(1 + Y1) ∼ 2n 2√π c A =1,C < 1. n1/2 · · 2 where c1, c2 are defined in Proposition 4.1. n Remark 4.6. The quantity Zn = E(1 + Y1) is closely related to the partition function in the literature, denoted by n for the discussion here, via (see explanation in [10, Remark 5]), Z 2n W V = Z h | i, Zn n (1 q)n − where the factor W V depends on the choice of the vectors in (2.7) and needs to be taken intoh account| i when comparing Lemma 4.5 with the literature. Par- tial results on asymptotics of partition function, including also low/high den- sity phases, have been known. See for example [20, (52),(53) and (55)] for the case q = 0,γ = δ = 0 (with W V = 1) and [5, (56)] for A < 1,C < 1,B = D = 0 with W V =h 1/|(ACi ; q) . In more generality, Uchiyama h | i ∞ et al. [44, (6.6) and (6.9)] compute n for A > 1,A > C and A,C < 1 (with W V = (ABCD; q) /(q,AB,AC,AD,BC,BD,CDZ ; q) ). We do not find gen- h | i ∞ ∞ eral results on asymptotic of Zn for A =1,C < 1 in the literature. Proof. It follows from (4.10) that, taking s = 0 therein, as n , for u> 0, → ∞ u 2c A< 1,C < 1 1 n u r (4.20) π 1 1 − 2n ∼ n 2c A =1,C < 1. 2 u r Consider first the case A = 1. Then 1 n 1/2 n n 1/2 n (4.21) 2− n E(1 + Y1) =2− n (1 + y) π1(y)dy 1 Z− u n u du 1 = u 4n 1 π1 1 1/2 . R+ { ≤ } − 4n − 2n 2n Z 22 WLODZIMIERZ BRYC AND YIZAO WANG In view of (4.20), the integrand on the right-hand side of (4.21) converges to u/4 ∞ u/4 c2e− /√u. Since 0 e− /√udu = 2√π, to conclude the proof by the domi- nated convergence theorem, we now give an integrable bound for the integrands. R Recalling (4.11), and noting that π1 (1 u/(2n)) = 0 for u > 2n, we have, for n large enough, − n u u u/4 n 1 π 1 ce− , for all u> 0 − 4n 1 − 2n ≤ u r for some constant c that depends on C and q< 1. The proof for A< 1,C < 1 is similar, starting from (4.21) and eventually by the dominated convergence theorem. In this process, the above bound is replaced by, for n large enough, n u u u/4 u 1 π 1 ce− , for all u> 0. − 4n 1 − 2n ≤ n r Details are omitted. Next, we look at the numerator in (4.18). We write, recalling Y in (4.6), d+1 nk nk 1 n sk − − b(n) (n) E 2s /√n E (4.22) 2− cosh + Ye k = Gn Ys1 ,..., Ysd+1 " √n − # kY=1 b b with d+1 nk nk 1 n sk uk − − Gn(u)=2− cosh +1 1 u 4n √n − 2n { ≤ } kY=1 d+1 nk nk 1 2 sk uk − − = 1 + sinh 1 u 4n . 2√n − 4n { ≤ } kY=1 Recall that Y takes values from [0, 4n]. We introduce the indicator function above for the convenience in later analysis (u 4n stands for max u 4n ≤ k=1,...,d+1 k ≤ here). Sinceb (nk nk 1)/n xk xk 1, we have − − → − − d+1 u 1 2 u (4.23) lim Gn( ) = exp (sk uk)(xk xk 1) =: G( ). n − →∞ 4 − − ! Xk=1 The key step is to show the following. Recall that s > >s >s = 0. 1 ··· d d+1 Proposition 4.7. With the notations above and E ( )= E( Z = u), u · · | 0 E (n) (n) (n) Gn Ys1 ,..., Ysd , Y0 c 1 E Z Z 1/2 b b b uG( s1 ,..., sd ,u)u du A< 1,C < 1 n3/2 R Z + ∼ c 1 2 E Z Z uG( s1 ,..., sd ,u) du A =1,C < 1. n1/2 R u1/2 Z + LIMITFLUCTUATIONSOFOPENASEP 23 m Proof. We start with some properties of Gn as preparations. First, since (1+x) exp(mx) for m N, x 1, we get an exponential bound on G : ≤ ∈ ≥− n d+1 nk nk 1 2 sk Gn(u) exp − − 4n sinh uk 1 u 4n ≤ 4n 2√n − { ≤ } kY=1 d+1 nk nk 1 (4.24) c exp − − u . ≤ − 4n k kY=1 Here and below, c denotes a constant that does not depend on u and n, but may vary from line to line. This inequality also shows that Gn(u) is uniformly bounded u Rd+1 N for all + and n . ∈ ∈ n n n Next, by the inequality a b M n a b provided that | k=1 k − k=1 k|≤ k=1 | k − k| ak , bk M, we have, for all u, u′ [0, 4n], | | | |≤ Q ∈Q P n d+1 2 s1 uk uk′ G (u) G (u′) 1 + sinh (n n )| − | | n − n |≤ 2√n k − k+1 4n kX=1 d+1 c u u′ . ≤ | k − k| Xk=1 It follows that for all u , u Rd+1, n ∈ + u u u u Rd+1 (4.25) lim Gn( n)= G( ), if n + as n . n →∞ → ∈ → ∞ In the remaining part of the proof with a little abuse of notation we write E ( ) u · E (n) (n) which is to be understood as ( Y0 = u) when dealing with Y and as E Z Z· | (n) (n) ( 0 = u) when dealing with . Conditioning on the value of Ysd+1 = Y0 , we write· | b b b b E (n) (n) ∞ E (n) (n) (n) (4.26) Gn Ys1 ,..., Ysd+1 = uGn Ys1 ,..., Ysd ,u π0 (u)du. 0 Z Because of the uniformb boundedness,b (4.25) andb the weakb convergenceb of the tan- gent process established in Proposition 4.1, we have E (n) (n) E Z Z (4.27) lim uGn Y ,..., Y ,u = uG( s1 ,..., s ,u) for all u> 0 n s1 sd d →∞ (see e.g., [3, Exercise 6.6]).b b By (4.8) the integrands on the right hand side of (4.26) converge pointwise under appropriate normalization. That is, when A< 1,C < 1 we have 3/2 n E (n) (n) (n) E Z Z 1/2 (4.28) lim uGn Ys1 ,..., Ysd ,u π0 (u)= uG( s1 ,..., sd ,u)u n c1 →∞ and when A =1,C < 1 we have b b b n1/2 1 (4.29) lim E G Y (n),..., Y (n),u π(n)(u)= E G(Z ,..., Z ,u) . u n s1 sd 0 u s1 sd 1/2 n c2 u →∞ To conclude the proof we now apply the dominated convergence theorem. We see b b b that, if xd < 1, (4.24) yields Gn(u) c exp( (1 xd)ud+1/4). Recall upper bounds (n) ≤ − − on π0 in (4.9). Therefore, the functions of u that appear on the left-hand side of (4.28) and (4.29) are bounded by the integrable functions cu1/2 exp( (1 x )u/4) − − d b 24 WLODZIMIERZ BRYC AND YIZAO WANG 1/2 and cu− exp( (1 x )u/4), respectively, for n large enough. This proves the − − d case xd < 1. To prove the case xd = 1, we need to work a little harder, although the approach is very similar. Notice that when xd = 1, both Gn(u) and its limit G(u) do not depend on the last coordinate ud+1. Therefore we introduce Gn∗ (u1,...,ud)= Gn (u1,...,ud, 0) and G∗ (u1,...,ud)= G (u1,...,ud, 0) (the choice of 0 in the last variable of Gn and G is irrelevant for the definitions), (n) and write E∗ ( ) as either E( Ys = u) or E( Z = u), depending on whether u · · | d · | sd (n) Z (n) the conditional expectation is for Y or . By conditioning on the value of Ysd , we write b b b (n) (n) (n) (n) EGn Y ,..., Y EG∗ Y ,..., Y s1 sd+1 ≡ n s1 sd ∞ E (n) (n) (n) b b = u∗bGn∗ Ys1 b,..., Ysd 1 ,u πsd (u)du. 0 − Z The same argument in the proof of Propositionb 4.1 forb the convergeb nce in law of tangent process leads to (n) (n) (n) Z Z Z Ys1 ,..., Ysd 1 Ysd = u ( s1 ,..., sd 1 ) sd = u . L − ⇒ L − | Therefore, byb the sameb argument b for (4.27), we have E (n) (n) E Z Z lim u∗ Gn∗ Ys1 ,..., Ysd 1 ,u = u∗ G∗ s1 ,..., sd 1 ,u , n − →∞ − and instead of (4.28) andb (4.29), web have, for A< 1,C < 1, 3/2 n E (n) (n) (n) E Z Z 1/2 lim u∗ Gn∗ Ys1 ,..., Ysd 1 ,u πsd (u)= u∗ G∗( s1 ,..., sd 1 ,u)u n c1 − − →∞ and for A =1,C < 1,b b b 1/2 1/2 n E (n) (n) (n) E Z Z u lim u∗ Gn∗ Ys1 ,..., Ysd 1 ,u πsd (u)= u∗ G∗( s1 ,..., sd 1 ,u) . n − 2 c2 − u + sd →∞ b b b (n) (n) One can show that for n large enough, πsd has the same upper bounds as π0 (n) in (4.9) (independent of s ), by applying (4.11) to πs (u) = π (1 u/(2n))/2n. d d sd − By (4.24), for n large enough, Gn(u) bc exp( (xd xd 1)ud/8) (we cannotb use ≤ − − − (xd xd 1)/4 as before, because of the rounding issueb caused by nk = nxk ). So the− dominated− convergence theorem applies, and we arrive at ⌊ ⌋ E (n) (n) Gn∗ Ys1 ,..., Ysd c 1 E Z Z 1/2 b b u∗ G∗( s1 ,..., sd 1 ,u)u du A< 1,C < 1 n3/2 R − Z + ∼ 1/2 c2 E Z Z u u∗ G∗( s1 ,..., sd 1 ,u) 2 du A =1,C < 1. n1/2 R − u + s Z + d LIMITFLUCTUATIONSOFOPENASEP 25 These are however not the same expressions as desired yet, and we need to rewrite them. For the case A < 1,C < 1, since Z is stationary with respect to the distri- bution u1/2du, we have (recalling that E ( )= E( Z = u)) u · · | sd+1 E Z Z 1/2 E Z Z Z 1/2 u∗ G∗( s1 ,..., sd 1 ,u)u du = uG∗( s1 ,..., sd 1 , sd )u du. R − R − Z + Z + For the case A =1,C < 1, because of the fact that ∞ 1 √y (4.30) g (x, y)dx = , x1/2 0,t y + t2 Z0 (recall (4.1)) we have 1/2 E Z Z u E Z Z Z 1 u∗ G∗( s1 ,..., sd 1 ,u) 2 du = uG∗( s1 ,..., sd 1 , sd ) du. R − u + s R − u1/2 Z + d Z + So to complete the proof it remains to show (4.30). For this purpose, introduce change of variables x = u2,y = v2, and we get elementary integrals ∞ 1 ∞ g (x, y)dx =2 g (u2, v2)du x1/2 0,t 0,t Z0 Z0 ∞ 1 t t = 2 2 2 2 du 0 πu t + (u v) − t + (u + v) Z − u= t t2 + (u + v)2 ∞ = log 2π (t2 + v2) t2 + (u v)2 u=0 − u= v u v u + v ∞ + arctan − + arctan π (t2 + v2) t t u=0 v √y = = . t2 + v2 y + t2 Proof of Theorem 4.4. By Lemma 4.5 and Proposition 4.7, we have (4.31) 1 E Z Z 1/2 uG( s1 ,..., sd ,u)u du A< 1,C < 1 4√π R c Z + lim ϕx,n = n √ →∞ n 1 du E Z Z uG( s1 ,..., sd ,u) 1/2 A =1,C < 1. 2√π R+ u Z We first prove the case of A< 1,C < 1. Observe that E Z Z 1/2 (4.32) uG( s1 ,..., sd ,u)u du R Z + d+1 d+1 1 2 E 1 Z 1/2 = exp sk(xk xk 1) uexp sk (xk xk 1) u du. − R − 4 − ! × + −4 − ! Xk=1 Z Xk=1 The first exponential function on the right-hand side above corresponds to the 1/2 Laplace transform of the scaled Brownian motion 8− B (see (3.1) and recall that in (3.1) we used sk = ck + +cd, but we have been using sk = (ck + +cd)/2 since (4.17)). The integral on the··· right-hand side of (4.32) corresponds··· to the Laplace 26 WLODZIMIERZ BRYC AND YIZAO WANG transform of a Brownian excursion. Indeed, by self-similarity of the process Z, it equals d+1 1/2 E 1 1Z u du u exp sk (xk xk 1) R − + −2 2 − ! Z Xk=1 d+1 1/2 E 1 Z = u du u/2 exp s /√2(xk xk 1) R k − + −2 − ! Z kX=1 d+1 1 √ 1/2 E Z = 8 u du u exp s /√2(xk xk 1) . R k − + −2 − ! Z Xk=1 By the duality expression in (4.2), this becomes d d s s c √ E k k+1 Bex √ E k Bex 4 π exp − xk =4 π exp xk . − √2 ! − 2√2 ! Xk=1 Xk=1 Combining all the identities together we have proved (4.19) for A< 1,C < 1. Now we prove the case of A =1,C < 1. This time we have E Z Z 1 (4.33) uG( s1 ,..., sd ,u) du Rd+1 u1/2 Z + d d+1 E ck B 1 E 1 Z = exp xk du u exp sk (xk xk 1) . R 1/2 − − 2√2 ! × + u −4 − ! Xk=1 Z kX=1 E Z Z du uG( s1 ,..., sd ,u) R u1/2 Z + d d+1 E ck B E 1 Z du = exp xk uexp sk (xk xk 1) . R − 1/2 − 2√2 ! × + −4 − !u Xk=1 Z kX=1 Again by self-similarity and duality (4.3), we rewrite the integral on the right-hand side above as 1 d+1 du √ E Z 2 uexp s /√2(xk xk 1) R k − 1/2 + −2 − !u Z Xk=1 d s s d c √ E k k+1 Bme √ E k Bme =2 π exp − xk =2 π exp xk . − √2 ! − 2√2 ! Xk=1 Xk=1 This completes the proof for the case A = 1 and C < 1. 4.4. Proof for the case A< 1, C =1. This case of Theorem 1.2 follows from the case A = 1,C < 1 again by the particle-hole duality, which we have explained in Section 3.2. Write εj =1 τn j+1, and define − − nx ⌊ ⌋ 1 h (x)= ε . n j − 2 j=1 X b LIMITFLUCTUATIONSOFOPENASEP 27 Similarly as in (3.4), this time we have 0 nx = nx ⌊ ⌋ h (x)+ h (1) h (1 x) = , n n − n − 1 ε n(1 x) +1 nx = nx − ⌊ − ⌋ − 2 6 ⌊ ⌋ b b A,B,C,D and hn(x) x [0,1] with respect to πn has the same law as hn(x) x [0,1] { } ∈ C,D,A,B { } ∈ with respect to πn . Therefore, finite-dimensional distributions of the process 1/2 n− bhn(x) x [0,1] have the same limit as the finite-dimensional distributions of { } ∈ 1/2 the process n− hn(1 x) hn(1) x [0,1]. We get { − − } ∈ 1 f.d.d. 1 B Bme B Bme hn(x) bx [0,1] = b 1 x + 1 x ( 1 + 1 ) x [0,1] √n { } ∈ ⇒ 2√2 − − − ∈ d 1 B Bme Bme = x + 1 x 1 x [0,1] . 2√2 − − ∈ This completes the proof of Theorem 1.2. Acknowledgement. The authors benefited from references provided by G´erard Letac, Jacek Wesolowski, and Przemyslaw Matula. WB’s research was supported in part by the Charles Phelps Taft Research Center at the University of Cincinnati. YW’s research was supported in part by NSA grant H98230-16-1-0322 and Army Research Laboratory grant W911NF-17-1-0006. Appendix A. Weak Convergence from convergence of Laplace transforms It is well known that convergence of Laplace transforms in a neighborhood of 0 Rd implies weak convergence, but it is less known under what conditions convergence∈ on an open set away from the origin suffices. Hoffmann-Jørgensen [31, Section 5.14, page 378, (5.14.8)] and Mukherjea, Rao, and Suen [36, Theorem 2] independently discovered the pertinent result in the univariate case, and the argument in [36] generalizes to the multivariate setting, compare also [29, Theorem 2.1]. X(n) (n) (n) (n) Let = (X1 ,X2 ,...,Xd ) be a sequence of random vectors with Laplace transform d L (z)= L (z ,...,z )= E exp z X(n) . n n 1 d j j j=1 X Theorem A.1. Suppose that Ln(z) are finite and converge pointwise to a function L(z) for all z from an open set in Rd. If on this open set L(z) is the Laplace trans- (n) form of a random variable Y = (Y1,...,Yd), then X converges in distribution to Y . Proof. The proof is a modification of [36, proof of Theorem 2]. Denote by Pn the law of X(n) and by P the law of Y . Choose ε > 0 and c Rd such that ∞ ∈ Ln(z) L(z) for all z from the open ball z c <ε. Set Cn = Ln(c), C = L(c) and consider→ probability measures k − k ∞ 1 x c Qn(dx)= e · Pn(dx), Cn 28 WLODZIMIERZ BRYC AND YIZAO WANG with n = standing for the limiting measure. Due to our choice of c, ∞ t x Ln(t + c) L(t + c) t x e · Qn(dx)= = e · Q (dx) Rd Cn → C Rd ∞ Z ∞ Z for all t < ε. By the “usual form” of the Laplace criterion (for example, by the Cram´er–Woldk k device and Curtiss [13, Theorem 3]), this implies weak convergence of Qn to Q , i.e. for every bounded continuous function h, ∞ h(x)Qn(dx) h(x)Q (dx). Rd → Rd ∞ Z Z x c Mimicking [36] we take 0 h 1, ε > 0 and note that f(x) = e · is strictly positive, so ≤ ≤ h(x)f(x) h (x)= h(x) for all x as ε 0. ε f(x)+ ε ր ց Therefore, lim inf h(x)Pn(dx) lim inf hε(x)Pn(dx) n ≥ n →∞ Z →∞ Z h(x) h(x) = lim Cn Qn(dx)= C Q (dx) n f(x)+ ε ∞ f(x)+ ε ∞ →∞ Z Z = hε(x)P (dx). ∞ Z Taking the limit as ε 0, we get → lim inf h(x)Pn(dx) h(x)P (dx). n ≥ ∞ →∞ Z Z Applying the above to 1 h, we see that − lim h(x)Pn(dx)= h(x)P (dx) n ∞ →∞ Z Z for all continuous functions 0 h 1, which ends the proof. ≤ ≤ References [1] Askey, R. and Wilson, J. (1985). Some basic hypergeometric orthogonal polyno- mials that generalize Jacobi polynomials. Mem. Amer. Math. Soc., 54(319):iv+55. [2] Biane, P. (1998). Processes with free increments. Math. Z., 227(1):143–174. [3] Billingsley, P. (1999). Convergence of probability measures. Wiley Series in Probability and Statistics: Probability and Statistics. John Wiley & Sons Inc., New York, second edition. A Wiley-Interscience Publication. [4] Blythe, R. A. and Evans, M. R. (2007). Nonequilibrium steady states of matrix- product form: a solver’s guide. J. Phys. A, 40(46):R333–R441. [5] Blythe, R. A., Evans, M. R., Colaiori, F., and Essler, F. H. L. (2000). Exact solution of a partially asymmetric exclusion model using a deformed oscillator algebra. J. Phys. A, 33(12):2313–2332. [6] Bryc, W., Matysiak, W., and Wesolowski, J. (2007). Quadratic harnesses, q- commutations, and orthogonal martingale polynomials. Trans. Amer. Math. Soc., 359(11):5449–5483. [7] Bryc, W. and Wang, Y. (2016). The local structure of q-Gaussian processes. Probability and Mathematical Statistics, 36(2):335–252. LIMITFLUCTUATIONSOFOPENASEP 29 [8] Bryc, W. and Wang, Y. (2018). Dual representations of Laplace transforms of Brownian excursion and generalized meanders. Statist. Probab. Lett., 140:77–83. [9] Bryc, W. and Wesolowski, J. (2010). Askey–Wilson polynomials, quadratic harnesses and martingales. Ann. Probab., 38(3):1221–1262. [10] Bryc, W. and Wesolowski, J. (2017). Asymmetric simple exclusion process with open boundaries and quadratic harnesses. J. Stat. Phys., 167(2):383–415. [11] Corteel, S., Josuat-Verg`es, M., and Williams, L. K. (2011). The matrix ansatz, orthogonal polynomials, and permutations. Advances in Applied Mathematics, 46(1-4):209–225. [12] Corwin, I. and Shen, H. (2018). Open ASEP in the weakly asymmetric regime. Communications on Pure and Applied Mathematics, to appear. [13] Curtiss, J. H. (1942). A note on the theory of moment generating functions. The Annals of Mathematical Statistics, 13(4):430–433. [14] de Gier, J. and Essler, F. H. (2011). Large deviation function for the current in the open asymmetric simple exclusion process. Phys. Rev. Lett., 107(1):010602. [15] Derrida, B. (2006). Matrix ansatz large deviations of the density in exclusion processes. In International Congress of Mathematicians. Vol. III, pages 367–382. Eur. Math. Soc., Z¨urich. [16] Derrida, B. (2007). Non-equilibrium steady states: fluctuations and large deviations of the density and of the current. J. Stat. Mech. Theory Exp., 2007(7):P07023, 45. [17] Derrida, B., Domany, E., and Mukamel, D. (1992). An exact solution of a one-dimensional asymmetric exclusion model with open boundaries. J. Statist. Phys., 69(3-4):667–687. [18] Derrida, B., Enaud, C., Landim, C., and Olla, S. (2005). Fluctuations in the weakly asymmetric exclusion process with open boundary conditions. J. Stat. Phys., 118(5-6):795–811. [19] Derrida, B., Enaud, C., and Lebowitz, J. L. (2004). The asymmetric exclusion process and Brownian excursions. J. Statist. Phys., 115(1-2):365–382. [20] Derrida, B., Evans, M. R., Hakim, V., and Pasquier, V. (1993). Exact solution of a 1D asymmetric exclusion model using a matrix formulation. J. Phys. A, 26(7):1493–1517. [21] Derrida, B., Lebowitz, J. L., and Speer, E. R. (2002). Exact free energy functional for a driven diffusive open stationary nonequilibrium system. Phys. Rev. Lett., 89(3):030601, 4. [22] Derrida, B., Lebowitz, J. L., and Speer, E. R. (2003). Exact large deviation functional of a stationary open driven diffusive system: the asymmetric exclusion process. J. Statist. Phys., 110(3-6):775–810. Special issue in honor of Michael E. Fisher’s 70th birthday (Piscataway, NJ, 2001). [23] Durrett, R. T., Iglehart, D. L., and Miller, D. R. (1977). Weak convergence to Brownian meander and Brownian excursion. Ann. Probability, 5(1):117–129. [24] Enaud, C. and Derrida, B. (2004). Large deviation functional of the weakly asymmetric exclusion process. J. Statist. Phys., 114(3-4):537–562. [25] Essler, F. H. L. and Rittenberg, V. (1996). Representations of the quadratic algebra and partially asymmetric diffusion with open boundaries. J. Phys. A, 29(13):3375–3407. [26] Eyink, G., Lebowitz, J. L., and Spohn, H. (1990). Hydrodynamics of stationary nonequilibrium states for some stochastic lattice gas models. Comm. Math. Phys., 30 WLODZIMIERZ BRYC AND YIZAO WANG 132(1):253–283. [27] Eyink, G., Lebowitz, J. L., and Spohn, H. (1991). Lattice gas models in contact with stochastic reservoirs: local equilibrium and relaxation to the steady state. Comm. Math. Phys., 140(1):119–131. [28] Falconer, K. J. (2003). The local structure of random processes. J. London Math. Soc. (2), 67(3):657–672. [29] Farrell, R. H. (1976). Techniques of Multivariate Calculation, volume 520 of Lecture Notes in Mathematics. Springer. [30] Gon¸calves, P. c., Landim, C., and Milan´es, A. (2017). Nonequilibrium fluctua- tions of one-dimensional boundary driven weakly asymmetric exclusion processes. Ann. Appl. Probab., 27(1):140–177. [31] Hoffmann-Jørgensen, J. (1994). Probability with a view toward statistics. Vol. 1. Chapman & Hall. [32] Janson, S. (2007). Brownian excursion area, Wright’s constants in graph enu- meration, and other Brownian areas. Probab. Surv., 4:80–145. [33] Koekoek, R. and Swarttouw, R. F. (1996). The Askey-scheme of hypergeomet- ric orthogonal polynomials and its q-analogue. arXiv preprint math/9602214. [34] Liggett, T. M. (1985). Interacting particle systems, volume 276 of Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer-Verlag, New York. [35] MacDonald, C. T., Gibbs, J. H., and Pipkin, A. C. (1968). Kinetics of biopoly- merization on nucleic acid templates. Biopolymers, 6(1):1–25. [36] Mukherjea, A., Rao, M., and Suen, S. (2006). A note on moment generating functions. Statistics & Probability Letters, 76(11):1185–1189. [37] Pitman, J. (1999). Brownian motion, bridge, excursion, and meander charac- terized by sampling at independent uniform times. Electron. J. Probab., 4:no. 11, 33. [38] Pitman, J. (2006). Combinatorial stochastic processes, volume 1875 of Lecture Notes in Mathematics. Springer-Verlag, Berlin. Lectures from the 32nd Sum- mer School on Probability Theory held in Saint-Flour, July 7–24, 2002, With a foreword by Jean Picard. [39] Revuz, D. and Yor, M. (1999). Continuous martingales and Brownian motion, volume 293 of Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer-Verlag, Berlin, third edition. [40] Sandow, S. (1994). Partially asymmetric exclusion process with open bound- aries. Physical Review E, 50(4):2660–2667. [41] Sasamoto, T. (1999). One-dimensional partially asymmetric simple exclusion process with open boundaries: orthogonal polynomials approach. J. Phys. A, 32(41):7109–7131. [42] Sch¨utz, G. and Domany, E. (1993). Phase transitions in an exactly soluble one- dimensional exclusion process. Journal of Statistical Physics, 72(1-2):277–296. [43] Spitzer, F. (1970). Interaction of Markov processes. Advances in Math., 5:246– 290 (1970). [44] Uchiyama, M., Sasamoto, T., and Wadati, M. (2004). Asymmetric simple exclusion process with open boundaries and Askey-Wilson polynomials. J. Phys. A, 37(18):4985–5002. LIMITFLUCTUATIONSOFOPENASEP 31 [45] Uchiyama, M. and Wadati, M. (2005). Correlation function of asymmetric simple exclusion process with open boundaries. J. Nonlinear Math. Phys., 12:676– 688. [46] Wang, Y. (2018). Extremes of q-Ornstein–Uhlenbeck processes. Stochastic Process. Appl., 128(9):2979–3005. [47] Yen, J.-Y. and Yor, M. (2013). Local times and excursion theory for Brownian motion, volume 2088 of Lecture Notes in Mathematics. Springer, Cham. A tale of Wiener and Itˆomeasures. W lodzimierz Bryc, Department of Mathematical Sciences, University of Cincinnati, 2815 Commons Way, Cincinnati, OH, 45221-0025, USA. E-mail address: [email protected] Yizao Wang, Department of Mathematical Sciences, University of Cincinnati, 2815 Commons Way, Cincinnati, OH, 45221-0025, USA. E-mail address: [email protected]