dynamics of a nonequilibrium system is one of the most stimulating topics in statistical mechanics .
the dynamical property is determined by the interplay among various elements such as interaction , initial and boundary conditions , and so on . but little is known about how it depends on these elements and realizes universality in a situation where the elements are intertwined .
the aim of this article is to approach such questions by analyzing the dynamics of the one - dimensional asymmetric simple exclusion process ( asep ) .
the asep is one of the most typical models of interacting particle processes where particles diffuse to a preferred direction with hard core repulsive interaction  @xcite . in spite of the simplicity of the model
, it has been known that the asep shows various interesting phenomena caused by the collaboration between the diffusion and exclusion effect .
furthermore , in the one - dimensional case , it has an integrable mathematical structure which allows us to analyze some physical quantities exactly .
for example , in the open boundary condition a steady state can be constructed using the matrix product method and @xmath0-orthogonal polynomials , and the boundary - induced phase transition can be discussed  @xcite .    on the other hand , the dynamical ( non - stationary ) properties of the one - dimensional asep
have been also analyzed recently . among various studies on the topic , we focus on the problem of the diffusion of a particular particle in the one - dimensional asep .
( we call the particle a _ tagged _ particle . )
this is a fundamental and elementary problem for the understanding of the dynamics of the asep .
we are interested in how initial conditions and the exclusion effect through hard - core repulsive interaction affect the diffusion property and how it is different from the ( normal ) diffusion of the brownian particle .
the scaling exponent of the diffusion provides insights into these problems .
it was found that it depends on the initial conditions . in a steady state with a given density , the position fluctuation of a tagged particle grows as @xmath1  @xcite .
the exponent 1/2 is the same as that of the brownian particle . on the other hand , for the fixed initial conditions where initially we create particles with a given density and fix the initial configuration , the exponent changes from 1/2 to 1/3  @xcite .
such a diffusion is called an anomalous diffusion .
this exponent 1/3 can be understood from the fact that the time evolution of the one - dimensional asep can be described by the one - dimensional kardar - parisi - zhang universality class  @xcite . in the study of the current in the one - dimensional asep on a ring ,
the exponent 1/3 was calculated by the bethe ansatz technique  @xcite .
note that the two cases above have a common property that the density is invariant under both time and space translation and a tagged particle moves with a fixed velocity on average . on the other hand , let us now consider another initial condition where all sites on the left of a certain bond are occupied while all right sites are empty .
we call this initial condition a _ step _ initial condition .
this is the typical initial condition where the density and the average velocity of a tagged particle change with time . under the initial condition
, we focus on the @xmath2th particle from the right as a tagged particle and consider the time evolution of the particle when @xmath2 is large ( but finite ) . in order to study this problem
, we investigate the limiting behavior of the diffusion of the tagged particle as both @xmath2 and time go to infinity . in the totally asymmetric simple exclusion process ( tasep ) , where particles move only to the right , the average position of a tagged particle has been given in the scaling limit  @xcite .
on the other hand , we are interested in the limiting process itself which the tagged particle obeys , as well as the average position and the scaling exponent .
such an attempt as getting more sophisticated information than the scaling exponent has been recently discussed in the study on the current fluctuations in the tasep . not only the exponent but also the scaling function were obtained and it was revealed that the function is equivalent to the largest eigenvalue distribution in the random matrix theory with various universality classes .
this development is based on the results on the longest increasing subsequence in a random permutation  @xcite .
remarkable is that the scaling function detects a difference in initial conditions even if the scaling exponent does not .
for instance , in the case of the step initial condition , the scaling function is equivalent to the gue tracy - widom distribution  @xcite , which is the largest eigenvalue distribution of the gaussian unitary ensemble ( gue )  @xcite . in other initial conditions , on the other hand ,
these are described as the largest eigenvalue distributions with other universality classes in rmt such as gaussian orthogonal ensemble ( goe )  @xcite , two independent goe s which is denoted by goe@xmath3  @xcite , and so on  @xcite .
in addition , the equal - time multipoint distribution of the current fluctuations is obtained in the step and periodic initial conditions  @xcite by use of the green function of the tasep  @xcite .    in this article
, we calculate the multi - time distribution function of position fluctuations of a tagged particle in the tasep with the step initial condition based on the techniques of the directed polymer problem in a 01 matrix  @xcite and stochastic growth of young diagram characterized by the schur process  @xcite .
we express the function as the fredholm determinant and discuss the scaling limit by the saddle point analysis .
we remark that what will be considered in this paper is the pointwise convergence of the kernel .
the convergence of the fredholm determinant itself is certainly expected to hold from previous works on related problems @xcite , but more elaborate asymptotic analysis is necessary for its proof . in the study on the asep ,
such correlations between different times have not been studied yet although the green function  @xcite and the equal - time joint distributions of particle positions  @xcite have been studied recently . in the one - time case ,
on the other hand , the position fluctuation is essentially the same as the current fluctuation which has already been studied in  @xcite .
but it has not been stated as the position fluctuation .
furthermore , we also consider the dependence of the distribution function on defect particles located in front of the tagged particle . here
defect particle means a particle with different hopping rate from other normal particles and we assume the number of them is finite .
such a situation as the tasep with particle - dependent hopping rates has been recently discussed in  @xcite .
we clarify that when their hopping rate is smaller than that of normal particles , the limiting behavior of the multi - time distribution changes due to their presence although their number is finite .    the paper is organized as follows . in the next section ,
we describe the definition of the model .
the main result of this article is given in section  [ main ] .
in section  [ ftts ] , we express the multi - time distribution of position fluctuations in the tasep as fredholm determinant . for this purpose , we first show that the directed polymer problem on a 01 matrix is related to the time evolution of a tagged particle in the tasep with the step initial condition .
next , by applying the dual robinson - schensted - knuth algorithm to the 01 matrix , we map the tagged particle problem in the tasep to the stochastic growth of young diagram described by the schur process . in section  [ asymp ]
, we discuss the two types of scaling limit for the multi - time distribution function , the case where both time @xmath4 and the label of the tagged particle @xmath2 tend to infinity , and the one where @xmath4 goes to infinity with @xmath2 fixed .
some discussions and concluding remark are given in section  [ discussion ] and section  [ conclusion ] respectively .
in this section , we define the model and the quantity which we study in this article .
let us consider the one - dimensional infinite lattice and particles as illustrated in fig .  1 .
each site can be occupied by at most one particle .
suppose all particles are labeled an integer @xmath5 from the right
. a particle moves stochastically obeying the following rules . during each time step between @xmath6 and @xmath7 ,
the particle labeled @xmath5 hops to the right neighboring site with probability @xmath8 ( @xmath9 ) and stays at the same site with probability @xmath10 if the right neighboring site is empty . on the other hand , if the site is occupied by the particle labeled @xmath11 , the @xmath5th particle stays at the same site with probability 1 .
this incorporates the exclusion effect which describes the hard core repulsive interaction .
when hopping rates of all particles are the same , i.e. , @xmath12 , the system is the usual totally asymmetric simple exclusion process ( tasep ) with parallel update .
we allow the particle - dependent hopping rates and study the effects of defect particles .    in this article , we consider the step initial condition in which all sites on the left of a certain bond between two sites are occupied by particles and all sites on the right of the bond are empty .
1(b ) illustrates the initial condition .    under these settings , we label the particles at @xmath13 from the right ( see fig .  1(b ) ) and consider the dynamics of the particle labeled @xmath2 .
( this particle which we focus on is called the tagged particle . ) for this purpose , we set the position coordinate such that the tagged particle is at the origin at time 0 as depicted in fig .
let us define @xmath14 as the position of the @xmath2th particle at time @xmath4 .
in other words , @xmath14 represents the distance travelled by the @xmath2th particle from time @xmath15 through @xmath4 .
the quantity which we will investigate is the multi - time joint distribution function for @xmath14 defined as @xmath16 in particular , we are interested in the asymptotic behavior of   as @xmath4 goes to infinity . in this article , two types of scaling limit will be discussed .
the first one is the limit where both time @xmath4 and the label of the tagged particle @xmath2 go to infinity with the ratio @xmath17 fixed . in the second one ,
on the other hand , we take the @xmath18 limit with @xmath2 fixed .    in order to see how the whole particles from the first to @xmath2th travel
, we performed the numerical simulation of which the result is shown in fig .  2 .
we considered the time evolution of the particles from the first to 100th and assumed the hopping rate of all particles is 0.9 except the first , 25th , 50th and 75th particles whose hopping rate is 0.8 .
thus these four particles represent the defect particles .
the typical example of the time evolution of the 100 particles from time 0 through 3000 is illustrated in fig .
2(a ) . in this figure ,
the vertical axis represents the position coordinate introduced in fig .
1(b ) , whereas the horizontal axis is time axis , and each particle is expressed as @xmath19 . from the macroscopic point of view such as fig .
2(a ) , we only find that they travel in a group .
however , if we zoom the configuration of the particles in a certain time zone , we can see the microscopic pattern of the configuration . here
we show two characteristic configurations which are depicted in fig .
2(b ) and 2(c ) .
2(b ) ( resp.(c ) ) represents the configuration of the particles when time is around 200 ( resp .
3000 ) . when the time zone which we focus on is early enough as in fig .
2(b ) , the particles move forming only one group and the four defect particles are included in the group .
on the other hand , when time has passed sufficiently as in fig .
2(c ) , they form the four groups automatically and the top of each group is the defect particle .
the similar platoon structure was also discussed in the study of the tasep with disordered hopping rates  @xcite .
in fact , as will be explained in the next section , the position fluctuation of a tagged particle also changes with time corresponding to the change of the configuration .
our whole discussions in this article are based on a fact that the tagged particle problem for the step initial condition has a nice combinatorial structure which admits us to obtain a closed expression for the multi - time joint distribution in the form of the fredholm determinant .
the results are summarized in the following theorem .
the proof will be presented in section  [ ftts ] .
note that due to the rule of the tasep , the particle labeled @xmath2 can not move when @xmath20 .
thus , throughout the article , we consider the multi - time function   under the condition @xmath21    * theorem 1 *    _ _    when holds , @xmath22 here @xmath23 is the fredholm determinant defined as @xmath24 where @xmath25 and the kernel @xmath26 is given by @xmath27 here @xmath28 in   and   denotes a contour enclosing the origin anticlockwise with radius @xmath29 and @xmath30 in   satisfy the conditions @xmath31 and @xmath32 .
using the results in theorem 1 , one can study the asymptotics of the joint distribution . in this subsection
, we consider the scaling limit such that both @xmath4 and @xmath2 go to infinity with their ratio @xmath34 fixed .
@xmath35 represents the scaled time .
notice that we also take the limit of the label of the tagged particle @xmath2 .
this may sound strange if the particle we are focusing on is varied as @xmath2 goes to infinity .
rather we think that the results below give the asymptotic behaviors of the path statistics of the tagged particle for large but finite @xmath2 which is fixed . alternatively one might be able to interpret them as giving the asymptotic behaviors of the path statistics of the tagged particle near a large but finite @xmath4 .
we also discuss the effect of defect particles .
in particular , we consider the situation where there are finite @xmath36 defect particles with stay rates @xmath37 in front of the tagged particle . for this purpose , we set the stay rate @xmath10 of @xmath5th particle as follows . for some set @xmath38 with
@xmath36 fixed , we assume @xmath39 note that the particle with the rate @xmath40 represents the defect particle whose label is @xmath41 whereas the remaining particles with the rate @xmath0 are the normal particles .      the average position of the tagged particle divided by @xmath2 has a deterministic limit as @xmath33 .
let us call it @xmath42 , @xmath43 we also define @xmath44 .
( note that @xmath45 is the stay rate of the slowest defect particle . )
if @xmath46 , which is the case where the slowest defect particle is slower than the normal particles , we obtain @xmath47 where @xmath48 and @xmath49 @xmath50 can be understood from theorem 1.1 . in  @xcite .
if @xmath51 , the average position is represented as only the first and second cases in  .
( @xmath52 and @xmath53 do not appear in the case . ) note that @xmath50 does not depend on the stay rate @xmath45 of the slowest defect particle whereas @xmath53 does .
@xmath42 is illustrated in fig .
next we consider the position fluctuations of the tagged particle around the average position .
there are four typical regions shown in fig .
3 . in each region , we can obtain the result on the position fluctuations by focusing on a point within the region and taking a proper scaling around the point .
note that in fig .  3 , the regions 2 and 4 spread in the form of a line whereas the regions 1 and 3 are point - like . thus , in the region 1(resp .
3 ) , we always consider the fluctuation property around the point @xmath54 ( resp .
@xmath55 ) .
this is the region where the tagged particle just starts to move .
this corresponds to @xmath56 .
we set @xmath57 where @xmath58 .
note that the time @xmath59   is approximately equal to the time on which some leftmost holes arrive at the site occupied by the tagged particle .
one easily finds that the arrival time can be scaled as in   which is the same scaling as the central limit theorem .    applying the scaling to theorem 1
, we have the theorem as follows . the proof will be given in section  [ pot21 ] .
* theorem 2 - 1 * _ @xmath60 here the fredholm determinant in the right hand side is defined in where @xmath61 .
@xmath62 is defined in  .
the kernel is given by @xmath63 here @xmath64 can be represented by use of the parabolic cylinder function @xmath65 in  @xcite as @xmath66 where @xmath67 is taken to be positive , and @xmath68 where @xmath69 in   represents the contour enclosing the origin anticlockwise with radius 1 , and @xmath70 is the hermite polynomial with degree @xmath36  @xcite . _    the special case where @xmath71 has appeared in the distribution of the height fluctuation in the `` critical regime '' of the oriented digital boiling model  @xcite .
recently the distribution in one time case ( @xmath72 ) has been also given in  @xcite .
our formula above gives a generalization to the multi - time version .
in particular , the limiting distribution in the case @xmath72 can be described as @xmath73 where @xmath74 note that in   the summation on @xmath75 where @xmath76 can be omitted .
thus in this case , we can obtain exactly the limiting probability by calculating the determinants with finite rank .
another expression of the kernel   has recently appeared as the `` discrete hermite kernel '' in  @xcite . in  @xcite ,
the specific values of the probability are given in the case @xmath77 .
this is the region where the effect of the defect particles does not affect the time evolution of the tagged particle and the dynamics of the ordinary tasep is dominant .
the typical picture of the time evolution of the whole particles from 1st to @xmath2th in this region is illustrated in fig .
2(b ) . in this figure , we are interested in the position fluctuations of the bottommost particle .
let us scale as @xmath79 where @xmath80 and @xmath81 .
@xmath50 is defined in   and @xmath82 note that the scaling exponent 1/3 in   also appears in an anomalous diffusion of the tagged particle  @xcite .
the scaling exponents @xmath83 and @xmath84 in   and   are characteristic of the one - dimensional kpz universality class .
they are expected to be universal for all models belonging to the kpz universality class .    under the scaling
defined above , we get the following theorem for the scaling function .
this function is also universal but may depend on various elements such as initial conditions and the effect of the defect particles .
for example , as we will show in theorem 2 - 3 , the scaling function in the region 3 is different from that in region 2 although we take the same scaling   and   in both regions .
thus , the function provides more detailed viewpoint in the kpz universality than the scaling exponents .
the proof of this theorem will be given in section  [ pot223 ] .
* theorem 2 - 2 * _ @xmath85 here the right hand side is the fredholm determinant defined as @xmath86 where @xmath87 are defined in terms of @xmath62  , as @xmath88 and the kernel @xmath89 is given by @xmath90 _
the kernel @xmath91 is called the extended airy kernel  @xcite . the process characterized by the fredholm determinant with this kernel
is called the airy process  @xcite .
this process appears as the limiting process of the largest eigenvalue in dyson s brownian motion model  @xcite of the unitary class .
the model is described as @xmath92 hermitian matrix where each independent element of @xmath93 obeys the ornstein - uhlenbeck process .
the transition probability density @xmath94 from matrix @xmath95 to @xmath96 during @xmath97 is given by @xmath98 where @xmath99 is the normalization constant .
if we choose the initial matrix @xmath100 to be gue random matrix , the joint density function of the probability that matrix @xmath96 is at time @xmath101 is represented as the hermitian multi - matrix model , @xmath102 where @xmath103 is the normalization constant .
let @xmath104 be the largest eigenvalue of @xmath95 and we consider the quantity , @xmath105 when we set @xmath106 the limiting distribution   is described as @xmath107 in the case @xmath72 , the distribution is the gue tracy - widom distribution  @xcite .
hence our theorem 2 - 2 says that in the appropriate scaling limit the dynamics of the tagged particle in this region is equivalent to the dynamics of the largest eigenvalue of dyson s brownian motion of the unitary type . in the context of the one dimensional kpz universality class
, the airy process has already appeared in the study of the polynuclear growth ( png ) model . in  @xcite
, the process has first appeared as the equal - time multi - point height fluctuation .
recently in  @xcite , it has been shown that the process also describes the correlations among any `` space - like '' points in the space - time plane .
our result is closely related to the latter situation .
this region is the border between the region 2 where the normal particles in the tasep are dominant , and the region 4 where only the finite number of defect particles are dominant .
first we consider the case where the stay rates @xmath109 of the defect particles are distinct . taking the same scaling as   with @xmath108 , we obtain the following theorem .
* theorem 2 - 3 * _ @xmath110 where the fredholm determinant and @xmath111 in the right hand side are defined in   and   respectively .
the kernel is given by @xmath112 where @xmath91 is the extended airy kernel  .
_    this kernel has appeared in  @xcite in the study of the height fluctuation property of the png model with external sources .
when @xmath71 , it has been known that the fredholm determinant of this kernel is described as the distribution of the larger of the largest eigenvalues in two independent goes , which is denoted as goe@xmath3  @xcite . on the other hand ,
when @xmath113 , the fredholm determinant becomes the gue tracy - widom distribution .
thus describes the transition of the largest eigenvalue distribution between gue and goe@xmath114 .
next we consider the situation where the hopping rates @xmath115 are more or less the same .
in addition to the scaling   , we set the hopping rate of the defect particles as @xmath116 where @xmath117 is defined in  .
note that the parameters @xmath118 characterize the inhomogeneity of the hopping rates and the case @xmath119 corresponds to the situation where the hopping rates degenerate completely .    under these settings , we obtain the following theorem .    * theorem 2 - 3 * _ @xmath120 where the fredholm determinant and @xmath111 in the right hand side are defined in   and   respectively and the kernel is given by @xmath121 where the contour @xmath122 runs from @xmath123 to @xmath124 passing the down side of the points @xmath125 .
_    in the case @xmath72 , this fredholm determinant has appeared in  @xcite .
note that if we set @xmath126 for some @xmath127 and @xmath128 for @xmath129 , we can realize the situation where the maximum of @xmath130 is unique such as the former case .
indeed in the setting above , the kernel @xmath131   is reduced to @xmath132  . in that sense ,
the theorem is the generalization of theorem 2 - 3 .
later we will give only the proof of the theorem in section  [ pot223 ] .
the fredholm determinant   also describes the limiting largest eigenvalue distribution in the following multi - matrix model .
let @xmath133 be the @xmath92 hermitian matrices .
analogous to , the joint density function of the model is defined as @xmath134 here @xmath135 is the normalization constant and @xmath136 are parameters of the model such that @xmath137 gauges the strength of the connection between matrix @xmath138 and @xmath95 .
@xmath139 represents the external source with rank @xmath36 of the model . under the scalings   and @xmath140 the limiting joint distribution of the largest eigenvalues @xmath104 of @xmath141
is described as the fredholm determinant  , @xmath142 in the case @xmath143 , was shown in  @xcite .
we can also show this equation for the general @xmath36 case in a similar manner .
this is the region where the effect of the defect particles has become dominant for the dynamics of the tagged particle .    as in the region 3 , we first consider the case @xmath145 . in this case
, one can expect that the dynamics is effectively the same as that of the slowest defect particle with stay rate @xmath146 . to see nontrivial correlations , times @xmath147 should be macroscopically separated .
we set @xmath148 where @xmath81 , @xmath149 is given in  , and @xmath150 under this setting and introducing the parameter @xmath151 such that @xmath152 we get the theorem as follows .
* theorem 2 - 4 * _ @xmath153 where the fredholm determinant and @xmath111 are defined in   and   respectively and @xmath154 _
this kernel represents the propagation of a brownian particle ( a particle obeying the ornstein - uhlenbeck process ) .
actually the two point distribution is calculated as @xmath155    next we consider the case @xmath156 .
the typical picture of the time evolution of the particles is depicted in fig .
2(c ) . in this figure
, it seems that four groups of which each top particle is the defect particle are formed and each group behaves like one particle in the tasep .
thus in the region , one can expect that the system is effectively the same as @xmath36 tasep particles .
in addition to the scaling    , we also set @xmath157 note that the scaling exponent 1/2 in   is the same as that in the case of the central limit theorem .
however the following theorem indicates that the limiting process of @xmath14 depends on the number @xmath36 of the defect particles . note that this is not the information which the scaling exponent can detect .
* theorem 2 - 4 * _ @xmath158 the fredholm determinant and @xmath111 in the right hand side are defined in and   respectively and the kernel is given by @xmath159 where @xmath160 in , @xmath161 represents the contour enclosing @xmath162 anticlockwise and @xmath163 represents the arbitrary path running from @xmath164 to @xmath165 .
_    the kernel above has been given in  @xcite .
as in theorem 2 - 3 , it is reduced to @xmath166   if we set @xmath167 for some @xmath127 and @xmath168 for @xmath129 .
thus , later we will give only the proof of this theorem in section  [ pot24 ] .
the fredholm determinant also describes the joint distribution function of the largest eigenvalue in the hermitian multi - matrix model  , where we assume the rank of @xmath169 is @xmath36 and @xmath170 .
the joint distribution of the largest eigenvalues @xmath104 of @xmath95 can be described as the fredholm determinant in the above theorem , @xmath171 note that in the equation above , we do not take the limit as the rank @xmath36 goes to infinity .
here we consider another scaling limit such that @xmath2 is fixed and @xmath4 goes to infinity .
we set @xmath172 and consider the asymptotic behavior as @xmath173 goes to infinity .
we get the following theorem .
the proof of the theorem will be given in section  [ pot3 ] .    * theorem 3 * _
@xmath174 where the right hand side is the same as that of   with @xmath175 .
_      in this subsection , we consider the continuous time version of the tasep which is also usually discussed in the study of the asep . in the case without defect particles , the rule is defined as follows .
let @xmath176 be the time variable which can take any positive real number . between time @xmath176 and @xmath177 , a particle can hop to its right neighboring site with probability @xmath178 if the site is empty .
if the site is occupied , the particle stays at the same site with probability 1 .
we easily find that the continuous time version can be realized from the discrete time version defined in section  [ model ] by taking the limit , @xmath179 with @xmath180 fixed .
in addition , the scaling limit in section  [ sclim1 ] can also be studied by taking the @xmath181 limit with their ratio @xmath182 fixed . applying   and   to the theorems in section  [ sclim1 ] , we can rewrite the theorems for the continuous time tasep . here
we show only the result which corresponds to theorem 2 - 2 .
analogous to   and  , we scale time @xmath183 and the position @xmath184 as @xmath185 where @xmath186 and @xmath187 then we have @xmath188 here the right hand side is defined in    .
in this section , we discuss a relationship between the time evolution of a tagged particle in the tasep with the step initial condition described in section  [ model ] and the stochastic growth of a young diagram  @xcite which turns out to be a special case of the schur process  @xcite . in the last part of this section , the proof of theorem 1 in section  [ multitimed ]
is obtained .    in the study of the tasep ,
the technique using the enumeration of young diagrams was used in  @xcite to study the current fluctuation .
the discussion in  @xcite was based on the relation of the tasep to the directed polymer problem on a matrix of which the elements obey the geometric distribution . in this section , on the other hand , we consider the correspondence between the directed polymer problem on a 01 matrix  @xcite and the tasep , which is a similar but different mapping . for the analysis of the tagged particle problem ,
the mapping is more natural and convenient .
let us consider an ensemble of @xmath189 01 matrices @xmath190 , whose matrix elements are either 0 or 1 .
all matrix elements are independent random variables and each element obeys the bernoulli distribution , @xmath191 here @xmath192 is taken to be the same value as the hopping rate of the @xmath127th particle of the tasep defined in section  [ model ] .
each realization of the 01 matrix @xmath193 can be translated to a time evolution of the tasep as follows .    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ @xmath194(resp .
1 ) : between time @xmath195 and @xmath196 , the @xmath127th particle tries to hop to the right neighboring site ( resp
. stays at the same site ) .
when the right neighboring site is occupied at time @xmath195 , particle can not move due to the exclusion effect even if @xmath194 .
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    since @xmath194 could mean both hopping or stay depending on whether the target site is occupied or empty , the mapping from a @xmath197 matrix to the time evolution of the tasep is not one - to - one . in fig .  4 , we depict an example of the time evolution in the tasep and the corresponding two 01 matrices .
note that in the step initial condition , the @xmath127th particle can not move until time @xmath198 and that the time evolution of the @xmath127th particle at time step between @xmath5 and @xmath199 does not depend on that of the @xmath198th particle at the time step .
thus one finds that the 01 matrix defined above has all informations needed for the time evolution of the @xmath2th particle from time 0 through @xmath200 .
let us define a left - down path @xmath201 of a @xmath189 01 matrix @xmath190 by @xmath202 note that the row indices @xmath203 are _ strictly _ increasing while the column indices @xmath204 are _ weakly _ decreasing . in @xmath205 , the position of @xmath206
is always located on the down side or left - down side of the position of @xmath207 in the 01 matrix .
we also define the quantity @xmath208 by @xmath209 where @xmath210 denotes the number of elements in the path @xmath205 .
if we regard the 01 matrix as the random media in a plane , @xmath201 as a spatial configuration of a polymer chain , and @xmath211 as its energy , we can interpret   as a statistical mechanical problem of a directed polymer .
this kind of problem is studied in  @xcite .
now we find that this maximum length of the left - down path is directly related to the position of a tagged particle as @xmath212 where @xmath213 represents the number of times the @xmath2th particle stays at the same site from time 0 through @xmath214 .
the equation ( [ prop1 ] ) enables us to study the dynamics of the tasep in terms of the 01 matrices and plays a fundamental role in our subsequent analysis .
we can show   in the following way .
for the case @xmath215 ( the case where we tag the first particle ) , we can easily check this equation , since @xmath216 defined in   is simply the total number of the value 1 in the sequence @xmath217 and it is clear that this quantity represents the number of times the particle stays from @xmath218 through @xmath219 .    for the case @xmath220
, we can prove   by mathematical induction about the row @xmath219 .
here , we mainly consider only the case @xmath221 . for the case
@xmath222 , we can check   in a similar fashion to this case .    when @xmath223 , we can check   easily for all four cases ( @xmath224 and @xmath225 ) .
we assume that   holds for @xmath226 , and @xmath215 and @xmath227 , i.e. , @xmath228 under this assumption , we will show that @xmath229 as follows . from
, one easily finds @xmath230 where @xmath231 , and @xmath232 and @xmath227 . considering these properties
, we classify the problem into two cases ,    ( i ) : :    @xmath233 ( ii ) : :    @xmath234    the case ( i ) corresponds to the situation where the two particles occupy the neighboring sites at @xmath235 .
thus the second particle ( the tagged particle ) must stay at time step between @xmath236 and @xmath237 .
hence one finds @xmath238 regardless of the value of @xmath239 . on the other hand , due to  ,   and
, we have @xmath240 thus from  ,   and  , we have for the case ( i ) , @xmath241    the case ( ii ) corresponds to the situation where at @xmath235 , the distance between the first and the second particles is at least one site .
note that in this case we do not need to consider the exclusion effect of the first particle and thus the situation is essentially the same as the case @xmath215 .
we have @xmath242 on the other hand , from   and   we can also check that @xmath243 thus we find that   holds also for the case ( ii ) and therefore   holds for @xmath221 .    in the case
@xmath244 , we can also classify the situation into two cases as in ( i ) and ( ii ) for @xmath221 , i.e. , the case ( i) where all @xmath75 particles are packed and ( ii) where the particles form several groups which are separated by some successive empty sites . for the case ( i)
, we can check   in a similar manner to the case @xmath221 . for the case ( ii)
, we can easily show that the case is essentially the same as the case ( i) for @xmath245 which is smaller than @xmath75 .
the quantity we would like to discuss is the multi - time distribution of @xmath14 , @xmath246 where @xmath247 represents the probability measure of the tasep defined in section  [ model ] . due to the result in section  [ 01table](especially   and  )
, one finds @xmath248 where @xmath249 represents the probability measure of 01 table defined in   and @xmath250 is defined in  .    in this section
, we further restate our problem in terms of the combinatorics of young tableaux .
for the definitions and basic properties of young tableaux and related subjects , we refer the readers to @xcite . there is a bijective mapping between an @xmath189 01 matrix and a pair @xmath251 where @xmath252 ( the transpose of @xmath253 ) and @xmath254 are the semistandard young tableaux ( ssyt ) with the condition that the shape of @xmath253 is the same as that of @xmath254 .
the mapping is called the dual robinson - schensted - knuth ( rsk ) algorithm . in order to get @xmath251
, we first construct a two - line array ( or generalized permutation ) , @xmath255 where @xmath256 , from the 01 matrix by listing @xmath257 s for which @xmath258 .
next , we construct @xmath253 ( resp .
@xmath254 ) by arranging on a plane the figures in the second row ( resp . the first row ) . note that each positive integer constructing @xmath253(resp .
@xmath254 ) is less than @xmath2(resp .
@xmath219 ) . for details of the algorithm ,
see  @xcite .
the generalized permutation and @xmath251 corresponding to the 01 table in fig .
5(b ) are given in fig .
now we express @xmath208   as the quantity related to a pair @xmath251 .
let us denote the shape of @xmath253 ( or @xmath254 ) by the young diagram , @xmath259 where @xmath260 represents the length of @xmath5th row of the tableaux @xmath253 ( or @xmath254 ) . in the example of fig .
6(b ) , we see @xmath261 .
then we have @xmath262 here @xmath263 means the transpose of the young diagram @xmath264 .
thus @xmath265 is equal to the length of the first column of @xmath266 .
.   may be shown by examining the dual rsk algorithm directly but is also understood as follows . for the sequence of the numbers @xmath267
where @xmath268 , we define a nondecreasing ( resp .
nonincreasing ) subsequence @xmath269 such that @xmath270 ( resp .
@xmath271 ) . from the construction
it is not difficult to see that @xmath208 is equivalent to the length of the longest nonincreasing subsequence of the second row @xmath272 in a generalized permutation  .
thus , @xmath208 is also regarded as the length of the longest nondecreasing subsequence of the opposite sequence @xmath273 of the second row in  .
we can apply the ( normal ) rsk algorithm  @xcite to this reverse sequence to obtain a ssyt .
let us denote the ssyt by @xmath29 .
then ( [ dl ] ) is a consequence of the facts that the length of the longest nondecreasing sequence in the reverse sequence is equal to the length of the first row of @xmath29 ( this is known as a property of the normal rsk algorithm  @xcite ) and the symmetry property , @xmath274 whose proof can be found in appendix a of  @xcite .
notice that in the example of fig .
5(a ) , @xmath275 thus @xmath276 while in fig .
6(b ) , @xmath277 .
thus from  , we get @xmath278 here @xmath279 in the right hand side represents the probability measure of the set of young diagrams obtained by the dual rsk algorithm from the 01 matrix defined in and @xmath280 means the young diagram obtained by applying the dual rsk algorithm to the @xmath281 01 submatrix @xmath282 of the @xmath189 01 matrix  .    from and
, we managed to interpret the time evolution of a tagged particle in the tasep as the time evolution of @xmath283 , the length of the first column , regarding @xmath284 as a time parameter .      in order to investigate the right hand side of
, we consider the following probability , @xmath285 for a given set of young diagrams , @xmath286 .
this equation describes the growth of a young diagram .
.   represents the joint distribution of them .
7 illustrates the growth of young diagram which corresponds to fig .
notice that from   and  , the growth of @xmath287 describes the time evolution of @xmath2th particle in the tasep .
let us denote the @xmath288 as the part of the @xmath254 tableau obtained by picking up only the symbols @xmath289 . from the rule of the dual rsk algorithm  @xcite
, we find that @xmath288 is equivalent to the @xmath254 tableau obtained by applying the dual rsk algorithm to the the submatrix @xmath290 .
thus one has @xmath291 where @xmath292 denotes the shape of @xmath293 .
hence the @xmath254 tableau records the growth of young diagram , that is , it one - to - one corresponds to the set of the young diagrams @xmath294 .
compare the tableau @xmath254 in fig .
6(b ) with @xmath295 in fig .  7
. we can recognize that @xmath292 constructed from the tableau @xmath254 in fig .
6(b ) is equivalent to @xmath296 in fig .
hence we have @xmath297 here @xmath298 represents the probability measure on the set of the pairs @xmath299 obtained by the dual rsk algorithm from the @xmath189 01 matrices defined in .
hence it is found that the growth process of the young diagram is characterized by the @xmath254 tableau in the dual rsk correspondence .
this type of growth process has been also discussed in  @xcite .
remembering @xmath253 is constructed from the second row of a generalized permutation , we find that the symbol @xmath127 in @xmath253 represents the column index of the element of the 01 matrix where @xmath193=1 .
thus the total number of @xmath127s in @xmath253 is equal to the total number of the figure 1 in the @xmath127th column of the 01 matrix .
hence , for a given @xmath300 , one finds @xmath301 where @xmath302 means the total number of the symbol @xmath5 in @xmath300 . applying the above equation to the combinatorial definition of the schur function , @xmath303 where the summation is taken to all ssyt @xmath304 with shape @xmath305
, we have @xmath306 where @xmath307 . the factor
@xmath308 ensures that the probability is zero unless @xmath309 and @xmath310 have no two squares in the same column . thus combining   with
, we finally find @xmath311 this measure is the special case of the joint distribution function in the schur process  @xcite .
since the schur function is expressed as a determinant , the function   is described as a product of determinants .
furthermore , it is known that the process has a remarkable mathematical structure that the correlation function of @xmath312 can be expressed as a determinant and its scaling limit can be discussed exactly .
the schur process also appears in other fields of physics and mathematics such as random growth process  @xcite , melting problem of a three - dimensional crystal  @xcite , random tiling model  @xcite , and so on .      from  , and , we find @xmath313 here @xmath314 represents the probabilistic measure of the schur process  .    in general , the schur process is defined as follows . for the set of the young diagrams such that @xmath315 the following weight is assigned @xmath316 here
@xmath317 is the schur function with the specialization of algebra @xmath318 . in our case
, we notice only two cases @xmath319 where @xmath320 is defined in  . substituting to  , @xmath321 where @xmath322 and applying theorem 1 in  @xcite or theorem 2.2 in  @xcite to this case
, we finally obtain the fredholm determinant representation(- ) of .
in this section , we discuss the scaling limit of the multi - time distribution function   by applying the saddle point method to the kernel  - in theorem 1 .
we consider two types of scaling limit .
the first case is explained in section  [ sclim1 ] .
we focus on the situation where both time @xmath4 and the label of the tagged particle @xmath2 tend to infinity when there are @xmath36 defect particles with the stay rates @xmath323 in front of the tagged particle .
the results on the limiting distributions are summarized as theorem 2 .
their proofs are given in section 5.1 , 5.2 and 5.3 .
in the second case , we take @xmath4 to be infinite with @xmath2 fixed .
the limiting distribution is described in theorem 3 in section  [ sclim2 ] , whose proof is given in section 5.4 .      in this section ,
we consider the scaling limit in the region 1 in fig .
we scale time @xmath59 as  . in order to calculate the limiting distribution , we rewrite the kernel  - with the condition   as @xmath324 where @xmath325 and @xmath326 .
@xmath327 denotes the contour with radius @xmath29 enclosing the origin anticlockwise and @xmath328 satisfy the condition , @xmath329 for @xmath330 and @xmath331 for @xmath332 . using   and the relation @xmath333 one finds @xmath334 here the function @xmath335 and @xmath336 are given by @xmath337 where the radius @xmath338 of the contour @xmath339 in   is taken such that @xmath340 . in the following discussion
, we evaluate the asymptotics of the kernel   by applying the saddle point method to the functions @xmath341 .
first we consider the asymptotics of @xmath342 .
substituting   with  , we set @xmath343 where @xmath344 from the equation @xmath345 , we find the critical point @xmath346 is given by @xmath347 scaling the variable @xmath348 around the critical point as @xmath349 where @xmath350 , we evaluate the asymptotics of the term @xmath351 in   by the saddle point method , @xmath352 as @xmath353 .
next we consider the asymptotics of other terms in  . from  , one obtains @xmath354 transforming @xmath355 into @xmath356 , @xmath357 and considering  , one gets @xmath358 note that by the transformation  , the function @xmath359   in the fredholm determinant   is changed to @xmath360 . from  - ,
we obtain the asymptotic form of @xmath342 , @xmath361 where @xmath362 .
next we consider the asymptotics of @xmath363  .
similar to @xmath342 , we set @xmath364 here @xmath365 is defined in   and @xmath69 denotes the contour enclosing the origin anticlockwise with radius 1 . applying the saddle point method to this equation in the same way as @xmath342
, we get @xmath366 thus from  ,  , and  , we find @xmath367 where @xmath368 introducing the function @xmath369 by @xmath370
one finds that the function satisfies weber s equation  @xcite , @xmath371 the initial condition was given in  @xcite as @xmath372 where @xmath373 represents the gamma function .
the function satisfying   is known as the parabolic cylinder function  @xcite .
in addition , one easily finds @xmath374 can be represented as   from the integral representation of the hermite polynomial @xmath375  @xcite , @xmath376    hence , noticing the prefactor @xmath377 does not contribute to the determinant in  , we get the limiting kernel  - .      in this section
we consider the scaling limit of the regions 2 and 3 in fig .
we fix the scaled time @xmath35   such that @xmath378 for the region 2 and @xmath108 for the region 3 , where @xmath379 .
we scale the time @xmath101 , the position @xmath184 of the tagged particle at @xmath101 , and the stay rates of the defect particles @xmath380 as  ,   and   respectively .
first we analyze the asymptotics of @xmath381 in   by the saddle point method . using the variable @xmath382 which will be fixed later , we set @xmath383 here @xmath384 , @xmath325 , and @xmath326 . the function @xmath385 is defined as @xmath386    we fix the value of @xmath382 in such a way that two saddle points of @xmath385 merge to one point .
we have @xmath387 and the double critical point @xmath388 of @xmath385is obtained as @xmath389 note that @xmath390 is satisfied with @xmath382 fixed as  .
thanks to the property that @xmath391 when @xmath392 which is the condition of regions 2 and 3 , it is found that we can deform the contours @xmath393 around the double saddle point  .
thus we scale @xmath394 and @xmath395 as @xmath396 where @xmath117 is defined in  . noting   and the relation @xmath397 we have @xmath398 hence we had the asymptotic form of the term @xmath399 in  .
next we consider that of other terms in  .
combining   with the relation @xmath400 we find @xmath401 furthermore , associated with  , we scale @xmath402 as @xmath403 here @xmath50 is defined in  . using  , and  , we get @xmath404 note that in  , the asymptotic form is different between regions 2 and 3 .
this leads to the difference of the limiting distribution between theorems 2 - 2 and 2 - 3.    thus we obtain the asymptotic form of @xmath405 , @xmath406 here @xmath407 in  , we used the integral representation of the airy function @xmath408 and the relation @xmath409    next we consider the asymptotics of @xmath410  . using @xmath385   and @xmath382  ,
the function is rewritten as for @xmath332 , @xmath411 let the variable @xmath348 scale as @xmath412 where @xmath388 defined in   is the double saddle point of  .
due to  , we have for @xmath413 @xmath414 from the above two equations ,   and  , we get @xmath415 thus we obtain for @xmath332 , @xmath416 where @xmath417 is given in  .
from   we finally find @xmath418 hence from   and   and noting that the term @xmath417   does not affect the determinant , we finally obtain the desired kernels for the regions 2 and 3 .      in this section , we discuss the asymptotics of the region 4 in fig .  3 . in this region
we take time @xmath101 in such a way that they are macroscopically separated and the scaled time @xmath81 is taken as @xmath419 . the position @xmath184 of the tagged particle at @xmath101 and the stay rates @xmath380 of the defect particles are also scaled as   and   respectively .
deforming the contour @xmath339 in the kernel @xmath420  , we divide it into two parts , @xmath421 here for @xmath413 @xmath422 where @xmath325 and @xmath326 . in the equation above , the contour @xmath423 encloses
@xmath424 anticlockwise while the contour @xmath425 is chosen in such a way that it encloses @xmath426 anticlockwise .
first , we discuss the asymptotic form of @xmath427 .
this is rewritten as @xmath428 where @xmath429 is a variable which will be fixed later and @xmath430 we choose the value of @xmath429 in such a way that the saddle point of @xmath431 is @xmath432 where @xmath433 .
from the condition , @xmath434 we have @xmath435 scaling @xmath436 around the saddle point   as @xmath437 and noting @xmath438 we can obtain the asymptotic form of @xmath439 in   by the saddle point method , @xmath440 associated with  , we set @xmath441 for @xmath413 . from this equation ,   and  , we find @xmath442 where we used the parameter @xmath443 defined in  .
thus we eventually get @xmath444 here the contour @xmath161 encloses @xmath445 anticlockwise and @xmath163 is an arbitrary path running from @xmath164 to @xmath165 .    next we consider the asymptotics of @xmath410  . by use of @xmath446   and
@xmath429  , we set for @xmath332 , @xmath447 scaling @xmath348 as @xmath448 from this equation and  , we get @xmath449 thus we obtain for @xmath332 , @xmath450    at last , we discuss that @xmath451 in does not contribute the asymptotic form of the kernel
. the derivation can be done in the similar fashion to the one in theorem 3.1 . in  @xcite . here
we give only its outline .
first we consider the case of equal time , @xmath452 . scaling @xmath436 as
, we have @xmath453 where @xmath454 . noticing that @xmath455 for the region 4 ( @xmath456 ) and the asymptotic form of the airy function @xmath457 as @xmath458 , we find @xmath459 this indicates that under the scaling  , the kernel @xmath460 vanishes as @xmath2 goes to infinity .
next we consider the case for arbitrary @xmath461 .
we scale @xmath394 around the double saddle point @xmath462 of @xmath463 and @xmath395 around the saddle point @xmath432 of @xmath464 as @xmath465 under these scalings , we find @xmath466 where @xmath467 considering the case @xmath468 in   and the former result  , one easily finds @xmath469 hence we finally get @xmath470    from  ,   and   and noting the factor @xmath471 does not affect the determinant , we obtain the desired expression for the limiting kernel .      in this section
we discuss the scaling limit explained in section  [ sclim2 ] where time @xmath4 goes to infinity with @xmath2 fixed . in this limit
, we scale the time @xmath101 , the particle position @xmath184 , and the stay rates @xmath472 as  ,  , and   respectively .
first we consider the scaling form of @xmath381  .
changing the variables @xmath436 to @xmath473 , we have @xmath474 where the contour @xmath475 encloses @xmath476 anticlockwise , and @xmath339 encloses the origin @xmath477 and @xmath476 . introducing @xmath478 , which will be chosen later , we set @xmath479 where the parameter @xmath97 is defined in   and @xmath480 we choose @xmath481 in a way that @xmath482 has the saddle point at @xmath483 where @xmath0 is defined in  .
thus from the condition , @xmath484 , we have @xmath485 we scale the variable @xmath486 around the saddle point @xmath487 as @xmath488 considering  , we also scale @xmath402 as @xmath489 from these equations and  , we find @xmath490 from these equations , we eventually obtain @xmath491 here the contour @xmath161 encloses @xmath492 anticlockwise and @xmath163 is arbitrary path from @xmath164 to @xmath165 .
next we consider the scaling limit of @xmath410  . changing the variable @xmath348 to @xmath493 , we set @xmath494 where @xmath482 and @xmath481 are given in   and .
we set @xmath495 from this equation and   @xmath496 combining these equations , we finally find @xmath497 thus from   and   and noting that the prefactor @xmath498 does not affect the determinant , we get the desired expression for the kernel .
in this section , we give the result of monte - carlo simulations about the position fluctuation of a tagged particle in the tasep in order to check our analysis discussed in the preceding sections .
we performed the simulations in two situations , the case without defect particle and that with one defect particle .
8 shows the data of the scaled position of a tagged particle obtained by the monte - carlo simulation and the probability distribution functions which must fit them .
8(a ) corresponds to the time region @xmath499 where @xmath0 ( resp .
@xmath45 ) is the stay rate of the normal particle ( resp .
the defect particle ) . note that in the region , the both situations belong to the region 2 and due to theorem 2 - 2 , the gue tracy - widom distribution must fit the data .
8(b ) represents the case @xmath108 . in this case
the first situation comes under the region 2 . on the other hand
, the second one belongs to the region 3 where the position fluctuation are described by the limiting largest eigenvalue distribution of goe@xmath3 as explained in theorem 2 - 3 . in fig .
8(c ) , only the data for the second situation are shown for the case @xmath500 .
this case is classified as the region 4 and from theorem 2 - 4 , the fluctuation is supposed to be gaussian . in all figures , we see a good agreement between the data and the distribution functions which they must obey .      in this article , we have discussed the multi - time distribution of the tagged particle  . in order to analyze it ,
we have introduced the directed polymer problem of the 01 matrix in section  [ ftts ] .
here we consider the directed polymer problem of another random matrix where each element is geometric random variable . by this analysis
we can discuss the correlations of other quantities in the tasep with the step initial condition .
let @xmath501 be the @xmath189 matrix and the element @xmath502 is the geometric random variable , @xmath503 where @xmath504 is the parameter of the geometric distribution and we identified this with the stay rate of the @xmath127th particle of the tasep with the step initial condition defined in section  [ model ] .
for the matrix , we introduce the quantity @xmath505 analogous to @xmath208   as @xmath506 here @xmath507 is the set of right / down paths from @xmath508 to @xmath509 , @xmath510 in this setting , we consider the following quantity , @xmath511    it is known that @xmath505 can be interpreted as a quantity in the tasep as follows  @xcite ,    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ @xmath512 : in the tasep with the step initial condition , the time until which the @xmath2th particle moves @xmath219 sites to its right is @xmath4 . _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    thus the quantity   represents the correlations of the arrival times of the @xmath2th particle at the site @xmath513 .
( note that we set the site coordinate as in fig .
furthermore we easily find that @xmath514 here @xmath515 represents the number of particles which passed the site @xmath219 until time @xmath4 .
thus the probability   also means the correlations of currents between different times and different sites .    applying the similar technique in section  [ ftts ] to
, we can also represent it as the growth process of young diagram characterized by the schur process .
the result is as follows , @xmath516 where @xmath280 is the young diagram obtained by applying the normal rsk algorithm to the submatrix @xmath517 and @xmath518 is the length of its first row .
the probability measure in the right hand side of this equation is characterized by the following joint distribution function , @xmath519 where @xmath520 is the schur function .
this is also a special case of the schur process and we can get the fredholm determinant representation of  .
although in one point case , the asymptotics of this equation was discussed in  @xcite , the multi - point distribution has not been discussed yet .
in this article , we have studied the multi - time distribution function   of position fluctuations of a tagged particle in the tasep with the step initial condition .
the main results are summarized as theorems 1 , 2 and 3 in section  [ main ] .
first , we have obtained the fredholm determinant expression of   in theorem 1 . for this purpose
, we have mapped the time evolution of a tagged particle in the tasep to the growth process of young diagram which is related to the special case of the schur process .
next , using the fredholm determinant in theorem 1 , we have studied the two types of scaling limit .
the first one is the case where both time @xmath4 and the label of a tagged particle @xmath2 go to infinity . in the second one , we take the @xmath18 limit with @xmath2 fixed .
the results for the first and second ones are shown in theorems 2 and 3 respectively in section  [ main ] .    in the first scaling limit ,
if the hopping rates of all particles are the same , we can divide the scaled time into two characteristic regions according to the limiting behavior of  , the region 1 where a tagged particle begins to move and the region 2 which is after the region 1 . in the region 1 , the limiting process of a tagged particle converges to the spatially discrete process which is described in theorem 2 - 1 .
the process reflects on the discreteness of the model . in the region 2 , on the other hand ,
we have shown in theorem 2 - 2 that the limiting process becomes the airy process , which is characteristic of the one - dimensional kpz universality class .    if there are @xmath36 defect particles with small hoping rates in front of a tagged particle , the limiting distribution changes at the scaled time @xmath52 which is determined by the hopping rate of the slowest defect particle .
the limiting process around @xmath52 ( region 3 ) is equivalent to that of the largest eigenvalue in gue dyson s brownian motion with rank one external source as described in theorem 2 - 3 .
theorem 2 - 3 indicates that the rank becomes @xmath36 if the hopping rates of the defect particles are the same .
when the scaled time is after @xmath52 ( region 4 ) , we have found in theorem 2 - 4 that the process is equivalent to the one dimensional brownian motion .
if the hopping rates are degenerate , it is equivalent to the process of the largest eigenvalue of @xmath521 gue dyson s brownian motion model as described in theorem 2 - 4. this indicates that in the region 4 , the effect of the defect particles is dominant whereas that of an infinite number of normal particles is irrelevant .
theorem 3 shows the result for the second scaling limit .
the limiting distribution is also described as the fredholm determinant whose kernel is the same as that in theorem 2 - 4.
the authors would like to thank a. borodin for drawing out attention to the reference  @xcite .
the work of t.i . is supported by core research for evolutional science and technology of japan science and technology agency .
the work of t.s . is supported by the grant - in - aid for young scientists ( b ) , the ministry of education , culture , sports , science and technology , japan .
fig .  1 : definition of the tasep with the step initial condition . ( a ) during each time step
, the @xmath5th particle can hop to the right with probability @xmath8 .
however , @xmath522th particle can not hop to the right neighboring site since the site is occupied by @xmath199th particle .
( b ) step initial condition .
fig .  2 : typical time evolution of the particles from the first to 100th in the tasep with the step initial condition .
the hopping rate of all particles is 0.9 ( @xmath523 ) except the four defect particles ( the first , 25th , 50th and 75th particles ) whose hopping rate is 0.8 ( @xmath524 ) .
( a ) illustrates the whole time evolution from @xmath218 through 3000 .
its closeup around @xmath525 and @xmath526 is shown in ( b ) and ( c ) respectively .
fig .  3 : average position of the tagged particle .
the thick line shows  .
we divide this line into four regions , where @xmath35 is near @xmath54 , @xmath378 , @xmath35 is near @xmath52 , and @xmath456 . we denote them as region 1 , 2 , 3 , and 4 respectively .    fig .  4 : example of the time evolution of the tasep and corresponding 01 matrices . in this case ,
two tables is assigned for one example of the time evolution .
since the second particle can not hop to its right neighboring site between @xmath527 and @xmath528 , there are two possibilities of assigning both 0 and 1 in the @xmath529 element of the 01 matrix .        fig .  7 : growth of young diagram obtained from fig .  5(b ) by the dual rsk algorithm . the diagrams @xmath531 are obtained from the submatrices @xmath532 of fig .
the set of the diagrams @xmath533 corresponds to the @xmath254 tableaux in fig .
fig .  8 : probability distributions of the scaled position of the 100th particle from the right ( @xmath534 ) for @xmath525 ( a ) , @xmath535 ( b ) and @xmath526 ( c ) .
in these figures , @xmath536 represents the data for the case without defect particles and we set the stay rate @xmath523 . on the other hand , @xmath19 corresponds to the case where the first particle is a defect particle with @xmath524 while remaining particles are normal ones with @xmath523 .
the number of samples are 10000 for each case . in ( a )
, both cases belong to the region 2 and are fitted into the gue tracy - widom distribution shown as the dashed line . in ( b ) , the second case belongs to the region 3 since @xmath537 where @xmath538 and it is described by the distribution denoted as goe@xmath3 while the first case remains in region 2 . in ( c ) , the second cases comes under the region 4 where the distribution is described by gaussian ( dashed line ) .
( 400,171 ) ( 0,160)(a ) ( 65,171)@xmath218 ( 160,171)@xmath539 ( 257,171)@xmath527 ( 65,147)@xmath540 ( 160,147)@xmath543 ( 257,147)@xmath544 ( 65,123)@xmath545 ( 160,123)@xmath546 ( 257,123)@xmath547 ( 65,98)@xmath548 ( 50,85 ) ( 0,70)(b ) ( 70,55)@xmath232 ( 70,40)@xmath541 ( 70,26)@xmath549 ( 70,11)@xmath550 ( 70,-3)@xmath551 ( 70,-18)@xmath552 ( 82,70)@xmath542 ( 118,70)@xmath227 ( 133,70)@xmath528 ( 147,70)@xmath553 ( 100,-20 )