Search

The size of a protoplanetary disk is a fundamental property, yet most remain unresolved, even in nearby star-forming regions (d $\sim$ 140-200 pc). We present the complete continuum size distribution for the $105$ brightest protoplanetary disks (M$_{\text{dust}}$ $\gtrsim$ 2 M$_{\oplus}$) in the Ophiuchus cloud, obtained from ALMA Band 8 (410 GHz) observations at 0.05$^{\prime\prime}$ (7 au) to 0.15$^{\prime\prime}$ (21 au) resolution. This sample includes 54 Class II and 51 Class I and Flat Spectrum sources, providing a comprehensive distribution across evolutionary stages. We measure the Half Width at Half Maximum (HWHM) and the radius encircling $68\%$ of the flux ($R_{68\%}$) for most non-binary disks, yielding the largest flux-limited sample of resolved disks in any star-forming region. The distribution is log-normal with a median value of $\sim$14 au and a logarithmic standard deviation $\sigma_{\log} = 0.46$ (factor of 2.9 in linear scale). Disks in close binary systems ($<$ 200 au separation) have smaller radii, with median value of $\sim$5 au, indicating efficient radial drift as predicted by dust evolution models. The size distribution for young embedded objects (SED Class I and Flat Spectrum, age $\lesssim$ 1 Myr) is similar to that of Class II objects (age $\sim$ a few Myr), implying that pressure bumps must be common at early disk stages to prevent mm-sized particle migration at au scales., Comment: 12 pages, 6 figures

The last decade has witnessed remarkable advances in the characterization of the (sub-)millimeter emission from planet-forming disks. Instead, the study of the (sub-)centimeter emission has made more limited progress, to the point that only a few exceptional disk-bearing objects have been characterized in the centimeter regime. This work takes a broad view of the centimeter emission from a large sample with VLA observations that is selected from previous ALMA surveys of more representative disks in brightness and extent. We report on the detection and characterization of flux at centimeter wavelengths from 21 sources in the Taurus star-forming region. Complemented by literature and archival data, the entire photometry from 0.85 mm to 6 cm is fitted by a two-component model that determines the ubiquitous presence of free-free emission entangled with the dust emission. The flux density of the free-free emission is found to scale with the accretion rate but is independent of the outer disk morphology depicted by ALMA. The dust emission at 2 cm is still appreciable, and offers the possibility to extract an unprecedented large set of dust spectral indices in the centimeter regime. A pronounced change between the median millimeter indices (2.3) and centimeter indices (2.8) suggests that a large portion of the disk emission is optically thick up to 3 mm. The comparison of both indices and fluxes with the ALMA disk extent indicates that this portion can be as large as 40 au, and suggests that the grain population within this disk region that emits the observed centimeter emission is similar in disks with different size and morphology. All these results await confirmation and dedicated dust modeling once facilities like ngVLA or SKA-mid are able to resolve the centimeter emission from planet-forming disks and disentangle the various components., Comment: Accepted for publication by A&A

We present a JWST MIRI/MRS spectrum of the inner disk of WISE J044634.16$-$262756.1B (hereafter J0446B), an old ($\sim$34 Myr) M4.5 star but with hints of ongoing accretion. The spectrum is molecule-rich and dominated by hydrocarbons. We detect 14 molecular species (H$_2$, CH$_3$, CH$_4$, C$_2$H$_2$, $^{13}$CCH$_2$, C$_2$H$_4$, C$_2$H$_6$, C$_3$H$_4$, C$_4$H$_2$, C$_6$H$_6$, HCN, HC$_3$N, CO$_2$ and $^{13}$CO$_2$) and 2 atomic lines ([Ne II] and [Ar II]), all observed for the first time in a disk at this age. The detection of spatially unresolved H$_2$ and Ne gas strongly supports that J0446B hosts a long-lived primordial disk, rather than a debris disk. The marginal H$_2$O detection and the high C$_2$H$_2$/CO$_2$ column density ratio indicate that the inner disk of J0446B has a very carbon-rich chemistry, with a gas-phase C/O ratio $\gtrsim$2, consistent with what have been found in most primordial disks around similarly low-mass stars. In the absence of significant outer disk dust substructures, inner disks are expected to first become water-rich due to the rapid inward drift of icy pebbles, and evolve into carbon-rich as outer disk gas flows inward on longer timescales. The faint millimeter emission in such low-mass star disks implies that they may have depleted their outer icy pebble reservoir early and already passed the water-rich phase. Models with pebble drift and volatile transport suggest that maintaining a carbon-rich chemistry for tens of Myr likely requires a slowly evolving disk with $\alpha-$viscosity $\lesssim10^{-4}$. This study represents the first detailed characterization of disk gas at $\sim$30 Myr, strongly motivating further studies into the final stages of disk evolution., Comment: Accepted by ApJL. 8 figures in the main text

We present an investigation of Complex Organic Molecules (COMs) in the spatially resolved Keplerian disk around V883 Ori, an eruptive young star, based on a spectral survey carried out with ALMA in Band 6 (220.7$-$274.9 GHz). We identified about 3,700 molecular emission lines and discovered 23 COMs in the disk. We estimated the column densities of COMs detected through the iterative LTE line fitting method. According to our analyses, using only optically thin lines is critical to deriving the reliable column densities of COMs. Therefore, covering a large frequency range is important for the studies of COMs. The most distinct phenomenon found from the spectra of the V883 Ori disk is that nitrogen-bearing COMs other than CH$_{3}$CN are missing, whereas various oxygen-bearing COMs, except for the CH$_2$OH-bearing molecules, are detected. The missing CH$_2$OH-bearing COMs may indicate the warm water-ice dominant environment for forming COMs. We compared our results with various objects in different evolutionary stages, from Class 0 hot corinos to a Solar System comet 67P/Churyumov-Gerasimenko, to examine the effect of evolution on the COM compositions. In general, the COMs abundances relative to methanol in V883 Ori are higher than in the hot corinos and hot cores, while they are comparable to the cometary values. This may indicate the planet-forming material chemically evolves in the disk midplane after being accreted from the envelope. In addition, as found in the comet 67P/Churyumov-Gerasimenko, nitrogen might also be trapped as ammonium salt within the dust grains in the V883 Ori disk., Comment: Accepted for publication in ApJS

We present a study of the abundance of calcium in the innermost disk of 70 T Tauri stars in the star-forming regions of Chamaeleon I, Lupus and Orion OB1b. We use calcium as a proxy for the refractory material that reaches the inner disk. We used magnetospheric accretion models to analyze the Ca II emission lines and estimate abundances in the accretion flows of the stars, which feed from the inner disks. We find Ca depletion in disks of all three star-forming regions, with 57% of the sample having [Ca/H] < -0.30 relative to the solar abundance. All disks with cavities and/or substructures show depletion, consistent with trapping of refractories in pressure bumps. Significant Ca depletion ([Ca/H] < -0.30) is also measured in 60% of full disks, although some of those disks may have hidden substructures or cavities. We find no correlation between Ca abundance and stellar or disk parameters except for the mass accretion rate onto the star. This could suggest that the inner and outer disks are decoupled, and that the mass accretion rate is related to a mass reservoir in the inner disk, while refractory depletion reflects phenomena in the outer disk related to the presence of structure and forming planets. Our results of refractory depletion and timescales for depletion are qualitatively consistent with expectations of dust growth and radial drift including partitioning of elements and constitute direct evidence that radial drift of solids locked in pebbles takes place in disks., Comment: Accepted for publication in The Astrophysical Journal

The eruptive YSO FU Ori is the eponym of its variable class. FU Ori stars are known to undergo outbursts with amplitudes of $>4$ magnitudes in the $V$ band and durations of several decades. Interaction with a binary companion is one proposed outburst trigger, so understanding both components of the FU Ori system is crucial. A recent HST/STIS observation of the FU Ori system clearly resolves its North and South components. We report here on the spectrum of FU Ori South. We detect NUV continuum emission but no FUV continuum, although several bright emission lines consistent with those seen in T Tauri stars are present. The presence of the C II] 2325 multiplet and many H$_2$ lines indicate active accretion. We estimate the extinction to the source and find that the UV spectrum favors $A_V < 4$, contrary to past estimates based on the NIR spectrum., Comment: 3 pages, 1 figure, accepted by RNAAS

FU Ori objects are the most extreme eruptive young stars known. Their 4 to 5 magnitude photometric outbursts last for decades and are attributed to a factor of up to 10,000 increase in the stellar accretion rate. The nature of the accretion disk-to-star interface in FU Ori objects has remained a mystery for decades. To date, attempts to directly observe a shock or boundary layer have been thwarted by the apparent lack of emission in excess of the accretion disk photosphere down to $\lambda = 2300$ \AA. We present a new NUV and the first high-sensitivity FUV spectrum of FU Ori. The FUV continuum is detected for the first time and, at $\lambda = 1400$ \AA, is more than $10^4$ times brighter than predicted by a viscous accretion disk. We interpret the excess as arising from a shock at the boundary between the disk and the stellar surface. We model the shock emission as a blackbody and find that the temperature of the shocked material is $T_\mathrm{FUV} \approx 16,000 \pm 2000$ K. The shock temperature corresponds to an accretion flow along the surface of the disk that reaches a velocity of 40 km s$^{-1}$ at the boundary, consistent with predictions from simulations., Comment: 16 pages, 8 figures, accepted to ApJL

We present a comprehensive analysis of young stellar object (YSO) variability within the M17 Southwest Extension (M17 SWex), using 3.5 years of monitoring data from the JCMT Transient Survey at sub-millimeter (sub-mm) and 9 years from the NEOWISE mission at mid-infrared (mid-IR). Our study encompasses observations of 147 bright sub-mm peaks identified within our deep JCMT co-added map as well as 156 YSOs in NEOWISE W1 and 179 in W2 that were previously identified in Spitzer surveys. We find three robust sub-mm variables: two are candidate YSOs and one is a likely extragalactic source. At mid-IR wavelengths, our analysis reveals secular and stochastic variability in 47 YSOs, with the highest fraction of secular variability occurring at the earliest evolutionary stage. This is similar to what has previously been observed for low-mass YSO variability within the Gould Belt. However, we observe less overall variability in M17SWex at both the sub-mm and mid-IR. We suspect that this lower fraction is due to the greater distance to M17 SWex. Our findings showcase the utility of multi-wavelength observations to better capture the complex variability phenomena inherent to star formation processes and demonstrate the importance of years-long monitoring of a diverse selection of star-forming environments., Comment: Accepted for Publication in The Astronomical Journal

Constraints on the binary fraction of young massive stellar objects (mYSOs) are important for binary and massive star formation theory. Here, we present speckle imaging of 34 mYSOs located in the Large (1/2 $Z_{\odot}$) and Small Magellanic Clouds ($\sim$1/5 $Z_{\odot}$), probing projected separations between the 2000-20000 au (at angular scales of 0.02-0.2") range, for stars above 8 $M_{\odot}$. We find two wide binaries in the Large Magellanic Cloud (from a sample of 23 targets), but none in a sample of 11 in the Small Magellanic Cloud, leading us to adopt a wide binary fraction of 9$\pm$5%, and $<$5%, respectively. We rule out a wide binary fraction greater than 35% in the Large, and 38% in the Small Magellanic Cloud at the 99% confidence level. This is in contrast to the wide binary fraction of mYSOs in the Milky Way (presumed $Z_{\odot}$), which within the physical parameter space probed by this study is $\sim$15-60% from the literature. We argue that while selection effects could be responsible for the lower binary fraction observed; it is more likely that there are underlying physical mechanisms responsible for the observed properties. This indicates that metallicity and environmental effects may influence the formation of wide binaries among massive stars. Future larger, statistically more significant samples of high-mass systems in low-metallicity environments, and for comparison in the Milky Way, are essential to confirm or repudiate our claim., Comment: Accepted for publication in ApJ

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., Comment: 37 pages, 14 figures

An unbiased spectral survey of V883 Ori, an eruptive young star, was carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 6. The detected line emission from various molecules reveals morphological/kinematical features in both the Keplerian disk and the infalling envelope. A direct infall signature, red-shifted absorption against continuum, has been detected in CO, HCO$^+$, HCN, HNC, and H$_2$CO. HCO$^+$ and SO show large arm-like structures that probably connect from an infalling envelope to the disk. HCN and H$_2$CO reveal a distinct boundary between the inner and outer disk and reveal tentative spiral structures connecting the outer disk to the inner disk. HNC shows a large central emission hole (r $\sim$0.3\arcsec) due to its chemical conversion to HCN at high temperatures. The HDO emission, a direct tracer of the water sublimation region, has been detected in the disk. Molecular emission from complex organic molecules (COMs) is confined within the HDO emission boundary, and HCO$^+$ has an emission hole in its distribution due to its destruction by water. Together, these features suggest that the current episode of eruption in V883 Ori may be triggered by the infall from the envelope to the outer disk, generating a spiral wave, which propagates inward and greatly enhances the accretion onto the central star., Comment: Accepted for publication in ApJ

High-angular resolution imaging by ALMA has revealed the near-universality and diversity of substructures in protoplanetary disks. However, disks around M-type pre-main-sequence stars are still poorly sampled, despite the prevalence of M-dwarfs in the galaxy. Here we present high-resolution (~50 mas, 8 au) ALMA Band 6 observations of six disks around mid-M stars in Taurus. We detect dust continuum emission in all six disks, 12CO in five disks, and 13CO line in two disks. The size ratios between gas and dust disks range from 1.6 to 5.1. The ratio of about 5 for 2M0436 and 2M0450 indicates efficient dust radial drift. Four disks show rings and cavities and two disks are smooth. The cavity sizes occupy a wide range: 60 au for 2M0412, and ~10 au for 2M0434, 2M0436 and 2M0508. Detailed visibility modeling indicates that small cavities of 1.7 and 5.7 au may hide in the two smooth disks 2M0450 and CIDA 12. We perform radiative transfer fitting of the infrared SEDs to constrain the cavity sizes, finding that micron-sized dust grains may have smaller cavities than millimeter grains. Planet-disk interactions are the preferred explanation to produce the large 60 au cavity, while other physics could be responsible for the three ~10 au cavities under current observations and theories. Currently, disks around mid-to-late M stars in Taurus show a higher detection frequency of cavities than earlier type stars, although a more complete sample is needed to evaluate any dependence of substructure on stellar mass., Comment: 25 pages, 6 figures. Accepted for publication in ApJ

We present ${\sim}0.2$ arcsec ($\sim$80 au) resolution observations of the CO (2-1) and SiO (5-4) lines made with the Atacama large millimeter/submillimeter array toward an extremely young intermediate-mass protostellar source (t$_{\rm dyn}<$1000 years), MMS 1 located in the Orion Molecular Cloud-3 region. We have successfully imaged a very compact CO molecular outflow associated with MMS 1, having deprojected lobe sizes of $\sim$18000 au (red-shifted lobe) and $\sim$35000 au (blue-shifted lobe). We have also detected an extremely compact ($\lesssim$1000 au) and collimated SiO protostellar jet within the CO outflow. The maximum deprojected jet speed is measured to be as high as 93 km s$^{-1}$. The SiO jet wiggles and displays a chain of knots. Our detection of the molecular outflow and jet is the first direct evidence that MMS 1 already hosts a protostar. The position-velocity diagram obtained from the SiO emission shows two distinct structures: (i) bow-shocks associated with the tips of the outflow, and (ii) a collimated jet, showing the jet velocities linearly increasing with the distance from the driving source. Comparisons between the observations and numerical simulations quantitatively share similarities such as multiple-mass ejection events within the jet and Hubble-like flow associated with each mass ejection event. Finally, while there is a weak flux decline seen in the 850 $\mu$m light curve obtained with JCMT/SCUBA 2 toward MMS 1, no dramatic flux change events are detected. This suggests that there has not been a clear burst event within the last 8 years., Comment: 19 pages, 9 figures, Accepted for publication in ApJ

We present the outflows detected in HOPS 373SW, a protostar undergoing a modest $30\%$ brightness increase at 850 $\mu$m. Atacama Large Millimeter/submillimeter Array (ALMA) observations of shock tracers, including SiO 8--7, CH$_3$OH 7$_{\rm k}$--6$_{\rm k}$, and $^{12}$CO 3--2 emission, reveal several outflow features around HOPS 373SW. The knots in the extremely high-velocity SiO emission reveal the wiggle of the jet, for which a simple model derives a 37$^\circ$ inclination angle of the jet to the plane of the sky, a jet velocity of 90 km s$^{-1}$, and a period of 50 years. The slow SiO and CH$_3$OH emission traces U-shaped bow shocks surrounding the two CO outflows. One outflow is associated with the high-velocity jets, while the other is observed to be close to the plane of the sky. The misaligned outflows imply that previous episodic accretion events have either reoriented HOPS 373SW or that it is an unresolved protostellar binary system with misaligned outflows., Comment: Accepted for publication in ApJ

The JCMT Transient Survey has been monitoring eight Gould Belt low-mass star-forming regions since December 2015 and six somewhat more distant intermediate-mass star-forming regions since February 2020 with SCUBA-2 on the JCMT at \ShortS and \LongS and with an approximately monthly cadence. We introduce our Pipeline v2 relative calibration procedures for image alignment and flux calibration across epochs, improving on our previous Pipeline v1 by decreasing measurement uncertainties and providing additional robustness. These new techniques work at both \LongS and \ShortNS, where v1 only allowed investigation of the \LongS data. Pipeline v2 achieves better than $0.5^{\prime\prime}$ relative image alignment, less than a tenth of the submillimeter beam widths. The v2 relative flux calibration is found to be 1\% at \LongS and $<5$\% at \ShortNS. The improvement in the calibration is demonstrated by comparing the two pipelines over the first four years of the survey and recovering additional robust variables with v2. Using the full six years of the Gould Belt survey the number of robust variables increases by 50\,\%, and at \ShortS we identify four robust variables, all of which are also robust at \LongNS. The multi-wavelength light curves for these sources are investigated and found to be consistent with the variability being due to dust heating within the envelope in response to accretion luminosity changes from the central source., Comment: Accepted for Publication in the The Astrophysical Journal. DOI link to data will become public after the proof stage is complete

Strong accretion outbursts onto protostars are associated with emission dominated by a viscously heated disk, which is characterized by high luminosities. We report the discovery and characterization of a strong mid-IR (3.4, 4.6 $\mu$m) outburst in the embedded protostar SSTgbs J21470601+4739394 (hereafter SSTgbsJ214706). SSTgbsJ214706 has steadily brightened in the mid-infrared by $\sim2$ magnitudes over the past decade, as observed by NEOWISE. Follow-up investigations with the Gemini near-IR spectrograph reveal that SSTgbsJ214706 is a binary system with a spatially extended outflow. The outburst is occurring on the more embedded southeast (SE) component, which dominates the mid- and far-infrared emission from the source. The outbursting component exhibits a spectrum consistent with an FU Ori-type outburst, including the presence of enhanced absorption observed in the molecular bands of CO. The luminosity of the SE component is estimated to be $\sim 0.23\,$ L$_\odot$ before the outburst and $\sim 0.95\,$ L$_\odot$ during the outburst, which is 1 to 2 orders of magnitude fainter than bonafide FU Ori outbursts. We interpret this eruption as an FU Ori-type outburst, although the possibility of brightening following an extinction episode cannot be ruled out. We discuss the implications and potential explanations for such a low-luminosity eruption.

Protostellar outflows often present a knotty appearance, providing evidence of sporadic accretion in stellar mass growth. To understand the direct relation between mass accretion and ejection, we analyze the contemporaneous accretion activity and associated ejection components in B335. B335 has brightened in the mid-IR by 2.5 mag since 2010, indicating increased luminosity, presumably due to increased mass accretion rate onto the protostar. ALMA observations of 12CO emission in the outflow reveal high-velocity emission, estimated to have been ejected 4.6 - 2 years before the ALMA observation and consistent with the jump in mid-IR brightness. The consistency in timing suggests that the detected high-velocity ejection components are directly linked to the most recent accretion activity. We calculated the kinetic energy, momentum, and force for the ejection component associated with the most recent accretion activity and found that at least, about 1.0% of accreted mass has been ejected. More accurate information on the jet inclination and the temperature of the ejected gas components will better constrain the ejected mass induced by the recently enhanced accretion event., Comment: 11 pages, 5 figures, and 1 table

Variability in the brightness of Young Stellar Objects (YSOs) is a common phenomenon that can be caused by changes in various factors, including accretion, extinction, disk morphology, interactions between the disk and the stellar photosphere, and the rotation of hot or cold magnetic spots on the stellar photosphere. Analyzing the variability on different timescales gives insight into the mechanisms driving the brightness changes in YSOs. We investigate the variability of YSOs on both long and short timescales using two mid-infrared datasets: the NEOWISE 7.5-year W2 (4.6$\mu$m) data and the YSOVAR 40-day Spitzer/IRAC2 (4.5$\mu$m) data, respectively. We classify the variability types in each timescale following Park et al (2021). We find a higher detection rate of variable sources in the short-term (77.6%) compared to the long-term (43.0%) due to the higher sensitivity of the Spitzer observations. In addition, the higher cadence of the YSOVAR data results in the weeks-long short-term variability being mostly secular, while the years-long long-term variability explored with the coarsely sampled NEOWISE data is mostly stochastic. By cross-matching the two catalogs, we statistically analyze the variability types exhibited by YSOs across both timescales. The long-term variability amplitude is mostly three times (up to ten times) greater than the short-term variability. Furthermore, we evaluate variability on very short (1-2 days) timescales and recover a trend of the increasing amplitude of variability as the timescales increase. By comprehensively analyzing the variability of YSOs over various timescales, we contribute to a deeper understanding of the underlying mechanisms driving their variability., Comment: Accepted for publication in ApJ

Lyman-$\alpha$ (Ly$\alpha$) is the strongest emission line in the accretion-generated UV spectra from T-Tauri stars and, as such, plays a critical role in regulating chemistry within the surrounding protoplanetary disks. Due to its resonant nature, the scattering of Ly$\alpha$ photons along the line-of-sight encodes information about the physical properties of the intervening H I medium. In this work, we present the first spatially resolved spectral images of Ly$\alpha$ emission across a protoplanetary disk in the iconic face-on T-Tauri star TW Hya, observed with HST-STIS at spatial offsets 0$''$, $\pm 0.2''$, and $\pm 0.4''$. To comprehensively interpret these Ly$\alpha$ spectra, we utilize a 3D Monte-Carlo Ly$\alpha$ radiative transfer simulation considering the H I wind and protoplanetary disk. From the simulation, we constrain the wind's properties: the H I column density $\sim 10^{20}\, \rm cm^{-2}$ and the outflow velocity $\sim 200\, \rm km\, s^{-1}$. Our findings indicate that successfully interpreting the observed spectra necessitates scattering contributions in the H I layer within the disk. Furthermore, to explore the effect of Ly$\alpha$ radiative transfer on protoplanetary disk chemistry, we compute the radiation field within the scattering medium and reveal that the wind reflection causes more Ly$\alpha$ photons to penetrate the disk. Our results show the necessity of spatially resolved Ly$\alpha$ observations of a broad range of targets, which will decode the complex interactions between the winds, protoplanetary disks, and surrounding environments., Comment: 17 pages, 18 figures, accepted for publication in MNRAS

We present an analysis of the JWST/MIRI-MRS spectrum of Sz114, an accreting M5 star surrounded by a large dust disk with a gap at ~39au. The spectrum is molecular-rich: we report the detection of water, CO, CO2, HCN, C2H2, and H2. The only identified atomic/ionic transition is from [NeII] at 12.81 micron. A distinct feature of this spectrum is the forest of water lines with the 17.22 micron emission surpassing that of most late M-star disks by an order of magnitude in flux and aligning instead with disks of earlier-type stars. Moreover, flux rations of C2H2/H2O and HCN/H2O in Sz114 also resemble those of earlier-type disks, with a slightly elevated CO2/H2O ratio. While accretion heating can boost all infrared lines, the unusual properties of Sz114 could be explained by the young age of the source, its formation under unusual initial conditions (a large massive disk), and/or the presence of dust substructures. The latter delay the inward drift of icy pebbles and thus preserve a lower C/O ratio over an extended period. In contrast, late M-star disks-which are typically faint, small in size, and likely lack significant substructures-may have more quickly depleted the outer icy reservoir and already evolved out of a water-rich inner disk phase. Our findings underscore the unexpected diversity within mid-infrared spectra of late M-star disks, highlighting the need to expand the observational sample for a comprehensive understanding of their variations and thoroughly test pebble drift and planet formation models., Comment: 15 pages, 7 figures, submitted to ApJL

While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data. To characterize the asymmetric structure, a parametric model is used to fit the observed visibilities. We report a tentative azimuthal shift of the dust emission peaks between the observations at the two wavelengths. This shift is consistent with a dust trap caused by a vortex, which may be formed by an embedded protoplanet or other hydrodynamical instabilities, such as a dead zone. Deep high-spatial observations of dust and molecular gas are needed to constrain the mechanisms that formed the observed millimeter cavity and dust asymmetry in the protoplanetary disk around CIDA 9A., Comment: Accepted for publication in ApJ, 12 pages, 7 figures with appendix

The proposed protoplanet AB Aur b is a spatially concentrated emission source imaged in the mm-wavelength disk gap of the Herbig Ae/Be star AB Aur. Its near-infrared spectrum and absence of strong polarized light have been interpreted as evidence supporting the protoplanet interpretation. However, the complex scattered light structures in the AB Aur disk pose challenges in resolving the emission source and interpreting the true nature of AB Aur b. We present new images of the AB Aur system obtained using the Hubble Space Telescope Wide Field Camera 3 in the ultraviolet (UV) and optical bands. AB Aur b and the known disk spirals are recovered in the F336W, F410M, and F645N bands. The spectral energy distribution of AB Aur b shows absorption in the Balmer jump, mimicking those of early-type stars. By comparing the colors of AB Aur b to those of the host star, the disk spirals, and predictions from scattered light and self-luminous models, we find that the emission from AB Aur b is inconsistent with planetary photospheric or accretion shock models. Instead, it is consistent with those measured in the circumstellar disks that trace scattered light. We conclude that the UV and visible emission from AB Aur b does not necessitate the presence of a protoplanet. We synthesize observational constraints on AB Aur b and discuss inconsistent interpretations of AB Aur b among different datasets. Considering the significance of the AB Aur b discovery, we advocate for further observational evidence to verify its planetary nature., Comment: 11 pages, 6 figures. Accepted for publication in AJ

Accretion plays a central role in the physics that governs the evolution and dispersal of protoplanetary disks. The primary goal of this paper is to analyze the stability over time of the mass accretion rate onto TW Hya, the nearest accreting solar-mass young star. We measure veiling across the optical spectrum in 1169 archival high-resolution spectra of TW Hya, obtained from 1998--2022. The veiling is then converted to accretion rate using 26 flux-calibrated spectra that cover the Balmer jump. The accretion rate measured from the excess continuum has an average of $2.51\times10^{-9}$~M$_\odot$~yr$^{-1}$ and a Gaussian distribution with a FWHM of 0.22 dex. This accretion rate may be underestimated by a factor of up to 1.5 because of uncertainty in the bolometric correction and another factor of 1.7 because of excluding the fraction of accretion energy that escapes in lines, especially Ly$\alpha$. The accretion luminosities are well correlated with He line luminosities but poorly correlated with H$\alpha$ and H$\beta$ luminosity. The accretion rate is always flickering over hours but on longer timescales has been stable over 25 years. This level of variability is consistent with previous measurements for most, but not all, accreting young stars., Comment: Accepted by ApJ. 31 pages

EX Lup is the archetype for the class of young stars that undergoes repeated accretion outbursts of $\sim 5$ mag at optical wavelengths and that last for months. Despite extensive monitoring that dates back 130 years, the accretion history of EX Lup remains mostly qualitative and has large uncertainties. We assess historical accretion rates of EX Lup by applying correlations between optical brightness and accretion, developed on multi-band magnitude photometry of the $\sim 2$ mag optical burst in 2022. Two distinct classes of bursts occur: major outbursts ($\Delta V\sim5$ mag) have year-long durations, are rare, reach accretion rates of $\dot{M}_{\rm acc}\sim10^{-7}~M_\odot~{\rm yr^{-1}}$ at peak, and have a total accreted mass of around 0.1 Earth masses. The characteristic bursts ($\Delta V\sim2$ mag) have durations of $\sim 2-3$ months, are more common, reach accretion rates of $\dot{M}_{\rm acc}\sim10^{-8}~M_\odot~{\rm yr^{-1}}$ at peak, and have a total accreted mass of around $10^{-3}$ Earth masses. The distribution of total accreted mass in the full set of bursts is poorly described by a power law, which suggests different driving causes behind the major outburst and characteristic bursts. The total mass accreted during two classes of bursts is around two times the masses accreted during quiescence. Our analysis of the light curves reveals a color-dependent time lag in the 2022 post-burst light curve, attributed to the presence of both hot and cool spots on the stellar surface., Comment: Accepted by ApJ; 29 pages,19 figures, 6 tables

Previous analyses of mid-infrared water spectra from young protoplanetary disks observed with the Spitzer-IRS found an anti-correlation between water luminosity and the millimeter dust disk radius observed with ALMA. This trend was suggested to be evidence for a fundamental process of inner disk water enrichment, used to explain properties of the Solar System 40 years ago, in which icy pebbles drift inward from the outer disk and sublimate after crossing the snowline. Previous analyses of IRS water spectra, however, were uncertain due to the low spectral resolution that blended lines together. We present new JWST-MIRI spectra of four disks, two compact and two large with multiple radial gaps, selected to test the scenario that water vapor inside the snowline is regulated by pebble drift. The higher spectral resolving power of MIRI-MRS now yields water spectra that separate individual lines, tracing upper level energies from 900 K to 10,000 K. These spectra clearly reveal excess emission in the low-energy lines in compact disks, compared to the large disks, demonstrating an enhanced cool component with $T \approx$ 170-400 K and equivalent emitting radius $R_{\rm{eq}}\approx$ 1-10 au. We interpret the cool water emission as ice sublimation and vapor diffusion near the snowline, suggesting that there is indeed a higher inwards mass flux of icy pebbles in compact disks. Observation of this process opens up multiple exciting prospects to study planet formation chemistry in inner disks with JWST., Comment: Accepted for publication in ApJ Letters

The spatial distribution and evolution of gas in the inner 10 au of protoplanetary disks form the basis for estimating the initial conditions of planet formation. Among the most important constraints derived from spectroscopic observations of the inner disk are the radial distributions of the major gas phase constituents, how the properties of the gas change with inner disk dust evolution, and how chemical abundances and excitation conditions are influenced by the high-energy radiation from the central star. We present a survey of the radial distribution, excitation, and evolution of inner disk molecular hydrogen (H$_{2}$) obtained as part of the $HST$/ULLYSES program. We analyze far-ultraviolet spectroscopy of 71 (63 accreting) pre-main sequence systems in the ULLYSES DR5 release to characterize the H$_{2}$ emission lines, H$_{2}$ dissociation continuum emission, and major photochemical/disk evolution driving UV emissions (Ly$\alpha$, UV continuum, and C IV). We use the widths of the H$_{2}$ emission lines to show that most fluorescent H$_{2}$ arises between 0.1 - 1.4 au from the parent star, and show positive correlations of the average emitting radius with the accretion luminosity and with the dust disk mass. We find a strong correlation between H$_{2}$ dissociation emission and both the accretion-dominated Ly$\alpha$ luminosity and the inner disk dust clearing, painting a picture where water molecules in the inner 3 au are exposed to and dissociated by strong Ly$\alpha$ emission as the opacity of the inner disk declines with time., Comment: 19 pages, 11 figures. Accepted to the Astronomical Journal

Substructures in protoplanetary disks can act as dust traps that shape the radial distribution of pebbles. By blocking the passage of pebbles, the presence of gaps in disks may have a profound effect on pebble delivery into the inner disk, crucial for the formation of inner planets via pebble accretion. This process can also affect the delivery of volatiles (such as H$_2$O) and their abundance within the water snow line region (within a few au). In this study, we aim to understand what effect the presence of gaps in the outer gas disk may have on water vapor enrichment in the inner disk. Building on previous work, we employ a volatile-inclusive multi-Myr disk evolution model that considers an evolving ice-bearing drifting dust population, sensitive to dust-traps, which loses its icy content to sublimation upon reaching the snow line. We find that vapor abundance in the inner disk is strongly affected by fragmentation velocity (v$_{\rm f}$) and turbulence, which control how intense vapor enrichment from pebble delivery is, if present, and how long it may last. Generally, for disks with low to moderate turbulence ($\alpha$ $\le$ 1 $\times$ 10$^{-3}$) and for a range of v$_{\rm f}$, radial location, and gap depth (especially that of the innermost gaps), can significantly alter enrichment. Shallow inner gaps may continuously leak material from beyond it, despite the presence of additional deep outer gaps. We finally find that the for realistic v$_{\rm f}$ ($\le$ 10 m s$^{-1}$), presence of gaps is more important than planetesimal formation beyond the snow line in regulating pebble and volatile delivery into the inner disk., Comment: 24 Pages, 15 Figures, Accepted to ApJ

We present the spectra of Complex Organic Molecules (COMs) detected in HOPS 373SW with the Atacama Large Millimeter/submillimeter Array (ALMA). HOPS 373SW, which is a component of a protostellar binary with a separation of 1500 au, has been discovered as a variable protostar by the JCMT Transient monitoring survey with a modest ~30% brightness increase at submillimeter wavelengths. Our ALMA Target of Opportunity (ToO) observation at ~345 GHz for HOPS 373SW revealed extremely young chemical characteristics with strong deuteration of methanol. The dust continuum opacity is very high toward the source center, obscuring line emission from within 0.03 arcsec. The other binary component, HOPS 373NE, was detected only in C17O in our observation, implying a cold and quiescent environment. We compare the COMs abundances relative to CH3OH in HOPS 373SW with those of V883 Ori, which is an eruptive disk object, as well as other hot corinos, to demonstrate the chemical evolution from envelope to disk. High abundances of singly, doubly, and triply deuterated methanol (CH2DOH, CHD2OH, and CD3OH) and a low CH3CN abundance in HOPS 373SW compared to other hot corinos suggest a very early evolutionary stage of HOPS 373SW in the hot corino phase. Since the COMs detected in HOPS 373SW would have been sublimated very recently from grain surfaces, HOPS 373SW is a promising place to study the surface chemistry of COMs in the cold prestellar phase, before sublimation., Comment: Accepted for publication in ApJ

The effects of metallicity on the evolution of protoplanetary disks may be studied in the outer Galaxy where the metallicity is lower than in the solar neighbourhood. We present the VLT/KMOS integral field spectroscopy in the near-infrared of $\sim$120 candidate young stellar objects (YSOs) in the CMa-$\ell$224 star-forming region located at a Galactocentric distance of 9.1 kpc. We characterise the YSO accretion luminosities and accretion rates using the hydrogen Br$\gamma$ emission and find the median accretion luminosity of $\log{(L_{\rm acc})} = -0.82^{+0.80}_{-0.82} L_\odot$. Based on the measured accretion luminosities, we investigate the hypothesis of star formation history in the CMa-$\ell$224. Their median values suggest that Cluster C, where most of YSO candidates have been identified, might be the most evolved part of the region. The accretion luminosities are similar to those observed toward low-mass YSOs in the Perseus and Orion molecular clouds, and do not reveal the impact of lower metallicity. Similar studies in other outer Galaxy clouds covering a wide range of metallicities are critical to gain a complete picture of star formation in the Galaxy., Comment: Accepted for publication in APJS, 51 pages, 37 figures, 6 tables

Episodic accretion is a low-mass pre-main sequence phenomenon characterized by sudden outbursts of enhanced accretion. These objects are classified into two: protostars with elevated levels of accretion that lasts for decades or more, called FUors, and protostars with shorter and repetitive bursts, called EXors. HBC 494 is a FUor object embedded in the Orion Molecular Cloud. Earlier Atacama Large (sub-)Millimeter Array (ALMA) continuum observations showed an asymmetry in the disk at 0.''2 resolution. Here, we present follow-up observations at ~0.''03, resolving the system into two components: HBC 494 N (primary) and HBC 494 S (secondary). No circumbinary disk was detected. Both disks are resolved with a projected separation of ~0.''18 (75 au). Their projected dimensions are 84+/-1.8 x 66.9+/-1.5 mas for HBC 494 N and 64.6+/-2.5 x 46.0+/-1.9 mas for HBC 494 S. The disks are almost aligned and with similar inclinations. The observations show that the primary is ~5 times brighter/more massive and ~2 times bigger than the secondary. We notice that the northern component has a similar mass to the FUors, while the southern has to EXors. The HBC 494 disks show individual sizes that are smaller than single eruptive YSOs. In this work, we also report 12CO, 13CO, and C18O molecular line observations. At large scale, the 12CO emission shows bipolar outflows, while the 13CO and C18O maps show a rotating and infalling envelope around the system. At a smaller scale, the 12CO and 13CO moment zero maps show cavities within the continuum disks' area, which may indicate continuum over-subtraction or slow-moving jets and chemical destruction along the line-of-sight., Comment: Published in MNRAS. Main text (11 pages, 8 figures). Appendix (11 pages, 9 figures). Total (22 pages, 17 figures)

Disk winds are thought to play a critical role in the evolution and dispersal of protoplanetary disks. A primary diagnostic of this physics is emission from the wind, especially in the low-velocity component of the [O I] $\lambda6300$ line. However, the interpretation of the line is usually based on spectroscopy alone, which leads to confusion between magnetohydrodynamic winds and photoevaporative winds. Here we report that in high-resolution spectral mapping of TW Hya by the multi-unit spectroscopic explorer at the Very Large Telescope, 80% of the [O I] emission is confined to within 1 AU radially from the star. A generic model of a magnetothermal wind produces [O I] emission at the base of the wind that broadly matches the flux and the observed spatial and spectral profiles. The emission at large radii is much fainter that predicted from models of photoevaporation, perhaps because the magnetothermal wind partially shields the outer disk from energetic radiation from the central star. This result calls into question the previously assessed importance of photoevaporation in disk dispersal predicted by models., Comment: Published in Nature Astronomy

We present an analysis of 288 young stellar objects (YSOs) in the Perseus Molecular Cloud that have well defined $g$ and $r$-band lightcurves from the Zwicky Transient Facility. Of the 288 YSOs, 238 sources (83% of our working sample) are identified as variables based on the normalized peak-to-peak variability metric, with variability fraction of 92% for stars with disks and 77% for the diskless populations. These variables are classified into different categories using the quasiperiodicity ($Q$) and flux asymmetry ($M$) metrics. Fifty-three variables are classified as strictly periodic objects that are well phased and can be attributed to spot modulated stellar rotation. We also identify 22 bursters and 25 dippers, which can be attributed to accretion burst and variable extinction, respectively. YSOs with disks tend to have asymmetric and non-repeatable lightcurves, while the YSOs without disks tend to have (quasi)periodic lightcurves. The periodic variables have the steepest change in $g$ versus $g-r$, while bursters have much flatter changes than dippers in $g$ versus $g-r$. Periodic and quasiperiodic variables display the lowest variability amplitude. Simple models suggest that the variability amplitudes of periodic variables correspond to changes of the spot coverage of 30% to 40%, burster variables are attributed to accretion luminosity changes in the range of $L_{\rm acc}/L_{\star}=0.1-0.3$, and dippers are due to variable extinction with $A_{V}$ changes in the range of $0.5-1.3\;$mag., Comment: 26 pages, 16 figures, 3 tables, accepted for publication in RAA

The process of accretion through circumstellar disks in young stellar objects is an integral part of star formation and the $H\alpha$ emission line is a prominent signature of accretion in low-mass stars. We present the detection and characterization of $H\alpha$ emission line sources in the central region of a distant, low-metallicity young stellar cluster - Dolidze 25 (at $\sim$ 4.5 kpc) - using medium-resolution optical spectra (4750-9350 \r{A} ) obtained with the Multi-Unit Spectroscopic Explorer (MUSE) at the VLT. We have identified 14 potential accreting sources within a rectangular region of (2$'$ x 1$'$) towards the center of the cluster based on the detection of strong and broad emissions in $H\alpha$ as well as the presence of other emission lines such as [OI] and $H\beta$. Based on their positions in both photometric color-magnitude and color-color diagrams, we have also confirmed that these objects belong to the pre-main sequence phase of star formation. Our results were compared with the disk and diskless members of the cluster previously identified by Guarcello et al. (2021) using near-IR colors, and all sources they had identified as disks were confirmed to be accreting based on the spectroscopic characteristics.

Low-mass stars and sub-stellar objects are essential in tracing the initial mass function (IMF). We study the nearby young $\sigma$ Orionis cluster (d$\sim$408 pc; age$\sim$1.8 Myr) using deep NIR photometric data in J, W and H-bands from WIRCam on the Canada-France-Hawaii Telescope. We use the water absorption feature to photometrically select the brown dwarfs and confirm their nature spectroscopically with the IRTF-SpeX. Additionally we select candidate low-mass stars for spectroscopy and analyze their membership and that of literature sources using astrometry from Gaia DR3. We obtain the near-IR spectra for 28 very low-mass stars and brown dwarfs and estimate their spectral type between M3-M8.5 (mass ranging between 0.3-0.01 M$_{\odot}$). Apart from these, we also identify 5 new planetary mass candidates which require further spectroscopic confirmation of youth. We compile the comprehensive catalog of 170 spectroscopically confirmed members in the central region of the cluster, for a wide mass range of $\sim$19-0.004 M$_{\odot}$. We estimate the star/BD ratio to be $\sim$4, within the range reported for other nearby star forming regions. With the updated catalog of members we trace the IMF down to 4 M$_\mathrm{Jup}$ and we find that a two-segment power-law fits the sub-stellar IMF better than the log-normal distribution., Comment: Accepted for publication in The Astrophysical Journal (ApJ). 27 pages, 9 figures, 2 tables

Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcsec ($\sim$70 and 116 au), with Gaussian $\sigma$ widths of 5.6 and 8.5 au, respectively. The width of the outer ring is smaller than the estimated pressure scale height by $\sim25\%$, suggesting dust trapping in a radial pressure bump. The dust disk size, set by the location of the outermost ring, is significantly larger (by $3\sigma$) than other disks with similar millimeter luminosity, which can be explained by an early formation of local pressure bump to stop radial drift of millimeter dust grains. After considering the disk's physical structure and accretion properties, we prefer planet--disk interaction over dead zone or photoevaporation models to explain the observed dust disk morphology. We carry out high-contrast imaging at $L'$ band using Keck/NIRC2 to search for potential young planets, but do not identify any source above $5\sigma$. Within the dust gap between the two rings, we reach a contrast level of $\sim$7 mag, constraining the possible planet below $\sim$2--4 $M_{\rm Jup}$. Analyses of the gap/ring properties suggest a $\sim$Saturn mass planet at $\sim$90 au is likely responsible for the formation of the outer ring, which can be potentially revealed with JWST., Comment: 15 pages, 5 figures. Accepted for publication in ApJ

Mid-infrared (mid-IR) variability in young stellar objects (YSOs) is driven by several physical mechanisms, which produce a variety of amplitudes and light curve shapes. One of these mechanisms, variable disk accretion is predicted by models of episodic accretion to drive secular variability, including in the mid-IR. Because the largest accretion bursts are rare, adding new objects to the YSO eruptive variable class aids our understanding of the episodic accretion phenomenon and its possible impact on stellar and planetary formation. A previous analysis of 6.5 yr of NeoWISE light curves (3-5 $\mu$m) of ~7000 nearby YSOs found an increase in the fraction of variability and variability amplitude for objects at younger stages of evolution. To help interpret these light curves, we have obtained low- and high-resolution near-IR spectra of 78 objects from this sample of YSOs. In this work, we present the analysis of nine nearby YSOs (d$<$1 kpc) that show the characteristics of known classes of eruptive variable YSOs. We find one FUor-like source, one EX Lupi-type object, and six YSOs with mixed characteristics, or V1647 Ori-like objects. The varied characteristics observed in our sample are consistent with recent discoveries of eruptive YSOs. We discuss how a wide range in YSO outburst parameters (central mass, maximum accretion rate during outburst, evolutionary stage and/or instability leading to the outburst) may play a significant role in the observed spectro-photometric properties of YSO outbursts., Comment: 17 pages, 7 figures, accepted for publication in MNRAS on 14 March 2023

T Tauri stars produce broad Lyman-alpha emission lines that contribute $\sim$88% of the total UV flux incident on the inner circumstellar disks. Lyman-alpha photons are generated at the accretion shocks and in the protostellar chromospheres and must travel through accretion flows, winds and jets, the protoplanetary disks, and the interstellar medium before reaching the observer. This trajectory produces asymmetric, double-peaked features that carry kinematic and opacity signatures of the disk environments. To understand the link between the evolution of Lyman-alpha emission lines and the disks themselves, we model HST-COS spectra from targets included in Data Release 3 of the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We find that resonant scattering in a simple spherical expanding shell is able to reproduce the high velocity emission line wings, providing estimates of the average velocities within the bulk intervening H I. The model velocities are significantly correlated with the K band veiling, indicating a turnover from Lyman-alpha profiles absorbed by outflowing winds to emission lines suppressed by accretion flows as the hot inner disk is depleted. Just 30% of targets in our sample have profiles with red-shifted absorption from accretion flows, many of which have resolved dust gaps. At this stage, Lyman-alpha photons may no longer intersect with disk winds along the path to the observer. Our results point to a significant evolution of Lyman-alpha irradiation within the gas disks over time, which may lead to chemical differences that are observable with ALMA and JWST., Comment: accepted to ApJ

Protostars likely accrete material at a highly time variable rate, however, measurements of accretion variability from the youngest protostars are rare, as they are still deeply embedded within their envelopes. Sub-mm/mm observations can trace the thermal response of dust in the envelope to accretion luminosity changes, allowing variations in the accretion rate to be quantified. In this paper, we present contemporaneous sub-mm/mm light curves of variable protostars in Serpens Main, as observed by the ALMA ACA, SMA, and JCMT. The most recent outburst of EC 53 (V371 Ser), an $\sim 18$ month periodic variable, is well-sampled in the SMA and JCMT observations. The SMA light curve of EC 53 is observed to peak weeks earlier and exhibit a stronger amplitude than at the JCMT. Stochastic variations in the ACA observations are detected for SMM 10 IR with a factor $\sim 2$ greater amplitude than as seen by the JCMT. We develop a toy model of the envelope response to accretion outbursts to show EC 53's light curves are plausibly explained by the delay associated with the light travel time across the envelope and the additional dilution of the JCMT response by the incorporation of cold envelope material in the beam. The larger JCMT beam can also wash out the response to rapid variations, which may be occurring for SMM 10 IR. Our work thus provides a valuable proof of concept for the usage of sub-mm/mm observations as a probe of both the underlying accretion luminosity variations and the protostellar environment., Comment: 29 pages, 14 figures. Accepted for publication in the Astrophysical Journal

Short-duration flares at millimeter wavelengths provide unique insights into the strongest magnetic reconnection events in stellar coronae, and combine with longer-term variability to introduce complications to next-generation cosmology surveys. We analyze 5.5 years of JCMT Transient Survey 850 micron submillimeter monitoring observations toward eight Gould Belt star-forming regions to search for evidence of transient events or long-duration variability from faint sources. The eight regions (30 arcmin diameter fields), including ~1200 infrared-selected YSOs, have been observed on average 47 times with integrations of approximately half an hour, or one day total spread over 5.5 years. Within this large data set, only two robust faint source detections are recovered: JW 566 in OMC 2/3 and MGM12 2864 in NGC 2023. JW 566, a Class II TTauri binary system previously identified as an extraordinary submillimeter flare, remains unique, the only clear single-epoch transient detection in this sample with a flare eight times bright than our ~4.5 sigma detection threshold of 55 mJy/beam. The lack of additional recovered flares intermediate between JW 566 and our detection limit is puzzling, if smaller events are more common than larger events. In contrast, the other submillimeter variable identified in our analysis, Source 2864, is highly variable on all observed timescales. Although Source 2864 is occasionally classified as a YSO, the source is most likely a blazar. The degree of variability across the electromagnetic spectrum may be used to aid source classification., Comment: Accepted for publication in the Astrophysicall Journal. 17 pages, 6 figures, 2 tables

The coevolution of T Tauri stars and their surrounding protoplanetary disks dictates the timescales of planet formation. In this paper, we present magnetospheric accretion and inner disk wall model fits to NUV-NIR spectra of nine classical T Tauri stars in Orion OB1b as part of the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey. Using NUV-optical spectra from the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program and optical-NIR spectra from the PENELLOPE VLT Large Programme, we find that the accretion rates of these targets are relatively high for the region's intermediate age of 5.0 Myr; rates range from $0.5-17.2 \times 10^{-8}$ M$_{\odot}$/yr, with a median value of $1.2\times 10^{-8}$ M$_{\odot}$/yr. The NIR excesses can be fit with 1200-1800 K inner disk walls located at 0.05-0.10 AU from the host stars. We discuss the significance of the choice in extinction law, as the measured accretion rate depends strongly on the adopted extinction value. This analysis will be extended to the complete sample of T Tauri stars being observed through ULLYSES to characterize accretion and inner disks in star-forming regions of different ages and stellar populations., Comment: 22 pages, 3 figures, 8 tables. Accepted for publication in AJ

Identifying the young optically visible population in a star-forming region is essential for fully understanding the star formation event. In this paper, We identify 211 candidate members of the Perseus molecular cloud based on Gaia astronomy. We use LAMOST spectra to confirm that 51 of these candidates are new members, bringing the total census of known members to 856. The newly confirmed members are less extincted than previously known members. Two new stellar aggregates are identified in our updated census. With the updated member list, we obtain a statistically significant distance gradient of $\rm 4.84\;pc\;deg^{-1}$ from west to east. Distances and extinction corrected color-magnitude diagrams indicate that NGC 1333 is significantly younger than IC 348 and the remaining cloud regions. The disk fraction in NGC 1333 is higher than elsewhere, consistent with its youngest age. The star formation scenario in the Perseus molecular cloud is investigated and the bulk motion of the distributed population is consistent with the cloud being swept away by the Per-Tau Shell., Comment: 27 pages, 14 figures, 1 table; Accepted for publication in ApJ

FU Ori type objects (FUors) are decades-long outbursts of accretion onto young stars that are strong enough to viscously heat disks so that the disk outshines the central star. We construct models for FUor objects by calculating emission components from a steady-state viscous accretion disk, a passively-heated dusty disk, magnetospheric accretion columns, and the stellar photosphere. We explore the parameter space of the accretion rate and stellar mass to investigate implications on the optical and near-infrared spectral energy distribution and spectral lines. The models are validated by fitting to multi-wavelength photometry of three confirmed FUor objects, FU Ori, V883 Ori and HBC 722 and then comparing the predicted spectrum to observed optical and infrared spectra. The brightness ratio between the viscous disk and the stellar photosphere, $\eta$, provides an important guide for identifying viscous accretion disks, with $\eta=1$ ("transition line") and $\eta=5$ ("sufficient dominance line") marking turning points in diagnostics, evaluated here in the near-infrared. These turning points indicate the emergence and complete development of FUor-characteristic strong CO absorption, weak metallic absorption, the triangular spectral continuum shape in the $H$-band, and location in color-magnitude diagrams. Lower stellar mass $M_*$ and higher accretion rate $\dot{M}$ lead to larger $\eta$; for $M_*=0.3~{\rm M_\odot}$, $\eta=1$ corresponds to $\dot{M}=2\times10^{-7}~{\rm M_\odot}/$yr and $\eta=5$ to $\dot{M}=6\times10^{-7}~{\rm M_\odot}/$yr. The sufficient dominance line also coincides with the expected accretion rate where accreting material directly reaches the star. We discuss implications of the models on extinction diagnostics, FUor brightening timescales, viscous disks during initial protostellar growth, and eruptive young stellar objects (YSO) associated with FUors., Comment: 32 pages, 23 figures, accepted for publication in ApJ

Molecular hydrogen is the most abundant molecule in the Galaxy and plays important roles for planets, their circumstellar environments, and many of their host stars. We have confirmed the presence of molecular hydrogen in the AU Mic system using high-resolution FUV spectra from HST-STIS during both quiescence and a flare. AU Mic is a $\sim$23 Myr M dwarf which hosts a debris disk and at least two planets. We estimate the temperature of the gas at 1000 to 2000 K, consistent with previous detections. Based on the radial velocities and widths of the H$_2$ line profiles and the response of the H$_2$ lines to a stellar flare, the H$_2$ line emission is likely produced in the star, rather than in the disk or the planet. However, the temperature of this gas is significantly below the temperature of the photosphere ($\sim$3650 K) and the predicted temperature of its star spots ($\gtrsim$2650 K). We discuss the possibility of colder star spots or a cold layer in the photosphere of a pre-main sequence M dwarf., Comment: accepted to ApJ, 20 pages, many figures

The Serpens Molecular Cloud is one of the most active star-forming regions within 500 pc, with over one thousand of YSOs at different evolutionary stages. The ages of the member stars inform us about the star formation history of the cloud. In this paper, we develop a spectral energy distribution (SED) fitting method for nearby evolved (diskless) young stars from members of the Pleiades to estimate their ages, with a temperature scale adopted from APOGEE spectra. When compared with literature temperatures of selected YSOs in Orion, the SED fits to cool (<5000 K) stars have temperatures that differ by an average of <~ 50 K and have a scatter of ~ 210 K for both disk-hosting and diskless stars. We then apply this method to YSOs in the Serpens Molecular Cloud to estimate ages of optical members previously identified from Gaia DR2 astrometry data. The optical members in Serpens are concentrated in different subgroups with ages from ~4 Myr to ~22 Myr; the youngest clusters, W40 and Serpens South, are dusty regions that lack enough optical members to be included in this analysis. These ages establish that the Serpens Molecular Cloud has been forming stars for much longer than has been inferred from infrared surveys., Comment: 28 pages, 22 figures, accepted for publication in The Astrophysical Journal

Variable accretion in young stellar objects reveals itself photometrically and spectroscopically over a continuum of timescales and amplitudes. Most dramatic are the large outbursts (e.g., FU Ori, V1647 Ori, and EX Lup type events), but more frequent are the less coherent, smaller burst-like variations in accretion rate. Improving our understanding of time-variable accretion directly addresses the fundamental question of how stars gain their masses. We review variability phenomena, as characterized from observations across the wavelength spectrum, and how those observations probe underlying physical conditions. The diversity of observed lightcurves and spectra at optical and infrared wavelengths defies a simple classification of outbursts and bursts into well-defined categories. Mid-infrared and submillimeter wavelengths are sensitive to lower-temperature phenomena and more embedded, younger sources, and it is currently unclear if observed flux variations probe similar or distinct physics relative to the shorter wavelengths. We highlight unresolved issues and emphasize the value of spectroscopy, multiwavelength studies, and ultimately patience in using variable accretion to understand stellar mass assembly., Comment: 22 pages, 7 figures. Review chapter for Protostars and Planets VII. Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Final accepted version

Observed changes in protostellar brightness can be complicated to interpret. In our JCMT~Transient monitoring survey, we discovered that a young binary protostar, HOPS 373, is undergoing a modest $30\%$ brightness increase at 850 $\mu$m, caused by a factor of 1.8$-$3.3 enhancement in the accretion rate. The initial burst occurred over a few months, with a sharp rise and then shallower decay. A second rise occurred soon after the decay, and the source is still bright one year later. The mid-IR emission, the small-scale CO outflow mapped with ALMA, and the location of variable maser emission indicate that the variability is associated with the SW component. The near-infrared and NEOWISE $W1$ and $W2$ emission is located along the blueshifted CO outflow, spatially offset by $\sim3$ to $4^{\prime\prime}$ from the SW component. The $K$-band emission imaged by UKIRT shows a compact H$_2$ emission source at the edge of the outflow, with a tail tracing the outflow back to the source. The $W1$ emission, likely dominated by scattered light, brightens by 0.7 mag, consistent with expectations based on the sub-mm lightcurve. The signal of continuum variability in $K$-band and $W2$ is masked by stable H$_2$ emission, as seen in our Gemini/GNIRS spectrum, and perhaps by CO emission. These differences in emission sources complicate infrared searches for variability of the youngest protostars., Comment: 29 pages, 14 figures, accepted for publication in ApJ

Low mass star formation inside massive clusters is crucial to understand the effect of cluster environment on processes like circumstellar disk evolution, planet and brown dwarf formation. The young massive association of Cygnus OB2, with a strong feedback from massive stars, is an ideal target to study the effect of extreme environmental conditions on its extensive low-mass population. We aim to perform deep multi-wavelength studies to understand the role of stellar feedback on the IMF, brown dwarf fraction and circumstellar disk properties in the region. We introduce here, the deepest and widest optical photometry of 1.5$^\circ$ diameter region centred at Cygnus OB2 in r$_{2}$, i$_{2}$, z and Y-filters using Subaru Hyper Suprime-Cam (HSC). This work presents the data reduction, source catalog generation, data quality checks and preliminary results about the pre-main sequence sources. We obtain 713,529 sources in total, with detection down to $\sim$ 28 mag, 27 mag, 25.5 mag and 24.5 mag in r$_{2}$, i$_{2}$, z and Y-band respectively, which is $\sim$ 3 - 5 mag deeper than the existing Pan-STARRS and GTC/OSIRIS photometry. We confirm the presence of a distinct pre-main sequence branch by statistical field subtraction of the central 18$^\prime$ region. We find the median age of the region as $\sim$ 5 $\pm$ 2 Myrs with an average disk fraction of $\sim$ 9$\%$. At this age, combined with A$_V$ $\sim$ 6 - 8 mag, we detect sources down to a mass range $\sim$ 0.01 - 0.17 M$_\odot$. The deep HSC catalog will serve as the groundwork for further studies on this prominent active young cluster., Comment: Accepted for publication in Monthly Notices of Royal Astronomical Society (MNRAS), 20 pages, 21 figures, 4 Tables, typos corrected

Using HST-COS FUV spectra, we have discovered warm molecular hydrogen in the TWA 7 system. TWA 7, a $\sim$9 Myr old M2.5 star, has a cold debris disk and has previously shown no signs of accretion. Molecular hydrogen is expected to be extremely rare in a debris disk. While molecular hydrogen can be produced in star spots or the lower chromospheres of cool stars such as TWA 7, fluxes from progressions that get pumped by the wings of Ly$\alpha$ indicate that this molecular hydrogen could be circumstellar and thus that TWA 7 is accreting at very low levels and may retain a reservoir of gas in the near circumstellar environment., Comment: accepted to ApJ, 18 pages, 9 figures

We present an analysis of wind absorption in the C II ${\lambda}1335$ doublet towards 40 classical T Tauri stars with archival far-ultraviolet (FUV) spectra obtained by the Hubble Space Telescope. Absorption features produced by fast or slow winds are commonly detected (36 out of 40 targets) in our sample. The wind velocity of the fast wind decreases with disk inclination, consistent with expectations for a collimated jet. Slow wind absorption is detected mostly in disks with intermediate or high inclination, without a significant dependence of wind velocity on disk inclination. Both the fast and slow wind absorption are preferentially detected in FUV lines of neutral or singly ionized atoms. The Mg II ${\lambda}{\lambda}2796,2804$ lines show wind absorption consistent with the absorption in the C II lines. We develop simplified semi-analytical disk/wind models to interpret the observational disk wind absorption. Both fast and slow winds are consistent with expectations from a thermal-magnetized disk wind model and are generally inconsistent with a purely thermal wind. Both the models and the observational analysis indicate that wind absorption occurs preferentially from the inner disk, offering a wind diagnostic in complement to optical forbidden line emission that traces the wind in larger volumes., Comment: 27 pages, 16 figures, plus appendix. Accepted for publication in ApJ

Variability in young stellar objects (YSOs) can be caused by various time-dependent phenomena associated with star formation, including accretion rates, geometric changes in the circumstellar disks, stochastic hydromagnetic interactions between stellar surfaces and inner disk edges, reconnections within the stellar magnetosphere, and hot/cold spots on stellar surfaces. We uncover and characterize $\sim$1700 variables from a sample of $\sim$5400 YSOs in nearby low-mass star-forming regions using mid-IR light curves obtained from the 6.5-years NEOWISE All Sky Survey. The mid-IR variability traces a wide range of dynamical, physical, and geometrical phenomenon. We classify six types of YSO mid-IR variability based on their light curves: secular variability ($Linear, Curved, Periodic$) and stochastic variability ($Burst, Drop, Irregular$). YSOs in earlier evolutionary stages have higher fractions of variables and higher amplitudes for the variability, with the recurrence timescale of FUor-type outbursts (defined here as $\Delta$W1 or $\Delta$W2 $>1$ mag followed by inspection of candidates) of $\sim$1000 years in the early embedded protostellar phase. Known eruptive young stars and subluminous objects show fractions of variables similar to the fraction ($\sim55\%$) found in typical protostars, suggesting that these two distinct types are not distinct in variability over the 6.5-year timescale. Along with brightness variability, we also find a diverse range of secular color variations, which can be attributed to a competitive interplay between the variable accretion luminosity of the central source and the variable extinction by material associated with the accretion process., Comment: Accepted for publication in ApJ