|
|
Set up n particles randomly distributed in a 1/r^2 profile (so that the circular velocity is constant). Truncate the halo at radius of 1Mpc. These particles are all fixed in place.
|
|
|
|
|
|
Put particle in a circular orbit at 0.5 Mpc and evolve for around 30 orbits.
|
|
|
|
|
|
![orbit_plot_1](/uploads/e80734543b568b3f41a845465182ff7c/orbit_plot_1.png)
|
|
|
|
|
|
Decrease timestep:
|
|
|
|
|
|
![orbit_plot_2](/uploads/9b7c0b14224f284766588339ec7f7b79/orbit_plot_2.png)
|
|
|
|
|
|
![orbit_plot_3](/uploads/31288711ec049dbdb665f864b78e6da9/orbit_plot_3.png)
|
|
|
|
|
|
Change softening:
|
|
|
|
|
|
![orbit_plot_1](/uploads/88e9891ed5f46291e58e3d48281ec147/orbit_plot_1.png)
|
|
|
|
|
|
![orbit_plot_3](/uploads/e6a697213ed301a5aab7e3418031fa3d/orbit_plot_3.png)
|
|
|
|
|
|
Increase order of multipoles:
|
|
|
|
|
|
![orbit_plot_2](/uploads/ef9667f92a5e227c8ba9bf92ee962cce/orbit_plot_2.png) |
|
|
\ No newline at end of file |