At this time, the orbital condition code is 4
Orbital Elements at Epoch 2457800.5 (2017-Feb-16.0) TDB
Reference: JPL 13 (heliocentric ecliptic J2000)
This asteroid was simulated with the Mercury6 orbit simulator together with 100 virtual clones generated with the package R.
These 100 clones were generated so that their orbital parameters are normally distributed around the nominal value of asteroid 2016 WF9 and their standard deviation is almost equal to the uncertainty shown above.
Virtual Asteroids: summary
- period simulated: past 1e8 days
- time step 0.1 days
- ejection distance = 100 au
- N-body algorithm: Conservative Bulirsch-Stoer
- 1 out of 101 virtual asteroid was discarded because it would have collided with Jupiter
- 44 out of 101 virtual asteroids came from the outskirt of the solar system (i.e., there was a time in the past when their distance was more than 100 au - ejection distance from the simulator point of view), so it makes sense to think they were comets.
This is the density distribution of arrival time in the solar system:
Probability of being a comet
I built a table with three columns:
- virtual asteroid id
- year: time of arrival into solar system, or end-time of the simulation (right censored data)
- event: in case of arrival from the outskirt of the solar system, this flag is TRUE, otherwise it is FALSE coherently with column year.
|Year(in the past)||n.risk||n.event||survival||std.err||lower 95% CI||upper 95% CI|
we started with 101 virtual asteroids, then at Year 30084 in the past we see one of them coming from the outskirt of the solar system, so we do not count it as an asteroid ...then, at Year 67281 in the past, this event occurs again ... and so on till we arrive at the last virtual asteroid being counted as a comet at Year 273009 in the past (last row not shown).
Based on this, we can draw a plot showing the probability of a virtual asteroid being a comet as a function of time:
SurvminerMercury Simulator - Mercury6
Alboukadel Kassambara and Marcin Kosinski (2016). survminer: Drawing Survival Curves using 'ggplot2'
J.E.Chambers (1999) "A Hybrid
Symplectic Integrator that Permits Close Encounters between
Massive Bodies''. Monthly Notices of the Royal Astronomical
Society, vol 304, pp793-799.