System Data:
Simulations:
Simulation data:
System Data:
Simulations:
Simulation data:
System Data:
Simulations:
Simulation data:
System Data:
Simulations:
Simulation data:
System Data:
Simulations:
Simulation data:
System Data:
Simulations:
Simulation data:
System Data:
Simulations:
Simulation data:
And sma vs time: zoom
System Data:
Simulations:
Simulation data:
And comparing 4 different
runs with the same initial parameters...
System Data:
Simulations:
Simulation data:
And sma vs time: zoom and zoom2 and zoom3!
HD169830
HD169830
Mstar = 1.4 Msun
M1 = 2.88 Mjup M2 = 4.04 Mjup
e1 = 0.31 e2 = 0.33
a1 = 0.81 AU a2 = 3.6 AU
HZin = 1.372 AU
HZout = 3.038 AU
TP data:
Results:
NTP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [1.372,3.038] AU
Ran for 106 years with dt = 0.031 yrs = 11.241 days
After 5x105 years there were 0 TPs left. System very unstable!
And comparing 4 different
runs with the same initial parameters...
position inc vs time ecc vs time sma vs time HD160691
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
TP data:
Results:
NTP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [0.816,1.808] AU
Ran for 106 years with dt = 0.035 yrs = 12.965 days
After just 5x104 years there were 0 TP left. Extremely unstable - even the planets get kicked!
And comparing 4 different runs with the same
initial parameters...
(final eccentricity of the second planet is 20!)
position inc vs time ecc vs time sma vs time Eps Eridani
Eps Eridani
Mstar = 0.80 Msun
M1 = 0.86 Mjup M2 = 0.10 Mjup
e1 = 0.608 e2 = 0.30
a1 = 3.30 AU a2 = 40 AU
HZin = 0.448 AU
HZout = 0.992 AU
TP data:
Results:
NTP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [0.448,0.992] AU
Ran for 106 years with dt = 0.017 yrs = 6.123 days
After 106 years there were 911 TPs left. System pretty stable.
position inc vs time ecc vs time sma vs time HD37124-new
HD37124 - new data
Mstar = 0.91 Msun
M1 = 0.75 Mjup M2 = 1.20 Mjup
e1 = 0.10 e2 = 0.69
a1 = 0.54 AU a2 = 2.50 AU
HZin = 0.580 AU
HZout = 1.284 AU
TP data:
Results:
NTP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [0.580,1.284] AU
Ran for 106 years with dt = 0.021 yrs = 7.597 days
After only 104 years there were 0 TPs left. Very unstable...
And comparing 3 different runs with the same
initial parameters...
position inc vs time ecc vs time sma vs time HD37124-old
HD37124 - old data
Mstar = 0.91 Msun
M1 = 0.86 Mjup M2 = 1.01 Mjup
e1 = 0.10 e2 = 0.40
a1 = 0.54 AU a2 = 2.95 AU
HZin = 0.580 AU
HZout = 1.284 AU
TP data:
Results:
NTP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [0.580,1.284] AU
Ran for 106 years with dt = 0.021 yrs = 7.597 days
After 106 years there were 369 TPs left. The most stable system yet!
And comparing 3 different runs with the same
initial parameters...
position inc vs time ecc vs time sma vs time HD74156-new
HD74156 - new data
Mstar = 1.05 Msun
M1 = 1.86 Mjup M2 = 6.17 Mjup
e1 = 0.636 e2 = 0.583
a1 = 0.294 AU a2 = 3.4 AU
HZin = 0.772 AU
HZout = 1.709 AU
TP data:
Results:
NTP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [0.772,1.709] AU
Ran for 106 years with dt = 0.008 yrs = 2.839 days
After 106 years there was only 1 TP left. Unstable system.
And comparing 3 different runs with the same
initial parameters...
position inc vs time ecc vs time sma vs time HD74156-old
HD74156 - old data
Mstar = 1.05 Msun
M1 = 1.56 Mjup M2 = 7.5 Mjup
e1 = 0.649 e2 = 0.395
a1 = 0.28 AU a2 = 3.47 AU
HZin = 0.772 AU
HZout = 1.709 AU
TP data:
Results:
NTP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [0.772,1.709] AU
Ran for 106 years with dt = 0.007 yrs = 2.639 days
After 106 years there was only 1 TP left. Unstable system.
position inc vs time ecc vs time sma vs time HD38529-new
HD38529 - new data
Mstar = 1.39 Msun
M1 = 0.78 Mjup M2 = 12.7 Mjup
e1 = 0.29 e2 = 0.36
a1 = 0.129 AU a2 = 3.68 AU
HZin = 1.352 AU
HZout = 2.995 AU
TP data:
Results:
NTP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [1.352,2.995] AU
Ran for 106 years with dt = 0.002 yrs = 0.717 days
After 106 years there were only 11 TP left. Possible a thin stable zone??
And sma vs time: zoom1 and zoom2
position inc vs time ecc vs time sma vs time HD38529-old
HD38529 - old data
Mstar = 1.39 Msun
M1 = 0.77 Mjup M2 = 11.2 Mjup
e1 = 0.27 e2 = 0.34
a1 = 0.129 AU a2 = 3.51 AU
HZin = 1.352 AU
HZout = 2.995 AU
TP data:
Results:
NTP = 1000, ecc = [0.0,0.1], inc = [0.0,0.2], sma = [1.352,2.995] AU
Ran for 106 years with dt = 0.002 yrs = 0.720 days
After 106 years there was only 5 TP left. Unstable system.
position inc vs time ecc vs time sma vs time Last Updated: 5 June 2004