In recent years evidence has been building that planet formation starts early, in the first $\sim$ 0.5 Myr. Studying the dust masses available in young disks enables understanding the origin of planetary systems since mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka-band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities. Using the DIANA opacity model including large grains, with a dust opacity value of $\kappa_{\rm 9\ mm}$ = 0.28 cm$^{2}$ g$^{-1}$, the median dust masses of the embedded disks in Perseus are 158 M$_\oplus$ for Class 0 and 52 M$_\oplus$ for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47 M$_\oplus$ and 12 M$_\oplus$ for Class 0 and Class I, respectively, obtained using the maximum dust opacity value $\kappa_{\rm 1.3mm}$ = 2.3 cm$^{2}$ g$^{-1}$. The dust masses of young Class 0 and I disks are larger by at least a factor of 10 and 3, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions. The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of $\sim$ 15%. Higher efficiency of $\sim$ 30% is necessary if the planet formation is set to start in Class I disks.

We present ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature.

Episodic accretion has been used to explain the wide range of protostellar luminosities, but its origin and influence on the star forming process are not yet fully understood. We present an ALMA survey of N$_2$H$^+$ ($1-0$) and HCO$^+$ ($3-2$) toward 39 Class 0 and Class I sources in the Perseus molecular cloud. N$_2$H$^+$ and HCO$^+$ are destroyed via gas-phase reactions with CO and H$_2$O, respectively, thus tracing the CO and H$_2$O snowline locations. A snowline location at a much larger radius than that expected from the current luminosity suggests that an accretion burst has occurred in the past which has shifted the snowline outward. We identified 18/18 Class 0 and 9/10 Class I post-burst sources from N$_2$H$^+$, and 7/17 Class 0 and 1/8 Class I post-burst sources from HCO$^+$.The accretion luminosities during the past bursts are found to be $\sim10-100~L_\odot$. This result can be interpreted as either evolution of burst frequency or disk evolution. In the former case, assuming that refreeze-out timescales are 1000 yr for \ce{H2O} and 10,000 yr for CO, we found that the intervals between bursts increases from 2400 yr in the Class 0 to 8000 yr in the Class I stage. This decrease in the burst frequency may reflect that fragmentation is more likely to occur at an earlier evolutionary stage when the young stellar object is more prone to instability.

We present ALMA long-baseline observations toward the Class 0 protostar IRAS 16253-2429 (hereafter IRAS 16253) with a resolution down to 0.12" (~15 au). The 1.3 mm dust continuum emission has a deconvolved Gaussian size of 0.16" x 0. 07" (20 au x 8.8 au), likely tracing an inclined dusty disk. Interestingly, the position of the 1.38 mm emission is offset from that of the 0.87 mm emission along the disk minor axis. Such an offset may come from a torus-like disk with very different optical depths between these two wavelengths. Furthermore, through CO (2 - 1) and C18O (2 - 1) observations, we study rotation and infall motions in this disk-envelope system and infer the presence of a Keplerian disk with a radius of 8 - 32 au. This result suggests that the disk could have formed by directly evolving from a first core, because IRAS16253 is too young to gradually grow a disk to such a size considering the low rotation rate of its envelope. In addition, we find a quadruple pattern in the CO emission at low velocity, which may originate from CO freeze out at the disk/envelope midplane. This suggests that the "cold disk" may appear in the early stage, implying a chemical evolution for the disk around this proto-brown dwarf (or very low-mass protostar) different from that of low-mass stars.

Episodic accretion has been proposed as a solution to the long-standing luminosity problem in star formation; however, the process remains poorly understood. We present observations of line emission from N2H+ and CO isotopologues using the Atacama Large Millimeter/submillimeter Array (ALMA) in the envelopes of eight Very Low Luminosity Objects (VeLLOs). In five of the sources the spatial distribution of emission from N2H+ and CO isotopologues shows a clear anti-correlation. It is proposed that this is tracing the CO snow line in the envelopes: N2H+ emission is depleted toward the center of these sources in contrast to the CO isotopologue emission which exhibits a peak. The positions of the CO snow lines traced by the N2H+ emission are located at much larger radii than those calculated using the current luminosities of the central sources. This implies that these five sources have experienced a recent accretion burst because the CO snow line would have been pushed outwards during the burst due to the increased luminosity of the central star. The N2H+ and CO isotopologue emission from DCE161, one of the other three sources, is most likely tracing a transition disk at a later evolutionary stage. Excluding DCE161, five out of seven sources (i.e., ~70%) show signatures of a recent accretion burst. This fraction is larger than that of the Class 0/I sources studied by J{\o}rgensen et al. (2015) and Frimann et al. (2016) suggesting that the interval between accretion episodes in VeLLOs is shorter than that in Class 0/I sources.

We present an outflow survey toward 20 Low Luminosity Objects (LLOs), namely protostars with an internal luminosity lower than 0.2 Lsun. Although a number of studies have reported the properties of individual LLOs, the reasons for their low luminosity remain uncertain. To answer this question, we need to know the evolutionary status of LLOs. Protostellar outflows are found to widen as their parent cores evolve, and therefore, the outflow opening angle could be used as an evolutionary indicator. The infrared scattered light escapes out through the outflow cavity and highlights the cavity wall, giving us the opportunity to measure the outflow opening angle. Using the Canada-France-Hawaii Telescope, we detected outflows toward eight LLOs out of 20 at Ks band, and based on archival Spitzer IRAC1 images, we added four outflow-driving sources from the remaining 12 sources. By fitting these images with radiative transfer models, we derive the outflow opening angles and inclination angles. To study the widening of outflow cavities, we compare our sample with the young stellar objects from Arce & Sargent 2006 and Velusamy et al. 2014 in the plot of opening angle versus bolometric temperature taken as an evolutionary indicator.Our LLO targets match well the trend of increasing opening angle with bolometric temperature reported by Arce & Sargent and are broadly consistent with that reported by Velusamy et al., suggesting that the opening angle could be a good evolutionary indicator for LLOs. Accordingly, we conclude that at least 40% of the outflow-driving LLOs in our sample are young Class 0 objects.

The formation mechanism of brown dwarfs (BDs) is one of the long-standing problems in star formation because the typical Jeans mass in molecular clouds is too large to form these substellar objects. To answer this question, it is crucial to study a BD at the embedded phase. IRAS 16253-2429 is classified as a very low luminosity object (VeLLO) with internal luminosity 0.1 Lsun. VeLLOs are believed to be very low-mass protostars or even proto-BDs. We observed the jet/outflow driven by IRAS 16253-2429 in CO (2-1), (6-5), and (7-6) using the IRAM 30 m and APEX telescopes and the SMA in order to study its dynamical features and physical properties. Our SMA map reveals two protostellar jets, indicating the existence of a proto-binary system as implied by the precessing jet detected in H2 emission. We detect a wiggling pattern in the position-velocity diagrams along the jet axes, which is likely due to the binary orbital motion. Based on this, we derive the current mass of the binary as ~0.032 Msun. Given the low envelope mass, IRAS 16253-2429 will form a binary that probably consist of one or two BDs. Furthermore, we found that the outflow force as well as the mass accretion rate are very low based on the multi-transition CO observations, which suggests that the final masses of the binary components are at the stellar/substellar boundary. Since IRAS 16253 is located in an isolated environment, we suggest that BDs can form through fragmentation and collapse like low-mass stars.

We present a survey toward 16 Low Luminosity Objects (LLOs with an internal luminosity, Lint, lower than 0.2 Lsun) with N2H+ (1-0), N2H+ (3-2), N2D+ (3-2), HCO+ (3-2) and HCN (3-2) using the Arizona Radio Observatory Kitt Peak 12m Telescope and Submillimeter Telescope. Our goal is to probe the nature of these faint protostars which are believed to be either very low mass or extremely young protostars. We find that the N2D+/N2H+ column density ratios of LLOs are similar to those of typical starless cores and Class 0 objects. The N2D+/N2H+ column density ratios are relatively high (> 0.05) for LLOs with kinetic temperatures less than 10 K in our sample. The distribution of N2H+ (1-0) line widths spreads between that of starless cores and young Class 0 objects. If we use the line width as a dynamic evolutionary indicator, LLOs are likely young Class 0 protostellar sources. We further use the optically thick tracers, HCO+ (3-2) and HCN (3-2), to probe the infall signatures of our targets. We derive the asymmetry parameters from both lines and estimate the infall velocities by fitting the HCO+ (3-2) spectra with two-layer models. As a result, we identify eight infall candidates based on the infall velocities and seven candidates have infall signatures supported by asymmetry parameters from at least one of HCO+ (3-2) and HCN (3-2).

We develop a new method to identify YSOs from star-forming regions using the photometry data from Spitzer's c2d Legacy Project. The aim is to obtain YSO lists as complete as possible for studying the statistical properties, such as Star Formation Rate (SFR) and lifetimes of YSOs in different evolutionary stages. The largest obstacle for identifying YSOs comes from background galaxies with similar SEDs to YSOs. Traditionally, selected color-color and color-magnitude criteria are used to separate YSOs and galaxies. However, since there is no obvious boundary between YSOs and galaxies in Color-Color Diagrams (CCDs) and Color-Magnitude Diagrams (CMDs), those criteria may exclude faint YSOs near the boundary. In this paper, we separate the YSOs and galaxies in multi-dimensional (Multi-D) magnitude space, which is equivalent to using all variations of CMDs simultaneously. Comparing sources from molecular clouds to Spitzer's SWIRE data, which have negligible amount of YSOs, we can naturally identify YSO candidates locating outside of the galaxy populated regions in the Multi-D space. In the five c2d-surveyed clouds, we select 322 new YSO candidates (YSOc), miss/exclude 33 YSOc compared to Evans et al. 2009 and result in 1313 YSOc in total. As a result, SFR increases 28% correspondingly, but the lifetimes of YSOs in different evolutionary stages remain unchanged. Comparing to theories Krumholz & McKee 2005, our derived SFR suggests that star formation in large scale is dominated by supersonic turbulence rather than magnetic fields. Furthermore, we identify 7 new Very Low Luminosity Objects.