O-type Bloated Star Candidate

• O-type bloated star candidates are massive young stellar objects characterized by inflated envelopes, lower effective temperatures, and pulsational variability due to high accretion rates.
• Multiwavelength observations, including optical, infrared, ALMA, and EVLA data, directly test the period–luminosity relation and reveal detailed disk–jet–outflow structures.
• Empirical data from IRAS 19520+2759 validate theoretical models, offering constraints on stellar mass, radius, and accretion dynamics in early high-mass star formation.

An O-type bloated star candidate refers to a massive, young, optically visible stellar object that exhibits characteristics predicted by theoretical models of protostellar evolution under high accretion rates. These characteristics include a luminosity exceeding $10^5 L_\odot$, effective temperature significantly below main-sequence O stars, a markedly inflated stellar radius, and distinctive photometric and spectroscopic variability. The massive young stellar object (MYSO) IRAS 19520+2759 currently stands as the archetype for such a candidate, integrating multiwavelength observational evidence with direct tests of the predicted period–luminosity (P–L) relation for bloated stars in the O-type mass regime (Pandey et al., 20 Mar 2025, Pandey et al., 18 Jan 2026).

1. Theoretical Foundation for O-type Bloated Stars

High-mass star formation models predict that massive young stellar objects accreting at high rates $$(\dot{M} \gtrsim 10^{-4} M_\odot\,\mathrm{yr}^{-1})$$ undergo a substantial swelling of their envelopes prior to reaching the zero-age main sequence (ZAMS). This "bloating" increases the stellar radius ($R_*$) and decreases the effective temperature ($T_\mathrm{eff}$) and ultraviolet output, as discussed by Hosokawa et al. (2009) and Kuiper et al. (2013). Inayoshi et al. (2013) demonstrated, via linear stability analysis, that such bloated envelopes become unstable to $\kappa$-mechanism pulsations, particularly in the He$^+$ ionization zone.

A key result of these models is a Period–Luminosity relation for pulsationally unstable, bloated MYSOs in the spherical accretion regime: $$\log\left(\frac{L}{L_\odot}\right) = 4.62 + 0.98 \log \left( \frac{P}{100\,\mathrm{days}} \right)$$ where $L$ denotes the stellar luminosity and $P$ the pulsation period (Pandey et al., 20 Mar 2025). For a given $L \sim 10^5 L_\odot$, this predicts long periodic variability on the order of several hundred days, directly testable by multi-epoch photometry.

2. Observational Diagnostics and Periodicity

Direct observational tests have focused on IRAS 19520+2759, utilizing time-series data across optical (TJO Rc, Ic; Gaia G, Rp, Bp), space-based (TESS), and mid-infrared (NEOWISE) bands. The measured periods are:

• Rc (TJO): $P_1 = 270 \pm 40$ days
• Ic (TJO): $P_2 = 272 \pm 50$ days
• G (Gaia): $P_3 = 460 \pm 60$ days (primary), $250 \pm 30$ days (secondary)
• Rp (Gaia): $P_4 = 440 \pm 70$ days (primary), $240 \pm 20$ days (secondary)
• TESS: $P_5 \simeq 6.02$ days (short, likely rotational)

Photometric amplitudes are $\Delta \mathrm{Rc} \simeq 0.30$ mag and $\Delta \mathrm{Ic} \simeq 0.25$ mag for the TJO bands. Color–magnitude diagrams from Gaia exhibit a negative slope ("bluer when fainter"), inconsistent with variable extinction or spots and supporting pulsational or accretion-driven intrinsic color change. The NEOWISE bands show large mid-IR amplitude ($\Delta \mathrm{W1}\sim1$ mag) but lack evidence for periodicity matching the optical bands (Pandey et al., 20 Mar 2025).

The observed periods align closely with the P–L relation in the theoretical models for $L \sim 10^5\,L_\odot$, providing empirical support for the existence of the predicted pulsational instability in a bloated O-type protostar.

3. Stellar and Circumstellar Parameters

Physical parameters have been inferred using both spherical and disc accretion models. For a representative period of 270–460 days in the spherical accretion scenario:

Parameter	Value (P=270 days)	Value (P=460 days)
Stellar Mass ($M_*$)	$24 \pm 1.3\,M_\odot$	$28 \pm 2\,M_\odot$
Radius ($R_*$)	$650 \pm 75\,R_\odot$	$900 \pm 115\,R_\odot$
Accretion Rate ($\dot{M}_*$)	$$(6.4 \pm 0.9) \times 10^{-3}\,M_\odot\,\mathrm{yr}^{-1}$$	$$(9 \pm 1.4) \times 10^{-3}\,M_\odot\,\mathrm{yr}^{-1}$$

However, the parameters are highly sensitive to the adopted $T_\mathrm{eff}$ and to the accretion geometry. Spectroscopic estimates suggest $T_\mathrm{eff} \gtrsim 17,000\,\mathrm{K}$, which, using the Stefan–Boltzmann law with $L=10^5\,L_\odot$, yields $R_* \simeq 36\,R_\odot$ and $M_* \simeq 15\,M_\odot$. Disc accretion models predict even smaller radii ($R_* \sim 30$–$40\,R_\odot$) and lower accretion rates than spherical estimates, with $M_* \lesssim 15\,M_\odot$, as discussed in Hosokawa & Yorke (2010) (Pandey et al., 20 Mar 2025).

4. High-Resolution Imaging: Disk–Jet–Outflow System

ALMA and EVLA observations have resolved the circumstellar environment of IRAS 19520+2759 in detail (Pandey et al., 18 Jan 2026). EVLA continuum imaging in C, K, and Q bands reveals a thermal jet characterized by a spectral index $\alpha=0.54 \pm 0.04$ (where $S_\nu \propto \nu^\alpha$), demonstrating consistency with collimated ionized jets from MYSOs. This emission extends along the east–west (P.A. ≈ 78°) axis, with a weaker component in the northeast–southwest (P.A. ≈ 28°) direction. ALMA 1.3 mm continuum maps reveal a compact core (MM1) spatially coincident (within 0.05″) with the optical counterpart. The core mass, calculated under optically thin assumptions, is $\sim 50\,M_\odot$ for $T_\mathrm{dust}=100$ K, establishing the presence of a massive reservoir for ongoing accretion.

Molecular line mapping in $^{13}$CO exposes two bipolar flows: a strong east–west outflow ($M_A \approx 30\,M_\odot$) and a weaker northeast–southwest outflow ($M_B \approx 10\,M_\odot$). The collimation factors, dynamical timescales ($t_\mathrm{dyn} \sim 1.5 \times 10^4$ yr), and energetics are congruent with other known MYSO jets (Pandey et al., 18 Jan 2026).

5. Disk Structure and Dynamical Mass

Multiple $\mathrm{SO}_2$ transitions, including high excitation lines ($E_u$ up to $\sim 200$ K), trace a compact hot core and a velocity gradient perpendicular to the dominant outflow, interpreted as a Keplerian disk. Position–velocity analysis yields a dynamical mass $M_\mathrm{dyn}=10$–$15\,M_\odot$ for the central object, assuming an edge-on orientation. Variations of the inclination by $\pm10^\circ$ alter the mass estimate by less than 30%, reinforcing placement of the star within the O-type regime.

Optical and near-IR spectroscopies indicate $T_\mathrm{eff} \approx 18,000$ K (B0–B1 type), $L_\mathrm{bol} > 10^5 L_\odot$, and $R_* \approx 25$–$35 R_\odot$. The surface gravity is $0.5$–$1.0$ dex lower than ZAMS O stars. These traits necessitate a bloated photosphere: SED modeling shows that only a swollen envelope (not a ZAMS model) can simultaneously account for the observed $L_\mathrm{bol}$ and $T_\mathrm{eff}$ (Pandey et al., 18 Jan 2026).

6. Distinction from Other High-Mass Protostars

IRAS 19520+2759 satisfies all theoretically motivated observational criteria for O-type bloated stars: high $L_\mathrm{bol}$ at relatively low $T_\mathrm{eff}$, absence of free–free emission, collimated molecular outflows, robust periodic variability, and a circumsystem disk–jet morphology. Gravity-sensitive spectral features and an inflated radius distinguish the bloated phase from both pre-main-sequence contraction and evolved stages. The negative color–magnitude trend ("bluer when fainter") in the optical strongly disfavors alternative explanations such as dust extinction variations or starspot modulations, supporting a pulsational or accretion-driven scenario.

The direct imaging of the disk–jet–outflow system in conjunction with characteristic periodicity and spectral properties represents the first unambiguous identification of an O-type bloated star candidate. The object serves as a critical anchor for empirical calibration of theoretical high-mass protostellar evolution models, particularly concerning accretion feedback and early pre-main-sequence development (Pandey et al., 20 Mar 2025, Pandey et al., 18 Jan 2026).

7. Significance and Future Prospects

The observation and characterization of IRAS 19520+2759 as an O-type bloated star candidate inaugurate the direct empirical testing of the period–luminosity relation for pulsating, highly accreting massive protostars. It establishes a benchmark for verifying models of high-mass star formation under extreme accretion, including the effects of bloating, pulsational instability, and feedback via jets and outflows.

Further multi-epoch spectroscopy, particularly targeting hydrogen recombination lines and higher-$T_\mathrm{eff}$ diagnostics, combined with high-resolution imaging of the disk and outflow kinematics, will refine estimates of stellar radius, mass, accretion rate, and pulsational properties. The distinction between spherical and disc-dominated accretion regimes remains crucial for accurately constraining the evolutionary status and internal structure of such objects. This system thus stands as the principal laboratory for advancing understanding of the transient, high-accretion bloated phase in O-star formation (Pandey et al., 20 Mar 2025, Pandey et al., 18 Jan 2026).