Skip to main content
1 of 4
Joseph O'Rourke
  • 150.9k
  • 36
  • 358
  • 958

Random N-body problem

Suppose there are $n$ unit-mass particles whose initial positions are uniformly distributed in a unit-radius disk. Each particle is assigned a randomly oriented initial velocity vector of length $1$ (green below). Then the particles act upon one another via inverse-square gravity.


          [![NBodyRand8][1]][1]
      Dots show initial positions inside unit disk. Green vectors: initial velocity; red vectors: final velocity.

Q. What is the probability that $k \ge 1$ of the $n$ particles remain within a disk of some radius $R \ge 1$ forever?

In the illustration above, $n=8$ and $R=3$.

I am wondering if the answer is: zero, independent of $k$ and $R$ and the gravitational constant? Perhaps the answer differs for points in $\mathbb{R}^2$ vs. points in $\mathbb{R}^3$ (or $\mathbb{R}^d$, $d > 3$)?

Joseph O'Rourke
  • 150.9k
  • 36
  • 358
  • 958