Pith. sign in

REVIEW 3 major objections 5 minor 2 cited by

A systematic flux excess at 1.5–2.5 microns, strongest in young (≤6 Myr) low-mass (≤3000 solar masses) clusters, is not captured by current stellar population models; stochastic IMF sampling plus pre-main-sequence stars shrinks but does not

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Observed 1.5-2.5 micron excess in embedded young star clusters in the FEAST galaxies that CIGALE stellar population models fail to fit, strongest for the youngest and lowest-mass clusters.

T0 review reviewed 2026-08-05 challenge →

load-bearing objection Solid multi-galaxy measurement of a real NIR excess in young embedded clusters; the 'missing ingredients' interpretation is plausible but the nebular grid coverage caveat is not fully closed. the 3 major comments →

arxiv 2509.01670 v1 pith:5O3FOLJI submitted 2025-09-01 astro-ph.GA

The near infrared SED of young star clusters in the FEAST galaxies: Missing ingredients at 1-5 $\mu$m

classification astro-ph.GA
keywords Star clustersSpectral energy distributionNear-infrared excessEmbedded young star clustersPre-main-sequence starsStochastic IMF samplingPolycyclic aromatic hydrocarbonsJWST NIRCam
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The reading

This paper uses combined HST and JWST photometry from 0.2 to 5 microns to fit the spectral energy distributions of about 3,800 emerging young star clusters in four nearby galaxies and claims that current stellar population models systematically underestimate their 1.5–2.5 micron flux. The excess is largest in the youngest (≤6 Myr) and least massive (≤3,000 solar masses) clusters, appears in the F150W and F200W (and F277W) bands, and survives tests with different apertures, different dust grids, and fits with dust emission removed. Stochastic sampling of the initial mass function with pre-main-sequence stars in the slug code moves model colors toward the observed ones but cannot close the gap in the metal-rich spirals. If the claim holds, SED-fit ages and masses of embedded young clusters are biased, and a fraction of strong Pa-alpha emitters are misclassified as older than 6 Myr. The paper concludes that emission from young stellar objects and more realistic IMF sampling are the missing ingredients.

Core claim

The central discovery is a systematic near-infrared excess: observed fluxes in the 1.5 and 2.0 micron JWST bands exceed the best-fit model fluxes by roughly 0.2–0.6 mag for the young, low-mass subpopulation, while neighboring blue and red bands are fit well. The excess is strongest for clusters fitted with ages ≤3 Myr and masses ≤3,000 solar masses, and in every galaxy it persists when dust emission is removed from the fit and when smaller apertures are used. As a consequence, the fitting code assigns ages ≥6 Myr to a subset of clusters whose Pa-alpha equivalent widths indicate much younger ages, and it prefers dust parameters that are not those expected from mid/far-infrared studies of star

What carries the argument

The load-bearing quantity is the per-filter residual Δm = m_model − m_observed, computed for thousands of eYSCs and summarized as medians in age/mass bins; positive Δm in F150W/F200W defines the NIR excess. The argument is carried by two SED tools: CIGALE, a deterministic fitting code using instantaneous-burst stellar populations, CLOUDY nebular emission, an attenuation law, and dust templates, and slug, a stochastic stellar-population synthesis code that samples the IMF and includes pre-main-sequence stars. The paired Pa-alpha equivalent-width age comparison provides an independent clock that exposes the age bias.

Load-bearing premise

The load-bearing premise is that the grid of ionized-gas models used in the fitting covers all realistically strong nebular continuum at 1.5–2.5 microns; if the true nebular light is stronger than any grid entry, the reported excess is a limitation of the grid, not a missing stellar or dust component.

What would settle it

Take spectra of about ten low-mass (≤3,000 solar masses) clusters with SED-fitted ages ≤3 Myr in M51 over 1.4–2.6 µm with JWST/NIRSpec. If the observed continuum matches a stellar-plus-nebular model within a few percent band by band, the excess is a photometric or modeling artifact; if a broad smooth excess remains above the model and is not Pa-alpha, the missing component is physical.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Ages from CIGALE for young, low-mass eYSCs are biased: a fraction of strong Pa-alpha emitters with equivalent widths pointing to ages well below 6 Myr are fitted with ages ≥6 Myr.
  • Any SED-fitting pipeline that uses deterministic, fully sampled IMFs with no pre-main-sequence stars will underestimate the 1.5–2.5 µm flux of embedded clusters and will push dust parameters to unrealistic values to compensate.
  • The Pa-alpha equivalent-width age relation is unreliable for clusters below roughly 5,000 solar masses, where stochastic sampling makes a given EW consistent with both young and old ages.
  • Adding young stellar object SEDs and stochastically sampled IMFs with pre-main-sequence emission to fitting codes is necessary to recover credible ages, masses, and dust parameters for emerging clusters.

Where Pith is reading between the lines

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

  • If optically selected young clusters already show ~0.2 mag residuals, legacy HST-based cluster catalogs may carry a milder version of the same NIR bias, systematically underestimating the NIR flux of their youngest objects.
  • The paper does not test whether the excess tracks accretion; correlating the 1.5–2.5 µm residual with Br-gamma or Pa-alpha line width within the ≤3 Myr bin would separate a YSO-disk origin from a purely photospheric PMS origin.
  • Adding YSO SED templates to the fitting grid is a direct falsification test: the excess should disappear and the fitted dust parameters should return to the Umin~1–10, gamma~0.5–1 ranges expected from MIR-FIR studies.
  • Because the 3.3 µm PAH continuum subtraction assumes stellar plus nebular plus dust model continua, a missing YSO component would bias measured PAH equivalent widths, with consequences for PAH-based star-formation-rate calibrations in embedded regions.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper combines HST and JWST/NIRCam photometry to fit 0.2–5 μm SEDs of ~3800 emerging young star clusters (eYSCs) in four FEAST galaxies (M51, M83, NGC 628, NGC 4449) with CIGALE. The central empirical claim is a systematic model-observation residual at 1.5–2.5 μm (F150W and F200W), where CIGALE best-fit fluxes underestimate the observed fluxes; the residual is reported for all four galaxies and appears strongest for clusters with CIGALE-fitted ages ≤6 Myr and masses ≤3000 M⊙. The authors test and reject aperture size and dust-model choice as the origin, show that a dust-free fit up to F200W preserves the excess, and use slug stochastic IMF simulations to show that PMS stars plus extinction can partially but not fully reproduce the scatter. They conclude that current SSP-based models are missing one or more NIR ingredients, possibly YSO emission, stochastic IMF/PMS contributions, or non-homogeneous dust geometries. The paper also reports that CIGALE assigns ages ≥6 Myr to a fraction of strong Paα emitters whose Paα equivalent widths suggest younger ages.

Significance. If the central claim holds, this is an important empirical falsification test for current SED models of embedded young clusters: a reproducible 1.5–2.5 μm excess across four galaxies and hundreds to thousands of clusters would affect age/mass derivations and motivate inclusion of YSO SEDs, PMS emission, and stochastic IMF sampling in cluster SED fitting. The paper has clear strengths: the residual is quantified consistently across multiple environments; the aperture test (App. B), the dust-free test (Sec. 5.2), and the slug comparison (Sec. 5.4) are good robustness checks; and the public data products are a useful community resource. However, the headline claim that the excess is 'not accounted for by current models' is stronger than what the tested CIGALE grid demonstrates, because only a restricted nebular parameter cube is explored. In addition, the demographic statement that the excess is strongest in young, low-mass clusters is made using age and mass bins from the very fits that fail in the NIR, which introduces a partial circularity. These issues are addressable, but they are load-bearing for the paper's main conclusions.

major comments (3)
  1. [Sec. 3.1.1 / Tab. 3] The claim that the F150W/F200W excess is 'not accounted for by current stellar population models' is too strong relative to the experiments. The tested nebular grid is log U = -3.5…-2, n_e = 10, 100 cm^-3, f_esc = 0.01…0.6, f_dust = 0.01…0.3. The insensitivity test in Sec. 3.1.1 varies parameters only within these ranges; it does not establish that the ranges bracket the true nebular emission of compact, embedded, very young HII regions. At 1.5–2.5 μm, free-free and free-bound nebular continuum can be significant, and a stronger continuum than any grid point (e.g., higher n_e, very low f_esc, or different covering factor) would produce exactly this residual pattern. Please either extend the grid to physically motivated extremes (e.g., n_e up to 10^4 cm^-3, f_esc below 0.01, larger f_dust), quantify the maximal nebular contribution at F150W/F200W, or temper the Abstract/Sec. 6 wording to
  2. [Sec. 4 / Fig. 6] The statement that the NIR excess is 'most prominent in low-mass (≤3000 M⊙) and young (≤6 Myr) clusters' is based on bins of best-fit age and stellar mass from the same CIGALE fits that fail in the NIR. A missing NIR component can bias both fitted ages and masses, so the demographic conclusion is partially circular. What is robust is the direct model-observation residual itself; the dependence of that residual on fitted parameters is a secondary, model-dependent result. The Paα EW analysis in Sec. 5.6 addresses only the misclassification of the >6 Myr bin and is itself affected by excess continuum in F150W/F200W. To support the headline trend, use an independent age or mass indicator (e.g., Paα EW, HST/optical colors, or slug prior predictions), or explicitly rescale the claim to 'clusters that CIGALE classifies as young and low mass'.
  3. [Sec. 3.1.1 / Sec. 5.2] The choice of the 'adopted' dust grid is justified by lower reduced χ² rather than by physical priors, and the paper interprets the preference for low U_min / high γ as evidence of model shortcomings at 3–5 μm. Because the dust parameters are essentially unconstrained without MIR/FIR data, this preference could be a consequence of the same missing NIR component being absorbed by flexible dust templates rather than a separate physical failure. The dust-free fit in Sec. 5.2 is a good test for the 1.5–2.5 μm excess and supports that result, but it does not establish the separate claim about the dust parameter grid. I recommend presenting the dust-grid preference as a suggestive byproduct of the NIR excess, not as an independent demonstration of missing model components.
minor comments (5)
  1. [Sec. 2.2] Typographical issues to correct: 'FULLBOX 4TIGHT' appears to be a malformed dither pattern name, and 'V arun' in the author list should be 'Varun'. These do not affect the science.
  2. [Sec. 5.6 / Eq. (1)] The Paα EW formula mixes line flux and bandpass notation. Please clarify whether Fλ,line and Fλ,cont are in the same units and state how the F187N filter width (240 Å) is applied. Propagating photometric uncertainties into EW uncertainties would make the age comparison in Fig. 10 more informative.
  3. [Appendix A] The test of nebular-grid insensitivity in Appendix A is performed on optically selected YSCs in NGC 4449 using Hα, not on eYSCs in the NIR bands where the excess is found. This is a useful consistency check, but it is not equivalent to testing grid coverage at 1.5–2.5 μm for embedded clusters. The text should state this limitation explicitly.
  4. [Fig. 5 / Sec. 4] The F814W residual distribution shows a negative tail that is described only as 'the opposite behavior'. A brief quantitative statement (e.g., median Δm and fraction of objects beyond ±0.5 mag) would help the reader assess whether this is a separate model deficiency or a compensating effect of the NIR excess.
  5. [Fig. 9 / Fig. 15] The yggdrasil evolutionary track is shown as a gray line, but the caption does not state all model assumptions (e.g., f_cov = 0.5, E(B-V)=0). Adding these values to the caption would improve reproducibility.

Circularity Check

0 steps flagged

No significant circularity: the NIR excess is a direct data–model residual, and the paper's central claims are supported by independent checks.

full rationale

The central claim is an observational residual: m_CIG − m_obs in F150W/F200W (Sec. 4, Figs. 5–6). This residual is not an input to CIGALE; it is the output of a χ² minimization over the full 0.2–5 μm SED. The conclusion that the excess is 'not accounted for by current stellar population models' is a statement about the tested CIGALE grid, not a quantity defined to equal an input. Binning residuals by CIGALE-recovered age and mass (Fig. 6) introduces a statistical covariance between the residual and fitted parameters, and the paper acknowledges the resulting age biases; however, no equation defines the residual in terms of age/mass, and the independent Paα EW age analysis (Sec. 5.6) provides an external check that the qualitative age trend is not purely an artifact. The nebular-grid coverage (Sec. 3.1.1) is a model-completeness assumption: the paper shows insensitivity within the tested log U, n_e, f_esc, f_dust ranges but cannot prove the true eYSC nebular emission is bracketed. This is a correctness risk, not circularity. Self-citations to FEAST companion papers (Gregg et al. 2024; Pedrini et al. 2024; Knutas et al. 2025; Linden et al. in prep; Adamo et al. in prep) are used for data reduction, catalogs, and grid details that are also specified in this paper; they are not used to force the central conclusion. The slug comparison (Sec. 5.4) is a forward model, not a fit to the residuals. No load-bearing step reduces by construction to its input.

Axiom & Free-Parameter Ledger

2 free parameters · 4 axioms · 0 invented entities

No new physical entities are introduced. The candidate explanations (YSO disks, PMS stars, Be/WR stars, clumpy dust, multiple bursts) are existing astrophysical components invoked from the literature (Sec. 5.4-5.5). The ledger's main debt is the adequacy of the physical model grid rather than any invented ingredient.

free parameters (2)
  • Dust grid ranges (umin, gamma) = umin: 0.3-1.0; gamma: 0.01-0.5 (AD grid)
    Chosen because they minimize reduced chi2 in CIGALE; the authors state explicitly there is no physical motivation for these ranges (Sec. 3.1.1, Fig. 2). The excess persists when dust is excluded (Sec. 5.2), so this parameter choice affects the 'model shortcoming' interpretation but not the existence of the residual.
  • Per-cluster fitted age, stellar mass, E(B-V) = Age 1-10 Myr, M* ~ 10^2-10^4 M_sun, E(B-V) medians ~0.6-0.8 mag (Fig. 3)
    CIGALE outputs used to define the age and mass bins that characterize the excess (Sec. 4). These estimates are themselves contaminated by the excess, adding a circular component to the 'young, low-mass' claim.
axioms (4)
  • domain assumption Each eYSCI is treated as a single stellar population from one instantaneous burst of star formation.
    Sec. 3.1: 'We consider an eYSC SED to be generated by an instantaneous burst of star formation.' If many objects are complexes of sub-clusters or contain age spreads of a few Myr, the residual interpretation changes; the authors acknowledge this in Sec. 5.5.
  • domain assumption CIGALE's Bruzual & Charlot (2003) templates, with no PMS or YSO emission, adequately represent the stellar continuum.
    The excess is defined against these templates. The 'missing ingredients' conclusion depends on these templates being otherwise complete; slug with Dotter (2016) PMS tracks is the only test of PMS contribution (Sec. 5.4).
  • domain assumption The nebular emission grid (CLOUDY, Inoue 2011) with tested logU, ne, f_esc, f_dust ranges brackets the true ionized gas contribution.
    Sec. 3.1.1 asserts results are insensitive to the nebular grid, but only within the tested ranges. If the true nebular continuum at 1.5-2.5 microns lies outside this grid, the residual would be a grid-coverage artifact rather than a missing stellar ingredient.
  • domain assumption Pa-alpha equivalent width derived from F150W-F200W continuum interpolation, interpreted with yggdrasil tracks, is a valid age indicator.
    Used to argue that CIGALE wrongly assigns ages >=6 Myr to a fraction of strong Pa-alpha emitters (Sec. 5.6). The paper notes the EW is itself affected by the NIR excess and that the EW-age relation breaks down below 3000 M_sun (Fig. 11, Sec. 5.6), so this assumption carries its own caveats that the authors explicitly flag.

reviewed 2026-08-05 · how reviews work

0 comments
Cite this review

Pith. "Pith review of The near infrared SED of young star clusters in the FEAST galaxies: Missing ingredients at 1-5 $\mu$m." pith.science (2026). https://pith.science/paper/5O3FOLJI

@misc{pith2026250901670,
  author       = {Pith},
  title        = {Pith review of: The near infrared SED of young star clusters in the FEAST galaxies: Missing ingredients at 1-5 $\mu$m},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5O3FOLJI}},
  note         = {Machine review of arXiv:2509.01670}
}
Share X Bluesky LinkedIn Reddit HN
abstract

We present a combined HST and JWST 0.2 - to - 5 $\mu$m analysis of the spectral energy distributions (SEDs) of emerging young star clusters (eYSCs) in four nearby galaxies from the Feedback in Emerging extrAgalactic Star clusTers (FEAST) survey: M51, M83, NGC 628, and NGC 4449. These clusters, selected for their bright Pa$\alpha$ and 3.3 $\mu$m polycyclic aromatic hydrocarbon (PAH) emission, are still associated to their natal gas cloud and have been largely missed in previous HST optical campaigns. We modeled their SEDs using the CIGALE fitting code and identified: i) a systematic flux excess at 1.5 - 2.5 $\mu$m that is not accounted for by current stellar population models; ii) the preference for a set of dust model parameters that is not aligned with expectations from self-consistent analyses of star-forming regions, suggesting model shortcomings also in the 3 - 5 $\mu$m. The near-infrared (NIR) excess is most prominent in low-mass ($\leq 3000$ M$_\odot$) and young ($\leq 6$ Myr) clusters. Additionally, we see that the SED fitting analysis wrongly assigns ages $\geq 6$ Myr to a fraction of strong Pa$\alpha$ emitters with equivalent widths suggestive of significantly younger ages. A parallel analysis with the slug code suggests that stochastic initial mass function (IMF) sampling of pre-main-sequence stars combined with extinction might partially reduce the gap. We conclude that the inclusion of young stellar object SEDs, along with more realistic sampling of the cluster IMF, might be needed to fully account for the stellar population and dust properties of eYSCs.

Figures

Figures reproduced from arXiv: 2509.01670 by Alex Pedrini, Ana Duarte-Cabral, Angela Adamo, Anne S. M. Buckner, Arjan Bik, Benjamin Gregg, Bruce G. Elmegreen, Daniela Calzetti, Drew Lapeer, Giacomo Bortolini, G\"oran \"Ostlin, Helena Faustino Vieira, Jenna E. Ryon, John S. Gallagher, Kathryn Grasha, Kelsey E. Johnson, Linda J. Smith, Linn Roos, Mark R. Krumholz, Matteo Correnti, Matteo Messa, Michele Cignoni, Monica Tosi, Sean T. Linden, Thomas S. -Y. Lai, Varun Bajaj.

Figure 1
Figure 1. Figure 1: NIRCam detected eYSC populations in M51 (left panel), M83 (top right panel), NGC 628 (middle right panel) and NGC 4449 (bottom right panel). For each galaxy, the RGB color scheme is the following: F335M continuum subtracted map (red), F187N continuum subtracted map (green), F115W map (blue). For NGC 628, we replaced F187N with F405N, due to better visual quality of the map. We highlight the position of the… view at source ↗
Figure 2
Figure 2. Figure 2 [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Normalized distributions of age, stellar mass, and E(B-V) for eYSCs in our sample of galaxies, obtained with CIGALE. We use different colors and line styles to distinguish the distributions of each galaxy. In the middle panel, we drew a dashed line at M∗ = 1000 M⊙ to guide in the plot visualization. In the right panel, we added median values of the E(B-V) distributions as vertical lines. (see Sec. 2.2). If… view at source ↗
Figure 4
Figure 4. Figure 4: Observed SEDs, best fit models and recovered best fluxes for four representative eYSCs in M51. For each panel, the title of the plot describes the fitted properties of the specific eYSC (represented by a unique ID): χ 2 reduced, Age and stellar mass. The observed fluxes in HST and JWST filters are shown with blue squares with error bars, while we show the CIGALE best fitted fluxes convolved to each filter … view at source ↗
Figure 5
Figure 5. Figure 5: Distributions of the residuals mCIG−mobs in eight selected filters (F555W, F814W, F115W, F150W, F187N, F200W, F300M, F444W) for the sample of eYSCI in M51. For each filter, mCIG−mobs represents the difference in magnitude between the best fitted flux from the CIGALE fit and the observed flux, respectively. For each panel, we show the total number of eYSCI, which varies in HST filters due to the additional … view at source ↗
Figure 6
Figure 6. Figure 6: Median values of the residual distributions mCIG−mobs (see [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Median values of the residual distributions mCIG−mobs for eYSCI in M51, as shown in the top left panel of [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Median values of the residual distributions mCIG−mobs for the optically selected YSC population in M51 younger than 10 Myr. The layout of the figure is as presented in [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: NIR color-color diagrams in ABmag for simulated slug cluster libraries with different fixed stellar mass and color￾coded by cluster age. The grey line represents a stellar evolutionary track computed with the yggdrasil code. Each point of the track is color-coded by age as the slug cluster libraries. For an observational comparison, we plotted median values and 16th-84th percentiles for the young cluster p… view at source ↗
Figure 10
Figure 10. Figure 10: Distribution of age differences for eYSCI in M51, M83, NGC 628 and NGC 4449, obtained using two different methods. For each panel, on the x-axis, AgeCIGALE represents the best fit values of age from CIGALE, while AgeEW is estimated from Paα EWs, as described in Sec. 5.6. The unfilled histogram shows the entire population of eYSCI in the four FEAST galaxies, while we use filled histograms to compare the ag… view at source ↗
Figure 11
Figure 11. Figure 11: Paα EW as a function of age for simulated slug cluster libraries with different fixed stellar masses of 500, 1000, 3000 5000 and 10000 M⊙, as defined in Sec. 5.4. Each panel corresponds to a specific stellar mass, as indicated by the label. The darker lines represent the median tracks of the simulated populations. we modeled their SEDs using the CIGALE fitting code. The main findings of this analysis are … view at source ↗
Figure 12
Figure 12. Figure 12: Comparison between the CIGALE-estimated ages of the optical YSC population in NGC 4449 (see Sec. 2.3) and the age estimates from the LEGUS survey (Whitmore et al. 2020). In the left panel, LEGUS ages are derived using only broadband filters, while in the right panel, the Hα recombination line is included. The color bar indicates the Hα morphology (from 1: strong Hα emission, to 4: no Hα), as described in … view at source ↗
Figure 13
Figure 13. Figure 13: Same as the top right panel of [PITH_FULL_IMAGE:figures/full_fig_p022_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Same as [PITH_FULL_IMAGE:figures/full_fig_p023_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Same as [PITH_FULL_IMAGE:figures/full_fig_p023_15.png] view at source ↗

discussion (0)

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

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. The emerging timescale of young star clusters regulated by cluster stellar mass

    astro-ph.GA 2026-03 unverdicted novelty 6.0

    Massive young star clusters clear their natal gas faster than lower-mass clusters, based on HST and JWST imaging of four galaxies.

  2. Active Measurement of Two-Point Correlations

    cs.CV 2026-04 unverdicted novelty 5.0

    A classifier-guided active sampling framework produces unbiased, lower-variance estimates of two-point correlations for rare target sources with far fewer human annotations than standard Monte Carlo sampling.

Reference graph

Works this paper leans on

102 extracted references · 15 canonical work pages · cited by 2 Pith papers · 2 internal anchors

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

  3. [3]

    n7eٲe1jԨ =ztL:5֮][ b w?>>ҥKFkD

    thebibliography [1] 20pt to REFERENCES 6pt =0pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key on reference command E...

  4. [4]

    E., Messa , M., et al

    Adamo , A., Ryon , J. E., Messa , M., et al. 2017, , 841, 131, 10.3847/1538-4357/aa7132

  5. [5]

    Adamo , A., Zeidler , P., Kruijssen , J. M. D., et al. 2020, , 216, 69, 10.1007/s11214-020-00690-x

  6. [6]

    2016, The Journal of Open Source Software, 1, 58, 10.21105/joss.00058

    Barbary , K. 2016, The Journal of Open Source Software, 1, 58, 10.21105/joss.00058

  7. [7]

    2014, , 444, 3829, 10.1093/mnras/stu1695

    Bastian , N., Adamo , A., Schirmer , M., et al. 2014, , 444, 3829, 10.1093/mnras/stu1695

  8. [8]

    A., Pogge , R

    Berg , D. A., Pogge , R. W., Skillman , E. D., et al. 2020, , 893, 96, 10.3847/1538-4357/ab7eab

  9. [9]

    2019, , 622, A103, 10.1051/0004-6361/201834156

    Boquien , M., Burgarella , D., Roehlly , Y., et al. 2019, , 622, A103, 10.1051/0004-6361/201834156

  10. [10]

    A., et al

    Bresolin , F., Kudritzki , R.-P., Urbaneja , M. A., et al. 2016, , 830, 64, 10.3847/0004-637X/830/2/64

  11. [11]

    2003, , 344, 1000, 10.1046/j.1365-8711.2003.06897.x

    Bruzual , G., & Charlot , S. 2003, , 344, 1000, 10.1046/j.1365-8711.2003.06897.x

  12. [12]

    2013, in Secular Evolution of Galaxies, ed

    Calzetti , D. 2013, in Secular Evolution of Galaxies, ed. J. Falc \'o n-Barroso & J. H. Knapen , 419, 10.48550/arXiv.1208.2997

  13. [13]

    C., et al

    Calzetti , D., Armus , L., Bohlin , R. C., et al. 2000, , 533, 682, 10.1086/308692

  14. [14]

    C., et al

    Calzetti , D., Chandar , R., Lee , J. C., et al. 2010, , 719, L158, 10.1088/2041-8205/719/2/L158

  15. [15]

    C., Sabbi , E., et al

    Calzetti , D., Lee , J. C., Sabbi , E., et al. 2015, , 149, 51, 10.1088/0004-6256/149/2/51

  16. [16]

    T., et al

    Calzetti , D., Adamo , A., Linden , S. T., et al. 2024, , 971, 118, 10.3847/1538-4357/ad53c0

  17. [17]

    A., Clayton , G

    Cardelli , J. A., Clayton , G. C., & Mathis , J. S. 1989, , 345, 245, 10.1086/167900

  18. [18]

    2004, , 413, 145, 10.1051/0004-6361:20031454

    Cervi \ n o , M., & Luridiana , V. 2004, , 413, 145, 10.1051/0004-6361:20031454

  19. [19]

    2003, , 115, 763, 10.1086/376392

    Chabrier , G. 2003, , 115, 763, 10.1086/376392

  20. [20]

    C., Calzetti , D., & O'Connell , R

    Chandar , R., Whitmore , B. C., Calzetti , D., & O'Connell , R. 2014, , 787, 17, 10.1088/0004-637X/787/1/17

  21. [21]

    2019, , 876, 62, 10.3847/1538-4357/ab16cf

    Chastenet , J., Sandstrom , K., Chiang , I.-D., et al. 2019, , 876, 62, 10.3847/1538-4357/ab16cf

  22. [22]

    PDRs4All XIII. Empirical prescriptions for the interpretation of JWST imaging observations of star-forming regions

    Chown , R., Okada , Y., Peeters , E., et al. 2024, arXiv e-prints, arXiv:2411.06061, 10.48550/arXiv.2411.06061

  23. [23]

    A., Hadfield , L

    Crowther , P. A., Hadfield , L. J., Clark , J. S., Negueruela , I., & Vacca , W. D. 2006, , 372, 1407, 10.1111/j.1365-2966.2006.10952.x

  24. [24]

    I., Vogl , C., et al

    Cs \"o rnyei , G., Anderson , R. I., Vogl , C., et al. 2023, , 678, A44, 10.1051/0004-6361/202346971

  25. [25]

    L., Fumagalli , M., & Krumholz , M

    da Silva , R. L., Fumagalli , M., & Krumholz , M. 2012, , 745, 145, 10.1088/0004-637X/745/2/145

  26. [26]

    A., Helou , G., Magdis , G

    Dale , D. A., Helou , G., Magdis , G. E., et al. 2014, , 784, 83, 10.1088/0004-637X/784/1/83

  27. [27]

    2024, , 977, 214, 10.3847/1538-4357/ad7a63

    De Marchi , G., Giardino , G., Biazzo , K., et al. 2024, , 977, 214, 10.3847/1538-4357/ad7a63

  28. [28]

    2023, Frontiers in Astronomy and Space Sciences, 10, 1272771, 10.3389/fspas.2023.1272771

    Dobbs , C. 2023, Frontiers in Astronomy and Space Sciences, 10, 1272771, 10.3389/fspas.2023.1272771

  29. [29]

    2016, , 222, 8, 10.3847/0067-0049/222/1/8

    Dotter , A. 2016, , 222, 8, 10.3847/0067-0049/222/1/8

  30. [30]

    T., & Li , A

    Draine , B. T., & Li , A. 2007, , 657, 810, 10.1086/511055

  31. [31]

    T., Dale , D

    Draine , B. T., Dale , D. A., Bendo , G., et al. 2007, , 663, 866, 10.1086/518306

  32. [32]

    T., Aniano , G., Krause , O., et al

    Draine , B. T., Aniano , G., Krause , O., et al. 2014, , 780, 172, 10.1088/0004-637X/780/2/172

  33. [33]

    2001, , 365, 110, 10.1051/0004-6361:20000029

    Eiroa , C., Garz \'o n , F., Alberdi , A., et al. 2001, , 365, 110, 10.1051/0004-6361:20000029

  34. [34]

    G., & Efremov , Y

    Elmegreen , B. G., & Efremov , Y. N. 1996, , 466, 802, 10.1086/177554

  35. [35]

    2023, , 671, L14, 10.1051/0004-6361/202346240

    Fahrion , K., & De Marchi , G. 2023, , 671, L14, 10.1051/0004-6361/202346240

  36. [36]

    2019, , 628, A60, 10.1051/0004-6361/201834403

    Fensch , J., Duc , P.-A., Boquien , M., et al. 2019, , 628, A60, 10.1051/0004-6361/201834403

  37. [37]

    J., Porter , R

    Ferland , G. J., Porter , R. L., van Hoof , P. A. M., et al. 2013, , 49, 137, 10.48550/arXiv.1302.4485

  38. [38]

    Gaia Collaboration , Vallenari , A., Brown , A. G. A., et al. 2023, , 674, A1, 10.1051/0004-6361/202243940

  39. [39]

    2024, , 971, 115, 10.3847/1538-4357/ad54b4

    Gregg , B., Calzetti , D., Adamo , A., et al. 2024, , 971, 115, 10.3847/1538-4357/ad54b4

  40. [40]

    A., Pipher , J

    Gutermuth , R. A., Pipher , J. L., Megeath , S. T., et al. 2011, , 739, 84, 10.1088/0004-637X/739/2/84

  41. [41]

    2019, , 624, A137, 10.1051/0004-6361/201935051

    Haemmerl \'e , L., Eggenberger , P., Ekstr \"o m , S., et al. 2019, , 624, A137, 10.1051/0004-6361/201935051

  42. [42]

    H., Ladd , E

    Heyer , M. H., Ladd , E. F., Myers , P. C., & Campbell , B. 1990, , 99, 1585, 10.1086/115441

  43. [43]

    A., Strom , S

    Hillenbrand , L. A., Strom , S. E., Calvet , N., et al. 1998, , 116, 1816, 10.1086/300536

  44. [44]

    J., & Tielens , A

    Hollenbach , D. J., & Tielens , A. G. G. M. 1999, Reviews of Modern Physics, 71, 173, 10.1103/RevModPhys.71.173

  45. [45]

    Hosokawa , T., Offner , S. S. R., & Krumholz , M. R. 2011, , 738, 140, 10.1088/0004-637X/738/2/140

  46. [46]

    R., Brogan , C

    Hunter , T. R., Brogan , C. L., MacLeod , G., et al. 2017, , 837, L29, 10.3847/2041-8213/aa5d0e

  47. [47]

    A., Johnson , K

    Indebetouw , R., Whitney , B. A., Johnson , K. E., & Wood , K. 2006, , 636, 362, 10.1086/497886

  48. [48]

    Inoue , A. K. 2011, , 415, 2920, 10.1111/j.1365-2966.2011.18906.x

  49. [49]

    S., & Haugb lle , T

    Jensen , S. S., & Haugb lle , T. 2018, , 474, 1176, 10.1093/mnras/stx2844

  50. [50]

    E., Indebetouw , R., & Pisano , D

    Johnson , K. E., Indebetouw , R., & Pisano , D. J. 2003, , 126, 101, 10.1086/375459

  51. [51]

    E., Indebetouw , R., Watson , C., & Kobulnicky , H

    Johnson , K. E., Indebetouw , R., Watson , C., & Kobulnicky , H. A. 2004, , 128, 610, 10.1086/422017

  52. [52]

    E., Leroy , A

    Johnson , K. E., Leroy , A. K., Indebetouw , R., et al. 2015, , 806, 35, 10.1088/0004-637X/806/1/35

  53. [53]

    C., Seth , A

    Johnson , L. C., Seth , A. C., Dalcanton , J. J., et al. 2016, , 827, 33, 10.3847/0004-637X/827/1/33

  54. [54]

    J., & Hartmann , L

    Kenyon , S. J., & Hartmann , L. 1995, , 101, 117, 10.1086/192235

  55. [55]

    Kim , J.-G., Kim , W.-T., & Ostriker , E. C. 2018, , 859, 68, 10.3847/1538-4357/aabe27

  56. [56]

    2025, arXiv e-prints, arXiv:2505.08874, 10.48550/arXiv.2505.08874

    Knutas , A., Adamo , A., Pedrini , A., et al. 2025, arXiv e-prints, arXiv:2505.08874, 10.48550/arXiv.2505.08874

  57. [57]

    Kruijssen , J. M. D. 2015, , 454, 1658, 10.1093/mnras/stv2026

  58. [58]

    R., Fumagalli , M., da Silva , R

    Krumholz , M. R., Fumagalli , M., da Silva , R. L., Rendahl , T., & Parra , J. 2015, , 452, 1447, 10.1093/mnras/stv1374

  59. [59]

    R., McKee , C

    Krumholz , M. R., McKee , C. F., & Bland-Hawthorn , J. 2019, , 57, 227, 10.1146/annurev-astro-091918-104430

  60. [60]

    Lada , C. J. 1987, in IAU Symposium, Vol. 115, Star Forming Regions, ed. M. Peimbert & J. Jugaku , 1

  61. [61]

    D., et al

    Leitherer , C., Schaerer , D., Goldader , J. D., et al. 1999, , 123, 3, 10.1086/313233

  62. [62]

    K., Schinnerer , E., Hughes , A., et al

    Leroy , A. K., Schinnerer , E., Hughes , A., et al. 2021, , 257, 43, 10.3847/1538-4365/ac17f3

  63. [63]

    C., Bolatto , A

    Levy , R. C., Bolatto , A. D., Mayya , D., et al. 2024, , 973, L55, 10.3847/2041-8213/ad7af3

  64. [64]

    T., Lai , T., Evans , A

    Linden , S. T., Lai , T., Evans , A. S., et al. 2024, , 974, L27, 10.3847/2041-8213/ad7eae

  65. [65]

    2009, , 699, 1938, 10.1088/0004-637X/699/2/1938

    Ma \' z Apell \'a niz , J. 2009, , 699, 1938, 10.1088/0004-637X/699/2/1938

  66. [66]

    C., Thilker , D

    Maschmann , D., Lee , J. C., Thilker , D. A., et al. 2024, , 273, 14, 10.3847/1538-4365/ad3cd3

  67. [67]

    I., & Hosokawa , T

    Matsukoba , R., Vorobyov , E. I., & Hosokawa , T. 2024, arXiv e-prints, arXiv:2404.05805, 10.48550/arXiv.2404.05805

  68. [68]

    C., Van Dyk , S

    Mauerhan , J. C., Van Dyk , S. D., & Morris , P. W. 2009, , 121, 591, 10.1086/603544

  69. [69]

    2011, , 142, 40, 10.1088/0004-6256/142/2/40

    ---. 2011, , 142, 40, 10.1088/0004-6256/142/2/40

  70. [70]

    2012, , 538, A110, 10.1051/0004-6361/201117955

    Meilland , A., Millour , F., Kanaan , S., et al. 2012, , 538, A110, 10.1051/0004-6361/201117955

  71. [71]

    H., Federrath , C., & Krumholz , M

    Menon , S. H., Federrath , C., & Krumholz , M. R. 2022, , 517, 1313, 10.1093/mnras/stac2702

  72. [72]

    2018, , 473, 996, 10.1093/mnras/stx2403

    Messa , M., Adamo , A., \"O stlin , G., et al. 2018, , 473, 996, 10.1093/mnras/stx2403

  73. [73]

    2021, , 909, 121, 10.3847/1538-4357/abe0b5

    Messa , M., Calzetti , D., Adamo , A., et al. 2021, , 909, 121, 10.3847/1538-4357/abe0b5

  74. [74]

    Linking stellar populations to HII regions across nearby galaxies. II. Infrared Reprocessed and UV Direct Radiation Pressure in HII Regions

    Pathak , D., Leroy , A., Thompson , T., et al. 2025, arXiv e-prints, arXiv:2502.00165, 10.48550/arXiv.2502.00165

  75. [75]

    Peck , M. L. 1980, , 238, 79, 10.1086/157959

  76. [76]

    2024, , 971, 32, 10.3847/1538-4357/ad534d

    Pedrini , A., Adamo , A., Calzetti , D., et al. 2024, , 971, 32, 10.3847/1538-4357/ad534d

  77. [77]

    2023, arXiv e-prints, arXiv:2310.08720, 10.48550/arXiv.2310.08720

    Peeters , E., Habart , E., Berne , O., et al. 2023, arXiv e-prints, arXiv:2310.08720, 10.48550/arXiv.2310.08720

  78. [78]

    S., Grebel , E

    Pilyugin , L. S., Grebel , E. K., & Zinchenko , I. A. 2015, , 450, 3254, 10.1093/mnras/stv932

  79. [79]

    Richert , A. J. W., Feigelson , E. D., Getman , K. V., & Kuhn , M. A. 2015, , 811, 10, 10.1088/0004-637X/811/1/10

  80. [80]

    Richert , A. J. W., Getman , K. V., Feigelson , E. D., et al. 2018, , 477, 5191, 10.1093/mnras/sty949

Showing first 80 references.

This paper was first reviewed by deepseek-v4-flash on August 5, 2026.