arxiv: 2604.16731 · v1 · submitted 2026-04-17 · 🌌 astro-ph.EP · astro-ph.IM

Recognition: unknown

heliostack: A Novel Approach to Minor Planet Discovery

Kevin J Napier , Matthew J Holman , Hsing-Wen Lin , David W Gerdes , Thomas R Ruch

Authors on Pith no claims yet

Pith reviewed 2026-05-10 06:43 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.IM

keywords heliostackKuiper Belt Objectsshift-and-stackminor planet discoveryHubble Space Telescopenonlinear motionfaint solar system objects

0 comments

The pith

Heliostack is an algorithm that discovers faint minor planets by nonlinearly stacking images over time spans longer than one day.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces heliostack as a method to search for faint solar system objects by stacking images while accounting for their nonlinear motion across longer periods. This matters because conventional surveys struggle with faint objects and existing stacking methods were limited to short time spans where signals would smear. By applying it to Hubble Space Telescope data spanning 15 days, the authors recover previously known objects below detection thresholds and identify two new ones. These represent the first discoveries made from such extended stacks. The approach promises to unlock more discoveries from archival and future survey data.

Core claim

Heliostack is introduced as an algorithm for nonlinear shift-and-stack searches for solar system objects. When applied to Hubble Space Telescope images spanning 15 days, it successfully recovers both known sub-threshold objects and adds two new discoveries. These two new objects are the first to be discovered in stacks of images taken over a time span longer than about one day.

What carries the argument

heliostack, an algorithm for nonlinear shift-and-stack image combination that models the curved paths of solar system objects to align and stack images over extended time baselines.

If this is right

Existing archival images from surveys can be reanalyzed for additional faint object detections.
Searches can now use longer time baselines, increasing sensitivity to slow-moving or faint objects.
Forthcoming surveys can apply the method to maximize their potential for minor planet discoveries.
The technique supports exhaustive searches for specific populations such as Cold Classical Kuiper Belt Objects.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The method could be adapted to ground-based surveys that face similar long-baseline motion challenges.
Routine use in data pipelines for large surveys would allow routine detection of sub-threshold sources.
A larger sample of such faint objects could refine models of solar system formation.

Load-bearing premise

The nonlinear motion model and detection thresholds accurately distinguish real objects from noise and artifacts in the 15-day HST data without significant false positives or missed detections.

What would settle it

Independent follow-up observations that fail to confirm the two new objects, or blind tests on similar data that produce many unconfirmed candidates, would show the thresholds do not reliably separate real objects from noise.

Figures

Figures reproduced from arXiv: 2604.16731 by David W Gerdes, Hsing-Wen Lin, Kevin J Napier, Matthew J Holman, Thomas R Ruch.

**Figure 1.** Figure 1: (Left): The canvas (approximately 14,000 by 13,000 pixels) after stacking all 1118 images for some initial condition. The inset plot shows the lone peak above 7σ in the stack, which turns out to be a 21σ detection of 2003 BH91. Since the asserted parameters are not quite correct, the detection is somewhat smeared. (Right): A heatmap of the number of images that contributed to each pixel in the significance… view at source ↗

**Figure 2.** Figure 2: (Top) Distribution of the signal-to-noise ratio of all of the unique candidate trajectories resulting from the heliostack search, with ν ≥ 7. The green histogram represents the implanted synthetic sources, and the gray histogram represents real candidate moving objects. The dashed black line shows the shape of a Gaussian distribution with mean 0 and width 1. It is clear that the distribution of ν is very n… view at source ↗

**Figure 3.** Figure 3: (Top): Histogram of the mean apparent magnitude of the implanted synthetic sources, with an accounting of whether or not the sources were recovered. (Bottom): Best fit of our efficiency to Equation (20), including 68% and 95% confidence intervals. The best-fit parameters are η0 = 1.0, m50 = 29.21 ± 0.03, and σ = 0.11 ± 0.02 [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗

**Figure 4.** Figure 4: From left to right, stacked cutouts of 2003 BF91, 2003 BH91, 2003 ABCD, and 2003 WXYZ. (2004) reported 28.15 ± 0.04. Similarly, we measure the mean magnitude of 2003 BH91 to be 28.16 ± 0.04, while G. M. Bernstein et al. (2004) reported 28.38 ± 0.05. Second, we have confident detections of both 2003 ABCD and 2003 WXYZ, while G. M. Bernstein et al. (2004) did not find either object. One possible reason for t… view at source ↗

**Figure 5.** Figure 5: Cutouts of each of the recovered sources. Objects that are not implanted synthetic sources are highlighted in red [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗

**Figure 6.** Figure 6: Trajectory of 2003 BF91 on the sky as seen by HST over the course of the survey. The inset plot emphasizes how nonlinear the trajectory is, due to HST’s low Earth orbit; the proper motion in the sky plane remains linear for at most a few minutes at a time. The red square shows the size of one pixel; it is clear that the parallax loops are too large to simply ignore. The heliostack algorithm automatically c… view at source ↗

read the original abstract

The study of faint solar system objects is a promising avenue for understanding the origin and evolution of planetary systems. However, such objects are difficult to detect in conventional surveys. Here we introduce heliostack, an algorithm for nonlinear shift-and-stack searches for solar system objects, which enables us to combine images taken over longer time spans than was previously possible. Applying this algorithm to a number of existing archival and forthcoming surveys will allow us to maximize their potential for discovering faint solar system objects. In this work, we apply heliostack to archival Hubble Space Telescope (HST) data, completing an exhaustive search for Cold Classical Kuiper Belt Objects in a set of HST images taken over a 15-day time span in 2003. We successfully recover both of the known sub-threshold objects in the data, and add two new discoveries. These two new objects are the first to ever be discovered in stacks of images taken over a time span longer than about one day.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Heliostack extends shift-and-stack to nonlinear motion over 15-day baselines and pulls two new Cold Classical KBOs from 2003 HST frames.

read the letter

The main point is that this gives a practical way to stack images over much longer times than the usual one-night limit, and the authors recover the two known sub-threshold objects while reporting two new discoveries as the first ever found that way. They apply the method exhaustively to the archival HST set and position it as something that can be dropped onto other surveys to boost faint-object yields. That is the concrete advance: the longer baseline without losing the ability to combine signal. The recoveries of the known objects supply direct empirical support for the approach working on real data. The new finds are presented as independent outcomes of the search rather than fitted results. Soft spots sit in the implementation details. The abstract gives no numbers on false-positive rates, threshold choices, or how the nonlinear motion model was validated against noise and artifacts, so it is hard to judge how cleanly the new objects separate from background. If the full paper has only limited simulations or no code release, that would be the section needing the most referee attention. The argument itself is not circular; the discoveries follow from applying the algorithm to the data. This is for people running or reprocessing minor-planet surveys who want to squeeze more out of existing image sets. A reader focused on detection techniques or Kuiper Belt population studies would get direct value from the method description and the example results. The work is grounded enough and addresses a clear practical gap that it deserves a serious referee rather than a desk reject.

Referee Report

3 major / 2 minor

Summary. The manuscript introduces heliostack, a novel algorithm for nonlinear shift-and-stack searches that enables detection of faint solar system objects (such as Cold Classical Kuiper Belt Objects) in image stacks spanning longer time baselines than previously feasible (up to 15 days). The authors apply the method to archival HST data from 2003, claiming successful recovery of both known sub-threshold objects plus two new discoveries, which are presented as the first minor planets found in stacks exceeding ~1-day baselines.

Significance. If the central claims hold, heliostack could meaningfully expand the discovery space for sub-threshold solar system bodies by leveraging existing archival surveys and forthcoming data sets with extended temporal coverage. The empirical recovery of known objects supplies direct (if limited) support for the approach; however, the absence of detailed validation leaves the significance conditional on further demonstration of robustness against noise and artifacts.

major comments (3)

[Abstract and §3] Abstract and §3 (Algorithm): The nonlinear motion model and detection thresholds are described at a high level but without explicit equations, parameter definitions, or optimization procedure. This omission is load-bearing because the weakest assumption (accurate separation of real objects from noise over 15-day baselines) cannot be evaluated without the mathematical formulation.
[§4 and §5] §4 (Application to HST data) and §5 (Results): No false-positive rate, completeness simulations, or error analysis is reported despite the claim of two new discoveries. The recovery of known objects is stated but not quantified (e.g., no S/N values, positional residuals, or comparison to independent detections), undermining the assertion that the new objects are reliable.
[§5] §5 (Results): The statement that the two new objects are 'the first to ever be discovered in stacks... longer than about one day' requires a clear literature comparison or citation to prior shift-and-stack limits; without it, the novelty claim rests on an unverified assertion.

minor comments (2)

[§3] Notation for the heliocentric motion parameters is introduced without a consistent symbol table or reference to prior work on nonlinear stacking.
[§5] Figure captions and axis labels in the results section should explicitly indicate which panels show the recovered known objects versus the new candidates.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed comments, which have identified important areas for improving the clarity, rigor, and validation of our manuscript on heliostack. We address each major comment point by point below and will incorporate revisions to strengthen the paper.

read point-by-point responses

Referee: [Abstract and §3] Abstract and §3 (Algorithm): The nonlinear motion model and detection thresholds are described at a high level but without explicit equations, parameter definitions, or optimization procedure. This omission is load-bearing because the weakest assumption (accurate separation of real objects from noise over 15-day baselines) cannot be evaluated without the mathematical formulation.

Authors: We agree that the absence of explicit mathematical details limits the ability to fully evaluate the method. In the revised manuscript, we will expand §3 to include the complete set of equations for the nonlinear motion model, which uses a quadratic approximation to the apparent heliocentric motion over the 15-day baseline. All parameters will be defined, including the grid search ranges for velocity and acceleration components, and the optimization procedure will be described as a maximization of the stacked signal-to-noise ratio subject to the detection threshold derived from the noise statistics in the data. revision: yes
Referee: [§4 and §5] §4 (Application to HST data) and §5 (Results): No false-positive rate, completeness simulations, or error analysis is reported despite the claim of two new discoveries. The recovery of known objects is stated but not quantified (e.g., no S/N values, positional residuals, or comparison to independent detections), undermining the assertion that the new objects are reliable.

Authors: We acknowledge that quantitative validation metrics are essential to support the discovery claims. In the revision, we will add to §§4 and 5: (i) completeness simulations via injection of synthetic objects with realistic Cold Classical KBO motions and magnitudes into the HST images, reporting recovery rates versus S/N; (ii) an estimate of the false-positive rate based on the distribution of noise peaks in the stacked images; and (iii) for both known and new objects, tabulated S/N values, positional residuals relative to predicted ephemerides, and cross-checks against independent detections where available. This will provide a clearer assessment of reliability. revision: yes
Referee: [§5] §5 (Results): The statement that the two new objects are 'the first to ever be discovered in stacks... longer than about one day' requires a clear literature comparison or citation to prior shift-and-stack limits; without it, the novelty claim rests on an unverified assertion.

Authors: We will revise the statement in §5 (and update the abstract) to include a direct literature comparison. We will cite prior shift-and-stack studies for Kuiper Belt objects, which have generally been restricted to baselines of approximately one day or shorter due to the effects of nonlinear motion, and explain how the heliostack nonlinear model extends this capability to 15 days. This will properly contextualize the novelty of the two new discoveries. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper introduces heliostack as a novel algorithmic extension for nonlinear shift-and-stack detection over multi-day baselines and demonstrates it empirically by recovering two known sub-threshold objects plus two new discoveries in 15-day HST archival data. No derivation chain, equations, or first-principles results are presented that reduce by construction to fitted inputs, self-definitions, or self-citation load-bearing premises. The central claims rest on the empirical outcomes of applying the method to external data rather than on any internal renaming, ansatz smuggling, or uniqueness theorem imported from prior author work. The recovery of known objects supplies direct, non-circular validation of the detection thresholds and motion model.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No specific free parameters, axioms, or invented entities are identifiable from the abstract alone. The algorithm description implies motion modeling parameters but provides no details.

pith-pipeline@v0.9.0 · 5476 in / 1243 out tokens · 38954 ms · 2026-05-10T06:43:01.785972+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

28 extracted references · 21 canonical work pages

[1]

L., Bernstein, G

Allen, R. L., Bernstein, G. M., & Malhotra, R. 2002, AJ, 124, 2949, doi: 10.1086/343773

work page doi:10.1086/343773 2002
[2]

2023, The ACS/WFC Focus-Diverse ePSF Webtool,, Instrument Science Report ACS 2023-06, 12 pages

Bellini, A. 2023, The ACS/WFC Focus-Diverse ePSF Webtool,, Instrument Science Report ACS 2023-06, 12 pages

2023
[3]

F., Durda, D

Batygin, K., Adams, F. C., Brown, M. E., & Becker, J. C. 2019, PhR, 805, 1, doi: 10.1016/j.physrep.2019.01.009

work page doi:10.1016/j.physrep.2019.01.009 2019
[4]

M., Trilling , D

Bernstein, G. M., Trilling, D. E., Allen, R. L., et al. 2004, AJ, 128, 1364, doi: 10.1086/422919

work page doi:10.1086/422919 2004
[5]

M., Trilling, D

Bernstein, G. M., Trilling, D. E., Allen, R. L., et al. 2006, AJ, 131, 2364, doi: 10.1086/503194

work page doi:10.1086/503194 2006
[6]

W., Spencer, J

Buie, M. W., Spencer, J. R., Porter, S. B., et al. 2024, PSJ, 5, 196, doi: 10.3847/PSJ/ad676d

work page doi:10.3847/psj/ad676d 2024
[7]

Danby, J. M. A. 1992, Fundamentals of Celestial Mechanics (Richmond, VA: Willmann-Bell, Inc.) Dark Energy Survey Collaboration, Abbott, T., Abdalla, F. B., et al. 2016, MNRAS, 460, 1270, doi: 10.1093/mnras/stw641

work page doi:10.1093/mnras/stw641 1992
[8]

J., Nicholson, P

Gladman, B., Kavelaars, J. J., Nicholson, P. D., Loredo, T. J., & Burns, J. A. 1998, AJ, 116, 2042, doi: 10.1086/300573

work page doi:10.1086/300573 1998
[9]

2025, Astronomy and Computing, 53, 100987, doi: 10.1016/j.ascom.2025.100987

Golovich, N., Steil, T., Geringer-Sameth, A., et al. 2025, Astronomy and Computing, 53, 100987, doi: 10.1016/j.ascom.2025.100987

work page doi:10.1016/j.ascom.2025.100987 2025
[10]

2004, Multiple View Geometry in Computer Vision, 2nd edn

Hartley, R., & Zisserman, A. 2004, Multiple View Geometry in Computer Vision, 2nd edn. (Cambridge University Press)

2004
[11]

2018, AJ, 156, 135, doi: 10.3847/1538-3881/aad69a

Kuindersma, S. 2018, AJ, 156, 135, doi: 10.3847/1538-3881/aad69a

work page doi:10.3847/1538-3881/aad69a 2018
[12]

Research Notes of the American Astronomical Society , keywords =

Holman, M. J., Payne, M. J., & P´ al, A. 2019, Research Notes of the American Astronomical Society, 3, 160, doi: 10.3847/2515-5172/ab4ea6

work page doi:10.3847/2515-5172/ab4ea6 2019
[13]

J., Kavelaars, J

Holman, M. J., Kavelaars, J. J., Grav, T., et al. 2004, Nature, 430, 865, doi: 10.1038/nature02832

work page doi:10.1038/nature02832 2004
[14]

2022, ApJS, 259, 54, doi: 10.3847/1538-4365/ac559a

Huang, Y., Gladman, B., & Volk, K. 2022, ApJS, 259, 54, doi: 10.3847/1538-4365/ac559a

work page doi:10.3847/1538-4365/ac559a 2022
[15]

Huber, P. J. 1964, The Annals of Mathematical Statistics, 35, 73 , doi: 10.1214/aoms/1177703732 Juri´ c, M., Jones, R. L., Bryce Kalmbach, J., et al. 2019, arXiv e-prints, arXiv:1901.08549, doi: 10.48550/arXiv.1901.08549

work page doi:10.1214/aoms/1177703732 1964
[16]

J., Petit, J.-M., Gladman, B., et al

Kavelaars, J. J., Petit, J.-M., Gladman, B., et al. 2021, ApJL, 920, L28, doi: 10.3847/2041-8213/ac2c72

work page doi:10.3847/2041-8213/ac2c72 2021
[17]

Lang, D., & Hogg, D. W. 2026, AJ, 171, 231, doi: 10.3847/1538-3881/ae4505

work page doi:10.3847/1538-3881/ae4505 2026
[18]

F., et al

Levison, H. F., et al. 2021, PSJ, 2

2021
[19]

Manzano, L. E. S., Gerdes, D. W., Napier, K. J., et al. 2025, PSJ, 6, 213, doi: 10.3847/PSJ/adf5c2

work page doi:10.3847/psj/adf5c2 2025
[20]

J., & Holman, M

Napier, K. J., & Holman, M. J. 2024, PSJ, 5, 270, doi: 10.3847/PSJ/ad938f

work page doi:10.3847/psj/ad938f 2024
[21]

J., Lin, H

Napier, K. J., Lin, H. W., Gerdes, D. W., et al. 2024, PSJ, 5, 50, doi: 10.3847/PSJ/ad1528

work page doi:10.3847/psj/ad1528 2024
[22]

R., Napier, K

Ruch, T. R., Napier, K. J., et al. 2026, in prep 14

2026
[23]

J., Kalmbach, J

Smotherman, H., Connolly, A. J., Kalmbach, J. B., et al. 2021, AJ, 162, 245, doi: 10.3847/1538-3881/ac22ff

work page doi:10.3847/1538-3881/ac22ff 2021
[24]

A., Fraser, W

Stansberry, J. A., Fraser, W. C., Trilling, D. E., et al. 2021, An Ultra-Sensitive Pencil Beam Search for 10 km Trans-Neptunian Objects,, JWST Proposal. Cycle 1, ID. #1568

2021
[26]

E., Gerdes, D

Trilling, D. E., Gerdes, D. W., Juri´ c, M., et al. 2024b, AJ, 167, 132, doi: 10.3847/1538-3881/ad1529

work page doi:10.3847/1538-3881/ad1529
[27]

E., Fraser, W

Trilling, D. E., Fraser, W. C., Bernstein, G., et al. 2025, The supremely deep trans-Neptunian object survey: A critical test of planet formation models,, JWST Proposal. Cycle 4, ID. #7700

2025
[28]

A., Fuentes, C., Gerdes, D

Trujillo, C. A., Fuentes, C., Gerdes, D. W., et al. 2024, AJ, 167, 133, doi: 10.3847/1538-3881/ad1523

work page doi:10.3847/1538-3881/ad1523 2024
[29]

A., Guhathakurta, P., Bernstein, G

Tyson, J. A., Guhathakurta, P., Bernstein, G. M., & Hut, P. 1992, in American Astronomical Society Meeting

1992