Critical Percolation in High Dimensions

We present Monte Carlo estimates for site and bond percolation thresholds in simple hypercubic lattices with 4 to 13 dimensions. For d<6 they are preliminary, for d >= 6 they are between 20 to 10^4 times more precise than the best previous estimates. This was achieved by three ingredients: (i) simple and fast hashing which allowed us to simulate clusters of millions of sites on computers with less than 500 MB memory; (ii) a histogram method which allowed us to obtain information for several p values from a single simulation; and (iii) a new variance reduction technique which is especially efficient at high dimensions where it reduces error bars by a factor up to approximately 30 and more. Based on these data we propose a new scaling law for finite cluster size corrections.Comment: 5 pages including figures, RevTe

Universal Formulae for Percolation Thresholds

A power law is postulated for both site and bond percolation thresholds. The formula writes $p_c=p_0[(d-1)(q-1)]^{-a}d^{\ b}$, where $d$ is the space dimension and $q$ the coordination number. All thresholds up to $d\rightarrow \infty$ are found to belong to only three universality classes. For first two classes $b=0$ for site dilution while $b=a$ for bond dilution. The last one associated to high dimensions is characterized by $b=2a-1$ for both sites and bonds. Classes are defined by a set of value for $\{p_0; \ a\}$. Deviations from available numerical estimates at $d \leq 7$ are within $\pm 0.008$ and $\pm 0.0004$ for high dimensional hypercubic expansions at $d \geq 8$. The formula is found to be also valid for Ising critical temperatures.Comment: 11 pages, latex, 3 figures not include

Critical Exponent for the Density of Percolating Flux

This paper is a study of some of the critical properties of a simple model for flux. The model is motivated by gauge theory and is equivalent to the Ising model in three dimensions. The phase with condensed flux is studied. This is the ordered phase of the Ising model and the high temperature, deconfined phase of the gauge theory. The flux picture will be used in this phase. Near the transition, the density is low enough so that flux variables remain useful. There is a finite density of finite flux clusters on both sides of the phase transition. In the deconfined phase, there is also an infinite, percolating network of flux with a density that vanishes as $T \rightarrow T_{c}^{+}$. On both sides of the critical point, the nonanalyticity in the total flux density is characterized by the exponent $(1-\alpha)$. The main result of this paper is a calculation of the critical exponent for the percolating network. The exponent for the density of the percolating cluster is $\zeta = (1-\alpha) - (\varphi-1)$. The specific heat exponent $\alpha$ and the crossover exponent $\varphi$ can be computed in the $\epsilon$-expansion. Since $\zeta < (1-\alpha)$, the variation in the separate densities is much more rapid than that of the total. Flux is moving from the infinite cluster to the finite clusters much more rapidly than the total density is decreasing.Comment: 20 pages, no figures, Latex/Revtex 3, UCD-93-2

Site percolation and random walks on d-dimensional Kagome lattices

The site percolation problem is studied on d-dimensional generalisations of the Kagome' lattice. These lattices are isotropic and have the same coordination number q as the hyper-cubic lattices in d dimensions, namely q=2d. The site percolation thresholds are calculated numerically for d= 3, 4, 5, and 6. The scaling of these thresholds as a function of dimension d, or alternatively q, is different than for hypercubic lattices: p_c ~ 2/q instead of p_c ~ 1/(q-1). The latter is the Bethe approximation, which is usually assumed to hold for all lattices in high dimensions. A series expansion is calculated, in order to understand the different behaviour of the Kagome' lattice. The return probability of a random walker on these lattices is also shown to scale as 2/q. For bond percolation on d-dimensional diamond lattices these results imply p_c ~ 1/(q-1).Comment: 11 pages, LaTeX, 8 figures (EPS format), submitted to J. Phys.

Hard squares with negative activity

We show that the hard-square lattice gas with activity z= -1 has a number of remarkable properties. We conjecture that all the eigenvalues of the transfer matrix are roots of unity. They fall into groups (``strings'') evenly spaced around the unit circle, which have interesting number-theoretic properties. For example, the partition function on an M by N lattice with periodic boundary condition is identically 1 when M and N are coprime. We provide evidence for these conjectures from analytical and numerical arguments.Comment: 8 page

Square lattice site percolation at increasing ranges of neighbor interactions

We report site percolation thresholds for square lattice with neighbor interactions at various increasing ranges. Using Monte Carlo techniques we found that nearest neighbors (N$^2$), next nearest neighbors (N$^3$), next next nearest neighbors (N$^4$) and fifth nearest neighbors (N$^6$) yield the same $p_c=0.592...$. At odds, fourth nearest neighbors (N$^5$) give $p_c=0.298...$. These results are given an explanation in terms of symmetry arguments. We then consider combinations of various ranges of interactions with (N$^2$+N$^3$), (N$^2$+N$^4$), (N$^2$+N$^3$+N$^4$) and (N$^2$+N$^5$). The calculated associated thresholds are respectively $p_c=0.407..., 0.337..., 0.288..., 0.234...$. The existing Galam--Mauger universal formula for percolation thresholds does not reproduce the data showing dimension and coordination number are not sufficient to build a universal law which extends to complex lattices.Comment: 4 pages, revtex

Relaxation properties in a lattice gas model with asymmetrical particles

We study the relaxation process in a two-dimensional lattice gas model, where the interactions come from the excluded volume. In this model particles have three arms with an asymmetrical shape, which results in geometrical frustration that inhibits full packing. A dynamical crossover is found at the arm percolation of the particles, from a dynamical behavior characterized by a single step relaxation above the transition, to a two-step decay below it. Relaxation functions of the self-part of density fluctuations are well fitted by a stretched exponential form, with a $\beta$ exponent decreasing when the temperature is lowered until the percolation transition is reached, and constant below it. The structural arrest of the model seems to happen only at the maximum density of the model, where both the inverse diffusivity and the relaxation time of density fluctuations diverge with a power law. The dynamical non linear susceptibility, defined as the fluctuations of the self-overlap autocorrelation, exhibits a peak at some characteristic time, which seems to diverge at the maximum density as well.Comment: 7 pages and 9 figure

Complex-Temperature Singularities in the d=2d=2 Ising Model. III. Honeycomb Lattice

We study complex-temperature properties of the uniform and staggered susceptibilities $\chi$ and $\chi^{(a)}$ of the Ising model on the honeycomb lattice. From an analysis of low-temperature series expansions, we find evidence that $\chi$ and $\chi^{(a)}$ both have divergent singularities at the point $z=-1 \equiv z_{\ell}$ (where $z=e^{-2K}$), with exponents $\gamma_{\ell}'= \gamma_{\ell,a}'=5/2$. The critical amplitudes at this singularity are calculated. Using exact results, we extract the behaviour of the magnetisation $M$ and specific heat $C$ at complex-temperature singularities. We find that, in addition to its zero at the physical critical point, $M$ diverges at $z=-1$ with exponent $\beta_{\ell}=-1/4$, vanishes continuously at $z=\pm i$ with exponent $\beta_s=3/8$, and vanishes discontinuously elsewhere along the boundary of the complex-temperature ferromagnetic phase. $C$ diverges at $z=-1$ with exponent $\alpha_{\ell}'=2$ and at $v=\pm i/\sqrt{3}$ (where $v = \tanh K$) with exponent $\alpha_e=1$, and diverges logarithmically at $z=\pm i$. We find that the exponent relation $\alpha'+2\beta+\gamma'=2$ is violated at $z=-1$; the right-hand side is 4 rather than 2. The connections of these results with complex-temperature properties of the Ising model on the triangular lattice are discussed.Comment: 22 pages, latex, figures appended after the end of the text as a compressed, uuencoded postscript fil

Trypanosoma cruzi IIc: phylogenetic and phylogeographic insights from sequence and microsatellite analysis and potential impact on emergent Chagas disease.

Trypanosoma cruzi, the etiological agent of Chagas disease, is highly genetically diverse. Numerous lines of evidence point to the existence of six stable genetic lineages or DTUs: TcI, TcIIa, TcIIb, TcIIc, TcIId, and TcIIe. Molecular dating suggests that T. cruzi is likely to have been an endemic infection of neotropical mammalian fauna for many millions of years. Here we have applied a panel of 49 polymorphic microsatellite markers developed from the online T. cruzi genome to document genetic diversity among 53 isolates belonging to TcIIc, a lineage so far recorded almost exclusively in silvatic transmission cycles but increasingly a potential source of human infection. These data are complemented by parallel analysis of sequence variation in a fragment of the glucose-6-phosphate isomerase gene. New isolates confirm that TcIIc is associated with terrestrial transmission cycles and armadillo reservoir hosts, and demonstrate that TcIIc is far more widespread than previously thought, with a distribution at least from Western Venezuela to the Argentine Chaco. We show that TcIIc is truly a discrete T. cruzi lineage, that it could have an ancient origin and that diversity occurs within the terrestrial niche independently of the host species. We also show that spatial structure among TcIIc isolates from its principal host, the armadillo Dasypus novemcinctus, is greater than that among TcI from Didelphis spp. opossums and link this observation to differences in ecology of their respective niches. Homozygosity in TcIIc populations and some linkage indices indicate the possibility of recombination but cannot yet be effectively discriminated from a high genome-wide frequency of gene conversion. Finally, we suggest that the derived TcIIc population genetic data have a vital role in determining the origin of the epidemiologically important hybrid lineages TcIId and TcIIe

Influence of extended dynamics on phase transitions in a driven lattice gas

Monte Carlo simulations and dynamical mean-field approximations are performed to study the phase transition in a driven lattice gas with nearest-neighbor exclusion on a square lattice. A slight extension of the microscopic dynamics with allowing the next-nearest-neighbor hops results in dramatic changes. Instead of the phase separation into high- and low-density regions in the stationary state the system exhibits a continuous transition belonging to the Ising universality class for any driving. The relevant features of phase diagram are reproduced by an improved mean-field analysis.Comment: 3 pages, 3 figure