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

Christian R. Scullard

Christian R. Scullard 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.MPphysics.comp-phphysics.plasm-ph

Catalog footprint

What is connected

15works
6topics
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

15 published item(s)

preprint2021arXiv

Critical points of the random cluster model with Newman-Ziff sampling

We present a method for computing transition points of the random cluster model using a generalization of the Newman-Ziff algorithm, a celebrated technique in numerical percolation, to the random cluster model. The new method is straightforward to implement and works for real cluster weight $q>0$. Furthermore, results for an arbitrary number of values of $q$ can be found at once within a single simulation. Because the algorithm used to sweep through bond configurations is identical to that of Newman and Ziff, which was conceived for percolation, the method loses accuracy for large lattices when $q>1$. However, by sampling the critical polynomial, accurate estimates of critical points in two dimensions can be found using relatively small lattice sizes, which we demonstrate here by computing critical points for non-integer values of $q$ on the square lattice, to compare with the exact solution, and on the unsolved non-planar square matching lattice. The latter results would be much more difficult to obtain using other techniques.

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
preprint2019arXiv

Bond percolation thresholds on Archimedean lattices from critical polynomial roots

We present percolation thresholds calculated numerically with the eigenvalue formulation of the method of critical polynomials; developed in the last few years, it has already proven to be orders of magnitude more accurate than traditional techniques. Here we report the result of large parallel calculations to produce what we believe may become the reference values of bond percolation thresholds on the Archimedean lattices for years to come. For example, for the kagome lattice we find $p_{\rm c}=0.524\,404\,999\,167\,448\,20 (1)$, whereas the best estimate using standard techniques is $p_{\rm c}=0.524\,404\,99(2)$. We further provide strong evidence that there are two classes of lattices: one for which the first three scaling exponents characterizing the finite-size corrections to $p_{\rm c}$ are $Δ=6,7,8$, and another for which $Δ=4,6,8$. We discuss the open questions related to the method, such as the full scaling law, as well as its potential for determining critical points of other models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical solution of the quantum Lenard-Balescu equation for a one-component plasma

We present a numerical solution of the quantum Lenard-Balescu equation using a spectral method, namely an expansion in Laguerre polynomials. This method exactly conserves both particles and energy and facilitates the integration over the dielectric function. To demonstrate the method, we solve the equilibration problem for a spatially homogeneous one-component plasma with various initial conditions. Unlike the more usual Landau/Fokker-Planck system, this method requires no input Coulomb logarithm; the logarithmic terms in the collision integral arise naturally from the equation along with the non-logarithmic order-unity terms. The spectral method can also be used to solve the Landau equation and a quantum version of the Landau equation in which the integration over the wavenumber requires only a lower cutoff. We solve these problems as well and compare them with the full Lenard-Balescu solution in the weak-coupling limit. Finally, we discuss the possible generalization of this method to include spatial inhomogeneity and velocity anisotropy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the growth constant for square-lattice self-avoiding walks

The growth constant for two-dimensional self-avoiding walks on the honeycomb lattice was conjectured by Nienhuis in 1982, and since that time the corresponding results for the square and triangular lattices have been sought. For the square lattice, a possible conjecture was advanced by one of us (AJG) more than 20 years ago, based on the six significant digit estimate available at the time. This estimate has improved by a further six digits over the intervening decades, and the conjectured value continued to agree with the increasingly precise estimates. We discuss the three most successful methods for estimating the growth constant, including the most recently developed Topological Transfer-Matrix method, due to another of us (JLJ). We show this to be the most computationally efficient of the three methods, and by parallelising the algorithm we have estimated the growth constant significantly more precisely, incidentally ruling out the conjecture, which fails in the 12th digit. Our new estimate of the growth constant is $$μ(\mathrm{square}) = 2.63815853032790\, (3).$$

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Potts-model critical manifolds revisited

We compute the critical polymials for the q-state Potts model on all Archimedean lattices, using a parallel implementation of the algorithm of (Jacobsen, J. Phys. A: Math. Theor. 47 135001) that gives us access to larger sizes than previously possible. The exact polynomials are computed for bases of size $6 \times 6$ unit cells, and the root in the temperature variable $v=e^K-1$ is determined numerically at $q=1$ for bases of size $8 \times 8$. This leads to improved results for bond percolation thresholds, and for the Potts-model critical manifolds in the real $(q,v)$ plane. In the two most favourable cases, we find now the kagome-lattice threshold to eleven digits and that of the $(3,12^2)$ lattice to thirteen. Our critical manifolds reveal many interesting features in the antiferromagnetic region of the Potts model, and determine accurately the extent of the Berker-Kadanoff phase for the lattices studied.

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
preprint2013arXiv

The percolation critical polynomial as a graph invariant

Every lattice for which the bond percolation critical probability can be found exactly possesses a critical polynomial, with the root in [0,1] providing the threshold. Recent work has demonstrated that this polynomial may be generalized through a definition that can be applied on any periodic lattice. The polynomial depends on the lattice and on its decomposition into identical finite subgraphs, but once these are specified, the polynomial is essentially unique. On lattices for which the exact percolation threshold is unknown, the polynomials provide approximations for the critical probability with the estimates appearing to converge to the exact answer with increasing subgraph size. In this paper, I show how this generalized critical polynomial can be viewed as a graph invariant, similar to the Tutte polynomial. In particular, the critical polynomial is computed on a finite graph and may be found using the recursive deletion-contraction algorithm. This allows calculation on a computer, and I present such results for the kagome lattice using subgraphs of up to 36 bonds. For one of these, I find the prediction p_c=0.52440572..., which differs from the numerical value, p_c=0.52440503(5), by only 6.9 x 10^{-7}.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Critical manifold of the kagome-lattice Potts model

Any two-dimensional infinite regular lattice G can be produced by tiling the plane with a finite subgraph B of G; we call B a basis of G. We introduce a two-parameter graph polynomial P_B(q,v) that depends on B and its embedding in G. The algebraic curve P_B(q,v) = 0 is shown to provide an approximation to the critical manifold of the q-state Potts model, with coupling v = exp(K)-1, defined on G. This curve predicts the phase diagram both in the ferromagnetic (v>0) and antiferromagnetic (v<0) regions. For larger bases B the approximations become increasingly accurate, and we conjecture that P_B(q,v) = 0 provides the exact critical manifold in the limit of infinite B. Furthermore, for some lattices G, or for the Ising model (q=2) on any G, P_B(q,v) factorises for any choice of B: the zero set of the recurrent factor then provides the exact critical manifold. In this sense, the computation of P_B(q,v) can be used to detect exact solvability of the Potts model on G. We illustrate the method for the square lattice, where the Potts model has been exactly solved, and the kagome lattice, where it has not. For the square lattice we correctly reproduce the known phase diagram, including the antiferromagnetic transition and the singularities in the Berker-Kadanoff phase. For the kagome lattice, taking the smallest basis with six edges we recover a well-known (but now refuted) conjecture of F.Y. Wu. Larger bases provide successive improvements on this formula, giving a natural extension of Wu's approach. The polynomial predictions are in excellent agreement with numerical computations. For v>0 the accuracy of the predicted critical coupling v_c is of the order 10^{-4} or 10^{-5} for the 6-edge basis, and improves to 10^{-6} or 10^{-7} for the largest basis studied (with 36 edges).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Potts critical frontiers of inhomogeneous and asymmetric bow-tie lattices

We study the critical frontiers of the Potts model on two-dimensional bow-tie lattices with fully inhomogeneous coupling constants. Generally, for the Potts critical frontier to be found exactly, the underlying lattice must be a 3-uniform hypergraph. A more general class of lattices are the 4-uniform ones, with unit cells contained within four boundary vertices. We demonstrate that in some cases, such lattices can be decomposed into triangular cells, and solved using a modification of standard techniques. This leads to the exact inhomogeneous Potts critical frontiers on various lattices, such as the bow-tie lattice with five different couplings, and critical points for asymmetric bow-tie lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The computation of generalized percolation critical polynomials by the deletion-contraction algorithm

Although every exactly known bond percolation critical threshold is the root in $[0,1]$ of a lattice-dependent polynomial, it has recently been shown that the notion of a critical polynomial can be extended to any periodic lattice. The polynomial is computed on a finite subgraph, called the base, of an infinite lattice. For any problem with exactly known solution, the prediction of the bond threshold is always correct, regardless of the base chosen. For unsolved problems, the polynomial is referred to as the generalized critical polynomial and provides an approximation that becomes more accurate with increasing number of bonds in the base, appearing to approach the exact answer. The polynomials are computed using the deletion-contraction algorithm, which quickly becomes intractable by hand for more than about 18 bonds. Here, I present generalized critical polynomials calculated with a computer program for bases of up to 36 bonds for all the Archimedean lattices, except the kagome which was considered in an earlier work. The polynomial estimates are generally within $10^{-5}$ to $10^{-7}$ of the numerical values, but the prediction for the $(4,8^2)$ lattice, though not exact, is not ruled out by simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transfer matrix computation of critical polynomials for two-dimensional Potts models

In our previous work we have shown that critical manifolds of the q-state Potts model can be studied by means of a graph polynomial P_B(q,v), henceforth referred to as the critical polynomial. This polynomial may be defined on any periodic two-dimensional lattice. It depends on a finite subgraph B, called the basis, and the manner in which B is tiled to construct the lattice. The real roots v = e^K - 1 of P_B(q,v) either give the exact critical points for the lattice, or provide approximations that, in principle, can be made arbitrarily accurate by increasing the size of B in an appropriate way. In earlier work, P_B(q,v) was defined by a contraction-deletion identity, similar to that satisfied by the Tutte polynomial. Here, we give a probabilistic definition of P_B(q,v), which facilitates its computation, using the transfer matrix, on much larger B than was previously possible. We present results for the critical polynomial on the (4,8^2), kagome, and (3,12^2) lattices for bases of up to respectively 96, 162, and 243 edges, compared to the limit of 36 edges with contraction-deletion. We discuss in detail the role of the symmetries and the embedding of B. The critical temperatures v_c obtained for ferromagnetic (v>0) Potts models are at least as precise as the best available results from Monte Carlo simulations or series expansions. For instance, with q=3 we obtain v_c(4,8^2) = 3.742 489(4), v_c(kagome) = 1.876 459 7(2), and v_c(3,12^2) = 5.033 078 49(4), the precision being comparable or superior to the best simulation results. More generally, we trace the critical manifolds in the real (q,v) plane and discuss the intricate structure of the phase diagram in the antiferromagnetic (v<0) region.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transfer matrix computation of generalised critical polynomials in percolation

Percolation thresholds have recently been studied by means of a graph polynomial $P_B(p)$, henceforth referred to as the critical polynomial, that may be defined on any periodic lattice. The polynomial depends on a finite subgraph $B$, called the basis, and the way in which the basis is tiled to form the lattice. The unique root of $P_B(p)$ in $[0,1]$ either gives the exact percolation threshold for the lattice, or provides an approximation that becomes more accurate with appropriately increasing size of $B$. Initially $P_B(p)$ was defined by a contraction-deletion identity, similar to that satisfied by the Tutte polynomial. Here, we give an alternative probabilistic definition of $P_B(p)$, which allows for much more efficient computations, by using the transfer matrix, than was previously possible with contraction-deletion. We present bond percolation polynomials for the $(4,8^2)$, kagome, and $(3,12^2)$ lattices for bases of up to respectively 96, 162, and 243 edges, much larger than the previous limit of 36 edges using contraction-deletion. We discuss in detail the role of the symmetries and the embedding of $B$. For the largest bases, we obtain the thresholds $p_c(4,8^2) = 0.676 803 329 ...$, $p_c(\mathrm{kagome}) = 0.524 404 998 ...$, $p_c(3,12^2) = 0.740 420 798 ...$, comparable to the best simulation results. We also show that the alternative definition of $P_B(p)$ can be applied to study site percolation problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Polynomial sequences for bond percolation critical thresholds

In this paper, I compute the inhomogeneous (multi-probability) bond critical surfaces for the (4,6,12) and (3^4,6) lattices using the linearity approximation described in (Scullard and Ziff, J. Stat. Mech. P03021), implemented as a branching process of lattices. I find the estimates for the bond percolation thresholds, p_c(4,6,12)=0.69377849... and p_c(3^4,6)=0.43437077..., compared with Parviainen's numerical results of p_c \approx 0.69373383 and p_c \approx 0.43430621 . These deviations are of the order 10^{-5}, as is standard for this method, although they are outside Parviainen's typical standard error of 10^{-7}. Deriving thresholds in this way for a given lattice leads to a polynomial with integer coefficients, the root in [0,1] of which gives the estimate for the bond threshold. I show how the method can be refined, leading to a sequence of higher order polynomials making predictions that likely converge to the exact answer. Finally, I discuss how this fact hints that for certain graphs, such as the kagome lattice, the exact bond threshold may not be the root of any polynomial with integer coefficients.

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