arxiv: 2604.16045 · v1 · submitted 2026-04-17 · 🌌 astro-ph.GA

Kinematically cold and warm planetary nebulae samples, HII regions and supernovae remnants in the disc of the face-on spiral galaxy NGC 628 (M74) -- The Planetary Nebulae Spectrograph with the Hα arm

Magda Arnaboldi , Ortwin Gerhard , Surya Aniyan , Kenneth C. Freeman , Anastasia Ponomareva , Lodovico Coccato , Johanna Hartke , Steven P. Bamford

show 9 more authors

Arianna Cortesi Nigel Douglas Crescenzo Tortora Michael Merrifield Konrad Kuijken Massimo Capaccioli Nicola R. Napolitano Claudia Pulsoni Aaron J. Romanowsky

This is my paper

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

classification 🌌 astro-ph.GA

keywords planetary nebulaeNGC 628velocity dispersionmaximal discrotation curveHII regionssupernova remnantsgalactic kinematics

0 comments

The pith

Two planetary nebula populations in NGC 628 show a maximal disc where baryons supply 78 percent of the rotational velocity.

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

The paper uses the Planetary Nebulae Spectrograph with an added H-alpha arm to observe the face-on spiral NGC 628 and separate planetary nebulae from HII regions and supernova remnants via narrow-band colors. It identifies two distinct kinematic groups among the 442 planetary nebulae: a cold population with low velocity dispersion perpendicular to the disc that dominates the bright end of the luminosity function, and a warmer population with higher dispersion that dominates at fainter magnitudes. By matching the warm component's dispersion to the expected vertical scale height of older stars, the authors decompose the galaxy's rotation curve. The resulting model shows the baryonic component accounts for 78 percent of the total rotational velocity across the studied radii.

Core claim

The study finds two kinematically distinct planetary nebula populations in the disc of NGC 628. The cold population has a velocity dispersion orthogonal to the plane of sigma_z,cold = 8.8 km/s and dominates the planetary nebula luminosity function near the bright cut-off, while the warm population has sigma_z,warm = 26.1 km/s and increasingly dominates at fainter magnitudes. These populations contribute 46 percent and 54 percent respectively. Matching the velocity dispersion of the old warm component with the stellar population's scale height allows a rotation-curve decomposition in which the baryonic component supplies 78 percent of the total rotational velocity, indicating a maximal disc.

What carries the argument

Separation of planetary nebulae into cold and warm populations by their line-of-sight velocity dispersions perpendicular to the galactic plane, used to trace distinct stellar age groups and to constrain the rotation-curve decomposition.

If this is right

The cold planetary-nebula population traces younger, more massive progenitors that set the bright cut-off of the planetary nebula luminosity function.
The warm planetary-nebula population traces older stars whose greater vertical thickness is reflected in their higher velocity dispersion.
The rotation-curve decomposition yields a maximal-disc model in which baryons dominate the inner rotation of NGC 628.
The two populations together allow a direct kinematic probe of the vertical structure of the stellar disc without photometric assumptions alone.

Where Pith is reading between the lines

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

The same PN.S plus H-alpha technique could be applied to other face-on spirals to test whether maximal discs are typical.
If the scale-height matching is reliable, dark-matter contributions in the inner disc of NGC 628 are limited to roughly 22 percent of the rotational support.
Kinematic separation of planetary nebulae offers a way to isolate age-dependent vertical structures in discs where photometric age dating is difficult.

Load-bearing premise

The measured velocity dispersion of the warm planetary-nebula population can be directly equated to the vertical scale height of the older stellar population to enable the rotation-curve decomposition.

What would settle it

An independent measurement of the vertical scale height of the older stellar population in NGC 628 that differs from the height implied by the warm component's 26.1 km/s dispersion would invalidate the maximal-disc result.

Figures

Figures reproduced from arXiv: 2604.16045 by Aaron J. Romanowsky, Anastasia Ponomareva, Arianna Cortesi, Claudia Pulsoni, Crescenzo Tortora, Johanna Hartke, Kenneth C. Freeman, Konrad Kuijken, Lodovico Coccato, Magda Arnaboldi, Massimo Capaccioli, Michael Merrifield, Nicola R. Napolitano, Nigel Douglas, Ortwin Gerhard, Steven P. Bamford, Surya Aniyan.

**Figure 1.** Figure 1: (a): Left panel. The optical layout of the PN.S showing the [OIII] left and right arms. The pair of diffraction gratings on the right split the light beam into two identical spectrograph arms, each with opposite dispersion directions. (b): Right panel. The optical arrangement of the PN.S showing the Hα arm. The “B/S first surface” label indicates the beam splitter. The optics to the left of the Beam splitt… view at source ↗

**Figure 2.** Figure 2: RMS plot from the PN.S alignspots calibration, before and after the higher order geometric distortion correction. Y0 is the y coordinate of the emission lines in the image of the pin-hole mask before and after the correction of the residual distortion by the larger order polynomial. The goal of this step is to correct for the spatial distortions introduced by the spectrograph [PITH_FULL_IMAGE:figures/full… view at source ↗

**Figure 3.** Figure 3: Velocity error (pairs) as function of m5007. Distribution of the ∆V/1.805 pairs vs m5007. The normalization factor is introduced as the two sets of images are registered to the same high image quality reference frame. Photometric error and completeness – The PN.S magnitudes were extensively compared with literature values for PN samples in ETGs and disc galaxies, e.g. M31, giving typical dispersion values … view at source ↗

**Figure 4.** Figure 4: a): Left [OIII] 5007Å Narrow band dispersed image from the PN.S L arm. Monochromatic spatially unresolved sources with red circles are bright HII regions, with no continuum; East is down and North is on the Left. b): Hα image of NGC 628 registered on the [OIII] PN.S L image. This narrow band Hα image is not continuum subtracted. 3.1 Selection criteria with PN.S CDI+Hα for PNe in NGC 628 – Matching sources … view at source ↗

**Figure 5.** Figure 5: Color magnitude diagram ([OIII] - Hα) vs. m5007 for the spatially unresolved [OIII] 5007 Å monochromatic emitters with no continuum selected in the PN.S L/R arms + Hα. The continuum black line shows the selection criteria for the PNs vs HII regions expressed in Eq. 5. In the lower panel, blue open circles indicate the PN.S [OIII] emitters matched with the Herrmann & Ciardullo (2009a) spectroscopically conf… view at source ↗

**Figure 6.** Figure 6: Number counts of the isolated, spatially unresolved HII regions (dashed line histograms) and the PNLF of the selected PN candidates in three radial bins (continuum line histogram). The radial bins are selected to have more than 130 PN candidates in each bin to allow for the kinematical decomposition of the LOSVD as in Aniyan et al. (2018). The PNLFs show a clear rise at m5007 ~ 25 and a decline at magnitud… view at source ↗

**Figure 8.** Figure 8: (a) Left panel - Histogram of the line-of-sight velocity VLOS distribution for the color selected PN candidates (green) and the HII regions (red), plus the PN candidates with no Hα emission (cyan). The VLOS is measured from the Doppler shifted λobs with heliocentric correction measured from the unresolved [OIII] emission with the PN.S. (b) Right panel - Histogram of the difference between the line-of-sight… view at source ↗

**Figure 9.** Figure 9: Histogram of the ∆V = VLOS – VHI for the PN candidates in the inner radial bin 80” < r < 170” (red continuum line) vs the histogram of the ∆V of the SNR candidates (black continuum line). There are only 5 matched SNRs in this radial interval. The majority of SNRs in Kreckel et al [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗

**Figure 10.** Figure 10: ∆VPN,LOS vs. m5007 distribution for the PN candidates in NGC 628, with detected Hα emission (green full dots) and those with none (red full squares). The two samples have the same kinematics as function of m5007, see [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗

**Figure 12.** Figure 12: Total RMS homogenized ∆VPN,LOS of the PN populations in NGC 628 covering the radial range 80” ≤ R ≤ 425”. This histogram is obtained by merging the PN’s ∆V distributions in the three radial bins, once they have been rescaled by the RMS value of ∆VPN,LOS distribution in the middle radial bin. The RMS of the combined Gaussian is 21.15 kms–1 , and kurtosis = 0.783. The GMM applied to the merged histogram giv… view at source ↗

**Figure 13.** Figure 13: Cumulative PNLFs in the three radial bins for the kinematically cold-blue lines- and warmer-red lines- PN population. Kinematical decomposition from Aniyan et al. (2018). The green line depicts the cumulative PNLF in the outermost radial bin where the GMM decomposition assigns all PNe to the cold component in the outermost radial bin (240’’ - 425” ). In the inner bin (radial range: 80” – 170”, light conti… view at source ↗

read the original abstract

We present the results for the galaxy NGC 628 observed with the Planetary Nebulae Spectrograph (PN.S) equipped with the H$\alpha$ arm. With the third PN.S arm, the H$\alpha$ arm, we measure the H$\alpha$ fluxes, in addition to fluxes and line-of-sight velocities (LOSV) of monochromatic spatially unresolved [OIII] 5007{\AA} sources. The narrow band color ([OIII] 5007{\AA}-H$\alpha$) vs m5007 magnitude diagram separates planetary nebulae (PNe) from single compact ionized HII regions and supernovae remnants (SNRs), which also emit in [OIII]5007 {\AA}. The goals are to detect bona-fide PNe in the face-on spiral galaxy NGC 628 (M74) so that we can measure the velocity dispersion of the stars perpendicular to the main plane of the disc. This study validates the empirical selection criteria for PNe with the PN.S in star forming discs. We classified 442 PNe and 251 spatially isolated, unresolved HII regions: the PN.S with the H$\alpha$ arm increased the number of known PNe by a factor 4. We find evidence for two kinematically distinct PN populations in the NGC 628 disc. The kinematically cold PN population dominates the PN luminosity function close to the bright cut-off magnitude, indicating that the PN massive, short-lived progenitors dominate the PNLF bright cut-off in NGC 628. The warmer PN component increasingly dominates at fainter magnitudes. The velocity dispersion orthogonal to the disc plane are {\sigma}z,cold = 8.8 kms-1 and {\sigma}z,warm =26.1 kms-1 respectively, over a range of radii 80 to 425 arcsec. These components contribute with the ratio 46% (cold) and 54% (warm). Once the velocity dispersion of the old component is matched with the population's scale height, the decomposition of the rotation curve for NGC 628 leads to a maximal disc, with the rotation of the baryonic component accounting for 78% of the total rotational velocity in NGC 628.

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

New PN data in NGC 628 shows two clear kinematic populations and supports a maximal disc, but the 78% baryon fraction rests on an untested mapping from dispersion to scale height.

read the letter

The paper's real contribution is the expanded sample: 442 PNe in NGC 628, four times previous counts, plus clean separation from HII regions and SNRs using the new Hα arm on PN.S. They measure line-of-sight velocities and fluxes, then split the PNe into cold (σz = 8.8 km/s) and warm (σz = 26.1 km/s) groups with a 46/54 split that holds across radii. The cold group dominates the bright end of the PNLF, which fits expectations for younger progenitors. These are concrete, new numbers from a face-on spiral where vertical motions are directly accessible in the LOS velocities. That part of the work stands on its own and will be useful for anyone modeling disc kinematics with PNe as tracers. The maximal-disc claim is the part that needs more scrutiny. They take the warm dispersion, match it to an old stellar scale height, and then decompose the rotation curve to get baryons supplying 78% of Vrot. This step assumes the warm PNe trace a single old population in vertical equilibrium, with no age mixing from PNLF selection and with the same scale height applying radially. The abstract does not spell out the exact vertical equilibrium formula or how they propagate uncertainties from the population split into the final fraction. In a face-on galaxy those links are not automatic, so the 78% number could shift with different assumptions. The data themselves look solid and the two-population finding is new for this galaxy. This is the kind of observational paper that belongs in a journal like MNRAS or A&A. It gives people working on nearby spirals and on PN kinematics something concrete to cite or test against simulations. I would send it to referees rather than desk-reject; the measurements are worth checking even if the dynamical interpretation gets revised.

Referee Report

2 major / 3 minor

Summary. The paper reports PN.S observations of the face-on spiral NGC 628 with a new Hα arm, yielding 442 classified PNe and 251 isolated HII regions (increasing the known PN sample by a factor of four). It identifies two kinematically distinct PN populations with line-of-sight velocity dispersions σz,cold = 8.8 km s^{-1} (46% contribution) and σz,warm = 26.1 km s^{-1} (54% contribution) over 80–425 arcsec. The warm component is associated with the old stellar population; matching its dispersion to the population scale height allows a standard rotation-curve decomposition in which the baryonic disc supplies 78% of the total circular velocity, implying a maximal disc.

Significance. If the population-to-scale-height mapping and vertical equilibrium assumptions hold, the work supplies new, direct kinematic constraints on the stellar contribution to the rotation curve of a face-on spiral, strengthening the case for maximal discs and limiting the inner dark-matter density. The fourfold increase in PN detections and the empirical validation of the narrow-band selection criteria against HII regions and SNRs constitute a clear observational advance for future PN.S studies of star-forming discs.

major comments (2)

[rotation-curve decomposition paragraph] The rotation-curve decomposition (abstract and the paragraph linking σz,warm to scale height) states that the warm PN dispersion is 'matched with the population's scale height' to obtain the stellar surface-density profile, yet no explicit vertical-equilibrium equation (e.g., σz² = π G Σ⋆ hz for an isothermal sheet or the sech² equivalent), no adopted value of the constant, and no radial dependence of hz are provided. This step is load-bearing for the quoted 78% baryonic fraction.
[kinematic population analysis] The claim that the two kinematic PN populations map one-to-one onto distinct age cohorts (cold = young, warm = old) with different vertical structures is used to assign the 54% warm fraction to the old thin-disc tracer. The text shows velocity histograms but does not quantify possible PNLF selection biases or progenitor-age mixing that could alter the warm-component weight entering the decomposition.

minor comments (3)

[abstract] Abstract: LaTeX artifacts remain (e.g., 'kms-1', 'H$alpha$', '5007{AA}'); units should be written consistently as km s^{-1}.
[classification section] The narrow-band colour–magnitude diagram used to separate PNe from HII regions and SNRs should display the adopted selection boundaries and any magnitude-dependent completeness corrections.
[results] Error propagation on the 78% baryonic fraction, including uncertainties in the adopted scale height and in the exclusion of the cold component, is not reported.

Simulated Author's Rebuttal

2 responses · 0 unresolved

Thank you for the referee's insightful comments on our paper. We address the major comments below and have revised the manuscript accordingly to improve clarity and completeness.

read point-by-point responses

Referee: [rotation-curve decomposition paragraph] The rotation-curve decomposition (abstract and the paragraph linking σz,warm to scale height) states that the warm PN dispersion is 'matched with the population's scale height' to obtain the stellar surface-density profile, yet no explicit vertical-equilibrium equation (e.g., σz² = π G Σ⋆ hz for an isothermal sheet or the sech² equivalent), no adopted value of the constant, and no radial dependence of hz are provided. This step is load-bearing for the quoted 78% baryonic fraction.

Authors: We agree with the referee that the vertical equilibrium relation should be made explicit for clarity and reproducibility. In the revised manuscript, we will state the adopted equation (the isothermal sheet approximation σ_z² = π G Σ_* h_z) and specify that h_z is assumed constant with radius, as is conventional for thin-disc populations in such analyses. This will allow readers to follow how the stellar surface density is derived from the warm component dispersion and how it contributes to the 78% baryonic fraction in the rotation curve decomposition. revision: yes
Referee: [kinematic population analysis] The claim that the two kinematic PN populations map one-to-one onto distinct age cohorts (cold = young, warm = old) with different vertical structures is used to assign the 54% warm fraction to the old thin-disc tracer. The text shows velocity histograms but does not quantify possible PNLF selection biases or progenitor-age mixing that could alter the warm-component weight entering the decomposition.

Authors: The velocity histograms presented in the paper exhibit a clear bimodal distribution, justifying the separation into cold and warm populations, with the cold one dominating the bright PNLF as expected for younger, more massive progenitors. We recognize that quantifying selection biases in the PNLF and potential age mixing would provide additional robustness. We will add a discussion section addressing these issues qualitatively, based on the observed trends and the validation of our selection criteria, while noting that a detailed statistical modeling lies outside the scope of this work. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation uses independent kinematic measurements and standard methods

full rationale

The paper reports direct observational results: classification of 442 PNe and 251 HII regions, identification of two kinematically distinct populations with measured σz,cold = 8.8 km/s and σz,warm = 26.1 km/s, and their radial contributions. The rotation-curve decomposition applies these σz values to match an old stellar population scale height and then uses standard baryonic mass modeling to conclude that baryons supply 78% of V_rot. No equation in the provided chain reduces a claimed prediction or result to a fitted parameter or definition taken from the same dataset by construction. No load-bearing self-citation chain or ansatz smuggling is evident that would make the maximal-disc conclusion tautological. The vertical equilibrium step is a standard external relation applied to new data, not derived internally.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis rests on standard domain assumptions about planetary nebulae as kinematic tracers of the stellar disc and on the validity of matching observed velocity dispersions to scale heights for rotation-curve decomposition; no new free parameters or invented entities are introduced beyond conventional astronomical modeling choices.

axioms (2)

domain assumption Planetary nebulae trace the kinematics of the underlying stellar population
Invoked when interpreting the two PN components as cold and warm stellar populations with different scale heights.
standard math The distance and inclination of NGC 628 are known sufficiently well to convert angular scales to physical radii and velocities
Required for the reported radial range (80–425 arcsec) and velocity dispersions in km/s.

pith-pipeline@v0.9.0 · 5830 in / 1556 out tokens · 48889 ms · 2026-05-10T08:08:53.144816+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

3 extracted references

[1]

C., Gerhard, O

Aniyan, S., Freeman, K. C., Gerhard, O. E., Arnaboldi, M., & Flynn, C. 2016, MNRAS, 456, 1484 Aniyan, S., Freeman, K. C., Arnaboldi, M., et al. 2018, MNRAS, 476, 1909 Aniyan, S., Ponomareva, A. A., Freeman, K. C., et al. 2021, MNRAS, 500, 3579 Arnaboldi, M., Aguerri, J. A. L., Napolitano, N. R., et al. 2002, AJ, 123, 760 Arnaboldi, M., Freeman, K. C., Oka...

2016
[2]

the [OIII] – Hα color ID RA (J2000) DEC (J2000) ∆V kms–1 VLOS kms–1 m5007 – 24.66 mag [OIII] – Hα mag HII N628 J013623.9+154513.7 1:36:23.9 15:45:13.7 2.5 664.4 -0.1 2.8 HII N628 J013624.4+154540.1 1:36:24.4 15:45:40.1 2.9 665.2 0.0 2.5 HII N628 J013625.5+154848.9 1:36:25.5 15:48:48.9 -7.5 669.6 0.1 3.9 HII N628 J013625.5+154858.0 1:36:25.5 15:48:58.0 -3....

2007
[3]

2002), hence m*5007 = 25.3 or m0 = 0.6

The distance modulus of NGC 628 is 29.8 and the bright end of the PNLF is nominally at an absolute 24 magnitude M*5007 = –4.54±0.05 (Ciardullo et al. 2002), hence m*5007 = 25.3 or m0 = 0.6. 25 26 For the Hα arm we carry out independent calibration using a similar set of equations. From the measurements of the 27 spectrophotometric standard stars observed ...

2002