I was trying to simulate the precession of Mercury based on the perturbed solution, and my questions about its implementation in python can be seen here: https://scicomp.stackexchange.com/questions/19778/precession-of-mercury-python-simulation.

And as I was searching on-line, I saw many papers were using post-Newtonian method to do it, but the formula involved in PN always have very complicated expression and I didn't really find any code implementation, why don't they just use the perturbed solution(the solution in the link above) coming from Lagrangian(the geodesic of time-like particle) instead of using those complicated expression? Can someone explain?

(EDIT: here is a blog that talks about the post-Newtonian method approach: https://astrokode.wordpress.com/2014/05/03/the-precession-of-mercurys-perihelion-simulation/)

(Update: Eventually I used Leapfrog method and adding some correction to the Newtonian gravity to do the simulation in Vpython, the code can be found in the first link. And I assume the reason why people use PN method is simply because in real life you don't get static and nice spherically symmetric objects, and just use the perturbed solution in schwarzschild metric would not be accurate in a long time simulation.)

  • 2
    $\begingroup$ Might Computational Science be better suited for the programming aspect of this question? $\endgroup$ – Kyle Kanos May 29 '15 at 1:02
  • $\begingroup$ @KyleKanos I didn't know that exists, just posted one there. Thanks for the info $\endgroup$ – Sam May 29 '15 at 1:15
  • 2
    $\begingroup$ Our comp-phys is somewhat limited in scope (not about implementations). Typically, SE does not allow for cross-posting; your options would be either deleting the question & re-asking it on Computational Science or flagging it for moderator attention. I personally think that the last question (about PN corrections) should be on-topic here and can be left alone (it really stands alone as a separate question from the first bit anyways), so maybe you'd want to remove the first bit & post that as a separate question on Computational Science? $\endgroup$ – Kyle Kanos May 29 '15 at 1:15
  • $\begingroup$ @KyleKanos hope it's better now. $\endgroup$ – Sam May 29 '15 at 1:25
  • 1
    $\begingroup$ Mercury's relativistic precession is 43 arc seconds per century. That is tiny. I truly doubt you'll see such a subtle effect using leapfrog. The blog you cited uses IAS15, a 15th order integrator specifically designed for gravitational dynamics. $\endgroup$ – David Hammen Jun 6 '15 at 11:51

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.