2- and 3-body problems when gravity is not inverse-square - MathOverflow most recent 30 from http://mathoverflow.net 2013-05-21T12:47:37Z http://mathoverflow.net/feeds/question/35980 http://www.creativecommons.org/licenses/by-nc/2.5/rdf http://mathoverflow.net/questions/35980/2-and-3-body-problems-when-gravity-is-not-inverse-square 2- and 3-body problems when gravity is not inverse-square Joseph O'Rourke 2010-08-18T16:16:37Z 2012-09-04T16:01:55Z <p>Suppose that gravity did not follow an inverse-square law, but was instead a central force diminishing as $1/d^p$ for distance separation $d$ and some power $p$. Two questions:</p> <ol> <li><p>Presumably the 2-body problem still factors into two independent 1-body problems, results in planar motion, and can be solved. Have the orbits (the equivalents of elliptical and parabolic orbits for $p=2$) been worked out for other (perhaps specific) values of $p$?</p></li> <li><p>In some sense, the 3-body problem for $p=2$ cannot be solved. Most systems are choatic; see <a href="http://faculty.ifmo.ru/butikov/Projects/Collection.html" rel="nofollow">this interesting collection</a> of Eugene Butikov. Only a few periodic solutions are known; see the nice <a href="http://www.ams.org/samplings/feature-column/fcarc-orbits1" rel="nofollow">article</a> by Bill Casselman on the discovery of "choreographies." Is the situation simpler for other values of $p$? Perhaps $p=1$?</p></li> </ol> <p>References and pointers would be appreciated. Thanks!</p> <p><b>Edit</b>. Thanks to Agol, Ken, and José. I've now looked at Arnolʹd's <em>Huygens</em> and Needham (but not yet Arnolʹd's <em>Classical Mechanics</em>). Indeed, as the commenters say, there is a remarkable 2-body result for $p=1$: the orbits for a linearly attractive force are ellipses, for a linearly repulsive force, hyperbolae. This depends on the Kasner-Arnolʹd theorem stating that for each power law, there is a <em>dual</em> power law that maps orbits of one to orbits of the other. Newton proved in <em>Principia</em> that elliptical orbits result if and only if the force is inverse-linear or inverse-square. The Kasner-Arnolʹd theorem explains why.</p> http://mathoverflow.net/questions/35980/2-and-3-body-problems-when-gravity-is-not-inverse-square/35998#35998 Answer by José Figueroa-O'Farrill for 2- and 3-body problems when gravity is not inverse-square José Figueroa-O'Farrill 2010-08-18T18:04:24Z 2010-08-18T18:04:24Z <p>Not a complete answer, but there is a theorem (proved, for instance, in Arnold's <em>Mathematical Methods in Classical Mechanics</em>) that says that the only radial potentials for which all bound orbits are closed are the Kepler potential $\propto 1/r$ and the harmonic oscillator $\propto r^2$.</p> http://mathoverflow.net/questions/35980/2-and-3-body-problems-when-gravity-is-not-inverse-square/36007#36007 Answer by Andreas Rüdinger for 2- and 3-body problems when gravity is not inverse-square Andreas Rüdinger 2010-08-18T18:49:31Z 2010-08-18T19:01:28Z <p>Regarding question 1):</p> <p>a) In addition to the potentials for which all bound orbits are closed there are further central force potentials with analytic solutions, please cf.<br> <a href="http://en.wikipedia.org/wiki/Classical_central-force_problem#Central_forces_with_exact_solutions" rel="nofollow">http://en.wikipedia.org/wiki/Classical_central-force_problem#Central_forces_with_exact_solutions</a> (Sorry, not quite what you were asking for [I read too fast], but perhaps interesting nevertheless.) </p> <p>b) Perhaps this paper on "Multiparticle systems with a particle-interaction potential homogeneous of degree α = −2" by Borisov et al. is interesting for you: <a href="http://ics.org.ru/doc?pdf=1359&amp;dir=e" rel="nofollow">http://ics.org.ru/doc?pdf=1359&amp;dir=e</a>. ("Some new integrable and superintegrable systems generalizing the classical ones are also described."). </p> http://mathoverflow.net/questions/35980/2-and-3-body-problems-when-gravity-is-not-inverse-square/37413#37413 Answer by Richard Montgomery for 2- and 3-body problems when gravity is not inverse-square Richard Montgomery 2010-09-01T19:42:58Z 2010-09-01T19:42:58Z <p>The answers to question (1) for the 2 body problem are fine, and complete enough.</p> <p>Regarding (2). The 3 body problem (and N-body) with p =3 is significantly simpler than with $p \ne 3$. The added simplicity is due to the occurenc of an additional integral which comes out of the Lagrange Jacobi identity for the evolution of the total moment of inertia $I$. This identity asserts that $d^2 I/ dt^2 = 4 H + (4 - 2(p-1)) U$ where $H = K - U$ is the total energy, with $K$ the kinetic energy and $U$ the NEGATIVE of the total potential energy, a function which is homogeneous of degree $p-1$. When $p =3$ we get $d^2 I/ dt^2 = 4H = const.$!. (The total moment of inertia is a the squared norm relative to the mass inner product' and as such measures the total size of the system. )</p> <p>For details on this Lagrange-Jacobi identity and its use see the first sections of my paper Hyperbolic Pants fit a three-body problem', Ergodic Theory and Dynamical Systems, Volume 25, - June 2005, 921-947, which you can also find on my web site <a href="http://count.ucsc.edu/~rmont/papers/list.html" rel="nofollow">http://count.ucsc.edu/~rmont/papers/list.html</a> or on the arXivs. Also see the references there.</p> <p>For a study of choreographies with various $p$- potentials see the paper by Fujiwara et al. Choreographic Three Bodies on the Lemniscate':J. Phys. A: Math. Gen. 36 (21 March 2003) 2791-2800, available on his web site ( or the ArXivs). <a href="http://www.clas.kitasato-u.ac.jp/~fujiwara/nBody/nbody.html" rel="nofollow">http://www.clas.kitasato-u.ac.jp/~fujiwara/nBody/nbody.html</a></p> <p>The discoverer of the figure eight, Cris Moore, in his beautiful 2 page paper Braids and Classical Gravity' (which Casselman should have a ref. to) found numerically, and argues convincingly that as one increases $p$ more and more braid types'' (and hence choreographies) appear. It is known that all possible braid types (and so choreography types) occur as soon as $p =2$. </p> <p>The cases $p \ge 2$ are often called `strong-force potentials'' and from the variational perspective are much simpler than $p &lt; 2$ for the reason that with the strong force potentials all collision paths have infinite action. This fact regarding action is surprising, since with the strong force condition in force it seems that almost all bounded solutions end in collision. This "seems" is a theorem for the 2-body problem, and for the negative energy three body problem when $p=2$. </p> http://mathoverflow.net/questions/35980/2-and-3-body-problems-when-gravity-is-not-inverse-square/106354#106354 Answer by jbc for 2- and 3-body problems when gravity is not inverse-square jbc 2012-09-04T15:21:50Z 2012-09-04T15:21:50Z <p>Not really an answer but too long for a comment.</p> <p>In addition to the cases already mentioned, there is a third one where the orbits can be described simply in analytical terms. This is where the force is proportional to $\frac 1{r^3}$---the curves are then Cotes spirals. These three cases were known to Newton.</p> <p>Maclaurin then discovered a two parameter family of orbits for each power law. They are obtained from the Maclaurin spirals (as they are now known), i.e., the curves $r^n=\cos (n \theta)$---by dilation and rotation. If $n$ is a whole number, then these orbits are closed---this is not in contradition to Bertrand's theorem, of course.</p> <p>The reason for this is that the functions $(\cos (d \theta))^{\frac 1 d}$ are such that $f+f''$ is proportional to a power of $f$.</p> <p>This fact lies behind further remarkable properties of these curves (for example, their curvatures are proportional to a power of the distance from the origin and they are catenaries for suitable power laws (the index of the power law depends in a simple fashion on $d$).</p> http://mathoverflow.net/questions/35980/2-and-3-body-problems-when-gravity-is-not-inverse-square/106361#106361 Answer by Bill Satzer for 2- and 3-body problems when gravity is not inverse-square Bill Satzer 2012-09-04T16:01:55Z 2012-09-04T16:01:55Z <p>A couple of Don Saari's papers from the 1970's era look at p = 2 + epsilon and show that there is no dramatic change in dynamics for a small perturbation of the exponent.</p>