Asteroid Ephemerides / New XEPH Files - 2023 Release

Juan Conejero

PixInsight Staff
Staff member
A new set of 13 asteroid ephemerides files in XEPH format is now available on our distribution server. The new files cover the time span from 2019-12-01 to 2024-02-02 and include the first 620,000 numbered asteroids. These XEPH files can be used to compute high-precision ephemerides with our solar system ephemerides engine, as well as to generate asteroid annotations for plate-solved images with our standard AnnotateImage script.

The new XEPH files provide Chebyshev polynomial expansions for position in barycentric rectangular coordinates coherent with fundamental JPL planetary and lunar DE440 ephemerides (ICRS), with maximum truncation errors below 1e–10 au (about 15 meters) for the vast majority of objects in the main asteroid belt, and proportionally smaller errors for objects closer to the Sun.

These XEPH files have been generated by numerical integration of perturbed orbital motion with our standard EphemerisGenerator process. Initial conditions have been obtained from JPL Small–Body Element databases, latest ELEMENTS.NUMBR file for epoch 2023-02-25. The implemented physical model includes:
  • Perturbations from Mercury, Venus, Earth, Moon, Mars, Jupiter, Saturn, Uranus, Neptune and Pluto. State vectors and constants are provided by XEPH fundamental ephemerides (currently JPL DE440).

  • Perturbations from the 343 most massive asteroids used in DE440 numerical integrations. State vectors and constants are provided by XEPH asteroid ephemerides (currently JPL DE441 asteroid ephemerides).

  • Perturbations from the 30 most massive Kuiper belt objects (KBOs). State vectors and constants are provided by XEPH asteroid ephemerides (currently JPL DE441 asteroid ephemerides).

  • Relativistic perturbation terms (DE440's parametrized post-Newtonian (PPN) formulation) up to O(1/c^2) for the point–mass accelerations induced by all planets, Moon, Pluto and the five most massive asteroids.
The applied numerical integration process uses the Bulirsch–Stoer method for solution of ordinary differential equations. Our current implementation of the Bulirsch–Stoer algorithm is based on Boost's ODEINT C++ library.

Contrarily to previous editions, acceleration terms due to oblateness (zonal harmonics J2) have not been included in this release, since their contribution is negligible to the accuracy of computed positional ephemerides. This has allowed us to save a significant amount of computation time. Also in contrast to previous releases, where we used initial conditions obtained from the Asteroid Orbital Elements Database by Lowell Observatory, on this occasion we have preferred JPL orbital elements since they lead to more consistent results in all of our tests.

Please note that these files don't include the 343 most massive asteroids and the 30 largest Kuiper belt objects, which are already available directly on the PixInsight platform for a much larger time span (from 1550 to 2650) from DE440 ephemerides. So if you are going to generate asteroid annotations, remember to include an Asteroids annotation layer besides custom XEPH layers.

The entire numerical integration executed to generate these 13 XEPH files (1.3 GiB) has required about 60 hours on a dedicated workstation with an AMD Ryzen 9 5950X 16–core processor and 128 GiB of RAM running Kubuntu 20.04 Linux.

We hope you enjoy this new release. It is dedicated to our users who are interested in astronomy (hopefully there are still some). Merry Christmas to all!
 
Back
Top