Orbit Fit and Astrometric record for 16GR391

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: 16GR391   
# Created Sun Oct 19 01:11:28 2025
# Orbit generated from Bernstein formalism
# Fitting     22 observations of     22
# Arc:   1.08y
# First observation: 2015/03/18
#  Last observation: 2016/04/15
Preliminary a, adot, b, bdot, g, gdot:
   0.000020   0.036379  -0.000000   0.007239   0.027967   0.000000
# Chi-squared of fit:     9.49 DOF:     38 RMS:  0.12
# Min/Max residuals:    -0.37    0.24
# Exact a, adot, b, bdot, g, gdot:
  1.985950E-05  3.620641E-02  3.873324E-06  7.210231E-03  2.771449E-02 -4.439604E-03
# Covariance matrix:
  4.9391E-13  4.8727E-12 -1.6191E-14  1.1318E-12  3.5062E-12  1.4576E-10
  4.8727E-12  1.3116E-10 -4.2604E-13  2.9826E-11  7.7253E-11  3.8436E-09
 -1.6191E-14 -4.2604E-13  2.5275E-13 -3.3636E-13 -2.4638E-13 -1.2560E-11
  1.1318E-12  2.9826E-11 -3.3636E-13  7.1328E-12  1.7230E-11  8.7948E-10
  3.5062E-12  7.7253E-11 -2.4638E-13  1.7230E-11  5.1172E-11  2.2192E-09
  1.4576E-10  3.8436E-09 -1.2560E-11  8.7948E-10  2.2192E-09  1.1334E-07
#      lat0       lon0       xBary       yBary       zBary        JD0
    1.503317 -178.252443    0.075743    0.026079   -0.988918  2457099.995558
# Heliocentric elements and errors
Epoch:              2457090.5000  =  2015/03/09
Mean Anomaly:          357.85562 +/-     0.151
Argument of Peri:       24.27042 +/-     1.196
Long of Asc Node:      174.31250 +/-     0.001
Inclination:            11.33987 +/-     0.001
Eccentricity:         0.70054404 +/-    0.0045
Semi-Major Axis:    121.62674414 +/-    1.5418
Time of Perihelion: 2460008.8696 +/-     197.9
Perihelion:          36.42185384 +/-    0.7158
Aphelion:           206.83163444 +/-    2.6783
Period (y)             1341.3800 +/-     25.51
# Ecliptic coordinates at JD0 (AU and AU/d)
Ecliptic X          -37.04722902 +/-    0.0093
Ecliptic Y           -1.05472302 +/-    0.0003
Ecliptic Z            0.94676824 +/-    0.0002
Ecliptic XDOT         0.00056256 +/-    0.0000
Ecliptic YDOT        -0.00356692 +/-    0.0000
Ecliptic ZDOT         0.00070062 +/-    0.0000
# Distances at JD0 (AU)
Heliocenter to KBO   37.07433059 +/-    0.0093
Geocenter to KBO     36.08221392 +/-    0.0093
# Hcoef:  8.66

The following table shows the complete astrometric record for 16GR391. 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 (16GR391) followed by the observatory code and reference code for the source of the astrometry.

2015 03  18.49478  12 08 48.28   +00 41 03.6   24.2r 16GR391   T09  C~94sd      
2015 03  18.50552  12 08 48.21   +00 41 03.9   24.1r 16GR391   T09  C~94sd      
2015 03  18.57964  12 08 47.86   +00 41 06.6   24.1r 16GR391   T09  C~94sd      
2015 03  29.26630  12 07 58.01   +00 47 18.2   23.8y 16GR391   T09  C~94sd      
2015 03  29.27727  12 07 57.95   +00 47 18.6   23.9y 16GR391   T09  C~94sd      
2015 03  29.36170  12 07 57.55   +00 47 21.4   24.0y 16GR391   T09  C~94sd      
2016 04  02.38656  12 15 54.47   +00 23 07.9   23.7i 16GR391   T09  C~8HFu      
2016 04  02.41462  12 15 54.34   +00 23 08.7   23.9i 16GR391   T09  C~8HFu      
2016 04  02.43920  12 15 54.22   +00 23 09.7   23.9i 16GR391   T09  C~8HFu      
2016 04  02.46633  12 15 54.10   +00 23 10.8   24.0i 16GR391   T09  C~8HFu      
2016 04  02.47448  12 15 54.06   +00 23 10.9   24.0i 16GR391   T09  C~8HFu      
2016 04  02.52059  12 15 53.84   +00 23 12.4   23.8i 16GR391   T09  C~8HFu      
2016 04  04.34312  12 15 45.47   +00 24 13.9   25.2g 16GR391   T09  C~8HFu      
2016 04  04.36588  12 15 45.37   +00 24 14.6   25.3g 16GR391   T09  C~8HFu      
2016 04  04.39376  12 15 45.22   +00 24 15.6   25.4g 16GR391   T09  C~8HFu      
2016 04  04.41084  12 15 45.14   +00 24 16.3   25.1g 16GR391   T09  C~8HFu      
2016 04  15.35150  12 14 56.56   +00 30 09.2   24.4y 16GR391   T09  C~8HFu      
2016 04  15.36307  12 14 56.50   +00 30 09.8   23.7y 16GR391   T09  C~8HFu      
2016 04  15.37454  12 14 56.43   +00 30 10.4   23.9y 16GR391   T09  C~8HFu      
2016 04  15.41093  12 14 56.29   +00 30 11.1   23.7y 16GR391   T09  C~8HFu      
2016 04  15.45396  12 14 56.10   +00 30 12.8   23.5y 16GR391   T09  C~8HFu      
2016 04  15.45934  12 14 56.09   +00 30 12.6   24.4y 16GR391   T09  C~8HFu      

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.16       0.00     0.05
     2   0.0000     -1.08    -0.07      -0.14    -0.13
     3   0.0002     -6.97    -0.10       0.25    -0.03
     4   0.0295   -840.77     0.08      43.93     0.09
     5   0.0295   -841.76    -0.05      43.94     0.06
     6   0.0298   -848.38    -0.03      44.13    -0.04
     7   1.0428   6295.04    -0.04    1556.31    -0.03
     8   1.0429   6292.93     0.05    1556.27    -0.15
     9   1.0430   6290.88    -0.07    1556.47    -0.02
    10   1.0430   6288.79    -0.03    1556.77     0.19
    11   1.0431   6288.20     0.02    1556.62     0.02
    12   1.0432   6284.57     0.00    1556.68    -0.06
    13   1.0482   6144.79     0.09    1563.14     0.00
    14   1.0482   6143.14     0.21    1563.18    -0.02
    15   1.0483   6140.67    -0.09    1563.20    -0.08
    16   1.0484   6139.29    -0.14    1563.37     0.04
    17   1.0783   5329.72     0.24    1597.14    -0.13
    18   1.0784   5328.65     0.01    1597.33     0.04
    19   1.0784   5327.45    -0.37    1597.46     0.15
    20   1.0785   5325.24     0.05    1597.27    -0.13
    21   1.0786   5321.95    -0.14    1597.69     0.20
    22   1.0786   5321.89     0.19    1597.45    -0.05

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.

16GR391    quality flag:3

Type:       SCATNEAR  SCATNEAR  SCATNEAR

axisobj       117.458   120.134   115.593
ecceobj         0.690     0.698     0.683
incobj         11.380    11.379    11.382
qmin           34.947    34.679    35.394
qmax          214.145   232.140   209.444
amean         119.202   127.268   117.622
amin          114.784   120.133   114.467
amax          125.300   134.215   123.044
emean           0.696     0.716     0.688
emin            0.679     0.696     0.677
emax            0.711     0.731     0.702
imean          10.052    10.227    10.999
imin            9.237     9.306    10.376
imax           10.948    11.126    12.048
excite_mean     0.717     0.738     0.714
fracstop        1.000     1.000     1.000
cjmean          3.067     3.059     3.070

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       145.9     143.9      82.3
kozaiamp        122.0     103.5      76.9