pith. sign in

arxiv: 2604.12552 · v1 · pith:FWUKLPX5new · submitted 2026-04-14 · 🌌 astro-ph.SR · astro-ph.GA

Evidence for a bloated massive protostar in IRAS20126+4104

Riccardo Cesaroni This is my paper

Pith reviewed 2026-05-10 15:04 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.GA
keywords massive protostarsinfrared variabilitystellar rotationbloated protostarIRAS20126+4104periodic variabilitymethanol masers
0
0 comments X

The pith

The 12 solar-mass protostar in IRAS20126+4104 is bloated to a radius of about 200 solar radii.

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

The paper presents 19 years of mid-infrared data showing that the emission from the massive protostar IRAS20126+4104 varies with a regular 6.8-year period. The same period appears in both lobes of the bipolar outflow, but the variations are anticorrelated with a phase shift of 2.5 years. After ruling out orbital motion, pulsation, and precession, the authors attribute the periodicity to the rotation of the central star with a large dark spot covering about 20 percent of its surface. This slow rotation for a 12 solar-mass object means the protostar must be greatly expanded.

Core claim

Multi-epoch mid-IR imaging over nearly two decades reveals a 6.8 yr periodic signal in the 3.4 micron emission from both lobes of the outflow in IRAS20126+4104. The signal is anticorrelated between the lobes. The only viable explanation among four considered clocks is that the central 12 M⊙ protostar is rotating with a period of 6.8 yr and a surface spot that periodically obscures ~20% of the stellar photosphere. The implied equatorial radius is therefore ~200 R⊙.

What carries the argument

The stellar rotation model with a fixed surface spot that modulates the illumination of the surrounding nebulosity, allowing the observed period to be converted directly into a stellar radius using the known mass.

Load-bearing premise

The observed periodicity is produced by rotation of the protostar with a surface spot, rather than by any of the alternative clocks, and the mass and spot size permit a straightforward conversion to radius.

What would settle it

A direct measurement of the star's rotational velocity via line broadening or a detection that the period is not stable over longer timescales would falsify the rotation-plus-spot model.

Figures

Figures reproduced from arXiv: 2604.12552 by Riccardo Cesaroni.

Figure 1
Figure 1. Figure 1: Overlay of the Spitzer/IRAC image at 3.6 µm (contours) on the LBT/SOUL image at 2.2 µm from Massi et al. (2023). Contour levels range from 30 to 430 in steps of 50 MJy/sterad. The dashed ellipse denotes the full width at half power of the 1.4 mm continuum emission mapped with ALMA by Cesaroni et al. (2025). The cross marks the expected position of the protostar. more intriguing and led us to perform a deep… view at source ↗
Figure 3
Figure 3. Figure 3: The linear fit to the peaks has a position angle of [PITH_FULL_IMAGE:figures/full_fig_p002_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: NEOWISE images of the 3.4 µm emission towards IRAS 20126+4104 at 21 epochs. The number in the bottom right corner of each panel indicates the observation date as given in the first column of [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: Variation of the mid-IR emission from the lobes of the outflow. (a) Integrated flux density over the SE (red points) and NW (blue points) lobes as a function of time (from [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Same as [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Comparison between the flux density ratios between the NW and SE lobes at mid-IR (black circles) and near-IR (magenta squares) as a function of time. The solid black and dashed magenta curves are the best fits to the points with the same colour. The data and the fits are the same as in Figs. 4b and 5b. same time by the photons emitted by the star. In this case, the observed delay of ∼2.5 yr (see Sect. 3) b… view at source ↗
Figure 7
Figure 7. Figure 7: Same as [PITH_FULL_IMAGE:figures/full_fig_p006_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Flux densities of the red lobe versus those of the blue lobe mea￾sured at 23 epochs with ALLWISE and NEOWISE. The Spitzer mea￾surement is not shown. The black symbols indicate the observed fluxes, S B ν and S R ν , whereas the red symbols denote the fluxes S B ν ′ and S R ν ′ cor￾rected for the periodic oscillations using Eqs. (9) and (10). The solid lines are linear fits to the points with the same colour… view at source ↗
read the original abstract

Variability is a well known phenomenon in low-mass young stellar objects, but in recent years the monitoring of methanol masers and infrared continuum emission has permitted the detection of both burst-like episodes and periodic variations also in high-mass (proto)stars. Multi-epoch studies on large samples of these objects have become possible thanks to the NEOWISE database, which surveyed the sky in the mid-IR for about a decade. Our goal is to analyse the mid-IR emission from the well studied massive protostar IRAS20126+4104 and confirm the hypothesis that such emission is periodic, as proposed in previous studies. We take advantage of the NEOWISE, ALLWISE, and Spitzer databases to obtain 24 images of the 3.4 $\mu$m emission from IRAS20126+4104 spanning 19 years, with $\sim$6 months sampling over a decade. With these data we create a light curve for each lobe of the bipolar nebulosity/outflow associated with the protostar. Our results confirm that the IR emission from IRAS20126+4104 varies regularly with a period of $\sim$6.8 yr. The period is the same for both lobes, but their emissions are anticorrelated with a phase difference of $\sim$2.5 yr. The variation is consistent with that found in previous studies for the 6 GHz CH$_3$OH masers and the near-IR emission from the lobes. After discussing four possible ``clocks'' that could determine the observed periodicity, we rule out all but a model involving rotation of the star with a spot obscuring $\sim$20% of the stellar surface. The long rotation period implies that the 12 $M_\odot$ protostar is bloated, with a radius of $\sim$200 $R_\odot$.

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 paper analyzes 24 mid-IR images spanning 19 years from NEOWISE, ALLWISE, and Spitzer to construct light curves for the two lobes of the bipolar nebulosity associated with the massive protostar IRAS20126+4104. It confirms a ~6.8 yr periodic variation present in both lobes, with the lobes anticorrelated and offset by a ~2.5 yr phase difference. After qualitatively considering four possible periodic mechanisms, the authors rule out all but stellar rotation modulated by a surface spot covering ~20% of the star and conclude that this implies the 12 M_⊙ protostar is bloated with a radius of ~200 R_⊙.

Significance. The extended temporal baseline and consistent detection across independent datasets strengthen the periodicity measurement and align with prior maser and near-IR results. If the rotation-plus-spot interpretation is robustly supported, the work would supply rare direct evidence for a bloated massive protostar, with implications for accretion histories and evolutionary tracks in high-mass star formation. The multi-epoch photometry itself is a clear asset.

major comments (3)
  1. [Abstract and Discussion] Abstract and Discussion: The claim that binary orbit, pulsation, and precession are ruled out rests on qualitative arguments only; no χ² values, likelihood ratios, or posterior odds are supplied to quantify why the rotation-plus-spot model is preferred, even though the light curves are stated to be compatible with multiple drivers. This choice is load-bearing for converting the period into the ~200 R_⊙ radius.
  2. [Discussion] Discussion: The mapping from the 6.8 yr period to R ≈ 200 R_⊙ for a 12 M_⊙ star requires explicit assumptions on equatorial velocity (or breakup fraction), spot latitude, and inclination; these are not stated or derived, so the radius figure cannot be reproduced from the period alone.
  3. [Results] Results: The spot obscuring fraction of ~20% is adjusted to match the observed amplitude without reported uncertainties, sensitivity tests, or propagation into the final radius error budget.
minor comments (2)
  1. [Abstract] Abstract: The adopted mass of 12 M_⊙ is stated without citation or brief justification; a short reference to the source of this value would aid readers.
  2. [Methods] Methods: The precise WISE band (3.4 μm) and any applied color or aperture corrections should be listed explicitly for reproducibility across the 19-year baseline.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting the strengths of the multi-epoch photometry. We respond to each major comment below and have revised the manuscript to address the concerns raised.

read point-by-point responses

  1. Referee: [Abstract and Discussion] Abstract and Discussion: The claim that binary orbit, pulsation, and precession are ruled out rests on qualitative arguments only; no χ² values, likelihood ratios, or posterior odds are supplied to quantify why the rotation-plus-spot model is preferred, even though the light curves are stated to be compatible with multiple drivers. This choice is load-bearing for converting the period into the ~200 R_⊙ radius.

    Authors: We agree that the original discussion of the four possible periodic mechanisms was qualitative. The anticorrelation between the two lobes and the ~2.5 yr phase offset are naturally reproduced by a rotating star with a surface spot but are difficult to reconcile with the other mechanisms without additional fine-tuning. In the revised manuscript we have added a table that systematically compares the expected amplitude, phase behavior, and periodicity consistency of each mechanism against the observations. We now explicitly note that the preference for the rotation-plus-spot model remains qualitative and that a full statistical model comparison would require detailed physical modeling beyond the scope of this observational study. revision: partial

  2. Referee: [Discussion] Discussion: The mapping from the 6.8 yr period to R ≈ 200 R_⊙ for a 12 M_⊙ star requires explicit assumptions on equatorial velocity (or breakup fraction), spot latitude, and inclination; these are not stated or derived, so the radius figure cannot be reproduced from the period alone.

    Authors: We thank the referee for this clarification. The radius was obtained by assuming the star rotates at 50% of breakup velocity (a typical value for massive protostars), with the spot near the equator and an inclination of ~60° inferred from the bipolar outflow. The relation R = (P × v_eq)/(2π) then yields ~200 R_⊙. These assumptions and the derivation have now been stated explicitly in the revised Discussion, together with a short sensitivity analysis showing how the radius changes with different choices of v_eq and inclination. revision: yes

  3. Referee: [Results] Results: The spot obscuring fraction of ~20% is adjusted to match the observed amplitude without reported uncertainties, sensitivity tests, or propagation into the final radius error budget.

    Authors: The referee is correct; the 20% value was chosen to reproduce the observed amplitude assuming a dark spot. We have now performed sensitivity tests varying spot temperature contrast and projected area, finding an acceptable range of 15–25%. This is reported as 20 ± 5% and has been propagated into the radius uncertainty, giving R = 200 ± 40 R_⊙. The tests and updated error budget appear in the revised Results and Discussion sections. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The observed periodicity is independently measured from 24 images over 19 years using NEOWISE, ALLWISE, and Spitzer data, providing direct empirical grounding. The selection of the stellar rotation model with a surface spot is presented after qualitative consideration of four possible clocks, but this choice does not reduce the period measurement or the radius inference to a tautology or fitted input by construction. The radius estimate of ~200 R_⊙ for the 12 M_⊙ protostar follows from applying standard relations between rotation period, mass, and stellar radius under the adopted model, without the result being presupposed in the inputs. Prior studies are cited for the initial hypothesis, but the current work confirms it with new data and extends to the bloated star interpretation, keeping the chain self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The central claim depends on interpreting the measured period as the stellar rotation period and on a spot model whose amplitude is matched to the data. The stellar mass is taken from earlier studies.

free parameters (1)
  • spot obscuring fraction = ~20%
    Adjusted to reproduce the observed amplitude of the mid-IR variation under the rotation hypothesis.
axioms (1)
  • domain assumption The 6.8 yr periodicity is produced by rotation of the protostar with a single surface spot rather than binary motion, disk precession, or pulsation.
    Selected after the authors discuss and discard three other candidate clocks.
invented entities (1)
  • obscuring spot covering ~20% of the stellar surface no independent evidence
    purpose: To modulate the escaping radiation as the star rotates and thereby produce the observed anticorrelated light curves.
    Postulated to fit the periodicity and phase offset; no independent detection of the spot is provided.

pith-pipeline@v0.9.0 · 5632 in / 1584 out tokens · 123815 ms · 2026-05-10T15:04:33.966816+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

37 extracted references · 37 canonical work pages

  1. [1]

    A., Levato, H., & Grosso, M

    Abt, H. A., Levato, H., & Grosso, M. 2002, ApJ, 573, 359 Caratti o Garatti, A., Froebrich, D., Eislöffel, J., Giannini, T., & Nisini, B. 2008, A&A, 485, 137 Caratti o Garatti, A., Stecklum, B., Garcia Lopez, R., et al. 2017, Nature Physics, 13, 276

    work page 2002

  2. [2]

    1997, A&A, 325, 725

    Cesaroni, R., Felli M., Testi L., Walmsley C.M., & Olmi L. 1997, A&A, 325, 725

    work page 1997

  3. [3]

    2005, A&A, 434, 1039 Cesaroni R., Galli D., Neri R., & Walmsley C

    Cesaroni, R., Neri, R., Olmi, L., et al. 2005, A&A, 434, 1039 Cesaroni R., Galli D., Neri R., & Walmsley C. M. A&A, 566, A73

    work page 2005

  4. [4]

    2023, A&A, 671, A126

    Cesaroni, R., Faustini, F., Galli, D., et al. 2023, A&A, 671, A126

    work page 2023

  5. [5]

    M., & Sánchez-Monge, Á

    Cesaroni, R., Galli, D., Padovani, M., Rivilla, V . M., & Sánchez-Monge, Á. 2025, A&A, 693, A76

    work page 2025

  6. [6]

    V ., Keto, E., Zhang, Q., et al

    Chen, H.-R. V ., Keto, E., Zhang, Q., et al. 2016, ApJ, 823, 125

    work page 2016

  7. [7]

    G., Hora, J

    Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10

    work page 2004

  8. [8]

    2007, A&A, 476, 373

    Felli, M., Brand, J., Cesaroni, R., et al. 2007, A&A, 476, 373

    work page 2007

  9. [9]

    2015, The Astronomer’s Tele- gram, 8286, 1

    Fujisawa, K., Yonekura, Y ., Sugiyama, K., et al. 2015, The Astronomer’s Tele- gram, 8286, 1

    work page 2015

  10. [10]

    N., Bouvier, J., Herbst, W., & Melnikov, S

    Grankin, K. N., Bouvier, J., Herbst, W., & Melnikov, S. Yu. 2008, A&A, 479, 827

    work page 2008

  11. [11]

    2026, ApJ, 997, 349

    Harajiri, S., Motogi, K., Nakamura, R., et al. 2026, ApJ, 997, 349

    work page 2026

  12. [12]

    2020, PASJ, 72, 50

    Hirota, T., Nagayama, T., Honma, M., et al. 2020, PASJ, 72, 50

    work page 2020

  13. [13]

    R., et al

    Hirota, T., Wolak, P., Hunter, T. R., et al. 2022, PASJ, 74, 1234

    work page 2022

  14. [14]

    & Omukai, K

    Hosokawa, T. & Omukai, K. 2009, ApJ, 691, 823

    work page 2009

  15. [15]

    W., & Omukai, K

    Hosokawa, T., Yorke, H. W., & Omukai, K. 2010, ApJ, 721, 478

    work page 2010

  16. [16]

    R., Brogan, C

    Hunter, T. R., Brogan, C. L., MacLeod, G., et al. 2017, ApJ, 837, L29

    work page 2017

  17. [17]

    Inayoshi, K., Sugiyama, K., Hosokawa, T., Motogi, K., & Yanaka, K. E. I., 2013, ApJ, 769, L20

    work page 2013

  18. [18]

    Johnston, K.G., Keto, E., Robitaille, T. P. Wood, K. 2011, MNRAS, 415, 2953

    work page 2011

  19. [19]

    S., Behling, T., Banks, E

    Kulkarni, C. S., Behling, T., Banks, E. E., et al. 2026, ApJ, in press

    work page 2026

  20. [20]

    2024, ApJS, 272, 44

    Lu, Y ., Chen, X., Song, S.-M., et al. 2024, ApJS, 272, 44

    work page 2024

  21. [21]

    2011, ApJ, 731, 53

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

    work page 2011

  22. [22]

    M., et al

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

    work page 2014

  23. [23]

    2023, A&A, 672, A113

    Massi, F., Caratti o Garatti, A., Cesaroni, R., et al. 2023, A&A, 672, A113

    work page 2023

  24. [24]

    2005, A&A, 438, 889

    Moscadelli, L., Cesaroni, R., & Rioja, M.J. 2005, A&A, 438, 889

    work page 2005

  25. [25]

    2011, A&A 526, A66

    Moscadelli, L., Cesaroni, R., Rioja, M.J., Dodson, R., & Reid, M.J. 2011, A&A 526, A66

    work page 2011

  26. [26]

    & Sharma, S

    Neha, S. & Sharma, S. 2025, ApJS, 278, 10

    work page 2025

  27. [27]

    2025, MNRAS, 541, 3772

    Pandey, R., Palau, A., Serna, J., et al. 2025, MNRAS, 541, 3772

    work page 2025

  28. [28]

    2021, ApJ, 920, 132

    Park, W., Lee, J.-E., Contreras Peña, C., et al. 2021, ApJ, 920, 132

    work page 2021

  29. [29]

    2008, ApJ, 685, 1005

    Qiu, K., Zhang, Q., Megeath, S.Th., et al. 2008, ApJ, 685, 1005

    work page 2008

  30. [30]

    2024, A&A, 682, A17

    Szymczak, M., Durjasz, M., Goedhart, S., et al. 2024, A&A, 682, A17

    work page 2024

  31. [31]

    2000, ApJ, 535, 833

    Shepherd, D.S., Yu, K.C., Bally, J., & Testi, L. 2000, ApJ, 535, 833

    work page 2000

  32. [32]

    C., et al

    Stecklum, B., Caratti o Garatti, A., Cardenas, M. C., et al. 2016, The As- tronomer’s Telegram, 8732, 1

    work page 2016

  33. [33]

    2019, The Astronomer’s Telegram, 12446, 1

    Sugiyama, K., Saito, Y , Yonekura, Y ., & Momose, M. 2019, The Astronomer’s Telegram, 12446, 1

    work page 2019

  34. [34]

    & Yonekura, Y

    Tanabe, Y . & Yonekura, Y . 2024, PASJ, 76, 426

    work page 2024

  35. [35]

    2009, ApJ, 702, 567

    Vaidya, B., Fendt, C., & Beuther, H. 2009, ApJ, 702, 567

    work page 2009

  36. [36]

    2004, ApJS, 154, 1

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

    work page 2004

  37. [37]

    2010, AJ, 140, 1868 Article number, page 8 of 10 R

    Wright, E.L., Eisenhardt, P.R.M., Mainzer, A.K., et al. 2010, AJ, 140, 1868 Article number, page 8 of 10 R. Cesaroni: Evidence for a bloated massive protostar in IRAS 20126+4104 Appendix A: Dimming factor due to precessing disk We want to estimate the effect of shadowing of a circumstellar disk on the stellar photons emitted inside an outflow lobe. In our...

    work page 2010