System Data:

HD160691   Mstar = 1.08 Msun M1 = 1.70 Mjup    M2 = 1.00 Mjup e1 = 0.31               e2 = 0.80 a1 = 1.50 AU        a2 = 2.3 AU HZin = 0.816 AU HZout = 1.808 AU

Simulations:

TP data:
N_TP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [0.816,1.808] AU Ran for 10⁶ years with dt = 0.035 yrs = 12.965 days

After just 10⁴ years there were 0, 611, 40 and 922 TP left. But in all cases a planet was ejected! Extremely unstable system...

So let's look at this system is a bit more detail and compare 4 different runs to check for sensitivity to initial conditions. If this is correct, then one could assume that the observational parameters are actually not correct since they would not result in a stable system for long enough for us to observe it.

Final planet parameters for the four simulations was:

• V1:
  planet1 - a=0.07, e=20.3, i=180
  planet2 - a=25.2, e=0.91, i=0.0
  so P1 flung into star & P2 eccentricity pumped
• V2:
  planet1 - a=1.07, e=0.26, i=0.0
  planet2 - a=10.4, e=1.54, i=180
• V3:
  planet1 - a=0.81, e=0.86, i=0.0
  planet2 - a=0.77, e=129.6, i=0.0
• V4:
  planet1 - a=2.93, e=0.62, i=0.0
  planet2 - a=12.7, e=1.12, i=180.0

position - V1	position - V2	position - V3	position - V4
ecc vs time - V1	ecc vs time - V2	ecc vs time - V3	ecc vs time - V4
inc vs time - V1	inc vs time - V2	inc vs time - V3	inc vs time - V4
sma vs time - V1	sma vs time - V2	sma vs time - V3	sma vs time - V4