I recently read about the Breakthrough Initiative to launch "StarShot", a nano-probe that is designed to travel to Alpha Centauri at $0.2c$. One of the challenges to be solved involves the precise determination of the orbital ephemerides of any exoplanets in the system.
My question is related to how we might go about simulating not only the orbital ephemerides of the planets but also the probe itself, given that it is travelling at relativistic speeds. In essence, what do we need to consider to simulate the entire physical system?
Presumably, the usual method for speeds $v \ll c$ is to consider a Newtonian n-body approach and simply integrate the equations of motion with the gravitational inverse square law, using a numerical scheme such as Runge-Kutta or, something more sophisticated, like a symplectic integrator?
To me, this seems appropriate only for the planets themselves and not the probe. Given that speeds near $c$ are occurring it should be necessary to include at least a special relativistic approach. However, after more consideration, given that gravitation is involved (even without strong-field regions such as black holes) it would seem to be necessary to involve general relativity.
If the latter is the case, does that mean we would need to bring in the heavy machinery from numerical relativity, such as a foliating curved space-time and solving the Einstein field equations numerically, with a matter distribution of planets as the initial condition?
To me this seems to be "using a sledgehammer to crack a nut" when all we really need is the trajectory of the probe from the Alpha Centauri reference frame and the positions of the planets from the probe reference frame.
In addition, does anybody know of any available software (open source or otherwise) that might help with this simulation?