The following information shows the result of the orbit fit based on Gary Bernstein's method. Most of the information should be self-explanatory. Take special note that while the original Bernstein software works with barycentric coordinates, we convert these results into a heliocentric coordinate system.
# Object: 16WY85 # Created Wed Jul 17 01:11:26 2024 # Orbit generated from Bernstein formalism # Fitting 14 observations of 14 # Arc: 65.02d # First observation: 2016/11/28 # Last observation: 2017/02/01 Preliminary a, adot, b, bdot, g, gdot: 0.000001 0.031270 0.000001 0.002307 0.028029 0.000000 # WARNING Fitting with energy constraint # Chi-squared of fit: 3.25 DOF: 23 RMS: 0.09 # Min/Max residuals: -0.20 0.15 # Exact a, adot, b, bdot, g, gdot: 1.826568E-05 3.050723E-02 9.384020E-07 2.241411E-03 2.741710E-02 -4.459714E-03 # Covariance matrix: 1.9127E-13 -6.2481E-11 -1.8584E-15 -1.6165E-11 -7.5522E-12 1.7932E-09 -6.2481E-11 1.7683E-07 7.3312E-11 3.9434E-08 4.8758E-08 2.3100E-06 -1.8584E-15 7.3312E-11 2.0521E-13 1.2812E-11 3.0714E-11 3.8466E-09 -1.6165E-11 3.9434E-08 1.2812E-11 9.0852E-09 9.8020E-09 2.2022E-07 -7.5522E-12 4.8758E-08 3.0714E-11 9.8020E-09 1.7629E-08 1.7882E-06 1.7932E-09 2.3100E-06 3.8466E-09 2.2022E-07 1.7882E-06 3.4696E-04 # lat0 lon0 xBary yBary zBary JD0 -13.780563 144.365283 0.969584 -0.048978 -0.200531 2457721.022749 # Heliocentric elements and errors Epoch: 2457720.5000 = 2016/11/28 Mean Anomaly: 327.31747 +/- 45.965 Argument of Peri: 338.48524 +/- 102.724 Long of Asc Node: 217.02454 +/- 0.643 Inclination: 14.21456 +/- 0.043 Eccentricity: 0.25042197 +/- 0.5400 Semi-Major Axis: 45.01993789 +/- 11.8471 Time of Perihelion: 2467737.0664 +/- 13521.1 Perihelion: 33.74595643 +/- 25.8802 Aphelion: 56.29391936 +/- 28.4671 Period (y) 302.0755 +/- 119.24 # Ecliptic coordinates at JD0 (AU and AU/d) Ecliptic X -28.39727576 +/- 0.1394 Ecliptic Y 21.54244270 +/- 0.0999 Ecliptic Z -8.68816181 +/- 0.0421 Ecliptic XDOT -0.00146434 +/- 0.0015 Ecliptic YDOT -0.00270365 +/- 0.0011 Ecliptic ZDOT 0.00032342 +/- 0.0004 # Distances at JD0 (AU) Heliocenter to KBO 36.68741288 +/- 0.1232 Geocenter to KBO 36.47358472 +/- 0.1766 # Hcoef: 7.86
The following table shows the complete astrometric record for 16WY85. The first three columns show the date of observation. The next six columns are RA and DEC. The next column (when provided) is the observed magnitude and filter. The next column is the object name (16WY85) followed by the observatory code and reference code for the source of the astrometry.
2016 11 28.52196 09 28 30.40 +00 22 27.4 23.4r 16WY85 T09 C~85Ik 2016 11 28.55492 09 28 30.38 +00 22 26.8 23.7r 16WY85 T09 C~85Ik 2016 11 28.56332 09 28 30.37 +00 22 26.3 23.9r 16WY85 T09 C~85Ik 2016 11 28.61822 09 28 30.33 +00 22 25.4 23.9r 16WY85 T09 C~85Ik 2016 11 28.63097 09 28 30.32 +00 22 25.1 24.0r 16WY85 T09 C~85Ik 2016 11 29.50349 09 28 29.63 +00 22 08.8 23.1i 16WY85 T09 C~85Ik 2016 11 29.62468 09 28 29.52 +00 22 06.6 23.6z 16WY85 T09 C~85Ik 2016 11 29.62736 09 28 29.51 +00 22 06.6 23.5z 16WY85 T09 C~85Ik 2016 11 29.63275 09 28 29.51 +00 22 06.6 23.5z 16WY85 T09 C~85Ik 2016 11 29.64147 09 28 29.49 +00 22 06.3 23.0z 16WY85 T09 C~85Ik 2017 02 01.49516 09 24 43.81 +00 24 54.7 23.0i 16WY85 T09 C~85Ik 2017 02 01.50344 09 24 43.77 +00 24 55.1 23.1i 16WY85 T09 C~85Ik 2017 02 01.53643 09 24 43.59 +00 24 55.7 22.9i 16WY85 T09 C~85Ik 2017 02 01.54461 09 24 43.54 +00 24 55.9 23.1i 16WY85 T09 C~85Ik
The following table shows the residuals to the orbit fit. The first coumn is the point number. The second column is the time, in years, measured from the first observation. The third and fifth columns are the regularized positions used in the orbit fit. The fourth and sixth columns are the residuals, in arc seconds, for RA and Dec respectively.
1 0.0000 0.00 -0.14 0.00 0.07
2 0.0001 -0.09 -0.05 -0.66 0.13
3 0.0001 -0.07 0.02 -1.19 -0.20
4 0.0003 -0.35 0.06 -2.23 -0.03
5 0.0003 -0.39 0.09 -2.56 -0.08
6 0.0027 -4.92 0.05 -21.34 -0.03
7 0.0030 -5.77 0.09 -23.95 0.02
8 0.0030 -5.91 -0.03 -24.00 0.03
9 0.0030 -5.91 0.01 -24.00 0.15
10 0.0031 -6.10 -0.11 -24.38 -0.04
11 0.1779 -3264.22 0.05 -959.34 -0.08
12 0.1779 -3264.91 0.03 -959.15 0.14
13 0.1780 -3267.66 0.01 -959.46 -0.02
14 0.1780 -3268.44 -0.09 -959.51 -0.04
The following table comes from a 10My integration of the orbit of the object. Three columns are shown. The first column is the result of integrating the nominal orbit. The other two columns are based on clones of the nominal orbit that are +/- 3 sigma from the nominal orbit. If all three types agree then the classificiation is deemed secure. The basis for these calculations is described in more detail in AJ, 129, 1117 (2005). Any use made of these calculations should refer to and credit this publication and the Deep Ecliptic Survey Team.