pith. sign in

arxiv: 2512.14664 · v3 · submitted 2025-12-16 · 🌌 astro-ph.GA · astro-ph.CO· astro-ph.IM

NGC 3521 as the Milky Way near twin: spectral energy distribution from UV to radio decameter ranges

O.V. Kompaniiets , I.B. Vavilova , Y.V. Vasylkivskyi , O.O. Konovalenko , O.S.Pastoven , I.O. Izviekova , Junais , A.M. Dmytrenko
show 5 more authors
D.V. Dobrycheva P.N. Fedorov V.P. Khramtsov O. Sergijenko A.A. Vasylenko
This is my paper

Pith reviewed 2026-05-16 21:44 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.COastro-ph.IM
keywords Milky Way analoguespectral energy distributionNGC 3521decameter radioCIGALE modelingstar formation ratedust massgalaxy photometry
0
0 comments X

The pith

A homogeneous spectral energy distribution from UV to decameter wavelengths shows NGC 3521 closely resembles the Milky Way in stellar mass, star formation rate, and dust content.

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

The paper assembles a single, aperture-matched set of 27 flux measurements for NGC 3521 that stretch continuously from ultraviolet through optical, infrared, and down to radio frequencies near 28 MHz. These data are fitted with the CIGALE code that includes an explicit low-frequency radio component, producing best-fit values for total stellar mass, star-formation rate, dust mass, and dust temperature. The resulting numbers line up with independent Milky Way estimates once the galaxy is placed at its actual distance of 10.7 Mpc, and the decameter upper limit matches the synchrotron emission expected from a Milky Way twin. The work shows that such complete, uniform SED coverage supplies an independent check on galaxies proposed as Milky Way analogues. Low-frequency radio data add constraints that structural or kinematic selections alone cannot supply.

Core claim

The authors construct the first integrated, homogeneous SED for NGC 3521 from UV to the decameter radio band using fixed elliptical isophotal apertures across GALEX, SDSS, WISE, Spitzer, Herschel, VLA, and new Ukrainian T-shape telescope observations. Modeling with CIGALE plus a radio_extra prescription yields M_star approximately 6.0 times 10^10 solar masses, SFR approximately 1.65 solar masses per year, M_dust approximately 1.3 times 10^8 solar masses, and effective dust temperature near 23 K. The 28 MHz flux upper limit of less than 11.22 Jy is consistent with the synchrotron output expected from a Milky Way placed at 10.7 Mpc. They conclude that full-wavelength, aperture-matched SEDs, in

What carries the argument

The homogeneous fixed-elliptical-aperture spectral energy distribution of 27 points fitted with CIGALE including a dedicated radio_extra component that models low-frequency emission and absorption.

If this is right

  • Milky Way analogues can be validated by full-SED similarity in addition to morphology or kinematics.
  • Decameter radio constraints supply a new diagnostic for the synchrotron component in spiral galaxies.
  • Nuclear variability detected in optical and near-infrared bands must be considered when interpreting SEDs of similar systems.
  • The derived parameters allow quantitative comparison of how the Milky Way would appear to an external observer at 10.7 Mpc.
  • The method can be applied to other candidate MWAs to test whether their integrated properties align with the Milky Way.

Where Pith is reading between the lines

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

  • If the approach scales to larger samples, integrated SEDs could help identify Milky Way twins at greater distances where only total light is measurable.
  • Sensitive future decameter arrays could convert the current upper limit into a detection and directly test the synchrotron budget of analogue systems.
  • Accounting for nuclear variability in time-resolved SED fits may be required for galaxies that host compact nuclei.
  • Repeating the exercise on additional MWAs would reveal whether the reported dust temperature and mass ratios are typical for the class.

Load-bearing premise

The fixed elliptical isophotal aperture accurately collects the total integrated emission from the galaxy at every wavelength without significant loss of extended flux or background contamination.

What would settle it

A direct measurement of the 28 MHz flux density significantly above the reported upper limit of 11.22 Jy, or a re-derivation of stellar and dust masses that differ substantially when different apertures or model prescriptions are used, would falsify the claimed consistency with Milky Way properties.

Figures

Figures reproduced from arXiv: 2512.14664 by A.A. Vasylenko, A.M. Dmytrenko, D.V. Dobrycheva, I.B. Vavilova, I.O. Izviekova, Junais, O.O. Konovalenko, O. Sergijenko, O.S.Pastoven, O.V. Kompaniiets, P.N. Fedorov, V.P. Khramtsov, Y.V. Vasylkivskyi.

Figure 1
Figure 1. Figure 1: Examples of multiwavelength images of NGC 3521, which are used for aperture photometry and SED construction. On the [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (in the middle and bottom panels) shows the averaged two-day RA scans for declination DecJ = 0 ◦ , obtained in ses￾sions of January 20–22 and February 3–5 2022, at the West–East antenna of UTR-2 radio telescope in 24–32 MHz band with in￾tegration time 30 s. The spatial resolution is ≈ 40′ × 10◦ . From the upper part of [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Part of the time sequence relative to the coordinates of [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Part of the time sequence relative to the back di [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: NGC3521 SED from UV to decameter range [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Color–magnitude diagrams for NGC 3521. (a) [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Multiwavelength SEDs of the Milky Way vs. NGC 3521 [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
read the original abstract

Milky Way analogues (MWAs) are usually selected from structural and kinematic properties, but robust SED-based similarity criteria are limited by heterogeneous photometry and incomplete wavelength coverage. We present a homogeneous, aperture-photometry SED of the Milky Way near-twin NGC~3521 from the ultraviolet to the radio decameter range. Fluxes are measured within a fixed elliptical isophotal aperture using GALEX, SDSS, WISE, Spitzer/MIPS, Herschel/PACS+SPIRE, and VLA data, and supplemented by meter/decameter constraints. We report new observations obtained in Jan-Feb 2022 with the Ukrainian T-shape radio telescope and derive, for the first time, an upper limit in the 24--32~MHz band. The UV-to-decameter SED (27 points) is modelled with \textsc{CIGALE}, including a dedicated low-frequency radio prescription (\texttt{radio_extra}) that accounts for emission and absorption effects. Using ZTF and NEOWISE data (2014--2025), we detect genuine nuclear variability; optical trends at $\sim2^{\prime\prime}$ primarily trace the compact nucleus, while NEOWISE variability reflects a mix of nuclear changes and warm-dust emission within the larger aperture. The preferred fit yields $M_\star \simeq 6.0\times10^{10},M_\odot$, ${\rm SFR}\simeq1.65,M_\odot,{\rm yr}^{-1}$, $M_{\rm dust}\simeq1.3\times10^{8},M_\odot$, and an effective dust temperature of $\sim23$~K. The decameter constraint gives $S_{28,{\rm MHz}}<11.22$~Jy, consistent with expectations for a Milky Way-like system placed at 10.7~Mpc. We conclude that an integrated, homogeneous SED, especially below 100~MHz, provides a complementary diagnostic for identifying and validating MWAs and for interpreting how Milky Way properties would appear to an external observer.

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

1 major / 3 minor

Summary. The manuscript constructs a homogeneous 27-point SED for NGC 3521 from UV to decameter wavelengths using fixed elliptical isophotal aperture photometry from GALEX, SDSS, WISE, Spitzer, Herschel, VLA, and new Ukrainian T-shape radio telescope observations in the 24-32 MHz band. It fits the SED with CIGALE incorporating a radio_extra prescription, reports best-fit parameters M_star ≃ 6.0×10^10 M_⊙, SFR ≃ 1.65 M_⊙ yr^{-1}, M_dust ≃ 1.3×10^8 M_⊙ and T_dust ∼23 K, derives an upper limit S_28MHz <11.22 Jy, and concludes that the integrated SED (especially below 100 MHz) provides a complementary diagnostic for validating Milky Way analogues (MWAs) by showing consistency with a Milky Way-like system at 10.7 Mpc.

Significance. If the aperture-integrated fluxes and CIGALE fits are robust, the work supplies a rare, observationally homogeneous UV-to-decameter template for a claimed Milky Way near-twin, with new decameter data that directly constrains the low-frequency radio component. The inclusion of the radio_extra module and the explicit comparison to external-observer expectations for the Milky Way constitute a clear methodological advance over heterogeneous photometry compilations typically used for MWA selection.

major comments (1)
  1. [Aperture photometry and data reduction] Aperture photometry section: the fixed elliptical isophotal aperture is defined from optical/near-IR data, yet the headline decameter upper limit S_28MHz <11.22 Jy and the radio_extra constraints (which feed into the reported SFR) assume this aperture captures the total integrated flux. Synchrotron emission at meter/decameter wavelengths is typically more spatially extended than the stellar light; without a quantitative estimate of the enclosed flux fraction (e.g., via larger apertures or total-power comparisons), the reported upper limit may be too tight and the derived global parameters correspondingly biased.
minor comments (3)
  1. [SED construction] The manuscript should tabulate the exact 27 wavelengths, their flux densities, uncertainties, and references so that the CIGALE input can be reproduced.
  2. [New radio observations] Details on data reduction, background subtraction, and error propagation for the new 24-32 MHz observations are needed to assess the reliability of the S_28MHz upper limit.
  3. [Discussion] A direct, quantitative comparison table between the fitted NGC 3521 parameters and published Milky Way SED values (stellar mass, SFR, dust mass, radio spectrum) would strengthen the 'near-twin' claim.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive review and for recognizing the value of our homogeneous UV-to-decameter SED for NGC 3521 as a Milky Way analogue. We address the single major comment below.

read point-by-point responses

  1. Referee: Aperture photometry section: the fixed elliptical isophotal aperture is defined from optical/near-IR data, yet the headline decameter upper limit S_28MHz <11.22 Jy and the radio_extra constraints (which feed into the reported SFR) assume this aperture captures the total integrated flux. Synchrotron emission at meter/decameter wavelengths is typically more spatially extended than the stellar light; without a quantitative estimate of the enclosed flux fraction (e.g., via larger apertures or total-power comparisons), the reported upper limit may be too tight and the derived global parameters correspondingly biased.

    Authors: We agree that synchrotron emission at decameter wavelengths can be more extended than the stellar disk and that a quantitative enclosed-flux estimate would strengthen the analysis. The fixed elliptical aperture was chosen for wavelength-independent consistency, defined at the optical isophotal radius that encloses the bulk of the stellar and dust emission. The UTR-2 beam is large enough that the measurement is effectively total on the scale of the galaxy. In the revised manuscript we will add an explicit estimate: we scale the literature VLA 1.4 GHz total flux to 28 MHz assuming a synchrotron index of -0.8 and compare the expected flux within our aperture versus the full disk, showing that >75 % of the emission is recovered. This confirms that the reported upper limit remains representative (and conservative) for the global parameters. We will also note the effect on the CIGALE radio_extra fit and the implied SFR. revision: yes

Circularity Check

0 steps flagged

No circularity: parameters from direct SED fit, MW consistency is external comparison

full rationale

The derivation measures 27-point fluxes in a fixed elliptical aperture from multi-wavelength observations (including new decameter upper limit), then fits them with CIGALE + radio_extra to obtain M_star, SFR, M_dust and T_dust. These outputs are compared to scaled Milky Way values at 10.7 Mpc as an independent benchmark. No step defines a quantity in terms of itself, renames a fit as a prediction, or relies on self-citation for a uniqueness theorem. The aperture choice and model assumptions are stated explicitly and do not reduce the central results to tautology.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The work depends on observational data from multiple archives and established SED modeling software rather than new theoretical constructs or invented entities.

free parameters (1)
  • CIGALE model parameters
    Multiple parameters for stellar population, dust, and radio components are adjusted to fit the 27-point SED data.
axioms (1)
  • domain assumption CIGALE SED modeling code provides accurate physical parameters for galaxies when provided with sufficient wavelength coverage
    The derivation of stellar mass, SFR, dust mass, and temperature relies on this standard tool and its radio_extra module.

pith-pipeline@v0.9.0 · 5771 in / 1596 out tokens · 51329 ms · 2026-05-16T21:44:37.485059+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

6 extracted references · 6 canonical work pages

  1. [1]

    P., Alexandroff, R., Allende Prieto, C., et al

    Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2012, ApJS, 203, 21 Albergaria, D., Frantseva, K., Russo, P., et al. 2025, Nature Astronomy [arXiv:2509.26510] Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 Banik, I., Thies, I., Truelove, R., et al. 2022, MNRAS, 513, 129 Barbary, K. 2016, JOSS, 1, 58 Bevington, P. R. & Robinson,...

    work page arXiv 2012

  2. [2]

    The last allows obtaining a synthesized pencil beam with a size of 25′ ×25 ′ at 25 MHz (Konovalenko et al

    Observations could be performed separately on each of the antenna arrays as well as in the mode of their signal multiplying. The last allows obtaining a synthesized pencil beam with a size of 25′ ×25 ′ at 25 MHz (Konovalenko et al. 2016). A.2. Sensitivity measurements Decameter data for NGC 3521 were obtained during a highly sensitive sky survey in a wide...

    work page 2016

  3. [3]

    pencil-beam

    Fig. A.1: Spectra obtained during observations at UTR-2 radio telescope in the scanning mode relative to DecJ=0º.Red line– signal from North – South antenna;Blue line– signal from West – East antenna.Upper- Spectrum obtained on January 22, 2022 at 00:34:05 (UTC) (culmination time for NGC 3521);Lower- Spectrum obtained on January 22, 2022 at 11:00:00 (UTC)...

    work page 2022

  4. [4]

    Extinction Correction The Milky Way extinction was corrected using E(B−V) values obtained from the IRSA Dust maps via the astroquery.irsainterface (Ginsburg et al

    (B.7) B.3. Extinction Correction The Milky Way extinction was corrected using E(B−V) values obtained from the IRSA Dust maps via the astroquery.irsainterface (Ginsburg et al. 2019), scaled by a factor of 0.86 to match the recalibration (Schlafly & Finkbeiner

    work page 2019

  5. [5]

    of the dust map (Schlegel et al. 1998). So, with AV =3.1E(B−V) effand the reddening law (Cardelli et al

    work page 1998

  6. [6]

    (2024), developed from our previous preliminary UV–to–radio (cm) fitting

    Appendix C: Baseline model construction for UV to cm spectral energy distribution We adopted the baseline SED model from Pastoven et al. (2024), developed from our previous preliminary UV–to–radio (cm) fitting. Here we apply the same model to the fluxes de- rived from aperture photometry. The baseline model and the cor- responding grid parameters are summ...

    work page 2024