Finite-size scaling study of the d=4 site-diluted Ising

We study the four dimensional site-diluted Ising model using finite-size scaling techniques. We explore the whole parameter space (density-coupling) in order to determine the Universality Class of the transition line. Our data are compatible with Mean Field behavior plus logarithmic corrections.Comment: Contribution to LATTICE 9

A multisite microcanonical updating method

We have made a study of several update algorithms using the XY model. We find that sequential local overrelaxation is not ergodic at the scale of typical Monte Carlo simulation time. We have introduced a new multisite microcanonical update method, which yields results compatible with those of random overrelaxation and the microcanonical demon algorithm, which are very much slower, all being incompatible with the sequential overrelaxation results.Comment: 13 pages, 4 figure

The four dimensional site-diluted Ising model: a finite-size scaling study

Using finite-size scaling techniques, we study the critical properties of the site-diluted Ising model in four dimensions. We carry out a high statistics Monte Carlo simulation for several values of the dilution. The results support the perturbative scenario: there is only the Ising fixed point with large logarithmic scaling corrections. We obtain, using the Perturbative Renormalization Group, functional forms for the scaling of several observables that are in agreement with the numerical data.Comment: 30 pages, 8 postscript figure

Phase diagram of d=4 Ising Model with two couplings

We study the phase diagram of the four dimensional Ising model with first and second neighbour couplings, specially in the antiferromagnetic region, by using Mean Field and Monte Carlo methods. From the later, all the transition lines seem to be first order except that between ferromagnetic and disordered phases in a region including the first-neighbour Ising transition point.Comment: Latex file and 4 figures (epsfig required). It replaces the preprint entitled "Non-classical exponents in the d=4 Ising Model with two couplings". New analysis with more statistical data is performed. Final version to appear in Phys. Lett.

An Extended Variational Principle for the SK Spin-Glass Model

The recent proof by F. Guerra that the Parisi ansatz provides a lower bound on the free energy of the SK spin-glass model could have been taken as offering some support to the validity of the purported solution. In this work we present a broader variational principle, in which the lower bound, as well as the actual value, are obtained through an optimization procedure for which ultrametic/hierarchal structures form only a subset of the variational class. The validity of Parisi's ansatz for the SK model is still in question. The new variational principle may be of help in critical review of the issue.Comment: 4 pages, Revtex

Experimental proof of faster-is-slower in systems of frictional particles flowing through constrictions

The "faster-is-slower" (FIS) effect was first predicted by computer simulations of the egress of pedestrians through a narrow exit D. Helbing, I. J. Farkas, and T. Vicsek, Nature (London) 407, 487 (2000)NATUAS0028-083610.1038/35035023]. FIS refers to the finding that, under certain conditions, an excess of the individuals'' vigor in the attempt to exit causes a decrease in the flow rate. In general, this effect is identified by the appearance of a minimum when plotting the total evacuation time of a crowd as a function of the pedestrian desired velocity. Here, we experimentally show that the FIS effect indeed occurs in three different systems of discrete particles flowing through a constriction: (a) humans evacuating a room, (b) a herd of sheep entering a barn, and (c) grains flowing out a 2D hopper over a vibrated incline. This finding suggests that FIS is a universal phenomenon for active matter passing through a narrowing

From finite geometry exact quantities to (elliptic) scattering amplitudes for spin chains: the 1/2-XYZ

Initially, we derive a nonlinear integral equation for the vacuum counting function of the spin 1/2-XYZ chain in the {\it disordered regime}, thus paralleling similar results by Kl\"umper \cite{KLU}, achieved through a different technique in the {\it antiferroelectric regime}. In terms of the counting function we obtain the usual physical quantities, like the energy and the transfer matrix (eigenvalues). Then, we introduce a double scaling limit which appears to describe the sine-Gordon theory on cylindrical geometry, so generalising famous results in the plane by Luther \cite{LUT} and Johnson et al. \cite{JKM}. Furthermore, after extending the nonlinear integral equation to excitations, we derive scattering amplitudes involving solitons/antisolitons first, and bound states later. The latter case comes out as manifestly related to the Deformed Virasoro Algebra of Shiraishi et al. \cite{SKAO}. Although this nonlinear integral equations framework was contrived to deal with finite geometries, we prove it to be effective for discovering or rediscovering S-matrices. As a particular example, we prove that this unique model furnishes explicitly two S-matrices, proposed respectively by Zamolodchikov \cite{ZAMe} and Lukyanov-Mussardo-Penati \cite{LUK, MP} as plausible scattering description of unknown integrable field theories.Comment: Article, 41 pages, Late

Hybrid Monte Carlo algorithm for the Double Exchange Model

The Hybrid Monte Carlo algorithm is adapted to the simulation of a system of classical degrees of freedom coupled to non self-interacting lattices fermions. The diagonalization of the Hamiltonian matrix is avoided by introducing a path-integral formulation of the problem, in $d+1$ Euclidean space-time. A perfect action formulation allows to work on the continuum euclidean time, without need for a Trotter-Suzuki extrapolation. To demonstrate the feasibility of the method we study the Double Exchange Model in three dimensions. The complexity of the algorithm grows only as the system volume, allowing to simulate in lattices as large as $16^3$ on a personal computer. We conclude that the second order paramagnetic-ferromagnetic phase transition of Double Exchange Materials close to half-filling belongs to the Universality Class of the three-dimensional classical Heisenberg model.Comment: 20 pages plus 4 postscript figure

Nature of the spin-glass phase at experimental length scales

We present a massive equilibrium simulation of the three-dimensional Ising spin glass at low temperatures. The Janus special-purpose computer has allowed us to equilibrate, using parallel tempering, L=32 lattices down to T=0.64 Tc. We demonstrate the relevance of equilibrium finite-size simulations to understand experimental non-equilibrium spin glasses in the thermodynamical limit by establishing a time-length dictionary. We conclude that non-equilibrium experiments performed on a time scale of one hour can be matched with equilibrium results on L=110 lattices. A detailed investigation of the probability distribution functions of the spin and link overlap, as well as of their correlation functions, shows that Replica Symmetry Breaking is the appropriate theoretical framework for the physically relevant length scales. Besides, we improve over existing methodologies to ensure equilibration in parallel tempering simulations.Comment: 48 pages, 19 postscript figures, 9 tables. Version accepted for publication in the Journal of Statistical Mechanic

Simulating spin systems on IANUS, an FPGA-based computer

We describe the hardwired implementation of algorithms for Monte Carlo simulations of a large class of spin models. We have implemented these algorithms as VHDL codes and we have mapped them onto a dedicated processor based on a large FPGA device. The measured performance on one such processor is comparable to O(100) carefully programmed high-end PCs: it turns out to be even better for some selected spin models. We describe here codes that we are currently executing on the IANUS massively parallel FPGA-based system.Comment: 19 pages, 8 figures; submitted to Computer Physics Communication