Generalized Green Functions and Current Correlations in the TASEP

We study correlation functions of the totally asymmetric simple exclusion process (TASEP) in discrete time with backward sequential update. We prove a determinantal formula for the generalized Green function which describes transitions between positions of particles at different individual time moments. In particular, the generalized Green function defines a probability measure at staircase lines on the space-time plane. The marginals of this measure are the TASEP correlation functions in the space-time region not covered by the standard Green function approach. As an example, we calculate the current correlation function that is the joint probability distribution of times taken by selected particles to travel given distance. An asymptotic analysis shows that current fluctuations converge to the Airy₂ process.     

Bethe Ansatz Solution of Discrete Time Stochastic Processes with Fully Parallel Update

We present the Bethe ansatz solution for the discrete time zero range and asymmetric exclusion processes with fully parallel dynamics. The model depends on two parameters: p, the probability of single particle hopping, and q, the deformation parameter, which in the general case, |q| < 1, is responsible for long range interaction between particles. The particular case q = 0 corresponds to the Nagel-Schreckenberg traffic model with v _max = 1. As a result, we obtain the largest eigenvalue of the equation for the generating function of the distance travelled by particles. For the case q = 0 the result is obtained for arbitrary size of the lattice and number of particles. In the general case we study the model in the scaling limit and obtain the universal form specific for the Kardar-Parisi-Zhang universality class. We describe the phase transition occurring in the limit p→ 1 when q < 0.     

Dynein-Inspired Multilane Exclusion Process with Open Boundary Conditions

Motivated by the sidewise motions of dynein motors shown in experiments, we use a variant of the exclusion process to model the multistep dynamics of dyneins on a cylinder with open ends. Due to the varied step sizes of the particles in a quasi-two-dimensional topology, we observe the emergence of a novel phase diagram depending on the various load conditions. Under high-load conditions, our numerical findings yield results similar to the TASEP model with the presence of all three standard TASEP phases, namely the low-density (LD), high-density (HD), and maximal-current (MC) phases. However, for medium- to low-load conditions, for all chosen influx and outflux rates, we only observe the LD and HD phases, and the maximal-current phase disappears. Further, we also measure the dynamics for a single dynein particle which is logarithmically slower than a TASEP particle with a shorter waiting time. Our results also confirm experimental observations of the dwell time distribution: The dwell time distribution for dyneins is exponential in less crowded conditions, whereas a double exponential emerges under overcrowded conditions.     

Scale Invariance of the PNG Droplet and the Airy Process

We establish that the static height fluctuations of a particular growth model, the PNG droplet, converges upon proper rescaling to a limit process, which we call the Airy process A(y). The Airy process is stationary, it has continuous sample paths, its single time (fixed y) distribution is the Tracy–Widom distribution of the largest eigenvalue of a GUE random matrix, and the Airy process has a slow decay of correlations as y ^–2. Roughly the Airy process describes the last line of Dyson's Brownian motion model for random matrices. Our construction uses a multi-layer version of the PNG model, which can be analyzed through fermionic techniques. Specializing our result to a fixed value of y, one reobtains the celebrated result of Baik, Deift, and Johansson on the length of the longest increasing subsequence of a random permutation.     

Two Speed TASEP

We consider the TASEP on ? with two blocks of particles having different jump rates. We study the large time behavior of particles’ positions. It depends both on the jump rates and the region we focus on, and we determine the complete process diagram. In particular, we discover a new transition process in the region where the influence of the random and deterministic parts of the initial condition interact. Slow particles may create a shock, where the particle density is discontinuous and the distribution of a particle’s position is asymptotically singular. We determine the diffusion coefficient of the shock without using second class particles. We also analyze the case where particles are effectively blocked by a wall moving with speed equal to their intrinsic jump rate.     

Instability transitions and ensemble equivalence in diffusive flow

We investigate the critical behavior of one-dimensional (1D) stochastic flow with competing nonlocal and local hopping events, in context of the totally asymmetric simple exclusion process (TASEP) with a defect site in a 1D closed chain. The defect site can effectively generate various boundary conditions, controlling the total number of particles in the system. Both open and periodic-like setups exhibit dynamic instability transitions from a populated finite density phase to an empty road (ER) phase as the nonlocal hopping rate increases. In the stationary populated phase, strong clustering promoted by nonlocal skids drives such transitions and determines their scaling properties. By static and dynamic simulations, we locate such transition points, and discuss their nature and scaling properties. In the open TASEP variant, we numerically establish that the instability transition into the ER phase is second order in the regime where the entry point reservoir controls the current, while it is first order in the regime where the bulk controls the current. Since it is well known that such transitions are absent in the periodic TASEP variant, we compare our results in the open setup with those in the periodic-like setup, and discuss the issue of the ensemble equivalence. Finally, the same discussion is extended to the symmetric cases.      

Fluctuation Properties of the TASEP with Periodic Initial Configuration

We consider the joint distributions of particle positions for the continuous time totally asymmetric simple exclusion process (TASEP). They are expressed as Fredholm determinants with a kernel defining a signed determinantal point process. We then consider certain periodic initial conditions and determine the kernel in the scaling limit. This result has been announced first in a letter by one of us (Sasamoto in J. Phys. A 38:L549–L556, 2005) and here we provide a self-contained derivation. Connections to last passage directed percolation and random matrices are also briefly discussed.     

Current Distribution and Random Matrix Ensembles for an Integrable Asymmetric Fragmentation Process

We calculate the time-evolution of a discrete-time fragmentation process in which clusters of particles break up and reassemble and move stochastically with size-dependent rates. In the continuous-time limit the process turns into the totally asymmetric simple exclusion process (only pieces of size 1 break off a given cluster). We express the exact solution of the master equation for the process in terms of a determinant which can be derived using the Bethe ansatz. From this determinant we compute the distribution of the current across an arbitrary bond which after appropriate scaling is given by the distribution of the largest eigenvalue of the Gaussian unitary ensemble of random matrices. This result confirms universality of the scaling form of the current distribution in the KPZ universality class and suggests that there is a link between integrable particle systems and random matrix ensembles.     

Relaxation Times in the ASEP Model Using a DMRG Method

We compute the largest relaxation times for the totally asymmetric exclusion process (TASEP) with open boundary conditions with a DMRG method. This allows us to reach much larger system sizes than in previous numerical studies. We are then able to show that the phenomenological theory of the domain wall indeed predicts correctly the largest relaxation time for large systems. Besides, we can obtain results even when the domain wall approach breaks down, and show that the KPZ dynamical exponent z=3/2 is recovered in the whole maximal current phase.     

Determinant Representation for Some Transition Probabilities in the TASEP with Second Class Particles

We study the transition probabilities for the totally asymmetric simple exclusion process (TASEP) on the infinite integer lattice with a finite, but arbitrary number of first and second class particles. Using the Bethe ansatz we present an explicit expression of these quantities in terms of the Bethe wave function. In a next step it is proved rigorously that this expression can be written in a compact determinantal form for the case where the order of the first and second class particles does not change in time. An independent geometrical approach provides insight into these results and enables us to generalize the determinantal solution to the multi-class TASEP.     

An exclusion process for modelling fungal hyphal growth

A simple model for mass transport within a growing fungal filament is reviewed. Inspired by the role of microtubule-transported vesicles, we embody the dynamics of mass along a quasi-one-dimensional hypha with mutually excluding particles hopping on a growing one-dimensional lattice. The model is a generalisation of the totally asymmetric simple exclusion process (TASEP) to a dynamically extending lattice. We discuss mean-field and improved mean-field equations and present a phase diagram of the model's steady-state behaviour which generalises that of the TASEP. In particular we identify a region in which a shock in the density travels forward more slowly than the tip of the lattice and thus moves away from both the boundaries.     

Asymptotics in ASEP with Step Initial Condition

In previous work the authors considered the asymmetric simple exclusion process on the integer lattice in the case of step initial condition, particles beginning at the positive integers. There it was shown that the probability distribution for the position of an individual particle is given by an integral whose integrand involves a Fredholm determinant. Here we use this formula to obtain three asymptotic results for the positions of these particles. In one an apparently new distribution function arises and in another the distribution function F ₂ arises. The latter extends a result of Johansson on TASEP to ASEP, and hence proves KPZ universality for ASEP with step initial condition.     

Equilibrium time correlation functions in the low-density limit

We consider a system of hard spheres in thermal equilibrium. Using Lanford's result about the convergence of the solutions of the BBGKY hierarchy to the solutions of the Boltzmann hierarchy, we show that in the low-density limit (Boltzmann-Grad limit): (i) the total time correlation function is governed by the linearized Boltzmann equation (proved to be valid for short times), (ii) the self time correlation function, equivalently the distribution of a tagged particle in an equilibrium fluid, is governed by the Rayleigh-Boltzmann equation (proved to be valid for all times). In the latter case the fluid (not including the tagged particle) is to zeroth order in thermal equilibrium and to first order its distribution is governed by a combination of the Rayleigh-Boltzmann equation and the linearized Boltzmann equation (proved to be valid for short times).Supported in part by NSF Grant PHY 78-22302.     

Zero-Range Process with Open Boundaries

We calculate the exact stationary distribution of the one-dimensional zero-range process with open boundaries for arbitrary bulk and boundary hopping rates. When such a distribution exists, the steady state has no correlations between sites and is uniquely characterized by a space-dependent fugacity which is a function of the boundary rates and the hopping asymmetry. For strong boundary drive the system has no stationary distribution. In systems which on a ring geometry allow for a condensation transition, a condensate develops at one or both boundary sites. On all other sites the particle distribution approaches a product measure with the finite critical density ρ_c. In systems which do not support condensation on a ring, strong boundary drive leads to a condensate at the boundary. However, in this case the local particle density in the interior exhibits a complex algebraic growth in time. We calculate the bulk and boundary growth exponents as a function of the system parameters.     

Large Time Asymptotics of Growth Models on Space-like Paths II: PNG and Parallel TASEP

We consider the polynuclear growth (PNG) model in 1+1 dimension with flat initial condition and no extra constraints. The joint distributions of surface height at finitely many points at a fixed time moment are given as marginals of a signed determinantal point process. The long time scaling limit of the surface height is shown to coincide with the Airy₁ process. This result holds more generally for the observation points located along any space-like path in the space-time plane. We also obtain the corresponding results for the discrete time TASEP (totally asymmetric simple exclusion process) with parallel update.     

The Largest Eigenvalue of Real Symmetric,Hermitian and Hermitian Self-dual Random Matrix Models with Rank One External Source,Part I

We consider the limiting location and limiting distribution of the largest eigenvalue in real symmetric (β=1), Hermitian (β=2), and Hermitian self-dual (β=4) random matrix models with rank 1 external source. They are analyzed in a uniform way by a contour integral representation of the joint probability density function of eigenvalues. Assuming the “one-band” condition and certain regularities of the potential function, we obtain the limiting location of the largest eigenvalue when the nonzero eigenvalue of the external source matrix is not the critical value, and further obtain the limiting distribution of the largest eigenvalue when the nonzero eigenvalue of the external source matrix is greater than the critical value. When the nonzero eigenvalue of the external source matrix is less than or equal to the critical value, the limiting distribution of the largest eigenvalue will be analyzed in a subsequent paper. In this paper we also give a definition of the external source model for all β>0.     

Limit Processes for TASEP with Shocks and Rarefaction Fans

We consider the totally asymmetric simple exclusion process (TASEP) with two-sided Bernoulli initial condition, i.e., with left density ρ ₋ and right density ρ ₊. We study the associated height function, whose discrete gradient is given by the particle occurrences. Macroscopically one has a deterministic limit shape with a shock or a rarefaction fan depending on the values of ρ _±. We characterize the large time scaling limit of the multipoint fluctuations as a function of the densities ρ _± and of the different macroscopic regions. Moreover, using a slow decorrelation phenomena, the results are extended from fixed time to the whole space-time, except along the some directions (the characteristic solutions of the related Burgers equation) where the problem is still open.     

Dynamics of a stochastically driven Brownian particle in one dimension

We present a study on the dynamics of a system consisting of a pair of hardcore particles diffusing with different rates. We solved the drift-diffusion equation for this model in the case when one particle, labeled F, drifts and diffuses slowly toward the second particle, labeled M. The displacements of particle M exhibits a crossover from diffusion to drift at a characteristic time which depends on the rate constants. We show that the positional fluctuation of M exhibits an intermediate crossover regime of subdiffusion separating initial and asymptotic diffusive behavior; this is in agreement with the complete set of Master Equations that describe the stochastic evolution of the model. The intermediate crossover regime can be considerably large depending on the hopping probabilities of the two particles. This is in contrast to the known crossover from diffusive to subdiffusive behavior of a tagged particle that is in the interior of a large single-file system on an unbound real line. We discuss our model with respect to the biological phenomena of membrane protrusions, where polymerizing actin filaments (F) push the cell membrane (M).     

Derivation of Some Translation-Invariant Lindblad Equations for a Quantum Brownian Particle

We study the dynamics of a Brownian quantum particle hopping on an infinite lattice with a spin degree of freedom. This particle is coupled to free boson gases via a translation-invariant Hamiltonian which is linear in the creation and annihilation operators of the bosons. We derive the time evolution of the reduced density matrix of the particle in the van Hove limit in which we also rescale the hopping rate. This corresponds to a situation in which both the system-bath interactions and the hopping between neighboring sites are small and they are effective on the same time scale. The reduced evolution is given by a translation-invariant Lindblad master equation which is derived explicitly.     

Driven Tracer Particle in One Dimensional Symmetric Simple Exclusion

Consider an infinite system of particles evolving in a one dimensional lattice according to symmetric random walks with hard core interaction. We investigate the behavior of a tagged particle under the action of an external constant driving force. We prove that the diffusively rescaled position of the test particle εX(ε^-2 t), t > 0, converges in probability, as ε→ 0, to a deterministic function v(t). The function v(⋅) depends on the initial distribution of the random environment through a non-linear parabolic equation. This law of large numbers for the position of the tracer particle is deduced from the hydrodynamical limit of an inhomogeneous one dimensional symmetric zero range process with an asymmetry at the origin. An Einstein relation is satisfied asymptotically when the external force is small. Received: 5 December 1996 / Accepted: 30 June 1997