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

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

Jamming probabilities for a vacancy in the dimer model

Following the recent proposal made by Bouttier et al [Phys. Rev. E 76, 041140 (2007)], we study analytically the mobility properties of a single vacancy in the close-packed dimer model on the square lattice. Using the spanning web representation, we find determinantal expressions for various observable quantities. In the limiting case of large lattices, they can be reduced to the calculation of Toeplitz determinants and minors thereof. The probability for the vacancy to be strictly jammed and other diffusion characteristics are computed exactly.Comment: 19 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

The problem of predecessors on spanning trees

We consider the equiprobable distribution of spanning trees on the square lattice. All bonds of each tree can be oriented uniquely with respect to an arbitrary chosen site called the root. The problem of predecessors is finding the probability that a path along the oriented bonds passes sequentially fixed sites $i$ and $j$. The conformal field theory for the Potts model predicts the fractal dimension of the path to be 5/4. Using this result, we show that the probability in the predecessors problem for two sites separated by large distance $r$ decreases as $P(r) \sim r^{-3/4}$. If sites $i$ and $j$ are nearest neighbors on the square lattice, the probability $P(1)=5/16$ can be found from the analytical theory developed for the sandpile model. The known equivalence between the loop erased random walk (LERW) and the directed path on the spanning tree says that $P(1)$ is the probability for the LERW started at $i$ to reach the neighboring site $j$. By analogy with the self-avoiding walk, $P(1)$ can be called the return probability. Extensive Monte-Carlo simulations confirm the theoretical predictions.Comment: 7 pages, 2 figure

Role of gluons in soft and semi-hard multiple hadron production in pp collisions at LHC

Hadron inclusive spectra in pp collisions are analyzed within the modified quark-gluon string model including both the longitudinal and transverse motion of quarks in the proton in the wide region of initial energies. The self-consistent analysis shows that the experimental data on the inclusive spectra of light hadrons like pions and kaons at ISR energies can be satisfactorily described at transverse momenta not larger than 1-2 GeV/c. We discuss some difficulties to apply this model at energies above the ISR and suggest to include the distribution of gluons in the proton unintegrated over the internal transverse momentum. It leads to an increase in the inclusive spectra of hadrons and allows us to extend the satisfactory description of the data in the central rapidity region at energies higher than ISR.Comment: 19 pages, 20 figure