Eulerian Walkers as a model of Self-Organised Criticality

We propose a new model of self-organized criticality. A particle is dropped at random on a lattice and moves along directions specified by arrows at each site. As it moves, it changes the direction of the arrows according to fixed rules. On closed graphs these walks generate Euler circuits. On open graphs, the particle eventually leaves the system, and a new particle is then added. The operators corresponding to particle addition generate an abelian group, same as the group for the Abelian Sandpile model on the graph. We determine the critical steady state and some critical exponents exactly, using this equivalence.Comment: 4 pages, RevTex, 4 figure

The Oslo model, hyperuniformity, and the quenched Edwards-Wilkinson model

We present simulations of the 1-dimensional Oslo rice pile model in which the critical height at each site is randomly reset after each toppling. We use the fact that the stationary state of this sandpile model is hyperuniform to reach system of sizes $> 10^7$. Most previous simulations were seriously flawed by important finite size corrections. We find that all critical exponents have values consistent with simple rationals: $\nu=4/3$ for the correlation length exponent, $D =9/4$ for the fractal dimension of avalanche clusters, and $z=10/7$ for the dynamical exponent. In addition we relate the hyperuniformity exponent to the correlation length exponent $\nu$. Finally we discuss the relationship with the quenched Edwards-Wilkinson (qEW) model, where we find in particular that the local roughness exponent is $\alpha_{\rm loc} = 1$.Comment: 20 pages, 26 figure

Quasiadiabatic dynamics of ultracold bosonic atoms in a one-dimensional optical superlattice

We study the quasiadiabatic dynamics of a one-dimensional system of ultracold bosonic atoms loaded in an optical superlattice. Focusing on a slow linear variation in time of the superlattice potential, the system is driven from a conventional Mott insulator phase to a superlattice-induced Mott insulator, crossing in between a gapless critical superfluid region. Due to the presence of a gapless region, a number of defects depending on the velocity of the quench appear. Our findings suggest a power-law dependence similar to the Kibble-Zurek mechanism for intermediate values of the quench rate. For the temporal ranges of the quench dynamics that we considered, the scaling of defects depends nontrivially on the width of the superfluid region.Comment: 6 Pages, 4 Figure

Zero-temperature Hysteresis in Random-field Ising Model on a Bethe Lattice

We consider the single-spin-flip dynamics of the random-field Ising model on a Bethe lattice at zero temperature in the presence of a uniform external field. We determine the average magnetization as the external field is varied from minus infinity to plus infinity by setting up the self-consistent field equations, which we show are exact in this case. We find that for a 3-coordinated Bethe lattice, there is no jump discontinuity in magnetization for arbitrarily small gaussian disorder, but the discontinuity is present for larger coordination numbers. We have checked our results by Monte Carlo simulations employing a technique for simulating classical interacting systems on the Bethe lattice which avoids surface effects altogether.Comment: latex file with 5 eps figures. This version is substantially revised with new material. Submitted to J. Phys.

Drift and trapping in biased diffusion on disordered lattices

We reexamine the theory of transition from drift to no-drift in biased diffusion on percolation networks. We argue that for the bias field B equal to the critical value B_c, the average velocity at large times t decreases to zero as 1/log(t). For B < B_c, the time required to reach the steady-state velocity diverges as exp(const/|B_c-B|). We propose an extrapolation form that describes the behavior of average velocity as a function of time at intermediate time scales. This form is found to have a very good agreement with the results of extensive Monte Carlo simulations on a 3-dimensional site-percolation network and moderate bias.Comment: 4 pages, RevTex, 3 figures, To appear in International Journal of Modern Physics C, vol.

Enhanced conduction band density of states in intermetallic EuTSi3_3 (T=Rh, Ir)

We report on the physical properties of single crystalline EuRhSi$_3$ and polycrystalline EuIrSi$_3$, inferred from magnetisation, electrical transport, heat capacity and $^{151}$Eu M\"ossbauer spectroscopy. These previously known compounds crystallise in the tetragonal BaNiSn$_3$-type structure. The single crystal magnetisation in EuRhSi$_3$ has a strongly anisotropic behaviour at 2 K with a spin-flop field of 13 T, and we present a model of these magnetic properties which allows the exchange constants to be determined. In both compounds, specific heat shows the presence of a cascade of two close transitions near 50 K, and the $^{151}$Eu M\"ossbauer spectra demonstrate that the intermediate phase has an incommensurate amplitude modulated structure. We find anomalously large values, with respect to other members of the series, for the RKKY N\'eel temperature, for the spin-flop field (13 T), for the spin-wave gap ($\simeq$ 20-25 K) inferred from both resistivity and specific heat data, for the spin-disorder resistivity in EuRhSi$_3$ ($\simeq 35$ $\mu$Ohm.cm) and for the saturated hyperfine field (52 T). We show that all these quantities depend on the electronic density of states at the Fermi level, implying that the latter must be strongly enhanced in these two materials. EuIrSi$_3$ exhibits a giant magnetoresistance ratio, with values exceeding 600 % at 2 K in a field of 14 T.Comment: 6 pages, 8 figure

The Irreducible String and an Infinity of Additional Constants of Motion in a Deposition-Evaporation Model on a Line

We study a model of stochastic deposition-evaporation with recombination, of three species of dimers on a line. This model is a generalization of the model recently introduced by Barma {\it et. al.} (1993 {\it Phys. Rev. Lett.} {\bf 70} 1033) to $q\ge 3$ states per site. It has an infinite number of constants of motion, in addition to the infinity of conservation laws of the original model which are encoded as the conservation of the irreducible string. We determine the number of dynamically disconnected sectors and their sizes in this model exactly. Using the additional symmetry we construct a class of exact eigenvectors of the stochastic matrix. The autocorrelation function decays with different powers of $t$ in different sectors. We find that the spatial correlation function has an algebraic decay with exponent 3/2, in the sector corresponding to the initial state in which all sites are in the same state. The dynamical exponent is nontrivial in this sector, and we estimate it numerically by exact diagonalization of the stochastic matrix for small sizes. We find that in this case $z=2.39\pm0.05$.Comment: Some minor errors in the first version has been correcte