Ferromagnetic phase transition for the spanning-forest model (q \to 0 limit of the Potts model) in three or more dimensions

We present Monte Carlo simulations of the spanning-forest model (q \to 0 limit of the ferromagnetic Potts model) in spatial dimensions d=3,4,5. We show that, in contrast to the two-dimensional case, the model has a "ferromagnetic" second-order phase transition at a finite positive value w_c. We present numerical estimates of w_c and of the thermal and magnetic critical exponents. We conjecture that the upper critical dimension is 6.Comment: LaTex2e, 4 pages; includes 6 Postscript figures; Version 2 has expanded title as published in PR

Cluster simulations of loop models on two-dimensional lattices

We develop cluster algorithms for a broad class of loop models on two-dimensional lattices, including several standard O(n) loop models at n \ge 1. We show that our algorithm has little or no critical slowing-down when 1 \le n \le 2. We use this algorithm to investigate the honeycomb-lattice O(n) loop model, for which we determine several new critical exponents, and a square-lattice O(n) loop model, for which we obtain new information on the phase diagram.Comment: LaTex2e, 4 pages; includes 1 table and 2 figures. Totally rewritten in version 2, with new theory and new data. Version 3 as published in PR

Dynamic critical behavior of the Chayes-Machta-Swendsen-Wang algorithm

We study the dynamic critical behavior of the Chayes-Machta dynamics for the Fortuin-Kasteleyn random-cluster model, which generalizes the Swendsen-Wang dynamics for the q-state Potts model to noninteger q, in two and three spatial dimensions, by Monte Carlo simulation. We show that the Li-Sokal bound z \ge \alpha/\nu is close to but probably not sharp in d=2, and is far from sharp in d=3, for all q. The conjecture z \ge \beta/\nu is false (for some values of q) in both d=2 and d=3.Comment: Revtex4, 4 pages including 4 figure

Self-avoiding walks on scale-free networks

Several kinds of walks on complex networks are currently used to analyze search and navigation in different systems. Many analytical and computational results are known for random walks on such networks. Self-avoiding walks (SAWs) are expected to be more suitable than unrestricted random walks to explore various kinds of real-life networks. Here we study long-range properties of random SAWs on scale-free networks, characterized by a degree distribution $P(k) \sim k^{-\gamma}$. In the limit of large networks (system size $N \to \infty$), the average number $s_n$ of SAWs starting from a generic site increases as $\mu^n$, with $\mu = / - 1$. For finite $N$, $s_n$ is reduced due to the presence of loops in the network, which causes the emergence of attrition of the paths. For kinetic growth walks, the average maximum length, $$, increases as a power of the system size: $\sim N^{\alpha}$, with an exponent $\alpha$ increasing as the parameter $\gamma$ is raised. We discuss the dependence of $\alpha$ on the minimum allowed degree in the network. A similar power-law dependence is found for the mean self-intersection length of non-reversal random walks. Simulation results support our approximate analytical calculations.Comment: 9 pages, 7 figure

Critical speeding-up in a local dynamics for the random-cluster model

We study the dynamic critical behavior of the local bond-update (Sweeny) dynamics for the Fortuin-Kasteleyn random-cluster model in dimensions d=2,3, by Monte Carlo simulation. We show that, for a suitable range of q values, the global observable S_2 exhibits "critical speeding-up": it decorrelates well on time scales much less than one sweep, so that the integrated autocorrelation time tends to zero as the critical point is approached. We also show that the dynamic critical exponent z_{exp} is very close (possibly equal) to the rigorous lower bound \alpha/\nu, and quite possibly smaller than the corresponding exponent for the Chayes-Machta-Swendsen-Wang cluster dynamics.Comment: LaTex2e/revtex4, 4 pages, includes 5 figure

Grassmann Integral Representation for Spanning Hyperforests

Given a hypergraph G, we introduce a Grassmann algebra over the vertex set, and show that a class of Grassmann integrals permits an expansion in terms of spanning hyperforests. Special cases provide the generating functions for rooted and unrooted spanning (hyper)forests and spanning (hyper)trees. All these results are generalizations of Kirchhoff's matrix-tree theorem. Furthermore, we show that the class of integrals describing unrooted spanning (hyper)forests is induced by a theory with an underlying OSP(1|2) supersymmetry.Comment: 50 pages, it uses some latex macros. Accepted for publication on J. Phys.

Failure of vaccination to prevent outbreaks of foot-and-mouth disease

Outbreaks of foot-and-mouth disease persist in dairy cattle herds in Saudi Arabia despite revaccination at intervals of 4-6 months. Vaccine trials provide data on antibody responses following vaccination. Using this information we developed a mathematical model of the decay of protective antibodies with which we estimated the fraction of susceptible animals at a given time after vaccination. The model describes the data well, suggesting over 95% take with an antibody half-life of 43 days. Farm records provided data on the time course of five outbreaks. We applied a 'SLIR' epidemiological model to these data, fitting a single parameter representing disease transmission rate. The analysis provides estimates of the basic reproduction number R(0), which may exceed 70 in some cases. We conclude that the critical intervaccination interval which would provide herd immunity against FMDV is unrealistically short, especially for heterologous challenge. We suggest that it may not be possible to prevent foot-and-mouth disease outbreaks on these farms using currently available vaccines

Projected single-spin flip dynamics in the Ising Model

We study transition matrices for projected dynamics in the energy-magnetization space, magnetization space and energy space. Several single spin flip dynamics are considered such as the Glauber and Metropolis canonical ensemble dynamics and the Metropolis dynamics for three multicanonical ensembles: the flat energy-magnetization histogram, the flat energy histogram and the flat magnetization histogram. From the numerical diagonalization of the matrices for the projected dynamics we obtain the sub-dominant eigenvalue and the largest relaxation times for systems of varying size. Although, the projected dynamics is an approximation to the full state space dynamics comparison with some available results, obtained by other authors, shows that projection in the magnetization space is a reasonably accurate method to study the scaling of relaxation times with system size. The transition matrices for arbitrary single-spin flip dynamics are obtained from a single Monte-Carlo estimate of the infinite temperature transition-matrix, for each system size, which makes the method an efficient tool to evaluate the relative performance of any arbitrary local spin-flip dynamics. We also present new results for appropriately defined average tunnelling times of magnetization and compute their finite-size scaling exponents that we compare with results of energy tunnelling exponents available for the flat energy histogram multicanonical ensemble.Comment: 23 pages and 6 figure

Conformations, Transverse Fluctuations and Crossover Dynamics of a Semi-Flexible Chain in Two Dimensions

We present a unified scaling description for the dynamics of monomers of a semiflexible chain under good solvent condition in the free draining limit. We consider both the cases where the contour length $L$ is comparable to the persistence length $\ell_p$ and the case $L\gg \ell_p$. Our theory captures the early time monomer dynamics of a stiff chain characterized by $t^{3/4}$ dependence for the mean square displacement(MSD) of the monomers, but predicts a first crossover to the Rouse regime of $t^{2\nu/{1+2\nu}}$ for $\tau_1 \sim \ell_p^3$, and a second crossover to the purely diffusive dynamics for the entire chain at $\tau_2 \sim L^{5/2}$. We confirm the predictions of this scaling description by studying monomer dynamics of dilute solution of semi-flexible chains under good solvent conditions obtained from our Brownian dynamics (BD) simulation studies for a large choice of chain lengths with number of monomers per chain N = 16 - 2048 and persistence length $\ell_p = 1 - 500$ Lennard-Jones (LJ) units. These BD simulation results further confirm the absence of Gaussian regime for a 2d swollen chain from the slope of the plot of $\langle R_N^2 \rangle/2L \ell_p \sim L/\ell_p$ which around $L/\ell_p \sim 1$ changes suddenly from $\left(L/\ell_p \right) \rightarrow \left(L/\ell_p \right)^{0.5}$, also manifested in the power law decay for the bond autocorrelation function disproving the validity of the WLC in 2d. We further observe that the normalized transverse fluctuations of the semiflexible chains for different stiffness $\sqrt{\langle l_{\bot}^2\rangle}/L$ as a function of renormalized contour length $L/\ell_p$ collapse on the same master plot and exhibits power law scaling $\sqrt{\langle l_{\bot}^2\rangle}/L \sim (L/\ell_p)^\eta$ at extreme limits, where $\eta = 0.5$ for extremely stiff chains ($L/\ell_p \gg 1$), and $\eta = -0.25$ for fully flexible chains.Comment: 14 pages, 18 figure