Numerical Study of the Correspondence Between the Dissipative and Fixed Energy Abelian Sandpile Models

We consider the Abelian sandpile model (ASM) on the large square lattice with a single dissipative site (sink). Particles are added by one per unit time at random sites and the resulting density of particles is calculated as a function of time. We observe different scenarios of evolution depending on the value of initial uniform density (height) $h_0=0,1,2,3$. During the first stage of the evolution, the density of particles increases linearly. Reaching a critical density $\rho_c(h_0)$, the system changes its behavior sharply and relaxes exponentially to the stationary state of the ASM with $\rho_s=25/8$. We found numerically that $\rho_c(0)=\rho_s$ and $\rho_c(h_0>0) \neq \rho_s$. Our observations suggest that the equality $\rho_c=\rho_s$ holds for more general initial conditions with non-positive heights. In parallel with the ASM, we consider the conservative fixed-energy Abelian sandpile model (FES). The extensive Monte-Carlo simulations for $h_0=0,1,2,3$ have confirmed that in the limit of large lattices $\rho_c(h_0)$ coincides with the threshold density $\rho_{th}(h_0)$ of FES. Therefore, $\rho_{th}(h_0)$ can be identified with $\rho_s$ if the FES starts its evolution with non-positive uniform height $h_0 \leq 0$.Comment: 6 pages, 8 figure

Geometric expansion of the log-partition function of the anisotropic Heisenberg model

We study the asymptotic expansion of the log-partition function of the anisotropic Heisenberg model in a bounded domain as this domain is dilated to infinity. Using the Ginibre's representation of the anisotropic Heisenberg model as a gas of interacting trajectories of a compound Poisson process we find all the non-decreasing terms of this expansion. They are given explicitly in terms of functional integrals. As the main technical tool we use the cluster expansion method.Comment: 38 page

Rotor-Router Walk on a Semi-infinite Cylinder

We study the rotor-router walk with the clockwise ordering of outgoing edges on the semi-infinite cylinder. Imposing uniform conditions on the boundary of the cylinder, we consider growth of the cluster of visited sites and its internal structure. The average width of the surface region of the cluster evolves with time to the stationary value by a scaling law whose parameters are close to the standard KPZ exponents. We introduce characteristic labels corresponding to closed clockwise contours formed by rotors and show that the sequence of labels has in average an ordered helix structure.Comment: 17 pages, 6 figure

Euler tours and unicycles in the rotor-router model

A recurrent state of the rotor-routing process on a finite sink-free graph can be represented by a unicycle that is a connected spanning subgraph containing a unique directed cycle. We distinguish between short cycles of length 2 called "dimers" and longer ones called "contours". Then the rotor-router walk performing an Euler tour on the graph generates a sequence of dimers and contours which exhibits both random and regular properties. Imposing initial conditions randomly chosen from the uniform distribution we calculate expected numbers of dimers and contours and correlation between them at two successive moments of time in the sequence. On the other hand, we prove that the excess of the number of contours over dimers is an invariant depending on planarity of the subgraph but not on initial conditions. In addition, we analyze the mean-square displacement of the rotor-router walker in the recurrent state.Comment: 17 pages, 4 figures. J. Stat. Mech. (2014

Green functions for the TASEP with sublattice parallel update

We consider the totally asymmetric simple exclusion process (TASEP) in discrete time with the sublattice parallel dynamics describing particles moving to the right on the one-dimensional infinite chain with equal hoping probabilities. Using sequentially two mappings, we show that the model is equivalent to the TASEP with the backward-ordered sequential update in the case when particles start and finish their motion not simultaneously. The Green functions are obtained exactly in a determinant form for different initial and final conditions.Comment: 11 pages, 4 figure

The totally asymmetric exclusion process with generalized update

We consider the totally asymmetric exclusion process in discrete time with generalized updating rules. We introduce a control parameter into the interaction between particles. Two particular values of the parameter correspond to known parallel and sequential updates. In the whole range of its values the interaction varies from repulsive to attractive. In the latter case the particle flow demonstrates an apparent jamming tendency not typical for the known updates. We solve the master equation for $N$ particles on the infinite lattice by the Bethe ansatz. The non-stationary solution for arbitrary initial conditions is obtained in a closed determinant form.Comment: 11 pages, 3 figure