Shape, Orientation and Colors Combined approach for Asteroids (SOCCA)

Technical Summary: Shape, Orientation and Colors Combined approach for Asteroids (SOCCA)

Problem Statement
Large photometric surveys (e.g., ZTF, LSST) provide sparse, multi-band photometry for millions of Solar System Objects (SSOs), offering a unique opportunity to constrain their physical and compositional properties simultaneously. However, current methodologies face significant limitations:

1. Phase Function Models: Standard models (e.g., $HG_1G_2$) do not account for rotational variability or the changing aspect angle of the target, leading to biased phase parameters and inaccurate absolute magnitudes ($H$).
2. Shape and Spin Inversion: Traditional light curve inversion (LCI) methods are computationally expensive and often require dense time-series data, making them prohibitive for large-scale sparse datasets.
3. Data Scarcity: While taxonomic classifications are increasing, physical parameters (spin state, shape, density) remain constrained for a small fraction of the known SSO population, hindering the study of long-term dynamical evolution (e.g., the Yarkovsky effect).

Methodology
The authors introduce SOCCA (Shape, Orientation and Colors Combined approach for Asteroids), a model designed to jointly retrieve absolute magnitudes, phase parameters, spin state, and shape proportions from sparse, multi-filter photometry.

• Model Formulation: SOCCA extends the $HG_1G_2$ formalism by incorporating the projected surface area of a rotating triaxial ellipsoid. The apparent magnitude $m$ is modeled as:
  $$m = H + f(r, \Delta) + g(\gamma) + s(\alpha, \delta, t)$$
  Where $f(r, \Delta)$ accounts for distance, $g(\gamma)$ is the phase function, and $s(\alpha, \delta, t)$ is a time-dependent shape function describing the projected area of a triaxial ellipsoid ($a > b > c$) rotating with a sidereal period $P_{sid}$ around a spin axis $(\alpha_0, \delta_0)$.
• Parameterization: The model utilizes a minimal set of free parameters:
  • 3 parameters for phase behavior ($H, G_1, G_2$) per filter.
  • 6 geometric parameters independent of wavelength: spin axis orientation ($\alpha_0, \delta_0$), rotation parameters ($P_{sid}, W_0$), and shape ratios ($a/b, a/c$).
• Initialization Strategy: To avoid computationally expensive grid searches, SOCCA employs a frugal initialization pipeline:
  1. Phase/Shape: Initial $H, G_1, G_2$ are derived from a standard $sHG_1G_2$ fit. Shape ratios ($a/b, a/c$) are estimated from the peak-to-peak amplitude of the $sHG_1G_2$ residuals.
  2. Spin Axis: A grid search over the celestial sphere identifies local minima in the $sHG_1G_2$ residuals to seed the spin axis coordinates.
  3. Rotational Period: A two-step process is used:
     • Synodic Period: Lomb-Scargle frequency analysis is applied to the $sHG_1G_2$ residuals. A bootstrap resampling technique assesses stability, and an alias-flagging scheme (based on the 24-hour observing cadence) distinguishes true periods from aliases.
     • Sidereal Period: The synodic period is converted to a sidereal period ($P_{sid}$) by accounting for the target's orbital motion and spin axis orientation. A search window is defined based on the semi-major axis and the survey baseline to resolve period aliases.
• Optimization: The final parameters are determined via gradient-descent least-squares minimization. Model selection between complexity levels (e.g., $HG_1G_2$ vs. SOCCA) is governed by an F-test to ensure statistical significance.

Key Results
The model was validated using a 10-year LSST simulation of 2,207 SSOs and real data from the asteroid (45) Eugenia observed by ZTF.

• Success Rate: SOCCA achieved a success rate of 52.8% for the simulated sample. The inclusion of shape and rotation increased the number of physically meaningful solutions by $\sim$10–20% per filter compared to previous models.
• Fit Quality: SOCCA significantly improved the fit to photometric data, reducing mean residuals to 0.07 mag, compared to 0.11 mag for $sHG_1G_2$ and 0.14 mag for $HG_1G_2$. The reduced chi-squared ($\chi^2_\nu$) dropped from 108.9 ($HG_1G_2$) to 8.81 (SOCCA).
• Parameter Accuracy:
  • Absolute Magnitude ($H$): Retrieved with a scatter approximately three times smaller than existing approaches.
  • Phase Parameters ($G_1, G_2$): Improved determination by a factor of three, resulting in a sharper distribution in the $G_1G_2$ plane that better separates taxonomic classes (e.g., S- and C-types).
  • Spin Axis: 68.5% of solutions were within $15^\circ$ of the true obliquity.
  • Shape: Axis ratios $a/b$ and $a/c$ were recovered with success rates of 83.3% and 67.4%, respectively (within 20% of true values).
  • Period: After alias correction and phase dispersion minimization, the purity of true period estimates increased to 75.8%.
• Computational Efficiency: The median runtime was ~157 seconds per object on a single core, scaling linearly with the number of observations. This makes the method scalable to datasets of several hundred thousand objects.
• Real Data Validation: For (45) Eugenia, SOCCA retrieved a sidereal period of 5.699,136 hours (vs. 5.699,152 hours from 3D reconstruction), spin axis coordinates within $4^\circ$ of established values, and shape ratios consistent with 3D models.

Significance and Claims
The authors claim that SOCCA provides a more complete physical description of SSOs from sparse photometry by unifying phase, shape, and rotational information into a single, computationally efficient model.

• Scalability: The method is explicitly designed for current and upcoming large surveys (ZTF, LSST), enabling population-level studies rather than individual object analysis.
• Physical Constraints: By accurately retrieving spin and shape, SOCCA provides the necessary physical parameters to study the Yarkovsky and YORP effects, which drive the long-term dynamical evolution and diffusion of asteroid families.
• Taxonomic Utility: The improved accuracy of absolute magnitudes and phase parameters offers a more robust basis for taxonomic classification of large datasets.
• Frugality: The approach balances physical realism (triaxial ellipsoid) with parameter economy, avoiding the over-parameterization of full 3D shape models while outperforming simple phase curve fits.

The paper concludes that SOCCA is a practical tool for shifting SSO characterization from individual case studies to a comprehensive, statistical description of the Solar System population.