VVV-WIT-13: Eruptive Infrared YSO Outburst

• VVV-WIT-13 is an eruptive young star characterized by a rapid infrared outburst, significant photometric variations, and unique molecular absorption features such as deep AlO bands.
• The study combines multi-survey photometry and spectroscopy to detail the outburst evolution, revealing a transition from cool molecular absorption to broad H2 emissions indicative of outflow dynamics.
• Hydrodynamic simulations and observed blueshifted absorption lines support a scenario where a tidally disrupted gas clump in the disk triggers accretion instabilities and drives a distinct outflow environment.

VVV-WIT-13 is an infrared eruptive source identified in the decade-long VISTA Variables in the Via Lactea survey and assigned to the VVV “WIT” class, where WIT stands for “What Is This?” because of its unusual photometric variation behaviour and the initially unclear origin of the outburst. The source is now interpreted as an eruptive young star with instability occurring in the accretion disk, distinguished by a large-amplitude infrared eruption, an intermediate-duration light curve, and unusually cool molecular absorption features—most notably deep AlO absorption during the outburst stage—followed by spectroscopic signatures of a wind or outflow during decline (Guo et al., 18 Sep 2025).

1. Identification, nomenclature, and survey setting

VVV-WIT-13 is located at RA $=16{:}53{:}44.38$, Dec $=-43{:}28{:}19.47$ and has appeared previously under the names VVVv746, WISEA J165344.39-432819.2, and L222_59 (Guo et al., 18 Sep 2025). The designation places it within the VVV WIT family, a naming convention used for rare VVV variables whose nature was initially uncertain. A closely related example is VVV-WIT-12, another unusual VVV source for which the term was explicitly glossed as “What Is This?”, and whose eventual interpretation as a deeply embedded variable YSO illustrates that the WIT label is a discovery category rather than a single physical class (Saito et al., 2023).

The observational basis for VVV-WIT-13 combines VVV and VVVX near-infrared monitoring with archival and follow-up mid-infrared and spectroscopic data. The photometric analysis uses VVV/VVVX, 2MASS, Spitzer/GLIMPSE, MIPSGAL, ALLWISE, and NEOWISE, together with follow-up imaging from SOFI/NTT, IRSF, REM, and a Magellan/FIRE $J$-band acquisition image (Guo et al., 18 Sep 2025). This multi-survey coverage is characteristic of VVV-based time-domain work in the obscured inner Milky Way, where the main variability information is typically carried by the near-infrared $K_s$ band (Hempel et al., 2014).

The source is interpreted as a low-mass embedded young star. Its pre-outburst properties place it between Class I and Class II / flat-spectrum evolutionary stages, and in the SPICY catalogue it is classified as a flat-spectrum YSO (Guo et al., 18 Sep 2025). The preferred distance adopted in the paper is 2 kpc, motivated by the surrounding YSO population and proper-motion constraints; the median parallax of nearby SPICY-group members is 0.487 mas, corresponding to 2.05 kpc (Guo et al., 18 Sep 2025).

2. Infrared outburst and time-domain evolution

Before the main eruption, VVV-WIT-13 was already variable. In the pre-outburst state the source showed $\Delta K_s = 1.5$ mag over a timescale of about 2000 d, together with $\Delta W1 = 1.6$ mag and $\Delta W2 = 0.5$ mag (Guo et al., 18 Sep 2025). A tentative Lomb–Scargle signal with period

$P = 1748 \pm 141~{\rm d}$

was extracted from the $K_s$-band data, although only one cycle is seen and the periodicity is therefore not secure (Guo et al., 18 Sep 2025).

The main outburst began after 25 August 2016. The rise was rapid but only sparsely sampled, so the rise time is constrained to

$180\text{--}360~{\rm d}.$

The abstract reports an outburst amplitude of 5.7 mag in the $=-43{:}28{:}19.47$0-band, a brightness plateau lasting 3.5 years, and a subsequent rapid fade to the pre-eruptive level (Guo et al., 18 Sep 2025). In the main text, the paper also reports context-dependent amplitude estimates of $=-43{:}28{:}19.47$1 mag, $=-43{:}28{:}19.47$2 mag, and $=-43{:}28{:}19.47$3 mag, while the colour analysis quotes $=-43{:}28{:}19.47$4 mag and $=-43{:}28{:}19.47$5 mag (Guo et al., 18 Sep 2025). The coexistence of these values reflects differing choices of reference epoch and dataset.

The plateau persisted for roughly 1000 days after reaching maximum, and the event as a whole is described as lasting about $=-43{:}28{:}19.47$6–3000 days, placing it between typical EXor and classical FUor timescales (Guo et al., 18 Sep 2025). The decline then became evident in 2021. The $=-43{:}28{:}19.47$7-band flux fell to below 14 mag in about two years, but the fade was not achromatic: the mid-IR light curves began to decline about 1 year earlier than the near-IR (Guo et al., 18 Sep 2025). By 2024 the source had reached a new lower plateau, still approximately 1 mag brighter than the pre-outburst state (Guo et al., 18 Sep 2025).

A notable photometric property is the strong mid-infrared response. The paper emphasizes that the mid-IR amplitudes are comparable to or larger than the near-IR amplitudes, which is unusual for classical FUor events (Guo et al., 18 Sep 2025). This is one of the reasons the outburst is regarded as atypical even within the broader class of accretion-driven YSO eruptions.

3. Colours, extinction, luminosity, and environment

The quiescent near-infrared colours are

$=-43{:}28{:}19.47$8

whereas during outburst they became

$=-43{:}28{:}19.47$9

Thus the source became brighter and bluer during the eruption (Guo et al., 18 Sep 2025). Although the colour change lies close to an extinction vector in the colour–colour plane, the observed amplitudes do not satisfy the expectation for pure extinction variation, since the measured

$J$0

do not match the relation $J$1 invoked in the paper’s discussion (Guo et al., 18 Sep 2025). The authors therefore argue that variable line-of-sight extinction is not the sole driver of the event.

Several extinction estimates are presented. Assuming intrinsic colours appropriate for a normal disk-bearing YSO, the near-infrared quiescent colours imply

$J$2

From the $J$3 diffuse interstellar band, the paper derives

$J$4

which is interpreted as tracing mainly the interstellar component. From the $J$5 line ratio $J$6 in 2023, the paper obtains

$J$7

The preferred interpretation is therefore that $J$8 mag is a lower limit and $J$9 mag is a plausible total extinction including circumstellar material (Guo et al., 18 Sep 2025).

Assuming $K_s$0 kpc and $K_s$1 mag, the pre-outburst bolometric luminosity is estimated as

$K_s$2

implying a stellar mass of roughly $K_s$3 for an age of 0.5–1 Myr (Guo et al., 18 Sep 2025). If the object were instead at 1 kpc, the luminosity would fall by 0.60 dex, implying a mass closer to $K_s$4 (Guo et al., 18 Sep 2025). During outburst, the bolometric luminosity rises to

$K_s$5

for $K_s$6 mag, and remains

$K_s$7

even for $K_s$8 mag (Guo et al., 18 Sep 2025).

The environmental context strongly supports a YSO interpretation. VVV-WIT-13 lies in a star-forming region containing molecular material, Spitzer dark clouds, Herschel Hi-GAL clumps, nearby H II regions, and a substantial SPICY YSO association (Guo et al., 18 Sep 2025). It is only 6″ from molecular cloud SDG G342.136+0.2045 and is projected against the Spitzer dark cloud G342.135+0.204 (Guo et al., 18 Sep 2025). The paper treats this setting as a major argument against a post-main-sequence interpretation.

4. Spectroscopy and cool molecular phenomenology

Near-infrared spectra were obtained during both the plateau and decay stages. The sequence consists of XSHOOTER/VLT spectra on 29 April 2021 and 3 May 2021, a FIRE/Magellan spectrum in July 2023, and a later XSHOOTER/VLT spectrum on 16 March 2024 (Guo et al., 18 Sep 2025). The 2021 plateau spectra were combined because no short-term spectral variability was detected.

During the outburst plateau, the spectrum is dominated by cool absorption features. The paper reports TiO absorption in the optical, deep AlO absorption bands in the $K_s$9 and $\Delta K_s = 1.5$0 bands, broad $\Delta K_s = 1.5$1 absorption through $\Delta K_s = 1.5$2, and strong $\Delta K_s = 1.5$3 first-overtone bandhead absorption beyond $\Delta K_s = 1.5$4 (Guo et al., 18 Sep 2025). There is no obvious detection of $\Delta K_s = 1.5$5 bandheads, which is used to argue against a chemically evolved progenitor (Guo et al., 18 Sep 2025). Hydrogen features include Pa$\Delta K_s = 1.5$6, Pa$\Delta K_s = 1.5$7, and Br$\Delta K_s = 1.5$8 absorption, together with broad He I $\Delta K_s = 1.5$9 absorption (Guo et al., 18 Sep 2025). Narrow, symmetric absorption lines of Mg II, Fe I, Si I, and Al I are also detected (Guo et al., 18 Sep 2025).

The kinematics are central to the paper’s interpretation. The expected system velocity at the preferred location is approximately $\Delta W1 = 1.6$0 km s$\Delta W1 = 1.6$1, corresponding to about $\Delta W1 = 1.6$2 km s$\Delta W1 = 1.6$3 heliocentric (Guo et al., 18 Sep 2025). By contrast, the molecular absorptions are strongly blueshifted:

• CO: $\Delta W1 = 1.6$4 km s$\Delta W1 = 1.6$5
• AlO: $\Delta W1 = 1.6$6 km s$\Delta W1 = 1.6$7

The hydrogen absorption lines are similarly blueshifted at about $\Delta W1 = 1.6$8 to $\Delta W1 = 1.6$9 km s$\Delta W2 = 0.5$0, whereas the narrow metal lines remain near the photospheric velocity (Guo et al., 18 Sep 2025). This separation between narrow metal-line velocities and molecular velocities is one of the strongest indications that the molecular absorbers are not simply photospheric.

The molecular bands were modeled using ExoMol line lists under the assumptions of a plane-parallel slab of gas in local thermal equilibrium and in the optically thin case (Guo et al., 18 Sep 2025). The most unusual result is the low temperature of the AlO absorber: $\Delta W2 = 0.5$1 The CO absorber is also cool by eruptive-YSO standards, with a best-fit temperature of about 1300 K in 2021 (Guo et al., 18 Sep 2025).

Epoch	Species	RV / temperature / column density
2021	CO	$\Delta W2 = 0.5$2 km s$\Delta W2 = 0.5$3, $\Delta W2 = 0.5$4 K, $\Delta W2 = 0.5$5
2021	AlO	$\Delta W2 = 0.5$6 km s$\Delta W2 = 0.5$7, $\Delta W2 = 0.5$8 K, $\Delta W2 = 0.5$9
2023	CO	$P = 1748 \pm 141~{\rm d}$0 km s$P = 1748 \pm 141~{\rm d}$1, $P = 1748 \pm 141~{\rm d}$2 K, $P = 1748 \pm 141~{\rm d}$3
2023	$P = 1748 \pm 141~{\rm d}$4	$P = 1748 \pm 141~{\rm d}$5 km s$P = 1748 \pm 141~{\rm d}$6, FWHM $P = 1748 \pm 141~{\rm d}$7 km s$P = 1748 \pm 141~{\rm d}$8, EW $P = 1748 \pm 141~{\rm d}$9
2024	$K_s$0	$K_s$1 km s$K_s$2, FWHM $K_s$3 km s$K_s$4, EW $K_s$5

The decay-stage spectra show a clear transformation. The AlO absorption disappeared, the CO absorption remained but became shallower, and broad blueshifted $K_s$6 1–0 S(1) $K_s$7 emission appeared (Guo et al., 18 Sep 2025). The paper identifies this $K_s$8 feature as a common outflow or wind diagnostic in YSOs. The observational evolution from AlO-dominated absorption to broad $K_s$9 emission is treated as direct evidence for a changing circumstellar flow.

5. Physical interpretation and theoretical models

The paper concludes that VVV-WIT-13 is best understood as an eruptive young star with instability occurring in the accretion disk (Guo et al., 18 Sep 2025). The classification is supported by the pre-outburst SED, embedded colours, star-forming environment, low inferred stellar mass, large infrared luminosity increase, and the spectroscopic transition from cool molecular absorption to wind/outflow tracers (Guo et al., 18 Sep 2025).

The object shows some similarities to FUor-like eruptions: a rapid rise of less than about a year, large near-infrared amplitude, broad molecular absorption, and a high outburst luminosity that can reach $180\text{--}360~{\rm d}.$0 for the preferred extinction (Guo et al., 18 Sep 2025). It also differs from classical FUors in several respects. The event duration is only multi-year rather than decades-long, the mid-IR amplitudes are as large as or larger than the near-IR amplitudes, the CO absorption is unusually cool, AlO absorption is unprecedented for eruptive YSOs, and broad $180\text{--}360~{\rm d}.$1 emission appears during the fading stage (Guo et al., 18 Sep 2025). These differences motivate the paper’s treatment of VVV-WIT-13 as an unusual, possibly new variant of accretion-driven YSO outburst behaviour.

Alternative interpretations are discussed and disfavoured. A red nova or merger-like scenario is weakened by the pre-outburst SED of an embedded YSO, the comparatively low luminosity unless a much larger distance is assumed, and the event timescale, which is longer than most red-nova plateaus though not without precedent (Guo et al., 18 Sep 2025). A pure extinction event is disfavoured by the colour-amplitude mismatch. A classical nova interpretation is not supported by the infrared SED or environment (Guo et al., 18 Sep 2025). The paper nevertheless remains cautious and notes that an unusual low-luminosity red-nova-like interpretation is not entirely excluded (Guo et al., 18 Sep 2025).

A central inference is that the CO, AlO, and probably TiO features arise in an outflow or wind environment rather than in a standard hot inner accretion disk (Guo et al., 18 Sep 2025). The argument is based on the strong blueshifts of the molecular absorptions relative to the narrow metal lines, the very low AlO temperature, and the later emergence of broad $180\text{--}360~{\rm d}.$2 emission (Guo et al., 18 Sep 2025). The paper explicitly states that the observational evidence suggests that the CO and TiO features originate from an outflow or a wind environment (Guo et al., 18 Sep 2025).

To explore the outburst mechanism, the authors performed 2D hydrodynamic simulations with FARGO3D using a time-dependent energy equation in the adiabatic approximation (Guo et al., 18 Sep 2025). The favoured scenario is the tidal disruption of a giant gas clump or planet embryo in the disk. In the abstract this is summarized as “a disrupted gas clump at a distance of 3 au from the source”, while the main text specifies a tidally disrupted $180\text{--}360~{\rm d}.$3 planet embryo (Guo et al., 18 Sep 2025). The possible pre-outburst period

$180\text{--}360~{\rm d}.$4

would correspond to a Keplerian radius of 2.8 au around a $180\text{--}360~{\rm d}.$5 star, or 1.8 au around a $180\text{--}360~{\rm d}.$6 star (Guo et al., 18 Sep 2025), which the paper treats as suggestive of the same few-au scale. If this interpretation is confirmed, the event would constitute the first such disrupted gas-clump event observed in real time (Guo et al., 18 Sep 2025).

The paper also gives an order-of-magnitude wind mass-loss estimate from the CO absorber: $180\text{--}360~{\rm d}.$7 with $180\text{--}360~{\rm d}.$8, $180\text{--}360~{\rm d}.$9, and $=-43{:}28{:}19.47$00, yielding

$=-43{:}28{:}19.47$01

for the adopted parameters (Guo et al., 18 Sep 2025). This value is described as smaller than red-nova ejecta rates but compatible with strong YSO outflows (Guo et al., 18 Sep 2025).

6. Relation to the wider VVV literature and open issues

The designation VVV-WIT-13 is established in the dedicated eruptive-source study (Guo et al., 18 Sep 2025). Earlier VVV papers on proper motions, clusters, survey status, or bulge windows do not identify the source or define that name [(Kurtev et al., 2016); (Hempel et al., 2014)]. In particular, the paper on VVV WIN 1733$=-43{:}28{:}19.47$023349—also called Oscar’s window—explicitly states that it does not use the designation VVV-WIT-13, and therefore does not provide evidence that the two refer to the same target (2002.04329). This distinction is important because the low-extinction window literature concerns Galactic structure, not an eruptive infrared star.

Within the WIT class, VVV-WIT-13 now occupies a distinctive position. VVV-WIT-12 established that WIT sources can include deeply embedded infrared variables with unusual nebular behaviour and uncertain early classification (Saito et al., 2023). VVV-WIT-13 extends that pattern to an eruptive object whose decisive anomaly is not only its light curve but also its molecular spectrum, especially the transient AlO absorption at 600 K (Guo et al., 18 Sep 2025). This suggests that the WIT family is best understood as a survey-discovery reservoir of astrophysically heterogeneous outliers.

Several uncertainties remain. The distance to VVV-WIT-13 is argued rather than directly measured; the extinction is uncertain between roughly $=-43{:}28{:}19.47$03 and $=-43{:}28{:}19.47$04; the proposed 1748 d pre-outburst period is based on only one cycle; the molecular modelling assumes LTE, plane-parallel geometry, and optically thin absorption; and the hydrodynamic explanation is explicitly presented as ad hoc rather than as a fully explored parameter survey (Guo et al., 18 Sep 2025). The paper therefore calls for continued photometric monitoring, longer-wavelength observations to constrain disk geometry and mass, and more detailed theoretical work (Guo et al., 18 Sep 2025).

The later observational trajectory appears to remain active. The paper notes newer spectroscopy not analyzed in detail there, in which the CO bandheads have reportedly turned from absorption to emission, and accretion indicators such as Br$=-43{:}28{:}19.47$05 and Na I have appeared (Guo et al., 18 Sep 2025). This suggests that VVV-WIT-13 may continue to illuminate the transition between embedded-disk accretion physics, molecular outflows, and transient circumstellar chemistry.