Abstract
The stick number of a knot is the minimum number of segments needed to build a polygonal version of the knot. Despite its elementary definition and relevance to physical knots, the stick number is poorly understood: for most knots we only know bounds on the stick number. We adopt a Monte Carlo approach to finding better bounds, producing very large ensembles of random polygons in tight confinement to look for new examples of knots constructed from few segments. We generated a total of 220 billion random polygons, yielding either the exact stick number or an improved upper bound for more than 40% of the knots with 10 or fewer crossings for which the stick number was not previously known. We summarize the current state of the art in Appendix A, which gives the best known bounds on stick number for all knots up to 10 crossings.
Acknowledgments
We are grateful to Colin Adams, Ryan Blair, Jason Cantarella, Harrison Chapman, Tetsuo Deguchi, Kate Hake, Gyo Taek Jin, Ken Millett, Eric Rawdon, Erica Uehara, and the anonymous referees for helpful comments, suggestions, and conversations about random polygons and stick knots.
Declaration of interest
No potential conflict of interest was reported by the author(s).
Notes
1 For knots up to 10 crossings we use Alexander–Briggs notation consistent with the (Perko-corrected) Rolfsen table [53]; for knots with 11 or more crossings, we use the Dowker–Thistlethwaite naming convention.
2 At least in high dimensions, a single hit-and-run step is unlikely to move very far, so it is preferable to do several steps at once, though not too many since hit-and-run steps are much more expensive than sampling from the torus. Our experience is that 10 steps provides a good balance.
3 The best theoretical result along these lines is due to Hake [26], who proved that the probability that a random equilateral hexagon is knotted is bounded above by ; empirically, the fraction of nontrivial knots among random hexagons is close to
.
4 As described above, we use , so that hit-and-run steps and torus samples are equally likely, and γ = 10, so that we iterate hit-and-run 10 times whenever we choose to do a hit-and-run step.
5 “Unfortunately, the diagrams of and
given by Rolfsen do not correspond to their Conway notation or Alexander polynomial—they are interchanged” [27]. This was only corrected in the 2003 edition of Rolfsen’s book—after Rawdon and Scharein’s article—so it is not surprising that our pipeline identified Rawdon and Scharein’s
as a
and identified their
as a
. This is an issue that researchers in the broader computational knot theory community should be aware of: for example, the invariants computed by SnapPy for
match the invariants in the Knot Atlas [38], pyknotid, and KnotInfo [18] for
, and conversely the SnapPy invariants for
agree with the Knot Atlas, pyknotid, and KnotInfo invariants for
, so whether you identify a given knot as
or
may depend on which software you use. We have used the Knot Atlas naming convention in this article.