Is it possible to reliably generate a particular approximation of an ideal knot via a simulated annealing approach? - MathOverflow most recent 30 from http://mathoverflow.net2013-05-24T00:51:48Zhttp://mathoverflow.net/feeds/question/76545http://www.creativecommons.org/licenses/by-nc/2.5/rdfhttp://mathoverflow.net/questions/76545/is-it-possible-to-reliably-generate-a-particular-approximation-of-an-ideal-knot-vIs it possible to reliably generate a particular approximation of an ideal knot via a simulated annealing approach?UltraBlue062011-09-27T19:17:55Z2012-07-30T19:38:49Z
<p>Say I take a cord, tie a loose knot in three-dimensional space, and pull tightly on the ends to generate an approximation of an ideal knot. If the cord has a fixed knot topology and a random initial configuration in space, are there any conditions under which I am guaranteed to generate a particular ideal knot approximation? In other words, if I perform something akin to simulated annealing to approximate the ideal knot, is it possible for me to always arrive at the same final state (whatever that may be) within some small error? </p>
<p>[Sept. 29th, 2011] - Are there any conditions/constraints under which I might be able to reliably achieve some local (or global) minima in the energy landscape of a knot? </p>
http://mathoverflow.net/questions/76545/is-it-possible-to-reliably-generate-a-particular-approximation-of-an-ideal-knot-v/76555#76555Answer by Joseph O'Rourke for Is it possible to reliably generate a particular approximation of an ideal knot via a simulated annealing approach?Joseph O'Rourke2011-09-27T20:31:21Z2011-10-17T16:40:57Z<p>I suspect what you are seeking is not known.
One of the experts in this area is Jason Cantarella, who coauthored (with R. Kusner and J. Sullivan) one of the definitive
papers on this topic, "On the minimum ropelength of knots and links," <em>Inventiones Mathematicae</em>,
Volume 150, Number 2, 257-286, 2002; <a href="http://www.springerlink.com/content/jjkq7d3bdw1qqkqd/" rel="nofollow">Springer link</a>).
Jason has developed software that he calls <em>RidgeRunner</em> which implements a rope-length minimization,
knot-tightening procedure. (I cannot speak to the details of his software, in particular,
I don't know if it is "akin to simulated annealing.") He maintains <a href="http://www.jasoncantarella.com/movs/" rel="nofollow">a wonderful web page</a>
with (Quicktime) movies of his minimizer approaching what you call the "ideal knot" for many (more
than 100!) knots and links. Here, for example, is the endpoint of his optimizer for $8_5$:
<br />
<img src="http://cs.smith.edu/~orourke/MathOverflow/Cantarella.8_5.jpg" alt="8_5"><br />
<b>Update</b>. See the just-released paper, "The Shapes of Tight Composite Knots," <a href="http://arxiv.org/abs/1110.3262" rel="nofollow">arXiv:1110.3262 (math.DG)</a>, by Jason Cantarella and Al LaPointe and Eric Rawdon, for a description
of the <em>RidgeRunner</em> software mentioned above. It "proceeds by constrained gradient descent," and "is designed to stop at local minima of the ropelength function." They reduce the probability of "false local minima" with several
strategies. Their dataset now contains almost 1000 knots and links.</p>
http://mathoverflow.net/questions/76545/is-it-possible-to-reliably-generate-a-particular-approximation-of-an-ideal-knot-v/103536#103536Answer by Jason Cantarella for Is it possible to reliably generate a particular approximation of an ideal knot via a simulated annealing approach?Jason Cantarella2012-07-30T19:38:49Z2012-07-30T19:38:49Z<p>The only conditions that I know of under which it's absolutely known that gradient flow will converge to a minimum energy state is when the initial configuration is planar. Zheng-Xu He proved this for his renormalized $1/r^2$ potential 'Mobius' energy in <em>The Euler-Lagrange equation and heat flow for the Möbius energy</em> (CPAM 53, 2000). We proved a similar theorem for 'repulsive' energies on planar polygons in <em>An energy-driven approach to linkage unfolding</em> (SOCG 2004). </p>
<p>For a nonplanar configuration, you're right to think that there's a general theorem that ``simulated annealing has a positive probability of converging to the global minimum" which applies to these problems, but it's ineffective in practice. The configuration space is in practice very high dimensional (3 x number of vertices) and reducing an energy functional such as ropelength generally requires a coordinated global motion of these vertices, so it's quite rare to generate such a move randomly. AFAIK, nobody knows how to estimate the probability of convergence under these circumstances, so there's no hard information on how long to run an annealer. </p>
<p>FWIW, <em>RidgeRunner</em> really isn't an annealer: it generates 'coordinated motions' for tightening deterministically using a linear algebra algorithm to deflect the gradient of length for a polygon into a ropelength-decreasing direction using the active distance and curvature constraints. The algorithm is described in <em>Knot Tightening By Constrained Gradient Descent</em> (Experimental Math, 2011). </p>
<p>You're welcome to play with <em>RidgeRunner</em> yourself if you'd like to try tightening: the software can be downloaded at <a href="http://www.jasoncantarella.com/webpage/index.php?title=Ridgerunner" rel="nofollow">link text</a>.</p>