Search

While protoplanetary disks (PPDs) are generally thought to dissipate within several Myr, recent observations have revealed gas in debris disks. The origin of this gas remains uncertain, with one possibility being the unexpectedly long survival of PPDs (the primordial-origin scenario). To explore the plausibility of this scenario, we conduct 1D disk evolution simulations, varying parameters like stellar mass, disk mass, turbulent stress, and magnetohydrodynamic winds, while incorporating stellar evolution to account for time-varying photoevaporation rates. Our focus is on disks where small grains are depleted, as these are potentially long-lived due to reduced far-ultraviolet photoevaporation. Our results show that gas in these disks can survive beyond 10 Myr regardless of the stellar mass, provided they are initially massive ($M_{\mathrm{disk}}\approx 0.1M_*$) with relatively weak turbulent stress ($\alpha \ll 10^{-2}$). The longest lifetimes are consistently found for $M_* = 2 M_{\odot}$ across a wide parameter space, with gas typically persisting at $\sim 10$--$10^3 \mathrm{au}$. Roughly estimated CO masses for these disks fall within the observed range for the most massive gas-rich debris disks around early A stars. These alignments support the plausibility of the primordial-origin scenario. Additionally, our model predicts that accretion persists for as long as the disk survives, which could explain the accretion signatures detected in old disks hosted by low-mass stars, including Peter Pan disks. Our finding also suggests that ongoing accretion may exist in gas-rich debris disks. Thus, searching for accretion signatures could be a key factor to identifying the origin of gas in debris disks., Comment: 20 pages, 16 figures

Ultraviolet flare emission can drive photochemistry in exoplanet atmospheres and even serve as the primary source of uncertainty in atmospheric retrievals. Additionally, flare energy budgets are not well-understood due to a paucity of simultaneous observations. We present new near-UV (NUV) and optical observations of flares from three M dwarfs obtained at 20 s cadence with Swift and TESS, along with a re-analysis of flares from two M dwarfs in order to explore the energy budget and timing of flares at NUV--optical wavelengths. We find a 9000 K blackbody underestimates the NUV flux by $\geq$2$\times$ for 54$\pm$14% of flares and 14.8$\times$ for one flare. We report time lags between the bands of 0.5--6.6 min and develop a method to predict the qualitative flare shape and time lag to 36$\pm$30% accuracy. The scatter present in optical-NUV relations is reduced by a factor of 2.0$\pm$0.6 when comparing the total NUV energy with the TESS energy during the FWHM duration due to the exclusion of the $T_\mathrm{eff}\approx$5000 K tail. We show the NUV light curve can be used to remove flares from the optical light curve and consistently detect planets with 20% smaller transits than is possible without flare detrending. Finally, we demonstrate a 10$\times$ increase in the literature number of multi-wavelength flares with the Early eVolution Explorer (EVE), an astrophysics Small Explorer concept to observe young clusters with simultaneous NUV and optical bands in order to detect young planets, assess their photochemical radiation environments, and observe accretion., Comment: 27 pages, 13 figures, 4 tables, accepted to The Astronomical Journal

The composition of protoplanetary disks, and hence the initial conditions of planet formation, may be strongly influenced by the infall and thermal processing of material during the protostellar phase. Composition of dust and ice in protostellar envelopes, shaped by energetic processes driven by the protostar, serves as the fundamental building material for planets and complex organic molecules. As part of the JWST GO program, "Investigating Protostellar Accretion" (IPA), we observed an intermediate-mass protostar HOPS 370 (OMC2-FIR3) using NIRSpec/IFU and MIRI/MRS. This study presents the gas and ice phase chemical inventory revealed with the JWST in the spectral range of $\sim$2.9 to 28 $\mu$m and explores the spatial variation of volatile ice species in the protostellar envelope. We find evidence for thermal processing of ice species throughout the inner envelope. We present the first high-spatial resolution ($\sim 80$ au) maps of key volatile ice species H$_{2}$O, CO$_{2}$, $^{13}$CO$_2$, CO, and OCN$^-$, which reveal a highly structured and inhomogeneous density distribution of the protostellar envelope, with a deficiency of ice column density that coincides with the jet/outflow shocked knots. Further, we observe high relative crystallinity of H$_{2}$O ice around the shocked knot seen in the H$_2$ and OH wind/outflow, which can be explained by a lack of outer colder material in the envelope along the line of sight due to the irregular structure of the envelope. These observations show clear evidence of thermal processing of the ices in the inner envelope, close to the outflow cavity walls, heated by the luminous protostar., Comment: Submitted to ApJ. Main text:16 pages with 11 figures

Emission lines from complex organic molecules in B335 were observed in four epochs, spanning a luminosity burst of about 10 years duration. The emission lines increased dramatically in intensity as the luminosity increased, but they have decreased only slightly as the luminosity has decreased. This behavior agrees with expectations of rapid sublimation as the dust temperature increases, but slower freeze-out after the dust temperature drops. Further monitoring of this source, along with detailed chemical models, will exploit this natural laboratory for astrochemistry., Comment: Accepted for publication in ApJL

HOPS 87 is a Class 0 protostellar core known to harbor an extremely young bipolar outflow and a hot corino. We report the discovery of localized, chemically rich regions near the bases of the two-lobe bipolar molecular outflow in HOPS 87 containing molecules such as H$_2$CO, $^{13}$CS, H$_2$S, OCS, and CH$_3$OH, the simplest complex organic molecule (COM). The locations and kinematics suggest that these localized features are due to jet-driven shocks rather than being part of the hot corino region encasing the protostar. The COM compositions of the molecular gas in these jet-localized regions are relatively simpler than those in the hot corino zone. We speculate that this simplicity is due to either the liberation of ice with a less complex chemical history or the effects of shock chemistry. Our study highlights the dynamic interplay between the protostellar bipolar outflow, disk, inner core environment, and the surrounding medium, contributing to our understanding of molecular complexity in solar-like young stellar objects., Comment: 16 pages, 6+2 figures, accepted by ApJ

Molecular clouds (MCs) are the birthplaces of new stars in galaxies. A key component of MCs are photodissociation regions (PDRs), where far-ultraviolet radiation plays a crucial role in determining the gas's physical and chemical state. Traditional PDR models assume chemical steady state (CSS), where the rates of H$_2$ formation and photodissociation are balanced. However, real MCs are dynamic and can be out of CSS. In this study, we demonstrate that combining H$_2$ emission lines observed in the far-ultraviolet or infrared with column density observations can be used to derive the rates of H$_2$ formation and photodissociation. We derive analytical formulae that relate these rates to observable quantities, which we validate using synthetic H$_2$ line emission maps derived from the SILCC-Zoom hydrodynamical simulation. Our method estimates integrated H$_2$ formation and dissociation rates to within 29\% accuracy. Our simulations cover a wide dynamic range in H$_2$ formation and photodissociation rates, showing significant deviations from CSS, with 74\% of the MC's mass deviating from CSS by a factor greater than 2. Our analytical formulae can effectively distinguish between regions in and out of CSS. When applied to actual H$_2$ line observations, our method can assess the chemical state of MCs, providing insights into their evolutionary stages and lifetimes., Comment: Submitted to ApJ. **Comments Welcome!** This paper is the first ("Paper-I") in 2-paper series. Paper II, entitled "The Molecular Cloud Lifecycle II: Formation and Destruction of Molecular Clouds Diagnosed via H 2 Fluorescent Emission Emission" by Burkhart et al has been recently accepted for publication in ApJ, and may be found here: arXiv:2402.01587

We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) observations of warm CO and H$_2$O gas in emission toward the low-mass protostar IRAS 15398-3359, observed as part of the CORINOS program. The CO is detected via the rovibrational fundamental band and hot band near 5 $\mu$m, whereas the H$_2$O is detected in the rovibrational bending mode at 6-8 $\mu$m. Rotational analysis indicates that the CO originates in a hot reservoir of $1598\pm118$ K, while the water is much cooler at $204\pm 7$ K. Neither the CO nor the H$_2$O line images are significantly spatially extended, constraining the emission to within $\sim$40 au of the protostar. The compactness and high temperature of the CO are consistent with an origin in the embedded protostellar disk, or a compact disk wind. In contrast, the water must arise from a cooler region and requires a larger emitting area (compared to CO) to produce the observed fluxes. The water may arise from a more extended part of the disk, or from the inner portion of the outflow cavity. Thus, the origin of the molecular emission observed with JWST remains ambiguous. Better constraints on the overall extinction, comparison with realistic disk models, and future kinematically-resolved observations may all help to pinpoint the true emitting reservoirs., Comment: 15 pages, 13 figures. Accepted for publication in ApJ

Photoevaporation driven by hydrogen-ionizing radiation, also known as extreme-ultraviolet (EUV), profoundly shapes the lives of diverse astrophysical objects. Focusing here mainly on the dispersal of protoplanetary disks, we construct an analytical model accounting for the finite timescales of photoheating and photoionization. The model offers improved estimates for the ionization, temperature, and velocity structures versus distance from the central source, for a given EUV emission rate and spectral hardness. Compared to the classical picture of fully-ionized and isothermal winds with temperatures $\approx 10^4{\rm \,K}$ and speeds $\approx 10{\rm \,km\,s^{-1}}$, our model unveils broader hydrodynamical and thermochemical states of photoevaporative winds. In contrast to the classical picture, T~Tauri stars with EUV luminosities around $10^{30}{\rm \,erg\,s^{-1}}$ have non-isothermal ionized winds at lower temperatures than the classical value if the spectrum is soft, with an average deposited energy per photoionization less than about 3.7\,eV. Conversely, if the spectrum is hard, the winds tend to be atomic and isothermal at most radii in the disk. For lower EUV intensities, even with soft spectra, atomic winds can emerge beyond $\sim 10{\, \rm au}$ through advection. We demonstrate that the analytical model's predictions are in general agreement with detailed radiation-hydrodynamics calculations. The model furthermore illustrates how the energy efficiency of photoevaporation varies with the intensity and spectral hardness of the EUV illumination, as well as addressing discrepancies in the literature around the effectiveness of X-ray photoevaporation. These findings highlight the importance of considering the finite timescales of photoheating and photoionization, both in modeling and in interpreting observational data., Comment: Submitted to ApJ

The role of metallicity in shaping protoplanetary disk evolution remains poorly comprehended. This study analyzes the disk fraction of 10 young (0.9-2.1 Myr) and low-metallicity (0.34-0.83 Z$_{\odot}$) clusters located in the outer Milky Way with Galactocentric distances between 10 and 13 kpc. Using $JHK$ data obtained from UKIDSS, the calculated disk fraction values for low-mass stars (0.2-2 M$_{\odot}$) ranged from 42% to 7%. To enhance the statistical reliability of our analysis, eight additional low-metallicity clusters are sourced from previous studies with metallicity range 0.25-0.85 Z$_{\odot}$ along with our sample, resulting in a total of 18 regions with low-metallicity. We find that low-metallicity clusters exhibit on average $2.6\pm0.2$ times lower disk fraction compared to solar-metallicity clusters in all the age bins we have. Within the age range we can probe, our study does not find evidence of faster disk decay in sub-solar metallicity regions compared to solar-metallicity regions. Furthermore, we observe a positive correlation between cluster disk fraction and metallicity for two different age groups of 0.3-1.4 and 1.4-2.5 Myr. We emphasize that both cluster age and metallicity significantly affect the fraction of stars with evidence of inner disks., Comment: 24 pages, 17 figures, 2 tables, accepted for publication in The Astrophysical Journal

Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative strengths of the observed OH lines provide a powerful diagnostic of the ultraviolet radiation field in a heavily-extinguished region ($A_V \sim 10 - 20$) where direct UV observations are impossible. To high precision, the OH line strengths are consistent with a picture in which the suprathermal OH states are populated following the photodissociation of water in its $\tilde B - X$ band by ultraviolet radiation produced by fast ($\sim 80\,\rm km\,s^{-1}$) shocks along the jet. The observed dominance of emission from symmetric ($A^\prime$) OH states over that from antisymmetric ($A^{\prime\prime}$) states provides a distinctive signature of this particular population mechanism. Moreover, the variation of intensity with rotational quantum number suggests specifically that Ly$\alpha$ radiation is responsible for the photodissociation of water, an alternative model with photodissociation by a 10$^4$ K blackbody being disfavored at a high level of significance. Using measurements of the Br$\alpha$ flux to estimate the Ly$\alpha$ production rate, we find that $\sim 4\%$ of the Ly$\alpha$ photons are absorbed by water. Combined with direct measurements of water emissions in the $\nu_2 = 1 -0$ band, the OH observations promise to provide key constraints on future models for the diffusion of Ly$\alpha$ photons in the vicinity of a shock front., Comment: 30 pages, 7 figures. Accepted for publication in ApJ Letters

Leveraging the high resolution, high sensitivity, and wide frequency coverage of the Atacama Large Millimeter/submillimeter Array (ALMA), the QUARKS survey, standing for "Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures", is observing 139 massive star-forming clumps at ALMA Band 6 ($\lambda\sim$ 1.3 mm). This paper introduces the Atacama Compact Array (ACA) 7-m data. Combining multi-wavelength data, we provide the first edition of QUARKS atlas, offering insights into the multiscale and multiphase interstellar medium in high-mass star formation. The ACA 1.3 mm catalog includes 207 continuum sources that are called ACA sources. Their gas kinetic temperatures are estimated using three formaldehyde (H$_2$CO) transitions with a non-LTE radiation transfer model, and the mass and density are derived from a dust emission model. The ACA sources are massive (16-84 percentile values of 6-160 $M_{\odot}$), gravity-dominated ($M\propto R^{1.1}$) fragments within massive clumps, with supersonic turbulence ($\mathcal{M}>1$) and embedded star-forming protoclusters. We find a linear correlation between the masses of the fragments and the massive clumps, with a ratio of 6% between the two. When considering the fragments as representative of dense gas, the ratio indicates a dense gas fraction (DGF) of 6%, although with a wide scatter ranging from 1% to 10%. If we consider the QUARKS massive clumps to be what is observed at various scales, then the size-independent DGF indicates a self-similar fragmentation or collapsing mode in protocluster formation. With the ACA data over four orders of magnitude of luminosity-to-mass ratio ($L/M$), we find that the DGF increases significantly with $L/M$, which indicates clump evolutionary stage. We observed a limited fragmentation at the subclump scale, which can be explained by dynamic global collapse process., Comment: 24 pages, 7 figures. Accepted for publication in Research in Astronomy and Astrophysics. QUARKS atlas link: https://drive.google.com/file/d/1KTqXxCDduYepvLd9kIvZVSSytK48OmfL/view?usp=sharing

The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $\mu$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viable alternative tracer is the weaker $^{13}$CO$_2$ isotopologue band at 4.39 $\mu$m which has now become accessible at high S/N with the $\textit{James Webb}$ Space Telescope (JWST). We present JWST NIRSpec observations of $^{13}$CO$_2$ ice in five deeply embedded Class 0 sources spanning a wide range in luminosities (0.2 - 10$^4$ L$_{\odot}$ ) taken as part of the Investigating Protostellar Accretion Across the Mass Spectrum (IPA) program. The band profiles vary significantly, with the most luminous sources showing a distinct narrow peak at 4.38 $\mu$m. We first apply a phenomenological approach and show that a minimum of 3-4 Gaussian profiles are needed to fit the $^{13}$CO$_2$ absorption feature. We then combine these findings with laboratory data and show that a 15.2 $\mu$m $^{12}$CO$_2$ band inspired five-component decomposition can be applied for the isotopologue band where each component is representative of CO$_2$ ice in a specific molecular environment. The final solution consists of cold mixtures of CO$_2$ with CH$_3$OH, H$_2$O and CO as well as segregated heated pure CO$_2$ ice. Our results are in agreement with previous studies of the $^{12}$CO$_2$ ice band, further confirming that $^{13}$CO$_{2}$ is a useful alternative tracer of protostellar heating events. We also propose an alternative solution consisting only of heated CO$_2$:CH$_3$OH and CO$_2$:H$_2$O ices and warm pure CO$_2$ ice for decomposing the ice profiles of the two most luminous sources in our sample.

Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ formation, stellar radiation, sink particles, and turbulence. We find that, at early times in the cloud evolution, H$_2$ formation rapidly outpaces dissociation and molecular clouds build their mass from the atomic reservoir in their environment. Rapid H$_2$ formation is also associated with a higher early star formation rate. For the clouds studied here, H$_2$ is strongly out of chemical equilibrium during the early stages of cloud formation but settles into a bursty chemical steady-state about 2 Myrs after the first stars form. At the latest stage of cloud evolution, dissociation outweighs formation and the clouds enter a dispersal phase. We discuss how theories for the molecular cloud lifecycle and the star formation efficiency may be distinguished with observational measurements of H$_2$ fluorescence with a space-based high-resolution FUV spectrometer, such as the proposed Hyperion and Eos NASA Explorer missions. Such missions would enable measurements of the H$_2$ dissociation and formation rates, which we demonstrate can be connected to different phases in a molecular cloud's star-forming life, including cloud building, rapidly star-forming, H$_2$ chemical equilibrium, and cloud destruction., Comment: Submitted to ApJ, comments welcome

We investigate the bright CO fundamental emission in the central regions of five protostars in their primary mass assembly phase using new observations from JWST's Near-Infrared Spectrograph (NIRSpec) and Mid-Infrared Instrument (MIRI). CO line emission images and fluxes are extracted for a forest of $\sim$150 ro-vibrational transitions from two vibrational bands, $v=1-0$ and $v=2-1$. However, ${}^{13}$CO is undetected, indicating that ${}^{12}$CO emission is optically thin. We use H$_2$ emission lines to correct fluxes for extinction and then construct rotation diagrams for the CO lines with the highest spectral resolution and sensitivity to estimate rotational temperatures and numbers of CO molecules. Two distinct rotational temperature components are required for $v=1$ ($\sim600$ to 1000 K and 2000 to $\sim 10^4$ K), while one hotter component is required for $v=2$ ($\gtrsim 3500$ K). ${}^{13}$CO is depleted compared to the abundances found in the ISM, indicating selective UV photodissociation of ${}^{13}$CO; therefore, UV radiative pumping may explain the higher rotational temperatures in $v=2$. The average vibrational temperature is $\sim 1000$ K for our sources and is similar to the lowest rotational temperature components. Using the measured rotational and vibrational temperatures to infer a total number of CO molecules, we find that the total gas masses range from lower limits of $\sim10^{22}$ g for the lowest mass protostars to $\sim 10^{26}$ g for the highest mass protostars. Our gas mass lower limits are compatible with those in more evolved systems, which suggest the lowest rotational temperature component comes from the inner disk, scattered into our line of sight, but we also cannot exclude the contribution to the CO emission from disk winds for higher mass targets., Comment: 31 pages, 7 figures, 4 tables, received to ApJ December 10 2023, accepted to ApJ August 4 2024

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

While most debris disks consist of dust with little or no gas, a fraction has significant amounts of gas detected via emission lines of CO, ionized carbon, and/or atomic oxygen. Almost all such gaseous debris disks known are around A-type stars with ages up to 50 Myr. We show, using semi-analytic disk evolution modeling, that this can be understood if the gaseous debris disks are remnant protoplanetary disks that have become depleted of small grains compared to the interstellar medium. Photoelectric heating by the A-stars' FUV radiation is then inefficient, while the stars' EUV and X-ray emissions are weak owing to a lack of surface convective zones capable of driving magnetic activity. In this picture, stars outside the range of spectral types from A through early B are relatively hard to have such long-lived gas disks. Less-massive stars have stronger magnetic activity in the chromosphere, transition region, and corona with resulting EUV and X-ray emission, while more-massive stars have photospheres hot enough to produce strong EUV radiation. In both cases, primordial disk gas is likely to photoevaporate well before 50 Myr. These results come from 0D disk evolution models where we incorporate internal accretion stresses, MHD winds, and photoevaporation by EUV and X-ray photons with luminosities that are functions of the stellar mass and age. A key issue this work leaves open is how some disks become depleted in small dust so that FUV photoevaporation slows. Candidates include grains' growth, settling, radial drift, radiation force, and incorporation into planetary systems., Comment: Submitted to ApJL

Investigating Protostellar Accretion (IPA) is a JWST Cycle~1 GO program that uses NIRSpec IFU and MIRI MRS to obtain 2.9--28~$\mu$m spectral cubes of young, deeply embedded protostars with luminosities of 0.2 to 10,000~L$_{\odot}$ and central masses of 0.15 to 12~M$_{\odot}$. In this Letter, we report the discovery of a highly collimated atomic jet from the Class~0 protostar IRAS~16253$-$2429, the lowest luminosity source ($L_\mathrm{bol}$ = 0.2 $L_\odot$) in the IPA program. The collimated jet is detected in multiple [Fe~II] lines, [Ne~II], [Ni~II], and H~I lines, but not in molecular emission. The atomic jet has a velocity of about 169~$\pm$~15~km\,s$^{-1}$, after correcting for inclination. The width of the jet increases with distance from the central protostar from 23 to~60 au, corresponding to an opening angle of 2.6~$\pm$~0.5\arcdeg. By comparing the measured flux ratios of various fine structure lines to those predicted by simple shock models, we derive a shock {speed} of 54~km\,s$^{-1}$ and a preshock density of 2.0$\times10^{3}$~cm$^{-3}$ at the base of the jet. {From these quantities and using a suite of jet models and extinction laws we compute a mass loss rate between $0.4 -1.1\times10^{-10}~M_{\odot}$~yr~$^{-1}$.} The low mass loss rate is consistent with simultaneous measurements of low mass accretion rate ($2.4~\pm~0.8~\times~10^{-9}~M_{\odot}$~yr$^{-1}$) for IRAS~16253$-$2429 from JWST observations (Watson et al. in prep), indicating that the protostar is in a quiescent accretion phase. Our results demonstrate that very low-mass protostars can drive highly collimated, atomic jets, even during the quiescent phase., Comment: Accepted to ApJL. Comments and feedback welcome

Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9--28 $\mu$m spectral cubes of five young protostars with luminosities of 0.2-10,000 L$_{\odot}$ in their primary accretion phase. This paper introduces the NIRSpec 2.9--5.3 $\mu$m data of the inner 840-9000 au with spatial resolutions from 28-300 au. The spectra show rising continuum emission; deep ice absorption; emission from H$_{2}$, H~I, and [Fe~II]; and the CO fundamental series in emission and absorption. Maps of the continuum emission show scattered light cavities for all five protostars. In the cavities, collimated jets are detected in [Fe~II] for the four $< 320$~L$_{\odot}$ protostars, two of which are additionally traced in Br-$\alpha$. Knots of [Fe~II] emission are detected toward the most luminous protostar, and knots of [FeII] emission with dynamical times of $< 30$~yrs are found in the jets of the others. While only one jet is traced in H$_2$, knots of H$_2$ and CO are detected in the jets of four protostars. H$_2$ is seen extending through the cavities, showing that they are filled by warm molecular gas. Bright H$_2$ emission is seen along the walls of a single cavity, while in three cavities narrow shells of H$_2$ emission are found, one of which has an [Fe~II] knot at its apex. These data show cavities containing collimated jets traced in atomic/ionic gas surrounded by warm molecular gas in a wide-angle wind and/or gas accelerated by bow shocks in the jets., Comment: 26 pages, 12 figures

The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au at 350 GHz (870 um) in continuum and line emission. 248 dense cores were identified, including 106 cores showing protostellar signatures and 142 prestellar core candidates. Compared to early-stage infrared dark clouds (IRDCs) by ASHES, the core mass and surface density within the ASSEMBLE clumps exhibited significant increment, suggesting concurrent core accretion during the evolution of the clumps. The maximum mass of prestellar cores was found to be 2 times larger than that in IRDCs, indicating evolved protoclusters have the potential to harbor massive prestellar cores. The mass relation between clumps and their most massive core (MMCs) is observed in ASSEMBLE but not in IRDCs, which is suggested to be regulated by multiscale mass accretion. The mass correlation between the core clusters and their MMCs has a steeper slope compared to that observed in stellar clusters, which can be due to fragmentation of the MMC and stellar multiplicity. We observe a decrease in core separation and an increase in central concentration as protoclusters evolve. We confirm primordial mass segregation in the ASSEMBLE protoclusters, possibly resulting from gravitational concentration and/or gas accretion., Comment: 37 pages, 13 figures, 5 tables; accepted for publication in ApJS

Modeling self-gravitating gas flows is essential to answering many fundamental questions in astrophysics. This spans many topics including planet-forming disks, star-forming clouds, galaxy formation, and the development of large-scale structures in the Universe. However, the nonlinear interaction between gravity and fluid dynamics offers a formidable challenge to solving the resulting time-dependent partial differential equations (PDEs) in three dimensions (3D). By leveraging the universal approximation capabilities of a neural network within a mesh-free framework, physics informed neural networks (PINNs) offer a new way of addressing this challenge. We introduce the gravity-informed neural network (GRINN), a PINN-based code, to simulate 3D self-gravitating hydrodynamic systems. Here, we specifically study gravitational instability and wave propagation in an isothermal gas. Our results match a linear analytic solution to within 1\% in the linear regime and a conventional grid code solution to within 5\% as the disturbance grows into the nonlinear regime. We find that the computation time of the GRINN does not scale with the number of dimensions. This is in contrast to the scaling of the grid-based code for the hydrodynamic and self-gravity calculations as the number of dimensions is increased. Our results show that the GRINN computation time is longer than the grid code in one- and two- dimensional calculations but is an order of magnitude lesser than the grid code in 3D with similar accuracy. Physics-informed neural networks like GRINN thus show promise for advancing our ability to model 3D astrophysical flows.

We report the detection of the millimeter CH$_3$OH masers including a new detection of class I (11$_{0,11}$-10$_{1,10}$A) and class II (6$_{1,5}$-5$_{2,4}$E) maser transitions toward the high-mass protostar S255IR NIRS3 in post-burst phase. The CH$_3$OH emissions were detected as a mixture of maser and thermal characteristics. We examine the detected transitions using an excitation diagram and LTE model spectra and compare the observed properties with those of thermal lines. Class II CH$_3$OH maser transitions showed distinctive intensity and velocity distributions from those of thermal transitions. Bright distinct emission components in addition to the fragmented and arc-shaped emissions are only detected in class I CH$_3$OH maser transitions toward southern and western directions from the protostellar position, implying the presence of the slow outflow shocks., Comment: 9 pages, 4 Figures, 1 Table, Accepted for publication in ApJL

Hot corinos are of great interest due to their richness in interstellar complex organic molecules (COMs) and the consequent potential prebiotic connection to solar-like planetary systems. Recent surveys have reported an increasing number of hot corino detections in Class 0/I protostars; however, the relationships between their physical properties and the hot-corino signatures remain elusive. In this study, our objective is to establish a general picture of the detectability of the hot corinos by identifying the origin of the hot-corino signatures in the sample of young stellar objects (YSOs) obtained from the Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project. We apply spectral energy distribution (SED) modeling to our sample and identify the physical parameters of the modeled YSOs directly, linking the detection of hot-corino signatures to the envelope properties of the YSOs. Imaging simulations of the methanol emission further support this scenario. We, therefore, posit that the observed COM emission originates from the warm inner envelopes of the sample YSOs, based on both the warm region size and the envelope density profile. The former is governed by the source luminosity and is additionally affected by the disk and cavity properties, while the latter is related to the evolutionary stages. This scenario provides a framework for detecting hot-corino signatures toward luminous Class 0 YSOs, with fewer detections observed toward similarly luminous Class I sources., Comment: 28 pages, 11 figures

Most stars form in multiple star systems. For a better understanding of their formation processes, it is important to resolve the individual protostellar components and the surrounding envelope and disk material at the earliest possible formation epoch because the formation history can be lost in a few orbital timescales. Here we present the ALMA observational results of a young multiple protostellar system, IRAS 04239+2436, where three well-developed large spiral arms were detected in the shocked SO emission. Along the most conspicuous arm, the accretion streamer was also detected in the SO$_2$ emission. The observational results are complemented by numerical magneto-hydrodynamic simulations, where those large arms only appear in magnetically weakened clouds. The numerical simulations also suggest that the large triple spiral arms are the result of gravitational interactions between compact triple protostars and the turbulent infalling envelope., Comment: Accepted for publication in ApJ

The solar system's Zodiacal Cloud is visible to the unaided eye, yet the origin of its constituent dust particles is not well understood, with a wide range of proposed divisions between sources in the asteroid belt and Jupiter Family comets. The amount of dust contributed by Oort Cloud comets is uncertain. Knowledge of the Zodiacal Cloud's structure and origins would help with NASA's aim of characterizing potentially Earth-like planets around nearby stars, since the exo-Earths must be studied against the light scattered from extrasolar analogs of our cloud. As the only example where the parent bodies can be tracked, our own cloud is critical for learning how planetary system architecture governs the interplanetary dust's distribution. Our cloud has been relatively little-studied in the near-ultraviolet, a wavelength range that is important for identifying potentially-habitable planets since it contains the broad Hartley absorption band of ozone. We show through radiative transfer modeling that our cloud's shape and size at near-UV wavelengths can be measured from Earth orbit by mapping the zodiacal light's flux and linear polarization across the sky. We quantify how well the cloud's geometric and optical properties can be retrieved from a set of simulated disk observations, using a Markov chain Monte Carlo analysis. The results demonstrate that observations with sufficient precision, covering a set of fields distributed along the ecliptic and up to the poles, can be used to determine the division between asteroidal, Jupiter Family, and Oort Cloud dust components, primarily via their differing orbital inclination distributions. We find that the observations must be repeated over a time span of several months in order to disentangle the zodiacal light from the Galactic background using the Milky Way's rotation across the sky., Comment: submitted to PASP

Starless cores represent the initial stage of evolution toward (proto)star formation, and a subset of them, known as prestellar cores, with high density (~ 10^6 cm^-3 or higher) and being centrally concentrated are expected to be embryos of (proto)stars. Determining the density profile of prestellar cores, therefore provides an important opportunity to gauge the initial conditions of star formation. In this work, we perform rigorous modeling to estimate the density profiles of three nearly spherical prestellar cores among a sample of five highly dense cores detected by our recent observations. We employed multi-scale observational data of the (sub)millimeter dust continuum emission including those obtained by SCUBA-2 on the JCMT with a resolution of ~5600 au and by multiple ALMA observations with a resolution as high as ~480 au. We are able to consistently reproduce the observed multi-scale dust continuum images of the cores with a simple prescribed density profile, which bears an inner region of flat density and a r^-2 profile toward the outer region. By utilizing the peak density and the size of the inner flat region as a proxy for the dynamical stage of the cores, we find that the three modeled cores are most likely unstable and prone to collapse. The sizes of the inner flat regions, as compact as ~500 au, signify them being the highly evolved prestellar cores rarely found to date., Comment: Accepted for publication in ApJ

We report the first detection of radio recombination lines (RRLs) of ions heavier than helium. In a highly sensitive multi-band (12--50 GHz) line survey toward Orion KL with the TianMa 65-m Radio Telescope (TMRT), we successfully detected more than fifteen unblended $\alpha$ lines of RRLs of singly ionized species (XII) recombined from XIII. The Ka-band (35--50 GHz) spectrum also shows tentative signals of $\beta$ lines of ions. The detected lines can be successfully crossmatched with the the rest frequencies of RRLs of CII and/or OII. This finding greatly expands the connotation of ion RRLs, since before this work only two blended lines (105$\alpha$ and 121$\alpha$) of HeII had been reported. Our detected lines can be fitted simultaneously under assumption of local thermodynamic equilibrium (LTE). An abundance of CIII and OIII of 8.8$\times$10$^{-4}$ is obtained, avoiding the complexities of optical/infrared observations and the blending of RRLs of atoms. It is consistent with but approaches the upper bound of the value (10$^{-4}$--$10^{-3}$) estimated from optical/infrared observations. The effects of dielectronic recombination may contribute to enhancing the level populations even at large $n$. We expect future observations using radio interferometers could break the degeneracy between C and O, and help to reveal the ionization structure and dynamical evolution of various ionized regions., Comment: 3 figures, 2 tables, accepted by A&A Letter

We present Hyperion, a mission concept recently proposed to the December 2021 NASA Medium Explorer announcement of opportunity. Hyperion explores the formation and destruction of molecular clouds and planet-forming disks in nearby star-forming regions of the Milky Way. It does this using long-slit, high-resolution spectroscopy of emission from fluorescing molecular hydrogen, which is a powerful far-ultraviolet (FUV) diagnostic. Molecular hydrogen (H2) is the most abundant molecule in the universe and a key ingredient for star and planet formation, but is typically not observed directly because its symmetric atomic structure and lack of a dipole moment mean there are no spectral lines at visible wavelengths and few in the infrared. Hyperion uses molecular hydrogen's wealth of FUV emission lines to achieve three science objectives: (1) determining how star formation is related to molecular hydrogen formation and destruction at the boundaries of molecular clouds; (2) determining how quickly and by what process massive star feedback disperses molecular clouds; and (3) determining the mechanism driving the evolution of planet-forming disks around young solar-analog stars. Hyperion conducts this science using a straightforward, highly-efficient, single-channel instrument design. Hyperion's instrument consists of a 48 cm primary mirror, with an f/5 focal ratio. The spectrometer has two modes, both covering 138.5-161.5 nm bandpasses. A low resolution mode has a spectral resolution of R>10,000 with a slit length of 65 arcmin, while the high resolution mode has a spectral resolution of R>50,000 over a slit length of 5 armin. Hyperion occupies a 2 week long, high-earth, Lunar resonance TESS-like orbit, and conducts 2 weeks of planned observations per orbit, with time for downlinks and calibrations. Hyperion was reviewed as Category I, which is the highest rating possible, but was not selected., Comment: Accepted to JATIS, 9 Figures

Models of the protostellar source, B335, are developed using axisymmetric three-dimensional models to resolve conflicts found in one-dimensional models. The models are constrained by a large number of observations, including ALMA, Herschel, and Spitzer data. Observations of the protostellar source B335 with ALMA show red-shifted absorption against a central continuum source indicative of infall in the HCO$^+$ and HCN $J = 4\rightarrow 3$ transitions. The data are combined with a new estimate of the distance to provide strong constraints to three-dimensional radiative transfer models including a rotating, infalling envelope, outflow cavities, and a very small disk. The models favor ages since the initiation of collapse between $3 \times 10^4$ and $4 \times 10^4$ yr for both the continuum and the lines, resolving a conflict found in one-dimensional models. The models under-predict the continuum emission seen by ALMA, suggesting an additional component such as a pseudo-disk. The best-fitting model is used to convert variations in the 4.5 $\mu m$ flux in recent years into a model for a variation of a factor of 5-7 in luminosity over the last 8 years., Comment: 22 pages, 20 figures, in press ApJ

We model substructure in the protoplanetary disks around DoAr 44 and HD 163296 in order to better understand the conditions under which planets may form. We match archival millimeter-wavelength thermal emission against models of the disks' structure that are in radiation balance with the starlight heating and in vertical hydrostatic equilibrium, and then compare to archival polarized scattered near-infrared images of the disks. The millimeter emission arises in the interior, while the scattered near-infrared radiation probes the disks' outer layers. Our best model of the HD 163296 disk has dust masses $81\pm 13$ $M_\oplus$ in the inner ring at 68 au and $82^{+26}_{-16}$ $M_\oplus$ in the outer ring at 102 au, both falling within the range of estimates from previous studies. Our DoAr 44 model has total dust mass $84^{+7.0}_{-3.5}$ $M_\oplus$. Unlike HD 163296, DoAr 44 as of yet has no detected planets. If the central cavity in the DoAr 44 disk is caused by a planet, the planet's mass must be at least 0.5 $M_J$ and is unlikely to be greater than 1.6 $M_J$. We demonstrate that the DoAr 44 disk's structure with a bright ring offset within a fainter skirt can be formed by dust particles drifting through a plausible distribution of gas., Comment: 18 pages, 19 figures

The origin of complex organic molecules (COMs) in young Class 0 protostars has been one of the major questions in astrochemistry and star formation. While COMs are thought to form on icy dust grains via gas-grain chemistry, observational constraints on their formation pathways have been limited to gas-phase detection. Sensitive mid-infrared spectroscopy with JWST enables unprecedented investigation of COM formation by measuring their ice absorption features. We present an overview of JWST/MIRI MRS spectroscopy and imaging of a young Class 0 protostar, IRAS 15398-3359, and identify several major solid-state absorption features in the 4.9-28 $\mu$m wavelength range. These can be attributed to common ice species, such as H$_2$O, CH$_3$OH, NH$_3$, and CH$_4$, and may have contributions from more complex organic species, such as C$_2$H$_5$OH and CH$_3$CHO. The MRS spectra show many weaker emission lines at 6-8 $\mu$m, which are due to warm CO gas and water vapor, possibly from a young embedded disk previously unseen. Finally, we detect emission lines from [Fe II], [Ne II], [S I], and H$_2$, tracing a bipolar jet and outflow cavities. MIRI imaging serendipitously covers the south-western (blue-shifted) outflow lobe of IRAS 15398-3359, showing four shell-like structures similar to the outflows traced by molecular emission at sub-mm wavelengths. This overview analysis highlights the vast potential of JWST/MIRI observations and previews scientific discoveries in the coming years., Comment: Accepted to ApJL. The methods of data reduction and spectral extraction are updated

We have mapped HCN and HCO$^{+}$ (J = 1 $\to$ 0) line emission toward a sample of seven star-forming regions (with 12 + log[O/H] range from 8.34 to 8.69) in the outer Milky Way (Galactocentric distance > 9.5 kpc), using the 14-meter radio telescope of the Taeduk Radio Astronomy Observatory (TRAO). We compare these two molecular lines with other conventional tracers of dense gas, millimeter-wave continuum emission from dust and extinction thresholds ($A_{V} \geq 8$ mag), inferred from the $^{13}$CO line data. HCN and HCO$^{+}$ correlate better with the millimeter emission than with the extinction criterion. A significant amount of luminosity comes from regions below the extinction criterion and outside the millimeter clump for all the clouds. The average fraction of HCN luminosity from within the regions with $A_{V} \geq 8$ mag is $0.343\pm0.225$; for the regions of millimeter emission, it is $0.478\pm0.149$. Based on a comparison with column density maps from Herschel, HCN and HCO$^{+}$ trace dense gas in high column density regions better than does $^{13}$CO. HCO$^{+}$ is less concentrated than HCN for outer Galaxy targets, in contrast with the inner Galaxy sample, suggesting that metallicity may affect the interpretation of tracers of dense gas. The conversion factor between the dense gas mass ($M_{dense}$) and line luminosities of HCN and HCO$^{+}$, when integrated over the whole cloud, is comparable with factors used in extragalactic studies., Comment: 20 pages, 8 figures, 5 tables, accepted for publication in The Astronomical Journal

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation., Comment: 16 pages

Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet., Comment: 16 pages, 7 figures. accepted for publication in ApJ

The observed star formation rate of the Milky Way can be explained by applying a metallicity-dependent factor to convert CO luminosity to molecular gas mass and a star formation efficiency per free-fall time that depends on the virial parameter of a molecular cloud. These procedures also predict the trend of star formation rate surface density with Galactocentric radius. The efficiency per free-fall time variation with virial parameter plays the major role in bringing theory into agreement with observations for the total star formation rate, while the metallicity dependence of the CO luminosity to mass conversion is most notable in the variation with Galactocentric radius. Application of these changes resolves a factor of over 100 discrepancy between observed and theoretical star formation rates that has been known for nearly 50 years., Comment: In press, Astrophysical Journal Letters 11 pages, 6 figures

In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $\delta V$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak velocities to the H$^{13}$CO$^+$ line width, was biased toward negative values, suggesting that inward motions were more dominant than outward motions. Three starless cores (10% of all starless cores surveyed) were identified as cores with blue-skewed line profiles (asymmetric profiles with more intense blue-shifted emission), and another two starless cores (7%) were identified as candidate blue-skewed line profiles. The peak velocity difference between HCO$^+$ and H$^{13}$CO$^+$ of them was up to 0.9 km s$^{-1}$, suggesting that some inward motions exceeded the speed of sound for the quiescent gas ($\sim10-17$ K). The mean of $\delta V$ of the five aforementioned starless cores was derived to be $-$0.5$\pm$0.3. One core, G211.16$-$19.33North3, observed using the ALMA ACA in DCO$^+$ $J = 3\rightarrow2$ exhibited blue-skewed features. Velocity offset in the blue-skewed line profile with a dip in the DCO$^+$ $J = 3\rightarrow2$ line was larger ($\sim 0.5$ km s$^{-1}$) than that in HCO$^+$ $J = 1\rightarrow0$ ($\sim 0.2$ km s$^{-1}$), which may represent gravitational acceleration of inward motions. It seems that this core is at the last stage in the starless phase, judging from the chemical evolution factor version 2.0 (CEF2.0)., Comment: 45 pages, 28 figures, ApJ in press

We present a comprehensive study of the gas kinematics associated with density structures at different spatial scales in the filamentary infrared dark cloud, G034.43+00.24 (G34). This study makes use of the H13CO+ (1-0) molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey, which has spatial and velocity resolution of 0.04 pc and 0.2 km/s, respectively. Several tens of dendrogram structures have been extracted in the position-position-velocity space of H13CO+, which include 21 small-scale leaves and 20 larger-scale branches. Overall, their gas motions are supersonic but they exhibit the interesting behavior where leaves tend to be less dynamically supersonic than the branches. For the larger-scale, branch structures, the observed velocity-size relation (i.e., velocity variation/dispersion versus size) are seen to follow the Larson scaling exponent while the smaller-scale, leaf structures show a systematic deviation and display a steeper slope. We argue that the origin of the observed kinematics of the branch structures is likely to be a combination of turbulence and gravity-driven ordered gas flows. In comparison, gravity-driven chaotic gas motion is likely at the level of small-scale leaf structures. The results presented in our previous paper and this current follow-up study suggest that the main driving mechanism for mass accretion/inflow observed in G34 varies at different spatial scales. We therefore conclude that a scale-dependent combined effect of turbulence and gravity is essential to explain the star-formation processes in G34., Comment: 11 pages, 6 figures, and 1 table. To appear in MNRAS

Understanding the thermal structure of protoplanetary disks is crucial for modeling planet formation and interpreting disk observations. We present a new two-layer radiative transfer model for computing the thermal structure of axisymmetric irradiated disks. Unlike the standard two-layer model, our model accounts for the radial as well as vertical transfer of the starlight reprocessed at the disk surface. The model thus allows us to compute the temperature below "shadowed" surfaces receiving no direct starlight. Thanks to the assumed axisymmetry, the reprocessed starlight flux is given in one-dimensional integral form that can be computed at a low cost. Furthermore, our model evolves the midplane temperature using a time-dependent energy equation and can therefore treat thermal instabilities. We apply our global two-layer model to disks with a planetary induced gap and confirm that the model reproduces the disks' temperature profiles obtained from more computationally expensive Monte Carlo radiative transfer calculations to an accuracy of less than 20%. We also apply the model to study the long-term behavior of the thermal wave instability in irradiated disks. Being simple and computationally efficient, the global two-layer model will be suitable for studying the interplay between disks' thermal evolution and dust evolution., Comment: 22 pages, 18 figures, accepted for publication in PASJ (Publications of the Astronomical Society of Japan)

Understanding the origin of accretion and dispersal of protoplanetary disks is fundamental for investigating planet formation. Recent numerical simulations show that launching winds are unavoidable when disks undergo magnetically driven accretion and/or are exposed to external UV radiation. Observations also hint that disk winds are common. We explore how the resulting wind mass loss rate can be used as a probe of both disk accretion and dispersal. As a proof-of-concept study, we focus on magnetocentrifugal winds, MRI (magnetorotational instability) turbulence, and external photoevapotaion. By developing a simple, yet physically motivated disk model and coupling it with simulation results available in the literature, we compute the mass loss rate as a function of external UV flux for each mechanism. We find that different mechanisms lead to different levels of mass loss rate, indicating that the origin of disk accretion and dispersal can be determined, by observing the wind mass loss rate resulting from each mechanism. This determination provides important implications for planet formation. This work thus shows that the ongoing and future observations of the wind mass loss rate for protoplanetary disks are paramount to reliably constrain how protoplanetary disks evolve with time and how planet formation takes place in the disks., Comment: 16 pages; 7 figures, 5 tables, accepted for publication in ApJL; no change in the conclusions; the additional Monte-Carlo based, parameter study was conducted, following the referee's comment

We present velocity resolved SOFIA/upGREAT observations of [OI] and [CII] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [OI] emission toward only the protostar and [CII] emission toward the protostar and the red-shifted outflow. The [OI] emission has a width of $\sim$100 km s$^{-1}$ only in the blue-shifted velocity, suggesting an origin in shocked gas. The [CII] lines are narrow, consistent with an origin in a photodissociation region. Differential dust extinction from the envelope due to the inclination of the outflows is the most likely cause of the missing red-shifted [OI] emission. Fitting the [OI] line profile with two Gaussian components, we find one component at the source velocity with a width of $\sim$20 km s$^{-1}$ and another extremely broad component at -30 km s$^{-1}$ with a width of 87.5 km s$^{-1}$, the latter of which has not been seen in L1551 IRS 5. The kinematics of these two components resemble cavity shocks in molecular outflows and spot shocks in jets. Radiative transfer calculations of the [OI], high-J CO, and H$_2$O lines in the cavity shocks indicate that [OI] dominates the oxygen budget, making up more than 70% of the total gaseous oxygen abundance and suggesting [O]/[H] of $\sim$1.5$\times$10$^{-4}$. Attributing the extremely broad [OI] component to atomic winds, we estimate the intrinsic mass loss rate of (1.3$\pm$0.8)$\times$ 10$^{-6}$ M$_{\odot}$ yr$^{-1}$. The intrinsic mass loss rates derived from low-J CO, [OI], and HI are similar, supporting the model of momentum-conserving outflows, where the atomic wind carries most momentum and drives the molecular outflows., Comment: 19 pages, 14 figures, accepted for publication in The Astrophysical Journal