# Package(s) for n-body simulations?

**URL:** <https://discourse.julialang.org/t/package-s-for-n-body-simulations/50294>\
**Category:** General Usage\
**Created:** [November 17, 2020, 12:59pm UTC](https://discourse.julialang.org/t/package-s-for-n-body-simulations/50294 "2020-11-17T12:59:21Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![oschulz](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/oschulz/32/2998_2.png) [@oschulz](https://discourse.julialang.org/u/oschulz)\
**Post date:** [November 17, 2020, 12:59pm UTC](https://discourse.julialang.org/t/package-s-for-n-body-simulations/50294/1 "2020-11-17T12:59:21Z")

</div>

What do we currently have in the ecosystem for n-body simulations (with electric-field-like interaction, point-like particles, no collisions)?

---

<div class="post-metadata">

**Author:** ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)\
**Post date:** [November 17, 2020, 2:24pm UTC](https://discourse.julialang.org/t/package-s-for-n-body-simulations/50294/2 "2020-11-17T14:24:21Z")

</div>

[https://github.com/SciML/NBodySimulator.jl](https://github.com/SciML/NBodySimulator.jl)

---

<div class="post-metadata">

**Author:** ![oschulz](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/oschulz/32/2998_2.png) [@oschulz](https://discourse.julialang.org/u/oschulz)\
**Post date:** [November 17, 2020, 2:31pm UTC](https://discourse.julialang.org/t/package-s-for-n-body-simulations/50294/3 "2020-11-17T14:31:53Z")

</div>

Can NBodySimulator handle up to say 1000 particles (via octree or so)?

I actually have these particles in a semiconductor, so I don’t need a fancy solver for the motion (momenta aren’t conserved for relevant time scales, particle speed is basically proportional to force, no acceleration either). But I need something that scales with N(log N), using Barnes–Hut or so.

---

<div class="post-metadata">

**Author:** ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)\
**Post date:** [November 17, 2020, 2:35pm UTC](https://discourse.julialang.org/t/package-s-for-n-body-simulations/50294/4 "2020-11-17T14:35:04Z")

</div>

It doesn’t have a multipole method yet, that’s still in development.

---

<div class="post-metadata">

**Author:** ![oschulz](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/oschulz/32/2998_2.png) [@oschulz](https://discourse.julialang.org/u/oschulz)\
**Post date:** [November 19, 2020, 12:20pm UTC](https://discourse.julialang.org/t/package-s-for-n-body-simulations/50294/5 "2020-11-19T12:20:33Z")

</div>

Thanks, I’ll definitely keep an eye on it!