pith. sign in

arxiv: 2604.09209 · v1 · submitted 2026-04-10 · 🌌 astro-ph.GA

An automated method for planetary nebula detection with SIGNALS: first applications to NGC 4214 and NGC 4449

Nancy Yang , Johanna Hartke , Martin Bureau , Chiara Spiniello , Louis-Simon Guit\'e , Guy Flint , Magda Arnaboldi , Ana In\'es Ennis
show 6 more authors
R. Pierre Martin Thomas Martin Carmelle Robert Laurie Rousseau-Nepton Lucas M. Valenzuela S\'ebastien Vicens-Mouret
This is my paper

Pith reviewed 2026-05-10 17:41 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords planetary nebulaeautomated detectionemission-line diagnosticsmorphology testsSIGNALS surveyNGC 4214NGC 4449planetary nebula luminosity function
0
0 comments X

The pith

An automated pipeline detects planetary nebulae in galaxy data with accuracy matching visual inspection methods.

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

The paper develops a new automated pipeline that identifies planetary nebulae by combining emission-line diagnostics with morphology tests on integral-field spectroscopy data. It validates the pipeline through completeness tests that insert mock planetary nebulae with full spectra into the actual data cubes from the SIGNALS survey. The method is then applied to two dwarf irregular galaxies, NGC 4214 and NGC 4449, yielding dozens of detections including many new ones. Distances are derived from the planetary nebula luminosity function and specific frequencies are calculated in both bolometric and V-band terms. A sympathetic reader would care because this automation makes large-scale planetary nebula surveys feasible without relying on time-intensive visual checks.

Core claim

The pipeline uses automated emission-line diagnostics and morphology tests to distinguish planetary nebulae from contaminants at an accuracy level comparable to past visual methods, as demonstrated by completeness and contamination rates measured from mock planetary nebulae inserted into the SITELLE data cubes of NGC 4214 and NGC 4449.

What carries the argument

The automated pipeline that applies emission-line diagnostics and morphology tests to SITELLE integral-field data cubes.

If this is right

  • The pipeline identifies 25 planetary nebulae in NGC 4214, including 6 new discoveries, and 23 in NGC 4449, including 13 new discoveries.
  • Planetary nebula luminosity function distances are measured as 3.09 Mpc for NGC 4214 and 3.91 Mpc for NGC 4449.
  • Bolometric and V-band planetary nebula specific frequencies are calculated for both galaxies based on their total luminosities.

Where Pith is reading between the lines

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

  • The method could be scaled to process the full set of 31 galaxies targeted by the SIGNALS survey without proportional increases in manual effort.
  • The newly defined V-band specific frequency offers a way to compare planetary nebula populations across galaxies using a single observational band.
  • Similar automated diagnostics might be tested on data from other integral-field spectrographs to expand planetary nebula searches beyond the current sample.

Load-bearing premise

The mock planetary nebulae inserted into the data cubes with full spectra accurately represent the spectral and morphological properties of real planetary nebulae.

What would settle it

A large discrepancy between the number of planetary nebulae found by the automated pipeline and the number independently confirmed by visual inspection in a well-studied galaxy would show the claimed accuracy does not hold.

Figures

Figures reproduced from arXiv: 2604.09209 by Ana In\'es Ennis, Carmelle Robert, Chiara Spiniello, Guy Flint, Johanna Hartke, Laurie Rousseau-Nepton, Louis-Simon Guit\'e, Lucas M. Valenzuela, Magda Arnaboldi, Martin Bureau, Nancy Yang, R. Pierre Martin, S\'ebastien Vicens-Mouret, Thomas Martin.

Figure 1
Figure 1. Figure 1: SITELLE SN2 deep images of NGC 4214 (left) and NGC 4449 (right). Blue rectangles are overlaid to indicate the FOV of the corresponding HST images, taken with WFC3 (left) and ACS (right). The red ellipses represent the apertures used by Cook et al. (2014), whose total 𝐵 and 𝑉 magnitudes we adopt. (a) NGC 4214 (b) NGC 4449 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Flux correction factors of NGC 4214 (bottom-left) and NGC 4449 (bottom-right), determined by taking the flux ratios of stars measured using the SITELLE SN2 data cube (integrated over a narrowband HST filter) and an HST image. The sigma-clipped medians are indicated by red dot-dashed vertical lines, ratios of 1 by grey dashed vertical lines. The apparent GAIA 𝐺-band magnitudes of the stars used are shown in… view at source ↗
Figure 3
Figure 3. Figure 3: Fitted map of velocity offsets across the SITELLE FOV, as measured from sky lines in the SN3 filter, for NGC 4214 (left) and NGC 4449 (right). The images cover the same FOV and have the same orientations as those in [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Schematic illustration of the SIGNALS PN identification pipeline. MNRAS 000, 1–18 (2025) [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Inner section of the deep image (left), [O iii] map (centre) and detection map (right) of NGC 4449. The single black pixels in the central panel show spaxels for which the emission-line fit failed. As expected, these are more common toward the edge of the map. and velocity dispersion of each line separately. We did not use this option, however, as we assume PN emission lines are from a single source (super… view at source ↗
Figure 6
Figure 6. Figure 6: Fitting accuracy of mock spectra as a function of initial velocity guess and signal-to-noise ratio. emission-line flux upper limits. For these, we consider whether the direction of the limit makes a source point away from, or toward, the PN-only region. In doing so, we consider the demarcation lines (and the regions they delineate) to be defined only as far as the emission-line flux ratio ranges plotted in… view at source ↗
Figure 7
Figure 7. Figure 7: Classification of the emission-line sources of NGC 4214 (top) and NGC 4449 (bottom), based on the [N ii]-BPT diagram (left) and the [S ii]-BPT diagram (right). Triangular markers indicate limits and their directions, in cases of non-detection of at least one of the emission lines. of a higher precision and can be calibrated using the skyline method described in Section 2.4. The exceptions to this are sourc… view at source ↗
Figure 8
Figure 8. Figure 8: [N ii]-BPT diagnostic diagram, illustrating how a source is classified based on its location and emission-line ratio lower/upper limits. Diagonally pointing arrows indicate an upper limit along two axes. The horizontal line at log( [O iii]/H𝛽) = 0.2 shows the intersection of the Kauffmann et al. (2003) and Sabin et al. (2013) demarcation lines. The horizontal line at log( [O iii]/H𝛽) = 1 and the vertical l… view at source ↗
Figure 9
Figure 9. Figure 9: Distortion maps of NGC 4214 (left) and NGC 4449 (right), fitted to stars. The stars used for the fit are plotted as filled circles, with a colour matching that of the background map and colour table. Lighter colours indicate higher asymmetric distortions. (a) NGC 4214 (b) NGC 4449 [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: 𝜎a maps of NGC 4214 (left) and NGC 4449 (right), fitted to stars. The stars used for the fit are plotted as filled circles, with a colour matching that of the background map and colour table. Lighter colours indicate larger 𝜎𝑎 and therefore more severe size distortions. bright-end cut-off 𝑀∗ 5007 = −4.47 (Ciardullo 2012). The only pa￾rameter we fit for is the apparent magnitude of the bright-end cut-off (… view at source ↗
Figure 11
Figure 11. Figure 11: Comparison of the SN2 (cross-matching) and SN3 (sky-line) velocity calibration methods for sources in NGC 4449. magnitudes from 𝑚 ∗ 5007. It should be noted that with a completeness threshold defined at 50%, the PN catalogue of NGC 4449 is only complete up to 0.5 magnitudes from the fitted bright-end cut-off. 4 RESULTS 4.1 NGC 4214 In NGC 4214, we identify 628 sources with reliably detected [O iii] emissi… view at source ↗
Figure 12
Figure 12. Figure 12: [N ii]-BPT diagrams of the mock PNe generated for NGC 4214 (top) and NGC 4449 (bottom). The mock PNe are colour-coded by input apparent magnitude (left) and galactocentric distance (right). The input emission-line ratio is indicated by the red star. MNRAS 000, 1–18 (2025) [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Recovery rates of mock PNe for NGC 4214 (top) and NGC 4449 (bottom). Left: PNLFs of the mock PNe, overlaid with the recovery fractions. Colours indicate the pipeline step at which the mock PNe fail to be recovered. From top to bottom in the legend: all generated mock PNe, recovered mock PNe with 𝑆/𝑁 ( [O iii] ) > 3, recovered mock PNe which also have the correct emission-line ratios on the BPT diagram(s) … view at source ↗
Figure 14
Figure 14. Figure 14: Left: bona fide and possible PNe of NGC 4214, overlaid on the [O iii] flux map. Right: PNLF of NGC 4214, computed using only the bona fide PNe. The black dotted line shows the best-fitting PNLF. The vertical red dashed line indicates the completeness limit. The likelihood of the apparent bright-end cut-off is plotted underneath [PITH_FULL_IMAGE:figures/full_fig_p015_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Comparison of our NGC 4214 PNLF distances to literature mea￾surements. The red dashed line indicates our fit to the bona fide PNe only, the blue dot-dashed line to both bona fide and possible PNe. Shaded areas represent 1𝜎 uncertainties. References are listed on the right, methods on the left (CMD: colour-magnitude diagram). Distance measurements with black markers were used in calculating the assumed dis… view at source ↗
Figure 16
Figure 16. Figure 16: Same as [PITH_FULL_IMAGE:figures/full_fig_p016_16.png] view at source ↗
Figure 18
Figure 18. Figure 18: Comparison of the 𝛼 parameters derived in this work to those of Buzzoni et al. (2006) and Ennis et al. (2026). The Buzzoni et al. (2006) elliptical galaxies and local group galaxies have been scaled using log(𝛼2.5 ) = log(𝛼8 ) − 1. bad fits can be difficult to identify, as the associated uncertainties are relatively small (𝑆/𝑁 ≈ 2). We have attempted to minimise the impact on our PN catalogue by only quot… view at source ↗
read the original abstract

Utilising the optical imaging Fourier transform spectrograph SITELLE, the Star-formation, Ionized Gas and Nebular Abundances Legacy Survey (SIGNALS) is designed to study the connection between star-forming regions and their environments. Targeting $31$ local star-forming galaxies, its data products also lend themselves to planetary nebula (PN) surveys. We present here a new pipeline to find PNe using automated emission-line diagnostics and morphology tests, that is able to distinguish PNe from contaminants with an accuracy similar to that of past visual methods. We also perform thorough completeness tests using mock PNe inserted into the data cubes with full spectra. We apply these tools to a pilot sample of two dwarf irregular galaxies from the SIGNALS survey, NGC 4214 and NGC 4449, with other galaxies to follow. For these two galaxies, we identify $25$ PNe (including $6$ new discoveries) and $23$ PNe (including $13$ new discoveries), respectively, and calculate PN luminosity function distances of $3.09^{+0.25}_{-0.46}$ and $3.91^{+0.33}_{-0.52}$ Mpc, respectively, the latter consistent with previous estimates. We also calculate the bolometric PN specific frequency of our galaxies ($\alpha_\mathrm{bol}$), as well as a newly defined $V$-band PN specific frequency ($\alpha_\mathrm{V}$) based solely on the galaxies' total luminosities in that band.

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.

Referee Report

2 major / 2 minor

Summary. The paper introduces an automated pipeline for detecting planetary nebulae (PNe) in SITELLE integral-field spectroscopic data from the SIGNALS survey. The method combines emission-line ratio diagnostics with morphological tests to separate PNe from contaminants such as H II regions and supernova remnants. Completeness and contamination rates are quantified via insertion of mock PNe with full spectra into the data cubes. The pipeline is applied to NGC 4214 and NGC 4449, yielding 25 and 23 PNe (including new discoveries), PNLF distances of 3.09^{+0.25}_{-0.46} Mpc and 3.91^{+0.33}_{-0.52} Mpc, and new measurements of bolometric and V-band PN specific frequencies.

Significance. If the automated method can be shown to match visual classification accuracy across a broader range of galaxy types and instrumental conditions, it would enable scalable PN surveys in large IFS datasets, improving distance estimates and constraints on stellar populations in star-forming galaxies. The provision of explicit completeness tests with full-spectrum mocks is a positive step toward reproducibility.

major comments (2)
  1. [Completeness tests and mock insertion procedure] The central performance claim (accuracy comparable to visual methods) rests on completeness and contamination rates measured from mock PNe. The manuscript must demonstrate that the inserted mocks reproduce the observed distributions of [O III]/Hβ, [N II]/Hα, spatial extents, and any SITELLE-specific artifacts (e.g., Fourier sidelobes or PSF variation) in the real data cubes for NGC 4214 and NGC 4449; otherwise the reported rates do not transfer to the actual detections.
  2. [PNLF distance calculation] The PNLF distances are derived using the standard bright-end cutoff; the paper should quantify how uncertainties in the automated classification (false positives/negatives) propagate into the distance errors, particularly for the smaller sample in NGC 4449.
minor comments (2)
  1. [Specific frequency calculations] Define α_bol and α_V explicitly in the text or a table, including the exact luminosity bands and any assumptions about extinction or bolometric corrections.
  2. [Application to NGC 4214 and NGC 4449] Clarify whether the 6 and 13 new discoveries were confirmed by independent visual inspection or solely by the automated pipeline.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful and constructive comments on our manuscript. We have carefully considered each major comment and provide point-by-point responses below. Where appropriate, we have revised the manuscript to incorporate the suggested improvements.

read point-by-point responses

  1. Referee: [Completeness tests and mock insertion procedure] The central performance claim (accuracy comparable to visual methods) rests on completeness and contamination rates measured from mock PNe. The manuscript must demonstrate that the inserted mocks reproduce the observed distributions of [O III]/Hβ, [N II]/Hα, spatial extents, and any SITELLE-specific artifacts (e.g., Fourier sidelobes or PSF variation) in the real data cubes for NGC 4214 and NGC 4449; otherwise the reported rates do not transfer to the actual detections.

    Authors: We agree that demonstrating the fidelity of the mock PNe is important for validating the completeness and contamination rates. Our mock insertion procedure draws line ratios and spatial properties from the observed distributions in the data and inserts them into the real cubes, thereby incorporating SITELLE-specific features like PSF variations and Fourier artifacts. To make this explicit, we have revised the manuscript to include direct comparisons of the key distributions between real and mock PNe, along with quantitative metrics of agreement. These are presented in an expanded Section 3 and a new figure. This addition strengthens the justification for applying the measured rates to our detections. revision: yes

  2. Referee: [PNLF distance calculation] The PNLF distances are derived using the standard bright-end cutoff; the paper should quantify how uncertainties in the automated classification (false positives/negatives) propagate into the distance errors, particularly for the smaller sample in NGC 4449.

    Authors: We thank the referee for highlighting the need to assess the impact of classification uncertainties on the PNLF distances. In the revised manuscript, we have added Monte Carlo simulations that account for potential false positives and negatives by resampling the PN catalogs according to the completeness and contamination estimates from our tests. We refit the PNLF for numerous realizations and derive the resulting spread in distance measurements. These additional uncertainties have been incorporated into the final error bars reported for both galaxies. Although the sample in NGC 4449 is slightly smaller, the effect remains comparable to the other sources of uncertainty and does not alter our conclusions. The details are now described in Section 4.2. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper introduces an automated PN detection pipeline based on standard emission-line diagnostics and morphology tests, with completeness quantified via mock PNe inserted into SITELLE cubes. No equations or steps reduce the reported accuracies, distances, or specific frequencies to quantities defined by the authors' own fitted parameters or self-citations. The PN luminosity function distances and bolometric/V-band frequencies follow established external methods without self-referential closure. Validation against visual methods and prior distance estimates remains independent of the pipeline's internal definitions.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the domain assumption that planetary nebulae occupy a distinct region in emission-line diagnostic space and on the standard use of the PN luminosity function as a distance indicator; no free parameters or invented entities are explicitly introduced in the abstract.

axioms (2)
  • domain assumption Planetary nebulae can be reliably separated from HII regions, supernova remnants, and background galaxies using a combination of emission-line ratios and morphological criteria.
    This assumption underpins the automated diagnostics and is required for the claimed accuracy to hold.
  • domain assumption The PN luminosity function provides a standard candle for distance estimation once a sufficient number of PNe are detected.
    Used directly to derive the reported distances of 3.09 and 3.91 Mpc.

pith-pipeline@v0.9.0 · 5641 in / 1451 out tokens · 27468 ms · 2026-05-10T17:41:58.369856+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

  1. [1]

    B., et al., 2017, A&A, 600, A90 Caminha G

    Annibali F., Aloisi A., Mack J., Tosi M., van der Marel R. P., Angeretti L., Leitherer C., Sirianni M., 2008, AJ, 135, 1900 Annibali F., et al., 2017, ApJ, 843, 20 Arnaboldi M., Freeman K. C., Gerhard O., Matthias M., Kudritzki R. P., Méndez R. H., Capaccioli M., Ford H., 1998, ApJ, 507, 759 Baldwin J. A., Phillips M. M., Terlevich R., 1981, Publications ...

    work page doi:10.5281/zenodo.6825092 2008