You can always add the binary companion, out-of-plane, to this system and see what happens. The difficulty you'll have is that you have to run the simulation very slowly since planet E is so close. So the binary companion will move very slowly. It will take a long time to simulate any significant amount of time.

You can always add the binary companion, out-of-plane, to this system and see what happens. The difficulty you'll have is that you have to run the simulation very slowly since planet E is so close. So the binary companion will move very slowly. It will take a long time to simulate any significant amount of time.

I really have some difficulty understanding the applied method in this paper . The role of the binary at >1000 Au is rather vague . Further the key feature of he paper is the proposed inclination of the primary star rotation axis and the plane of the planets which is estimated to be around 50°. I think it is not possible to account for the stars rotation axis in GravSim ( which can give tidal influence upon the plantes orbits ) . One method for compensating this in Gravsim could be treating the primary star is a very very close binary , which surely introduces tidal effects on the orbits as well as precession .