Evaluating star formation rates at z = 5

A. U. Kapoor; C. Accard; D. Ismail; F. Renaud; J. Freundlich; J.K. Jang; J. Rhee; K. Kraljic; M. B\'ethermin; S. Han

arxiv: 2601.05916 · v2 · submitted 2026-01-09 · 🌌 astro-ph.GA

Evaluating star formation rates at z = 5

D. Ismail , K. Kraljic , M. B\'ethermin , A. U. Kapoor , F. Renaud , C. Accard , J. Freundlich , Y. Dubois

show 6 more authors

S. Han J.K. Jang S. Jeon T. Kimm J. Rhee S. Yi

This is my paper

Pith reviewed 2026-05-16 15:56 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords star formation ratehigh-redshift galaxiesdust attenuationHalphainfrared luminosity[CII] linehydrodynamical simulationsradiative transfer

0 comments

The pith

Simulations show IR, hybrid and [CII] tracers recover star formation rates more reliably than Halpha at redshift 5.

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

The authors post-process two large hydrodynamical simulations of galaxy formation with radiative transfer to generate synthetic observations of z=5 galaxies. They compare how well four common tracers recover the true star-formation rate: the Halpha nebular line, the total infrared continuum, a hybrid infrared-plus-ultraviolet combination, and the [CII] 158-micron line. Halpha estimates display large scatter driven by dust attenuation that depends on viewing angle, dust-to-metal ratio, and the adopted attenuation curve. In contrast, infrared-based and hybrid methods track the 100-Myr-averaged star-formation rate with smaller scatter when sufficient wavelength coverage is available, while the [CII] line yields a usable but steeper calibration than earlier studies.

Core claim

Using synthetic observables from NewHorizon and NewCluster simulations post-processed with Skirt, the paper demonstrates that Halpha-inferred SFRs, averaged over 10 Myr, remain sensitive to calibration choice and suffer substantial scatter from dust attenuation, viewing angle, and dust-to-metal ratio. IR continuum SFRs trace the intrinsic 100-Myr-averaged SFR when full 8-1000 micron coverage exists but underestimate it under typical limited sampling; hybrid IR+UV SFRs reduce scatter without explicit attenuation corrections. The derived [CII]-SFR relation is steeper than previous calibrations yet shows scatter linked to gas density and metallicity. Overall, IR, hybrid, and [CII] tracers prove

What carries the argument

Synthetic observables generated by radiative-transfer post-processing of hydrodynamical simulations, used to compare recovery of intrinsic SFR by Halpha, IR continuum, hybrid IR+UV, and [CII] 158 micron tracers.

If this is right

Halpha SFRs require steeper attenuation curves and careful calibration to reduce but not eliminate systematic bias.
IR continuum SFRs track 100-Myr averages only when the full 8-1000 micron range is well sampled.
Hybrid IR+UV estimates reduce scatter relative to pure IR while avoiding explicit dust corrections.
The [CII] 158-micron line follows a steeper SFR relation at z=5 than calibrations derived from lower-redshift samples.

Where Pith is reading between the lines

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

Observers targeting z=5 galaxies with JWST or ALMA should prioritize hybrid or IR tracers whenever wavelength coverage allows.
Much of the observed discrepancy between optical and infrared SFR indicators at high redshift may trace dust geometry rather than differences in the underlying star-formation history.
Varying the dust-to-metal ratio or attenuation law across a wider suite of simulations would test whether the robustness ordering of tracers holds under different assumptions.

Load-bearing premise

The NewHorizon and NewCluster simulations plus Skirt post-processing faithfully reproduce the dust geometry, metallicity, and star-formation histories of real z=5 galaxies.

What would settle it

A large sample of spectroscopically confirmed z=5 galaxies in which Halpha-derived SFRs show markedly larger scatter than IR-derived SFRs after controlling for inclination and dust column would support the claim; comparable scatter across all tracers would falsify it.

read the original abstract

Inferring the star formation rates (SFR) in high redshift galaxies remains challenging, owing to observational limitations or uncertainties in calibration methods that link luminosities to SFRs. We utilize two state-of-the-art hydrodynamical simulations NewHorizon and NewCluster, post-processed with the radiative transfer code Skirt, to investigate the systematic uncertainties and biases in the inferred SFRs for z=5 galaxies; an epoch where galaxies build-up their stellar mass. We create synthetic observables for widely-used tracers: Halpha nebular line, [CII] 158 micron fine-structure line, total infrared (IR) continuum luminosity, and hybrid (IR + UV). We find that Halpha-inferred SFRs, time-averaged over 10 Myr, are sensitive to the choice of calibration and exhibit substantial scatter driven by dust attenuation, viewing angle, and dust-to-metal ratio. Adopting a steeper attenuation curve reduces this scatter significantly but does not fully eliminate systematic uncertainties. IR continuum-based SFRs trace intrinsic SFRs time-averaged over 100 Myr timescales when a well-sampled continuum emission between restframe 8 and 1000 micron is available and underestimate them with typical approaches when IR data are limited. Nevertheless, IR SFRs display a considerable scatter, largely due to UV photon leakage and strong variations in the star formation history. When UV data are available, hybrid (IR + UV) SFRs provide a more robust estimate, reducing scatter compared to IR-based SFRs while avoiding explicit attenuation corrections. Finally, we derive a [CII]-SFR relation finding a steeper relation than previous studies, however with significant scatter linked to gas density and metallicity. Overall, IR-, hybrid-, and [CII]-based tracers remain more robust than Halpha against variations in optical depth.

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

This sim study ranks IR, hybrid, and [CII] SFR tracers more robust than Halpha at z=5, but the ranking rests on how faithfully NewHorizon/NewCluster plus SKIRT capture real dust geometry and SFHs.

read the letter

The main point is that this paper runs synthetic observations from NewHorizon and NewCluster through SKIRT to test SFR tracers at z=5 and concludes that Halpha suffers more from dust attenuation, viewing angle, and calibration choices, while IR, hybrid IR+UV, and [CII] hold up better against optical depth variations. They also report a steeper [CII]-SFR relation than earlier studies, with scatter tied to gas density and metallicity. The direct match to intrinsic simulation SFRs (time-averaged over 10 or 100 Myr) is a clean way to quantify biases without observational circularity. The hybrid tracer looks especially practical because it lowers scatter without separate attenuation fixes. They test different attenuation curves and note how limited IR sampling affects results, which gives concrete guidance. The soft spot is the lack of external checks on whether the simulations reproduce the actual range of dust optical depths, dust-to-metal ratios, and burstiness seen in real z=5 galaxies. If the model dust geometry is too smooth or the SFHs too steady compared to observations, the claimed robustness ranking for IR and hybrid tracers could be an artifact. The abstract flags free parameters like the H-alpha calibration and dust-to-metal ratio, so sensitivity is explored, but anchoring to data is not described. This is useful for observers choosing indicators for JWST or ALMA data on high-redshift galaxies. It supplies numbers on scatter sources that can be folded into error budgets. I would send it for peer review because the controlled setup is reproducible and the question is directly relevant to cosmic star-formation measurements.

Referee Report

3 major / 2 minor

Summary. The manuscript uses the NewHorizon and NewCluster hydrodynamical simulations post-processed with SKIRT to generate synthetic observables for z=5 galaxies. It compares SFRs inferred from Hα (10 Myr average), [CII] 158 μm, total IR continuum, and hybrid (IR+UV) tracers against the simulations' intrinsic SFRs, concluding that IR-, hybrid-, and [CII]-based tracers are more robust than Hα against variations in optical depth, dust attenuation, viewing angle, and dust-to-metal ratio, while also reporting a steeper [CII]-SFR relation than prior work.

Significance. If the simulations accurately reproduce z=5 dust geometries, metallicities, and SFHs, the direct comparison of inferred versus intrinsic SFRs offers a valuable controlled test of tracer reliability. The multi-tracer synthetic approach and explicit treatment of time-averaging windows and viewing-angle effects are strengths that could inform observational strategies for high-redshift star formation.

major comments (3)

[Abstract and Methods] Abstract and Methods: The central claim that IR, hybrid, and [CII] tracers remain more robust than Hα against optical-depth variations is load-bearing on the assumption that NewHorizon/NewCluster + SKIRT produce realistic distributions of τ_dust, dust-to-metal ratios, and SFH burstiness at z=5; no external validation against observed attenuation curves or IR luminosity functions is described, so the reported scatter differences could be artifacts of the chosen dust model.
[Results on Hα] Results on Hα: The statement that a steeper attenuation curve reduces scatter but does not eliminate systematic uncertainties requires explicit quantification of the residual bias (e.g., median offset and its dependence on the specific curve parameters) and sensitivity tests to the dust-to-metal ratio variations.
[IR and hybrid SFRs] IR and hybrid SFRs: The distinction between 'well-sampled continuum emission between restframe 8 and 1000 μm' (which traces 100 Myr averages) and 'typical approaches' (which underestimate) needs concrete specification of the wavelength sampling and filter sets used in the synthetic photometry, as this directly affects the robustness ranking.

minor comments (2)

[Abstract] Abstract: Define the precise time-averaging windows (10 Myr vs. 100 Myr) and how they are extracted from the simulation star-formation histories for each tracer.
[Throughout] Notation: Ensure the dust-to-metal ratio is consistently defined and its range across the simulated sample is reported with a table or histogram.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed report. We address each major comment point by point below, providing the strongest honest defense of the manuscript while agreeing to revisions where the comments identify genuine gaps in quantification or specification.

read point-by-point responses

Referee: [Abstract and Methods] Abstract and Methods: The central claim that IR, hybrid, and [CII] tracers remain more robust than Hα against optical-depth variations is load-bearing on the assumption that NewHorizon/NewCluster + SKIRT produce realistic distributions of τ_dust, dust-to-metal ratios, and SFH burstiness at z=5; no external validation against observed attenuation curves or IR luminosity functions is described, so the reported scatter differences could be artifacts of the chosen dust model.

Authors: We acknowledge that our robustness conclusions are conditional on the fidelity of the dust and SFH modeling in these simulations. NewHorizon and NewCluster have been validated against a range of lower-redshift observables in prior publications, and SKIRT employs standard dust models calibrated to local galaxies; however, we did not include direct z=5 comparisons to observed attenuation curves or IR luminosity functions. In the revised manuscript we will add an explicit limitations subsection (in Methods or Discussion) that states the assumptions, cites the relevant validation papers, and notes that while relative tracer rankings are internally robust, absolute scatter values should be interpreted cautiously pending future observational constraints at z=5. revision: partial
Referee: [Results on Hα] Results on Hα: The statement that a steeper attenuation curve reduces scatter but does not eliminate systematic uncertainties requires explicit quantification of the residual bias (e.g., median offset and its dependence on the specific curve parameters) and sensitivity tests to the dust-to-metal ratio variations.

Authors: The referee correctly identifies that the current text is qualitative on this point. We will revise the Hα results subsection to report the median offset and 16–84 percentile range between inferred and intrinsic SFR for each attenuation curve tested, and we will add a short sensitivity analysis varying the dust-to-metal ratio by ±0.2 dex around the fiducial value. These numbers and the associated figure/table will be included in the revised manuscript. revision: yes
Referee: [IR and hybrid SFRs] IR and hybrid SFRs: The distinction between 'well-sampled continuum emission between restframe 8 and 1000 μm' (which traces 100 Myr averages) and 'typical approaches' (which underestimate) needs concrete specification of the wavelength sampling and filter sets used in the synthetic photometry, as this directly affects the robustness ranking.

Authors: We agree that the Methods section must be more explicit. In the revision we will list the exact rest-frame wavelength grid (or filter transmission curves) used for the well-sampled IR continuum integration and contrast it with the limited filter sets adopted to mimic typical observational coverage. This specification will be added to the synthetic photometry description and will be cross-referenced in the IR/hybrid results discussion. revision: yes

Circularity Check

0 steps flagged

No circularity: simulation outputs compared directly to intrinsic quantities

full rationale

The paper generates synthetic observables from NewHorizon/NewCluster + SKIRT simulations and directly compares inferred SFRs (from Halpha, IR, hybrid, [CII]) against the simulations' known intrinsic SFRs. No equations, calibrations, or self-citations reduce the reported robustness rankings to input parameters by construction. The central claim is an empirical ranking within the simulated dataset rather than a fitted or self-referential derivation. This is a standard, non-circular use of forward modeling.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The analysis rests on the assumption that the chosen hydrodynamical simulations and radiative-transfer code produce realistic dust and line emission at z=5; calibration constants for each tracer are taken from the literature and tested rather than re-derived.

free parameters (2)

H-alpha calibration factor
Multiple literature calibrations are applied and shown to affect the inferred SFR scatter.
dust-to-metal ratio
Varied in post-processing and shown to drive scatter in H-alpha estimates.

axioms (2)

domain assumption NewHorizon and NewCluster simulations accurately capture z=5 galaxy properties including dust geometry and star-formation histories
Invoked to justify the use of synthetic observables as ground truth.
domain assumption Skirt radiative transfer produces realistic nebular and continuum emission for the simulated galaxies
Required to generate the mock observations used for all comparisons.

pith-pipeline@v0.9.0 · 5681 in / 1430 out tokens · 55733 ms · 2026-05-16T15:56:38.411558+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

IR-, hybrid-, and [CII]-based tracers remain more robust than Halpha against variations in optical depth.

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.

Forward citations

Cited by 2 Pith papers

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

Cosmic evolution of the [CII]-to-molecular gas relation
astro-ph.GA 2026-04 unverdicted novelty 7.0

The [CII]-to-molecular gas conversion factor evolves from over 10,000 to about 10 solar masses per solar luminosity as galaxies enrich from very low metallicity at high redshift to higher values later, making a redshi...
An Updated Characterization of Luminous Ly{\alpha} emitters at the End of Reionization
astro-ph.GA 2026-05 unverdicted novelty 5.0

Luminous Lyα emitters at z≈6 are low-mass ultra-young dwarf starbursts with median Lyα escape fractions above 40 percent, driven by vigorous star formation and low dust content.