4 ms·

I'm the first author on this paper being discussed here. <https://arxiv.org/abs/2303.05567 https://arxiv.org/abs/2303.05567> This post has driven the most traf

by zyrxvo 4y ago

I'm the first author on this paper being discussed here. <https://arxiv.org/abs/2303.05567 https://arxiv.org/abs/2303.05567> This post has driven the most traffic to my site ever. Forgive me for not understanding the social rules here, but I made an account today to join the incredible discussion here. If you have any questions I will do my best to answer them. However, you've all been answering each other's questions as well as I ever could've already. Edit: Added link to the preprint.

qmarchi 4y ago
Welcome to HN! Guidelines and the FAQ exist at the bottom of the site. You noted that you ran over a thousand simulations for this. Can you explain a bit more into what you used to create the variations and any challenges that you encountered while running the simulations?
- zyrxvo 4y ago
  Thank you! It's exciting to be here. The Guidelines seemed brief, but informative. Thank you as well for pointing out the FAQ. We used REBOUND <https://rebound.readthedocs.io/ https://rebound.readthedocs.io/> to run our N-body simulations and a modified version of REBOUNDx <https://reboundx.readthedocs.io/ https://reboundx.readthedocs.io/> to add the first-order post-Newtonian correction from GR and adjust it over time <https://github.com/zyrxvo/reboundx/tree/GR_Sweep https://github.com/zyrxvo/reboundx/tree/GR_Sweep>. We used the exact same setup and initial conditions for each of the 1280 simulations. Thus, the only variation between them was our choice of how slowly to turn off GR. "Turning off GR" is similar too, but not the same as increasing the speed of light from ~3e8 m/s to infinity. We started with GR at the currently accepted value (equivalent to c=3e8) and then ramped it down linearly so that c=∞ at some time t=τ (again, not exactly the same as changing the speed of light, but related to it). We would expect that without any differences, we could run two N-body simulations with REBOUND and reproduce them exactly, bit for bit even though the system is chaotic. Therefore, with even the slightest changes, such as slowly altering the perihelion precession rate of Mercury, chaos would cause any two similar or nearby trajectories to diverge to different results. For the most part running the simulations themselves didn't pose any challenges. The most difficult aspect of running them was the walltime, or the real-world time it takes for a simulation to run. With a timestep of ~3 days, it took about 82 days to run one simulation to 12.5 billion years. Based on the ending times of all 1280 simulations, it required about 132 years of walltime to run. This is why we were very grateful for Compute Canada and the Niagara supercomputer, because we were able to run all of these simulations in less than 3-4 months.

by zyrxvo 4y ago

qmarchi 4y ago

Welcome to HN! Guidelines and the FAQ exist at the bottom of the site. You noted that you ran over a thousand simulations for this. Can you explain a bit more into what you used to create the variations and any challenges that you encountered while running the simulations?

zyrxvo 4y ago
Thank you! It's exciting to be here. The Guidelines seemed brief, but informative. Thank you as well for pointing out the FAQ. We used REBOUND <https://rebound.readthedocs.io/ https://rebound.readthedocs.io/> to run our N-body simulations and a modified version of REBOUNDx <https://reboundx.readthedocs.io/ https://reboundx.readthedocs.io/> to add the first-order post-Newtonian correction from GR and adjust it over time <https://github.com/zyrxvo/reboundx/tree/GR_Sweep https://github.com/zyrxvo/reboundx/tree/GR_Sweep>. We used the exact same setup and initial conditions for each of the 1280 simulations. Thus, the only variation between them was our choice of how slowly to turn off GR. "Turning off GR" is similar too, but not the same as increasing the speed of light from ~3e8 m/s to infinity. We started with GR at the currently accepted value (equivalent to c=3e8) and then ramped it down linearly so that c=∞ at some time t=τ (again, not exactly the same as changing the speed of light, but related to it). We would expect that without any differences, we could run two N-body simulations with REBOUND and reproduce them exactly, bit for bit even though the system is chaotic. Therefore, with even the slightest changes, such as slowly altering the perihelion precession rate of Mercury, chaos would cause any two similar or nearby trajectories to diverge to different results. For the most part running the simulations themselves didn't pose any challenges. The most difficult aspect of running them was the walltime, or the real-world time it takes for a simulation to run. With a timestep of ~3 days, it took about 82 days to run one simulation to 12.5 billion years. Based on the ending times of all 1280 simulations, it required about 132 years of walltime to run. This is why we were very grateful for Compute Canada and the Niagara supercomputer, because we were able to run all of these simulations in less than 3-4 months.