THOR: An Algorithm for Cadence-Independent Asteroid Discovery

Abstract

We present "Tracklet-less Heliocentric Orbit Recovery" (THOR), an algorithm for linking of observations of Solar System objects across multiple epochs that does not require intra-night tracklets or a predefined cadence of observations within a search window. By sparsely covering regions of interest in the phase space with "test orbits", transforming nearby observations over a few nights into the co-rotating frame of the test orbit at each epoch, and then performing a generalized Hough transform on the transformed detections followed by orbit determination (OD) filtering, candidate clusters of observations belonging to the same objects can be recovered at moderate computational cost and little to no constraints on cadence. We validate the effectiveness of this approach by running on simulations as well as on real data from the Zwicky Transient Facility (ZTF). Applied to a short, 2-week, slice of ZTF observations, we demonstrate THOR can recover 97.4% of all previously known and discoverable objects in the targeted ($a>1.7$ au) population with 5 or more observations and with purity between 97.7% and 100%. This includes 10 likely new discoveries, and a recovery of an $e\sim 1$ comet C/2018 U1 (the comet would have been a ZTF discovery had THOR been running in 2018 when the data were taken). The THOR package and demo Jupyter notebooks are open source and available at https://github.com/moeyensj/thor.

arXiv comment22 pages, 11 figures

References in corpus (16)

• Array Programming with NumPy
• Astropy: A Community Python Package for Astronomy
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• The Zwicky Transient Facility: System Overview, Performance, and First Results
• The Zwicky Transient Facility: Data Processing, Products, and Archive
• The Zwicky Transient Facility: Science Objectives
• A record of planet migration in the Main Asteroid Belt
• The Pan-STARRS Moving Object Processing System
• Orbit Fitting and Uncertainties for Kuiper Belt Objects
• The Zwicky Transient Facility Alert Distribution System
• The Large Synoptic Survey Telescope as a Near-Earth Object Discovery Machine
• Efficient intra- and inter-night linking of asteroid detections using kd-trees
• High-fidelity Simulations of the Near-Earth Object Search Performance of the Large Synoptic Survey Telescope
• Near-Earth Object Orbit Linking with the Large Synoptic Survey Telescope

Cited by in corpus (2)

• Give Me a Few Hours: Exploring Short Timescales in Rubin Observatory Cadence Simulations
• Exploring the Solar System with the NOIRLab Source Catalog I: Detecting Objects with CANFind