On Jan 5th, I downloaded the orbital parameters of asteroid 2016 WF9 and their uncertainty from the Horizons Web interface:

param value 1-sigma

a 2.8729806 0.00200340

e 0.6582656 0.00022291

i 14.9994893 0.00329270

w 342.4286190 0.00216310

om 125.4222393 0.00428050

M 353.6898676 0.00394660

In spite of the fact that the uncertainty is still very high, I tried to make a simulation generating 100 virtual asteroids to see how many of them in the past would seem to come from a distance greater than 100 AU (cometary origin).

First, I used the R programming environment to generate the virtual asteroids.

An obvious check is that the mean value of their orbital parameters must be almost equal to the nominal values of asteroid 2016 WF9 and their standard deviation should be also almost equal to the uncertainty shown above.

In fact, looking at the virtual asteroids, I got:

e mean 0.6583 sd 2e-04

i mean 14.9988 sd 0.0035

w mean 342.4285 sd 0.0023

om mean 125.4218 sd 0.0046

M mean 353.6903 sd 0.004

As a second step I used the Mercury6 simulator by John E. Chambers, to simulate the past 1e8 days, output interval 100 days (hybrid algorithm).

The program simulated the behaviour of 101 asteroids (100 virtual asteroid + the "real" asteroid):

The result was that about 60 virtual asteroids would seem to have a cometary orgin (i.e. there was a time in the past when their distance from the sun was grater than 100 AU).

The earliest date is about 37500 years ago, the oldest date is constrained by the simulation period (about 275000 years ago) and the median value is 155000 years ago.

(for curiosity: the nominal orbital values of 2016 WF9 itself tells us that this can be a comet arrived 76000 years ago).

This is the graph showing the arrival time of the 60 virtual comets:

In the next weeks, when the orbit of asteroid 2016 WF9 will be better defined, I would like to repeat the simulation to see if and how much the simulation results will be different.

Kind Regards,

Alessandro Odasso

Citations

Mercury Simulator - Mercury6

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.

## No comments:

## Post a Comment

Note: Only a member of this blog may post a comment.