pith. sign in

arxiv: 2604.19551 · v1 · submitted 2026-04-21 · 🌌 astro-ph.SR · astro-ph.GA

Characterizing and spectrally modeling embedded FUor eruptions in the near-infrared

Jiaxun Li , Tinggui Wang , Zheyu Lin This is my paper

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

classification 🌌 astro-ph.SR astro-ph.GA
keywords FUorsepisodic accretionyoung stellar objectsinfrared spectroscopyaccretion disksembedded starsextinctiontime-domain surveys
0
0 comments X

The pith

Near-infrared spectra and disk models confirm infrared-selected young stellar objects as embedded FUor eruptions with high extinction.

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

The paper seeks to demonstrate that episodic accretion outbursts in young stars, known as FUors, can be reliably identified even when deeply embedded in dust by combining infrared light-curve selection with near-infrared spectroscopy and spectral energy distribution modeling. Optical surveys miss these obscured sources, yet they are central to resolving how young stars build most of their mass and to explaining the luminosity problem in star formation. Using Gemini IGRINS high-resolution spectra of four candidates plus photometry, the authors detect CO and H2O absorption at 2.3 microns, place the sources in a diagnostic equivalent-width diagram, and fit a two-component disk model of an actively accreting inner disk plus passively irradiated outer disk. The fits reproduce the data, return extinctions of 10-20 magnitudes, and yield stellar mass times accretion rate values matching classical FUors such as FU Orionis. Of 18 candidates chosen from infrared light curves, six to seven of the eight with spectra show these FUor signatures, indicating the selection technique is efficient.

Core claim

All observed objects display CO and H2O absorption bands at 2.3 microns and occupy positions in the Na plus Ca versus CO equivalent width diagram that confirm their classification as FUors. The best-fitting disk models match both the spectral features and the overall continuum shape, returning high extinctions of A_V = 10-20 mag and M_star times M_dot values comparable to those of FU Orionis. Among the 18 sources initially selected by infrared light curves, 6-7 out of the 8 with available spectra exhibit these FUor characteristics.

What carries the argument

Two-component disk model of an actively accreting inner disk plus passively irradiated outer disk, combined with equivalent-width measurements of CO, Na, and Ca lines.

If this is right

  • High extinctions mean these FUors are invisible at optical wavelengths, so infrared time-domain surveys are required to find them.
  • Accretion rates similar to FU Orionis indicate that the same inner-disk outburst mechanism operates in embedded environments.
  • The high confirmation rate among infrared candidates shows that light-curve selection efficiently recovers FUor events in dense star-forming regions.
  • Episodic accretion confirmed in obscured sources supports its proposed role in solving the luminosity problem for young stellar objects.

Where Pith is reading between the lines

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

  • Repeating the infrared selection plus near-IR follow-up on larger samples could produce a more complete count of FUor frequency across the Galaxy.
  • The same modeling framework could be applied to other classes of eruptive young stellar objects to test whether they share the same inner-disk structure.
  • Future infrared facilities could catch the onset of these outbursts and obtain time-resolved spectra to track changes in the accretion flow.

Load-bearing premise

The observed absorption bands, diagram positions, and disk-model fits are assumed to identify these sources uniquely as FUors rather than other variable young stellar objects.

What would settle it

Detection of an infrared outburst light curve in a source that shows no CO or H2O absorption at 2.3 microns or whose spectrum cannot be fit by the accreting-inner plus irradiated-outer disk model.

Figures

Figures reproduced from arXiv: 2604.19551 by Jiaxun Li, Tinggui Wang, Zheyu Lin.

Figure 1
Figure 1. Figure 1: Near-infrared and mid-infrared light curves of the four FUor candidates. We show the lightcurves from K [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Near-infrared spectra of four FUor can￾didates. The spectra are vertically shifted for clarity and have been convolved to a resolv￾ing power of R = 1200. Prominent atomic ab￾sorption features are marked by dashed verti￾cal lines. We show the Brackett series (red), the CO bandheads ∆ν = 2 (orange), 3 (green) and the broad H2O absorption band (gray horizontal bar). ing a reliable flux calibration. In this ca… view at source ↗
Figure 3
Figure 3. Figure 3: Example of wavelength alignment of observed and template spectra for radial velocity determination. The observed spectrum [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Parameter distributions and corner plot of V10 assuming the far distance. The marginalized one-dimensional distributions are shown along the diagonal, where the dashed blue lines indicate the median and the 16th and 84th percentiles ±1σ. The two-dimensional panels show the joint posterior distributions, with contours enclosing 39%, 68%, 95%, and 99% credible regions. tematic offsets introduced by non-conte… view at source ↗
Figure 5
Figure 5. Figure 5: Comparison between the observed spectra and the corresponding best-fit model spectra for the sample objects. The gray [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Near-infrared H- and K-band spectra of V10 compared with best-fitting disk models. Top panel: H band. Second: K band. The [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Equivalent-width diagram of Na I + Ca I ver￾sus CO for the FUor template and the observed candidates. The dashed line represents the empirical criterion defined by Connelley & Reipurth (2018) to separate FUor-type objects from other young stellar objects. Objects lying below the dashed line show enhanced CO absorption relative to Na I and Ca I and are therefore consistent with FUor-like near-infrared spect… view at source ↗
read the original abstract

Context. Episodic accretion in young stellar objects (YSOs) is thought to play a critical role in addressing the "luminosity problem" associated with star formation. However, optical surveys tend to bias against sources that are heavily obscured. Infrared time-domain surveys, such as unTimely WISE, facilitate the identification of such sources within the dense star formation regions of our Galaxy. Aims. We aim to systematically identify and characterize FUor outbursts in infrared-selected YSOs using high-resolution spectroscopy and detailed disk modeling. Methods. We conducted follow-up high-resolution spectroscopy with Gemini South/IGRINS for four FUor candidates discovered in infrared time-domain surveys. Using a combination of photometric and spectroscopic observations, we constructed spectral energy distributions and fit them with a disk model that incorporates an actively accreting inner disk together with a passively irradiated outer disk. Results. All objects show CO and H$_2$O absorption bands at 2.3$\mu$m, and their positions in the Na + Ca versus CO equivalent width diagram further corroborate their classification as FUors. The best-fitting model spectra closely match both the observed spectral features and the overall continuum, providing additional confirmation of the FUor classification. The best-fit models reveal high extinction values ($A_V$ = 10-20 mag), with $M_*\dot{M}$ comparable to those of classical FUors such as FU Orionis. Among 18 sources initially selected via infrared light curves, $6-$7 out of 8 with available spectra exhibit FUor characteristics, implying a high selection efficiency.

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

3 major / 2 minor

Summary. The manuscript reports follow-up Gemini/IGRINS near-infrared spectroscopy and two-component disk modeling (actively accreting inner disk plus passively irradiated outer disk) for four infrared-selected FUor candidates from time-domain surveys. It finds that all objects exhibit CO and H2O absorption bands at 2.3 μm, occupy the appropriate region of the Na+Ca versus CO equivalent-width diagram, and are well-fit by models yielding A_V = 10–20 mag and M_*Ṁ values comparable to classical FUors such as FU Orionis. The authors conclude that 6–7 of the 8 sources with spectra (from an initial sample of 18) display FUor characteristics, implying high selection efficiency for embedded outbursts.

Significance. If the classifications prove robust, the work would be significant for episodic-accretion studies because it supplies a concrete demonstration that IR light-curve selection can efficiently recover heavily extincted FUors missed by optical surveys. The reported consistency of accretion-rate products with known FUors and the high yield among spectroscopically observed targets would strengthen the case that such events contribute substantially to resolving the luminosity problem in embedded star formation.

major comments (3)
  1. [Abstract and Results] Abstract, Results paragraph: The claim that CO/H2O bands, Na+Ca vs. CO EW diagram placement, and two-component disk fits 'further corroborate' and 'provide additional confirmation' of FUor classification is load-bearing but unsupported by any quantitative comparison to a control sample of non-outbursting Class I/II YSOs at comparable extinction. No equivalent-width statistics or model fits for such a sample are presented, so overlap cannot be excluded.
  2. [Results] Results (modeling description): The statement that 'best-fitting model spectra closely match' the data does not include reported χ² values, posterior widths on A_V and M_*Ṁ, or fits to alternative models (e.g., steady-state or truncated inner-disk geometries). Without these, degeneracies between extinction, accretion rate, and disk geometry remain unquantified.
  3. [Methods] Methods: No details are supplied on IGRINS data reduction, equivalent-width measurement procedures, or error propagation. These omissions directly affect the reliability of the diagnostic diagram and the claimed spectral matches.
minor comments (2)
  1. [Abstract] Abstract: The notation 'M_* Mdot' should be rendered as M_*Ṁ (or M_* dot{M}) for standard clarity.
  2. [Abstract] Abstract: The context sentence on infrared time-domain surveys would benefit from an explicit citation to the unTimely WISE survey.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed comments, which have identified areas where the manuscript can be strengthened. We address each major comment below and indicate the revisions we will implement.

read point-by-point responses

  1. Referee: [Abstract and Results] Abstract, Results paragraph: The claim that CO/H2O bands, Na+Ca vs. CO EW diagram placement, and two-component disk fits 'further corroborate' and 'provide additional confirmation' of FUor classification is load-bearing but unsupported by any quantitative comparison to a control sample of non-outbursting Class I/II YSOs at comparable extinction. No equivalent-width statistics or model fits for such a sample are presented, so overlap cannot be excluded.

    Authors: We acknowledge that a direct quantitative comparison against a control sample of non-outbursting Class I/II YSOs at similar extinctions would provide stronger evidence against potential overlap. The CO/H2O features and Na+Ca vs. CO EW diagram are established diagnostics in the FUor literature, and our infrared-variable sources align with known FUors. In the revised manuscript we will add a short discussion citing equivalent-width statistics from published non-FUor YSO samples and explicitly note the absence of a dedicated control sample in the present work as a limitation. This constitutes a partial revision. revision: partial

  2. Referee: [Results] Results (modeling description): The statement that 'best-fitting model spectra closely match' the data does not include reported χ² values, posterior widths on A_V and M_*Ṁ, or fits to alternative models (e.g., steady-state or truncated inner-disk geometries). Without these, degeneracies between extinction, accretion rate, and disk geometry remain unquantified.

    Authors: We agree that quantitative fit statistics and exploration of alternative models are needed to assess degeneracies. In the revised Results section we will report χ² values for the best-fit models, provide posterior widths or uncertainties on A_V and M_*Ṁ, and include a brief comparison to steady-state and truncated-disk geometries. These additions will be made in the next version of the manuscript. revision: yes

  3. Referee: [Methods] Methods: No details are supplied on IGRINS data reduction, equivalent-width measurement procedures, or error propagation. These omissions directly affect the reliability of the diagnostic diagram and the claimed spectral matches.

    Authors: We apologize for these omissions. The revised Methods section will describe the IGRINS data reduction steps, the exact procedures for equivalent-width measurements (including continuum definition and integration windows), and the error propagation methods used for the diagnostic diagram. These details will be added to ensure full transparency and reproducibility. revision: yes

Circularity Check

0 steps flagged

No significant circularity; classification rests on external spectral standards and new observations

full rationale

The paper reports new Gemini/IGRINS spectra and applies standard two-component disk modeling (active inner + passive outer) to photometric/spectroscopic data. FUor classification is corroborated by observed CO/H2O bands at 2.3 μm and positions in the Na+Ca vs. CO EW diagram, both drawn from prior literature on classical FUors rather than defined within this work. Model fits yield A_V and M_*Ṁ values but are not renamed as 'predictions' of the input data; they are direct fits. No self-citation chain is load-bearing for the central claim, and no equation reduces the result to its own inputs by construction. The derivation is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims rely on standard astrophysical assumptions for YSO classification and disk structure drawn from prior literature, plus fitted physical parameters for the specific sources. No new entities are postulated.

free parameters (2)
  • A_V (extinction)
    Fitted per source in the disk model to match observed SED and spectral continuum.
  • M_* Mdot (accretion rate product)
    Fitted to reproduce observed luminosity and spectral features in the actively accreting inner disk.
axioms (2)
  • domain assumption The system consists of an actively accreting inner disk plus a passively irradiated outer disk.
    Invoked as the basis for constructing and fitting the spectral energy distributions and model spectra.
  • domain assumption CO and H2O absorption at 2.3um combined with Na+Ca vs CO equivalent width diagram position uniquely identifies FUors.
    Used to corroborate classification, based on established criteria from prior YSO studies.

pith-pipeline@v0.9.0 · 5589 in / 1660 out tokens · 72787 ms · 2026-05-10T01:36:46.923360+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

53 extracted references · 53 canonical work pages

  1. [1]

    2014, in IAU Symposium, V ol

    Allard, F. 2014, in IAU Symposium, V ol. 299, Exploring the Formation and Evo- lution of Planetary Systems, ed. M. Booth, B. C. Matthews, & J. R. Graham, 271–272

    work page 2014

  2. [2]

    M., et al

    Audard, M., Ábrahám, P., Dunham, M. M., et al. 2014, in Protostars and Planets VI, ed. H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning, 387–410

    work page 2014

  3. [3]

    Bae, J., Hartmann, L., Zhu, Z., & Nelson, R. P. 2014, ApJ, 795, 61

    work page 2014

  4. [4]

    J., Aboobaker, A

    Bock, J. J., Aboobaker, A. M., Adamo, J., et al. 2026, ApJ, 999, 139

    work page 2026

  5. [5]

    S., Hillenbrand, L

    Carvalho, A. S., Hillenbrand, L. A., Hambsch, F.-J., et al. 2023, ApJ, 953, 86

    work page 2023

  6. [6]

    Connelley, M. S. & Reipurth, B. 2018, ApJ, 861, 145 Contreras Peña, C., Lucas, P. W., Minniti, D., et al. 2017, MNRAS, 465, 3011

    work page 2018

  7. [7]

    M., Evans, II, N

    Dunham, M. M., Evans, II, N. J., Terebey, S., Dullemond, C. P., & Young, C. H. 2010, ApJ, 710, 470

    work page 2010

  8. [8]

    J., Dunham, M

    Evans, II, N. J., Dunham, M. M., Jørgensen, J. K., et al. 2009, ApJS, 181, 321

    work page 2009

  9. [9]

    2013, ApJ, 768, 93

    Findeisen, K., Hillenbrand, L., Ofek, E., et al. 2013, ApJ, 768, 93

    work page 2013

  10. [10]

    J., Hillenbrand, L

    Fischer, W. J., Hillenbrand, L. A., Herczeg, G. J., et al. 2023, in Astronomical Society of the Pacific Conference Series, V ol. 534, Protostars and Planets VII, ed. S. Inutsuka, Y . Aikawa, T. Muto, K. Tomida, & M. Tamura, 355

    work page 2023

  11. [11]

    Fitzpatrick, E. L. 1999, PASP, 111, 63 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2021, A&A, 649, A1 Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1

    work page 1999

  12. [12]

    D., Misselt, K

    Gordon, K. D., Misselt, K. A., Bouwman, J., et al. 2021, ApJ, 916, 33

    work page 2021

  13. [13]

    W., Contreras Peña, C., et al

    Guo, Z., Lucas, P. W., Contreras Peña, C., et al. 2021, MNRAS, 504, 830

    work page 2021

  14. [14]

    2016, ARA&A, 54, 135

    Hartmann, L., Herczeg, G., & Calvet, N. 2016, ARA&A, 54, 135

    work page 2016

  15. [15]

    & Kenyon, S

    Hartmann, L. & Kenyon, S. J. 1996, ARA&A, 34, 207

    work page 1996

  16. [16]

    2011, arXiv e-prints, arXiv:1106.3343

    Hartmann, L., Zhu, Z., & Calvet, N. 2011, arXiv e-prints, arXiv:1106.3343

    work page arXiv 2011

  17. [17]

    A., Carvalho, A., van Roestel, J., & De, K

    Hillenbrand, L. A., Carvalho, A., van Roestel, J., & De, K. 2023, ApJ, 958, L27

    work page 2023

  18. [18]

    A., Carvalho, A

    Hillenbrand, L. A., Carvalho, A. S., Stern, D., et al. 2025, ApJ, 988, 77

    work page 2025

  19. [19]

    A., Contreras Peña, C., Morrell, S., et al

    Hillenbrand, L. A., Contreras Peña, C., Morrell, S., et al. 2018, ApJ, 869, 146

    work page 2018

  20. [20]

    2024, igrins/plp: v3.0.0

    Kaplan, K., Lee, J.-J., Sawczynec, E., & Kim, H.-J. 2024, igrins/plp: v3.0.0

    work page 2024

  21. [21]

    2018, AJ, 156, 84

    Kounkel, M., Covey, K., Suárez, G., et al. 2018, AJ, 156, 84

    work page 2018

  22. [22]

    & Wang, T

    Li, J. & Wang, T. 2024, MNRAS, 532, 2683

    work page 2024

  23. [23]

    J., Johnstone, D., et al

    Liu, H., Herczeg, G. J., Johnstone, D., et al. 2022, ApJ, 936, 152

    work page 2022

  24. [24]

    W., Elias, J., Points, S., et al

    Lucas, P. W., Elias, J., Points, S., et al. 2020, MNRAS, 499, 1805

    work page 2020

  25. [25]

    W., Smith, L

    Lucas, P. W., Smith, L. C., Guo, Z., et al. 2024, MNRAS, 528, 1789

    work page 2024

  26. [26]

    2018, in Society of Photo-Optical Instru- mentation Engineers (SPIE) Conference Series, V ol

    Mace, G., Sokal, K., Lee, J.-J., et al. 2018, in Society of Photo-Optical Instru- mentation Engineers (SPIE) Conference Series, V ol. 10702, Ground-based and Airborne Instrumentation for Astronomy VII, ed. C. J. Evans, L. Simard, & H. Takami, 107020Q

    work page 2018

  27. [27]

    M., et al

    Mainzer, A., Bauer, J., Cutri, R. M., et al. 2014, ApJ, 792, 30

    work page 2014

  28. [28]

    2011, ApJ, 731, 53

    Mainzer, A., Bauer, J., Grav, T., et al. 2011, ApJ, 731, 53

    work page 2011

  29. [29]

    M., Caselden, D., Schlafly, E

    Meisner, A. M., Caselden, D., Schlafly, E. F., & Kiwy, F. 2023, AJ, 165, 36

    work page 2023

  30. [30]

    W., Emerson, J

    Minniti, D., Lucas, P. W., Emerson, J. P., et al. 2010, New A, 15, 433 Morales-Calderón, M., Stauffer, J. R., Hillenbrand, L. A., et al. 2011, ApJ, 733, 50

    work page 2010

  31. [31]

    T., Yuk, I.-S., et al

    Park, C., Jaffe, D. T., Yuk, I.-S., et al. 2014, in Society of Photo-Optical Instru- mentation Engineers (SPIE) Conference Series, V ol. 9147, Ground-based and Airborne Instrumentation for Astronomy V , ed. S. K. Ramsay, I. S. McLean, & H. Takami, 91471D

    work page 2014

  32. [32]

    2021, ApJ, 920, 132 Pérez, S., Hales, A., Liu, H

    Park, W., Lee, J.-E., Contreras Peña, C., et al. 2021, ApJ, 920, 132 Pérez, S., Hales, A., Liu, H. B., et al. 2020, ApJ, 889, 59

    work page 2021

  33. [33]

    1996, MNRAS, 283, 1102

    Persic, M., Salucci, P., & Stel, F. 1996, MNRAS, 283, 1102

    work page 1996

  34. [34]

    J., Menten, K

    Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, ApJ, 783, 130

    work page 2014

  35. [35]

    Reipurth, B., Aspin, C., & Herbig, G. H. 2012, ApJ, 748, L5

    work page 2012

  36. [36]

    Robitaille, T. P. 2017, A&A, 600, A11

    work page 2017

  37. [37]

    P., Whitney, B

    Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 167, 256

    work page 2006

  38. [38]

    Rodriguez, A. C. & Hillenbrand, L. A. 2022, ApJ, 927, 144

    work page 2022

  39. [39]

    & Hillenbrand, L

    Roychowdhury, T. & Hillenbrand, L. A. 2024, Research Notes of the AAS, 8, 171

    work page 2024

  40. [40]

    2010, The Messenger, 141, 24

    Saito, R., Hempel, M., Alonso-García, J., et al. 2010, The Messenger, 141, 24

    work page 2010

  41. [41]

    F., Mace, G

    Sawczynec, E., Kaplan, K. F., Mace, G. N., et al. 2025, PASP, 137, 034505

    work page 2025

  42. [42]

    F., Meisner, A

    Schlafly, E. F., Meisner, A. M., & Green, G. M. 2019, ApJS, 240, 30

    work page 2019

  43. [43]

    Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337

    work page 1973

  44. [44]

    2025, A&A, 695, A130

    Siwak, M., Kóspál, Á., Ábrahám, P., et al. 2025, A&A, 695, A130

    work page 2025

  45. [45]

    F., Cutri, R

    Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

    work page 2006

  46. [46]

    Ulrich, R. K. 1976, ApJ, 210, 377 V orobyov, E. I. & Basu, S. 2006, ApJ, 650, 956

    work page 1976

  47. [47]

    Mace, G., et al

    Wang, T., Li, J., N. Mace, G., et al. 2023, ApJ, 957, 8

    work page 2023

  48. [48]

    2000, A&AS, 143, 9

    Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9

    work page 2000

  49. [49]

    V ., Balser, D

    Wenger, T. V ., Balser, D. S., Anderson, L. D., & Bania, T. M. 2018, ApJ, 856, 52

    work page 2018

  50. [50]

    W., Roellig, T

    Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1

    work page 2004

  51. [51]

    L., Eisenhardt, P

    Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

    work page 2010

  52. [52]

    2007, ApJ, 669, 483

    Zhu, Z., Hartmann, L., Calvet, N., et al. 2007, ApJ, 669, 483

    work page 2007

  53. [53]

    Zhu, Z., Jiang, Y .-F., & Stone, J. M. 2020, MNRAS, 495, 3494 Article number, page 11 A&A proofs:manuscript no. aa56647-25 Appendix A: Fitting procedure The fitting procedure is based on Bayesian inference. The obser- vational data used in the fitting consist of two components: (i) the ground-based near-infrared spectra and (ii) the SPHEREx low- resolutio...

    work page 2020