Asymmetric Simple Exclusion Process Research Articles

Totally asymmetric simple exclusion process on a dynamic lattice with local inhomogeneity

In this manuscript, we introduce a totally asymmetric simple exclusion process on a dynamic lattice with inhomogeneity to model motor-based long-range transport along microtubules (MTs) which exhibit cycles of growth, shortening and regrowth at their growing tips. By mean-field approximation and numerical simulations, we explore phase diagrams for the particle density near the dynamic end of the lattice as well as over the entire lattice. In particular, we find seven different phases for the density over the entire lattice and explore the corresponding phase diagram by analytical analysis. Interestingly, we find that the density may have a linear profile in part of the lattice. When the lattice has a defect, we separate the lattice into two subsystems: a right subsystem with fixed length and a defect in the middle and a left subsystem with dynamic length. When the right subsystem is in maximum-current phase, we calculate an analytical approximation of the density away from the boundaries and the defect site and the density is in association with the reduced hopping rate at the defect site; the density at the left side of the defect site serves as the effective entry rate of the left subsystem. We find that with defects, the linear profiles no longer appear at the left end of the left subsystem and only five possible phases could occur. This is a step-forward to a more comprehensive mathematical description of intracellular long-range transport of organelles along MTs with instability.

Approximating the Stationary Distribution of the ASEP with Open Boundaries

We investigate the stationary distribution of asymmetric and weakly asymmetric simple exclusion processes with open boundaries. We project the stationary distribution onto a subinterval, whose size is allowed to grow with the length of the underlying segment. Depending on the boundary parameters of the exclusion process, we provide conditions such that the stationary distribution projected onto a subinterval is close in total variation distance to a product measure.

Exact steady states in the asymmetric simple exclusion process beyond one dimension

The asymmetric simple exclusion process (ASEP) is a paradigmatic nonequilibrium many-body system that describes the asymmetric random walk of particles with exclusion interactions in a lattice. Although the ASEP is recognized as an exactly solvable model, most of the exact results obtained so far are limited to one-dimensional systems. Here, we construct the exact steady states of the ASEP with closed and periodic boundary conditions in arbitrary dimensions. This is achieved through the concept of transition decomposition, which enables the treatment of the multidimensional ASEP as a composite of the one-dimensional ASEPs. Published by the American Physical Society 2024

KPZ fluctuations in finite volume

These lecture notes, adapted from the habilitation thesis of the author, survey in a first part various exact results obtained in the past few decades about KPZ fluctuations in one dimension, with a special focus on finite volume effects describing the relaxation to its stationary state of a finite system starting from a given initial condition. The second part is more specifically devoted to an approach allowing to express in a simple way the statistics of the current in the totally asymmetric simple exclusion process in terms of a contour integral on a compact Riemann surface, whose infinite genus limit leads to KPZ fluctuations in finite volume.

Askey–Wilson Signed Measures and Open ASEP in the Shock Region

Abstract We introduce a family of multi-dimensional Askey–Wilson signed measures. We offer an explicit description of the stationary measure of the open asymmetric simple exclusion process (ASEP) in the full phase diagram, in terms of integrations with respect to these Askey–Wilson signed measures. Using our description, we provide a rigorous derivation of the density profile and limit fluctuations of open ASEP in the entire shock region, including the high and low density phases as well as the coexistence line. This in particular confirms the existing physics postulations of the density profile.

Combining Reinforcement Learning and Tensor Networks, with an Application to Dynamical Large Deviations.

We present a framework to integrate tensor network (TN) methods with reinforcement learning (RL) for solving dynamical optimization tasks. We consider the RL actor-critic method, a model-free approach for solving RL problems, and introduce TNs as the approximators for its policy and value functions. Our "actor-critic with tensor networks" (ACTeN) method is especially well suited to problems with large and factorizable state and action spaces. As an illustration of the applicability of ACTeN we solve the exponentially hard task of sampling rare trajectories in two paradigmatic stochastic models, the East model of glasses and the asymmetric simple exclusion process, the latter being particularly challenging to other methods due to the absence of detailed balance. With substantial potential for further integration with the vast array of existing RL methods, the approach introduced here is promising both for applications in physics and to multi-agent RL problems more generally.

Phase diagram of generalized totally asymmetric simple exclusion process on an open chain: Liggett-like boundary conditions.

The totally asymmetric simple exclusion process with generalized update is a version of the discrete time totally asymmetric exclusion process with an additional interparticle interaction that controls the degree of particle clustering. Though the model was shown to be integrable on the ring and on the infinite lattice, on the open chain it was studied mainly numerically, while no analytic results existed even for its phase diagram. In this paper, we introduce boundary conditions associated with the infinite translation invariant stationary states of the model, which allow us to obtain the exact phase diagram analytically. We discuss the phase diagram in detail and confirm the analytic predictions by extensive numerical simulations.

Open 2–TASEP with integrable boundaries

In this paper, we explore a two-species extension of the totally asymmetric simple exclusion process (TASEP) known as ‘2–TASEP’ with open boundaries. In this model, carriers on a one-dimensional lattice exhibit distinct behaviors: loaded carriers move right, empty carriers move left, and they exchange positions at a unit rate. At the boundaries carriers are loaded at the left and unloaded at the right. We focus on the stationary state of this model. Using the integrability of this model, we introduce spectral parameters to deform it. This process uncovers an underlying algebraic structure linked to the Weyl group CL . Exploiting this structure, we analytically compute the stationary state of the inhomogeneous model, which permits, after spectral parameter specialization, to derive the exact average particle current of the original 2–TASEP.

Exit times of totally asymmetric simple exclusion processes.

We address the question of the time needed by N particles, initially located on the first sites of a finite one-dimensional lattice of size L, to exit that lattice when they move according to a TASEP transport model. Using analytical calculations and numerical simulations, we show that when N≪L, the mean exit time of the particles is asymptotically given by T_{N}(L)∼L+β_{N}sqrt[L] for large lattices. Building upon exact results obtained for two particles, we devise an approximate continuous space and time description of the random motion of the particles that provides an analytical recursive relation for the coefficients β_{N}. The results are shown to be in very good agreement with numerical results. This approach sheds some light on the exit dynamics of N particles in the regime where N is finite while the lattice size L→∞. This complements previous asymptotic results obtained by Johansson [Commun. Math. Phys. 209, 437 (2000)0010-361610.1007/s002200050027] in the limit where both N and L tend to infinity while keeping the particle density N/L finite.

Current fluctuations in a partially asymmetric simple exclusion process with a defect particle.

We study an exclusion process on a ring comprising a free defect particle in a bath of normal particles. The model is one of the few integrable cases in which the bath particles are partially asymmetric. The presence of the free defect creates localized or shock phases according to parameter values. We use a functional approach to Bethe equationsresulting from a nested Bethe ansatz to calculate exactly the mean currents and diffusion constants. The results agree very well with Monte Carlo simulations and reveal the main modes of fluctuation in the different phases of the steady state.

The steady state of the boundary-driven multiparticle asymmetric diffusion model

We consider the multiparticle asymmetric diffusion model (MADM) introduced by Sasamoto and Wadati with integrability preserving reservoirs at the boundaries. In contrast to the open asymmetric simple exclusion process the number of particles allowed per site is unbounded in the MADM. Taking inspiration from the stationary measure in the symmetric case, i.e. the rational limit, we first obtain the length 1 solution and then show that the steady state can be expressed as an iterated product of Jackson q-integrals. In the proof of the stationarity condition, we observe a cancellation mechanism that closely resembles the one of the matrix product ansatz. To our knowledge, the occupation probabilities in the steady state of the boundary-driven MADM were not available before.

Totally asymmetric simple exclusion process with local resetting in a resource-constrained environment.

Inspired by the process of mRNA translation, in which the stochastic degradation of mRNA-ribosome machinery is modeled by the resetting dynamics, we study an open totally asymmetric simple exclusion process with local resetting at the entry site in a resource-constrained environment. The effect of constrained resources on the stationary properties of the system has been comprehended in the form of the filling factor. The mean-field approximations are utilized to obtain stationary state features, such as density profiles and phase diagrams. The phase diagram possesses pure phases as well as coexisting phases, including a low-density-high-density phase separation, which did not manifest under periodic boundary conditions despite the system being closed there as well. The role of the resetting rate has been investigated on the stationary properties of the system, depending on how the filling factor scales with the system size. In contrast to the resetting model for infinite resources, two distinct phase transitions are observed for the smaller values of the filling factor leading to a change in the topology of the phase diagram. The impact of the resetting rate along with the finite-size effect has also been examined on the shock dynamics. All the mean-field results are found in remarkable agreement with the Monte Carlo simulations.

Homogenization of a reaction-diffusion problem with large nonlinear drift and Robin boundary data

We study the periodic homogenization of a reaction-diffusion problem with large nonlinear drift and Robin boundary condition posed in an unbounded perforated domain. The nonlinear problem is associated with the hydrodynamic limit of a totally asymmetric simple exclusion process (TASEP) governing a population of interacting particles crossing a domain with obstacle. We are interested in deriving rigorously the upscaled model equations and the corresponding effective coefficients for the case when the microscopic dynamics are linked to a particular choice of characteristic length and time scales that lead to an exploding nonlinear drift. The main mathematical difficulty lies in proving the two-scale compactness and strong convergence results needed for the passage to the homogenization limit. To cope with the situation, we use the concept of two-scale compactness with drift, which is similar to the more classical two-scale compactness result but it is defined now in moving coordinates. We provide as well a strong convergence result for the corrector function, starting this way the search for the order of the convergence rate of the homogenization process for our target nonlinear drift problem.

Boundary current fluctuations for the half‐space ASEP and six‐vertex model

AbstractWe study fluctuations of the current at the boundary for the half‐space asymmetric simple exclusion process (ASEP) and the height function of the half‐space six‐vertex model at the boundary at large times. We establish a phase transition depending on the effective density of particles at the boundary, with Gaussian symplectic ensemble (GSE) and Gaussian orthogonal ensemble (GOE) limits as well as the Baik–Rains crossover distribution near the critical point. This was previously known for half‐space last‐passage percolation, and recently established for the half‐space log‐gamma polymer and Kardar–Parisi–Zhang equation in the groundbreaking work of Imamura, Mucciconi, and Sasamoto. The proof uses the underlying algebraic structure of these models in a crucial way to obtain exact formulas. In particular, we show a relationship between the half‐space six‐vertex model and a half‐space Hall–Littlewood measure with two boundary parameters, which is then matched to a free boundary Schur process via a new identity of symmetric functions. Fredholm Pfaffian formulas are established for the half‐space ASEP and six‐vertex model, indicating a hidden free fermionic structure.

Mixing time for the asymmetric simple exclusion process in a random environment

Mixing time for the asymmetric simple exclusion process in a random environment

The second class particle process at shocks

We consider the totally asymmetric simple exclusion process (TASEP) starting with a shock discontinuity at the origin, with asymptotic densities λ to the left of the origin and ρ to the right of it and λ<ρ. We find an exact identity for the distribution of a second class particle starting at the origin. Then we determine the limiting joint distributions of the second class particle. Bypassing the last passage percolation model, we work directly in TASEP, allowing us to extend previous one-point distribution results via a more direct and shorter ansatz.

General criterion applicable to exclusion processes with an Ising-like invariant measure

This paper explores extensions of the totally asymmetric simple exclusion process in one dimension, assuming periodic boundary conditions. Specifically, we investigate situations where a particle’s movement rate is influenced by the positions of its two nearest neighbors. We establish the necessary and sufficient conditions for these exclusion processes to maintain an Ising-like measure as their invariant measure. Our analysis covers several renowned models, such as the cooperative exclusion process (Antal and Schütz, 2000), its dual model (Gabel and Redner, 2011), the Katz–Lebowitz–Spohn model (Katz et al., 1984), and the facilitated exclusion process (Basu and Mohanty, 2009; Gabel et al., 2010). To compare the obtained models, we will examine the average particle velocity and flux. Upon doing so, we highlight that this research contributes to the comprehension of the collective pushing phenomenon observed, for example, in biological traffic models (Epshtein and Nudler, 2003; Epshtein et al., 2003).

Codebase release 1.1 for TASEPy

The totally asymmetric simple exclusion process (TASEP) is a paradigmatic lattice model for one-dimensional particle transport subject to excluded-volume interactions. Solving the inhomogeneous TASEP in which particles’ hopping rate vary across the lattice is a long-standing problem. In recent years, a power series approximation (PSA) has been developed to tackle this problem, however no computer algorithm currently exists that implements this approximation. This paper addresses this issue by providing a Python-based package TASEPy that finds the steady state solution of the inhomogeneous TASEP for any set of hopping rates using the PSA truncated at a user-defined order.

TASEPy: A Python-based package to iteratively solve the inhomogeneous exclusion process

The totally asymmetric simple exclusion process (TASEP) is a paradigmatic lattice model for one-dimensional particle transport subject to excluded-volume interactions. Solving the inhomogeneous TASEP in which particles’ hopping rate vary across the lattice is a long-standing problem. In recent years, a power series approximation (PSA) has been developed to tackle this problem, however no computer algorithm currently exists that implements this approximation. This paper addresses this issue by providing a Python-based package TASEPy that finds the steady state solution of the inhomogeneous TASEP for any set of hopping rates using the PSA truncated at a user-defined order.

Infinite order phase transition in the slow bond TASEP

AbstractIn the slow bond problem the rate of a single edge in the Totally Asymmetric Simple Exclusion Process (TASEP) is reduced from 1 to for some small . Janowsky and Lebowitz posed the well‐known question of whether such very small perturbations could affect the macroscopic current. Different groups of physicists, using a range of heuristics and numerical simulations reached opposing conclusions on whether the critical value of is 0. This was ultimately resolved rigorously in Basu‐Sidoravicius‐Sly which established that .Here we study the effect of the current as tends to 0 and in doing so explain why it was so challenging to predict on the basis of numerical simulations. In particular we show that the current has an infinite order phase transition at 0, with the effect of the perturbation tending to 0 faster than any polynomial. Our proof focuses on the Last Passage Percolation formulation of TASEP where a slow bond corresponds to reinforcing the diagonal. We give a multiscale analysis to show that when is small the effect of reinforcement remains small compared to the difference between optimal and near optimal geodesics. Since geodesics can be perturbed on many different scales, we inductively bound the tails of the effect of reinforcement by controlling the number of near optimal geodesics and giving new tail estimates for the local time of (near) geodesics along the diagonal.