The lifespan method as a tool to study criticality in absorbing-state phase transitions

In a recent work, a new numerical method (the lifespan method) has been introduced to study the critical properties of epidemic processes on complex networks [Phys. Rev. Lett. \textbf{111}, 068701 (2013)]. Here, we present a detailed analysis of the viability of this method for the study of the critical properties of generic absorbing-state phase transitions in lattices. Focusing on the well understood case of the contact process, we develop a finite-size scaling theory to measure the critical point and its associated critical exponents. We show the validity of the method by studying numerically the contact process on a one-dimensional lattice and comparing the findings of the lifespan method with the standard quasi-stationary method. We find that the lifespan method gives results that are perfectly compatible with those of quasi-stationary simulations and with analytical results. Our observations confirm that the lifespan method is a fully legitimate tool for the study of the critical properties of absorbing phase transitions in regular lattices

Critical properties of Ising model on Sierpinski fractals. A finite size scaling analysis approach

The present paper focuses on the order-disorder transition of an Ising model on a self-similar lattice. We present a detailed numerical study, based on the Monte Carlo method in conjunction with the finite size scaling method, of the critical properties of the Ising model on some two dimensional deterministic fractal lattices with different Hausdorff dimensions. Those with finite ramification order do not display ordered phases at any finite temperature, whereas the lattices with infinite connectivity show genuine critical behavior. In particular we considered two Sierpinski carpets constructed using different generators and characterized by Hausdorff dimensions d_H=log 8/log 3 = 1.8927.. and d_H=log 12/log 4 = 1.7924.., respectively. The data show in a clear way the existence of an order-disorder transition at finite temperature in both Sierpinski carpets. By performing several Monte Carlo simulations at different temperatures and on lattices of increasing size in conjunction with a finite size scaling analysis, we were able to determine numerically the critical exponents in each case and to provide an estimate of their errors. Finally we considered the hyperscaling relation and found indications that it holds, if one assumes that the relevant dimension in this case is the Hausdorff dimension of the lattice.Comment: 21 pages, 7 figures; a new section has been added with results for a second fractal; there are other minor change

Critical behaviour of the Rouse model for gelling polymers

It is shown that the traditionally accepted "Rouse values" for the critical exponents at the gelation transition do not arise from the Rouse model for gelling polymers. The true critical behaviour of the Rouse model for gelling polymers is obtained from spectral properties of the connectivity matrix of the fractal clusters that are formed by the molecules. The required spectral properties are related to the return probability of a "blind ant"-random walk on the critical percolating cluster. The resulting scaling relations express the critical exponents of the shear-stress-relaxation function, and hence those of the shear viscosity and of the first normal stress coefficient, in terms of the spectral dimension $d_{s}$ of the critical percolating cluster and the exponents $\sigma$ and $\tau$ of the cluster-size distribution.Comment: 9 pages, slightly extended version, to appear in J. Phys.

Theory of continuum percolation II. Mean field theory

I use a previously introduced mapping between the continuum percolation model and the Potts fluid to derive a mean field theory of continuum percolation systems. This is done by introducing a new variational principle, the basis of which has to be taken, for now, as heuristic. The critical exponents obtained are $\beta= 1$, $\gamma= 1$ and $\nu = 0.5$, which are identical with the mean field exponents of lattice percolation. The critical density in this approximation is \rho_c = 1/\ve where \ve = \int d \x \, p(\x) \{ \exp [- v(\x)/kT] - 1 \}. p(\x) is the binding probability of two particles separated by \x and v(\x) is their interaction potential.Comment: 25 pages, Late

Griffiths singularities in the two dimensional diluted Ising model

We study numerically the probability distribution of the Yang-Lee zeroes inside the Griffiths phase for the two dimensional site diluted Ising model and we check that the shape of this distribution is that predicted in previous analytical works. By studying the finite size scaling of the averaged smallest zero at the phase transition we extract, for two values of the dilution, the anomalous dimension, $\eta$, which agrees very well with the previous estimated values.Comment: 11 pages and 4 figures, some minor changes in Fig. 4, available at http://chimera.roma1.infn.it/index_papers_complex.htm

A quantum Monte Carlo algorithm realizing an intrinsic relaxation

We propose a new quantum Monte Carlo algorithm which realizes a relaxation intrinsic to the original quantum system. The Monte Carlo dynamics satisfies the dynamic scaling relation $\tau\sim \xi^z$ and is independent of the Trotter number. Finiteness of the Trotter number just appears as the finite-size effect. An infinite Trotter number version of the algorithm is also formulated, which enables us to observe a true relaxation of the original system. The strategy of the algorithm is a compromise between the conventional worldline local flip and the modern cluster loop flip. It is a local flip in the real-space direction and is a cluster flip in the Trotter direction. The new algorithm is tested by the transverse-field Ising model in two dimensions. An accurate phase diagram is obtained.Comment: 9 pages, 4 figure

Percolation and cluster Monte Carlo dynamics for spin models

A general scheme for devising efficient cluster dynamics proposed in a previous letter [Phys.Rev.Lett. 72, 1541 (1994)] is extensively discussed. In particular the strong connection among equilibrium properties of clusters and dynamic properties as the correlation time for magnetization is emphasized. The general scheme is applied to a number of frustrated spin model and the results discussed.Comment: 17 pages LaTeX + 16 figures; will appear in Phys. Rev.

Quantum Phase Transition of Randomly-Diluted Heisenberg Antiferromagnet on a Square Lattice

Ground-state magnetic properties of the diluted Heisenberg antiferromagnet on a square lattice are investigated by means of the quantum Monte Carlo method with the continuous-time loop algorithm. It is found that the critical concentration of magnetic sites is independent of the spin size S, and equal to the two-dimensional percolation threshold. However, the existence of quantum fluctuations makes the critical exponents deviate from those of the classical percolation transition. Furthermore, we found that the transition is not universal, i.e., the critical exponents significantly depend on S.Comment: RevTeX, 4 pages including 5 EPS figure

Monte Carlo computation of correlation times of independent relaxation modes at criticality

We investigate aspects of universality of Glauber critical dynamics in two dimensions. We compute the critical exponent $z$ and numerically corroborate its universality for three different models in the static Ising universality class and for five independent relaxation modes. We also present evidence for universality of amplitude ratios, which shows that, as far as dynamic behavior is concerned, each model in a given universality class is characterized by a single non-universal metric factor which determines the overall time scale. This paper also discusses in detail the variational and projection methods that are used to compute relaxation times with high accuracy