Pith. sign in

REVIEW 2 major objections 4 minor 2 cited by

ODIN: Star Formation Histories Reveal Formative Starbursts Experienced by Lyman Alpha Emitting Galaxies at Cosmic Noon

T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Most Lyman-alpha galaxies are experiencing their largest starburst yet, but a third have star-forming pasts.

desk verdict A careful first non-parametric SFH census of cosmic-noon LAEs, with a useful taxonomy whose exact fractions are softer than the headline implies. read the letter →

arxiv 2501.08568 v3 pith:7U474LE3 submitted 2025-01-15 astro-ph.GA

Nicole M. Firestone , Eric Gawiser , Kartheik G. Iyer , Kyoung-Soo Lee , Vandana Ramakrishnan , Francisco Valdes , Changbom Park , Yujin Yang
show 22 more authors
Anahita Alavi Robin Ciardullo Norman Grogin Caryl Gronwall Lucia Guaita Sungryong Hong Ho Seong Hwang Sang Hyeok Im Woong-Seob Jeong Seongjae Kim Anton M. Koekemoer Ankit Kumar Jaehyun Lee Vihang Mehta Gautam Nagaraj Julie Nantais Laura Prichard Marc Rafelski Hyunmi Song Ben Sunnquist Harry I. Teplitz Xin Wang
This is my paper · ORCID
classification astro-ph.GA
keywords LymanalphaemittersstarformationhistoriesDenseBasisnon-parametricSEDfittingcosmicnoongalaxyevolutionstarburstsODINsurvey
verification ladder T0 review T1 audit T2 compute T3 formal

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 tests a long-standing assumption about Lyman Alpha Emitting galaxies (LAEs): that they are observed during their very first major burst of star formation. Using the Dense Basis method to reconstruct non-parametric star formation histories from rest-UV-through-NIR photometry of 74 ODIN LAEs at redshifts 2.4, 3.1, and 4.5, the authors find that 67% of LAEs match the first-burst archetype, 28% have had significant past bursts with the current burst still largest, and 5% had their largest burst in the past. Combining the first two classes, roughly 95% of LAEs are experiencing their largest burst yet, which they call a formative burst. The authors also find that LAEs built a larger fraction of their stellar mass in the last 200 Myr than mass-matched Lyman Break Galaxies. If correct, strong Ly-alpha emission does not require a unique, first-time starburst, and multiple evolutionary paths can produce LAEs.

What carries the argument

The load-bearing tool is Dense Basis, a non-parametric, Gaussian-process-based method for reconstructing star formation histories from spectral energy distributions. Unlike parametric SFH models that assume a single functional form and tend to place the only peak at the time of observation, Dense Basis fits flexible histories with multiple possible bursts, allowing the authors to identify where the maximum star formation rate falls relative to the observation epoch and to separate first, dominant, and non-dominant burst classes. The classification itself uses a 200 Myr resolution at the time of observation and a 1 solar mass per year threshold for what counts as a significant past burst.

What would settle it

Re-run the same 74 LAEs through an independent non-parametric SFH code or vary the Dirichlet prior parameter alpha and the SFR prior, then recompute the first/dominant/non-dominant fractions; if the ~95% formative-burst result shifts by more than the quoted statistical uncertainties, the classification is a prior artifact rather than a property of the galaxies.

Watch

Extended reading notes

Core claim

The paper's central claim is that LAE star formation histories come in three archetypes, not one. From Dense Basis reconstructions, the authors classify 67% of their 74 LAEs as first-burst galaxies, with only modest star formation before observation; 28% as dominant-burst galaxies, where an earlier significant burst occurred but the current burst has the highest star formation rate; and 5% as non-dominant-burst galaxies, where the peak star formation rate lies in the past. Because the first and dominant classes together make up about 95% of the sample, the authors argue that almost all LAEs are undergoing a formative burst, a burst that plays a major role in building the galaxy. They further report that LAEs have a median fraction of stellar mass formed in the last 200 Myr about 1.3 times higher than mass-matched LBGs, and that most LBG controls show dominant-burst histories. The conclusion is that being a strong Ly-alpha emitter is not evidence of being a first-time starburst, and galaxies can become LAEs through several mass-assembly pathways.

Load-bearing premise

The classification assumes that the Dense Basis reconstruction, with its chosen priors, can reliably tell whether a galaxy's past star formation rate ever rose above about one solar mass per year.

Editorial extensions

If this is right

  • A significant minority of LAEs (about one-third) have had substantial past star formation, so first-burst is not a necessary condition for strong Ly-alpha emission.
  • About 95% of LAEs are experiencing their largest burst of star formation to date, meaning a formative burst is the norm even when the history is not a first burst.
  • LAEs assemble a larger fraction of their stellar mass in the most recent 200 Myr than mass-matched LBGs, confirming that LAEs have less evolved stellar populations.
  • A majority of LBG controls show dominant-burst rather than first-burst histories, pointing to a systematic difference in mass assembly between LAEs and other star-forming galaxies.
  • Because roughly a third of LAEs once had SFRs above the current sample minimum, LAE status may be temporary and repeatable within a galaxy's lifetime.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper, if LAEs can be in dominant bursts, then LAE-based estimates of star formation rate density or clustering may need to account for a fraction of galaxies that are not pristine first-time starbursts.
  • The non-dominant LAEs, though rare, offer a direct test of whether Ly-alpha escape depends more on dust and ISM geometry than on current star formation rate.
  • A testable extension would be to compare Dense Basis classifications with rest-frame optical line ratios or resolved stellar population ages for the same galaxies to see if past bursts leave independent signatures.
  • The paper's definition of a 'significant' past burst uses a fixed 1 solar mass per year threshold tied to the sample's own minimum SFR; an external threshold from Ly-alpha luminosity or stellar mass could change the archetype fractions.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

Summary. The manuscript uses the Dense Basis non-parametric SED-fitting method to reconstruct star formation histories for 74 ODIN narrowband-selected LAEs at z=2.4, 3.1, and 4.5 with UVCANDELS rest-UV-through-NIR photometry. The authors define three SFH archetypes: 'first burst' (no significant past SFR), 'dominant burst' (significant past SFR but current SFR largest), and 'non-dominant burst' (largest SFR in the past). They report 67% first-burst, 28% dominant-burst, and 5% non-dominant-burst LAEs, combine the first two into a 'formative burst' label (~95%), and compare the fraction of stellar mass formed in the last 200 Myr with mass-matched LBG control samples, finding LAEs assemble ~1.3 times more of their mass recently. The central claim is that multiple stellar mass assembly histories can produce strong Ly-alpha emission, and that the simple young-first-burst picture of LAEs is incomplete.

Significance. If the archetype fractions are robust, this is a valuable observational result that moves beyond parametric SFH assumptions and quantifies the diversity of LAE mass assembly at cosmic noon. The comparison with LBG controls, including KS tests, is a concrete and falsifiable check on the claim that LAEs differ from other star-forming galaxies. The paper also contains useful self-checks, such as examining the magnitude and EW dependence of the chi-squared cut and explicitly reporting that including chi-squared-rejected objects shifts the archetype distribution by about 10 percentage points while preserving the ~95% formative-burst fraction. These strengths make the paper a serious candidate for publication, but the missing prior-sensitivity analysis and lack of propagated uncertainties on the archetype fractions currently leave the headline percentages insufficiently protected.

major comments (2)
  1. [Section 4.1 (with §3.1)] The boundary between 'first' and 'dominant' bursts is defined by a threshold of SFR ≥ 1 M☉/yr (§4.1), which is set equal to the sample's minimum time-of-observation SFR. The classification is based on the median reconstructed SFH, yet the paper does not test how the 67/28/5 archetype fractions respond to the adopted Dense Basis priors. The only sensitivity test reported in §3.1 varies the number of t_x time bins; the Dirichlet concentration α=3.0 and the logSFR prior range [−1, 2] are held fixed. Because the reconstruction scatter quoted from Iyer et al. (2019) is 0.2 dex, individual objects can plausibly cross the 1 M☉/yr boundary if α or the SFR bounds change. I request a sensitivity analysis varying α (e.g., 1, 2, 4) and the SFR prior bounds, with the resulting first/dominant/non-dominant fractions and the combined ~95% formative-burst fraction reported for each setup.
  2. [Section 4.1] The archetype percentages (77%, 57%, 59%, etc.) are quoted as point values with no propagated uncertainties, even though each object's SFH is a posterior distribution and the 25–75% intervals are shown in Figures 1–4. Since the three archetypes are mutually exclusive, Poisson binomial counting uncertainties are straightforward to compute. A more complete treatment would classify each object by the fraction of posterior SFH samples satisfying each archetype condition (e.g., max(past SFR) ≥ 1 M☉/yr and max(past SFR) < current SFR) and propagate those per-object probabilities into the reported fractions. Without this, the reader cannot tell whether the 28% dominant-burst fraction is a robust population signal or a threshold artifact of the median SFH.
minor comments (4)
  1. [Abstract and Section 5 (conclusion 4)] The 'formative burst' fraction is written as '~95' and '∼95' without the percent sign; it should read '~95%' for clarity.
  2. [Section 3.1] In the sentence describing priors, 'the log of the instantaneous star formation rateSF R' has a missing space and the formatting of 'SF R' should be corrected to 'SFR'.
  3. [Section 4.2] The notation 'M ed(M∗‘NOW’/M∗)' is awkward and contains typographic quotes; define a clean symbol such as M∗,now/M∗ and use it consistently throughout the text and figures.
  4. [Section 3.3] The choice of a 200 Myr resolution for identifying the SFR peak is attributed to the near-UV SFR timescale of Broussard et al. (2019), but the connection between that timescale and the peak-search bin width could be stated more explicitly.

Circularity Check

1 steps flagged · score 1.0 of 10

No significant circularity: the 'formative burst' label is a definitional union of the measured first/dominant archetypes, so the ~95% figure is a summary of the classification rather than an independent prediction.

  1. self definitional [Section 4.2, 'Formative Bursts of Star Formation']
    "We designate this significant subset of LAEs (∼95) by combining the galaxies with first burst and dominant burst SFH archetypes."

    The 'formative burst' category is defined as the union of the first-burst and dominant-burst archetypes, so the statement that ~95% of LAEs are experiencing a formative burst is true by construction given the reported 67% + 28% archetype fractions. It is a labeling of the already-measured classification rather than an independent derivation. The empirical content resides in the archetype fractions themselves, which are direct readouts of the Dense Basis SFH reconstructions; the 95% figure is therefore a definitional summary, not a separate prediction. This does not undermine the paper's main measurement of the archetype distribution.

full rationale

The paper's central measurement is the archetype distribution (first/dominant/non-dominant) obtained by applying the Dense Basis SFH reconstruction to ODIN/UVCANDELS photometry. That measurement is not circular: the method is validated on simulations in Iyer et al. (2019), a self-cited but independently testable code-based result, and the priors (Dirichlet alpha=3, logSFR bounds, mass range) are stated assumptions, not fitted outputs of this paper. The LAE/LBG comparison and the M*/M* ratio are similarly direct outputs. The one definitional element is the 'formative burst' label: Section 4.2 defines it as the union of first- and dominant-burst archetypes, so the ~95% figure is an arithmetic restatement of the measured 67%+28% fractions rather than an independent derivation. This is a harmless labeling convention, not a circular prediction, because the empirical content is in the archetype fractions themselves. The SFR >= 1 Msun/yr burst threshold is set with reference to the sample's own minimum time-of-observation SFR and is not varied in a full prior sensitivity analysis; that is a robustness/correctness concern, not a circularity, since the threshold is a stated classification choice rather than a parameter fitted to produce the headline fractions. No other step reduces to its inputs or to a self-citation chain.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The central claim rests on the fidelity of Dense Basis SFH reconstruction and on chosen priors and thresholds. No new physical entities are introduced; the formative-burst category is a definitional classification, not a new physical ingredient.

free parameters (5)
  • Dirichlet prior concentration alpha for SFH time bins = 3.0
    Chosen in Section 3.1 from Leja et al. 2017/2019. Controls smoothness and number of SFH peaks, directly affecting first versus dominant burst classification. No sensitivity analysis for alpha is shown.
  • Significant past burst SFR threshold = 1 Msun/yr
    Defined in Section 4.1 as the minimum time-of-observation SFR of the sample. Separates first-burst from dominant-burst SFHs. The threshold is data-dependent and not varied.
  • Recent stellar mass window = 200 Myr
    Defines the time of observation for formative burst and M*_NOW/M* measurements, based on Broussard et al. 2019. If the timescale is wrong, the archetype fractions and recent-mass ratios shift.
  • Photometric error floor fractional uncertainty U = 0.03 for space-based, 0.10 for ground-based and IRAC
    Added to flux errors in Section 3.1. Inflates uncertainties and can smooth SFH reconstructions, affecting classification.
  • Stellar mass prior bounds = log M*/Msun uniform in [7,12]
    Uniform prior in Section 3.1. The paper notes posteriors do not push against the bounds, but the prior width still shapes the posterior and the derived SFH normalization.
assumptions (5)
  • domain assumption Dense Basis non-parametric SED fitting recovers true SFH shapes, including multiple bursts, with about 0.2 dex scatter out to lookback times of about 5 Gyr.
    Relied on for all archetype classifications. Validated in Iyer et al. 2019 Figures 6 and 7, but not proven for this specific photometry, redshift range, and prior choices.
  • domain assumption FSPS stellar population templates, Chabrier IMF, Calzetti dust law, and assumed IGM absorption correctly model LAE SEDs.
    Adopted in Section 3.1 without alternative template or dust law tests.
  • domain assumption Redshifts are confined to ODIN narrowband FWHM with a flat distribution, and UVCANDELS photometric redshifts and DESI spec-z agreement support this.
    Section 3.1. Incorrect redshifts would distort rest-frame wavelengths and the inferred SFH timescales.
  • domain assumption Ly-alpha emission traces recent star formation and AGN contamination is negligible after the chi-squared cut.
    Section 3.2. The paper excludes poor SED fits partly to remove AGN, but no direct AGN diagnostics are shown for the retained sample.
  • domain assumption The 200 Myr timescale approximates the SFR-sensitive stellar population age for defining the current burst.
    Taken from Broussard et al. 2019 and used to define the time of observation and the recent stellar mass fraction.

how reviews work

0 comments
Cite this review

Pith. "Pith review of ODIN: Star Formation Histories Reveal Formative Starbursts Experienced by Lyman Alpha Emitting Galaxies at Cosmic Noon." pith.science (2026). https://pith.science/paper/7U474LE3

@misc{pith2026250108568,
  author       = {Pith},
  title        = {Pith review of: ODIN: Star Formation Histories Reveal Formative Starbursts Experienced by Lyman Alpha Emitting Galaxies at Cosmic Noon},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7U474LE3}},
  note         = {Machine review of arXiv:2501.08568}
}
abstract

In this work, we test the frequent assumption that Lyman Alpha Emitting galaxies (LAEs) are experiencing their first major burst of star formation at the time of observation. To this end, we identify 74 LAEs from the ODIN Survey with rest-UV-through-NIR photometry from UVCANDELS. For each LAE, we perform non-parametric star formation history (SFH) reconstruction using the Dense Basis Gaussian process-based method of spectral energy distribution fitting. We find that a strong majority (67%) of our LAE SFHs align with the frequently assumed archetype of a first major star formation burst, with at most modest star formation rates (SFRs) in the past. However, the rest of our LAE SFHs have significant amounts of star formation in the past, with 28% exhibiting earlier bursts of star formation with the ongoing burst having the highest SFR (dominant bursts), and the final 5% having experienced their highest SFR in the past (non-dominant bursts). Combining the SFHs indicating first and dominant bursts, ~95% of LAEs are experiencing their largest burst yet -- a formative burst. We also find that the fraction of total stellar mass created in the last 200 Myr is ~1.3 times higher in LAEs than in mass-matched Lyman Break Galaxy (LBG) samples, and that a majority of LBGs are experiencing dominant bursts, reaffirming that LAEs differ from other star forming galaxies. Overall, our results suggest that multiple evolutionary paths can produce galaxies with strong observed Ly$\alpha$ emission.

Figures

Figures reproduced from arXiv: 2501.08568 by the authors.

Figure 1
Figure 1. Examples of three classes of observed LAE star formation histories and their SED fits: first burst, dominant burst, and non-dominant burst. Each LAE used in this figure has a confirmed spectroscopic redshift. In each star formation history (top), the x-axis represents the lookback time in Gyrs and the y-axis represents the star formation rate as a function of time in solar masses per year. The black solid line repre… view at source ↗
Figure 2
Figure 2. Star formation histories of z = 2.4 LAEs. The specifications in this figure match those of the star formation histories in [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Star formation histories of z = 3.1 LAEs. The specifications in this figure match those of the star formation histories in [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Star formation histories of z = 4.5 LAEs. The specifications in this figure match those of the star formation histories in [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]
Figure 5
Figure 5. Figure 5: Total stellar mass histograms for three star formation history archetypes of LAEs and LBGs at each redshift. The columns represent z = 2.4 (left), 3.1 (middle), and 4.5 (right). The rows represent the first burst (upper), dominant burst (middle), and non-dominant burst…
Figure 6
Figure 6. Figure 6: Histograms showing the percentage of stellar mass created at the time of observation for LAEs compared to LBGs. The x-axis represents the percentage of stellar mass created at the time of observation (the last 200 Myr) M∗‘NOW’/M∗. The y-axis represents the normalized f…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

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

  1. Subaru meets JWST: A Direct Measurement of Ly$\boldsymbol{\alpha}$ Escape Fraction at $\boldsymbol{z\simeq6.2}$ with Dual Narrow-Band Imaging

    astro-ph.GA 2026-07 conditional novelty 6.5 of 10

    Completeness-weighted stacking of 56 HAEs at z≃6.2 gives median f_esc^Lyα = 0.106^{+0.066}_{-0.044} with no strong Hα-luminosity dependence and UV-linked galaxy-to-galaxy trends.

  2. Quasar Negative Feedback to Surrounding Galaxies Probed with Ly$\alpha$ Emitters and Continuum-Selected Galaxies

    astro-ph.GA 2025-05 conditional novelty 6.0 of 10

    Within the quasar proximity region at z~2.2, Ly-alpha emitters are about 5 sigma less dense than continuum-selected galaxies, implying quasar photoevaporation suppresses low-mass galaxy formation.

Reference graph

Works this paper leans on

24 extracted references · 22 canonical work pages · cited by 2 Pith papers

  1. [1]

    2012, ApJL, 751, L26, doi: 10.1088/2041-8205/751/2/L26 Arrabal Haro, P., Rodr ´ ıguez Espinosa, J., Mu˜ noz-Tu˜ n´ on, C., et al

    Acquaviva, V., Vargas, C., Gawiser, E., & Guaita, L. 2012, ApJL, 751, L26, doi: 10.1088/2041-8205/751/2/L26 Arrabal Haro, P., Rodr ´ ıguez Espinosa, J., Mu˜ noz-Tu˜ n´ on, C., et al. 2020, MNRAS, 495, 1807, doi: 10.1093/mnras/staa1196 Ashby, M., Willner, S., Fazio, G., et al. 2013, ApJ, 769, 80, doi: 10.1088/0004-637X/769/1/80 Berger, V. W., & Zhou, Y. 20...

    arXiv 2012

  2. [4]

    Subaru/Suprime-CamIA427 426.51 20.5 Taniguchi et al. (2007,

    work page 2007

  3. [5]

    Subaru/Suprime-CamIA464 463.59 21.75 Taniguchi et al. (2007,

    work page 2007

  4. [6]

    Subaru/Suprime-CamIA484 485.0 22.8 Taniguchi et al. (2007,

    work page 2007

  5. [7]

    Subaru/Suprime-CamIA505 506.31 23.05 Taniguchi et al. (2007,

    work page 2007

  6. [8]

    Subaru/Suprime-CamIA527 526.11 24.2 Taniguchi et al. (2007,

    work page 2007

  7. [9]

    Subaru/Suprime-CamIA574 576.48 27.25 Taniguchi et al. (2007,

    work page 2007

  8. [10]

    Subaru/Suprime-CamIA624 623.44 29.95 Taniguchi et al. (2007,

    work page 2007

Show all 24 references
  1. [11]

    Subaru/Suprime-CamIA679 678.1 33.55 Taniguchi et al. (2007,

    2007

  2. [12]

    Subaru/Suprime-CamIA709 707.36 31.6 Taniguchi et al. (2007,

    2007

  3. [13]

    Subaru/Suprime-CamIA738 736.19 32.35 Taniguchi et al. (2007,

    2007

  4. [14]

    Subaru/Suprime-CamIA767 768.53 36.4 Taniguchi et al. (2007,

    2007

  5. [15]

    Subaru/Suprime-CamIA827 824.13 34.05 Taniguchi et al. (2007,

    2007

  6. [17]

    Subaru/Suprime-Camg + 474.21 133.5 Taniguchi et al. (2007,

    2007

  7. [19]

    Subaru/Suprime-Camr + 627.06 139.99 Taniguchi et al. (2007,

    2007

  8. [20]

    Subaru/Suprime-Cami + 795.23 136.88 Taniguchi et al. (2007,

    2007

  9. [21]

    Subaru/Suprime-Camz + 906.84 117.0 Taniguchi et al. (2007,

    2007

  10. [23]

    2011; Koekemoer et al

    HST/WFC3F160W1539.23 287.4 CANDELS (Grogin et al. 2011; Koekemoer et al

    2011

  11. [2007]

    H., Hwang, H

    15 Im, S. H., Hwang, H. S., Park, J., et al. 2024, ApJ, 972, 196, doi: 10.3847/1538-4357/ad67d2 Iyer, K., & Gawiser, E. 2017, ApJ, 838, 127, doi: 10.3847/1538-4357/aa63f0 Iyer, K. G. 2019, PhD thesis, Rutgers The State University of New Jersey, School of Graduate Studies Iyer,...

    1998 doi arXiv

  12. [2011]

    2011; Koekemoer et al

    HST/WFC3F125W1251.62 300.4 CANDELS (Grogin et al. 2011; Koekemoer et al

    2011

  13. [2012]

    Subaru/Suprime-CamV545.65 97.48 Taniguchi et al. (2007,

    2007

  14. [2015]

    Subaru/Suprime-CamN B816 815.01 11.8 Taniguchi et al. (2007,

    2007

  15. [2024]

    Subaru/Suprime-CamN B711 711.98 7.2 Taniguchi et al. (2007,

    2007

  16. [2025]

    Subaru/Suprime-CamB441.95 98.0 Taniguchi et al. (2007,

    2007

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.