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: 23PY4
# Created Tue Sep 30 01:12:14 2025
# Orbit generated from Bernstein formalism
# Fitting 12 observations of 12
# Arc: 2.10y
# First observation: 2023/08/13
# Last observation: 2025/09/17
Preliminary a, adot, b, bdot, g, gdot:
0.000391 0.077961 0.000095 -0.010997 0.046596 0.000000
# Chi-squared of fit: 6.16 DOF: 18 RMS: 0.13
# Min/Max residuals: -0.32 0.28
# Exact a, adot, b, bdot, g, gdot:
2.590507E-05 7.368670E-02 -3.332404E-06 -1.033090E-02 4.429144E-02 -2.521059E-02
# Covariance matrix:
6.5887E-13 1.9947E-13 5.4779E-15 -4.5858E-14 9.1290E-13 1.9177E-12
1.9947E-13 2.9956E-12 6.1916E-14 -4.5605E-13 1.6647E-12 1.9577E-11
5.4779E-15 6.1916E-14 2.4907E-13 -1.4486E-13 3.2655E-14 4.2592E-13
-4.5858E-14 -4.5605E-13 -1.4486E-13 2.1064E-13 -2.5116E-13 -3.1316E-12
9.1290E-13 1.6647E-12 3.2655E-14 -2.5116E-13 2.3558E-12 1.0684E-11
1.9177E-12 1.9577E-11 4.2592E-13 -3.1316E-12 1.0684E-11 1.3449E-10
# lat0 lon0 xBary yBary zBary JD0
4.854959 -14.864175 0.431442 0.076836 -0.907656 2460170.005101
# Heliocentric elements and errors
Epoch: 2460170.5000 = 2023/08/14
Mean Anomaly: 358.84477 +/- 0.004
Argument of Peri: 191.57519 +/- 0.013
Long of Asc Node: 194.16911 +/- 0.001
Inclination: 9.24660 +/- 0.000
Eccentricity: 0.89367272 +/- 0.0002
Semi-Major Axis: 194.22695730 +/- 0.4296
Time of Perihelion: 2463343.2075 +/- 0.7
Perihelion: 20.65162476 +/- 0.0652
Aphelion: 367.80228984 +/- 0.8149
Period (y) 2706.9016 +/- 8.98
# Ecliptic coordinates at JD0 (AU and AU/d)
Ecliptic X 22.52236502 +/- 0.0008
Ecliptic Y -6.41933354 +/- 0.0002
Ecliptic Z 1.91080426 +/- 0.0001
Ecliptic XDOT -0.00028352 +/- 0.0000
Ecliptic YDOT 0.00479511 +/- 0.0000
Ecliptic ZDOT -0.00076819 +/- 0.0000
# Distances at JD0 (AU)
Heliocenter to KBO 23.49714753 +/- 0.0007
Geocenter to KBO 22.57772706 +/- 0.0008
# Hcoef: 8.58
The following table shows the complete astrometric record for 23PY4. 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 (23PY4) followed by the observatory code and reference code for the source of the astrometry.
2023 08 13.50430 22 57 45.439 -01 22 36.03 22.5r 23PY4 T09 C~8nlW 2023 08 13.57965 22 57 44.971 -01 22 39.18 22.4r 23PY4 T09 C~8nlW 2023 09 14.19207 22 54 12.282 -01 46 57.57 22.5r 23PY4 269 C~8nlW 2023 09 15.16258 22 54 05.674 -01 47 46.88 22.4r 23PY4 269 C~8nlW 2023 10 06.08037 22 51 53.953 -02 05 11.27 22.3r 23PY4 269 C~8nlW 2023 10 07.09587 22 51 48.348 -02 05 59.71 22.4r 23PY4 269 C~8nlW 2024 10 03.23446 23 08 57.621 -00 57 48.90 23PY4 G37 C~8nlW 2024 10 04.19796 23 08 51.708 -00 58 37.95 23PY4 G37 C~8nlW 2025 09 17.392612 23 28 32.425 +00 26 05.85 21.8w 23PY4 F51 C~98UL 2025 09 17.404689 23 28 32.345 +00 26 05.09 21.9w 23PY4 F51 C~98UL 2025 09 17.416806 23 28 32.255 +00 26 04.49 21.8w 23PY4 F51 C~98UL 2025 09 17.428913 23 28 32.178 +00 26 03.85 21.8w 23PY4 F51 C~98UL
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.08 0.00 -0.07
2 0.0002 -7.69 0.06 -0.21 -0.32
3 0.0868 -3512.72 -0.25 -120.71 0.28
4 0.0894 -3623.16 -0.10 -128.17 0.10
5 0.1467 -5848.70 0.03 -334.25 0.15
6 0.1495 -5944.95 0.19 -346.74 -0.06
7 1.1409 9883.97 0.02 -2511.28 -0.04
8 1.1436 9783.00 -0.03 -2522.33 -0.06
9 2.0969 28255.65 -0.09 -4688.78 0.10
10 2.0969 28254.23 -0.02 -4689.02 -0.05
11 2.0970 28252.72 -0.02 -4689.05 0.01
12 2.0970 28251.38 0.14 -4689.20 -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.
23PY4 quality flag:3 Type: CENTAURR CENTAURR CENTAURR axisobj 184.687 185.126 184.249 ecceobj 0.889 0.889 0.888 incobj 9.256 9.256 9.256 qmin 19.483 19.448 19.258 qmax 404.951 350.866 544.858 amean 178.370 156.811 221.022 amin 154.984 138.143 175.745 amax 212.872 185.800 282.924 emean 0.883 0.869 0.906 emin 0.864 0.850 0.880 emax 0.902 0.889 0.926 imean 9.474 9.132 9.150 imin 8.728 8.677 8.556 imax 10.754 9.583 9.847 excite_mean 0.899 0.883 0.920 fracstop 0.253 0.161 0.322 cjmean 2.415 2.422 2.399 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 226.0 218.8 222.0 kozaiamp 31.5 27.2 28.0