🔭 astrophysics

Evidence for a bloated massive protostar in IRAS20126+4104

By analyzing 19 years of mid-infrared data from the NEOWISE, ALLWISE, and Spitzer databases, researchers confirmed that the massive protostar IRAS20126+4104 exhibits a regular 6.8-year emission cycle caused by a rotating star with a large spot, implying the 12 $M_\odot$ object is bloated with a radius of approximately 200 $R_\odot$ .

Original authors: Riccardo Cesaroni

Published 2026-04-15

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Riccardo Cesaroni

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a massive, newborn star named IRAS 20126+4104. It's a "baby" star, but not a tiny one—it's already about 12 times heavier than our Sun. It's currently in a messy, chaotic phase of growing up, surrounded by a swirling disk of gas and dust, and shooting out powerful jets of material in opposite directions (like a cosmic sprinkler).

For years, astronomers have watched this star, but they noticed something strange. The star isn't just glowing steadily; it's pulsing in brightness, but not in the way a heartbeat works. Instead, it's behaving like a cosmic lighthouse with a very specific, weird rhythm.

Here is the story of what the paper discovered, broken down simply:

1. The Mystery of the "Dancing" Lobes

The star is surrounded by two giant clouds of gas (lobes) shooting out in opposite directions: one pointing Southeast (the "Red" lobe) and one pointing Northwest (the "Blue" lobe).

Astronomers used a space telescope called NEOWISE (which acts like a giant, long-term camera) to take pictures of this star every six months for nearly 20 years. They found something fascinating:

When the Red lobe gets bright, the Blue lobe gets dim.
When the Blue lobe gets bright, the Red lobe gets dim.

It's like a cosmic seesaw. They are perfectly out of sync. This happens on a strict schedule: the whole cycle takes about 6.8 years to complete.

2. The Detective Work: Ruling Out the Suspects

The astronomers asked: What is causing this 6.8-year seesaw? They tested four different theories, like a detective eliminating suspects:

Suspect #1: The Star is Pulsing (Breathing).
- The Idea: Maybe the star is expanding and contracting like a lung.
- The Verdict: Guilty of being impossible. If the star just breathed, both lobes would get brighter and dimmer at the exact same time. Since they are opposite, this can't be it. Also, the math says the star would have to be impossibly huge and bright to do this.
Suspect #2: The Disk is Wobbling (Precession).
- The Idea: Imagine a spinning top that starts to wobble. Maybe the disk of dust around the star is wobbling, casting a shadow that moves back and forth.
- The Verdict: Not the culprit. For this to work, the disk would have to be wobbling so violently that it would break apart. The physics just doesn't add up.
Suspect #3: An Orbiting Body (Revolution).
- The Idea: Maybe a giant, dusty cloud or a second object is orbiting the star and blocking the light, like a planet passing in front of the Sun.
- The Verdict: Unlikely. To orbit that close and take 6.8 years, the object would have to be right next to the star. But the star is so hot that it would vaporize any dust there instantly. No dust, no shadow.
Suspect #4: A Giant "Sunspot" (Rotation).
- The Idea: Our Sun has dark spots (sunspots) that rotate with it. Maybe this baby star has a giant, dark spot covering a huge chunk of its surface. As the star spins, this dark spot moves in and out of view, blocking the light from reaching one lobe, then the other.
- The Verdict: This is the winner! It fits all the clues. The spot blocks light to one side, then the star spins, and the spot blocks the other side.

3. The Big Reveal: A "Bloated" Star

If the "Giant Sunspot" theory is true, there's a catch. We know how fast stars usually spin. A normal star of this size should spin once every 1.5 days. But this star is taking 6.8 years to spin once!

How can a star spin that slowly?

The Analogy: Imagine an ice skater spinning with their arms tucked in. They spin fast. If they stretch their arms out wide, they slow down because their mass is spread out.
The Reality: This star is bloated. It has expanded to be about 200 times wider than our Sun. Because it is so huge and puffy, it spins very slowly, just like the ice skater with arms wide open.

This "bloated" state is exactly what theory predicts for massive stars that are still actively eating (accreting) gas from their surroundings. They puff up before they settle down into their final, compact shape.

4. The "Outburst" Bonus

The astronomers also noticed that one of the early photos (from the Spitzer telescope) showed the star was twice as bright as usual. This suggests that while the star is slowly spinning with its giant spot, it also has occasional "burps" or outbursts where it gobbles up a lot of gas at once, making it flash brightly. This is like the star having a snack in the middle of its slow dance.

Summary

The paper concludes that IRAS 20126+4104 is a massive, puffy, "bloated" baby star. It is spinning very slowly (once every 6.8 years) because it is so huge. It has a giant dark spot on its surface that acts like a cosmic lighthouse, alternately dimming the two jets of gas shooting out from its sides.

It's a rare glimpse into the chaotic, puffy, and slow-spinning teenage years of a massive star, proving that even the biggest stars can have "bad hair days" (or in this case, giant spots) that we can actually see from light-years away.

1. Problem Statement

While variability is well-documented in low-mass young stellar objects (YSOs), detecting periodic variability in high-mass protostars has been challenging. Previous studies identified burst-like episodes and some periodic variations in methanol masers and infrared (IR) emission, but the physical mechanisms driving long-term periodicity in massive protostars remain debated. The specific object of study, IRAS 20126+4104, is a well-studied massive protostar ( $\sim12 M_\odot$ ) driving a precessing bipolar jet and outflow. Recent studies (Szymczak et al. 2024) suggested a correlation between 6 GHz Class II methanol ( $CH_3OH$ ) maser variability and mid-IR emission, hinting at a periodic mechanism. However, the nature of this periodicity (pulsation, precession, rotation, or revolution) and the physical state of the protostar required further investigation to confirm the hypothesis and constrain the star's physical properties.

2. Methodology

The authors conducted a multi-epoch analysis of the mid-IR emission from IRAS 20126+4104 over a 19-year baseline (2005–2024).

Data Sources:
- NEOWISE: 21 high-quality single-exposure images at 3.4 $\mu$ m (W1 band) spanning 2014–2024, with $\sim$ 6-month sampling.
- ALLWISE: Two additional images at 3.4 $\mu$ m from 2010.
- Spitzer/IRAC: One image at 3.6 $\mu$ m from 2005 (Qiu et al. 2008).
- Near-IR Data: Re-analysis of 2.2 $\mu$ m data from Massi et al. (2023) and Cesaroni et al. (1997) to extend the timeline.
Data Processing:
- Images within a few days were averaged to improve signal-to-noise.
- Flux densities were measured within two specific polygons defining the North-West (NW/blue-shifted) and South-East (SE/red-shifted) lobes of the bipolar outflow.
- Peak positions of the emission were determined via 2-D Gaussian fitting.
- Spitzer data were smoothed to match the resolution of WISE/NEOWISE data.
Analysis:
- Light curves were constructed for each lobe.
- Sinusoidal fitting was applied to determine periods, phase differences, and amplitude variations.
- Four physical models (Pulsation, Precession, Rotation, Revolution) were tested against the observational constraints.

3. Key Results

The analysis confirmed and refined the periodic nature of the source:

Periodicity: The mid-IR emission varies regularly with a period of $T \approx 6.8$ years. This is consistent with the period derived for 6 GHz methanol masers ( $6.9 \pm 0.2$ yr) and previous near-IR studies.
Anticorrelation: The emissions from the two lobes are anticorrelated. When the SE lobe brightens, the NW lobe dims, and vice versa.
Phase Difference: The two light curves exhibit a phase difference of $\sim2.5$ years (approx. 0.37 of the period), rather than the perfect 0.5 (180°) expected for a simple binary or symmetric shadowing.
Amplitude: The variation (dimming factor) is approximately 40% of the peak flux.
Light Curve Shape: The variations are sinusoidal, with peak widths similar to dip widths.
Outburst Event: The Spitzer data point (2005) shows a flux $\sim$ 100% higher than the WISE/NEOWISE baseline, suggesting a separate episodic outburst (likely accretion-driven) superimposed on the periodic cycle.
Spatial Motion: The peak position of the mid-IR emission oscillates along a line parallel to the outflow axis, consistent with the movement of the emission source relative to the central protostar.

4. Discussion and Model Elimination

The authors evaluated four potential mechanisms ("clocks") to explain the 6.8-year periodicity and the specific constraints (periodicity, anticorrelation, 40% dimming, sinusoidal shape):

Stellar Pulsations: Rejected. Pulsations would affect both lobes simultaneously. The observed 2.5-year delay between lobes would require a light-travel time delay implying an outflow inclination angle ( $\alpha > 78^\circ$ ) that contradicts existing kinematic data ( $\alpha \approx 3^\circ-8^\circ$ ). Furthermore, the required luminosity for a 6.8-year pulsation period would be $10^6 L_\odot$ , far exceeding the source's actual luminosity ( $\sim10^4 L_\odot$ ).
Disk Precession: Rejected. For a precessing disk to cause anticorrelated variability, it must be asymmetric (opaque on one side, transparent on the other). Quantitative modeling showed this would require a disk radius $R_0 \le 0.28$ au to remain coherent. This is smaller than the dust sublimation radius for a $12 M_\odot$ star, and the rotation period at this radius ( $\sim16$ days) is far too short to maintain a static asymmetric sector over a 6.8-year precession cycle.
Orbital Revolution (Shadowing): Rejected. A body orbiting at $\sim8$ au could theoretically cause the shadowing. However, dust at this distance from a $12 M_\odot$ star would be sublimated (sublimation radius $>10$ au). Additionally, a compact body would produce sharp flux drops, inconsistent with the observed smooth sinusoidal light curves.
Stellar Rotation with a Spot: Accepted. This is the only viable model.
- Mechanism: A large, dark spot (or region of lower temperature) on the protostar's surface rotates with the star. As the spot rotates, it alternately obscures the view of the star from the two lobes, causing the anticorrelated brightness variations.
- Geometry: The spot covers $\sim20\%$ of the stellar surface ( $\Delta\theta \approx 24^\circ$ ). The star's rotation axis is inclined relative to the disk/outflow axis.
- Implication for Stellar Radius: Massive protostars are expected to be "bloated" (expanded) during active accretion. A standard $12 M_\odot$ main-sequence star rotates in $\sim1.5$ days. To have a rotation period of 6.8 years, the star must be significantly larger. Using angular momentum conservation ( $P \propto R^2$ ), the authors calculate a stellar radius of $R \approx 200 R_\odot$ . This size is consistent with theoretical predictions for massive protostars undergoing disk accretion (Hosokawa et al. 2010).

5. Significance and Conclusions

Confirmation of Bloated Protostars: The study provides robust observational evidence for the existence of "bloated" massive protostars. The derived radius of $\sim200 R_\odot$ supports theoretical models where accreting massive stars expand significantly before settling onto the main sequence.
Physical State: The large radius explains the low effective temperature and faint Lyman continuum of IRAS 20126+4104, which accounts for the lack of a detectable HII region despite the star's high mass.
Methodology: The paper demonstrates that long-term monitoring of IR variability, combined with multi-wavelength data (masers, near-IR, mid-IR), is a powerful tool for characterizing the physical properties (radius, rotation, surface features) of embedded massive protostars that cannot be resolved directly.
Surface Convection: The presence of a large spot on a massive star suggests active surface convection, a phenomenon predicted by models of bloated protostars but rarely observed directly.

In summary, the paper concludes that IRAS 20126+4104 is a slowly rotating ( $P \approx 6.8$ yr), bloated massive protostar ( $R \approx 200 R_\odot$ ) with a large surface spot, and that the observed periodic variability is caused by the rotation of this spot obscuring the stellar flux alternately toward the two outflow lobes.