$\it COD:$ An Algorithm for Shape Reconstruction of Transiting Celestial Bodies through Topological Optimization

TL;DR

The paper presents $ extit{COD}$, a linear programming-based algorithm for reconstructing the shape of occulting celestial bodies from light curves, demonstrated through simulations and a real astronomical event.

Contribution

Introduces a novel topological optimization algorithm for shape reconstruction of occulting celestial bodies using light-curve data, without assuming specific shapes.

Findings

High accuracy in shape and velocity reconstruction in tests

Successful application to a real astronomical occultation event

Reconstruction of an elliptical opacity distribution for VVV-WIT-08

Abstract

We introduce a novel algorithm, $\textit{COD}$ -- Compact Opacity Distribution, for shape reconstruction of a celestial body that has been observed to occult a star, using the photometric time-series observations of the occultation. $\textit{COD}$ finds a solution to the light-curve inversion problem for an optically thick occulter having an approximately convex shape, together with an estimate of its size, impact parameter and velocity, relative to the occulted star. The algorithm is based on an optimization scheme that uses topological constraints and an objective function for the geometry of the occulter. The constraints of the problem follow linear relations, which enable the use of linear programming optimization as the mathematical framework. Multiple tests of the algorithm were performed, all of which resulted in high correlations between the simulated and obtained shapes of the occulting objects, with errors within $5\%$ in their projected velocities and horizontal sizes, and within $0.1$ in their impact parameters. These tests include a video of a solar eclipse by Phobos, as seen by NASA's Curiosity rover, which was collapsed into its corresponding light curve and reconstructed afterwards. We applied $\textit{COD}$ to the mysterious case of VVV-WIT-08 -- a single deep occultation ($\sim 96 \%$) of a giant star lasting for over 200 days. The analysis, which did not assume any specific shape of the occulter, suggested an object with a projected opacity distribution resembling an ellipse with an eccentricity of $\sim 0.5$, tilted at $\sim 30$ degrees relative to the direction of motion, with a semi-minor axis similar to the stellar radius.