Genoid Algorithm for Orbit Determination

• Genoid algorithm is a meta-heuristic orbit-fitting tool that uses genetic search methods to derive orbital elements and system mass for binary asteroid systems.
• It applies a global search strategy over Keplerian orbital parameters using relative astrometric inputs from instruments like SPHERE and HST for sparse and degenerate datasets.
• Key outputs include orbital elements, system mass, and bulk density, which help infer physical properties such as ice-rich compositions of asteroids.

Searching arXiv for papers relevant to the "Genoid" orbital-fitting algorithm and its documented uses. Astrophysical usage of the term Genoid algorithm refers to an orbit-determination method used for binary or multiple asteroid systems. In the material documented for asteroid satellites, Genoid is described as a meta-heuristic orbit-determination algorithm and also as a genetic algorithm that searches for the dynamical parameters of the orbital motion in a binary or multiple system through successive generations of randomly matched solutions (Yang et al., 2020, Minker et al., 19 Mar 2025). Its role is to fit relative astrometry of a satellite with respect to a primary, recover a best-fit orbit in the EQJ2000 reference frame, and derive physically important quantities such as the system mass, the orbital pole, and, when combined with independent size or shape information, the bulk density of the primary (Yang et al., 2020, Minker et al., 19 Mar 2025).

1. Definition and scope

Within the documented asteroid-satellite literature, Genoid is an orbital-fitting tool for sparse astrometric datasets. In the study of the satellite S/2019 (31) 1 around (31) Euphrosyne, it is identified as the meta-heuristic orbit solver used to fit measured relative positions and derive the mass of the system from the satellite orbit (Yang et al., 2020). In the study of very wide binary asteroids, it is described more specifically as a genetic algorithm operating through successive generations of candidate solutions (Minker et al., 19 Mar 2025).

These descriptions indicate a global-search strategy rather than a purely local differential correction. This suggests that Genoid is intended for regimes in which orbital inversion is underconstrained, aliased, or strongly degenerate because of limited cadence or incomplete orbital coverage. That interpretation is consistent with both documented use cases: a newly detected close satellite observed at only a few epochs (Yang et al., 2020), and very wide binary asteroids with long periods and sparse archival observations (Minker et al., 19 Mar 2025).

The same corpus also makes clear what Genoid is not. It is not the categorical notion of a genoid used in lambda calculus and first-order logic (0712.3088), nor is it related to Genocchi type polynomials (Kurt et al., 2010). In the asteroid literature, the word denotes an orbit-fitting algorithm.

2. Observational inputs and orbital model

The direct inputs to Genoid are relative astrometric positions of the satellite with respect to the primary. For Euphrosyne, the algorithm was applied directly to the measured sky-plane offsets derived from SPHERE/ZIMPOL imaging, where the observed positions were measured by fitting two 2D Gaussians to the primary and satellite images (Yang et al., 2020). For very wide binary asteroids, the input astrometry was measured from archival high-angular-resolution imaging, using 2D Gaussian fitting for ground-based AO data and DOLPHOT with pre-computed ACS PSFs for HST data (Minker et al., 19 Mar 2025).

The fitted orbit is parameterized by the seven classical Keplerian elements:

• orbital period $P$
• semi-major axis $a$
• eccentricity $e$
• inclination $i$
• longitude of ascending node $\Omega$
• argument of pericenter $\omega$
• time of pericenter passage $t_p$

In both documented applications, the baseline dynamical model is a Keplerian orbit in the EQJ2000 reference frame (Yang et al., 2020, Minker et al., 19 Mar 2025). For the very wide binary asteroid sample, the authors explicitly state that, for most systems, they restrict the solutions to simple Keplerian orbits, neglecting external perturbation of the Sun and planets and influences due to the non-spherical nature of the systems' primaries because of limited dataset size (Minker et al., 19 Mar 2025). Exceptions are noted for (379) Huenna and (3548) Eurybates, where larger datasets allowed inclusion of solar perturbations or tides (Minker et al., 19 Mar 2025).

A further modeling assumption appears in the wide-binary study: for triple systems, the inner pair was approximated as a single central object for the outer-satellite fit, and the photocenter of the primary was used as the center-of-mass approximation because the position uncertainty was large enough that the difference was negligible for that dataset (Minker et al., 19 Mar 2025).

3. Optimization strategy and fit diagnostics

The algorithmic description given in the literature is concise but consistent. In the Euphrosyne study, Genoid is called a meta-heuristic, implying a global search over orbital parameter space (Yang et al., 2020). In the very distant asteroid satellite study, the mechanics are stated more explicitly: the first generation is drawn randomly over a large range of parameter values, each candidate orbit is propagated to each observing date using Eproc, and fitness is based on $\chi^2$ (Minker et al., 19 Mar 2025). The search used a high number of trial solutions (from 300,000–2,000,000 depending on the complexity of the system) combined with numerous successive generations (1000) (Minker et al., 19 Mar 2025).

The reported fitting goal is to minimize the residuals between observed and computed positions. The papers do not provide a full internal specification of the solver, but the residual-based structure is explicit. A reasonable mathematical representation, stated in the source material for the Euphrosyne application, is

$$\chi^2 = \sum_{k=1}^{N} \left[ \frac{(X_{o,k}-X_{c,k})^2}{\sigma_k^2} + \frac{(Y_{o,k}-Y_{c,k})^2}{\sigma_k^2} \right].$$

For the very wide binary asteroid study, the corresponding summary-level representation is given in terms of right ascension and declination residuals, or equivalently relative sky-plane coordinates:

$$\chi^2 = \sum_i \left[\left(\frac{x_{i,\mathrm{obs}} - x_{i,\mathrm{calc}}}{\sigma_{x,i}}\right)^2 + \left(\frac{y_{i,\mathrm{obs}} - y_{i,\mathrm{calc}}}{\sigma_{y,i}}\right)^2\right].$$

The principal reported quality metric is the root mean square (RMS) of the astrometric residuals in milliarcseconds (Yang et al., 2020). In the wide-binary study, uncertainties are described as statistical, not formal, representing the range of parameter values among solutions whose RMS residuals lie below the average $a$0 threshold (Minker et al., 19 Mar 2025). The authors also caution that such uncertainties may be overestimated because many fitted orbits have RMS values well below that threshold (Minker et al., 19 Mar 2025).

4. Derived quantities and scientific outputs

Genoid is used not only to recover orbital elements but also to derive secondary dynamical and physical quantities. In the Euphrosyne case, the fitted solution yields:

• the mass of Euphrosyne
• the orbital pole in both ecliptic and equatorial coordinates
• the orbit’s inclination relative to Euphrosyne’s equator

The Euphrosyne study reports that the mass inferred from the satellite orbit, combined with the volume from the ADAM 3D-shape model, yields a bulk density of

$a$1

with the corresponding mass

$a$2

(Yang et al., 2020). That density is then used to support the conclusion that Euphrosyne is ice-rich and likely contains a substantial fraction of water ice internally (Yang et al., 2020).

In the wide-binary asteroid study, system mass is likewise derived from orbital period and semimajor axis through Kepler’s third law,

$a$3

although the paper notes that the equation is standard rather than presenting it as part of the algorithm itself (Minker et al., 19 Mar 2025). The orbital fits are also combined with other data types. For Christophedumas and Alconrad, NEATM/WISE/Spitzer photometry was used to estimate diameters, and those size estimates were combined with Genoid-derived masses to infer densities (Minker et al., 19 Mar 2025). For (2577) Litva, the final system characterization combines a Genoid-based orbit for the outer satellite with a lightcurve-based orbit for the inner satellite (Minker et al., 19 Mar 2025).

A plausible implication is that Genoid functions as the dynamical core of a broader inference pipeline: orbit fitting provides mass, while external shape, thermal, or photometric models supply volume or size, enabling density constraints.

5. Documented applications

Two arXiv-documented applications are especially informative because they cover different observational regimes: a newly discovered close satellite and a set of very wide binaries.

System or sample	Data regime	Reported role of Genoid
(31) Euphrosyne / S/2019 (31) 1	Five detections over 26 days from VLT/SPHERE/ZIMPOL in the R filter	Meta-heuristic orbit solver used to fit astrometry, derive orbit, mass, orbital pole, and density (Yang et al., 2020)
Very wide binary asteroids	Sparse archival astrometry from HST, VLT/NACO, Keck/NIRC2, LBT/PISCES	Genetic algorithm used to derive or update orbital solutions, masses, and poles for several systems (Minker et al., 19 Mar 2025)

For Euphrosyne, Genoid used five astrometric detections over a 26-day time span from 2019-03-15, 2019-03-20, 2019-03-25, 2019-03-27, and 2019-04-10 (Yang et al., 2020). The reported orbit is

• $a$4 day
• $a$5 km
• $a$6
• $a$7
• $a$8
• $a$9
• $e$0 JD

with an RMS residual of

$e$1

(Yang et al., 2020). The orbit is described as circular, prograde, and equatorial, with inclination relative to Euphrosyne’s equator

$e$2

(Yang et al., 2020).

For the very wide binary asteroid sample, the reported Genoid-modeled systems include (379) Huenna, (2577) Litva outer satellite, (4674) Pauling, (17246) Christophedumas, (22899) Alconrad, and (3548) Eurybates (Minker et al., 19 Mar 2025). The outputs include the following representative results:

• Huenna: $e$3, time span $e$4 d, RMS $e$5 mas, $e$6 d, $e$7 km, $e$8, $e$9, $i$0 kg (Minker et al., 19 Mar 2025).
• Litva outer satellite: $i$1, time span $i$2 d, RMS $i$3 mas, $i$4 d, $i$5 km, $i$6, $i$7, $i$8 kg (Minker et al., 19 Mar 2025).
• Pauling: $i$9, time span $\Omega$0 d, RMS $\Omega$1 mas, $\Omega$2 d, $\Omega$3 km, $\Omega$4, $\Omega$5, $\Omega$6 kg (Minker et al., 19 Mar 2025).
• Alconrad: $\Omega$7, time span $\Omega$8 d, RMS $\Omega$9 mas, $\omega$0 d, $\omega$1 km, $\omega$2, $\omega$3, $\omega$4 kg (Minker et al., 19 Mar 2025).
• Eurybates: $\omega$5, time span $\omega$6 d, RMS $\omega$7 mas, $\omega$8 d, $\omega$9 km, $t_p$0, $t_p$1, $t_p$2 kg (Minker et al., 19 Mar 2025).

The Christophedumas case is methodologically notable because the formally best residual solution was rejected as physically implausible. The study reports one extremely low-residual solution with RMS $t_p$3 mas but an unrealistically high density of $t_p$4, and a second, constrained, physically reasonable solution with RMS $t_p$5 mas, density $t_p$6, and eccentricity $t_p$7 (Minker et al., 19 Mar 2025). This is a clear example of orbit selection constrained not only by formal residual minimization but also by external physical plausibility.

6. Strengths, limitations, and interpretive cautions

The principal strengths attributed to Genoid are linked to global exploration of a large and potentially degenerate parameter space. The very wide binary asteroid study states that it is well suited for sparse, long-period, and often degenerate systems, and useful when local minimizers can be trapped by aliases or incomplete coverage (Minker et al., 19 Mar 2025). The Euphrosyne study likewise indicates that the method was chosen because the satellite was detected at only a small number of epochs with modest coverage, so a robust search over orbital solutions was needed (Yang et al., 2020).

Several limitations are stated explicitly:

• Sparse datasets: except for Huenna, most systems in the wide-binary sample have few measurements and long gaps, allowing aliasing and broad uncertainty ranges (Minker et al., 19 Mar 2025).
• Short-arc constraints: for Euphrosyne, only five astrometric points over a 26-day arc were available (Yang et al., 2020).
• Poorly constrained eccentricity: the Euphrosyne orbit favors near-circularity, but $t_p$8 shows that eccentricity is not tightly determined (Yang et al., 2020).
• Simplified dynamics: most wide-binary fits are Keplerian only, with solar perturbations, planetary perturbations, and non-spherical gravity neglected for most systems (Minker et al., 19 Mar 2025).
• Uncertainty model: uncertainties are statistical rather than formal covariance errors, and may be overestimated (Minker et al., 19 Mar 2025).
• Non-detections are difficult to encode: for Christophedumas, the authors state that Genoid cannot easily incorporate a “false point” for a non-detection without biasing the fit (Minker et al., 19 Mar 2025).
• Model dependence: the Eurybates solution depends somewhat on whether solar tides are included, and the Christophedumas solution is sensitive to physical plausibility constraints (Minker et al., 19 Mar 2025).

These caveats indicate that Genoid should be understood as a high-leverage search procedure whose outputs remain conditioned on observational completeness, dynamical assumptions, and auxiliary physical constraints.

7. Position within asteroid-dynamics methodology

In the available literature, Genoid appears as a standard workhorse for binary-asteroid orbit determination rather than a fully exposed methodological framework. The Euphrosyne paper cites it as the standard tool in the literature but does not provide algorithmic internals (Yang et al., 2020). The very wide binary asteroid paper gives a more explicit but still concise description centered on random initial populations, propagation with Eproc, $t_p$9-based fitness, and extensive search over candidate solutions (Minker et al., 19 Mar 2025).

Its scientific importance lies in converting relative astrometry into dynamically meaningful quantities. In the Euphrosyne case, Genoid enabled five detections of a $\chi^2$0 km satellite to support a near-circular equatorial orbit, a system mass, and ultimately a density estimate of $\chi^2$1, one of the main results of the study (Yang et al., 2020). In the very wide binary asteroid study, it provided new or updated orbital solutions for several systems and thereby informed hypotheses about the origins of very wide binaries, including arguments that BYORP expansion is unlikely to be solely responsible for their formation and that YORP spin-up followed by scattering, or collisional formation history, may be required for some cases (Minker et al., 19 Mar 2025).

A common misconception is to treat Genoid as synonymous with any genetic optimizer or as a general-purpose optimization framework. The documented usage is narrower: it is an orbit-fitting algorithm for binary and multiple asteroid systems, designed around astrometric residual minimization under orbital dynamics models (Yang et al., 2020, Minker et al., 19 Mar 2025). Another potential confusion arises from the unrelated mathematical and logical term genoid; the categorical construct developed for substitution, lambda calculus, and first-order logic belongs to a separate body of theory and is not the asteroid-orbit algorithm (0712.3088).

In that restricted but important sense, Genoid occupies a specific methodological niche: global orbit determination from sparse relative astrometry, especially where incomplete coverage and parameter degeneracy make purely local fitting unreliable.