Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Solution of the Kolmogorov equation for TASEP

Solution of the Kolmogorov equation for TASEP MIHAI NICA, JEREMY QUASTEL, AND DANIEL REMENIK ABSTRACT. We provide a direct and elementary proof that the formula obtained in [MQR17] for the TASEP transition probabilities for general (one-sided) initial data solves the Kolmogorov backward equation. The same method yields the solution for the related PushASEP particle system. 1. INTRODUCTION The totally asymmetric simple exclusion process (TASEP) consists of particles on the lattice Z performing totally asymmetric nearest neighbour random walks with exclusion: Each particle indepen- dently attempts jumps to the neighbouring site to the right at rate 1, the jump being allowed only if that site is unoccupied. We will always consider initial conditions in which there is a rightmost particle; this always remains so, and the positions of the particles can be denoted X (1) > X (2) >  . The t t dynamics of the first N particles X (1) > X (2) >  > X (N ) is independent of the rest, so the t t t infinite system clearly makes sense. If we let X (u) = minfk 2 Z : X (k)  ug denote the label of the rightmost particle which sits to the left of, or at, u at time t, then the TASEP height function associated to X can be defined, for z 2 Z, by 1 1 h (z) = 2 X (z 1) X (1) z; t 0 which fixes h (0) = 0. The height function itself is a simple random walk path h (z +1) = h (z)+ ^ (z) 0 t t t with  ^ (z) = 1 if there is a particle at z at time t and1 if there is no particle at z at time t. The dynamics of the height function is that local max’s become local min’s at rate 1; i.e. if h (z) = h (z  1) + 1 t t then h (z) 7! h (z) 2 at rate 1, the rest of the height function remaining unchanged. The rate of t t decrease is 1 + 21 = (r h)(r h) 1 + h where 1 is the indicator function of a local max. Hence the TASEP height function can be seen as a simple discretization of the Kardar-Parisi-Zhang (KPZ) equation, 1 2 1 2 @ h = (@ h) + @ h + : t x 2 2 Because of its amenability to computations, TASEP has become the most popular model in the KPZ universality class. N -particle TASEP was solved by Schütz [Sch97] using the coordinate Bethe ansatz. The transition probabilities are given by a determinant, P(X (1) = x ; : : : ; X (N ) = x ) = det(G (t; x X (N + 1 j))) (1.1) t 1 t N ij N +1i 0 1i;jN with n n (1) (1 w) t(w1) G (t; x) = dw e ; xn+1 2i w 0;1 where is any positively oriented simple loop which includes w = 0 and w = 1. A direct proof of 0;1 this formula is not difficult and can be obtained in a couple of pages (see, for example [MQ17]). On the other hand, one is generally less interested at a later time in the exact positions of the particles, or the exact height function, but in the joint distribution of the height function at a finite number of points x ; : : : ; x , where m is fixed, and N is large, or infinite. Or, what amounts to the same 1 m Date: January 20, 2020. arXiv:1906.01692v2 [math.PR] 17 Jan 2020 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 2 thing, one would like to compute, for some sequence n of m indices n < : : : < n and any vector 1 m a = (a ; : : : ; a ) 2 Z , the joint probability 1 m F (X ; a; n) = P (X (n ) > a ; j = 1; : : : ; m) ; (1.2) t 0 X t j j where the subscript in P denotes the initial condition. In Schütz’s formula, this would involve a sum over the positions of the other N m particles. The resulting formula is not useful, and, in particular, not conducive to the N ! 1 limit. This was overcome by [Sas05; BFPS07], who were able to rewrite the right hand side of (1.1) as a (signed) determinantal point process on a space of Gelfand-Tsetlin patterns. This allowed them to employ the Eynard-Mehta technology [EM98] to conclude that (1.2) can be written as the Fredholm TASEP determinant F (X ; a; n) = det I   K   which (in principle) could be obtained from the t 0 a a TASEP 2 operators   and K acting on ` (fn ; : : : ; n g Z) given by a 1 m f (n ; x) = f (n ; x)1 a j j xa and n j TASEP n n j i i K (n ; x ; n ; x ) = Q (x ; x )1 + (x ) (x ); (1.3) i i j j i j n <n i j t i j n k n k i j k=1 where Q(x; y) = 1 and, for k  n 1, xy x>y 1 (1 w) n t(w1) (x) = dw e ; (1.4) xX (nk) x+k+1X (nk) 0 0 2i 2 w where is any positively oriented simple loop including the pole at w = 0 but not the one at w = 1. The functions  (x), k = 0; : : : ; n 1, are defined implicitly by: I. The biorthogonality relation n n x n (x) (x) = 1 ; II. 2  (x) is a polynomial of degree at most n 1 in x for each k. k=` x2Z k ` k Except for a few very special choices of initial data, the solution for the  (x) was not discovered until [MQR17]. Without discussing their exact form, we state the final result after performing some 1 n1 (n) xy manipulations to get a nice formula. Let Q (x; y) = 2 (xyj) be the real analytic j=0 (n1)! extension of Q , let r f (x) = f (x) f (x 1) be the backwards discrete difference operator, and observe that Q is invertible (with inverse Q (x; y) = 2 1 1 ). Define x=y1 x=y 1 (1 w) r n  t(w1=2) S (z ; z ) = (e Q ) (z ; z ) = dw e ; (1.5) t;n 1 2 1 2 z z n+1+z z 2 1 2 1 2i 2 w z z +n1 2 1 t 1 (1 w) (n) r t(w1=2) S (z ; z ) = Q e (z ; z ) = dw e ; (1.6) t;n 1 2 1 2 z z n 1 2 2i 2 w the first one being defined for all n 2 N and the second one for n  1. Here (e ) is the t0 semigroup of a Poisson process with jumps to the left at rate , which we may think of as an integral t t xy 2 r t 2 2 operator on ` (Z) with kernel e (x; y) = e 1 ; this formula is actually valid for all xy xy 2 (xy)! t t r r 2 2 t 2 R and it defines the whole group of operators (e ) (in particular e is invertible, with t2R inverse e ). Define also, for n  0, epi(X ) S (z ; z ) = E S (B ; z )1 ; (1.7) 1 2 B =z t;n  2 <n t;n 0 1 where  is defined to be the hitting time of the strict epigraph of the ‘curve’ X (k + 1) k=0;:::;n1 by a discrete time random walk B with transition probabilities Q. Then, as proven in [MQR17], the kernel from (1.3) can be expressed as epi(X ) TASEP n n  0 j i K (n ;; n ;) = Q 1 + (S ) S : (1.8) i j n <n t;n t i j i t;n Now let L denote the generator of TASEP. It acts on bounded cylinder (i.e. depending on finitely many coordinates) functions f : W ! R, whereW = fX 2 Z : X (1) > X (2) > : : :g, as follows: Lf (X ) = 1 f (X (1); : : : ; X (k) + 1; : : : ) f (X (1); : : : ; X (k); : : : ) ; (1.9) X (k1)X (k)>1 k1 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 3 where we take the convention X (0) = 1 in the indicator 1 . Since f is a cylinder X (0)X (1)>1 function, the sum has only finitely many non-zero terms. If X is a Markov process with generator L, its transition probabilities p(s; X ; T; A) = P(X 2 A j X = X ) satisfy the Kolmogorov T s backward equation, @ p(s; X ; T; A) = L p(s; X ; T; A) with final condition lim p(s; X ; T; A) = s s"T 1 . The generator L acts on the X variable. In our case, we have a time homogeneous process X2A p(s; X ; T; A) = p ~(T s; X ; A) and a generating family of sets A = fX (n ) > a ; j = 1; : : : ; mg j j and F (X ; a; n) = p ~(t; X ; A), so the backward equation for F reads, F = LF (1.10) dt with initial condition F (X ; a; n) = 1 : (1.11) 0 X (n )>a ;:::X (n )>a 1 1 m m The solution of the Kolmogorov backward equation for a continuous time Markov chain on a countable state space S is unique under the condition that the supremum over all states of the rate of leaving that state is finite (see, for example, Thm. 2.26 and Cor. 2.34 of [Lig10], where it is written for the point-to-point probabilities p (X; Y ) = P (X (t) = Y ); S being countable, the probability t X measure p (X; A) = P (X (t) 2 A) = p (X; Y ) is therefore prescribed, and satisfies uniquely t X t Y2A the backward equation in t; X , with initial condition p (X; A) = 1 ). Our state space W is 0 X2A not countable; however, in TASEP (and PushASEP), the evolution of particles X (1); : : : ; X (n) is unaffected by particles X (m), m > n. Hence the backward equation (1.10), (1.11) actually takes place on the countable set W = f(X (1); X (2); : : : ; X (n )) 2 Z : X (1) > X (2) > : : : > X (n )g. n m m Therefore we have Theorem 1.1. The unique solution of the Kolmogorov backward equation (1.10), (1.11) for TASEP is TASEP given explicitly by the Fredholm determinant of the kernel K introduced in (1.3): TASEP F (X ; a; n) = det(I   K   ) 2 : (1.12) t 0 a a t ` (fn ;:::;n gZ) 1 m However, a detailed proof along the historical lines sketched above would run about 30 pages. A natural and important question is whether one could just prove directly that the determinant satisfies the Kolmogorov equation. The purpose of this article is to provide such a proof, for TASEP and its variant PushASEP. We comment that the direct proof is relatively short, but far from obvious. It is not known if it can be obtained directly from the biorthogonal representation (1.3) without knowledge of the special form of the  (x). It is also worth noting that there are examples (such as discrete time TASEP, see [MQR+]) with explicit  (x) for which the present short proof does not work. TASEP has a rich history, some of which can be found in the introduction of [Sch97]. Since that paper there has been significant progress in exact solvability, with much interest stemming from TASEP’s role as one of the fundamental models in the KPZ universality class. There are many results; we mention here only a few besides the ones which fit directly into our story above. The one-point distribution of TASEP with step initial data (X (i) = i, i  1) was first solved in [Joh00] by exploiting a Toeplitz structure which is only available for that choice of initial condition; the main goal of that paper was to prove the now-classic fact that the fluctuations of the position of an appropriately chosen particle converge, as time goes to infinity, to the Tracy-Widom GUE distribution [TW94]. This was later extended in [Joh03] to the multi-point distributions of TASEP, which were shown there to converge to the Airy process [PS02]. Much work was devoted during the last fifteen years to extending this type of results to two other choices of initial data: periodic (X (i) = 2i, i 2 Z, which was the subject of the papers [Sas05; BFPS07] discussed above) and stationary (X corresponding to placing particles on Z according to a product measure) [IS04; BFS08; BFP10]. In the past ten years there has been a huge effort to extend some of these results to the partially asymmetric simple exclusion process (see, for example, [TW08; IS11; BCS14]). These still depend on very particular initial data. In the case of TASEP on a ring there is also a huge literature with some recent breakthroughs: [BL18], where asymptotics is done for a type of step initial data, and [Pro13; Pro14], which actually computes the entire spectrum. SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 4 2. KOLMOGOROV EQUATION FOR TASEP In this section we provide a short self-contained proof of Thm. 1.1. The strategy is to compute the TASEP two sides of (1.10) with F given by (1.12) and K given by (1.8), and check that they are equal. We begin with the right hand side of (1.10). We will first consider the effect of moving a single particle, and then we will sum over all particles to obtain the effect of the generatorL. Fix a particle label k and consider the original initial condition X and the initial condition X where particle k is 0 0 moved to the right by one, i.e. X (k) = X (k) + 1 (for now it does not even matter whether or not this 0 0 particle can actually be moved without violating the strict particle order condition). We will compare TASEP the kernel K with these two initial conditions; for notational convenience we put tildes on top of TASEP all objects when they depend on the modified initial condition (e.g. K refers to the kernel started from the initial condition X ). By (1.8) we have epi(X ) epi(X ) TASEP TASEP  0 0 K K (n ;; n ;) = (S ) S S ; i j t;n t t i t;n t;n j j while by (1.7) we have epi(X ) epi(X ) 0 0 S S (z ; z ) = E S (B ; z )1 S (B ; z )1 : 1 2 B =z t;n ~  ~ 2  ~<n t;n  2 <n 0 1 t;n t;n Recall that  means the hitting time of the strict epigraph of X (` + 1) ;  ~ is the same but with the `0 modified X . It is clear from the definitions that  ~ =  unless  = k 1 < n and B = X (k) + 1 0 0 k1 (see Fig. 1 for a visual explanation), and thus the above expectation equals E S (B ; z )1 S (B ; z )1 1 1 ; B =z t;n ~  ~ 2  ~<n t;n  2 <n  =k1;B =X (k)+1 kn 0 1 0 k1 (n) which by the Markov property we may write as f (z )g (z ) with k 1 2 f (z) = P ( = k 1; B = X (k) + 1) (2.1) k B =z k1 0 and h i (n) g (z) = E S (B ; z)1 S (X (k) + 1; z) 1 ; (2.2) (k) (k) (k) t;nk+1 0 kn B =X (k)+1 k k1 0 t;n ~  ~  ~ <n where the subscript in expectation on the right hand side indicates that the random walk is started now (k) at X (k) + 1 at time k 1 and the superscript in  ~ indicates that the hitting time is now restricted to times infk; k + 1; : : :g (we will use the same notation for other hitting times below). This means that (n ) (k) TASEP TASEP  j := K K (n ;; n ;) = (S ) f g : (2.3) i j t;n k t t i i;j t;k Crucially, this is a rank one kernel (acting on ` (Z)), and then the same is true of the extended kernel TASEP TASEP 2 K K when thought of as acting on the extended space ` (fn ; : : : ; n g Z), so using 1 m t t the fact that det(I + A + B) det(I + A) = det(I + A) tr[(I + A) B] for B any rank one operator, we deduce that TASEP TASEP det(I   K   ) det(I   K   ) a a a a t t TASEP TASEP 1 (k) = det(I   K   ) tr[(I   K   )      ] a a a a a a t t (k) TASEP TASEP (here  = K K as in (2.3)). This difference corresponds to the bracket in the k-th term t t of the sum (1.9) defining the action of L on our function F (X ; a; n). The conclusion then is that if t 0 J  N is the set of labels of particles which can be moved (i.e. those so that X (k 1) X (k) > 1), 0 0 then LF (X ; a; n) equals the sum over k 2 J of the above difference, i.e. t 0 h i TASEP TASEP 1 (k) LF (X ; a; n) = det(I   K   ) tr (I   K   )      : t 0 a a a a a a t t k2J (k) Note that, as above, the sum only has finitely many non-zero terms, because  = 0 for all k > n (as is clear, for instance, from (2.2) and (2.3)). Finally we note if k 2= J (i.e. if X (k) + 1 = X (k 1)), 0 0 then since the random walk B has increments B B  1, we have fB = X (k) + 1g t1 t k1 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 5 X (1) X (2) X (3) X (4) X (4) X (5) 1 2 3 4 F IGURE 1. A sample path of the random walk B (circles) and the graph of the initial n1 n1 conditionsfX (i + 1)g andfX (i + 1)g which differ only at position k with 0 0 i=0 i=0 X (k) = X (k) + 1. The random walk has i.i.d. geometric decrements, and  and  ~ 0 0 are the first times the random walk B hits the strict epigraph of the initial condition curve with X and X respectively. Note that  ~ =  unless the hit happens exactly at 0 0 = k 1 < n and B = X (k) = X (k) + 1. In this example k = 4,  = 3 and k1 0 0 ~ > 4. fB  X (k 1) + 1g  f  k 2g and thus f (z) = 0. Thus summing over j 2= J has no k2 k effect, and we have then h i TASEP TASEP 1 (k) LF (X ; a; n) = det(I   K   ) tr (I   K   )      : (2.4) t 0 a a a a a a t t k1 Now we turn to the left hand side of (1.10). In general, if a kernel K depends smoothly on a d 1 d parameter h one has det(I K ) = det(I K ) tr[(I K ) K ]. Thus h h h h dh dh h i d TASEP TASEP 1 d TASEP F (X ; a; n) = det(I   K   ) tr (I   K   )   K t 0 a a a a a a t t t dt dt and then in view of (2.4) all we need in order to deduce (1.10) is to check that d TASEP (k) K =  : (2.5) t k1 dt We have from (1.8) that epi(X ) epi(X ) d TASEP d  0  d 0 := K (n ;; n ;) = (S ) S + (S ) ( S ): i;j i j t;n t;n t i t;n i t;n dt dt dt From (1.5) we get d  1 1 (S ) = r (S = (S r t;n t;n t;n dt 2 2 d 1 1 and from (1.6) we get S = S r = r S , which means t;n t;n t;n dt 2 2 epi(X ) d 0 1 S (z ; z ) = E r S (B ; z )1 1 2 B =z t;n  2 <n t;n 0 1 dt 2 (we have used here thatr commutes with (S ) andS since they all have Toeplitz kernels). t;n t;n This gives = (S ) H i;j t;n n i j with epi(X ) H (z ; z ) = r S (z ; z ) + E r S (B ; z )1 n 1 2 1 2 B =z t;n  2 <n t;n 0 1 epi(X ) = S (z 1; z ) E S (B 1; z )1 : 1 2 B =z t;n  2 <n t;n 0 1 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 6 Shifting the curve defined by X up by one, the first term on the second line can be expressed as E S (B 1; z )1 with  ^ the hitting time of the strict epigraph of X (m+1)+1 , B =z t;n ^  ^ 2  ^<n 0 0 1 m0 and thus H (z ; z ) = E S (B 1; z )1 S (B 1; z )1 : n 1 2 B =z t;n ^  ^ 2  ^<n t;n  2 <n 0 1 Note that  ^   , so the difference inside the brackets vanishes for   n, while if  < n and B  X ( + 1) + 2 we have  ^ =  , so the difference vanishes again. We deduce that n1 H (z ; z ) = P ( = k; B = X (k + 1) + 1) n 1 2 B =z k 0 0 1 k=0 h h i i E S (B 1; z )1 S (X (k + 1); z ) (k) (k) (k) 2 t;nk 0 2 B =X (k+1)+1 t;n ^  ^  ^ <n k 0 n1 (n) = f (z )gb (z ) k+1 1 2 k+1 k=0 (n) with f as in (2.1) and gb (z) given by the bracket in the middle line (see also the comment after k+1 (2.2)). Shifting now the curve defining  ^ back to X in the last expectation we get (n) gb (z) = E S (B ; z)1 S (X (k); z): (2.6) B =X (k) t;n  <n t;nk+1 0 k k1 0 We have then n (n ) d TASEP 1 j  j K (n ;; n ;) = (S ) f (z ) gb (z ): i j t;n k 1 2 t i dt 2 k=1 k (n) Comparing this formula and (2.3), and since g  0 for k  n by its definition (2.2), all that remains in order to prove (2.5) is that (n) (n) gb = 2g for all k  n: (2.7) k k t t r (n) r 2 2 But from (1.6) we have S e = S = Q , so applying e on the right in (2.7) and t;n 0;n (n) (n) using the definitions (2.2) and (2.6) of g and gb we see that (2.7) is equivalent to k k h i (k) (n ) (nk+1) E Q (B (k); z)1 (k) Q (X (k); z) B =X (k) k1 0   <n h i (k) (n ~ ) (nk+1) = 2E Q (B (k); z)1 (k) 2Q (X (k) + 1; z) (2.8) B =X (k)+1 k1 0  ~  ~ <n (n) z for k  n. By definition of Q , if we multiply the equation by 2 then both sides are polynomials in z, so it is enough to prove the equality for z < X (n). But it is easy to check that (n) n Q (x; y) = Q (x; y) for all x y  1: (2.9) (see [MQR17, Eqn. (2.23)]), so for such z, since X (i) is decreasing in i then z is also smaller than B , (`) ` B and X (k) in (2.8), and we deduce that all the Q ’s in the identity can be replaced by Q . As a ~ 0 consequence, what we need to prove is that h i (k) n nk+1 E Q (B ; z)1 Q (X (k); z) (k) (k) B =X (k) 0 k1 0   <n h i (k) n ~ nk+1 = 2E Q (B (k); z)1 (k) 2Q (X (k) + 1; z); B =X (k)+1 k1 0  ~  ~ <n or simply (k) (k) P B = z;   n = 2P B = z;  ~  n ; n n B =X (k) B =X (k)+1 k1 0 k1 0 which is easy to see since the walk takes Geom[ ] steps to the left: in fact, since X (k + 1) < X (k) 0 0 (k) (k) and  and  ~ may only differ on the event that the walk hits at time k 1, then for k < n, SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 7 decomposing according to the first step, (k) (k) P B = z;  ~  n = Q(X (k) + 1; y)P B = z;  ~  n n 0 B =y n B =X (k)+1 k1 0 yX (k+1) k 1 (k) 1 (k) = Q(X (k); y)P B = z;   n = P B = z;   n ; 0 B =y n B =X (k) n yX (k+1) k 0 0 2 2 k1 the case k = n is even simpler. This yields the desired equality (2.7). 3. INITIAL CONDITION Now we check that the initial condition (1.11) is satisfied. It is simpler to deal first with the one-point case (m = 1), so we start there. 3.1. One-point case. Let (n) TASEP K = K (n;; n;): t t We need to prove that (n) det I   K   = 1 : (3.1) a a X (n)>a 0 0 ` (Z) Remark 3.1. In the one point case, the initial condition (3.1) can be checked rather easily from the biorthogonal ensemble representation (1.3), as mentioned by one of the referees; the proof using this biorthogonalization is given below. However, since our goal in this paper is to prove (1.12) directly from the representation (1.8) (and thus avoid having to use the result of [MQR17] to connect this formula back to (1.3)), we provide the direct (but relatively more complicated) proof from (1.8) in this subsection. (n) n1 n n From (2.3) we have K (z ; z ) = (z ) (z ) with t = 0 in the definition (1.4). Let 1 2 1 2 0 k=0 k k n n ;  : Z  f0; : : : ; n 1g ! R be given by (z; k) = (z) and (z; k) =  , so that the left k k hand side of (3.1) equals det(I      ) = det(I    ); a a a ` (Z) where we have used the cyclic property of the Fredholm determinant. Note that on the determinant on the right hand side is finite dimensional, with the n n matrix inside being given by n n (`; k) = (z) (z); `; k = 0; : : : ; n 1: k ` za Assume first that a < X (n). Then a < X (n k) k for all k = 0; : : : ; n, which implies from (1.4) 0 0 that (z) = 0 for all z  a and all such k, so det(I    ) = det(I ) = 1 as desired. Assume otherwise that a  X (n). From (1.4) it is straightforward to check that (z) = 1 and 0 z=X (n) n n n thus the first column of    is given by    (`; 0) = (z) (z) =  (X (n)). But a a 0 za k ` ` n n n the biorthogonalization condition (z) (z) = 1 implies that  (X (n)) = 1 . We `=0 0 `=0 z2Z ` ` deduce that the first column of I    equals zero, so det(I    ) vanishes as desired. a a We turn now to the proof of (3.1) directly from (1.8). Consider first the case a  X (n) 1, which (n) implies that in (3.1) we need to evaluate K (z ; z ) only for z  X (n) 1. By definition of the 1 2 1 0 epi(X ) (n) hitting time  in (1.7) we haveS (z ; z ) = Q (z ; z ) for z > X (1), so we may write 1 2 1 2 1 0 0;n (n) epi(X ) n (n) n 0 K = Q  Q + Q   S : (3.2) X (1) X (1) 0 0 0 0;n Now n yx+n yx n Q (x; y) = (1) 2 1 (3.3) 0yxn yx so, in particular, Q (x; y) = 0 for y x > n. If y > X (1) and z  X (n) 1 we have 0 1 0 y z > X (1) X (n) + 1  n, so Q (z ; y) = 0 and the first term on the right hand side of (3.2) 1 0 0 1 epi(X ) 0 n1 vanishes. For the second term we consider S (y; z ) = E S (B ; z )1 with 2 B =y 0;nk k 2  =k 0;n k=0 y  X (1) and note that  = k implies B > X (k + 1), which cannot happen unless the starting 0 k 0 location satisfies y > X (k + 1) + k (because the walk takes downward steps of size at least 1). But 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 8 for such y (and z  X (n) 1) we have y z > X (k + 1) + k X (n) + 1  n, and thus again 1 0 1 0 0 we get Q (z ; y) = 0. We have proved that (n) K = 0 when a < X (n); (3.4) a 0 and then the left hand side of (3.1) equals det(I ) = 1 for a < X (n), as desired. Next we turn to the case a  X (n), and focus on the column corresponding to the index X (n) in 0 0 (n) (n) K   , that is K (z; X (n)) . We claim that, in this case, a a 0 0 0 za (n) K (z; X (n)) = 1 for all z 2 Z: (3.5) z=X (n) 0 0 (n) If that is the case then the matrix I   K   has a column which is identically 0, and thus its a a determinant vanishes as required. In order to prove (3.5), observe first that the argument we used to prove (3.4) actually shows that (n) K (z; X (n)) = 0 for all z < X (n). Consider next the case z = X (n). For the first term in (3.2), 0 0 0 n (n) a simple residue computation using the contour integral formulas for Q and Q which follow from (1.5) and (1.6) with t = 0 gives (remembering X (n) X (1)  1 n) 0 0 H H n X (n)X (1)+n1 0 0 (1w) (1v) n (n) 1 1 Q  Q (X (n); X (n)) = dw dv 0 0 2 X (1) X (n)X (1)+n n 0 (2i) 1vw w 0 0 v (3.6) = 1 : X (n)X (1)=1n 0 0 (n) Since we want to show that K (X (n); X (n)) = 1, then from (3.6) together with (3.2) we see that 0 0 we need to show that epi(X ) n 0 Q   S (X (n); X (n)) = 1 : (3.7) 0 0 X (1) X (n)X (1)<1n 0 0;n 0 0 We consider the two possibilities for X (n) X (1) separately. Assume first that X (1) X (n) = 0 0 0 0 n 1. This means that the first n particles are packed one next to the other, i.e. X (k) = X (1)k + 1, 0 0 and thus, since the walk starts below X (1) and takes downward steps of size at least 1, it cannot hit epi(X ) the strict epigraph of the curve by time n, so   S = 0 as desired. The other possibility X (1) 0 0;n epi(X ) (n) is that X (1) X (n) > n 1. Using S = Q we may express the kernel S (y; X (n)) 0 0 0;n 0 0;n n1 (nk) as E Q (B ; X (n))1 and then note that, since inside the expectation B is the B =y k 0  =k k k=0 location where the strict epigraph is hit, B X (n)  X (k + 1) + 1 X (n)  n k  1, k 0 0 0 n1 (nk) nk so by (2.9) we may replace Q by Q there to get P (B = X (n);  = k). But B =y n 0 k=0 B = X (n) implies B  X (n) + 1, which in turn implies   n 1, so the last sum equals n 0 n1 0 simply Q (y; X (n)), and then epi(X ) n n n Q   S (X (n); X (n)) = Q (X (n); y)Q (y; X (n)): 0 0 0 0 X (1) 0 0;n yX (1) And the last sum can be extended to all y’s because if y > X (1) then, under our assumption X (n) X (1) < 1 n and using (3.3), the first factor in the summand vanishes. This shows that, in 0 0 epi(X ) this case, Q   S (X (n); X (n)) = 1 as desired. X (1) 0 0 0 0;n Up to here we have proved (3.5) for z  X (n). Consider finally the case z > X (n). The exact 0 0 same argument used in the last paragraph to prove (3.7) also shows epi(X ) n 0 n (n) Q   S (z; X (n)) = Q  Q (z; X (n))1 : 0 0 X (1) X (1) X (1)X (n)>n1 0 0;n 0 0 0 (n) In particular, in the case X (1) X (n) > n 1 we get using (3.2) that K (z; X (n)) = 0 0 0 n (n) n (n) Q Q (z; X (n)) which equals 0 because in fact Q Q  0 (see [MQR17, Eqn. (2.24)]). It remains to prove that the first term in (3.2) vanishes when X (1) X (n) = n 1, but this is 0 0 straightforward, because in this case we have (using the contour integral formulas (1.5) and (1.6) as before) H H n H (1w) 1 1 1 1 Q  S (z; X (n)) = dw dv = dw = 0 0;n 0 2 X (1) zX (1)+n n zX (1)+n 0 1vw 2i (2i) 0 0 w v w since z X (1) + n > X (n) X (1) + n = 1. 0 0 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 9 TASEP TASEP 3.2. Multi-point case. We let L = K (n ;; n ;) and think of K as an operator-valued i;j i j matrix L with entries L . The initial condition (1.11) will follow from the arguments in the one-point i;j case and the identity (n ) n n n n i j i j i L = Q 1 + Q K ; (3.8) i;j n <n i j t which is (D.8) in [MQR17]. Consider first the case a  X (n ) 1 for each i. Here (3.4) gives i 0 i (n ) (n ) (n ) n n i n n i i j i j i K = 0, and then in the sum   Q K (z ; z ) = 1 Q (z ; y)K (y; z ) a a 1 2 z a 1 2 i i t 1 i y2Z 0 only terms with y > a survive. On the other hand, we have Q (z; y) = 0 for z < y (any m 2 Z) so n n n n j i j i Q (z ; y) = 0 for z  a < y. This shows that   L   = Q 1 , and since this 1 1 i a i;j a n <n i j i j holds for all i, we get F (X ; a; n) = det(I  L  ) with L a strictly upper triangular operator-valued 0 0 a a matrix, which yields F (X ; a; n) = 1 as desired. 0 0 Suppose next that a  X (n ) for some j, so we need to prove det(I   L  ) = 0. We may j 0 j a a assume without loss of generality that n is the largest particle label amongfn ; : : : ; n g for which j 1 m this inequality holds. We will focus on the j-th column of   L  , and more precisely on the subcolumn a a (j) (j) of this column of kernels corresponding to the index X (n ), that is v := v (z) with 0 j ` z2Z; i=1;:::;m (j) (n ) v (z) =   L (z; X (n )). In this case (3.5) gives us K (z; X (n )) = 1 for all z, and a i;j 0 j 0 j z=X (n ) i 0 0 ` j (n ) n n n n j i j i then Q K (z; X (n )) = Q (z; X (n )), so that using (3.8) we get 0 j 0 j (j) n n j i v (z) = Q (z; X (n ))1 1 : 0 j n n za i i j i (j) The case i = j is direct and yields v (z) = 1 . Otherwise, if n > n , then necessarily a < i j i z=X (n ) j 0 j X (n ) by our choice of j, and hence for z  a we have X (n )z > X (n )X (n )  n n 1 0 i i 0 j 0 j 0 i i j (j) (j) so as before v (z) = 0. The conclusion is that v (z) = 1 1 . But then I   L  has a i=j z=X (n ) a a i i 0 j (j) column (namely v ) which is identically 0, and thus F (X ; a; n) = det(I   L  ) = 0 as required. 0 0 a a 4. PUSHASEP We consider now PushASEP, a generalization of TASEP and the Toom model [DLSS91] introduced in [BF08]. Again we have particles on the lattice Z with positions X (1) > X (2) >  . Each particle t t attempts to jump one step to the right at rate r, with jumps being permitted only if the neighboring site is empty; this is the same TASEP dynamics considered above (except run at rate r). On the other hand, each particle has another (independent) exponential clock running at rate `, and when it rings the particle jumps to the nearest vacant site on its left. Relabeling the particles in order to keep the ordering after such a jump, we may think of the effect of a particle jumping left as pushing all its left neighbors one step to the left. TASEP and the Toom model (or PushTASEP) are recovered from PushASEP by setting ` = 0 and r = 0, respectively. The generatorL of this process can be written similary to the one for TASEP: letting b(k) denote r;` the length of the block of nearest neighbor occupied sites lying to the left of (and including) particle k, meaning that X (k + j) = X (k) j for j = 0; : : : ; b(k) 1 and X (m + b(k)) < X (k) b(k), 0 0 0 0 L f (X ) = r 1 f (: : : ; X (k) + 1; : : : ) f (: : : ; X (k); : : : ) r;` X (k1)X (k)>1 k1 + ` f (: : : ; X (k) 1; : : : ; X (k + b(k)) 1; : : : ) f (: : : ; X (k); : : : ; X (k + b(k)); : : : ) ; k1 where again we take X (0) = 1. As for the TASEP case we define F (X ; a; n) = P X (n ) > a ; : : : ; X (n ) > a ; t X t 1 1 t m m again we have that F (; a) satisfies the Kolmogorov backward equation F = L F; (4.1) r;` dt with initial condition F (X ; a; n) = 1 : (4.2) 0 X (1)>a ;:::X (n)>a 1 n SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 10 Theorem 4.1. The unique solution of (4.1), (4.2) is given by PushASEP F (X ; a; n) = det(I   K   ) 2 t 0 a a t ` (fn ;:::;n gZ) 1 m with epi(X ) PushASEP n n  0 j i K (n ;; n ;) = Q 1 + (S ) S ; i j n <n t;n i j i t;n where the operatorsS and S are now given by t;n t;n 1 + rtr +2`tr n S (z ; z ) = (e Q ) (z ; z ) (4.3) t;n 1 2 1 2 1 (1 w) t[r(w1=2)+`(1=w2)] = dw e ; z z n+1+z z 2 1 2 1 2i 2 w 1 + (n) rtr 2`tr S (z ; z ) = Q e (z ; z ) (4.4) t;n 1 2 1 2 z z +n1 2 1 1 (1 w) t[r(w1=2)+`(21=(1w))] = dw e ; z z n 1 2 2i 2 w 1 t rtr +2`tr r 2 2 herer f (x) = f (x + 1) f (x) and e is defined similarly to e after (1.6). This result, which is new, can be proved along the same lines as the solution of TASEP given in [MQR17]; such a proof will appear in [MQR+] (where convergence to the KPZ fixed point will also be proved). Our goal here, as for the Kolmorogov equation for TASEP in Thm. 1.1, is to prove directly that the Fredholm determinant provides a solution to the Kolmogorov backward equation (4.1), (4.2). Note first that for t = 0 the formulas forS andS given in Thm. 4.1 are the same as those t;n t;n for TASEP, so the initial condition (4.2) follows from the TASEP proof. Next we claim that in order to prove that the Fredholm determinant solves (4.1) it is enough to consider the case of PushTASEP, i.e. r = 0, since the TASEP part of the dynamics was already handled in Sec. 2. To see this it is convenient PushASEP to define F (X ; a) = det(I   K   ) 2 with t ;t 0 a a 1 2 t ;t ` (fn ;:::;n gZ) 1 2 1 m epi(X ) PushASEP n n  0 j i K (n ;; n ;) = Q 1 + (S ) S ; i j n <n t ;t ;n t ;t i j 1 2 i t ;t ;n 1 2 1 2 j 1 + 1 + epi(X ) rt r +2`t r n  (n) rt r 2`t r 0 1 2 1 2 2 2 S = (e Q ) , S = Q e , and S defined t ;t ;n t ;t ;n 1 2 1 2 t ;t ;n 1 2 analogously using S , so that (4.1) will follow (after setting t = t = t) if we prove that t ;t ;n 1 2 1 2 j @ @ ( + )F (X ; a) = L F (X ; a). But since L = L +L , it is actually enough to t ;t 0 `;r t ;t 0 `;r `;0 0;r 1 2 1 2 @t @t 1 2 prove that @ @ F (X ; a) = L F (X ; a) and F (X ; a) = L F (X ; a): (4.5) t ;t 0 0;r t ;t 0 t ;t 0 t ;t 0 `;0 1 2 1 2 1 2 1 2 @t @t 1 2 The first of these is a version of the Kolmogorov equation for TASEP (1.1) where in the contour integrals in (1.5) and (1.6) there are respectively additional factors (w) and (w) inside the integrands, both independent of t, and it is not hard to see that the proof in Sec. 2 works without any difficulty after this modification. We will prove below that (4.1) holds in the PushTASEP case r = 0, and inspecting the proof it is easy again to seen that this implies the second equation in (4.5). We conclude that it suffices to check the case r = 0, and we can also take ` = 1 as it is just a scaling factor. Given this reduction, one may naturally wonder why the proof works for PushASEP and not ASEP. It just turns out that when one tries to adapt Schütz’s formula [Sch97] to include jumps to the left (for example, as done in [BF08]) the resulting formula ends up solving the push dynamics instead of the leftward TASEP dynamics. Having made these comments, we begin the proof. Proceeding as for TASEP, let k 2 f1; : : : ; ng and consider the initial condition X and the one where particle k is moved one step to the left along with its neighboring block. Write b for the length of this block. As for TASEP, we need to compute epi(X ) epi(X ) 0 0 S S (z ; z ) = E S (B ; z )1 S (B ; z )1 ; 1 2 B =z t;n ~  ~ 2  ~<n t;n  2 <n t;n t;n 0 1 where now the tildes refer to the initial condition modified in this new way. We clearly have  =  ~ unless ~ 2 fk1; : : : ; k +b2g and B = X ( +1). But since the particles X (j) for j 2 fk; : : : ; k +b1g ~ 0 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 11 lie in a block,  ~ 2 fk 1; : : : ; k + b 2g implies  ~ = k 1 (because the walk always jumps down by at least one step). This means that the above expectation equals (n) E S (B ; z )1 S (B ; z )1 1 = f (z )g (z ) B =z t;n ~  ~ 2  ~<n t;n  2 <n  ~=k1;B =X (k) k 1 2 0 1 0 k1 k with f (z) = P ( ~ = k 1; B = X (k)) k B =z k1 0 (k) and (see also the comment after (2.2) for the meaning of  ) h i (n) g (z) = S (X (k); z) E S (k) (B (k); z)1 (k) 1 t;nk+1 0 kn B =X (k) k k1 0 t;n   <n (n) (k) (note that f does not depend on b but g does, through  ). Since this is a rank one kernel, the same arguments as in the TASEP case give h i PushTASEP PushTASEP 1 (k) L F (X ; a; n) = det(I   K   ) tr (I   K   ) 0;1 t 0 a a a a a a t t k=1 (n) (k) where  is defined exactly as in (2.3) but using the new versions of f and g (note that here L k 0;1 acts on all particles, since all particles are allowed to move in the PushTASEP dynamics). Now we consider the left hand side. As for TASEP, all we need to check in order to finish the proof is that d PushTASEP (k) K =  : t k1 dt Now we have epi(X ) epi(X ) d PushTASEP d  0  d 0 = K (n ;; n ;) = (S ) S + (S ) ( S ); i;j i j t;n t;n t i t;n i t;n dt dt dt j j d  + d  + and using (4.3) and (4.4) we have (S ) = 2(S r and (S ) = 2r (S . t;n t;n t;n t;n dt dt Therefore = 2S H i;j t;n n i j with epi(X ) + 0 + H (z ; z ) = r S (z ; z ) E r S (B ; z )1 n 1 2 1 2 B =z t;n  2 <n 0 1 t;n epi(X ) = S (z + 1; z ) E S (B + 1; z )1 : 1 2 B =z t;n  2 <n t;n 0 1 The first term on the second line can be expressed as E S (B + 1; z )1 with  ^ the B =z t;n ^  ^ 2  ^<n 0 1 hitting time of the strict epigraph of X (m + 1) 1 , and thus m0 H (z ; z ) = E S (B + 1; z )1 S (B + 1; z )1 : n 1 2 B =z 2 t;n  2 <n t;n ^  ^  ^<n 0 1 Now we have  ^   , so the difference inside the brackets vanishes for  ^  n, and it also vanishes when ^ < n and B  X ( ^ + 1) + 1. Thus H (z ; z ) equals ^ 0 n 1 2 n1 P ( ^ = k; B = X (k 1)) B =z k 0 0 1 k=0 h i (k) S (X (k + 1) + 1; z ) E S (B + 1; z )1 (k) t;nk 0 2 t;n 2 B =X (k+1) k 0  <n n1 (n) = f (z )gb (z ) k+1 1 2 k+1 k=0 with f (z) = P ( ^ = k 1; B = X (k)) k B =z k1 0 and h i (n) gb (z) = S (X (k) + 1; z) E S (B + 1; z)1 : (k) (k) (k) t;nk+1 0 B =X (k) t;n   <n k k1 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 12 We have then (n ) = 2 S f (z ) gb (z ); i;j t;n k 1 2 k=1 (n) (n) and so it remains to prove that f = f for all k and gb = g for k  n (using as for TASEP that k k k 2 k (n) g  0 for k > n). For fixed k let  ~ denote the hitting time of the initial condition X introduced above where the block starting at X (k) is moved one step to the left, so that f (z) = P ( ~ = k 1; B = X (k)). 0 k B =z k k1 0 We obviously have  ^   ~ , and thus  ^ = k 1 and B = X (k) implies  ~ = k 1, so that k k1 0 k f (z) = P ( ^ =  ~ = k 1; B = X (k)). This gives k B =z k k1 0 f (z) f (z) = P ( ~ >  ^ = k 1; B = X (k)) = 0: k k B =z k k1 0 (n) (n) 1 2tr Next we prove the identity gb = g for k  n. Proceeding as for TASEP, we apply e on the k 2 k right of this identity to see that it is equivalent to h i (k) (nk+1) (n ) Q (X (k) + 1; z) E Q (B + 1; z)1 (k) (k) 0 B =X (k) k1 0   <n h i (k) 1 (nk+1) 1 (n ) = Q (X (k); z) E Q (B (k); z)1 (k) : B =X (k) k1 0   <n 2 2 Following again the TASEP argument, the two sides are of the form 2 times a polynomial in z, so it is enough to prove the equality for z < X (n), in which case the identity becomes (k) 1 (k) P B = z 1;   n = P B = z;   n ; n n B =X (k) B =X (k) k1 0 2 k1 0 which is easy to prove again since the walk takes Geom[ ] steps to the left and z  X (n). Acknowledgements. MN and JQ were supported by the Natural Sciences and Engineering Research Council of Canada. DR was supported by Programa Iniciativa Científica Milenio grant number NC120062 through Nucleus Millenium Stochastic Models of Complex and Disordered Systems, by Conicyt Basal-CMM Proyecto/Grant PAI AFB-170001, and by Fondecyt Grant 1160174. R EFERENCES [BFP10] J. Baik, P. L. Ferrari, and S. Péché. Limit process of stationary TASEP near the characteristic line. Comm. Pure Appl. Math. 63.8 (2010), pp. 1017–1070. [BL18] J. Baik and Z. Liu. Fluctuations of TASEP on a ring in relaxation time scale. Comm. Pure Appl. Math. 71.4 (2018), pp. 747–813. [BCS14] A. Borodin, I. Corwin, and T. Sasamoto. From duality to determinants for q-TASEP and ASEP. Ann. Probab. 42.6 (2014), pp. 2314–2382. [BF08] A. Borodin and P. L. Ferrari. Large time asymptotics of growth models on space-like paths. I. PushASEP. Electron. J. Probab. 13 (2008), no. 50, 1380–1418. [BFPS07] A. Borodin, P. L. Ferrari, M. Prähofer, and T. Sasamoto. Fluctuation properties of the TASEP with periodic initial configuration. J. Stat. Phys. 129.5-6 (2007), pp. 1055–1080. [BFS08] A. Borodin, P. L. Ferrari, and T. Sasamoto. Transition between Airy and Airy processes 1 2 and TASEP fluctuations. Comm. Pure Appl. Math. 61.11 (2008), pp. 1603–1629. [DLSS91] B. Derrida, J. L. Lebowitz, E. R. Speer, and H. Spohn. Dynamics of an anchored Toom interface. J. Phys. A 24.20 (1991), pp. 4805–4834. [EM98] B. Eynard and M. L. Mehta. Matrices coupled in a chain. I. Eigenvalue correlations. J. Phys. A 31.19 (1998), pp. 4449–4456. [IS04] T. Imamura and T. Sasamoto. Fluctuations of the one-dimensional polynuclear growth model with external sources. Nuclear Phys. B 699.3 (2004), pp. 503–544. [IS11] T. Imamura and T. Sasamoto. Current moments of 1D ASEP by duality. J. Stat. Phys. 142.5 (2011), pp. 919–930. [Joh00] K. Johansson. Shape fluctuations and random matrices. Comm. Math. Phys. 209.2 (2000), pp. 437–476. SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 13 [Joh03] K. Johansson. Discrete polynuclear growth and determinantal processes. Comm. Math. Phys. 242.1-2 (2003), pp. 277–329. [Lig10] T. M. Liggett. Continuous time Markov processes. Vol. 113. Graduate Studies in Mathemat- ics. An introduction. American Mathematical Society, Providence, RI, 2010, pp. xii+271. [MQ17] K. Matetski and J. Quastel. From the totally asymmetric simple exclusion process to the KPZ fixed point. 2017. arXiv: 1710.02635[math.PR]. [MQR17] K. Matetski, J. Quastel, and D. Remenik. The KPZ fixed point. 2017. arXiv: 1701 . 00018[math.PR]. [MQR+] K. Matetski, J. Quastel, and D. Remenik. TASEP and generalizations: Method for exact solution. In preparation. 2019. [PS02] M. Prähofer and H. Spohn. Scale invariance of the PNG droplet and the Airy process. J. Stat. Phys. 108.5-6 (2002), pp. 1071–1106. [Pro13] S. Prolhac. Spectrum of the totally asymmetric simple exclusion process on a periodic lattice—bulk eigenvalues. J. Phys. A 46.41 (2013), pp. 415001, 36. [Pro14] S. Prolhac. Spectrum of the totally asymmetric simple exclusion process on a periodic lattice-first excited states. J. Phys. A 47.37 (2014), pp. 375001, 29. [Sas05] T. Sasamoto. Spatial correlations of the 1D KPZ surface on a flat substrate. Journal of Physics A: Mathematical and General 38.33 (2005), p. L549. [Sch97] G. M. Schütz. Exact solution of the master equation for the asymmetric exclusion process. J. Statist. Phys. 88.1-2 (1997), pp. 427–445. [TW94] C. A. Tracy and H. Widom. Level-spacing distributions and the Airy kernel. Comm. Math. Phys. 159.1 (1994), pp. 151–174. [TW08] C. A. Tracy and H. Widom. A Fredholm determinant representation in ASEP. J. Stat. Phys. 132.2 (2008), pp. 291–300. (M. Nica) DEPARTMENT OF M ATHEM ATICS, UNIVERSITY OF TORONTO, 40 S T. G EORGE STREET, TORONTO, ONTARIO, CANADA M5S 2E4 E-mail address: mnica@math.toronto.edu (J. Quastel) DEPARTMENT OF M ATHEM ATICS, UNIVERSITY OF TORONTO, 40 S T. GEORGE STREET, TORONTO, ONTARIO, CANADA M5S 2E4 E-mail address: quastel@math.toronto.edu (D. Remenik) DEPARTAMENTO DE I NGENIERÍA M ATEMÁTICA AND C ENTRO DE M ODELAMIENTO MATEMÁTICO (UMI-CNRS 2807), UNIVERSIDAD DE C HILE, AV. BEAUCHEF 851, TORRE NORTE, PISO 5, SANTIAGO, CHILE E-mail address: dremenik@dim.uchile.cl http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Mathematical Physics arXiv (Cornell University)

Solution of the Kolmogorov equation for TASEP

Nica, Mihai; Quastel, Jeremy; Remenik, Daniel
Mathematical Physics , Volume 2020 (1906) – Jun 4, 2019
Download PDF
13 pages

Loading next page...
 
/lp/arxiv-cornell-university/solution-of-the-kolmogorov-equation-for-tasep-7c3EptS3IY

References (27)

ISSN
0091-1798
eISSN
ARCH-3337
DOI
10.1214/20-AOP1425
Publisher site
See Article on Publisher Site

Abstract

MIHAI NICA, JEREMY QUASTEL, AND DANIEL REMENIK ABSTRACT. We provide a direct and elementary proof that the formula obtained in [MQR17] for the TASEP transition probabilities for general (one-sided) initial data solves the Kolmogorov backward equation. The same method yields the solution for the related PushASEP particle system. 1. INTRODUCTION The totally asymmetric simple exclusion process (TASEP) consists of particles on the lattice Z performing totally asymmetric nearest neighbour random walks with exclusion: Each particle indepen- dently attempts jumps to the neighbouring site to the right at rate 1, the jump being allowed only if that site is unoccupied. We will always consider initial conditions in which there is a rightmost particle; this always remains so, and the positions of the particles can be denoted X (1) > X (2) >  . The t t dynamics of the first N particles X (1) > X (2) >  > X (N ) is independent of the rest, so the t t t infinite system clearly makes sense. If we let X (u) = minfk 2 Z : X (k)  ug denote the label of the rightmost particle which sits to the left of, or at, u at time t, then the TASEP height function associated to X can be defined, for z 2 Z, by 1 1 h (z) = 2 X (z 1) X (1) z; t 0 which fixes h (0) = 0. The height function itself is a simple random walk path h (z +1) = h (z)+ ^ (z) 0 t t t with  ^ (z) = 1 if there is a particle at z at time t and1 if there is no particle at z at time t. The dynamics of the height function is that local max’s become local min’s at rate 1; i.e. if h (z) = h (z  1) + 1 t t then h (z) 7! h (z) 2 at rate 1, the rest of the height function remaining unchanged. The rate of t t decrease is 1 + 21 = (r h)(r h) 1 + h where 1 is the indicator function of a local max. Hence the TASEP height function can be seen as a simple discretization of the Kardar-Parisi-Zhang (KPZ) equation, 1 2 1 2 @ h = (@ h) + @ h + : t x 2 2 Because of its amenability to computations, TASEP has become the most popular model in the KPZ universality class. N -particle TASEP was solved by Schütz [Sch97] using the coordinate Bethe ansatz. The transition probabilities are given by a determinant, P(X (1) = x ; : : : ; X (N ) = x ) = det(G (t; x X (N + 1 j))) (1.1) t 1 t N ij N +1i 0 1i;jN with n n (1) (1 w) t(w1) G (t; x) = dw e ; xn+1 2i w 0;1 where is any positively oriented simple loop which includes w = 0 and w = 1. A direct proof of 0;1 this formula is not difficult and can be obtained in a couple of pages (see, for example [MQ17]). On the other hand, one is generally less interested at a later time in the exact positions of the particles, or the exact height function, but in the joint distribution of the height function at a finite number of points x ; : : : ; x , where m is fixed, and N is large, or infinite. Or, what amounts to the same 1 m Date: January 20, 2020. arXiv:1906.01692v2 [math.PR] 17 Jan 2020 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 2 thing, one would like to compute, for some sequence n of m indices n < : : : < n and any vector 1 m a = (a ; : : : ; a ) 2 Z , the joint probability 1 m F (X ; a; n) = P (X (n ) > a ; j = 1; : : : ; m) ; (1.2) t 0 X t j j where the subscript in P denotes the initial condition. In Schütz’s formula, this would involve a sum over the positions of the other N m particles. The resulting formula is not useful, and, in particular, not conducive to the N ! 1 limit. This was overcome by [Sas05; BFPS07], who were able to rewrite the right hand side of (1.1) as a (signed) determinantal point process on a space of Gelfand-Tsetlin patterns. This allowed them to employ the Eynard-Mehta technology [EM98] to conclude that (1.2) can be written as the Fredholm TASEP determinant F (X ; a; n) = det I   K   which (in principle) could be obtained from the t 0 a a TASEP 2 operators   and K acting on ` (fn ; : : : ; n g Z) given by a 1 m f (n ; x) = f (n ; x)1 a j j xa and n j TASEP n n j i i K (n ; x ; n ; x ) = Q (x ; x )1 + (x ) (x ); (1.3) i i j j i j n <n i j t i j n k n k i j k=1 where Q(x; y) = 1 and, for k  n 1, xy x>y 1 (1 w) n t(w1) (x) = dw e ; (1.4) xX (nk) x+k+1X (nk) 0 0 2i 2 w where is any positively oriented simple loop including the pole at w = 0 but not the one at w = 1. The functions  (x), k = 0; : : : ; n 1, are defined implicitly by: I. The biorthogonality relation n n x n (x) (x) = 1 ; II. 2  (x) is a polynomial of degree at most n 1 in x for each k. k=` x2Z k ` k Except for a few very special choices of initial data, the solution for the  (x) was not discovered until [MQR17]. Without discussing their exact form, we state the final result after performing some 1 n1 (n) xy manipulations to get a nice formula. Let Q (x; y) = 2 (xyj) be the real analytic j=0 (n1)! extension of Q , let r f (x) = f (x) f (x 1) be the backwards discrete difference operator, and observe that Q is invertible (with inverse Q (x; y) = 2 1 1 ). Define x=y1 x=y 1 (1 w) r n  t(w1=2) S (z ; z ) = (e Q ) (z ; z ) = dw e ; (1.5) t;n 1 2 1 2 z z n+1+z z 2 1 2 1 2i 2 w z z +n1 2 1 t 1 (1 w) (n) r t(w1=2) S (z ; z ) = Q e (z ; z ) = dw e ; (1.6) t;n 1 2 1 2 z z n 1 2 2i 2 w the first one being defined for all n 2 N and the second one for n  1. Here (e ) is the t0 semigroup of a Poisson process with jumps to the left at rate , which we may think of as an integral t t xy 2 r t 2 2 operator on ` (Z) with kernel e (x; y) = e 1 ; this formula is actually valid for all xy xy 2 (xy)! t t r r 2 2 t 2 R and it defines the whole group of operators (e ) (in particular e is invertible, with t2R inverse e ). Define also, for n  0, epi(X ) S (z ; z ) = E S (B ; z )1 ; (1.7) 1 2 B =z t;n  2 <n t;n 0 1 where  is defined to be the hitting time of the strict epigraph of the ‘curve’ X (k + 1) k=0;:::;n1 by a discrete time random walk B with transition probabilities Q. Then, as proven in [MQR17], the kernel from (1.3) can be expressed as epi(X ) TASEP n n  0 j i K (n ;; n ;) = Q 1 + (S ) S : (1.8) i j n <n t;n t i j i t;n Now let L denote the generator of TASEP. It acts on bounded cylinder (i.e. depending on finitely many coordinates) functions f : W ! R, whereW = fX 2 Z : X (1) > X (2) > : : :g, as follows: Lf (X ) = 1 f (X (1); : : : ; X (k) + 1; : : : ) f (X (1); : : : ; X (k); : : : ) ; (1.9) X (k1)X (k)>1 k1 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 3 where we take the convention X (0) = 1 in the indicator 1 . Since f is a cylinder X (0)X (1)>1 function, the sum has only finitely many non-zero terms. If X is a Markov process with generator L, its transition probabilities p(s; X ; T; A) = P(X 2 A j X = X ) satisfy the Kolmogorov T s backward equation, @ p(s; X ; T; A) = L p(s; X ; T; A) with final condition lim p(s; X ; T; A) = s s"T 1 . The generator L acts on the X variable. In our case, we have a time homogeneous process X2A p(s; X ; T; A) = p ~(T s; X ; A) and a generating family of sets A = fX (n ) > a ; j = 1; : : : ; mg j j and F (X ; a; n) = p ~(t; X ; A), so the backward equation for F reads, F = LF (1.10) dt with initial condition F (X ; a; n) = 1 : (1.11) 0 X (n )>a ;:::X (n )>a 1 1 m m The solution of the Kolmogorov backward equation for a continuous time Markov chain on a countable state space S is unique under the condition that the supremum over all states of the rate of leaving that state is finite (see, for example, Thm. 2.26 and Cor. 2.34 of [Lig10], where it is written for the point-to-point probabilities p (X; Y ) = P (X (t) = Y ); S being countable, the probability t X measure p (X; A) = P (X (t) 2 A) = p (X; Y ) is therefore prescribed, and satisfies uniquely t X t Y2A the backward equation in t; X , with initial condition p (X; A) = 1 ). Our state space W is 0 X2A not countable; however, in TASEP (and PushASEP), the evolution of particles X (1); : : : ; X (n) is unaffected by particles X (m), m > n. Hence the backward equation (1.10), (1.11) actually takes place on the countable set W = f(X (1); X (2); : : : ; X (n )) 2 Z : X (1) > X (2) > : : : > X (n )g. n m m Therefore we have Theorem 1.1. The unique solution of the Kolmogorov backward equation (1.10), (1.11) for TASEP is TASEP given explicitly by the Fredholm determinant of the kernel K introduced in (1.3): TASEP F (X ; a; n) = det(I   K   ) 2 : (1.12) t 0 a a t ` (fn ;:::;n gZ) 1 m However, a detailed proof along the historical lines sketched above would run about 30 pages. A natural and important question is whether one could just prove directly that the determinant satisfies the Kolmogorov equation. The purpose of this article is to provide such a proof, for TASEP and its variant PushASEP. We comment that the direct proof is relatively short, but far from obvious. It is not known if it can be obtained directly from the biorthogonal representation (1.3) without knowledge of the special form of the  (x). It is also worth noting that there are examples (such as discrete time TASEP, see [MQR+]) with explicit  (x) for which the present short proof does not work. TASEP has a rich history, some of which can be found in the introduction of [Sch97]. Since that paper there has been significant progress in exact solvability, with much interest stemming from TASEP’s role as one of the fundamental models in the KPZ universality class. There are many results; we mention here only a few besides the ones which fit directly into our story above. The one-point distribution of TASEP with step initial data (X (i) = i, i  1) was first solved in [Joh00] by exploiting a Toeplitz structure which is only available for that choice of initial condition; the main goal of that paper was to prove the now-classic fact that the fluctuations of the position of an appropriately chosen particle converge, as time goes to infinity, to the Tracy-Widom GUE distribution [TW94]. This was later extended in [Joh03] to the multi-point distributions of TASEP, which were shown there to converge to the Airy process [PS02]. Much work was devoted during the last fifteen years to extending this type of results to two other choices of initial data: periodic (X (i) = 2i, i 2 Z, which was the subject of the papers [Sas05; BFPS07] discussed above) and stationary (X corresponding to placing particles on Z according to a product measure) [IS04; BFS08; BFP10]. In the past ten years there has been a huge effort to extend some of these results to the partially asymmetric simple exclusion process (see, for example, [TW08; IS11; BCS14]). These still depend on very particular initial data. In the case of TASEP on a ring there is also a huge literature with some recent breakthroughs: [BL18], where asymptotics is done for a type of step initial data, and [Pro13; Pro14], which actually computes the entire spectrum. SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 4 2. KOLMOGOROV EQUATION FOR TASEP In this section we provide a short self-contained proof of Thm. 1.1. The strategy is to compute the TASEP two sides of (1.10) with F given by (1.12) and K given by (1.8), and check that they are equal. We begin with the right hand side of (1.10). We will first consider the effect of moving a single particle, and then we will sum over all particles to obtain the effect of the generatorL. Fix a particle label k and consider the original initial condition X and the initial condition X where particle k is 0 0 moved to the right by one, i.e. X (k) = X (k) + 1 (for now it does not even matter whether or not this 0 0 particle can actually be moved without violating the strict particle order condition). We will compare TASEP the kernel K with these two initial conditions; for notational convenience we put tildes on top of TASEP all objects when they depend on the modified initial condition (e.g. K refers to the kernel started from the initial condition X ). By (1.8) we have epi(X ) epi(X ) TASEP TASEP  0 0 K K (n ;; n ;) = (S ) S S ; i j t;n t t i t;n t;n j j while by (1.7) we have epi(X ) epi(X ) 0 0 S S (z ; z ) = E S (B ; z )1 S (B ; z )1 : 1 2 B =z t;n ~  ~ 2  ~<n t;n  2 <n 0 1 t;n t;n Recall that  means the hitting time of the strict epigraph of X (` + 1) ;  ~ is the same but with the `0 modified X . It is clear from the definitions that  ~ =  unless  = k 1 < n and B = X (k) + 1 0 0 k1 (see Fig. 1 for a visual explanation), and thus the above expectation equals E S (B ; z )1 S (B ; z )1 1 1 ; B =z t;n ~  ~ 2  ~<n t;n  2 <n  =k1;B =X (k)+1 kn 0 1 0 k1 (n) which by the Markov property we may write as f (z )g (z ) with k 1 2 f (z) = P ( = k 1; B = X (k) + 1) (2.1) k B =z k1 0 and h i (n) g (z) = E S (B ; z)1 S (X (k) + 1; z) 1 ; (2.2) (k) (k) (k) t;nk+1 0 kn B =X (k)+1 k k1 0 t;n ~  ~  ~ <n where the subscript in expectation on the right hand side indicates that the random walk is started now (k) at X (k) + 1 at time k 1 and the superscript in  ~ indicates that the hitting time is now restricted to times infk; k + 1; : : :g (we will use the same notation for other hitting times below). This means that (n ) (k) TASEP TASEP  j := K K (n ;; n ;) = (S ) f g : (2.3) i j t;n k t t i i;j t;k Crucially, this is a rank one kernel (acting on ` (Z)), and then the same is true of the extended kernel TASEP TASEP 2 K K when thought of as acting on the extended space ` (fn ; : : : ; n g Z), so using 1 m t t the fact that det(I + A + B) det(I + A) = det(I + A) tr[(I + A) B] for B any rank one operator, we deduce that TASEP TASEP det(I   K   ) det(I   K   ) a a a a t t TASEP TASEP 1 (k) = det(I   K   ) tr[(I   K   )      ] a a a a a a t t (k) TASEP TASEP (here  = K K as in (2.3)). This difference corresponds to the bracket in the k-th term t t of the sum (1.9) defining the action of L on our function F (X ; a; n). The conclusion then is that if t 0 J  N is the set of labels of particles which can be moved (i.e. those so that X (k 1) X (k) > 1), 0 0 then LF (X ; a; n) equals the sum over k 2 J of the above difference, i.e. t 0 h i TASEP TASEP 1 (k) LF (X ; a; n) = det(I   K   ) tr (I   K   )      : t 0 a a a a a a t t k2J (k) Note that, as above, the sum only has finitely many non-zero terms, because  = 0 for all k > n (as is clear, for instance, from (2.2) and (2.3)). Finally we note if k 2= J (i.e. if X (k) + 1 = X (k 1)), 0 0 then since the random walk B has increments B B  1, we have fB = X (k) + 1g t1 t k1 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 5 X (1) X (2) X (3) X (4) X (4) X (5) 1 2 3 4 F IGURE 1. A sample path of the random walk B (circles) and the graph of the initial n1 n1 conditionsfX (i + 1)g andfX (i + 1)g which differ only at position k with 0 0 i=0 i=0 X (k) = X (k) + 1. The random walk has i.i.d. geometric decrements, and  and  ~ 0 0 are the first times the random walk B hits the strict epigraph of the initial condition curve with X and X respectively. Note that  ~ =  unless the hit happens exactly at 0 0 = k 1 < n and B = X (k) = X (k) + 1. In this example k = 4,  = 3 and k1 0 0 ~ > 4. fB  X (k 1) + 1g  f  k 2g and thus f (z) = 0. Thus summing over j 2= J has no k2 k effect, and we have then h i TASEP TASEP 1 (k) LF (X ; a; n) = det(I   K   ) tr (I   K   )      : (2.4) t 0 a a a a a a t t k1 Now we turn to the left hand side of (1.10). In general, if a kernel K depends smoothly on a d 1 d parameter h one has det(I K ) = det(I K ) tr[(I K ) K ]. Thus h h h h dh dh h i d TASEP TASEP 1 d TASEP F (X ; a; n) = det(I   K   ) tr (I   K   )   K t 0 a a a a a a t t t dt dt and then in view of (2.4) all we need in order to deduce (1.10) is to check that d TASEP (k) K =  : (2.5) t k1 dt We have from (1.8) that epi(X ) epi(X ) d TASEP d  0  d 0 := K (n ;; n ;) = (S ) S + (S ) ( S ): i;j i j t;n t;n t i t;n i t;n dt dt dt From (1.5) we get d  1 1 (S ) = r (S = (S r t;n t;n t;n dt 2 2 d 1 1 and from (1.6) we get S = S r = r S , which means t;n t;n t;n dt 2 2 epi(X ) d 0 1 S (z ; z ) = E r S (B ; z )1 1 2 B =z t;n  2 <n t;n 0 1 dt 2 (we have used here thatr commutes with (S ) andS since they all have Toeplitz kernels). t;n t;n This gives = (S ) H i;j t;n n i j with epi(X ) H (z ; z ) = r S (z ; z ) + E r S (B ; z )1 n 1 2 1 2 B =z t;n  2 <n t;n 0 1 epi(X ) = S (z 1; z ) E S (B 1; z )1 : 1 2 B =z t;n  2 <n t;n 0 1 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 6 Shifting the curve defined by X up by one, the first term on the second line can be expressed as E S (B 1; z )1 with  ^ the hitting time of the strict epigraph of X (m+1)+1 , B =z t;n ^  ^ 2  ^<n 0 0 1 m0 and thus H (z ; z ) = E S (B 1; z )1 S (B 1; z )1 : n 1 2 B =z t;n ^  ^ 2  ^<n t;n  2 <n 0 1 Note that  ^   , so the difference inside the brackets vanishes for   n, while if  < n and B  X ( + 1) + 2 we have  ^ =  , so the difference vanishes again. We deduce that n1 H (z ; z ) = P ( = k; B = X (k + 1) + 1) n 1 2 B =z k 0 0 1 k=0 h h i i E S (B 1; z )1 S (X (k + 1); z ) (k) (k) (k) 2 t;nk 0 2 B =X (k+1)+1 t;n ^  ^  ^ <n k 0 n1 (n) = f (z )gb (z ) k+1 1 2 k+1 k=0 (n) with f as in (2.1) and gb (z) given by the bracket in the middle line (see also the comment after k+1 (2.2)). Shifting now the curve defining  ^ back to X in the last expectation we get (n) gb (z) = E S (B ; z)1 S (X (k); z): (2.6) B =X (k) t;n  <n t;nk+1 0 k k1 0 We have then n (n ) d TASEP 1 j  j K (n ;; n ;) = (S ) f (z ) gb (z ): i j t;n k 1 2 t i dt 2 k=1 k (n) Comparing this formula and (2.3), and since g  0 for k  n by its definition (2.2), all that remains in order to prove (2.5) is that (n) (n) gb = 2g for all k  n: (2.7) k k t t r (n) r 2 2 But from (1.6) we have S e = S = Q , so applying e on the right in (2.7) and t;n 0;n (n) (n) using the definitions (2.2) and (2.6) of g and gb we see that (2.7) is equivalent to k k h i (k) (n ) (nk+1) E Q (B (k); z)1 (k) Q (X (k); z) B =X (k) k1 0   <n h i (k) (n ~ ) (nk+1) = 2E Q (B (k); z)1 (k) 2Q (X (k) + 1; z) (2.8) B =X (k)+1 k1 0  ~  ~ <n (n) z for k  n. By definition of Q , if we multiply the equation by 2 then both sides are polynomials in z, so it is enough to prove the equality for z < X (n). But it is easy to check that (n) n Q (x; y) = Q (x; y) for all x y  1: (2.9) (see [MQR17, Eqn. (2.23)]), so for such z, since X (i) is decreasing in i then z is also smaller than B , (`) ` B and X (k) in (2.8), and we deduce that all the Q ’s in the identity can be replaced by Q . As a ~ 0 consequence, what we need to prove is that h i (k) n nk+1 E Q (B ; z)1 Q (X (k); z) (k) (k) B =X (k) 0 k1 0   <n h i (k) n ~ nk+1 = 2E Q (B (k); z)1 (k) 2Q (X (k) + 1; z); B =X (k)+1 k1 0  ~  ~ <n or simply (k) (k) P B = z;   n = 2P B = z;  ~  n ; n n B =X (k) B =X (k)+1 k1 0 k1 0 which is easy to see since the walk takes Geom[ ] steps to the left: in fact, since X (k + 1) < X (k) 0 0 (k) (k) and  and  ~ may only differ on the event that the walk hits at time k 1, then for k < n, SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 7 decomposing according to the first step, (k) (k) P B = z;  ~  n = Q(X (k) + 1; y)P B = z;  ~  n n 0 B =y n B =X (k)+1 k1 0 yX (k+1) k 1 (k) 1 (k) = Q(X (k); y)P B = z;   n = P B = z;   n ; 0 B =y n B =X (k) n yX (k+1) k 0 0 2 2 k1 the case k = n is even simpler. This yields the desired equality (2.7). 3. INITIAL CONDITION Now we check that the initial condition (1.11) is satisfied. It is simpler to deal first with the one-point case (m = 1), so we start there. 3.1. One-point case. Let (n) TASEP K = K (n;; n;): t t We need to prove that (n) det I   K   = 1 : (3.1) a a X (n)>a 0 0 ` (Z) Remark 3.1. In the one point case, the initial condition (3.1) can be checked rather easily from the biorthogonal ensemble representation (1.3), as mentioned by one of the referees; the proof using this biorthogonalization is given below. However, since our goal in this paper is to prove (1.12) directly from the representation (1.8) (and thus avoid having to use the result of [MQR17] to connect this formula back to (1.3)), we provide the direct (but relatively more complicated) proof from (1.8) in this subsection. (n) n1 n n From (2.3) we have K (z ; z ) = (z ) (z ) with t = 0 in the definition (1.4). Let 1 2 1 2 0 k=0 k k n n ;  : Z  f0; : : : ; n 1g ! R be given by (z; k) = (z) and (z; k) =  , so that the left k k hand side of (3.1) equals det(I      ) = det(I    ); a a a ` (Z) where we have used the cyclic property of the Fredholm determinant. Note that on the determinant on the right hand side is finite dimensional, with the n n matrix inside being given by n n (`; k) = (z) (z); `; k = 0; : : : ; n 1: k ` za Assume first that a < X (n). Then a < X (n k) k for all k = 0; : : : ; n, which implies from (1.4) 0 0 that (z) = 0 for all z  a and all such k, so det(I    ) = det(I ) = 1 as desired. Assume otherwise that a  X (n). From (1.4) it is straightforward to check that (z) = 1 and 0 z=X (n) n n n thus the first column of    is given by    (`; 0) = (z) (z) =  (X (n)). But a a 0 za k ` ` n n n the biorthogonalization condition (z) (z) = 1 implies that  (X (n)) = 1 . We `=0 0 `=0 z2Z ` ` deduce that the first column of I    equals zero, so det(I    ) vanishes as desired. a a We turn now to the proof of (3.1) directly from (1.8). Consider first the case a  X (n) 1, which (n) implies that in (3.1) we need to evaluate K (z ; z ) only for z  X (n) 1. By definition of the 1 2 1 0 epi(X ) (n) hitting time  in (1.7) we haveS (z ; z ) = Q (z ; z ) for z > X (1), so we may write 1 2 1 2 1 0 0;n (n) epi(X ) n (n) n 0 K = Q  Q + Q   S : (3.2) X (1) X (1) 0 0 0 0;n Now n yx+n yx n Q (x; y) = (1) 2 1 (3.3) 0yxn yx so, in particular, Q (x; y) = 0 for y x > n. If y > X (1) and z  X (n) 1 we have 0 1 0 y z > X (1) X (n) + 1  n, so Q (z ; y) = 0 and the first term on the right hand side of (3.2) 1 0 0 1 epi(X ) 0 n1 vanishes. For the second term we consider S (y; z ) = E S (B ; z )1 with 2 B =y 0;nk k 2  =k 0;n k=0 y  X (1) and note that  = k implies B > X (k + 1), which cannot happen unless the starting 0 k 0 location satisfies y > X (k + 1) + k (because the walk takes downward steps of size at least 1). But 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 8 for such y (and z  X (n) 1) we have y z > X (k + 1) + k X (n) + 1  n, and thus again 1 0 1 0 0 we get Q (z ; y) = 0. We have proved that (n) K = 0 when a < X (n); (3.4) a 0 and then the left hand side of (3.1) equals det(I ) = 1 for a < X (n), as desired. Next we turn to the case a  X (n), and focus on the column corresponding to the index X (n) in 0 0 (n) (n) K   , that is K (z; X (n)) . We claim that, in this case, a a 0 0 0 za (n) K (z; X (n)) = 1 for all z 2 Z: (3.5) z=X (n) 0 0 (n) If that is the case then the matrix I   K   has a column which is identically 0, and thus its a a determinant vanishes as required. In order to prove (3.5), observe first that the argument we used to prove (3.4) actually shows that (n) K (z; X (n)) = 0 for all z < X (n). Consider next the case z = X (n). For the first term in (3.2), 0 0 0 n (n) a simple residue computation using the contour integral formulas for Q and Q which follow from (1.5) and (1.6) with t = 0 gives (remembering X (n) X (1)  1 n) 0 0 H H n X (n)X (1)+n1 0 0 (1w) (1v) n (n) 1 1 Q  Q (X (n); X (n)) = dw dv 0 0 2 X (1) X (n)X (1)+n n 0 (2i) 1vw w 0 0 v (3.6) = 1 : X (n)X (1)=1n 0 0 (n) Since we want to show that K (X (n); X (n)) = 1, then from (3.6) together with (3.2) we see that 0 0 we need to show that epi(X ) n 0 Q   S (X (n); X (n)) = 1 : (3.7) 0 0 X (1) X (n)X (1)<1n 0 0;n 0 0 We consider the two possibilities for X (n) X (1) separately. Assume first that X (1) X (n) = 0 0 0 0 n 1. This means that the first n particles are packed one next to the other, i.e. X (k) = X (1)k + 1, 0 0 and thus, since the walk starts below X (1) and takes downward steps of size at least 1, it cannot hit epi(X ) the strict epigraph of the curve by time n, so   S = 0 as desired. The other possibility X (1) 0 0;n epi(X ) (n) is that X (1) X (n) > n 1. Using S = Q we may express the kernel S (y; X (n)) 0 0 0;n 0 0;n n1 (nk) as E Q (B ; X (n))1 and then note that, since inside the expectation B is the B =y k 0  =k k k=0 location where the strict epigraph is hit, B X (n)  X (k + 1) + 1 X (n)  n k  1, k 0 0 0 n1 (nk) nk so by (2.9) we may replace Q by Q there to get P (B = X (n);  = k). But B =y n 0 k=0 B = X (n) implies B  X (n) + 1, which in turn implies   n 1, so the last sum equals n 0 n1 0 simply Q (y; X (n)), and then epi(X ) n n n Q   S (X (n); X (n)) = Q (X (n); y)Q (y; X (n)): 0 0 0 0 X (1) 0 0;n yX (1) And the last sum can be extended to all y’s because if y > X (1) then, under our assumption X (n) X (1) < 1 n and using (3.3), the first factor in the summand vanishes. This shows that, in 0 0 epi(X ) this case, Q   S (X (n); X (n)) = 1 as desired. X (1) 0 0 0 0;n Up to here we have proved (3.5) for z  X (n). Consider finally the case z > X (n). The exact 0 0 same argument used in the last paragraph to prove (3.7) also shows epi(X ) n 0 n (n) Q   S (z; X (n)) = Q  Q (z; X (n))1 : 0 0 X (1) X (1) X (1)X (n)>n1 0 0;n 0 0 0 (n) In particular, in the case X (1) X (n) > n 1 we get using (3.2) that K (z; X (n)) = 0 0 0 n (n) n (n) Q Q (z; X (n)) which equals 0 because in fact Q Q  0 (see [MQR17, Eqn. (2.24)]). It remains to prove that the first term in (3.2) vanishes when X (1) X (n) = n 1, but this is 0 0 straightforward, because in this case we have (using the contour integral formulas (1.5) and (1.6) as before) H H n H (1w) 1 1 1 1 Q  S (z; X (n)) = dw dv = dw = 0 0;n 0 2 X (1) zX (1)+n n zX (1)+n 0 1vw 2i (2i) 0 0 w v w since z X (1) + n > X (n) X (1) + n = 1. 0 0 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 9 TASEP TASEP 3.2. Multi-point case. We let L = K (n ;; n ;) and think of K as an operator-valued i;j i j matrix L with entries L . The initial condition (1.11) will follow from the arguments in the one-point i;j case and the identity (n ) n n n n i j i j i L = Q 1 + Q K ; (3.8) i;j n <n i j t which is (D.8) in [MQR17]. Consider first the case a  X (n ) 1 for each i. Here (3.4) gives i 0 i (n ) (n ) (n ) n n i n n i i j i j i K = 0, and then in the sum   Q K (z ; z ) = 1 Q (z ; y)K (y; z ) a a 1 2 z a 1 2 i i t 1 i y2Z 0 only terms with y > a survive. On the other hand, we have Q (z; y) = 0 for z < y (any m 2 Z) so n n n n j i j i Q (z ; y) = 0 for z  a < y. This shows that   L   = Q 1 , and since this 1 1 i a i;j a n <n i j i j holds for all i, we get F (X ; a; n) = det(I  L  ) with L a strictly upper triangular operator-valued 0 0 a a matrix, which yields F (X ; a; n) = 1 as desired. 0 0 Suppose next that a  X (n ) for some j, so we need to prove det(I   L  ) = 0. We may j 0 j a a assume without loss of generality that n is the largest particle label amongfn ; : : : ; n g for which j 1 m this inequality holds. We will focus on the j-th column of   L  , and more precisely on the subcolumn a a (j) (j) of this column of kernels corresponding to the index X (n ), that is v := v (z) with 0 j ` z2Z; i=1;:::;m (j) (n ) v (z) =   L (z; X (n )). In this case (3.5) gives us K (z; X (n )) = 1 for all z, and a i;j 0 j 0 j z=X (n ) i 0 0 ` j (n ) n n n n j i j i then Q K (z; X (n )) = Q (z; X (n )), so that using (3.8) we get 0 j 0 j (j) n n j i v (z) = Q (z; X (n ))1 1 : 0 j n n za i i j i (j) The case i = j is direct and yields v (z) = 1 . Otherwise, if n > n , then necessarily a < i j i z=X (n ) j 0 j X (n ) by our choice of j, and hence for z  a we have X (n )z > X (n )X (n )  n n 1 0 i i 0 j 0 j 0 i i j (j) (j) so as before v (z) = 0. The conclusion is that v (z) = 1 1 . But then I   L  has a i=j z=X (n ) a a i i 0 j (j) column (namely v ) which is identically 0, and thus F (X ; a; n) = det(I   L  ) = 0 as required. 0 0 a a 4. PUSHASEP We consider now PushASEP, a generalization of TASEP and the Toom model [DLSS91] introduced in [BF08]. Again we have particles on the lattice Z with positions X (1) > X (2) >  . Each particle t t attempts to jump one step to the right at rate r, with jumps being permitted only if the neighboring site is empty; this is the same TASEP dynamics considered above (except run at rate r). On the other hand, each particle has another (independent) exponential clock running at rate `, and when it rings the particle jumps to the nearest vacant site on its left. Relabeling the particles in order to keep the ordering after such a jump, we may think of the effect of a particle jumping left as pushing all its left neighbors one step to the left. TASEP and the Toom model (or PushTASEP) are recovered from PushASEP by setting ` = 0 and r = 0, respectively. The generatorL of this process can be written similary to the one for TASEP: letting b(k) denote r;` the length of the block of nearest neighbor occupied sites lying to the left of (and including) particle k, meaning that X (k + j) = X (k) j for j = 0; : : : ; b(k) 1 and X (m + b(k)) < X (k) b(k), 0 0 0 0 L f (X ) = r 1 f (: : : ; X (k) + 1; : : : ) f (: : : ; X (k); : : : ) r;` X (k1)X (k)>1 k1 + ` f (: : : ; X (k) 1; : : : ; X (k + b(k)) 1; : : : ) f (: : : ; X (k); : : : ; X (k + b(k)); : : : ) ; k1 where again we take X (0) = 1. As for the TASEP case we define F (X ; a; n) = P X (n ) > a ; : : : ; X (n ) > a ; t X t 1 1 t m m again we have that F (; a) satisfies the Kolmogorov backward equation F = L F; (4.1) r;` dt with initial condition F (X ; a; n) = 1 : (4.2) 0 X (1)>a ;:::X (n)>a 1 n SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 10 Theorem 4.1. The unique solution of (4.1), (4.2) is given by PushASEP F (X ; a; n) = det(I   K   ) 2 t 0 a a t ` (fn ;:::;n gZ) 1 m with epi(X ) PushASEP n n  0 j i K (n ;; n ;) = Q 1 + (S ) S ; i j n <n t;n i j i t;n where the operatorsS and S are now given by t;n t;n 1 + rtr +2`tr n S (z ; z ) = (e Q ) (z ; z ) (4.3) t;n 1 2 1 2 1 (1 w) t[r(w1=2)+`(1=w2)] = dw e ; z z n+1+z z 2 1 2 1 2i 2 w 1 + (n) rtr 2`tr S (z ; z ) = Q e (z ; z ) (4.4) t;n 1 2 1 2 z z +n1 2 1 1 (1 w) t[r(w1=2)+`(21=(1w))] = dw e ; z z n 1 2 2i 2 w 1 t rtr +2`tr r 2 2 herer f (x) = f (x + 1) f (x) and e is defined similarly to e after (1.6). This result, which is new, can be proved along the same lines as the solution of TASEP given in [MQR17]; such a proof will appear in [MQR+] (where convergence to the KPZ fixed point will also be proved). Our goal here, as for the Kolmorogov equation for TASEP in Thm. 1.1, is to prove directly that the Fredholm determinant provides a solution to the Kolmogorov backward equation (4.1), (4.2). Note first that for t = 0 the formulas forS andS given in Thm. 4.1 are the same as those t;n t;n for TASEP, so the initial condition (4.2) follows from the TASEP proof. Next we claim that in order to prove that the Fredholm determinant solves (4.1) it is enough to consider the case of PushTASEP, i.e. r = 0, since the TASEP part of the dynamics was already handled in Sec. 2. To see this it is convenient PushASEP to define F (X ; a) = det(I   K   ) 2 with t ;t 0 a a 1 2 t ;t ` (fn ;:::;n gZ) 1 2 1 m epi(X ) PushASEP n n  0 j i K (n ;; n ;) = Q 1 + (S ) S ; i j n <n t ;t ;n t ;t i j 1 2 i t ;t ;n 1 2 1 2 j 1 + 1 + epi(X ) rt r +2`t r n  (n) rt r 2`t r 0 1 2 1 2 2 2 S = (e Q ) , S = Q e , and S defined t ;t ;n t ;t ;n 1 2 1 2 t ;t ;n 1 2 analogously using S , so that (4.1) will follow (after setting t = t = t) if we prove that t ;t ;n 1 2 1 2 j @ @ ( + )F (X ; a) = L F (X ; a). But since L = L +L , it is actually enough to t ;t 0 `;r t ;t 0 `;r `;0 0;r 1 2 1 2 @t @t 1 2 prove that @ @ F (X ; a) = L F (X ; a) and F (X ; a) = L F (X ; a): (4.5) t ;t 0 0;r t ;t 0 t ;t 0 t ;t 0 `;0 1 2 1 2 1 2 1 2 @t @t 1 2 The first of these is a version of the Kolmogorov equation for TASEP (1.1) where in the contour integrals in (1.5) and (1.6) there are respectively additional factors (w) and (w) inside the integrands, both independent of t, and it is not hard to see that the proof in Sec. 2 works without any difficulty after this modification. We will prove below that (4.1) holds in the PushTASEP case r = 0, and inspecting the proof it is easy again to seen that this implies the second equation in (4.5). We conclude that it suffices to check the case r = 0, and we can also take ` = 1 as it is just a scaling factor. Given this reduction, one may naturally wonder why the proof works for PushASEP and not ASEP. It just turns out that when one tries to adapt Schütz’s formula [Sch97] to include jumps to the left (for example, as done in [BF08]) the resulting formula ends up solving the push dynamics instead of the leftward TASEP dynamics. Having made these comments, we begin the proof. Proceeding as for TASEP, let k 2 f1; : : : ; ng and consider the initial condition X and the one where particle k is moved one step to the left along with its neighboring block. Write b for the length of this block. As for TASEP, we need to compute epi(X ) epi(X ) 0 0 S S (z ; z ) = E S (B ; z )1 S (B ; z )1 ; 1 2 B =z t;n ~  ~ 2  ~<n t;n  2 <n t;n t;n 0 1 where now the tildes refer to the initial condition modified in this new way. We clearly have  =  ~ unless ~ 2 fk1; : : : ; k +b2g and B = X ( +1). But since the particles X (j) for j 2 fk; : : : ; k +b1g ~ 0 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 11 lie in a block,  ~ 2 fk 1; : : : ; k + b 2g implies  ~ = k 1 (because the walk always jumps down by at least one step). This means that the above expectation equals (n) E S (B ; z )1 S (B ; z )1 1 = f (z )g (z ) B =z t;n ~  ~ 2  ~<n t;n  2 <n  ~=k1;B =X (k) k 1 2 0 1 0 k1 k with f (z) = P ( ~ = k 1; B = X (k)) k B =z k1 0 (k) and (see also the comment after (2.2) for the meaning of  ) h i (n) g (z) = S (X (k); z) E S (k) (B (k); z)1 (k) 1 t;nk+1 0 kn B =X (k) k k1 0 t;n   <n (n) (k) (note that f does not depend on b but g does, through  ). Since this is a rank one kernel, the same arguments as in the TASEP case give h i PushTASEP PushTASEP 1 (k) L F (X ; a; n) = det(I   K   ) tr (I   K   ) 0;1 t 0 a a a a a a t t k=1 (n) (k) where  is defined exactly as in (2.3) but using the new versions of f and g (note that here L k 0;1 acts on all particles, since all particles are allowed to move in the PushTASEP dynamics). Now we consider the left hand side. As for TASEP, all we need to check in order to finish the proof is that d PushTASEP (k) K =  : t k1 dt Now we have epi(X ) epi(X ) d PushTASEP d  0  d 0 = K (n ;; n ;) = (S ) S + (S ) ( S ); i;j i j t;n t;n t i t;n i t;n dt dt dt j j d  + d  + and using (4.3) and (4.4) we have (S ) = 2(S r and (S ) = 2r (S . t;n t;n t;n t;n dt dt Therefore = 2S H i;j t;n n i j with epi(X ) + 0 + H (z ; z ) = r S (z ; z ) E r S (B ; z )1 n 1 2 1 2 B =z t;n  2 <n 0 1 t;n epi(X ) = S (z + 1; z ) E S (B + 1; z )1 : 1 2 B =z t;n  2 <n t;n 0 1 The first term on the second line can be expressed as E S (B + 1; z )1 with  ^ the B =z t;n ^  ^ 2  ^<n 0 1 hitting time of the strict epigraph of X (m + 1) 1 , and thus m0 H (z ; z ) = E S (B + 1; z )1 S (B + 1; z )1 : n 1 2 B =z 2 t;n  2 <n t;n ^  ^  ^<n 0 1 Now we have  ^   , so the difference inside the brackets vanishes for  ^  n, and it also vanishes when ^ < n and B  X ( ^ + 1) + 1. Thus H (z ; z ) equals ^ 0 n 1 2 n1 P ( ^ = k; B = X (k 1)) B =z k 0 0 1 k=0 h i (k) S (X (k + 1) + 1; z ) E S (B + 1; z )1 (k) t;nk 0 2 t;n 2 B =X (k+1) k 0  <n n1 (n) = f (z )gb (z ) k+1 1 2 k+1 k=0 with f (z) = P ( ^ = k 1; B = X (k)) k B =z k1 0 and h i (n) gb (z) = S (X (k) + 1; z) E S (B + 1; z)1 : (k) (k) (k) t;nk+1 0 B =X (k) t;n   <n k k1 0 SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 12 We have then (n ) = 2 S f (z ) gb (z ); i;j t;n k 1 2 k=1 (n) (n) and so it remains to prove that f = f for all k and gb = g for k  n (using as for TASEP that k k k 2 k (n) g  0 for k > n). For fixed k let  ~ denote the hitting time of the initial condition X introduced above where the block starting at X (k) is moved one step to the left, so that f (z) = P ( ~ = k 1; B = X (k)). 0 k B =z k k1 0 We obviously have  ^   ~ , and thus  ^ = k 1 and B = X (k) implies  ~ = k 1, so that k k1 0 k f (z) = P ( ^ =  ~ = k 1; B = X (k)). This gives k B =z k k1 0 f (z) f (z) = P ( ~ >  ^ = k 1; B = X (k)) = 0: k k B =z k k1 0 (n) (n) 1 2tr Next we prove the identity gb = g for k  n. Proceeding as for TASEP, we apply e on the k 2 k right of this identity to see that it is equivalent to h i (k) (nk+1) (n ) Q (X (k) + 1; z) E Q (B + 1; z)1 (k) (k) 0 B =X (k) k1 0   <n h i (k) 1 (nk+1) 1 (n ) = Q (X (k); z) E Q (B (k); z)1 (k) : B =X (k) k1 0   <n 2 2 Following again the TASEP argument, the two sides are of the form 2 times a polynomial in z, so it is enough to prove the equality for z < X (n), in which case the identity becomes (k) 1 (k) P B = z 1;   n = P B = z;   n ; n n B =X (k) B =X (k) k1 0 2 k1 0 which is easy to prove again since the walk takes Geom[ ] steps to the left and z  X (n). Acknowledgements. MN and JQ were supported by the Natural Sciences and Engineering Research Council of Canada. DR was supported by Programa Iniciativa Científica Milenio grant number NC120062 through Nucleus Millenium Stochastic Models of Complex and Disordered Systems, by Conicyt Basal-CMM Proyecto/Grant PAI AFB-170001, and by Fondecyt Grant 1160174. R EFERENCES [BFP10] J. Baik, P. L. Ferrari, and S. Péché. Limit process of stationary TASEP near the characteristic line. Comm. Pure Appl. Math. 63.8 (2010), pp. 1017–1070. [BL18] J. Baik and Z. Liu. Fluctuations of TASEP on a ring in relaxation time scale. Comm. Pure Appl. Math. 71.4 (2018), pp. 747–813. [BCS14] A. Borodin, I. Corwin, and T. Sasamoto. From duality to determinants for q-TASEP and ASEP. Ann. Probab. 42.6 (2014), pp. 2314–2382. [BF08] A. Borodin and P. L. Ferrari. Large time asymptotics of growth models on space-like paths. I. PushASEP. Electron. J. Probab. 13 (2008), no. 50, 1380–1418. [BFPS07] A. Borodin, P. L. Ferrari, M. Prähofer, and T. Sasamoto. Fluctuation properties of the TASEP with periodic initial configuration. J. Stat. Phys. 129.5-6 (2007), pp. 1055–1080. [BFS08] A. Borodin, P. L. Ferrari, and T. Sasamoto. Transition between Airy and Airy processes 1 2 and TASEP fluctuations. Comm. Pure Appl. Math. 61.11 (2008), pp. 1603–1629. [DLSS91] B. Derrida, J. L. Lebowitz, E. R. Speer, and H. Spohn. Dynamics of an anchored Toom interface. J. Phys. A 24.20 (1991), pp. 4805–4834. [EM98] B. Eynard and M. L. Mehta. Matrices coupled in a chain. I. Eigenvalue correlations. J. Phys. A 31.19 (1998), pp. 4449–4456. [IS04] T. Imamura and T. Sasamoto. Fluctuations of the one-dimensional polynuclear growth model with external sources. Nuclear Phys. B 699.3 (2004), pp. 503–544. [IS11] T. Imamura and T. Sasamoto. Current moments of 1D ASEP by duality. J. Stat. Phys. 142.5 (2011), pp. 919–930. [Joh00] K. Johansson. Shape fluctuations and random matrices. Comm. Math. Phys. 209.2 (2000), pp. 437–476. SOLUTION OF THE KOLMOGOROV EQUATION FOR TASEP 13 [Joh03] K. Johansson. Discrete polynuclear growth and determinantal processes. Comm. Math. Phys. 242.1-2 (2003), pp. 277–329. [Lig10] T. M. Liggett. Continuous time Markov processes. Vol. 113. Graduate Studies in Mathemat- ics. An introduction. American Mathematical Society, Providence, RI, 2010, pp. xii+271. [MQ17] K. Matetski and J. Quastel. From the totally asymmetric simple exclusion process to the KPZ fixed point. 2017. arXiv: 1710.02635[math.PR]. [MQR17] K. Matetski, J. Quastel, and D. Remenik. The KPZ fixed point. 2017. arXiv: 1701 . 00018[math.PR]. [MQR+] K. Matetski, J. Quastel, and D. Remenik. TASEP and generalizations: Method for exact solution. In preparation. 2019. [PS02] M. Prähofer and H. Spohn. Scale invariance of the PNG droplet and the Airy process. J. Stat. Phys. 108.5-6 (2002), pp. 1071–1106. [Pro13] S. Prolhac. Spectrum of the totally asymmetric simple exclusion process on a periodic lattice—bulk eigenvalues. J. Phys. A 46.41 (2013), pp. 415001, 36. [Pro14] S. Prolhac. Spectrum of the totally asymmetric simple exclusion process on a periodic lattice-first excited states. J. Phys. A 47.37 (2014), pp. 375001, 29. [Sas05] T. Sasamoto. Spatial correlations of the 1D KPZ surface on a flat substrate. Journal of Physics A: Mathematical and General 38.33 (2005), p. L549. [Sch97] G. M. Schütz. Exact solution of the master equation for the asymmetric exclusion process. J. Statist. Phys. 88.1-2 (1997), pp. 427–445. [TW94] C. A. Tracy and H. Widom. Level-spacing distributions and the Airy kernel. Comm. Math. Phys. 159.1 (1994), pp. 151–174. [TW08] C. A. Tracy and H. Widom. A Fredholm determinant representation in ASEP. J. Stat. Phys. 132.2 (2008), pp. 291–300. (M. Nica) DEPARTMENT OF M ATHEM ATICS, UNIVERSITY OF TORONTO, 40 S T. G EORGE STREET, TORONTO, ONTARIO, CANADA M5S 2E4 E-mail address: mnica@math.toronto.edu (J. Quastel) DEPARTMENT OF M ATHEM ATICS, UNIVERSITY OF TORONTO, 40 S T. GEORGE STREET, TORONTO, ONTARIO, CANADA M5S 2E4 E-mail address: quastel@math.toronto.edu (D. Remenik) DEPARTAMENTO DE I NGENIERÍA M ATEMÁTICA AND C ENTRO DE M ODELAMIENTO MATEMÁTICO (UMI-CNRS 2807), UNIVERSIDAD DE C HILE, AV. BEAUCHEF 851, TORRE NORTE, PISO 5, SANTIAGO, CHILE E-mail address: dremenik@dim.uchile.cl

Journal

Mathematical Physics – arXiv (Cornell University)

Published: Jun 4, 2019

There are no references for this article.