Discrete analogue computing with rotor-routers

Rotor-routing is a procedure for routing tokens through a network that can implement certain kinds of computation. These computations are inherently asynchronous (the order in which tokens are routed makes no difference) and distributed (information is spread throughout the system). It is also possible to efficiently check that a computation has been carried out correctly in less time than the computation itself required, provided one has a certificate that can itself be computed by the rotor-router network. Rotor-router networks can be viewed as both discrete analogues of continuous linear systems and deterministic analogues of stochastic processes.Comment: To appear in Chaos Special Focus Issue on Intrinsic and Designed Computatio

Algebraic arctic curves in the domain-wall six-vertex model

The arctic curve, i.e. the spatial curve separating ordered (or `frozen') and disordered (or `temperate) regions, of the six-vertex model with domain wall boundary conditions is discussed for the root-of-unity vertex weights. In these cases the curve is described by algebraic equations which can be worked out explicitly from the parametric solution for this curve. Some interesting examples are discussed in detail. The upper bound on the maximal degree of the equation in a generic root-of-unity case is obtained.Comment: 15 pages, no figures; v2: metadata correcte

Functional relations for the six vertex model with domain wall boundary conditions

In this work we demonstrate that the Yang-Baxter algebra can also be employed in order to derive a functional relation for the partition function of the six vertex model with domain wall boundary conditions. The homogeneous limit is studied for small lattices and the properties determining the partition function are also discussed.Comment: 19 pages, v2: typos corrected, new section and appendix added. v3: minor corrections, to appear in J. Stat. Mech

The arctic curve of the domain-wall six-vertex model in its anti-ferroelectric regime

An explicit expression for the spatial curve separating the region of ferroelectric order (`frozen' zone) from the disordered one (`temperate' zone) in the six-vertex model with domain wall boundary conditions in its anti-ferroelectric regime is obtained.Comment: 12 pages, 1 figur

Refined Razumov-Stroganov conjectures for open boundaries

Recently it has been conjectured that the ground-state of a Markovian Hamiltonian, with one boundary operator, acting in a link pattern space is related to vertically and horizontally symmetric alternating-sign matrices (equivalently fully-packed loop configurations (FPL) on a grid with special boundaries).We extend this conjecture by introducing an arbitrary boundary parameter. We show that the parameter dependent ground state is related to refined vertically symmetric alternating-sign matrices i.e. with prescribed configurations (respectively, prescribed FPL configurations) in the next to central row. We also conjecture a relation between the ground-state of a Markovian Hamiltonian with two boundary operators and arbitrary coefficients and some doubly refined (dependence on two parameters) FPL configurations. Our conjectures might be useful in the study of ground-states of the O(1) and XXZ models, as well as the stationary states of Raise and Peel models.Comment: 11 pages LaTeX, 8 postscript figure

Interlaced particle systems and tilings of the Aztec diamond

Motivated by the problem of domino tilings of the Aztec diamond, a weighted particle system is defined on $N$ lines, with line $j$ containing $j$ particles. The particles are restricted to lattice points from 0 to $N$, and particles on successive lines are subject to an interlacing constraint. It is shown that marginal distributions for this particle system can be computed exactly. This in turn is used to give unified derivations of a number of fundamental properties of the tiling problem, for example the evaluation of the number of distinct configurations and the relation to the GUE minor process. An interlaced particle system associated with the domino tiling of a certain half Aztec diamond is similarly defined and analyzed.Comment: 17 pages, 4 figure

Exact sampling from non-attractive distributions using summary states

Propp and Wilson's method of coupling from the past allows one to efficiently generate exact samples from attractive statistical distributions (e.g., the ferromagnetic Ising model). This method may be generalized to non-attractive distributions by the use of summary states, as first described by Huber. Using this method, we present exact samples from a frustrated antiferromagnetic triangular Ising model and the antiferromagnetic q=3 Potts model. We discuss the advantages and limitations of the method of summary states for practical sampling, paying particular attention to the slowing down of the algorithm at low temperature. In particular, we show that such a slowing down can occur in the absence of a physical phase transition.Comment: 5 pages, 6 EPS figures, REVTeX; additional information at http://wol.ra.phy.cam.ac.uk/mackay/exac

Noisy Monte Carlo: Convergence of Markov chains with approximate transition kernels

Monte Carlo algorithms often aim to draw from a distribution $\pi$ by simulating a Markov chain with transition kernel $P$ such that $\pi$ is invariant under $P$. However, there are many situations for which it is impractical or impossible to draw from the transition kernel $P$. For instance, this is the case with massive datasets, where is it prohibitively expensive to calculate the likelihood and is also the case for intractable likelihood models arising from, for example, Gibbs random fields, such as those found in spatial statistics and network analysis. A natural approach in these cases is to replace $P$ by an approximation $\hat{P}$. Using theory from the stability of Markov chains we explore a variety of situations where it is possible to quantify how 'close' the chain given by the transition kernel $\hat{P}$ is to the chain given by $P$. We apply these results to several examples from spatial statistics and network analysis.Comment: This version: results extended to non-uniformly ergodic Markov chain

Criminal narrative experience: relating emotions to offence narrative roles during crime commission

A neglected area of research within criminality has been that of the experience of the offence for the offender. The present study investigates the emotions and narrative roles that are experienced by an offender while committing a broad range of crimes and proposes a model of Criminal Narrative Experience (CNE). Hypotheses were derived from the Circumplex of Emotions (Russell, 1997), Frye (1957), Narrative Theory (McAdams, 1988) and its link with Investigative Psychology (Canter, 1994). The analysis was based on 120 cases. Convicted for a variety of crimes, incarcerated criminals were interviewed and the data were subjected to Smallest Space Analysis (SSA). Four themes of Criminal Narrative Experience (CNE) were identified: Elated Hero, Calm Professional, Distressed Revenger and Depressed Victim in line with the recent theoretical framework posited for Narrative Offence Roles (Youngs & Canter, 2012). The theoretical implications for understanding crime on the basis of the Criminal Narrative Experience (CNE) as well as practical implications are discussed

The arctic curve of the domain-wall six-vertex model

The problem of the form of the `arctic' curve of the six-vertex model with domain wall boundary conditions in its disordered regime is addressed. It is well-known that in the scaling limit the model exhibits phase-separation, with regions of order and disorder sharply separated by a smooth curve, called the arctic curve. To find this curve, we study a multiple integral representation for the emptiness formation probability, a correlation function devised to detect spatial transition from order to disorder. We conjecture that the arctic curve, for arbitrary choice of the vertex weights, can be characterized by the condition of condensation of almost all roots of the corresponding saddle-point equations at the same, known, value. In explicit calculations we restrict to the disordered regime for which we have been able to compute the scaling limit of certain generating function entering the saddle-point equations. The arctic curve is obtained in parametric form and appears to be a non-algebraic curve in general; it turns into an algebraic one in the so-called root-of-unity cases. The arctic curve is also discussed in application to the limit shape of $q$-enumerated (with $0<q\leq 4$) large alternating sign matrices. In particular, as $q\to 0$ the limit shape tends to a nontrivial limiting curve, given by a relatively simple equation.Comment: 39 pages, 2 figures; minor correction