Episodic accretion plays an important role in the evolution of young stars. Although it has been under investigation for a long time, the origin of such episodic accretion events is not yet understood. We investigate the dust and gas emission of a sample of young outbursting sources in the infrared to get a better understanding of their properties and circumstellar material, and we use the results in a further work to model the objects. We used Herschel data, from our PI program of 12 objects and complemented with archival observations to obtain the spectral energy distributions (SEDs) and spectra of our targets. We report here the main characteristics of our sample, focussing on the SED properties and on the gas emission lines detected in the PACS and SPIRE spectra. The SEDs of our sample show the diversity of the outbursting sources, with several targets showing strong emission in the far-infrared from the embedded objects. Most of our targets reside in a complex environment, which we discuss in detail. We detected several atomic and molecular lines, in particular rotational CO emission from several transitions from J=38-37 to J=4-3. We constructed rotational diagrams for the CO lines, and derived in three domains of assumed local thermodynamic equilibrium (LTE) temperatures and column densities, ranging mainly between 0-100 K and 400-500K. We confirm correlation in our sample between intense CO $J=16-15$ emission and the column density of the warm domain of CO, N(warm). We notice a strong increase in luminosity of HH 381 IRS and a weaker increase for PP 13 S, which shows the beginning of an outburst.

The accretion-driven outbursts of young FU Orionis-type stars may be a common stage of pre-main sequence evolution and can have a significant impact on the circumstellar environment as it pertains to the growth of solids and eventually planets. This episodic accretion is thought to be sustained by additional gas in-falling from the circumstellar envelope and disk. We present APEX observations of the CO gas in the envelope around V883 Orionis, a young outbursting star. The observations mapped the $^{12}$CO(4-3), $^{12}$CO(3-2), and $^{13}$CO(3-2) lines with the FLASH$^{+}$ instrument and the $^{12}$CO(6-5) line with the SEPIA instrument. We detected high signal-to-noise emission extending out to radii $>10000$ au and calculated integrated fluxes of $1100~\rm Jy~km~s^{-1}$ for $^{12}$CO(6-5), $2400~\rm Jy~km~s^{-1}$ for $^{12}$CO(4-3), $1600~\rm Jy~km~s^{-1}$ for $^{12}$CO(3-2), and $450~\rm Jy~km~s^{-1}$ for $^{13}$CO(3-2). We used the thermo-chemical code P{\small RO}D{\small I}M{\small O} to test several models and find the data are best described by an envelope structure with $\rm M_{env}\approx 0.2-0.4\,M_{\odot}$ and a mass-infall rate of $\rm \dot{M}_{inf}=1-2\times10^{-6}\,M_{\odot}\,yr^{-1}$. We infer that the observed envelope and outflow structure around V883 Ori could be caused by multiple outbursts, consistent with episodic accretion.

Grids of stellar models are useful tools to derive the properties of stellar clusters, in particular young clusters hosting massive stars, and to provide information on the star formation process in various mass ranges. Because of their short evolutionary timescale, massive stars end their life while their low-mass siblings are still on the pre-main sequence (pre-MS) phase. Thus the study of young clusters requires consistent consideration of all the phases of stellar evolution. But despite the large number of grids that are available in the literature, a grid accounting for the evolution from the pre-MS accretion phase to the post-MS phase in the whole stellar mass range is still lacking. We build a grid of stellar models at solar metallicity with masses from 0.8 $M_\odot$ to 120 $M_\odot$, including pre-MS phase with accretion. We use the {\sc genec} code to run stellar models on this mass range. The accretion law is chosen to match the observations of pre-MS objects on the Hertzsprung-Russell diagram. We describe the evolutionary tracks and isochrones of our models. The grid is connected to previous MS and post-MS grids computed with the same numerical method and physical assumptions, which provides the widest grid in mass and age to date. Numerical tables of our models and corresponding isochrones are available online.

Episodic accretion may be a common occurrence in the evolution of young pre-main sequence stars and has important implications for our understanding of star and planet formation. Many fundamental aspects of what drives the accretion physics, however, are still unknown. The ngVLA will be a key tool in understanding the nature of these events. The high spatial resolution, broad spectral coverage, and unprecedented sensitivity will allow for the detailed analysis of outburst systems. The proposed frequency range of the ngVLA allows for observations of the gas, dust, and non-thermal emission from the star and disk.

Episodic accretion is an important process in the evolution of young stars and their environment. The observed strong luminosity bursts of young stellar objects likely have a long lasting impact on the chemical evolution of the disk and envelope structure. We want to investigate observational signatures of the chemical evolution in the post-burst phase for embedded sources. With such signatures it is possible to identify targets that experienced a recent luminosity burst. We present a new model for episodic accretion chemistry based on the 2D, radiation thermo-chemical disk code ProDiMo. We have extended ProDiMo with a proper treatment for envelope structures. For a representative Class I model, we calculated the chemical abundances in the post-burst phase and produced synthetic observables like intensity maps and radial profiles. During a burst many chemical species, like CO, sublimate from the dust surfaces. As the burst ends they freeze out again (post-burst phase). This freeze-out happens from inside-out due to the radial density gradient in the disk and envelope structure. This inside-out freeze-out produces clear observational signatures in spectral line emission, like rings and distinct features in the slope of radial intensity profiles. We fitted synthetic C18O J=2-1 observations with single and two component fits and find that post-burst images are much better matched by the latter. Comparing the quality of such fits allows identification of post-burst targets in a model-independent way. Our models confirm that it is possible to identify post-burst objects from spatially resolved CO observations. However, to derive proper statistics, like frequencies of bursts, from observations it is important to consider aspects like the inclination and structure of the target and also dust properties as those have a significant impact on the freeze-out timescale.