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: 16GS391
# Created Sun Oct 19 01:11:28 2025
# Orbit generated from Bernstein formalism
# Fitting 23 observations of 23
# Arc: 1.84y
# First observation: 2015/03/25
# Last observation: 2017/01/26
Preliminary a, adot, b, bdot, g, gdot:
-0.000006 0.023688 -0.000001 -0.000449 0.023893 0.000000
# Chi-squared of fit: 10.77 DOF: 40 RMS: 0.12
# Min/Max residuals: -0.28 0.32
# Exact a, adot, b, bdot, g, gdot:
1.546263E-05 2.378102E-02 -4.646521E-07 -4.512448E-04 2.431634E-02 -6.792019E-04
# Covariance matrix:
2.6660E-13 7.3550E-14 5.1641E-16 -4.4331E-15 7.4637E-13 8.6708E-12
7.3550E-14 7.0132E-12 2.2524E-14 -1.3049E-13 1.3511E-11 2.0605E-10
5.1641E-16 2.2524E-14 2.0478E-13 -1.5591E-13 4.3194E-14 6.7064E-13
-4.4331E-15 -1.3049E-13 -1.5591E-13 1.7437E-13 -2.5806E-13 -3.9324E-12
7.4637E-13 1.3511E-11 4.3194E-14 -2.5806E-13 2.8277E-11 4.1655E-10
8.6708E-12 2.0605E-10 6.7064E-13 -3.9324E-12 4.1655E-10 6.2641E-09
# lat0 lon0 xBary yBary zBary JD0
1.824651 -178.184705 -0.043430 0.031743 -0.992678 2457106.972828
# Heliocentric elements and errors
Epoch: 2457100.5000 = 2015/03/19
Mean Anomaly: 316.33569 +/- 2.826
Argument of Peri: 168.06354 +/- 3.156
Long of Asc Node: 60.44797 +/- 0.023
Inclination: 2.08751 +/- 0.001
Eccentricity: 0.03683472 +/- 0.0028
Semi-Major Axis: 43.24407796 +/- 0.0373
Time of Perihelion: 2469698.8170 +/- 815.1
Perihelion: 41.65119465 +/- 0.1260
Aphelion: 44.83696126 +/- 0.1268
Period (y) 284.3794 +/- 0.37
# Ecliptic coordinates at JD0 (AU and AU/d)
Ecliptic X -42.07756122 +/- 0.0090
Ecliptic Y -1.37717853 +/- 0.0003
Ecliptic Z 1.30944676 +/- 0.0003
Ecliptic XDOT 0.00015624 +/- 0.0000
Ecliptic YDOT -0.00267968 +/- 0.0000
Ecliptic ZDOT -0.00005313 +/- 0.0000
# Distances at JD0 (AU)
Heliocenter to KBO 42.12045145 +/- 0.0090
Geocenter to KBO 41.12460190 +/- 0.0090
# Hcoef: 8.03
The following table shows the complete astrometric record for 16GS391. 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 (16GS391) followed by the observatory code and reference code for the source of the astrometry.
2015 03 25.47205 12 09 33.86 +00 57 08.1 24.9g 16GS391 T09 C~94sd 2015 03 25.48280 12 09 33.82 +00 57 08.4 25.6g 16GS391 T09 C~94sd 2015 03 25.49994 12 09 33.73 +00 57 08.9 25.1g 16GS391 T09 C~94sd 2015 03 29.28550 12 09 16.83 +00 58 56.9 22.8y 16GS391 T09 C~94sd 2015 03 29.29365 12 09 16.81 +00 58 57.3 23.0y 16GS391 T09 C~94sd 2015 03 29.36170 12 09 16.48 +00 58 59.3 23.2y 16GS391 T09 C~94sd 2016 04 02.39692 12 13 53.73 +00 27 11.8 23.8i 16GS391 T09 C~8HFu 2016 04 02.49633 12 13 53.28 +00 27 14.5 23.7i 16GS391 T09 C~8HFu 2016 04 02.49902 12 13 53.29 +00 27 14.5 23.7i 16GS391 T09 C~8HFu 2016 04 02.52059 12 13 53.17 +00 27 15.4 23.7i 16GS391 T09 C~8HFu 2016 04 04.34312 12 13 45.19 +00 28 05.9 25.0g 16GS391 T09 C~8HFu 2016 04 04.35168 12 13 45.14 +00 28 06.1 25.3g 16GS391 T09 C~8HFu 2016 04 04.36588 12 13 45.08 +00 28 06.4 25.5g 16GS391 T09 C~8HFu 2016 04 04.39376 12 13 44.96 +00 28 07.3 25.2g 16GS391 T09 C~8HFu 2016 04 04.41084 12 13 44.89 +00 28 07.8 25.4g 16GS391 T09 C~8HFu 2016 04 15.35150 12 12 58.37 +00 32 57.4 23.7y 16GS391 T09 C~8HFu 2016 04 15.37454 12 12 58.28 +00 32 57.7 23.6y 16GS391 T09 C~8HFu 2016 04 15.40815 12 12 58.17 +00 32 58.7 23.5y 16GS391 T09 C~8HFu 2016 04 15.45127 12 12 57.96 +00 32 59.7 23.6y 16GS391 T09 C~8HFu 2017 01 26.62885 12 22 55.02 -00 34 02.5 24.5z 16GS391 T09 C~94sd 2017 01 26.63421 12 22 55.00 -00 34 02.3 24.1z 16GS391 T09 C~94sd 2017 01 26.64229 12 22 54.99 -00 34 02.3 23.9z 16GS391 T09 C~94sd 2017 01 26.64766 12 22 54.98 -00 34 02.2 24.6z 16GS391 T09 C~94sd
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.00 0.00 0.04
2 0.0000 -0.67 0.13 0.04 0.08
3 0.0001 -2.11 -0.04 -0.04 0.02
4 0.0104 -277.62 -0.03 -1.73 -0.16
5 0.0105 -278.05 0.14 -1.48 0.09
6 0.0106 -283.39 -0.20 -1.61 0.02
7 1.0238 4291.25 -0.11 -97.77 -0.03
8 1.0240 4283.98 -0.12 -97.98 -0.12
9 1.0240 4284.12 0.21 -97.92 -0.06
10 1.0241 4282.11 -0.23 -97.81 0.08
11 1.0291 4152.15 0.10 -99.10 0.13
12 1.0291 4151.38 -0.05 -99.22 0.03
13 1.0291 4150.44 0.03 -99.30 -0.04
14 1.0292 4148.43 0.04 -99.19 0.11
15 1.0293 4147.26 0.12 -99.15 0.18
16 1.0592 3391.62 -0.22 -111.04 -0.06
17 1.0593 3390.26 -0.03 -111.30 -0.28
18 1.0594 3388.35 0.32 -111.04 0.04
19 1.0595 3385.06 -0.06 -111.38 -0.22
20 1.8430 13219.52 0.08 -238.51 0.02
21 1.8430 13219.16 -0.10 -238.45 0.08
22 1.8430 13219.03 0.02 -238.51 0.01
23 1.8431 13218.85 0.01 -238.48 0.03
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.
16GS391 quality flag:3 Type: CLASSICAL CLASSICAL CLASSICAL axisobj 42.852 42.946 42.764 ecceobj 0.031 0.039 0.023 incobj 2.089 2.089 2.089 qmin 41.010 40.776 41.194 qmax 45.252 45.680 44.918 amean 43.052 43.148 42.961 amin 42.721 42.817 42.634 amax 43.397 43.484 43.293 emean 0.024 0.031 0.019 emin 0.005 0.013 0.000 emax 0.043 0.051 0.038 imean 1.006 0.996 1.023 imin 0.127 0.084 0.178 imax 1.618 1.615 1.619 excite_mean 0.030 0.036 0.027 fracstop 1.000 1.000 1.000 cjmean 3.089 3.090 3.088 libcent 0 -180.0 -180.0 -180.0 libamp 0 -180.0 -180.0 -180.0 libcent 1 -180.0 -180.0 -180.0 libamp 1 -180.0 -180.0 -180.0 libcent 2 -180.0 -180.0 -180.0 libamp 2 -180.0 -180.0 -180.0 libcent 3 -180.0 -180.0 -180.0 libamp 3 -180.0 -180.0 -180.0 libcent 4 -180.0 -180.0 -180.0 libamp 4 -180.0 -180.0 -180.0 kozaimean 179.7 184.5 169.5 kozaiamp 179.8 179.8 179.9