## It's knot what you know

Why is it that some knots seem to hold tight while others readily slip apart? Patil *et al.* develop a theoretical analysis of the stability of knots and find links between topological parameters (twist charge, crossing numbers, handedness) and mechanical stability. The theory is confirmed using simulations and experiments on color-changing fibers that optically show localized stress differences in different parts of the knot as the two strands are pulled apart. The authors show why some common knots slip easily and untie, whereas others hold tight.

*Science*, this issue p. 71

## Abstract

Knots play a fundamental role in the dynamics of biological and physical systems, from DNA to turbulent plasmas, as well as in climbing, weaving, sailing, and surgery. Despite having been studied for centuries, the subtle interplay between topology and mechanics in elastic knots remains poorly understood. Here, we combined optomechanical experiments with theory and simulations to analyze knotted fibers that change their color under mechanical deformations. Exploiting an analogy with long-range ferromagnetic spin systems, we identified simple topological counting rules to predict the relative mechanical stability of knots and tangles, in agreement with simulations and experiments for commonly used climbing and sailing bends. Our results highlight the importance of twist and writhe in unknotting processes, providing guidance for the control of systems with complex entanglements.

Knots are among the oldest, most enduring human technologies, as valuable to ancient builders (*1*) and mariners (*2*) as to modern engineers and surgeons (*3*). Thought to predate the wheel (*1*), knotted structures owe their extraordinary longevity and widespread usage to an inherent mechanical robustness that arises from the subtle interplay of topology, elasticity, and friction. Over the course of many centuries, sailors, weavers, climbers, and surgeons have acquired a wealth of knowledge about the benefits and drawbacks of various types of knots (*1*, *2*). Experience has taught us that certain knots are more stable than others, but we are still largely unable to predict the mechanical behavior of knots and tangles from basic topological observables (*4*), such as the number and relative ordering of crossings. Although recent experimental and theoretical research has revealed important insights into the competition between force transmission and friction in special classes of knots (*5*, *6*), hitches (*7*), and fabrics (*8*), there currently exists no comprehensive mathematical theory (*9*) linking the topological and mechanical properties of knotted elastic structures.

Physical knots and their topology first assumed a central role in science with the introduction of Kelvin’s vortex-atom model in the 1860s (*10*). Since then, the fundamental importance of entangled structures has become firmly established in a diverse range of disciplines and contexts (*9*). In physics, for example, interactions between knotted defect lines are essential to understanding and controlling dynamics and mixing in classical and complex fluids (*11*–*14*), including liquid crystals (*15*), plasmas (*16*), and quantum fluids (*17*). Whereas the energetic costs associated with topological transformations are typically low in liquids and gases (*17*), they tend to become prohibitively large in entangled solids (*5*, *6*). This fact has profound consequences for the stability and function of natural and engineered structures, from the microscopic knots in DNA (*18*, *19*), proteins (*20*–*22*), and polymers (*23*) to knitted clothes (*8*) and macroscopic meshworks (*24*). Achieving a unified understanding of these various systems requires taking into account not only their topological but also their elastic properties. Because key concepts from topology and elasticity theory remain applicable over a wide range of scales, deciphering the topological principles (*5*, *6*, *9*) that determine the mechanical stability of knots promises insights into a broad spectrum of physically entangled structures. Therefore, our main goal is to identify generic topological counting rules that enable us to estimate which members of a given family of elastic knots are the most robust against untying. To this end, we combined elements from mathematical and physical knot theory (*9*, *25*) with optomechanical experiments and quantitative continuum modeling (Fig. 1).

We were interested in tying two lines together so that they form a stable longer rope, a task known as “tying a bend” among sailors (*2*). Mathematically, this configuration describes an oriented 2-tangle, defined as the union of two oriented open curves embedded in space (*4*). Although an elegant mathematical formalism exists to describe certain simple families of 2-tangles (*26*), little is known theoretically about even the most basic bend knots used in practice. We constructed a topological phase diagram that explains the relative stability of a selection of bends that are commonly used in the sailing and climbing communities. To validate the underlying topological model, we compared its predictions with simulations of an optomechanically verified continuum theory and with quantitative measurements using laboratory “prisoner’s escape” experiments.

Our optomechanical experiments use recently developed color-changing photonic fibers (*27*) that allow for the imaging of strain in knots (Fig. 1, A and B). These fibers derive structural color from a multilayer cladding composed of alternating layers of transparent elastomers with distinct refractive indices wrapped around an elastic core. Their coloration varies with the thickness of layers in the periodic cladding, which changes upon elongation or bending. As is typical of macroscopic materials at room temperature, the persistence lengths of the fibers used in our experiments are several orders of magnitude larger than the diameters of the tightened knots (*28*), with empirical knowledge (*2*) indicating that essential knot properties are only weakly dependent on the elastic modulus. Theoretically, we describe knotted fibers using a damped Kirchhoff model (*5*, *28*, *29*) validated through comparison with photographs depicting the strain-induced color changes in mechanoresponsive photonic fibers (Fig. 1, A to C). Simulating the tightening process of a 1-tangle, corresponding to a single knotted fiber pulled at both ends (Fig. 1, A and B, and movie S1), reveals the relative strengths and localization of the bending and stretching strains (Fig. 1, D and E), which are not individually discernible in our experiments. Furthermore, the Kirchhoff model highlights why topological considerations (*4*) alone do not suffice to explain the mechanical behaviors (*2*) of real-world knots: Loosening or tightening a knot transforms any of its planar projections according to a sequence of three elementary topology-preserving Reidemeister moves, *R*_{1}, *R*_{2}, and *R*_{3} (Fig. 1F). Despite being topologically equivalent, the move *R*_{1} is energetically distinct as it involves substantial changes in strain, whereas moves *R*_{2} and *R*_{3} are energetically favored soft modes (Fig. 1F), implying that physical knots preferentially deform by *R*_{2} and *R*_{3}. Thus, to link the physical properties of tangled fibers to their topology, one must merge concepts from classical mathematical knot theory (*4*) with elasticity theory (*5*, *6*, *9*, *30*).

Continuum simulations provide guidance for how one can complement bare topological knot diagrams (*4*) with coarse-grained mechanical information that is essential for explaining why certain knots are more stable than others (Fig. 2). In contrast to a 1-tangle, which is tightened by pulling diametrically at its two ends (Fig. 2A), each strand of a bend knot has one pulled and one free end (Fig. 2B). Therefore, the local fiber velocity directions in the center-of-mass frame of the bend knot define natural fiber orientations on the underlying 2-tangle (Fig. 2B), thus establishing mapping between bends and oriented 2-tangles. At each contact crossing, the fibers mutually generate a frictional self-torque with well-defined handedness, depending on the relative velocity and ordering of the two fiber strands (Fig. 2, B and C). In analogy to the coarse-graining procedure underlying Ising-type spin models, we can associate a unit twist charge *q _{i}* = ±1 with each vertex

*i*in the planar 2-tangle diagram, where the sign of

*q*reflects the combined handedness of the torques acting on the two intersecting strands (Fig. 2C). The sum

_{i}A key puzzle of physical knot theory (*9*), the empirically observed mechanical difference between the visually similar reef and granny knots (Fig. 3, A and B), may be understood as a consequence of this torque–writhe correspondence in 2-tangles. The underlying mechanism becomes evident by considering a pair of crossings as shown in Fig. 3C. Whereas equally directed torques lead to rolling, opposite torques promote locking and thus stabilize a knot against untying. The overall stability of 2-tangles therefore depends on the self-torque distribution along the fibers, as encoded by the vertex twist charges *q _{i}* = ±1 in an untightened knot diagram (Fig. 3, D to F). The above argument suggests the following topological twist fluctuation energy per site:

*N*is the crossing number,

_{0}= 1 – 1/

*N*can be interpreted as a ground-state energy density (

*28*). Equation 1 has the form of a ferromagnetic energy for an Ising-type spin model with long-range interactions, emphasizing the concept of knots as strongly coupled systems.

In addition to twist locking for large values of τ, knots can be stabilized when their internal structure forces fiber segments to slide tangentially against each other. For example, the reef knot and the thief knot both have τ = 1, but because their pulled ends differ, friction makes the reef knot more stable (Fig. 3, D and G). At the coarse-grained level of planar knot diagrams, these friction effects correspond to edge-to-edge interactions dominated by pairs of edges sharing a face and pulled in opposite directions (Fig. 3, D to F). To formalize this notion, each edge around a face *F* is assigned a weight of +1 or –1 if it winds around *F* in the anticlockwise or clockwise direction, respectively. Each face then contributes a friction energy given by the net circulation of the edges around the face, *C _{F}*, normalized by the total number of edges

*e*. This yields the total circulation energy:

_{F}The topological parameters *N*, τ, and Γ allow us to rationalize the stability of a large class of popular knots used by sailors and climbers (Fig. 3G). These variables are easy to evaluate from knot diagrams (Fig. 3, D to F) and reflect topology-induced forces and torques throughout the knot. As such, the triplet (*N*, τ, Γ) captures both essential topological and mechanical structure hidden within knots. The (*N*, τ, Γ) phase diagram explains existing empirical knowledge for simple knots (*2*), as well as predictions of the Kirchhoff model about the relative strength and stability of more complex 2-tangles (Fig. 3G). We verified these predictions independently in experiments by mimicking the prisoner’s escape problem (Fig. 4A) with two thin Dyneema fibers made from ultra–high molecular weight polyethylene tied together (*28*). Of the two pulled ends for each knot, one is fixed in the experimental apparatus and the other is perturbed while suspending incrementally higher masses until the knot pulls through. Although the Kirchhoff model cannot account for surface contact details (*30*), the experimental data for the critical loads agree quantitatively with the simulations for simple knots and, more importantly, confirm the predicted qualitative stability differences between various commonly used knots (Fig. 4B). Notably, both theory and experiments indicate that the Zeppelin knot is more secure against untying than the popular alpine butterfly knot (Fig. 4, B and C).

To conclude, the above analysis shows how basic topological counting rules can be used to estimate the relative stability of frequently encountered knots and tangles. From a theoretical perspective, the parallels with long-range coupled spin systems suggest that the statistical mechanics (*4*, *21*) of general knotted structures can be formulated within this framework. Tangled vortices (*12*, *17*) in complex fluids and defect loops in liquid crystals (*15*) may permit similar statistical descriptions through reduction to topological crossing diagrams. In elastic systems, joint experimental and theoretical progress is needed to untangle longstanding puzzles regarding the statistics of knots in DNA (*18*) and proteins (*20*, *21*), where thermal effects induce a finite persistence length, and other macroscopic structures (*8*, *24*). In sailing, climbing, and many other applications, nontopological material parameters and contact geometry (*30*) also play important mechanical roles and must be included in more refined continuum models to quantitatively describe practically relevant knotting phenomena. From a broader conceptual and practical perspective, the above topological mechanics framework seems well suited for designing and exploring new classes of knots with desired behaviors under applied load.

## Supplementary Materials

science.sciencemag.org/content/367/6473/71/suppl/DC1

Materials and Methods

Supplementary Text

Figs. S1 to S6

Tables S1 and S2

Movies S1 and S2

This is an article distributed under the terms of the Science Journals Default License.

## References and Notes

**Acknowledgments:**J.D. thanks the Isaac Newton Institute for Mathematical Sciences for support and hospitality during the program “The Mathematical Design of New Materials” (supported by EPSRC grant EP/R014604/1) when work on this paper was undertaken. We thank J. Takagi for producing Fig. 4A.

**Funding:**This work was supported by an Alfred P. Sloan Research Fellowship (J.D.), a Complex Systems Scholar Award from the James S. McDonnell Foundation (J.D.), the Brigham and Women’s Hospital through a Stepping Strong Innovator Award (J.D.S. and M.K.), and the National Science Foundation through the “Designing Materials to Revolutionize and Engineer our Future” program (DMREF-1922321 to J.D.S. and M.K.).

**Author contributions:**V.P.P. and J.D. developed theory. V.P.P. performed simulations, for which J.D.S. and M.K. provided data and code for converting strain into perceived knot color. J.D.S. and M.K. designed color-changing fibers and conceived optomechanical experiments. J.D.S. conducted optomechanical experiments and provided the description of the experiments in the supplementary materials. J.D.S. and V.P.P. performed in-lab prisoner’s escape experiments. V.P.P. and J.D. wrote the first draft of the paper. All authors discussed and revised the manuscript.

**Competing interests:**The authors declare no competing interests.

**Data and materials availability:**The code used for numerical simulations is available on Zenodo (

*31*). All data are available in the main text or the supplementary materials.