Context. Gravitational fragmentation in massive protostellar disks can lead to the formation of bound gaseous clumps whose inward migration and disruption may trigger strong accretion outbursts. Yet the roles of numerical resolution and inner boundary treatment in shaping the burst properties remain poorly understood. Aims. We investigate how resolving the inner few astronomical units of a massive protostellar disk affects the migration, disruption, and accretion signatures of an inward-moving fragment. In particular, we aim to determine whether the predicted burst strength and duration depend on the adopted sink cell size. Methods. We present a new three-dimensional radiation-hydrodynamic simulation of a ∼5 M⊙ protostar surrounded by a self-gravitating disk that compares the original 30 AU sink model to a refined model with a 1 AU sink that resolves the inner disk. The resulting gas structures were post-processed with radiative transfer calculations to derive synthetic photometry and multiband images. Results. We find the both simulations produce a major accretion burst as a migrating fragment is tidally disrupted, but their detailed behavior differs markedly. The refined model shows faster migration, a complete tidal disruption of the fragment, and a shorter, sharper outburst (more consistent with observations) with nearly the same peak accretion rate as the 30 AU model, which yields a broader, smoother event. The refined run also produces much stronger near- and mid-infrared emission, reflecting the formation of a compact, hot inner disk. Conclusions. Resolving the inner few AU qualitatively changes the dynamics and observable appearance of fragment-driven bursts. Diffuse fragment disruption can reproduce decade-long events, but the much shorter (< 3 yr) bursts observed in some massive protostars likely require the tidal disruption of more compact objects, such as second Larson cores. Our trajectory analysis indicates that second Larson cores can migrate sufficiently close to the star to be tidally destroyed, offering a plausible mechanism for the fastest FU-Ori–like bursts observed in massive protostars.

The binary T Tauri system, KH 15D, is poised to provide unprecedented detail of gas and ices located between ∼3-5 AU of its circumbinary disk over the next few years. We analyze for the first time a complete set of ground-based time series photometry of the KH 15D system for which the 2022 October - 2025 April data have never been presented. By combining the latest I-band data with historical measurements, we show that the era of large amplitude photometric variability has ended. We define a revised photometric period of P = 48.365 ± 0.005 days, which accurately phases data over a time span of seventy-five years. We adopt magnitudes and colors for the stars that define the epochs when either star A (1995-2008) or star B (2010-2012) is the primary source of the out-of-eclipse brightness of the system. We determine for the first time the effective optical depth for the system between 2002 and 2025 assuming gray extinction. The current out-of-eclipse system color is consistent with a combination of the light from both stars, supporting the previously suggested model in which the trailing edge of the circumbinary disk is more transparent than the leading edge. We continue to see evidence for “clumps” and variable transparency along ingress and egress when the stars are probing the occulting material and suggest that observations made during these phases should permit studies of the physical and chemical nature of the planet-forming zone with the circumbinary disk.

We present a homogeneous analysis of rotational variability and spot properties in young stellar objects across multiple star-forming regions observed by the Hunting Outbursting Young Stars (HOYS) project. From over 2000 candidate members, we identify 144 YSOs with robust periodic signals and well-constrained multi-band amplitudes. The sample has a median age of $\sim$1~Myr, effective temperatures of 3500--6500~K (masses $\sim$0.6--2~M$_\odot$), and is dominated by Class~2 objects, one third of which exhibit inner disc dust emission. The rotation period distribution is strongly bimodal, with 55 percent fast rotators ($P<5.5$~d) and 45 percent slow rotators. Fast rotators are predominantly inner disc-less, whereas slow rotators include both disc-bearing and disc-free systems, indicating that disc braking alone cannot explain the observed rotational states. We derive spot properties from multi-band amplitudes and find that, after correcting for observational biases, the intrinsic cold-spot coverage distribution of fast rotators is well described by an exponential function. This implies that small spot coverages are intrinsically much more common than large ones, consistent with stochastic magnetic flux emergence governing spot formation. In contrast, slow rotators show a pronounced deficit of small cold spots. After considering observational biases and alternative physical explanations, we conclude that small spots on slowly rotating YSOs have significantly shorter lifetimes. These results provide new evidence that magnetic surface structure and its evolution depend on stellar rotation, placing new empirical constraints on models of magnetic activity and angular momentum evolution in young stars.

We investigate how resolving the inner few astronomical units of a massive protostellar disk affects the migration, disruption, and accretion signatures of an inward-moving fragment. In particular, we aim to determine whether the predicted burst strength and duration depend on the adopted sink cell size. We present a new three-dimensional radiation-hydrodynamic simulation of a ∼5M_⊙ protostar surrounded by a self-gravitating disk, comparing the original 30 AU sink model to a refined model with a 1 AU sink that resolves the inner disk. The resulting gas structures are post-processed with radiative transfer calculations to derive synthetic photometry and multi-band images. Both simulations produce a major accretion burst as a migrating fragment is tidally disrupted, but their detailed behavior differs markedly. The refined model shows faster migration, a complete tidal disruption of the fragment, and a shorter, sharper outburst (more consistent with observations) with nearly the same peak accretion rate as the 30 AU model, which yields a broader, smoother event. The refined run produces much stronger near- and mid-infrared emission, reflecting the formation of a compact, hot inner disk. Resolving the inner few AU qualitatively changes the dynamics and observable appearance of fragment-driven bursts. Diffuse fragment disruption can reproduce decade-long events, but the much shorter (<3 yr) bursts observed in some massive protostars likely require the tidal disruption of more compact objects such as second Larson cores. Our trajectory analysis indicates that second Larson cores can migrate sufficiently close to the star to be tidally destroyed, offering a plausible mechanism for the fastest FU-Ori-like bursts observed in massive protostars.

We analyze the H _2 emission observed in the HH 46 Class I system as part of PROtostellar JEts Cradle Tested with JWST (PROJECT-J), to investigate the origin and excitation of the warm molecular outflow. We used NIRSpec and MIRI spectral maps (1.6–27.9 μ m) to trace the structure and physical conditions of the outflow. By fitting the H _2 rotational diagrams with a multitemperature gas model, we derived key physical parameters, including temperature, extinction, column densities, and the ortho-to-para ratio. This information is combined with a detailed kinematical analysis and comparison with irradiated shock models. We find no evidence of H _2 temperature or velocity stratification from the axis to the edge of the outflow, as would be expected in MHD disk-wind models and as observed in other outflows. Instead, the observations suggest that the H _2 emission arises from shock interactions of jet bow shocks and/or wide-angle winds with the ambient medium and cavity walls. NIRSpec emission and velocity maps reveal expanding molecular shells, likely driven by the less luminous source in the binary system. We infer an accretion rate of ≲10 ^−9 M _⊙ yr ^−1 for the secondary source, approximately one order of magnitude lower than that of the primary. The H _2 emission is consistent with excitation by low-velocity (∼10 km s ^−1 ) J-type shocks, irradiated by an external UV field that may originate from strong dissociative shocks driven by the atomic jet. Future JWST observations will further constrain the evolution of the expanding shell and the mechanisms driving the outflow.

We use far-ultraviolet spectra of 36 T Tauri stars, predominately from the Hubble Space Telescope (HST) ULLYSES program, to examine the kinematic properties of fluorescent H _2 emission lines for evidence of disk outflows. Leveraging improvements to the HST Cosmic Origins Spectrograph wavelength solution, we coadd isolated lines within four fluorescent progressions ([ v’ , J’ ] = [1,4], [1,7], [0,2], and [3,16]) to improve signal-to-noise ratio (S/N), and we fit each coadded line profile with one or two Gaussian components. Of the high-S/N line profiles (S/N ≥ 12 at the peak of the profile), over half are best fit with a combination of a broad and a narrow Gaussian component. For profiles of the [1,4] and [1,7] progressions, we find a systematic blueshift of a few kilometers per second between the broad and narrow centroid velocities and stellar radial velocities. For the [0,2] progression, we find centroid velocities consistently blueshifted with respect to stellar radial velocities on the order of −5 km s ^−1 for the single and narrow components, and −10 km s ^−1 for the broad components. Overall, the blueshifts observed in our sample suggest that the molecular gas traces an outflow from a disk wind in some sources, and not solely disk gas in Keplerian rotation. The low-velocity systematic blueshifts, as well as emitting radii as inferred from line FWHMs, observed in our sample are similar to those observed with optical [O i ] surveys of T Tauri stars. We estimate H _2 mass-loss rates of 10 ^−9 to 10 ^−11 M _⊙ yr ^−1 , but incomplete knowledge of wind parameters limits comparisons to global models.

How a star forms is a fundamental question in astrophysics. In the earliest stages of protostellar evolution high extinction prevents a direct study of the accretion processes and their temporal evolution. Monitoring the variations of the accretion luminosity in a large protostar sample over decades is needed to reveal how protostars accrete – in major bursts or in a quasi-steady fashion. We here argue that a large ground based sub-millimeter single-dish facility with a wide FoV is required to fulfill this task.

Young stellar objects (YSOs) can display unpredictable and high-amplitude rises in brightness that can last from a few months to possibly over 100 years. These types of outbursts are explained by large changes in the mass accretion rate from the disk onto the central star. This type of variability has given support to a model of star formation (episodic accretion) where stars would spend most of their lifetimes accreting at low rates, and gain most of their mass through these short-lived accretion outbursts. The universality of episodic accretion, as well as its potential impact on stellar and planetary formation are still under debate. Improvement on the statistics of the members of the eruptive class is needed to better understand the episodic accretion phenomenon and its universality across different mass regimes and environments. In this paper we collect published information on the spectroscopic and photometric characteristics of 174 YSOs confirmed to belong to the eruptive variable class. We classify these objects into five different sub-classes (we find 49 FUor, 20 FUor-like, 16 EX Lupi-type, 81 Peculiar/V1647 Ori-like/MNors and 8 Periodic YSOs). The classification follows what has been done previously in the literature, and it is not an attempt to redefine these classes. In addition, we present a list of 18 embedded, and 6 massive YSOs, as additional categories of eruptive variable YSOs. Due to the complexity and/or faintness of these systems, it is hard to place them into the original classification scheme of this class of variable YSOs. Finally, we present a separate list of 355 candidate eruptive variable YSOs, which either lack spectroscopic information or the available spectroscopic data is not sufficient for an unambiguous classification. The online catalogue of confirmed and candidate eruptive YSOs will be maintained and updated in the future to serve as an important reference for the star formation community.

The main goal of this study is to screen the PENELLOPE/UVES targets for outflow activity and find microjets via spectro-astrometry in, e.g., the [OI]λ6300 line. In total, 34 T Tauri stars of the PENELLOPE survey have been observed with the high resolution slit spectrograph UVES in three different slit positions rotated by 120^o. Our spectro-astrometric analysis in the [OI]λ6300 wind line reveals two newly discovered microjets associated with Sz 103 and XX Cha. Both microjets have an extent of about 0.04 arcseconds, that is, <10 au, and we confined their orientation by the three slit observations. Furthermore, we confirm the binary nature of VW Cha and CVSO 109. We present (further) evidence that DK Tau B and CVSO 104 A are spectroscopic binaries. Sz 115 is tentatively a spectroscopic binary. The origin of the LVC, that is, MHD winds versus photoevaporative winds, of the Sz 103 and XX Cha microjets remains unclear.

The Classical T Tauri Star (CTTS) stage is a critical phase of the star and planet formation process. In an effort to better understand the mass accretion process, which can dictate further stellar evolution and planet formation, a multi-epoch, multi-wavelength photometric and spectroscopic monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was carried out in 2021 and 2022/2023 as part of the Outflows and Disks Around Young Stars: Synergies for the Exploration of ULYSSES Spectra (ODYSSEUS) program. Here we focus on the HST UV spectra obtained by the HST Director's Discretionary Time UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program. Using accretion shock modeling, we find that all targets exhibit accretion variability, varying from short increases in accretion rate by up to a factor of 3 within 48 hours, to longer decreases in accretion rate by a factor of 2.5 over the course of 1 year. This is despite the generally consistent accretion morphology within each target. Additionally, we test empirical relationships between accretion rate and UV luminosity and find stark differences, showing that these relationships should not be used to estimate the accretion rate for individual target. Our work reinforces that future multi-epoch and simultaneous multi-wavelength studies are critical in our understanding of the accretion process in low-mass star formation.

This paper is one in a series reporting results from small telescope observations of variable young stars. Here, we study the repeating outbursts of three likely Be stars based on long-term optical, near-infrared, and mid-infrared photometry for all three objects, along with follow-up spectra for two of the three. The sources are characterised as rare, truly regularly outbursting Be stars. We interpret the photometric data within a framework for modelling light curve morphology, and find that the models correctly predict the burst shapes, including their larger amplitudes and later peaks towards longer wavelengths. We are thus able to infer the start and end times of mass loading into the circumstellar disks of these stars. The disk sizes are typically 3-6 times the areas of the central star. The disk temperatures are 40 respectively. The available spectroscopy is consistent with inside-out evolution of the disk. Higher excitation lines have larger velocity widths in their double-horned shaped emission profiles. Our observations and analysis support the decretion disk model for outbursting Be stars.

Context. Recent studies have shown that 6.7 GHz methanol maser flares can be a powerful tool for verifying the mechanisms of maser production and even the specific signatures of accretion rate changes in the early stages of high-mass star formation. Aims. We characterize the spatial structure and evolution of methanol and water masers during a flare of methanol maser emission at 6.7 GHz in the high-mass young stellar object (HMYSO) G24.33+0.14. Methods. Very Long Baseline Array (VLBA) was used to image the 6.7 and 12.2 GHz methanol and 22.2 GHz water vapor masers at three epochs guided by monitoring the methanol line with the Torun 32m telescope. The 6.7 GHz maser maps were also obtained with the European VLBI Network (EVN) and Long Baseline Array (LBA) during the flare. The Wide-field Infrared Survey Explorer (WISE) data were used to find correlations between the 6.7 GHz maser and infrared (IR) fluxes. Results. The 6.7 GHz methanol maser cloudlets are distributed over ~3500 au, and the morphology of most of them is stable although their brightness varies following the course of the total flux density on a timescale of two months. The 12.2 GHz methanol maser cloudlets cover an area an order of magnitude smaller than that of 6.7 GHz emission, and both transitions emerge from the same masing gas. The 22.2 GHz maser cloudlets lie in the central region and show a systematic increase in brightness and moderate changes in size and orientation, together with the velocity drift of the strongest cloudlet during two months of the Very Long Baseline Interferometry (VLBI) observing period. Time lag estimates imply the propagation of changes in the physical conditions of the masing region with a subluminal speed (~0.3c). A tight correlation of IR (4.6 μm) and 6.7 GHz flux densities is found, supporting the radiative pumping model. Proper motion analysis does not reveal any signs of expansion or inflow of the methanol cloudlets within ~6 mas over ~10 yr. Comparison with the 230 GHz Atacama Large Millimeter Array (ALMA) data indicates that the methanol masers are distributed in the inner part of the rotating disk, whereas the 22.2 GHz emission traces the compact inner component of the bipolar outflow or a jet structure. Conclusions. The maser morphology in the target is remarkably stable over the course of the flare and is similar to the quiescent state, possibly due to less energetic accretion events that can repeat on a timescale of ~8 yr.

High-mass protostars (M$_{\star} >$ 8 M$_{\odot}$) are thought to gain the majority of their mass via short, intense bursts of growth. This episodic accretion is thought to be facilitated by gravitationally unstable and subsequently inhomogeneous accretion disks. Limitations of observational capabilities, paired with a lack of observed accretion burst events has withheld affirmative confirmation of the association between disk accretion, instability and the accretion burst phenomenon in high-mass protostars. Following its 2019 accretion burst, a heat-wave driven by a burst of radiation propagated outward from the high-mass protostar G358.93-0.03-MM1. Six VLBI (very long baseline interferometry) observations of the raditively pumped 6.7 GHz methanol maser were conducted during this period, tracing ever increasing disk radii as the heat-wave propagated outward. Concatenating the VLBI maps provided a sparsely sampled, milliarcsecond view of the spatio-kinematics of the accretion disk covering a physical range of $\sim$ 50 - 900 AU. We term this observational approach `heat-wave mapping'. We report the discovery of a Keplerian accretion disk with a spatially resolved four-arm spiral pattern around G358.93-0.03-MM1. This result positively implicates disk accretion and spiral arm instabilities into the episodic accretion high-mass star formation paradigm.

The Maser Monitoring Organisation (M2O) is a research community of telescope operators, astronomy researchers and maser theoreticians pursuing a joint goal of reaching a deeper understanding of maser emission and exploring its variety of uses as tracers of astrophysical events. These proceedings detail the origin, motivations and current status of the M2O, as was introduced at the 2021 EVN symposium.

AbstractRecently, remarkable progress has been made in understanding the formation of high mass stars. Observations provided direct evidence that massive young stellar objects (MYSOs), analogously to low-mass ones, form via disk-mediated accretion accompanied by episodic accretion bursts, possibly caused by disk fragmentation. In the case of MYSOs, the mechanism theoretically provides a means to overcome radiation pressure, but in practice it is poorly studied - only three accretion bursts in MYSOs have been caught in action to date. A significant contribution to the development of the theory has been made with the study of masers, which have proven to be a powerful tool for locating “bursting” MYSOs. This overview focuses on the exceptional role that masers play in the search and study of accretion bursts in massive protostars.

AbstractThe Maser Monitoring Organisation is a collection of researchers exploring the use of time-variable maser emission in the investigation of astrophysical phenomena. The forward directed aspects of research primarily involve using maser emission as a tool to investigate star formation. Simultaneously, these activities have deepened knowledge of maser emission itself in addition to uncovering previously unknown maser transitions. Thus a feedback loop is created where both the knowledge of astrophysical phenomena and the utilised tools of investigation themselves are iteratively sharpened. The project goals are open-ended and constantly evolving, however, the reliance on radio observatory maser monitoring campaigns persists as the fundamental enabler of research activities within the group.

The Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director’s Discretionary Program of low-mass pre-main-sequence stars, coupled with forthcoming data from Atacama Large Millimeter/submillimeter Array and James Webb Space Telescope, will provide the foundation to revolutionize our understanding of the relationship between young stars and their protoplanetary disks. A comprehensive evaluation of the physics of disk evolution and planet formation requires understanding the intricate relationships between mass accretion, mass outflow, and disk structure. Here we describe the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey and present initial results of the classical T Tauri Star CVSO 109 in Orion OB1b as a demonstration of the science that will result from the survey. ODYSSEUS will analyze the ULLYSES spectral database, ensuring a uniform and systematic approach in order to (1) measure how the accretion flow depends on the accretion rate and magnetic structures, (2) determine where winds and jets are launched and how mass-loss rates compare with accretion, and (3) establish the influence of FUV radiation on the chemistry of the warm inner regions of planet-forming disks. ODYSSEUS will also acquire and provide contemporaneous observations at X-ray, optical, near-IR, and millimeter wavelengths to enhance the impact of the ULLYSES data. Our goal is to provide a consistent framework to accurately measure the level and evolution of mass accretion in protoplanetary disks, the properties and magnitudes of inner-disk mass loss, and the influence of UV radiation fields that determine ionization levels and drive disk chemistry.

The Orion molecular cloud complex harbours the nearest GMCs and site of high-mass star formation. Its YSO populations are thoroughly characterized. The region is therefore a prime target for the study of star formation. Here, we verify the performance of the SuperCAM 64 pixel heterodyne array on APEX. We give a descriptive overview of a set of wide-field CO(3-2) spectral cubes obtained towards the Orion GMC complex, aimed at characterizing the dynamics and structure of the extended molecular gas in diverse regions of the clouds, ranging from very active sites of clustered star formation in Orion B to comparatively quiet regions in southern Orion A. We present a 2.7 square degree (130pc^2) mapping survey in the CO(3-2) transition, obtained using SuperCAM on APEX at an angular resolution of 19” (7600AU or 0.037pc at a distance of 400pc), covering L1622, NGC2071, NGC2068, OriB9, NGC2024, and NGC2023 in Orion B, and the southern part of the L1641 cloud in Orion A. We describe CO integrated emission and line moment maps and position-velocity diagrams and discuss a few sub-regions in some detail. Evidence for expanding bubbles is seen with lines splitting into double components, most prominently in NGC2024, where we argue that the bulk of the molecular gas is in the foreground of the HII region. High CO(3-2)/CO(1-0) line ratios reveal warm CO along the western edge of Orion B in the NGC2023/NGC2024 region facing the IC434 HII region. Multiple, well separated radial velocity components seen in L1641-S suggest that it consists of a sequence of clouds at increasingly larger distances. We find a small, spherical cloud - the 'Cow Nebula' globule - north of NGC2071. We trace high velocity line wings for the NGC2071-IR outflow and the NGC2024 CO jet. The protostellar dust core FIR4 (rather than FIR5) is the true driving source of the NGC2024 monopolar outflow.

The evolution of young stars and disks is driven by the interplay of several processes, notably the accretion and ejection of material. These processes, critical to correctly describe the conditions of planet formation, are best probed spectroscopically. Between 2020 and 2022, about 500orbits of the Hubble Space Telescope (HST) are being devoted in to the ULLYSES public survey of about 70 low-mass (M⋆ ≤ 2 M⊙) young (age < 10 Myr) stars at UV wavelengths. Here, we present the PENELLOPE Large Program carried out with the ESO Very Large Telescope (VLT) with the aim of acquiring, contemporaneously to the HST, optical ESPRESSO/UVES high-resolution spectra for the purpose of investigating the kinematics of the emitting gas, along with UV-to-NIR X-shooter medium-resolution flux-calibrated spectra to provide the fundamental parameters that HST data alone cannot provide, such as extinction and stellar properties. The data obtained by PENELLOPE have no proprietary time and the fully reduced spectra are being made available to the whole community. Here, we describe the data and the first scientific analysis of the accretion properties for the sample of 13 targets located in the Orion OB1 association and in the σ-Orionis cluster, observed in November–December 2020. We find that the accretion rates are in line with those observed previously in similarly young star-forming regions, with a variability on a timescale of days (≲3). The comparison of the fits to the continuum excess emission obtained with a slab model on the X-shooter spectra and the HST/STIS spectra shows a shortcoming in the X-shooter estimates of ≲10%, which is well within the assumed uncertainty. Its origin can be either due to an erroneous UV extinction curve or to the simplicity of the modeling and, thus, this question will form the basis of the investigation undertaken over the course of the PENELLOPE program. The combined ULLYSES and PENELLOPE data will be key in attaining a better understanding of the accretion and ejection mechanisms in young stars.

Aims. It is well known that low-mass young stellar objects (LMYSOs) gain a significant portion of their final mass through episodes of very rapid accretion, with mass accretion rates up to Ṁ∗ ∼ 10−4 M⊙ yr−1. Recent observations of high-mass young stellar objects (HMYSOs) with masses M∗ ≳ 10 M⊙ uncovered outbursts with accretion rates exceeding Ṁ∗ ∼ 10−3 M⊙ yr−1. Here, we examine which scenarios proposed in the literature so far to explain accretion bursts of LMYSOs can also apply to the episodic accretion in HMYSOs.Methods. We utilise 1D time-dependent models of protoplanetary discs around HMYSOs to study burst properties.Results. We find that discs around HMYSOs are much hotter than those around their low-mass cousins. As a result, a much more extended region of the disc is prone to the thermal hydrogen ionisation and magnetorotational activation instabilities. The former, in particular, is found to be ubiquitous in a very wide range of accretion rates and disc viscosity parameters. The outbursts triggered by these instabilities, however, always have too low of an Ṁ∗ and are one to several orders of magnitude too long compared to those observed from HMYSOs to date. On the other hand, bursts generated by tidal disruptions of gaseous giant planets formed by the gravitational instability of the protoplanetary discs yield properties commensurate with observations, provided that the clumps are in the post-collapse configuration with planet radius Rp ≳ 10 Jupiter radii. Furthermore, if observed bursts are caused by disc ionisation instabilities, then they should be periodic phenomena with the duration of the quiescent phase comparable to that of the bursts. This may yield potentially observable burst periodicity signatures in the jets, the outer disc, or the surrounding diffuse material of massive HMYSOs. Bursts produced by disruptions of planets or more massive objects are not expected to be periodic phenomena, although multiple bursts per protostar are possible.Conclusions. Observations and modelling of episodic accretion bursts across a wide range of young stellar object (YSO) masses is a new promising avenue to break the degeneracy between models of episodic accretion in YSOs.