## Transfer matrix computation of generalized critical polynomials in percolation

## Abstract

Percolation thresholds have recently been studied by means of a graph polynomial PB(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, 82), kagome, and (3, 122) 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, 82) = 0.676 803 329 · · ·, p _{c}(kagome) =more »

- Authors:

- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Ecole Normale Superieure, Paris (France). LPTENS; Univ. Pierre et Marie Curie, Paris (France)

- Publication Date:

- Research Org.:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

- Sponsoring Org.:
- USDOE

- OSTI Identifier:
- 1399712

- Report Number(s):
- LLNL-JRNL-593972

Journal ID: ISSN 1751-8113

- Grant/Contract Number:
- AC52-07NA27344

- Resource Type:
- Accepted Manuscript

- Journal Name:
- Journal of Physics. A, Mathematical and Theoretical

- Additional Journal Information:
- Journal Volume: 45; Journal Issue: 49; Journal ID: ISSN 1751-8113

- Publisher:
- IOP Publishing

- Country of Publication:
- United States

- Language:
- English

- Subject:
- 97 MATHEMATICS AND COMPUTING

### Citation Formats

```
Scullard, Christian R., and Jacobsen, Jesper Lykke. Transfer matrix computation of generalized critical polynomials in percolation. United States: N. p., 2012.
Web. doi:10.1088/1751-8113/45/49/494004.
```

```
Scullard, Christian R., & Jacobsen, Jesper Lykke. Transfer matrix computation of generalized critical polynomials in percolation. United States. doi:10.1088/1751-8113/45/49/494004.
```

```
Scullard, Christian R., and Jacobsen, Jesper Lykke. Thu .
"Transfer matrix computation of generalized critical polynomials in percolation". United States. doi:10.1088/1751-8113/45/49/494004. https://www.osti.gov/servlets/purl/1399712.
```

```
@article{osti_1399712,
```

title = {Transfer matrix computation of generalized critical polynomials in percolation},

author = {Scullard, Christian R. and Jacobsen, Jesper Lykke},

abstractNote = {Percolation thresholds have recently been studied by means of a graph polynomial PB(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 PB(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 PB(p) was defined by a contraction-deletion identity, similar to that satisfied by the Tutte polynomial. Here, we give an alternative probabilistic definition of PB(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, 82), kagome, and (3, 122) 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 pc(4, 82) = 0.676 803 329 · · ·, pc(kagome) = 0.524 404 998 · · ·, pc(3, 122) = 0.740 420 798 · · ·, comparable to the best simulation results. We also show that the alternative definition of PB(p) can be applied to study site percolation problems.},

doi = {10.1088/1751-8113/45/49/494004},

journal = {Journal of Physics. A, Mathematical and Theoretical},

number = 49,

volume = 45,

place = {United States},

year = {2012},

month = {9}

}

*Citation information provided by*

Web of Science

Web of Science