Use the 4S (Signal-Safe Speckle Subtraction): Explainable Machine Learning reveals the Giant Exoplanet AF Lep b in High-Contrast Imaging Data from 2011

Markus J. Bonse, Timothy D. Gebhard, Felix A. Dannert, Olivier Absil, Faustine Cantalloube, Valentin Christiaens, Gabriele Cugno, Emily O. Garvin, Jean Hayoz, Markus Kasper, Elisabeth Matthews, Bernhard Schölkopf, Sascha P. Quanz

The Astronomical Journal, 169 (194),

Preprint PDF BibTeX

Cite this paper

@article{Bonse_2025,
  author         = {{Bonse}, Markus J. and {Gebhard}, Timothy D. and {Dannert}, Felix A. and {Absil}, Olivier and {Cantalloube}, Faustine and {Christiaens}, Valentin and {Cugno},
Gabriele and {Garvin}, Emily O. and {Hayoz}, Jean and {Kasper}, Markus and {Matthews}, Elisabeth and {Schölkopf}, Bernhard and {Quanz}, Sascha P.},
  year           = 2025,
  month          = 4,
  title          = {{Use the 4S (Signal-Safe Speckle Subtraction): Explainable Machine Learning Reveals the Giant Exoplanet AF Lep b in High-contrast Imaging Data from 2011}},
  journal        = {\aj},
  day            = {1},
  doi            = {10.3847/1538-3881/adab79},
  number         = 4,
  pages          = {194},
  volume         = 169,
}
NASA/ADS Code DOI

Abstract:

The main challenge of exoplanet high-contrast imaging (HCI) is to separate the signal of exoplanets from their host stars, which are many orders of magnitude brighter. HCI for ground-based observations is further exacerbated by speckle noise originating from perturbations in Earth's atmosphere and imperfections in the telescope optics. Various data postprocessing techniques are used to remove this speckle noise and reveal the faint planet signal. Often, however, a significant part of the planet signal is accidentally subtracted together with the noise. In the present work, we use explainable machine learning to investigate the reason for the loss of the planet signal for one of the most used postprocessing methods: principal component analysis (PCA). We find that PCA learns the shape of the telescope point-spread function for high numbers of PCA components. This representation of the noise captures not only the speckle noise but also the characteristic shape of the planet signal. Building on these insights, we develop a new postprocessing method (4S) that constrains the noise model to minimize this signal loss. We apply our model to 11 archival HCI data sets from the Very Large Telescope NACO instrument in the L' band and find that our model consistently outperforms PCA. The improvement is largest at close separations to the star (≤ 4λ/D), providing up to 1.5 mag deeper contrast. This enhancement enables us to detect the exoplanet AF Lep b in data from 2011, 11 yr before its subsequent discovery. We present updated orbital parameters for this object.