試しに、総当たりで引力求めて動かす、って感じのをecsとijobforpalallelとburst使って作ってみました。エディタ上で500位までは60fpsでてましたが、800個とか1000位で40,30fps切ります。
計算は定数化出来るところは予め計算てはいますが、個数^2の計算量なので、もっと簡略化しないとこれ以上は無理なのかなと感じています。