BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Robert M. Ziff

Robert M. Ziff appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.dis-nnmath-phmath.MPcond-mat.mtrl-scicond-mat.softhep-latmath.PRphysics.comp-phBiological Physicsphysics.chem-phphysics.soc-phSocial and Information Networks

Catalog footprint

What is connected

44works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2026arXiv

Bond percolation in distorted square and triangular lattices

This article presents a Monte Carlo study on bond percolation in distorted square and triangular lattices. The distorted lattices are generated by dislocating the sites from their regular positions. The amount and direction of the dislocations are random, but can be tuned by the distortion parameter $α$. Once the sites are dislocated, the bond lengths $δ$ between the nearest neighbors change. A bond can only be occupied if its bond length is less than a threshold value called the connection threshold $d$. It is observed that when the connection threshold is greater than the lattice constant (assumed to be $1$), the bond percolation threshold $p_\mathrm{b}$ always increases with distortion. For $d\le 1$, no spanning configuration is found for the square lattice when the lattice is distorted, even very slightly. On the other hand, the triangular lattice not only spans for $d\le 1$, it also shows a decreasing trend for $p_\mathrm{b}$ in the low-$α$ range. These variation patterns have been linked with the average coordination numbers of the distorted lattices. A critical value $d_\mathrm{c}$ for the connection threshold has been defined as the value of $d$ below which no spanning configuration can be found even after occupying all the bonds satisfying the connection criterion $δ\le d$. The behavior of $d_\mathrm{c}(α)$ is markedly different for the two lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Site and bond percolation on four-dimensional simple hypercubic lattices with extended neighborhoods

The asymptotic behavior of the percolation threshold $p_c$ and its dependence upon coordination number $z$ is investigated for both site and bond percolation on four-dimensional lattices with compact extended neighborhoods. Simple hypercubic lattices with neighborhoods up to 9th nearest neighbors are studied to high precision by means of Monte-Carlo simulations based upon a single-cluster growth algorithm. For site percolation, an asymptotic analysis confirms the predicted behavior $zp_c \sim 16 η_c = 2.086$ for large $z$, and finite-size corrections are accounted for by forms $p_c \sim 16 η_c/(z+b)$ and $p_c \sim 1- \exp(-16 η_c/z)$ where $η_c \approx 0.1304$ is the continuum percolation threshold of four-dimensional hyperspheres. For bond percolation, the finite-$z$ correction is found to be consistent with the prediction of Frei and Perkins, $zp_{c} - 1 \sim a_{1} (\ln z)/z$, although the behavior $zp_{c} - 1 \sim a_1 z^{-3/4}$ cannot be ruled out.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Universal Critical Behavior of Percolation in Orientationally Ordered Janus Particles and Other Anisotropic Systems

We combine percolation theory and Monte Carlo simulation to study in two dimensions the connectivity of an equilibrium lattice model of interacting Janus disks which self-assemble into an orientationally ordered stripe phase at low temperature. As the patch size is increased or the temperature is lowered, clusters of patch-connected disks grow, and a percolating cluster emerges at a threshold. In the stripe phase, the critical clusters extend longer in the direction parallel to the stripes than in the perpendicular direction, and percolation is thus anisotropic. It is found that the critical behavior of percolation in the Janus system is consistent with that of standard isotropic percolation, when an appropriate spatial rescaling is made. The rescaling procedure can be applied to understand other anisotropic systems, such as the percolation of aligned rigid rods and of the $q$-state Potts model with anisotropic interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Universal behavior of site and bond percolation thresholds on regular lattices with compact extended-range neighborhoods in 2 and 3 dimensions

Extended-range percolation on various regular lattices, including all eleven Archimedean lattices in two dimensions, and the simple cubic (SC), body-centered cubic (BCC), and face-centered cubic (FCC) lattices in three dimensions, is investigated. In two dimensions, correlations between coordination number $z$ and site thresholds $p_c$ for Archimedean lattices up to 10th nearest neighbors (NN) are seen by plotting $z$ versus $1/p_{c}$ and $z$ versus $-1/\ln(1-p_c)$, using the data of d'Iribarne et al. [J. Phys. A 32:2611, 1999] and others. The results show that all the plots overlap on a line with a slope consistent with the theoretically predicted asymptotic value of $zp_{c} \sim 4 η_c = 4.51235$, where $η_c$ is the continuum threshold for disks. In three dimensions, precise site and bond thresholds for BCC and FCC lattices with 2nd and 3rd NN, and bond thresholds for the SC lattice with up to the 13th NN, are obtained by Monte-Carlo simulations, using an efficient single-cluster growth method. For site percolation, the values of thresholds for different types of lattices with compact neighborhoods also collapse together, and linear fitting is consistent with the predicted value of $zp_{c} \sim 8 η_c = 2.7351$, where $η_c$ is the continuum threshold for spheres. For bond percolation, Bethe-lattice behavior $p_c = 1/(z-1)$ is expected to hold for large $z$, and the finite-$z$ correction is confirmed to satisfy $zp_{c} - 1 \sim a_{1}z^{-x}$, with $x=2/3$ for three dimensions as predicted by Frei and Perkins [Electron. J. Probab. 21:56, 2016] and by Xu et al. [Phys. Rev. E, 103:022127, 2021]. Our analysis indicates that for compact neighborhoods, the asymptotic behavior of $zp_{c}$ is universal, depending only upon the dimension of the system and whether site or bond percolation, but not upon the type of lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bond percolation between $k$ separated points on a square lattice

We consider a percolation process in which $k$ points separated by a distance proportional to system size $L$ simultaneously connect together ($k>1$), or a single point at the center of a system connects to the boundary ($k=1$), through adjacent connected points of a single cluster. These processes yield new thresholds $\overline p_{ck}$ defined as the average value of $p$ at which the desired connections first occur. These thresholds are not sharp as the distribution of values of $p_{ck}$ for individual samples remains broad in the limit of $L \to \infty$. We study $\overline p_{ck}$ for bond percolation on the square lattice, and find that $\overline p_{ck}$ are above the normal percolation threshold $p_c = 1/2$ and represent specific supercritical states. The $\overline p_{ck}$ can be related to integrals over powers of the function $P_\infty(p)$ equal to the probability a point is connected to the infinite cluster; we find numerically from both direct simulations and from measurements of $P_\infty(p)$ on $L\times L$ systems that, for $L \to \infty$, $\overline p_{c1} = 0.51755(5)$, $\overline p_{c2} = 0.53219(5)$, $\overline p_{c3} = 0.54456(5)$, and $\overline p_{c4} = 0.55527(5).$ The percolation thresholds $\overline p_{ck}$ remain the same, even when the $k$ points are randomly selected within the lattice. We show that the finite-size corrections scale as $L^{-1/ν_k}$ where $ν_k = ν/(k β+1)$, with $β=5/36$ and $ν=4/3$ being the ordinary percolation critical exponents, so that $ν_1= 48/41$, $ν_2 = 24/23$, $ν_3 = 16/17$, $ν_4 = 6/7$, etc. We also study three-point correlations in the system, and show how for $p>p_c$, the correlation ratio goes to 1 (no net correlation) as $L \to \infty$, while at $p_c$ it reaches the known value of 1.022.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bond percolation on simple cubic lattices with extended neighborhoods

We study bond percolation on the simple cubic (SC) lattice with various combinations of first, second, third, and fourth nearest-neighbors by Monte Carlo simulation. Using a single-cluster growth algorithm, we find precise values of the bond thresholds. Correlations between percolation thresholds and lattice properties are discussed, and our results show that the percolation thresholds of these and other three-dimensional lattices decrease monotonically with the coordination number $z$ quite accurately according to a power law $p_{c} \sim z^{-a}$, with exponent $a = 1.111$. However, for large $z$, the threshold must approach the Bethe lattice result $p_c = 1/(z-1)$. Fitting our data and data for lattices with additional nearest neighbors, we find $p_c(z-1)=1+1.224 z^{-1/2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Critical percolation on the kagome hypergraph

We study the percolation critical surface of the kagome lattice in which each triangle is allowed an arbitrary connectivity. Using the method of critical polynomials, we find points along this critical surface to high precision. This kagome hypergraph contains many unsolved problems as special cases, including bond percolation on the kagome and $(3,12^2)$ lattices, and site percolation on the hexagonal, or honeycomb, lattice, as well as a single point for which there is an exact solution. We are able to compute enough points along the critical surface to find a very accurate fit, essentially a Taylor series about the exact point, that allows estimations of the critical point of any system that lies on the surface to precision rivaling Monte Carlo and traditional techniques of similar accuracy. We find also that this system sheds light on some of the surprising aspects of the method of critical polynomials, such as why it is so accurate for certain problems, like the kagome and $(3,12^2)$ lattices. The bond percolation critical points of these lattices can be found to 17 and 18 digits, respectively, because they are in close proximity, in a sense that can be made quantitative, to the exact point on the critical surface. We also discuss in detail a parallel implementation of the method which we use here for a few calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Precise bond percolation thresholds on several four-dimensional lattices

We study bond percolation on several four-dimensional (4D) lattices, including the simple (hyper) cubic (SC), the SC with combinations of nearest neighbors and second nearest neighbors (SC-NN+2NN), the body-centered cubic (BCC), and the face-centered cubic (FCC) lattices, using an efficient single-cluster growth algorithm. For the SC lattice, we find $p_c = 0.1601312(2)$, which confirms previous results (based on other methods), and find a new value $p_c=0.035827(1)$ for the SC-NN+2NN lattice, which was not studied previously for bond percolation. For the 4D BCC and FCC lattices, we obtain $p_c=0.074212(1)$ and 0.049517(1), which are substantially more precise than previous values. We also find critical exponents $τ= 2.3135(5)$ and $Ω= 0.40(3)$, consistent with previous numerical results and the recent four-loop series result of Gracey [Phys. Rev. D 92, 025012, (2015)].

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Renormalization group theory of percolation on pseudo-fractal simplicial and cell complexes

Simplicial complexes are gaining increasing scientific attention as they are generalized network structures that can represent the many-body interactions existing in complex systems raging from the brain to high-order social networks. Simplicial complexes are formed by simplicies, such as nodes, links, triangles and so on. Cell complexes further extend these generalized network structures as they are formed by regular polytopes such as squares, pentagons etc. Pseudo-fractal simplicial and cell complexes are a major example of generalized network structures and they can be obtained by gluing $2$-dimensional $m$-polygons ($m=2$ triangles, $m=4$ squares, $m=5$ pentagons, etc.) along their links according to a simple iterative rule. Here we investigate the interplay between the topology of pseudo-fractal simplicial and cell complexes and their dynamics by characterizing the critical properties of link percolation defined on these structures. By using the renormalization group we show that the pseudo-fractal simplicial and cell complexes have a continuous percolation threshold at $p_c=0$. When the pseudo-fractal structure is formed by polygons of the same size $m$, the transition is characterized by an exponential suppression of the order parameter $P_{\infty}$ that depends on the number of sides $m$ of the polygons forming the pseudo-fractal cell complex, i.e., $P_{\infty}\propto p\exp(-α/p^{m-2})$. Here these results are also generalized to random pseudo-fractal cell-complexes formed by polygons of different number of sides $m$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Site percolation on square and simple cubic lattices with extended neighborhoods and their continuum limit

By means of Monte Carlo simulations, we study long-range site percolation on square and simple cubic lattices with various combinations of nearest neighbors, up to the eighth neighbors for the square lattice and the ninth neighbors for the simple cubic lattice. We find precise thresholds for 23 systems using a single-cluster growth algorithm. Site percolation on lattices with compact neighborhoods can be mapped to problems of lattice percolation of extended shapes, such as disks and spheres, and the thresholds can be related to the continuum thresholds $η_c$ for objects of those shapes. This mapping implies $zp_{c} \sim 4 η_c = 4.51235$ in 2D and $zp_{c} \sim 8 η_c = 2.73512$ in 3D for large $z$ for circular and spherical neighborhoods respectively, where $z$ is the coordination number. Fitting our data to the form $p_c = c/(z+b)$ we find good agreement with $c = 2^d η_c$; the constant $b$ represents a finite-$z$ correction term. We also study power-law fits of the thresholds.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

The elastic and directed percolation backbone

We argue that the elastic backbone (EB) (union of shortest paths) on a cylindrical system, recently studied by Sampaio Filho et al. [Phys. Rev. Lett. 120, 175701 (2018)], is in fact the backbone of two-dimensional directed percolation (DP). We simulate the EB on the same system as considered by these authors, and also study the DP backbone directly using an algorithm that allows backbones to be generated in a completely periodic manner. We find that both the EB in the bulk and the DP backbone have a fractal dimension of $d_{b} = d_{B,\rm DP} = 1.681\,02(15)$ at the identical critical point $p_{c,\rm{DP}} \approx 0.705\,485\,22$. We also measure the fractal dimension at the edge of the EB system and for the full DP clusters, and find $d_e = d_{\rm DP} = 1.840 \, 54 (4)$. We argue that those two fractal dimensions follow from the DP exponents as $d_{B,\rm DP} = 2-2β/ν_\parallel = 1.681 \, 07 2 (12) $ and $d_{\rm DP} = 2-β/ν_\parallel = 1.840\, 536 (6)$. Our fractal dimensions differ from the value 1.750(3) found by Sampaio Filho et al., whose value may represent a crossover effect.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

In a search for a shape maximizing packing fraction for two-dimensional random sequential adsorption

Random sequential adsorption (RSA) of various two dimensional objects is studied in order to find a shape which maximizes the saturated packing fraction. This investigation was begun in our previous paper [Cieśla et al., Phys. Chem. Chem. Phys. 17, 24376 (2015)], where the densest packing was studied for smoothed dimers. Here this shape is compared with a smoothed $n$-mers, spherocylinders and ellipses. It is found that the highest packing fraction out of the studied shapes is $0.58405 \pm 0.0001$ and is obtained for ellipses having long-to-short axis ratio of $1.85$, which is also the largest anisotropy among the investigated shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

No-enclave percolation corresponds to holes in the cluster backbone

The no-enclave percolation (NEP) model introduced recently by Sheinman et al. can be mapped to a problem of holes within a standard percolation backbone, and numerical measurements of these holes gives the size-distribution exponent $τ= 1.82(1)$ of the NEP model. An argument is given that $τ=1 + d_B/2 \approx 1.822$ where $d_B$ is the backbone dimension. On the other hand, a model of simple holes within a percolation cluster implies $τ= 1 + d_f/2 = 187/96 \approx 1.948$, where $d_f$ is the fractal dimension of the cluster, and this value is consistent with Sheinman et al.'s experimental results of gel collapse which gives $τ= 1.91(6)$. Both models yield a discontinuous maximum hole size at $p_c$, signifying explosive percolation behavior. At $p_c$, the largest hole fills exactly half the system, due to symmetry. Extensive numerical simulations confirm our results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Percolation in Finite Matching Lattices

We derive an exact, simple relation between the average number of clusters and the wrapping probabilities for two-dimensional percolation. The relation holds for periodic lattices of any size. It generalizes a classical result of Sykes and Essam and it can be used to find exact or very accurate approximations of the critical density. The criterion that follows is related to the criterion Scullard and Jacobsen use to find precise approximate thresholds, and our work provides a new perspective on their approach.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dimer Covering and Percolation Frustration

Covering a graph or a lattice with non-overlapping dimers is a problem that has received considerable interest in areas such as discrete mathematics, statistical physics, chemistry and materials science. Yet, the problem of percolation on dimer-covered lattices has received little attention. In particular, percolation on lattices that are fully covered by non-overlapping dimers has not evidently been considered. Here, we propose a novel procedure for generating random dimer coverings of a given lattice. We then compute the bond percolation threshold on random and ordered coverings of the square and the triangular lattice, on the remaining bonds connecting the dimers. We obtain $p_c=0.367713(2)$ and $p_c=0.235340(1)$ for random coverings of the square and the triangular lattice, respectively. We observe that the percolation frustration induced as a result of dimer covering is larger in the low-coordination-number square lattice. There is also no relationship between the existence of long-range order in a covering of the square lattice, and its percolation threshold. In particular, an ordered covering of the square lattice, denoted by shifted covering in this work, has an unusually low percolation threshold, and is topologically identical to the triangular lattice. This is in contrast to the other ordered dimer coverings considered in this work, which have higher percolation thresholds than the random covering. In the case of the triangular lattice, the percolation thresholds of the ordered and random coverings are very close, suggesting the lack of sensitivity of the percolation threshold to microscopic details of the covering in highly-coordinated networks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Honeycomb Lattices with Defects

In this paper we introduce a variant of the honeycomb lattice in which we create defects by randomly exchanging adjacent bonds, producing a random tiling with a distribution of polygon edges. We study the percolation properties on these lattices as a function of the number of exchanged bonds using a novel computational method. We find the site and bond percolation thresholds are consistent with other three-coordinated lattices with the same standard deviation in the degree distribution of the dual; here we can produce a continuum of lattices with a range of standard deviations in the distribution. These lattices should be useful for modeling other properties of random systems as well as percolation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Percolation on Hypergraphs with Four-Edges

We study percolation on self-dual hypergraphs that contain hyperedges with four bounding vertices, or "four-edges", using three different generators, each containing bonds or sites with three distinct probabilities $p$, $r$, and $t$ connecting the four vertices. We find explicit values of these probabilities that satisfy the self-duality conditions discussed by Bollobás and Riordan. This demonstrates that explicit solutions of the self-duality conditions can be found using generators containing bonds and sites with independent probabilities. These solutions also provide new examples of lattices where exact percolation critical points are known. One of the generators exhibits three distinct criticality solutions ($p$, $r$, $t$). We carry out Monte-Carlo simulations of two of the generators on two different hypergraphs to confirm the critical values. For the case of the hypergraph and uniform generator studied by Wierman et al., we also determine the threshold $p = 0.441374\pm 0.000001$, which falls within the tight bounds that they derived. Furthermore, we consider a generator in which all or none of the vertices can connect, and find a soluble inhomogeneous percolation system that interpolates between site percolation on the union-jack lattice and bond percolation on the square lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shapes for maximal coverage for two-dimensional random sequential adsorption

The random sequential adsorption of various particle shapes is studied in order to determine the influence of particle anisotropy on the saturated random packing. For all tested particles there is an optimal level of anisotropy which maximizes the saturated packing fraction. It is found that a concave shape derived from a dimer of disks gives a packing fraction of 0.5833, which is comparable to the maximum packing fraction of ellipsoids and spherocylinders and higher than any other studied shape. Discussion why this shape is so beneficial for random sequential adsorption is given.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Short-range correlations in percolation at criticality

We derive the critical nearest-neighbor connectivity $g_n$ as $3/4$, $3(7-9p_c^{tri})/[4(5-4p_c^{tri})]$, and $3(2+7p_c^{tri})/[4(5-p_c^{tri})]$ for bond percolation on the square, honeycomb and triangular lattice respectively, where $p_c^{tri}=2\sin(π/18)$ is the percolation threshold for the triangular lattice; and confirm these values via Monte Carlo simulations. On the square lattice, we also numerically determine the critical next-nearest-neighbor connectivity as $g_{nn}=0.687\;500\;0(2)$, which confirms a conjecture by Mitra and Nienhuis in J. Stat. Mech. P10006 (2004), implying the exact value $g_{nn}=11/16$. We also determine the connectivity on a free surface as $g_n^{surf}=0.625\;000\;1(13)$ and conjecture that this value is exactly equal to $5/8$. In addition, we find that at criticality, the connectivities depend on the linear finite size L as $\sim L^{y_t-d}$, and the associated specific-heat-like quantities $C_n$ and $C_{nn}$ scale as $\sim L^{2y_t-d} \ln (L/L_0)$, where $d$ is the lattice dimensionality, $y_t=1/ν$ the thermal renormalization exponent, and $L_0$ a non-universal constant. We provide an explanation of this logarithmic factor in the theoretical framework reported recently by Vasseur et al. in J. Stat. Mech. L07001 (2012).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Percolation crossing probabilities in hexagons: a numerical study

In a recent article, one of the authors used $c=0$ logarithmic conformal field theory to predict crossing-probability formulas for percolation clusters inside a hexagon with free boundary conditions. In this article, we verify these predictions with high-precision computer simulations. Our simulations generate percolation-cluster perimeters with hull walks on a triangular lattice inside a hexagon. Each sample comprises two hull walks, and the order in which these walks strike the bottom and upper left/right sides of the hexagon determines the crossing configuration of the percolation sample. We compare our numerical results with the predicted crossing probabilities, finding excellent agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The critical manifolds of inhomogeneous bond percolation on bow-tie and checkerboard lattices

We give a conditional derivation of the inhomogeneous critical percolation manifold of the bow-tie lattice with five different probabilities, a problem that does not appear at first to fall into any known solvable class. Although our argument is mathematically rigorous only on a region of the manifold, we conjecture that the formula is correct over its entire domain, and we provide a non-rigorous argument for this that employs the negative probability regime of the triangular lattice critical surface. We discuss how the rigorous portion of our result substantially broadens the range of lattices in the solvable class to include certain inhomogeneous and asymmetric bow-tie lattices, and that, if it could be put on a firm foundation, the negative probability portion of our method would extend this class to many further systems, including F.Y. Wu's checkerboard formula for the square lattice. We conclude by showing that this latter problem can in fact be proved using a recent result of Grimmett and Manolescu for isoradial graphs, lending strong evidence in favour of our other conjectured results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cluster pinch-point densities in polygons

In a statistical cluster or loop model such as percolation, or more generally the Potts models or O(n) models, a pinch point is a single bulk point where several distinct clusters or loops touch. In a polygon P harboring such a model in its interior and with 2N sides exhibiting free/fixed side-alternating boundary conditions, "boundary" clusters anchor to the fixed sides of P. At the critical point and in the continuum limit, the density (i.e., frequency of occurrence) of pinch-point events between s distinct boundary clusters at a bulk point w in P is proportional to <psi_1^c(w_1)psi_1^c(w_2)...psi_1^c(w_{2N-1})psi_1^c(w_{2N})Psi_s(w)>_P. The w_i are the vertices of P, psi_1^c is a conformal field theory (CFT) corner one-leg operator, and Psi_s is a CFT bulk 2s-leg operator. In this article, we use the Coulomb gas formalism to construct explicit contour integral formulas for these correlation functions and thereby calculate the density of various pinch-point configurations at arbitrary points in the rectangle, in the hexagon, and for the case s=N, in the 2N-sided polygon at the system's critical point. Explicit formulas for these results are given in terms of algebraic functions or integrals of algebraic functions, particularly Lauricella functions. In critical percolation, the result for s=N=2 gives the density of red bonds between boundary clusters (in the continuum limit) inside a rectangle. We compare our results with high-precision simulations of critical percolation and Ising FK clusters in a rectangle of aspect ratio two and in a regular hexagon and find very good agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Crossing on hyperbolic lattices

We divide the circular boundary of a hyperbolic lattice into four equal intervals, and study the probability of a percolation crossing between an opposite pair, as a function of the bond occupation probability p. We consider the {7,3} (heptagonal), enhanced or extended binary tree (EBT), the EBT-dual, and {5,5} (pentagonal) lattices. We find that the crossing probability increases gradually from zero to one as p increases from the lower p_l to the upper p_u critical values. We find bounds and estimates for the values of p_ l and p_u for these lattices, and identify the self-duality point p* corresponding to where the crossing probability equals 1/2. Comparison is made with recent numerical and theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Crossover from Isotropic to Directed Percolation

We generalize the directed percolation (DP) model by relaxing the strict directionality of DP such that propagation can occur in either direction but with anisotropic probabilities. We denote the probabilities as $p_{\downarrow}= p \cdot p_d$ and $p_{\uparrow}=p \cdot (1-p_d)$, with $p $ representing the average occupation probability and $p_d$ controlling the anisotropy. The Leath-Alexandrowicz method is used to grow a cluster from an active seed site. We call this model with two main growth directions {\em biased directed percolation} (BDP). Standard isotropic percolation (IP) and DP are the two limiting cases of the BDP model, corresponding to $p_d=1/2$ and $p_d=0,1$ respectively. In this work, besides IP and DP, we also consider the $1/2<p_d<1$ region. Extensive Monte Carlo simulations are carried out on the square and the simple-cubic lattices, and the numerical data are analyzed by finite-size scaling. We locate the percolation thresholds of the BDP model for $p_d=0.6$ and 0.8, and determine various critical exponents. These exponents are found to be consistent with those for standard DP. We also determine the renormalization exponent associated with the asymmetric perturbation due to $p_d -1/2 \neq 0$ near IP, and confirm that such an asymmetric scaling field is relevant at IP.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Retention capacity of random surfaces

We introduce a "water retention" model for liquids captured on a random surface with open boundaries, and investigate it for both continuous and discrete surface heights 0, 1, ... n-1, on a square lattice with a square boundary. The model is found to have several intriguing features, including a non-monotonic dependence of the retention on the number of levels in the discrete case: for many n, the retention is counterintuitively greater than that of an n+1-level system. The behavior is explained using percolation theory, by mapping it to a 2-level system with variable probability. Results in 1-dimension are also found.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Shortest-Path Fractal Dimension for Percolation in Two and Three Dimensions

We carry out a high-precision Monte Carlo study of the shortest-path fractal dimension $\dm$ for percolation in two and three dimensions, using the Leath-Alexandrowicz method which grows a cluster from an active seed site. A variety of quantities are sampled as a function of the chemical distance, including the number of activated sites, a measure of the radius, and the survival probability. By finite-size scaling, we determine $\dm = 1.130 77(2)$ and $1.375 6(6)$ in two and three dimensions, respectively. The result in 2D rules out the recently conjectured value $\dm=217/192$ [Phys. Rev. E 81, 020102(R) (2010)].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cluster densities at 2-D critical points in rectangular geometries

Making use of a recent complete calculation of a chiral six-point correlation function C(z) in a rectangle we calculate various quantities of interest for percolation (SLE parameter κ= 6) and many other two-dimensional critical points. In particular, we specify the density at z of critical clusters conditioned to touch either or both vertical sides of the rectangle, with these sides 'wired,' i.e. constrained to be in a single cluster, and the horizontal sides free. These quantities probe the structure of various cluster configurations, including those that contribute to the crossing probability. We first examine the effects of boundary conditions on C for the critical O(n) loop models in both high and low density phases and for both Fortuin-Kasteleyn (FK) and spin clusters in the critical Q-state Potts models. A Coulomb gas analysis then allows us to calculate the cluster densities with various conditionings in terms of the known conformal blocks. Explicit formulas generalizing Cardy's horizontal crossing probability to these models (using previously known results) are also presented. These solutions are employed to generalize previous results demonstrating factorization of higher-order correlation functions to the critical systems mentioned. An explicit formula for the density of critical percolation clusters that cross a rectangle horizontally with free boundary conditions is also given. Simplifications of the hypergeometric functions in our solutions for various models are presented. High precision simulations verify these predictions for percolation and for the Q=2 and 3-state Potts models, including both FK and spin clusters. Our formula for the density of crossing clusters in percolation in open systems is also verified.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Correction-to-scaling exponent for two-dimensional percolation

We show that the correction-to-scaling exponents in two-dimensional percolation are bounded by Omega <= 72/91, omega = D Omega <= 3/2, and Delta_1 = nu omega <= 2, based upon Cardy's result for the critical crossing probability on an annulus. The upper bounds are consistent with many previous measurements of site percolation on square and triangular lattices, and new measurements for bond percolation presented here, suggesting this result is exact. A scaling form evidently applicable to site percolation is also found.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Factorization of correlations in two-dimensional percolation on the plane and torus

Recently, Delfino and Viti have examined the factorization of the three-point density correlation function P_3 at the percolation point in terms of the two-point density correlation functions P_2. According to conformal invariance, this factorization is exact on the infinite plane, such that the ratio R(z_1, z_2, z_3) = P_3(z_1, z_2, z_3) [P_2(z_1, z_2) P_2(z_1, z_3) P_2(z_2, z_3)]^{1/2} is not only universal but also a constant, independent of the z_i, and in fact an operator product expansion (OPE) coefficient. Delfino and Viti analytically calculate its value (1.022013...) for percolation, in agreement with the numerical value 1.022 found previously in a study of R on the conformally equivalent cylinder. In this paper we confirm the factorization on the plane numerically using periodic lattices (tori) of very large size, which locally approximate a plane. We also investigate the general behavior of R on the torus, and find a minimum value of R approx. 1.0132 when the three points are maximally separated. In addition, we present a simplified expression for R on the plane as a function of the SLE parameter kappa.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Results for a critical threshold, the correction-to-scaling exponent and susceptibility amplitude ratio for 2d percolation

We summarize several decades of work in finding values for the percolation threshold p_c for site percolation on the square lattice, the universal correction-to-scaling exponent Omega, and the susceptibility amplitude ratio C^+/C^-, in two dimensions. Recent studies have yielded the precise values p_c = 0.59274602(4), Omega = 72/91 = 0.791, and C^+/C^- = 161.5(1.5), resolving long-standing controversies about the last two quantities and verifying the widely used value p_c = 0.592746 for the first.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tricritical point in explosive percolation

The suitable interpolation between classical percolation and a special variant of explosive percolation enables the explicit realization of a tricritical percolation point. With high-precision simulations of the order parameter and the second moment of the cluster size distribution a fully consistent tricritical scaling scenario emerges yielding the tricritical crossover exponent $1/ϕ_t=1.8\pm0.1$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Computation of nucleation of a non-equilibrium first-order phase transition using a rare-event algorithm

We introduce a new Forward-Flux Sampling in Time (FFST) algorithm to efficiently measure transition times in rare-event processes in non-equilibrium systems, and apply it to study the first-order (discontinuous) kinetic transition in the Ziff-Gulari-Barshad model of catalytic surface reaction. The average time for the transition to take place, as well as both the spinodal and transition points, are clearly found by this method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Critical surfaces for general inhomogeneous bond percolation problems

We present a method of general applicability for finding exact or accurate approximations to bond percolation thresholds for a wide class of lattices. To every lattice we sytematically associate a polynomial, the root of which in $[0,1]$ is the conjectured critical point. The method makes the correct prediction for every exactly solved problem, and comparison with numerical results shows that it is very close, but not exact, for many others. We focus primarily on the Archimedean lattices, in which all vertices are equivalent, but this restriction is not crucial. Some results we find are kagome: $p_c=0.524430...$, $(3,12^2): p_c=0.740423...$, $(3^3,4^2): p_c=0.419615...$, $(3,4,6,4):p_c=0.524821...$, $(4,8^2):p_c=0.676835...$, $(3^2,4,3,4)$: $p_c=0.414120...$ . The results are generally within $10^{-5}$ of numerical estimates. For the inhomogeneous checkerboard and bowtie lattices, errors in the formulas (if they are not exact) are less than $10^{-6}$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Finite-size scaling of the stochastic susceptible-infected-recovered model

The critical behavior of the stochastic susceptible-infected-recovered model on a square lattice is obtained by numerical simulations and finite-size scaling. The order parameter as well as the distribution in the number of recovered individuals is determined as a function of the infection rate for several values of the system size. The analysis around criticality is obtained by exploring the close relationship between the present model and standard percolation theory. The quantity UP, equal to the ratio U between the second moment and the squared first moment of the size distribution multiplied by the order parameter P, is shown to have, for a square system, a universal value 1.0167(1) that is the same as for site and bond percolation, confirming further that the SIR model is also in the percolation class.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the critical behavior of the Susceptible-Infected-Recovered (SIR) model on a square lattice

By means of numerical simulations and epidemic analysis, the transition point of the stochastic, asynchronous Susceptible-Infected-Recovered (SIR) model on a square lattice is found to be c_0=0.1765005(10), where c is the probability a chosen infected site spontaneously recovers rather than tries to infect one neighbor. This point corresponds to an infection/recovery rate of lambda_c = (1-c_0)/c_0 = 4.66571(3) and a net transmissibility of (1-c_0)/(1 + 3 c_0) = 0.538410(2), which falls between the rigorous bounds of the site and bond thresholds. The critical behavior of the model is consistent with the 2-d percolation universality class, but local growth probabilities differ from those of dynamic percolation cluster growth, as is demonstrated explicitly.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Scaling behavior of explosive percolation on the square lattice

Clusters generated by the product-rule growth model of Achlioptas, D'Souza, and Spencer on a two-dimensional square lattice are shown to obey qualitatively different scaling behavior than standard (random growth) percolation. The threshold with unrestricted bond placement (allowing loops) is found precisely using several different criteria based upon both moments and wrapping probabilities, yielding p_c = 0.526565 +/- 0.000005, consistent with the recent result of Radicchi and Fortunato. The correlation-length exponent nu is found to be close to 1. The qualitative difference from regular percolation is shown dramatically in the behavior of the percolation probability P_(infinity) (size of largest cluster), the susceptibility, and of the second moment of finite clusters, where discontinuities appears at the threshold. The critical cluster-size distribution does not follow a consistent power-law for the range of system sizes we study L </- 8192) but may approach a power-law with tau > 2 for larger L.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Barrier Method: A Technique for Calculating Very Long Transition Times

In many dynamical systems there is a large separation of time scales between typical events and "rare" events which can be the cases of interest. Rare-event rates are quite difficult to compute numerically, but they are of considerable practical importance in many fields: for example transition times in chemical physics and extinction times in epidemiology can be very long, but are quite important. We present a very fast numerical technique that can be used to find long transition times (very small rates) in low-dimensional systems, even if they lack detailed balance. We illustrate the method for a bistable non-equilibrium system introduced by Maier and Stein and a two-dimensional (in parameter space) epidemiology model.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Capture of particles undergoing discrete random walks

It is shown that particles undergoing discrete-time jumps in 3D, starting at a distance r0 from the center of an adsorbing sphere of radius R, are captured with probability (R - c sigma)/r0 for r0 much greater than R, where c is related to the Fourier transform of the scaled jump distribution and sigma is the distribution's root-mean square jump length. For particles starting on the surface of the sphere, the asymptotic survival probability is non-zero (in contrast to the case of Brownian diffusion) and has a universal behavior sigma/(R sqrt(6)) depending only upon sigma/R. These results have applications to computer simulations of reaction and aggregation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fractal dimensions of the Q-state Potts model for the complete and external hulls

Fortuin-Kastelyn clusters in the critical $Q$-state Potts model are conformally invariant fractals. We obtain simulation results for the fractal dimension of the complete and external (accessible) hulls for Q=1, 2, 3, and 4, on clusters that wrap around a cylindrical system. We find excellent agreement between these results and theoretical predictions. We also obtain the probability distributions of the hull lengths and maximal heights of the clusters in this geometry and provide a conjecture for their form.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Factorization of percolation density correlation functions for clusters touching the sides of a rectangle

We consider the density at a point z = x + i y of critical percolation clusters that touch the left [P_L(z)], right [P_R(z)], or both [P_{LR}(z)] sides of a rectangular system, with open boundary conditions on the top and bottom. The ratio C(z) = P_{LR}(z) / sqrt[P_L(z) P_R(z) Pi_h], where Pi_h is the probability of horizontal crossing given by Cardy, is a universal function of z and goes to a constant value C_0 = 2^(7/2) 3^(-3/4) pi^(5/2) Gamma(1/3)^(-9/2) = 1.0299268... far from the ends. We observe numerically that C(z) depends upon x but not y for wired b.c., and this result leads to an explicit expression for C(z) via conformal field theory. For the semi-infinite strip we also derive explict expressions for P_L(z), P_R(z), and P_{LR}(z), for both wired and open b.c. Our results enable calculation of the finite-size corrections to the factorization near an isolated anchor point, for the case of clusters anchored at two boundary points. Finally, we present numerical results for a rectangle with periodic b.c. in the horizontal direction, and find that C(z) approaches a constant value C_1 = 1.022.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Re-examination of seven-dimensional site percolation thresholds

Monte Carlo simulations alone could not clarify the corrections to scaling for the size-dependent p_c(L) above the upper critical dimension. Including the previous series estimate for the bulk threshold $p_c(\infty)$ gives preference for the complicated corrections predicted by renormalization group and against the simple 1/L extrapolation. Additional Monte-Carlo simulations using the Leath method corroborate the series result for p_c.

0 citations0 reviewsNo code linkedOpen access