We study numerically the geometrical properties of minimally weighted paths that appear in the negative-weight percolation (NWP) model on two-dimensional lattices assuming a combination of periodic and free boundary conditions (BCs). Each realization of the disorder consists of a random fraction 1-rho of bonds with unit strength and a fraction rho of bond strengths drawn from a Gaussian distribution with zero mean and unit width. For each such sample, the path is forced to span the lattice along the direction with the free BCs. The path and a set of negatively weighted loops form a ground state (GS). A ground state on such a lattice can be determined performing a non-trivial transformation of the original graph and applying sophisticated matching algorithms. Here we examine whether the geometrical properties of the paths are in accordance with predictions of Schramm-Loewner evolution (SLE). Measuring the fractal dimension and reviewing Schramms left passage formula indicates that the paths cannot be described in terms of SLE.
We describe a percolation problem on lattices (graphs, networks), with edge weights drawn from disorder distributions that allow for weights (or distances) of either sign, i.e. including negative weights. We are interested whether there are spanning paths or loops of total negative weight. This kind of percolation problem is fundamentally different from conventional percolation problems, e.g. it does not exhibit transitivity, hence no simple definition of clusters, and several spanning paths/loops might coexist in the percolation regime at the same time. Furthermore, to study this percolation problem numerically, one has to perform a non-trivial transformation of the original graph and apply sophisticated matching algorithms. Using this approach, we study the corresponding percolation transitions on large square, hexagonal and cubic lattices for two types of disorder distributions and determine the critical exponents. The results show that negative-weight percolation is in a different universality class compared to conventional bond/site percolation. On the other hand, negative-weight percolation seems to be related to the ferromagnet/spin-glass transition of random-bond Ising systems, at least in two dimensions.
We investigate the geometric properties of loops on two-dimensional lattice graphs, where edge weights are drawn from a distribution that allows for positive and negative weights. We are interested in the appearance of spanning loops of total negative weight. The resulting percolation problem is fundamentally different from conventional percolation, as we have seen in a previous study of this model for the undiluted case. Here, we investigate how the percolation transition is affected by additional dilution. We consider two types of dilution: either a certain fraction of edges exhibit zero weight, or a fraction of edges is even absent. We study these systems numerically using exact combinatorial optimization techniques based on suitable transformations of the graphs and applying matching algorithms. We perform a finite-size scaling analysis to obtain the phase diagram and determine the critical properties of the phase boundary. We find that the first type of dilution does not change the universality class compared to the undiluted case whereas the second type of dilution leads to a change of the universality class.
We consider the negative weight percolation (NWP) problem on hypercubic lattice graphs with fully periodic boundary conditions in all relevant dimensions from d=2 to the upper critical dimension d=6. The problem exhibits edge weights drawn from disorder distributions that allow for weights of either sign. We are interested in in the full ensemble of loops with negative weight, i.e. non-trivial (system spanning) loops as well as topologically trivial (small) loops. The NWP phenomenon refers to the disorder driven proliferation of system spanning loops of total negative weight. While previous studies where focused on the latter loops, we here put under scrutiny the ensemble of small loops. Our aim is to characterize -using this extensive and exhaustive numerical study- the loop length distribution of the small loops right at and below the critical point of the hypercubic setups by means of two independent critical exponents. These can further be related to the results of previous finite-size scaling analyses carried out for the system spanning loops. For the numerical simulations we employed a mapping of the NWP model to a combinatorial optimization problem that can be solved exactly by using sophisticated matching algorithms. This allowed us to study here numerically exact very large systems with high statistics.
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}$.
We study the very long-range bond-percolation problem on a linear chain with both sites and bonds dilution. Very long range means that the probability $p_{ij}$ for a connection between two occupied sites $i,j$ at a distance $r_{ij}$ decays as a power law, i.e. $p_{ij} = rho/[r_{ij}^alpha N^{1-alpha}]$ when $ 0 le alpha < 1$, and $p_{ij} = rho/[r_{ij} ln(N)]$ when $alpha = 1$. Site dilution means that the occupancy probability of a site is $0 < p_s le 1$. The behavior of this model results from the competition between long-range connectivity, which enhances the percolation, and site dilution, which weakens percolation. The case $alpha=0$ with $p_s =1 $ is well-known, being the exactly solvable mean-field model. The percolation order parameter $P_infty$ is investigated numerically for different values of $alpha$, $p_s$ and $rho$. We show that in the ranges $ 0 le alpha le 1$ and $0 < p_s le 1$ the percolation order parameter $P_infty$ depends only on the average connectivity $gamma$ of sites, which can be explicitly computed in terms of the three parameters $alpha$, $p_s$ and $rho$.