Accurate measurements of statistical properties, such as the star formation rate and the lifetime of young stellar objects (YSOs) in different stages, is essential for constraining star formation theories. However, it is a difficult task to separate galaxies and YSOs based on spectral energy distributions (SEDs) alone, because they contain both thermal emission from stars and dust around them and no reliable theories can be applied to distinguish them. Here we develop a machine learning algorithm based on Convolutional Neural Network, named Spectrum Classifier of Astronomical Objects (SCAO), to classify regular stars, galaxies, and YSOs, solely based on their SEDs. Superior to previous classifiers, SCAO is solely trained by labeled data without a priori theoretical knowledge, and provides excellent results with high precision (>95%) and recall (>98%) for YSOs when data from only eight bands are included. We investigate the effects of observation errors and distance effects, and show that high accuracy performance is still maintained even when using fluxes of only three bands (IRAC 3, IRAC 4, and MIPS 1) in the long wavelengths regime. We apply SCAO to Spitzer Enhanced Imaging Products (SEIP), the most complete catalog of Spitzer observations, and found 136689 YSO candidates. Finally, based on results predicted by SCAO, we provide an intuitive contour plot for a direct identification of YSOs even without any calculation. The website from SCAO is available at http://scao.astr.nthu.edu.tw.

Mapping magnetic fields is the key to resolving what remains an unclear physical picture of circumstellar magnetic fields in late-type evolved stars. Observations of linearly polarized emission from thermal molecular line transitions due to the Goldreich-Kylafis (G-K) effect provides valuable insight into the magnetic field geometry in these sources that is complementary to other key studies. In this paper, we present the detection of spectral-line polarization from both the thermal $J=2-1$ CO line and the $v=1, J=5-4$ SiO maser line toward two thermal-pulsating (TP-) AGB stars, R Crt and R Leo. The observed fractional linear polarization in the CO emission is measured as $m_l\sim 3.1\%$ and $m_l\sim9.7\%$ for R Crt and R Leo respectively. A circumstellar envelope (CSE) model profile and the associated parameters are estimated and used as input to a more detailed modeling of the predicted linear polarization expected from the G-K effect. The observed thermal line polarization level is consistent with the predicted results from the G-K model for R Crt; additional effects need to be considered for R Leo.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

Polarized continuum emission from aligned grains in disks around young stellar objects can be used to probe the magnetic field, radiation anisotropy, or drift between dust and gas, depending on whether the non-spherical grains are aligned magnetically, radiatively or mechanically. We show that it can also be used to probe another key disk property -- the temperature gradient -- along sight lines that are optically thick, independent of the grain alignment mechanism. We first illustrate the technique analytically using a simple 1D slab model, which yields an approximate formula that relates the polarization fraction to the temperature gradient with respect to the optical depth tau at the tau=1 surface. The formula is then validated using models of stellar irradiated disks with and without accretion heating. The promises and challenges of the technique are illustrated with a number of Class 0 and I disks with ALMA dust polarization data, including NGC 1333 IRAS4A1, IRAS 16293B, BHB 07-11, L1527, HH 212 and HH 111. We find, in particular, that the sight lines passing through the near-side of a highly inclined disk trace different temperature gradient directions than those through the far-side, which can lead to a polarization orientation on the near-side that is orthogonal to that on the far-side, and that the HH 111 disk may be such a case. Our technique for probing the disk temperature gradient through dust polarization can complement other methods, particularly those using molecular lines.

We present the properties of magnetic fields (B-fields) in two clumps (clump 1 and clump 2), located at the waist of the bipolar H II region Sh2-201, based on JCMT SCUBA-2/POL-2 observations of 850 $\mu$m polarized dust emission. We find that B-fields in the direction of the clumps are bent and compressed, showing bow-like morphologies, which we attribute to the feedback effect of the H II region on the surface of the clumps. Using the modified Davis-Chandrasekhar-Fermi method we estimate B-fields strengths of 266 $\mu$G and 65 $\mu$G for clump 1 and clump 2, respectively. From virial analyses and critical mass ratio estimates, we argue that clump 1 is gravitationally bound and could be undergoing collapse, whereas clump 2 is unbound and stable. We hypothesize that the interplay between thermal pressure imparted by the H II region, B-field morphologies, and the various internal pressures of the clumps (such as magnetic, turbulent, and gas thermal pressure), has the following consequences: (a) formation of clumps at the waist of the H II region; (b) progressive compression and enhancement of the B-fields in the clumps; (c) stronger B-fields will shield the clumps from erosion by the H II region and cause pressure equilibrium between the clumps and the H II region, thereby allowing expanding I-fronts to blow away from the filament ridge, forming bipolar H II regions; and (d) stronger B-fields and turbulence will be able to stabilize the clumps. A study of a larger sample of bipolar H II regions would help to determine whether our hypotheses are widely applicable.

We present detailed analysis of high-resolution C18O (2-1), SO (88-77), CO (3-2) and DCO+ (3-2) data obtained by the Atacama Large Millimeter/sub-millimeter Array (ALMA) towards a Class 0 Keplerian circumbinary disk around VLA1623A, which represents one of the most complete analysis towards a Class 0 source. From the dendrogram analysis, we identified several accretion flows feeding the circumbinary disk in a highly anisotropic manner. Stream-like SO emission around the circumbinary disk reveals the complicated shocks caused by the interactions between the disk, accretion flows and outflows. A wall-like structure is discovered south of VLA1623B. The discovery of two outflow cavity walls at the same position traveling at different velocities suggests the two outflows from both VLA1623A and VLA1623B overlays on top of each other in the plane of sky. Our detailed flat and flared disk modeling shows that Cycle 2 C18O J = 2-1 data is inconsistent with the combined binary mass of 0.2 Msun as suggested by early Cycle 0 studies. The combined binary mass for VLA1623A should be modified to 0.3 ~ 0.5 Msun.

The deuterium fractionation in starless cores gives us a clue to estimate their lifetime scales, thus allowing us to distinguish between different dynamical theories of core formation. Cores also seem to be subject to a differential N2 and CO depletion which was not expected from models. We aim to make a survey of 10 cores to estimate their lifetime scales and depletion profiles in detail. After L183, in Serpens, we present the second cloud of the series, L1512 in Auriga. To constrain the lifetime scale, we perform chemical modeling of the deuteration profiles across L1512 based on dust extinction measurements from near-infrared observations and non-local thermal equilibrium radiative transfer with multiple line observations of N2H+, N2D+, DCO+, C18O, and 13CO, plus H2D+ (1$_{10}$--1$_{11}$). We find a peak density of 1.1$\times$10$^5$ cm$^{-3}$ and a central temperature of 7.5$\pm$1 K, which are respectively higher and lower compared with previous dust emission studies. The depletion factors of N2H+ and N2D+ are 27$^{+17}_{-13}$ and 4$^{+2}_{-1}$ in L1512, intermediate between the two other more advanced and denser starless core cases, L183 and L1544. These factors also indicate a similar freeze-out of N2 in L1512, compared to the two others despite a peak density one to two orders of magnitude lower. Retrieving CO and N2 abundance profiles with the chemical model, we find that CO has a depletion factor of $\sim$430-870 and the N2 profile is similar to that of CO unlike towards L183. Therefore, L1512 has probably been living long enough so that N2 chemistry has reached steady state. N2H+ modeling remains compulsory to assess the precise physical conditions in the center of cold starless cores, rather than dust emission. L1512 is presumably older than 1.4 Myr. Therefore, the dominating core formation mechanism should be ambipolar diffusion for this source.

The IC5146 cloud is a nearby star-forming region in Cygnus, consisting of molecular gas filaments in a variety of evolutionary stages. We used optical and near-infrared polarization data toward the IC5146 cloud, reported in the first paper of this series, to reveal the magnetic fields in this cloud. Using the newly released $Gaia$ data, we found that the IC5146 cloud may contain two separate clouds: a first cloud, including the densest main filament at a distance of $\sim$600 pc, and a second cloud, associated with the Cocoon Nebula at a distance of $\sim$800 pc. The spatially averaged H-band polarization map revealed a well-ordered magnetic field morphology, with the polarization segments perpendicular to the main filament but parallel to the nearby sub-filaments, consistent with models assuming that the magnetic field is regulating cloud evolution. We estimated the magnetic field strength using the Davis-Chandrasekhar-Fermi method, and found that the magnetic field strength scales with volume density with a power-law index of $\sim0.5$ in the density range from $N_{H_2}\sim$ 10 to 3000 cm$^{-3}$, which indicates an anisotropic cloud contraction with a preferred direction along the magnetic field. In addition, the mass-to-flux ratio of the cloud gradually changes from subcritical to supercritical from the cloud envelope to the deep regions. These features are consistent with strong magnetic field star-formation models and suggest that the magnetic field is important in regulating the evolution of the IC5146 cloud.

We report high angular resolution observations of linearly polarized dust emission towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A (hereafter, IRAS4A) using the Karl G. Jansky Very Large Array (JVLA) at K (11.5-16.7 mm), Ka (8.1-10.3 mm), and Q bands (6.3-7.9 mm), and using the Atacama Large Millimeter Array (ALMA) at Band 6 (1.2 mm) and Band 7 (0.85-0.89 mm). On 100-1000 AU scales, all of these observations consistently trace the hourglass shaped magnetic field topology as shown in the previous studies. In the innermost 100 AU region of IRAS4A1, the polarization position angles (E-field) detected at 6.3-16.7 mm are consistent, however, are nearly 90 degrees offset from those detected at 1.2 mm and 0.85-0.89 mm. Such a 90 degree offset may be explained by that the inner $\sim$100 AU area is optically thick at wavelengths shorter than $\sim$1.5 mm, whereby the observations probe the absorption of aligned dust against the weakly or unpolarized warm dust emission from the innermost region. This can also consistently explain why the highest angular resolution ALMA images at Band 7 show that the polarization percentage increases with dust brightness temperature in the inner $\sim$100 AU region of IRAS4A1. Following this interpretation and assuming that the dust grains are aligned with the magnetic fields, the inferred magnetic field position angle based on the 90$^{\circ}$ rotated at 6.3-7.9 mm in the central peak of IRAS4A1 is $\sim-22^{\circ}$, which is approximately consistent with the outflow direction $\sim-9^{\circ}$.

We present the results of mapping observations toward an outflow-shocked region, OMC-2 FIR 4 using the Nobeyama 45-m telescope. We observed the area in $^{13}$CO ($J=1-0$), C$^{18}$O ($J=1-0$), N$_2$H$^+$ ($J=1-0$), CCS ($J_N=8_7-7_6$), HCO$^+$ ($J=1-0$), H$^{13}$CO$^+$ ($J=1-0$), HN$^{13}$C ($J=1-0$), H$^{13}$CN ($J=1-0$), DNC ($J=1-0$), N$_2$D$^+$ ($J=1-0$), and DC$_3$N ($J=9-8$). We detected a dense molecular clump that contains FIR 4/5. We also detected in $^{13}$CO blueshifted and redshifted components driven presumably by protostellar outflows in this region. The axes of the FIR 3 and VLA 13 outflows, projected on the plane of the sky, appear to point toward the FIR 4 clump, suggesting that the clump may be compressed by protostellar outflows from Class I sources, FIR 3 and VLA 13. Applying the hyperfine fit of N$_2$H$^+$ lines, we estimated the excitation temperature to be $\sim$ 20 K. The high excitation temperature is consistent with the fact that the clump contains protostars. The CCS emission was detected in this region for the first time. Its abundance is estimated to be a few $\times 10^{-12}$, indicating that the region is chemically evolved at $\sim 10^5$ years, which is comparable to the typical lifetime of the Class I protostars. This timescale is consistent with the scenario that star formation in FIR 4 is triggered by dynamical compression of the protostellar outflows. The [HNC]/[HCN] ratio is evaluated to be $\sim 0.5$ in the dense clump and the outflow lobes, whereas it is somewhat larger in the envelope of the dense clump. The small [HNC]/[HCN] ratio indicates that the HNC formation was prevented due to high temperatures. Such high temperatures seem to be consistent with the scenario that either protostellar radiation or outflow compression, or both, affected the thermal properties of this region.

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.

Here we present new ALMA observations of polarized dust emission from six of the most massive clumps in W43-Main. The clumps MM2, MM3, MM4, MM6, MM7, and MM8, have been resolved into two populations of fragmented filaments. From these two populations we extracted 81 cores (96 with the MM1 cores) with masses between 0.9 \Msun\ to 425 \Msun\ and a mass sensitivity of 0.08 M$_{\odot}$. The MM6, MM7, and MM8 clumps show significant fragmentation, but the polarized intensity appears to be sparse and compact. The MM2, MM3, and MM4 population shows less fragmentation, but with a single proto-stellar core dominating the emission at each clump. Also, the polarized intensity is more extended and significantly stronger in this population. From the polarized emission, we derived detailed magnetic field patterns throughout the filaments which we used to estimate field strengths for 4 out of the 6 clumps. The average field strengths estimations were found between 500 $\mu$G to 1.8 mG. Additionally, we detected and modeled infalling motions towards MM2 and MM3 from single dish HCO$^{+}(J=4 \rightarrow 3)$ and HCN$(J=4 \rightarrow 3)$ data resulting in mass infall rates of $\dot{\mathrm{M}}_{\mathrm{MM2}} = 1.2 \times 10^{-2}$ \Msun\ yr$^{-1}$ and $\dot{\mathrm{M}}_{\mathrm{MM3}} = 6.3 \times 10^{-3}$ \Msun\ yr$^{-1}$. By using our estimations, we evaluated the dynamical equilibrium of our cores by computing the total virial parameter $\alpha_{\mathrm{total}}$. For the cores with reliable field estimations, we found that 71\% of them appear to be gravitationally bound while the remaining 29\% are not. We concluded that these unbound cores, also less massive, are still accreting and have not yet reached a critical mass. This also implies different evolutionary time-scales, which essentially suggests that star-formation in high mass filaments is not uniform.

The dependence of polarization fraction $p$ on total intensity $I$ in polarized submillimeter emission measurements is typically parameterized as $p\propto I^{-\alpha}$ $(\alpha \leq 1)$, and used to infer dust grain alignment efficiency in star-forming regions, with an index $\alpha=1$ indicating near-total lack of alignment of grains with the magnetic field. In this work we demonstrate that the non-Gaussian noise characteristics of polarization fraction may produce apparent measurements of $\alpha \sim 1$ even in data with significant signal-to-noise in Stokes $Q$, $U$ and $I$ emission, and so with robust measurements of polarization angle. We present a simple model demonstrating this behavior, and propose a criterion by which well-characterized measurements of polarization fraction may be identified. We demonstrate that where our model is applicable, $\alpha$ can be recovered by fitting the $p-I$ relationship with the mean of the Rice distribution, without statistical debiasing of polarization fraction. We apply our model to JCMT BISTRO Survey POL-2 850$\mu$m observations of three clumps in the Ophiuchus Molecular Cloud, finding that in the externally-illuminated Oph A region, $\alpha\approx 0.34$, while in the more isolated Oph B and C, despite their differing star formation histories, $\alpha \sim 0.6-0.7$. Our results thus suggest that dust grain alignment in dense gas is more strongly influenced by incident interstellar radiation field than by star formation history. We further find that grains may remain aligned with the magnetic field at significantly higher gas densities than has previously been believed, thus allowing investigation of magnetic field properties within star-forming clumps and cores.

The HH 211 protostellar system is currently the youngest Class 0 system found with a rotating disk. We have mapped it at ~ 50 au (0.16") resolution, studying its magnetic field morphology with dust polarization in continuum at 232 and 358 GHz and its kinematics in C18O J=2-1 line. A flattened envelope extending out to ~ 400 au from the disk is detected in the continuum and C18O, slightly misaligned with the disk by 8 degree. It is spiraling inwards and expected to transform into a rotating disk at ~ 20 au, consistent with the disk radius estimated before. It appears to have a constant specific angular momentum and itself can result from an inside-out collapse of an extended envelope detected before in NH$_3$. In the flattened envelope, the polarization is mainly due to the magnetically aligned dust grains, inferring a highly pinched poloidal field morphology there. Thus, both the kinematics and field morphology support that the flattened envelope is a pseudodisk formed as the infalling gas is guided by the field lines to the equatorial plane. Interestingly, a point symmetric polarization distribution is also seen in the flattened envelope, implying that the pinched field lines also have a significant toroidal component generated by the rotation. No significant loss of angular momentum and thus no clear magnetic braking are detected in the flattened envelope around the disk probably because of the large misalignment between the axis of the rotation and the axis of the magnetic field in the cloud core.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

We present the results based on the optical $R$-band polarization observations of 280 stars distributed towards the dark globule LDN\,1225. {\it Gaia} data release 2 parallaxes along with the polarization data of $\sim$200 stars have been used to (a) constrain the distance of LDN\,1225 as 830$\pm$83~pc, (b) determine the contribution of interstellar polarization (ISP), and (c) characterize the dust properties and delineate the magnetic field (B-field) morphology of LDN\,1225. We find that B-fields are more organized and exhibit a small dispersion of 12$\degr$. Using the $^{12}$CO molecular line data from the Purple Mountain Observatory (PMO), along with the column density, dispersion in B-fields, we estimate B-field strength to be $\sim$56\,$\pm$\,10\,$\mu$G, magnetic to turbulence pressure to be $\sim$3\,$\pm$\,2, and the mass-to-magnetic flux ratio (in units of critical value) to be~$<$\,1. These results indicate the dominant role of B-fields in comparison to turbulence and gravity in rendering the cloud support. B-fields are aligned parallel to the low-density parts (traced by $^{12}$CO map) of the cloud, in contrast they are neither parallel nor perpendicular to the high-density core structures (traced by $^{13}$CO and C$^{18}$O maps). LDN\,1225 hosts two 70\,$\mu$m sources which seem to be of low-mass Class 0 sources. The total-to-selective extinction derived using optical and near-infrared photometric data is found to be anomalous ($R_{V}$~$=$~3.4), suggesting dust grain growth in LDN\,1225. Polarization efficiency of dust grains follows a power-law index of $-$0.7 inferring that optical polarimetry traces B-fields in the outer parts of the cloud.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

Probing magnetic fields in self-gravitating molecular clouds are generally difficult even with the use of the polarimetry. Based on the properties of magneto-hydrodynamic (MHD) turbulence and turbulent reconnection, Velocity Gradient Technique (VGT) provides a new way in tracing magnetic field orientation and strength based on the spectroscopic data. Our study tests the applicability of VGT in various molecular tracers, e.g. 12CO, 13CO, and C18O. By inspecting synthetic molecular line maps of CO isotopologue generated through radiative transfer calculations, we show that the VGT method can be successfully applied in probing the magnetic field direction in the diffuse interstellar medium as well as in self-gravitating molecular clouds.

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.

We report Submillimeter Array molecular line observations in the 345 GHz band of five massive dense cores, Cyg-N38, Cyg-N43, Cyg-N48, Cyg-N51, and Cyg-N53 in the DR21 filament. The molecular line data reveal several dynamical features of the cores: (1) prominent outflows in all cores seen in the CO and SiO lines, (2) significant velocity gradients in Cyg-N43 and Cyg-N48 seen in the H13CN and H13CO+ lines suggesting 0.1-pc-scale rotational motions, and (3) possible infalls in Cyg-N48 found in the SiO and SO lines. Comparing the molecular line data and our dust polarization data in Ching et al. (2017), we find that the gradients of line-of-sight velocities appear to be randomly oriented relative to the plane-of-sky magnetic fields. Our simulations suggest that this random alignment implies parallel or random alignment between the velocity gradients and magnetic fields in the three dimensional space. The linewidths of H13CN emission are consistently wider than those of H13CO+ emission in the 3''-10'' detectable scales, which can be explained by the existence of ambipolar diffusion with maximum plane-of-the-sky magnetic field strengths of 1.9 mG and 5.1 mG in Cyg-N38 and Cyg-N48, respectively. Our results suggest that the gas dynamics may distort the magnetic fields of the cores of into complex structures and ambipolar diffusion could be important in dissipating the magnetic energies of the cores.

VLA 1623$-$2417 is a triple protostellar system deeply embedded in Ophiuchus A. Sources A and B have a separation of 1.1", making their study difficult beyond the submillimeter regime. Lack of circumstellar gas emission suggested that VLA 1623$-$2417 B has a very cold envelope and is much younger than source A, generally considered the prototypical Class 0 source. We explore the consequences of new ALMA Band 9 data on the spectral energy distribution (SED) of VLA 1623$-$2417 and their inferred nature. Using dust continuum observations spanning from centimeter to near-infrared wavelengths, the SED of each component in VLA 1623$-$2417 is constructed and analysed. The ALMA Band 9 data presented here show that the SED of VLA 1623$-$2417 B does not peak at 850 $\mu$m as previously expected, but instead presents the same shape as VLA 1623$-$2417 A at wavelengths shorter than 450 $\mu$m. The results presented here indicate that the previous assumption that the flux in $Herschel$ and Spitzer observations is solely dominated by VLA 1623$-$2417 A is not valid, and instead, VLA 1623$-$2417 B most likely contributes a significant fraction of the flux at $\lambda~<$ 450 $\mu$m. These results, however, do not explain the lack of circumstellar gas emission and puzzling nature of VLA 1623$-$2417 B.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present the first high-resolution, submillimeter-wavelength polarimetric observations of -- and thus direct observations of the magnetic field morphology within -- the dense gas of the Pillars of Creation in M16. These 850$\,\mu$m observations, taken as part of the BISTRO (B-Fields in Star-forming Region Observations) Survey using the POL-2 polarimeter on the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT), show that the magnetic field runs along the length of the pillars, perpendicular to, and decoupled from, the field in the surrounding photoionized cloud. Using the Chandrasekhar-Fermi method we estimate a plane-of-sky magnetic field strength of $170-320\,\mu$G in the Pillars, consistent with their having been formed through compression of gas with initially weak magnetization. The observed magnetic field strength and morphology suggests that the magnetic field may be slowing the pillars' evolution into cometary globules. We thus hypothesize that the evolution and lifetime of the Pillars may be strongly influenced by the strength of the coupling of their magnetic field to that of their parent photoionized cloud -- i.e. that the Pillars' longevity results from magnetic support

We performed the new JVLA full polarization observations at 40-48 GHz (6.3-7.5 mm) towards the nearby ($d$ $=$147$\pm$3.4 pc) Class 0 YSO IRAS 16293-2422, and compare with the previous SMA observations reported by Rao et al. (2009; 2014). We observed the quasar J1407+2827 which is weakly polarized and can be used as a leakage term calibrator for $<$9 GHz observations, to gauge the potential residual polarization leakage after calibration. We did not detect Stokes Q, U, and V intensities from the observations of J1407+2827, and constrain (3-$\sigma$) the residual polarization leakage after calibration to be $\lesssim$0.3\%. We detect linear polarization from one of the two binary components of our target source, IRAS\,16293-2422\,B. The derived polarization position angles from our observations are in excellent agreement with those detected from the previous observations of the SMA, implying that on the spatial scale we are probing ($\sim$50-1000 au), the physical mechanisms for polarizing the continuum emission do not vary significantly over the wavelength range of $\sim$0.88-7.5 mm. We hypothesize that the observed polarization position angles trace the magnetic field which converges from large scale to an approximately face-on rotating accretion flow. In this scenario, magnetic field is predominantly poloidal on $>$100 au scales, and becomes toroidal on smaller scales. However, this interpretation remains uncertain due to the high dust optical depths at the central region of IRAS\,16293-2422\,B and the uncertain temperature profile. We suggest that dust polarization at wavelengths comparable or longer than 7\,mm may still trace interstellar magnetic field. Future sensitive observations of dust polarization in the fully optically thin regime will have paramount importance for unambiguously resolving the magnetic field configuration.

We report results from analyzing the B612 asteroid observation data taken by the sCMOS cameras on board of Planet SkySat-3 using the synthetic tracking technique. The analysis demonstrates the expected sensitivity improvement in the signal-to-noise ratio of the asteroids from properly stacking up the the short exposure images in post-processing.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

Although, from a theoretical point of view, magnetic fields are believed to have a significant role during the early stages of star formation, especially during the main accretion phase, the magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. We carried out observations of the polarized dust continuum emission with the SMA interferometer at 0.87mm, in order to probe the structure of the magnetic field in a sample of 12 low-mass Class 0 envelopes, including both single protostars and multiple systems, in nearby clouds. Our SMA observations probe the envelope emission at scales of 600-5000 au with a spatial resolution ranging from 600 to 1500 au depending on the source distance. We report the detection of linearly polarized dust continuum emission in all of our targets, with average polarization fractions ranging from 2% to 10% in these protostellar envelopes. The polarization fraction decreases with the continuum flux density, which translates into a decrease with the H2 column density within an individual envelope. Our analysis show that the envelope-scale magnetic field is preferentially observed either aligned or perpendicular to the outflow direction. Interestingly, our results suggest for the first time a relation between the orientation of the magnetic field and the rotational energy of envelopes, with a larger occurrence of misalignment in sources where strong rotational motions are detected at hundreds to thousands of au scales. We also show that the best agreement between the magnetic field and outflow orientation is found in sources showing no small-scale multiplicity and no large disks at about 100 au scales.

The role of the magnetic field during protostellar collapse is poorly constrained from an observational point of view, although it could be significant if we believe state-of-the-art models of protostellar formation. We present polarimetric observations of the 233 GHz thermal dust continuum emission obtained with ALMA in the B335 Class 0 protostar. Linearly polarized dust emission arising from the circumstellar material in the envelope of B335 is detected at all scales probed by our observations, from radii of 50 to 1000 au. The magnetic field structure producing the dust polarization has a very ordered topology in the inner envelope, with a transition from a large-scale poloidal magnetic field, in the outflow direction, to strongly pinched in the equatorial direction. This is probably due to magnetic field lines being dragged along the dominating infall direction since B335 does not exhibit prominent rotation. Our data and their qualitative comparison to a family of magnetized protostellar collapse models show that, during the magnetized collapse in B335, the magnetic field is maintaining a high level of organization from scales 1000 au to 50 au: this suggests the field is dynamically relevant and capable of influencing the typical outcome of protostellar collapse, such as regulating the disk size in B335.

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.

Polarized emission is detected in two young nearly edge-on protostellar disks in 343 GHz continuum at ~ 50 au (~ 0.12") resolution with Atacama Large Millimeter/submillimeter Array. One disk is in HH 212 (Class 0) and the other in HH 111 (early Class I) protostellar system. Polarization fraction is ~ 1%. The disk in HH 212 has a radius of ~ 60 au. The emission is mainly detected from the nearside of the disk. The polarization orientations are almost perpendicular to the disk major axis, consistent with either self-scattering or emission by grains aligned with a poloidal field around the outer edge of the disk because of optical depth effect and temperature gradient; the presence of a poloidal field would facilitate the launching of a disk wind, for which there is already tentative evidence in the same source. The disk of HH 111 VLA 1 has a larger radius of ~ 220 au and is thus more resolved. The polarization orientations are almost perpendicular to the disk major axis in the nearside, but more along the major axis in the farside, forming roughly half of an elliptical pattern there. It appears that toroidal and poloidal magnetic field may explain the polarization on the near and far side of the disk, respectively. However, it is also possible that the polarization is due to self-scattering. In addition, alignment of dust grains by radiation flux may play a role in the farside. Our observations reveal a diversity of disk polarization patterns that should be taken into account in future modeling efforts.

We present optical and near-infrared stellar polarization observations toward the dark filamentary clouds associated with IC5146. The data allow us to investigate the dust properties (this paper) and the magnetic field structure (Paper II). A total of 2022 background stars were detected in $R_{c}$-, $i'$-, $H$-, and/or $K$-bands to $A_V \lesssim 25$ mag. The ratio of the polarization percentage at different wavelengths provides an estimate of $\lambda_{max}$, the wavelength of peak polarization, which is an indicator of the small-size cutoff of the grain size distribution. The grain size distribution seems to significantly change at $A_V \sim$ 3 mag, where both the average and dispersion of $P_{R_c}/P_{H}$ decrease. In addition, we found $\lambda_{max}$ $\sim$ 0.6-0.9 $\mu$m for $A_V>2.5$ mag, which is larger than the $\sim$ 0.55 $\mu$m in the general ISM, suggesting that grain growth has already started in low $A_V$ regions. Our data also reveal that polarization efficiency (PE $\equiv P_{\lambda}/A_V$) decreases with $A_V$ as a power-law in $R_c$-, $i'$-, and $K$-bands with indices of -0.71$\pm$0.10, -1.23$\pm$0.10 and -0.53$\pm$0.09. However, $H$-band data show a power index change; the PE varies with $A_V$ steeply (index of -0.95$\pm$0.30) when $A_V < 2.88\pm0.67$ mag but softly (index of -0.25$\pm$0.06) for greater $A_V$ values. The soft decay of PE in high $A_V$ regions is consistent with the Radiative Aligned Torque model, suggesting that our data trace the magnetic field to $A_V \sim 20$ mag. Furthermore, the breakpoint found in $H$-band is similar to the $A_V$ where we found the $P_{R_c}/P_{H}$ dispersion significantly decreased. Therefore, the flat PE-$A_V$ in high $A_V$ regions implies that the power index changes result from additional grain growth.

The influence of magnetic fields (B-fields) in the formation and evolution of bipolar bubbles, due to the expanding ionization fronts (I-fronts) driven by the Hii regions that are formed and embedded in filamentary molecular clouds, has not been well-studied yet. In addition to the anisotropic expansion of I-fronts into a filament, B-fields are expected to introduce an additional anisotropic pressure which might favor expansion and propagation of I-fronts to form a bipolar bubble. We present results based on near-infrared polarimetric observations towards the central $\sim$8'$\times$8' area of the star-forming region RCW57A which hosts an Hii region, a filament, and a bipolar bubble. Polarization measurements of 178 reddened background stars, out of the 919 detected sources in the JHKs-bands, reveal B-fields that thread perpendicular to the filament long axis. The B-fields exhibit an hour-glass morphology that closely follows the structure of the bipolar bubble. The mean B-field strength, estimated using the Chandrasekhar-Fermi method, is 91$\pm$8 {\mu}G. B-field pressure dominates over turbulent and thermal pressures. Thermal pressure might act in the same orientation as those of B-fields to accelerate the expansion of those I-fronts. The observed morphological correspondence among the B-fields, filament, and bipolar bubble demonstrate that the B-fields are important to the cloud contraction that formed the filament, gravitational collapse and star formation in it, and in feedback processes. The latter include the formation and evolution of mid-infrared bubbles by means of B-field supported propagation and expansion of I-fronts. These may shed light on preexisting conditions favoring the formation of the massive stellar cluster in RCW57A.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

We determine the magnetic field strength in the OMC 1 region of the Orion A filament via a new implementation of the Chandrasekhar-Fermi method using observations performed as part of the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-Forming Region Observations (BISTRO) survey with the POL-2 instrument. We combine BISTRO data with archival SCUBA-2 and HARP observations to find a plane-of-sky magnetic field strength in OMC 1 of $B_{\rm pos}=6.6\pm4.7$ mG, where $\delta B_{\rm pos}=4.7$ mG represents a predominantly systematic uncertainty. We develop a new method for measuring angular dispersion, analogous to unsharp masking. We find a magnetic energy density of $\sim1.7\times 10^{-7}$ Jm$^{-3}$ in OMC 1, comparable both to the gravitational potential energy density of OMC 1 ($\sim 10^{-7}$ Jm$^{-3}$), and to the energy density in the Orion BN/KL outflow ($\sim 10^{-7}$ Jm$^{-3}$). We find that neither the Alfv\'{e}n velocity in OMC 1 nor the velocity of the super-Alfv\'{e}nic outflow ejecta is sufficiently large for the BN/KL outflow to have caused large-scale distortion of the local magnetic field in the $\sim$500-year lifetime of the outflow. Hence, we propose that the hour-glass field morphology in OMC 1 is caused by the distortion of a primordial cylindrically-symmetric magnetic field by the gravitational fragmentation of the filament and/or the gravitational interaction of the BN/KL and S clumps. We find that OMC 1 is currently in or near magnetically-supported equilibrium, and that the current large-scale morphology of the BN/KL outflow is regulated by the geometry of the magnetic field in OMC 1, and not vice versa.

We present high angular resolution dust polarization and molecular line observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostar Serpens SMM1. By complementing these observations with new polarization observations from the Submillimeter Array (SMA) and archival data from the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the James Clerk Maxwell Telescopes (JCMT), we can compare the magnetic field orientations at different spatial scales. We find major changes in the magnetic field orientation between large (~0.1 pc) scales -- where the magnetic field is oriented E-W, perpendicular to the major axis of the dusty filament where SMM1 is embedded -- and the intermediate and small scales probed by CARMA (~1000 AU resolution), the SMA (~350 AU resolution), and ALMA (~140 AU resolution). The ALMA maps reveal that the redshifted lobe of the bipolar outflow is shaping the magnetic field in SMM1 on the southeast side of the source; however, on the northwestern side and elsewhere in the source, low velocity shocks may be causing the observed chaotic magnetic field pattern. High-spatial-resolution continuum and spectral-line observations also reveal a tight (~130 AU) protobinary system in SMM1-b, the eastern component of which is launching an extremely high-velocity, one-sided jet visible in both CO(2-1) and SiO(5-4); however, that jet does not appear to be shaping the magnetic field. These observations show that with the sensitivity and resolution of ALMA, we can now begin to understand the role that feedback (e.g., from protostellar outflows) plays in shaping the magnetic field in very young, star-forming sources like SMM1.

We report ALMA observations of polarized dust emission from the protostellar source Ser-emb 8 at a linear resolution of 140 au. Assuming models of dust-grain alignment hold, the observed polarization pattern gives a projected view of the magnetic field structure in this source. Contrary to expectations based on models of strongly magnetized star formation, the magnetic field in Ser-emb 8 does not exhibit an hourglass morphology. Combining the new ALMA data with previous observational studies, we can connect magnetic field structure from protostellar core (~80,000 au) to disk (~100 au) scales. We compare our observations with four magnetohydrodynamic gravo-turbulence simulations made with the AREPO code that have initial conditions ranging from super-Alfv\'enic (weakly magnetized) to sub-Alfv\'enic (strongly magnetized). These simulations achieve the spatial dynamic range necessary to resolve the collapse of protostars from the parsec scale of star-forming clouds down to the ~100 au scale probed by ALMA. Only in the very strongly magnetized simulation do we see both the preservation of the field direction from cloud to disk scales and an hourglass-shaped field at < 1000 au scales. We conduct an analysis of the relative orientation of the magnetic field and the density structure in both the Ser-emb 8 ALMA observations and the synthetic observations of the four AREPO simulations. We conclude that the Ser-emb 8 data are most similar to the weakly magnetized simulations, which exhibit random alignment, in contrast to the strongly magnetized simulation, where the magnetic field plays a role in shaping the density structure in the source. In the weak-field case, it is turbulence -- not the magnetic field -- that shapes the material that forms the protostar, highlighting the dominant role that turbulence can play across many orders of magnitude in spatial scale.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We present Submillimeter Array 880 $\mu$m dust polarization observations of six massive dense cores in the DR21 filament. The dust polarization shows complex magnetic field structures in the massive dense cores with sizes of 0.1 pc, in contrast to the ordered magnetic fields of the parsec-scale filament. The major axes of the massive dense cores appear to be aligned either parallel or perpendicular to the magnetic fields of the filament, indicating that the parsec-scale magnetic fields play an important role in the formation of the massive dense cores. However, the correlation between the major axes of the cores and the magnetic fields of the cores is less significant, suggesting that during the core formation, the magnetic fields below 0.1 pc scales become less important than the magnetic fields above 0.1 pc scales in supporting a core against gravity. Our analysis of the angular dispersion functions of the observed polarization segments yields the plane-of-sky magnetic field strengths of 0.4--1.7 mG of the massive dense cores. We estimate the kinematic, magnetic, and gravitational virial parameters of the filament and the cores. The virial parameters show that in the filament, the gravitational energy is dominant over magnetic and kinematic energies, while in the cores, the kinematic energy is dominant. Our work suggests that although magnetic fields may play an important role in a collapsing filament, the kinematics arising from gravitational collapse must become more important than magnetic fields during the evolution from filaments to massive dense cores.

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.

We present the first resolved observations of the 1.3mm polarized emission from the disk-like structure surrounding the high-mass protostar Cepheus A HW2. These CARMA data partially resolve the dust polarization, suggesting an uniform morphology of polarization vectors with an average position angle of 57 degrees and an average polarization fraction of 2.0%. The distribution of the polarization vectors can be attributed to (1) the direct emission of magnetically aligned grains of dust by a uniform magnetic field, or (2) the pattern produced by the scattering of an inclined disk. We show that both models can explain the observations, and perhaps a combination of the two mechanisms produce the polarized emission. A third model including a toroidal magnetic field does not match the observations. Assuming scattering is the polarization mechanism, these observations suggest that during the first few 10000 years of high-mass star formation, grain sizes can grow from 1 to several 10s micron.

Here we present the first results from ALMA observations of 1 mm polarized dust emission towards the W43-MM1 high mass star forming clump. We have detected a highly fragmented filament with source masses ranging from 14Msun to 312Msun, where the largest fragment, source A, is believed to be one of the most massive in our Galaxy. We found a smooth, ordered, and detailed polarization pattern throughout the filament which we used to derived magnetic field morphologies and strengths for 12 out of the 15 fragments detected ranging from 0.2 to 9 mG. The dynamical equilibrium of each fragment was evaluated finding that all the fragments are in a super-critical state which is consistent with previously detected infalling motions towards W43-MM1. Moreover, there are indications suggesting that the field is being dragged by gravity as the whole filament is collapsing.

L1551 is chosen because it is relatively isolated in the Taurus molecular cloud shielded from FUV photons, providing an ideal environment for studying the target properties. Our observations cover ~40'x40' with resolution ~30", which are the maps with highest spatial dynamical range to date. We derive the X(13CO)/X(C18O) value on the sub-parsec scales in the range of ~3-27 with a mean value of 8.0+-2.8. Comparing to the visual extinction map derived from the Herschel observations, we found that the abundance ratio reaches its maximum at low Av (i.e., Av ~ 1-4mag), and decreases to the typical solar system value of 5.5 inside L1551 MC. The high X(13CO)/X(C18O) value at the boundary of the cloud is most likely due to the selective FUV photodissociation of C18O. This is in contrast with Orion-A where its internal OB stars keep the abundance ratio at a high level greater than ~10. In addition, we explore the variation of the X-factor, because it is an uncertain but widely used quantity in extragalactic studies. We found that X-factor is proportional to N(H2)^1.0 which is consistent with previous simulations. Excluding the high density region, the average X-factor is similar to the Milky Way average value.

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.

(Abridged) The physical structure of deeply-embedded low-mass protostars (Class 0) on scales of less than 300 AU is still poorly constrained. Determining this is crucial for understanding the physical and chemical evolution from cores to disks. In this study two models of the emission, a Gaussian disk intensity distribution and a parametrized power-law disk model, are fitted to sub-arcsecond resolution interferometric continuum observations of five Class 0 sources, including one source with a confirmed Keplerian disk. For reference, a spherically symmetric single power-law envelope is fitted to the larger scale ($\sim$1000 AU) emission and investigated further for one of the sources on smaller scales. A thin disk model can approximate the emission and physical structure in the inner few 100 AU scales of the studied deeply-embedded low-mass protostars and paves the way for analysis of a larger sample with ALMA. While the disk radii agree with previous estimates the masses are different for some of the sources studied. Assuming a typical temperature distribution, the fractional amount of mass in the disk above 100 K varies in between 7% to 30%. Kinematic data are needed to determine the presence of any Keplerian disk. Using previous observations of p-H$_2^{18}$O, we estimate the relative gas phase water abundances roughly an order of magnitude higher than previously inferred when both warm and cold H$_2$ was used as reference. A spherically symmetric single power-law envelope model fails to simultaneously reproduce both the small and large scale emission.

We report new Karl G. Jansky Very Large Array (JVLA), 0$"$.5 angular resolution observations of linearly polarized continuum emission at 6.9 mm, towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A. This target source is a collapsing dense molecular core, which was resolved at short wavelengths to have hourglass shaped B-field configuration. We compare these 6.9 mm observations with previous polarization Submillimeter Array (SMA) observations at 0.88 mm, which have comparable angular resolution ($\sim$0$"$7). We found that at the same resolution, the observed polarization position angles at 6.9 mm are slightly deviated from those observed at 0.88 mm. Due to the lower optical depth of the emission at 6.9 mm, and the potential effect of dust grain growth, the new JVLA observations are likely probing B-field alignments in regions interior to those sampled by the previous polarization observations at higher frequencies. Our understanding can be improved by more sensitive observations, and observations for the more extended spatial scales.

We present Submillimeter Array (SMA) polarization observations of the CO $J$ = 3--2 line toward the NGC1333 IRAS 4A. The CO Stokes $I$ maps at an angular resolution of $\sim$1$\arcsec$ reveal two bipolar outflows from the binary sources of the NGC 1333 IRAS 4A. The kinematic features of the CO emission can be modeled by wind-driven outflows at $\sim$ 20$\arcdeg$ inclined from the plane of the sky. Close to the protostars the CO polarization, at an angular resolution of $\sim$$2\farcs3$, has a position angle approximately parallel to the magnetic field direction inferred from the dust polarizations. The CO polarization direction appears to vary smoothly from an hourglass field around the core to an arc-like morphology wrapping around the outflow, suggesting a helical structure of magnetic fields that inherits the poloidal fields at the launching point and consists of toroidal fields at a farther distance of outflow. The helical magnetic field is consistent with the theoretical expectations for launching and collimating outflows from a magnetized rotating disk. Considering that the CO polarized emission is mainly contributed from the low-velocity and low-resolution data, the helical magnetic field is likely a product of the wind-envelope interaction in the wind-driven outflows. The CO data reveal a PA of $\sim$ 30$\arcdeg$ deflection in the outflows. The variation in the CO polarization angle seems to correlate with the deflections. We speculate that the helical magnetic field contributes to $\sim$ 10$\arcdeg$ deflection of the outflows by means of Lorenz force.

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $\lambda$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

Fomalhaut plays an important role in the study of debris disks and small bodies in other planetary systems. The proximity and luminosity of the star make key features of its debris, like the water ice-line, accessible. Here we present ALMA cycle 1, 870 \mu m (345 GHz) observations targeted at the inner part of the Fomalhaut system with a synthesized beam of 0.45"x0.37" (~3 AU linear resolution at the distance of Fomalhaut) and a rms of 26 \mu Jy/beam. The high angular resolution and sensitivity of the ALMA data enable us to place strong constraints on the nature of the warm excess revealed by Spitzer and Herschel observations. We detect a point source at the star position with a total flux consistent with thermal emission from the stellar photosphere. No structures that are brighter than 3\sigma\ are detected in the central 15 AU x 15 AU region. Modeling the spectral energy distribution using parameters expected for a dust-producing planetesimal belt indicates a radial location in the range ~8-15 AU. This is consistent with the location where ice sublimates in Fomalhaut, i.e., an asteroid-belt analog. The 3\sigma\ upper limit for such a belt is <1.3 mJy at 870 \mu m. We also interpret the 2 and 8-13 \mu m interferometric measurements to reveal the structure in the inner 10 AU region as dust naturally connected to this proposed asteroid belt by Poynting-Robertson drag, dust sublimation, and magnetically trapped nano grains.

We highlight distinct and systematic observational features of magnetic field morphologies in polarized submm dust continuum. We illustrate this with specific examples and show statistical trends from a sample of 50 star-forming regions.

Context: High-resolution studies of class 0 protostars represent the key to constraining protostar formation models. VLA16234-2417 represents the prototype of class 0 protostars, and it has been recently identified as a triple non-coeval system. Aim: We aim at deriving the physical properties of the jets in VLA16234-2417 using tracers of shocked gas. Methods: ALMA Cycle 0 Early Science observations of CO(2-1) in the extended configuration are presented in comparison with previous SMA CO(3-2) and Herschel-PACS [OI}] 63 micron observations. Gas morphology and kinematics were analysed to constrain the physical structure and origin of the protostellar outflows. Results: We reveal a collimated jet component associated with the [OI] 63 micron emission at about 8'' (about 960 AU) from source B. This newly detected jet component is inversely oriented with respect to the large-scale outflow driven by source A, and it is aligned with compact and fast jet emission very close to source B (about 0.3'') rather than with the direction perpendicular to the A disk. We also detect a cavity-like structure at low projected velocities, which surrounds the [OI] 63 micron emission and is possibly associated with the outflow driven by source A. Finally, no compact outflow emission is associated with source W. Conclusions: Our high-resolution ALMA observations seem to suggest there is a fast and collimated jet component associated with source B. This scenario would confirm that source B is younger than A, that it is in a very early stage of evolution, and that it drives a faster, more collimated, and more compact jet with respect to the large-scale slower outflow driven by A. However, a different scenario of a precessing jet driven by A cannot be firmly excluded from the present observations.

The Spitzer Space Telescope mapped the Perseus molecular cloud complex with IRAC and MIPS as part of the c2d Spitzer Legacy project. This paper combines the observations from both instruments giving an overview of low-mass star formation across Perseus from 3.6 to 70 micron. We provide an updated list of young stellar objects with new classifications and source fluxes from previous works, identifying 369 YSOs in Perseus with the Spitzer dataset. By synthesizing the IRAC and MIPS maps of Perseus and building on the work of previous papers in this series (Jorgensen et al. 2006, Rebull et al. 2007), we present a current census of star formation across the cloud and within smaller regions. 67% of the YSOs are associated with the young clusters NGC 1333 and IC 348. The majority of the star formation activity in Perseus occurs in the regions around the clusters, to the eastern and western ends of the cloud complex. The middle of the cloud is nearly empty of YSOs despite containing regions of high visual extinction. The western half of Perseus contains three-quarters of the total number of embedded YSOs (Class 0+I and Flat SED sources) in the cloud and nearly as many embedded YSOs as Class II and III sources. Class II and III greatly outnumber Class 0+I objects in eastern Perseus and IC 348. These results are consistent with previous age estimates for the clusters. Across the cloud, 56% of YSOs and 91% of the Class 0+I and Flat sources are in areas where Av > 5 mag, indicating a possible extinction threshold for star formation.

Due to instrumental limitations and a lack of disk detections, the structure between the envelope and the rotationally supported disk has been poorly studied. This is now possible with ALMA through observations of CO isotopologs and tracers of freezeout. Class 0 sources are ideal for such studies given their almost intact envelope and young disk. The structure of the disk-envelope interface of the prototypical Class 0 source, VLA1623A which has a confirmed Keplerian disk, is constrained from ALMA observations of DCO+ 3-2 and C18O 2-1. The physical structure of VLA1623 is obtained from the large-scale SED and continuum radiative transfer. An analytic model using a simple network coupled with radial density and temperature profiles is used as input for a 2D line radiative transfer calculation for comparison with the ALMA Cycle 0 12m array and Cycle 2 ACA observations of VLA1623. DCO+ emission shows a clumpy structure bordering VLA1623A's Keplerian disk, suggesting a cold ring-like structure at the disk-envelope interface. The radial position of the observed DCO+ peak is reproduced in our model only if the region's temperature is between 11-16K, lower than expected from models constrained by continuum and SED. Altering the density has little effect on the DCO+ position, but increased density is needed to reproduce the disk traced in C18O. The DCO+ emission around VLA1623A is the product of shadowing of the envelope by the disk. Disk-shadowing causes a drop in the gas temperature outside of the disk on >200AU scales, encouraging deuterated molecule production. This indicates that the physical structure of the disk-envelope interface differs from the rest of the envelope, highlighting the drastic impact that the disk has on the envelope and temperature structure. The results presented here show that DCO+ is an excellent cold temperature tracer.

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).

Submillimeter dust polarization measurements of a sample of 50 star-forming regions, observed with the SMA and the CSO covering pc-scale clouds to mpc-scale cores, are analyzed in order to quantify the magnetic field importance. The magnetic field misalignment $\delta$ -- the local angle between magnetic field and dust emission gradient -- is found to be a prime observable, revealing distinct distributions for sources where the magnetic field is preferentially aligned with or perpendicular to the source minor axis. Source-averaged misalignment angles $\langle|\delta|\rangle$ fall into systematically different ranges, reflecting the different source-magnetic field configurations. Possible bimodal $\langle|\delta|\rangle$-distributions are found for the separate SMA and CSO samples. Combining both samples broadens the distribution with a wide maximum peak at small $\langle|\delta|\rangle$-values. Assuming the 50 sources to be representative, the prevailing source-magnetic field configuration is one that statistically prefers small magnetic field misalignments $|\delta|$. When interpreting $|\delta|$ together with an MHD force equation, as developed in the framework of the polarization-intensity gradient method, a sample-based scaling fits the magnetic field tension-to-gravity force ratio $\langle\Sigma_B\rangle$ versus $\langle|\delta|\rangle$ with $\langle\Sigma_B\rangle = 0.116 \cdot \exp(0.047\cdot \langle|\delta|\rangle)\pm 0.20$ (mean error), providing a way to estimate the relative importance of the magnetic field, only based on measurable field misalignments $|\delta|$. The force ratio $\Sigma_B$ discriminates systems that are collapsible on average ($\langle \Sigma_B\rangle <1$) from other molecular clouds where the magnetic field still provides enough resistance against gravitational collapse ($\langle \Sigma_B\rangle >1$) (abridged).

HH 211 is a young Class 0 protostellar system, with a flattened envelope, a possible rotating disk, and a collimated jet. We have mapped it with the Submillimeter Array in 341.6 GHz continuum and SiO J=8-7 at ~ 0.6 resolution. The continuum traces the thermal dust emission in the flattened envelope and the possible disk. Linear polarization is detected in the continuum in the flattened envelope. The field lines implied from the polarization have different orientations, but they are not incompatible with current gravitational collapse models, which predict different orientation depending on the region/distance. Also, we might have detected for the first time polarized SiO line emission in the jet due to the Goldreich-Kylafis effect. Observations at higher sensitivity are needed to determine the field morphology in the jet.

Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $\mu$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AU), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation.

We present subarcsecond resolution polarimetric observations of the 878 mum thermal dust continuum emission obtained with the Submillimeter Array (SMA) towards the IRAS 16293-2422 protostellar binary system. We report the detection of linearly polarized dust emission arising from the circumstellar disk associated with the IRAS 16293-2422 B protostar. The fractional polarization of ~1.4$% is only slightly lower than that expected from theoretical calculations in such disks. The magnetic field structure on the plane of the sky derived from the dust polarization suggests a complex magnetic field geometry in the disk, possibly associated with a rotating disk that is wrapping the field lines as expected from the simulations. The polarization around IRAS 16293-2422 A at subarcsecond angular resolution is only marginally detected.

Context: Rotationally supported disks are critical in the star formation process. The questions of when do they form and what factors influence or hinder their formation have been studied but are largely unanswered. Observations of early stage YSOs are needed to probe disk formation. Aims: VLA1623 is a triple non-coeval protostellar system, with a weak magnetic field perpendicular to the outflow, whose Class 0 component, VLA1623A, shows a disk-like structure in continuum with signatures of rotation in line emission. We aim to determine whether this structure is in part or in whole a rotationally supported disk, i.e. a Keplerian disk, and what are its characteristics. Methods: ALMA Cycle 0 Early Science 1.3 mm continuum and C$^{18}$O (2-1) observations in the extended configuration are presented here and used to perform an analysis of the disk-like structure using PV diagrams and thin disk modelling with the addition of foreground absorption. Results: The PV diagrams of the C$^{18}$O line emission suggest the presence of a rotationally supported component with a radius of at least 50 AU. Kinematical modelling of the line emission shows that the disk out to 180 AU is actually rotationally supported, with the rotation being well described by Keplerian rotation out to at least 150 AU, and the central source mass to be $\sim$0.2 M$_{sun}$ for an inclination of 55$^{\circ}$. Pure infall and conserved angular momentum rotation models are excluded. Conclusions: VLA1623A, a very young Class 0 source, presents a disk with an outer radius $R_{\rm out}$ = 180 AU with a Keplerian velocity structure out to at least 150 AU. The weak magnetic fields and recent fragmentation in this region of rho Ophiuchus may have played a lead role in the formation of the disk.

We present the various science cases for building Band 1 receivers as part of ALMA's ongoing Development Program. We describe the new frequency range for Band 1 of 35-52 GHz, a range chosen to maximize the receiver suite's scientific impact. We first describe two key science drivers: 1) the evolution of grains in protoplanetary disks and debris disks, and 2) molecular gas in galaxies during the era of re-ionization. Studies of these topics with Band 1 receivers will significantly expand ALMA's Level 1 Science Goals. In addition, we describe a host of other exciting continuum and line science cases that require ALMA's high sensitivity and angular resolution. For example, ALMA Band 1 continuum data will probe the Sunyaev-Zel'dovich Effect in galaxy clusters, Very Small Grains and spinning dust, ionized jets from young stars, spatial and flaring studies of Sgr A*, the acceleration sites of solar flares, pulsar wind nebulae, radio supernovae, and X-ray binaries. Furthermore, ALMA Band 1 line data will probe chemical differentiation in cloud cores, complex carbon chain molecules, extragalactic radio recombination lines, masers, magnetic fields through Zeeman effect measurements, molecular outflows from young stars, the co-evolution of star formation and active galactic nuclei, and the molecular content of galaxies at z ~ 3. ALMA provides similar to better sensitivities than the JVLA over 35-50 GHz, with differences increasing with frequency. ALMA's smaller antennas and shorter baselines, greater number of baselines, and single-dish capabilities, however, give it a significant edge for observing extended emission, making wide-field maps (mosaics), or attaining high image fidelity, as required by the described science cases.

We present high-angular-resolution observations of the massive star forming core DR21(OH) at 880 mum using the Submillimeter Array (SMA). The dense core exhibits an overall velocity gradient in a Keplerian-like pattern, which breaks at the center of the core where SMA 6 and SMA 7 are located. The dust polarization shows a complex magnetic field, compatible with a toroidal configuration. This is in contrast with the large, parsec-scale filament that surrounds the core, where there is a smooth magnetic field. The total magnetic field strengths in the filament and in the core are 0.9 and 2.1 mG, respectively. We found evidence of magnetic field diffusion at the core scales, far beyond the expected value for ambipolar diffusion. It is possible that the diffusion arises from fast magnetic reconnection in the presence of turbulence. The dynamics of the DR 21(OH) core appear to be controlled energetically in equal parts by the magnetic field, magneto-hydrodynamic (MHD) turbulence and the angular momentum. The effect of the angular momentum (this is a fast rotating core) is probably causing the observed toroidal field configuration. Yet, gravitation overwhelms all the forces, making this a clear supercritical core with a mass-to-flux ratio of ~6 times the critical value. However, simulations show that this is not enough for the high level of fragmentation observed at 1000 AU scales. Thus, rotation and outflow feedback is probably the main cause of the observed fragmentation.

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.

Commonplace at every evolutionary stage, Multiple Protostellar Systems (MPSs) are thought to be formed through fragmentation, but it is unclear when and how. The youngest MPSs, which have not yet undergone much evolution, provide important constraints to this question. It is then of interest to disentangle early stage MPSs. In this letter we present the results of our work on VLA1623 using our observations and archival data from the Submillimeter Array (SMA). Our continuum and line observations trace VLA1623's components, outflow and envelope, revealing unexpected characteristics. We construct the SED for each component using the results of our work and data from literature, as well as derive physical parameters from continuum and perform a simple kinematical analysis of the circumstellar material. Our results show VLA1623 to be a triple non-coeval system composed of VLA1623A, B & W, with each source driving its own outflow and unevenly distributed circumstellar material. From the SED, physical parameters and IR emission we conclude that VLA1623A & W are Class 0 and I protostars, respectively, and together drive the bulk of the observed outflow. Furthermore, we find two surprising results, first the presence of a rotating disk-like structure about VLA1623A with indications of pure Keplerian rotation, which, if real, would make it one of the first evidence of Keplerian disk structures around Class 0 protostars. Second, we find VLA1623B to be a bonafide extremely young protostellar object between the starless core and Class 0 stages.

We present a set of compelling science cases for the ALMA Band 1 receiver suite. For these cases, we assume in tandem the updated nominal Band 1 frequency range of 35-50 GHz with a likely extension up to 52 GHz; together these frequencies optimize the Band 1 science return. The scope of the science cases ranges from nearby stars to the re-ionization edge of the Universe. Two cases provide additional leverage on the present ALMA Level One Science Goals and are seen as particularly powerful motivations for building the Band 1 Receiver suite: (1) detailing the evolution of grains in protoplanetary disks, as a complement to the gas kinematics, requires continuum observations out to ~35 GHz (~9mm); and (2) detecting CO 3-2 line emission from galaxies like the Milky Way during the epoch of re-ionization, i.e., 6 < z < 10, also requires Band 1 receiver coverage. The range of Band 1 science is wide, however, and includes studies of very small dust grains in the ISM, pulsar wind nebulae, radio supernovae, X-ray binaries, the Galactic Center (i.e., Sgr A*), dense cloud cores, complex carbon-chain molecules, masers, magnetic fields in the dense ISM, jets and outflows from young stars, distant galaxies, and galaxy clusters (i.e., the Sunyaev-Zel'dovich Effect). A comparison of ALMA and the Jansky VLA (JVLA) at the same frequencies of Band 1 finds similar sensitivity performance at 40-50 GHz, with a slight edge for ALMA at higher frequencies (e.g., within a factor of 2 for continuum observations). With its larger number of instantaneous baselines, however, ALMA Band 1data will have greater fidelity than those from the JVLA at similar frequencies.

The J-Ks color distribution (CD) with a bin size of 0.05 magnitude for the entire Milky Way has been carried out by using the Two Micron All Sky Survey Point Source Catalog (2MASS PSC). The CDs are bimodal, which has a red peak at 0.8 < J-Ks < 0.85 and a blue peak at 0.3 < J-Ks < 0.4. The colors of the red peak are more or less the same for the whole sky, but that of the blue peak depend on Galactic latitude, (J-Ks ~ 0.35 at low Galactic latitudes and 0.35 < J-Ks < 0.4 for other sky areas). The blue peak dominates the bimodal CDs at low Galactic latitudes and becomes comparable with the red peak in other sky regions. In order to explain the bimodal distribution and the global trend shown by the all sky 2MASS CDs, we assemble an empirical HR diagram, which is composed by observational-based near infrared HR diagrams and color magnitude diagrams, and incorporate a Milky Way model. In the empirical HR diagram, the main sequence stars turnoff the thin disk is relatively bluer, (J-Ks)0 = 0.31, when we compare with the thick disk which is (J-Ks)0 = 0.39. The age of the thin/thick disk is roughly estimated to be around 4-5/8-9 Gyr according to the color-age relation of the main sequence turnoff. In general, the 2MASS CDs can be treated as a tool to census the age of stellar population of the Milky Way in a statistical manner and to our knowledge this is a first attempt to measure the age.

The Infrared Array Camera (IRAC) on the Spitzer Space Telescope has revealed that a number of high-mass protostars are associated with extended mid-infrared emission, particularly prominent at 4.5-micron. These are called "Green Fuzzy" emission or "Extended Green Objects". We present color analysis of this emission toward six nearby (d=2-3 kpc) well-studied high-mass protostars and three candidate high-mass protostars identified with the Spitzer GLIMPSE survey. In our color-color diagrams most of the sources show a positive correlation between the [3.6]-[4.5] and [3.5]-[5.8] colors along the extinction vector in all or part of the region. We compare the colors with those of scattered continuum associated with the low-mass protostar L 1527, modeled scattered continuum in cavities, shocked emission associated with low-mass protostars, modeled H2 emission for thermal and fluorescent cases, and modeled PAH emission. Of the emission mechanisms discussed above, scattered continuum provides the simplest explanation for the observed linear correlation. In this case, the color variation within each object is attributed to different foreground extinctions at different positions. Alternative possible emission mechanisms to explain this correlation may be a combination of thermal and fluorescent H2 emission in shocks, and a combination of scattered continuum and thermal H2 emission, but detailed models or spectroscopic follow-up are required to further investigate this possibility. Our color-color diagrams also show possible contributions from PAHs in two objects. However, none of our sample show clear evidence for PAH emission directly associated with the high-mass protostars, several of which should be associated with ionizing radiation. This suggests that those protostars are heavily embedded even at mid-infrared wavelengths.

We used the Submillimeter Array (SMA) to observe the thermal polarized dust emission from the protostellar source NGC 2024 FIR 5. The polarized emission outlines a partial hourglass morphology for the plane-of-sky component of the core magnetic field. Our data are consistent with previous BIMA maps, and the overall magnetic field geometries obtained with both instruments are similar. We resolve the main core into two components, FIR 5A and FIR 5B. A possible explanation for the asymmetrical field lies in depolarization effects due to the lack of internal heating from FIR 5B source, which may be in a prestellar evolutionary state. The field strength was estimated to be 2.2 mG, in agreement with previous BIMA data. We discuss the influence of a nearby H{\sc ii} region over the field lines at scales of $\sim 0.01$ pc. Although the hot component is probably compressing the molecular gas where the dust core is embedded, it is unlikely that the radiation pressure exceeds the magnetic tension. Finally, a complex outflow morphology is observed in CO (3 $\rightarrow$ 2) maps. Unlike previous maps, several features associated with dust condensations other than FIR 5 are detected.

Though flux freezing is a good approximation frequently assumed for molecular clouds, ambipolar diffusion (AD) is inevitable at certain scales. The scale at which AD sets in can be a crucial parameter for turbulence and the star formation process. However, both observation and simulation of AD are very challenging and our knowledge of it is very limited. We proposed earlier (Li and Houde 2008) that the difference between ion and neutral velocity spectra is a signature of turbulent AD and can be used to estimate the AD scales and magnetic field strength. Here we present observational evidence showing that this difference between the velocity dispersions from coexistent ions and neutrals is indeed correlated with magnetic field strength.

We present a near-infrared extinction study of nine dense cores at evolutionary stages between starless to Class I. Our results show that the density structure of all but one observed cores can be modeled with a single power law rho \propto r^p between ~ 0.2R-R of the cores. The starless cores in our sample show two different types of density structures, one follows p ~ -1.0 and the other follows p ~ -2.5, while the protostellar cores all have p ~ -2.5. The similarity between the prestellar cores with p ~ -2.5 and protostellar cores implies that those prestellar cores could be evolving towards the protostellar stage. The slope of p ~ -2.5 is steeper than that of an singular isothermal sphere, which may be interpreted with the evolutionary model of cores with finite mass.

We present polarization observations of the dust continuum emission from IRAS 16293 which is a Class 0 protostar and is known to have at least two components, source A and B. These measurements were conducted by the Submillimeter Array (SMA) at a frequency of ~341 GHz and with high angular resolution (~2-3 arcseconds). We find that the large scale global direction of the field, which is perpendicular to the observed polarization, appears to be along the dust ridge where the emission peaks. On smaller scales we find that the field structure is significantly different for the two components of the binary. The first component, source A, shows a magnetic field structure which is "hourglass" shaped as predicted from theoretical models of low mass star formation in the presence of strong magnetic fields. However, the other component, source B, shows a relatively ordered magnetic field with no evidence of any deformation. There is an observed decrease in polarization towards the center. Our calculations show that in IRAS 16293 the magnetic energy is stronger than the turbulent energy but is approximately similar to the centrifugal energy. Our results provide additional evidence to show that the two protostars appear to be in different stages during their evolution.

We report our observational results of 870 $\mu$m continuum emission and its linear polarization in the massive star formation site W51 e2/e8. Inferred from the linear polarization maps, the magnetic field in the plane of sky (B$_{\bot}$) is traced with an angular resolution of 0$\farcs$7 with the Submillimeter Array (SMA). Whereas previous BIMA observations with an angular resolution of 3$\arcsec$ (0.1 pc) showed a uniform B field, our revealed B$_{\bot}$ morphology is hourglass-like in the collapsing core near the Ultracompact H II region e2 and also possibly in e8. The decrease in polarization near the continuum peak seen at lower angular resolution is apparently due to the more complex structures at smaller scales. In e2, the pinched direction of the hourglass-like B field morphology is parallel to the plane of the ionized accretion flow traced by H53$\alpha$, suggesting that the massive stars are formed via processes similar to the low mass stars, i.e. accretion through a disk, except that the mass involved is much larger. Furthermore, our finding that the resolved collapsing cores in e2 and e8 lie within one subcritical 0.5 pc envelope supports the scenario of \textit{magnetic fragmentation} via ambipolar diffusion. We therefore suggest that magnetic fields control the dynamical evolution of the envelope and cores in W51 e2 and e8.

We report high angular resolution (3") Submillimeter Array (SMA) observations of the molecular cloud associated with the Ultra-Compact HII region G5.89-0.39. Imaged dust continuum emission at 870 micron reveals significant linear polarization. The position angles (PAs) of the polarization vary enormously but smoothly in a region of 2x10^4 AU. Based on the distribution of the PAs and the associated structures, the polarized emission can be separated roughly into two components. The component "x" is associated with a well defined dust ridge at 870 micron, and is likely tracing a compressed B field. The component "o" is located at the periphery of the dust ridge and is probably from the original B field associated with a pre-existing extended structure. The global B field morphology in G5.89, as inferred from the PAs, is clearly disturbed by the expansion of the HII region and the molecular outflows. Using the Chandrasekhar-Fermi method, we estimate from the smoothness of the field structures that the B field strength in the plane of sky can be no more than 2-3 mG. We then compare the energy densities in the radiation, the B field, and the mechanical motions as deduced from the C^17O 3-2 line emission. We conclude that the B field structures are already overwhelmed and dominated by the radiation, outflows, and turbulence from the newly formed massive stars.

We compute the mid-infrared extinction law from 3.6-24 microns in three molecular clouds: Ophiuchus, Perseus, and Serpens, by combining data from the "Cores to Disks" Spitzer Legacy Science program with deep JHKs imaging. Using a new technique, we are able to calculate the line-of-sight extinction law towards each background star in our fields. With these line-of-sight measurements, we create, for the first time, maps of the chi-squared deviation of the data from two extinction law models. Because our chi-squared maps have the same spatial resolution as our extinction maps, we can directly observe the changing extinction law as a function of the total column density. In the Spitzer IRAC bands, 3.6-8 microns, we see evidence for grain growth. Below $A_{K_s} = 0.5$, our extinction law is well-fit by the Weingartner & Draine (2001) $R_V = 3.1$ diffuse interstellar medium dust model. As the extinction increases, our law gradually flattens, and for $A_{K_s} >= 1$, the data are more consistent with the Weingartner & Draine $R_V = 5.5$ model that uses larger maximum dust grain sizes. At 24 microns, our extinction law is 2-4 times higher than the values predicted by theoretical dust models, but is more consistent with the observational results of Flaherty et al. (2007). Lastly, from our chi-squared maps we identify a region in Perseus where the IRAC extinction law is anomalously high considering its column density. A steeper near-infrared extinction law than the one we have assumed may partially explain the IRAC extinction law in this region.

We present maps of 14.4 deg^2 of the Ophiuchus dark clouds observed by the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). These high quality maps depict both numerous point sources as well as extended dust emission within the star-forming and non-star-forming portions of these clouds. Using PSF-fitting photometry, we detect 5779 sources at 24 um and 81 sources at 70 um at the 10 sigma level of significance. Three hundred twenty-three candidate young stellar objects (YSOs) were identified according to their positions on the MIPS/2MASS K versus K$-$[24] color-magnitude diagrams as compared to 24 um detections in the SWIRE extragalactic survey. We find that more than half of the YSO candidates, and almost all the ones with protostellar Class I spectral energy distributions, are confined to the known cluster and aggregates.

We present maps of 7.78 square degrees of the Lupus molecular cloud complex at 24, 70, and $160\:\mu$m. They were made with the Spitzer Space Telescope's Multiband Imaging Photometer for Spitzer (MIPS) instrument as part of the Spitzer Legacy Program, ``From Molecular Cores to Planet-Forming Disks'' (c2d). The maps cover three separate regions in Lupus, denoted I, III, and IV. We discuss the c2d pipeline and how our data processing differs from it. We compare source counts in the three regions with two other data sets and predicted star counts from the Wainscoat model. This comparison shows the contribution from background galaxies in Lupus I. We also create two color magnitude diagrams using the 2MASS and MIPS data. From these results, we can identify background galaxies and distinguish them from probable young stellar objects. The sources in our catalogs are classified based on their spectral energy distribution (SED) from 2MASS and Spitzer wavelengths to create a sample of young stellar object candidates. From 2MASS data, we create extinction maps for each region and note a strong corresponence between the extinction and the $160\:\mu$m emission. The masses we derived in each Lupus cloud from our extinction maps are compared to masses estimated from $^{13}$CO and C$^{18}$O and found to be similar to our extinction masses in some regions, but significantly different in others. Finally, based on our color-magnitude diagrams, we selected 12 of our reddest candidate young stellar objects for individual discussion. Five of the 12 appear to be newly-discovered YSOs.

One of the central goals of the Spitzer Legacy Project ``From Molecular Cores to Planet-forming Disks'' (c2d) is to determine the frequency of remnant circumstellar disks around weak-line T Tauri stars (wTTs) and to study the properties and evolutionary status of these disks. Here we present a census of disks for a sample of over 230 spectroscopically identified wTTs located in the c2d IRAC (3.6, 4.5, 4.8, and 8.0 um) and MIPS (24 um) maps of the Ophiuchus, Lupus, and Perseus Molecular Clouds. We find that ~20% of the wTTs in a magnitude limited subsample have noticeable IR-excesses at IRAC wavelengths indicating the presence of a circumstellar disk. The disk frequencies we find in these 3 regions are ~3-6 times larger than that recently found for a sample of 83 relatively isolated wTTs located, for the most part, outside the highest extinction regions covered by the c2d IRAC and MIPS maps. The disk fractions we find are more consistent with those obtained in recent Spitzer studies of wTTs in young clusters such as IC 348 and Tr 37. From their location in the H-R diagram, we find that, in our sample, the wTTs with excesses are among the younger part of the age distribution. Still, up to ~50% of the apparently youngest stars in the sample show no evidence of IR excess, suggesting that the circumstellar disks of a sizable fraction of pre-main-sequence stars dissipate in a timescale of ~1 Myr. We also find that none of the stars in our sample apparently older than ~10 Myrs have detectable circumstellar disks at wavelengths < 24 um. Also, we find that the wTTs disks in our sample exhibit a wide range of properties (SED morphology, inner radius, L_DISK/L*, etc) which bridge the gaps observed between the cTTs and the debris disk regimes.

We present maps of 1.5 square degrees of the Serpens dark cloud at 24, 70, and 160\micron observed with the Spitzer Space Telescope MIPS Camera. More than 2400 compact sources have been extracted at 24um, nearly 100 at 70um, and 4 at 160um. We estimate completeness limits for our 24um survey from Monte Carlo tests with artificial sources inserted into the Spitzer maps. We compare source counts, colors, and magnitudes in the Serpens cloud to two reference data sets, a 0.50 deg^2 set on a low-extinction region near the dark cloud, and a 5.3 deg^2 subset of the SWIRE ELAIS N1 data that was processed through our pipeline. These results show that there is an easily identifiable population of young stellar object candidates in the Serpens Cloud that is not present in either of the reference data sets. We also show a comparison of visual extinction and cool dust emission illustrating a close correlation between the two, and find that the most embedded YSO candidates are located in the areas of highest visual extinction.

We present observations of 10.6 square degrees of the Perseus molecular cloud at 24, 70, and 160 microns with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). The image mosaics show prominent, complex extended emission dominated by illuminating B stars on the East side of the cloud, and by cold filaments of 160 micron emission on the West side. Of 3950 point sources identified at 24 microns, 1141 have 2MASS counterparts. A quarter of these populate regions of the Ks vs. Ks-[24] diagram that are distinct from stellar photospheres and background galaxies, and thus are likely to be cloud members with infrared excess. Nearly half (46%) of these 24 micron excess sources are distributed outside the IC 348 and NGC 1333 clusters. NGC 1333 shows the highest fraction of stars with flat or rising spectral energy distributions (28%), while Class II SEDs are most common in IC 348. These results are consistent with previous relative age determinations for the two clusters. The intercluster region contains several tightly clumped (r~0.1 pc) young stellar aggregates whose members exhibit a wide variety of infrared spectral energy distributions characteristic of different circumstellar environments. One possible explanation is a significant age spread among the aggregate members, such that some have had time to evolve more than others. Alternatively, if the aggregate members all formed at roughly the same time, then remarkably rapid circumstellar evolution would be required to account for the association of Class I and Class III sources at ages <~1 Myr. We highlight important results for several other objects as well (full abstract in the paper).

Infrared images of the dark cloud core B59 were obtained with the Spitzer Space Telescope as part of the "Cores to Disks" Legacy Science project. Photometry from 3.6-70 microns indicates at least 20 candidate low-mass young stars near the core, more than doubling the previously known population. Out of this group, 13 are located within about 0.1 pc in projection of the molecular gas peak, where a new embedded source is detected. Spectral energy distributions span the range from small excesses above photospheric levels to rising in the mid-infrared. One other embedded object, probably associated with the millimeter source B59-MMS1, with a bolometric luminosity L(bol) roughly 2 L(sun), has extended structure at 3.6 and 4.5 microns, possibly tracing the edges of an outflow cavity. The measured extinction through the central part of the core is A(V) greater than of order 45 mag. The B59 core is producing young stars with a high efficiency.

We present a comprehensive study of a low-mass star-forming region,L1251B, at wavelengths from the near-infrared to the millimeter. L1251B, where only one protostar, IRAS 22376+7455, was known previously, is confirmed to be a small group of protostars based on observations with the Spitzer Space Telescope. The most luminous source of L1251B is located 5" north of the IRAS position. A near-infrared bipolar nebula, which is not associated with the brightest object and is located at the southeast corner of L1251B, has been detected in the IRAC bands. OVRO and SMA interferometric observations indicate that the brightest source and the bipolar nebula source in the IRAC bands are deeply embedded disk sources.Submillimeter continuum observations with single-dish telescopes and the SMA interferometric observations suggest two possible prestellar objects with very high column densities. Outside of the small group, many young stellar object candidates have been detected over a larger region of 12' x 12'. Extended emission to the east of L1251B has been detected at 850 micron; this "east core" may be a site for future star formation since no point source has been detected with IRAC or MIPS. This region is therefore a possible example of low-mass cluster formation, where a small group of pre- and protostellar objects (L1251B) is currently forming, alongside a large starless core (the east core).

We present observations of 3.86 sq. deg. of the Perseus molecular cloud complex with the Spitzer Space Telescope Infrared Array Camera (IRAC). The maps show strong extended emission arising from shocked H2 in outflows in the region and from polycyclic aromatic hydrocarbon features. More than 120,000 sources are extracted toward the cloud. Based on their IRAC colors and comparison to off-cloud and extragalactic fields, we identify 400 candidate young stellar objects. About two thirds of these are associated with the young clusters IC348 and NGC1333, while the last third is distributed over the remaining cloud. We classify the young stellar objects according to the traditional scheme based on the slope of their spectral energy distributions. Significant differences are found for the numbers of embedded Class I objects relative to more evolved Class II objects in IC348, NGC1333 and the remaining cloud with the embedded Class I and "flat spectrum" YSOs constituting 14%, 36% and 47% of the total number of YSOs identified in each of these regions. The high number of Class I objects in the extended cloud (61% of the Class I objects in the entire cloud) suggests that a significant fraction of the current star formation is occuring outside the two main clusters. Finally we discuss a number of outflows and identify their driving sources, including the known deeply embedded Class 0 sources outside the two major clusters. The Class 0 objects are found to be detected by Spitzer and have very red [3.6]-[4.5] colors but do not show similarly red [5.8]-[8.0] colors. The Class 0 objects are easily identifiable in color-color diagrams plotting these two colors but are problematic to extract automatically due to the extended emission from shocked gas or scattered light in cavities related to the associated outflows.

Using the Spitzer Space Telescope, we have observed 90 weak-line and classical T Tauri stars in the vicinity of the Ophiuchus, Lupus, Chamaeleon, and Taurus star-forming regions as part of the Cores to Disks (c2d) Spitzer Legacy project. In addition to the Spitzer data, we have obtained contemporaneous optical photometry to assist in constructing spectral energy distributions. These objects were specifically chosen as solar-type young stars with low levels of H alpha emission, strong X-ray emission, and lithium absorption i.e. weak-line T Tauri stars, most of which were undetected in the mid-to-far IR by the IRAS survey. Weak-line T Tauri stars are potentially extremely important objects in determining the timescale over which disk evolution may take place. Our objective is to determine whether these young stars are diskless or have remnant disks which are below the detection threshold of previous infrared missions. We find that only 5/83 weak-line T Tauri stars have detectable excess emission between 3.6 and 70 micron which would indicate the presence of dust from the inner few tenths of an AU out to the planet-forming regions a few tens of AU from the star. Of these sources, two have small excesses at 24 microns consistent with optically thin disks; the others have optically thick disks already detected by previous IR surveys. All of the seven classical T Tauri stars show excess emission at 24 and 70 micron, although their properties vary at the shorter wavelengths. Our initial results show that disks are rare among young stars selected for their weak H alpha emission.

We present maps of over 1.5 square degrees in Chamaeleon (Cha) II at 24, 70, and 160 micron observed with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS) and a 1.2 square degree millimeter map from SIMBA on the Swedish-ESO Submillimetre Telescope (SEST). The c2d Spitzer Legacy Team's data reduction pipeline is described in detail. Over 1500 24 micron sources and 41 70 micron sources were detected by MIPS with fluxes greater than 10-sigma. More than 40 potential YSOs are identified with a MIPS and 2MASS color-color diagram and by their spectral indices, including two previously unknown sources with 24 micron excesses. Our new SIMBA millimeter map of Cha II shows that only a small fraction of the gas is in compact structures with high column densities. The extended emission seen by MIPS is compared with previous CO observations. Some selected interesting sources, including two detected at 1 mm, associated with Cha II are discussed in detail and their SEDs presented. The classification of these sources using MIPS data is found to be consistent with previous studies.

We present observations of L1014, a dense core in the Cygnus region previously thought to be starless, but data from the Spitzer Space Telescope shows the presence of an embedded source. We propose a model for this source that includes a cold core, heated by the interstellar radiation field, and a low-luminosity internal source. The low luminosity of the internal source suggests a substellar object. If L1014 is representative, other "starless" cores may turn out to harbor central sources.

It has been suggested that under average interstellar field strengths the cyclotron interaction between ions and magnetic fields is strong enough to narrow the linewidth and suppress the line wings in the ion spectra. We present evidence for the cyclotron interaction effect at arcsec scale on the velocity dispersions in the spectra of ion/neutral molecular species in DR21(OH) observed with the OVRO-MMA. Using a spatial resolution $\sim$ 3 times higher than previous CSO observations by Houde et al. (2002), we show that H$^{13}$CO$^+$ and H$^{13}$CN are coexistent at the scale of our observations (6\arcsec). In the eastern parts of the DR21(OH) core where the dynamics is simple, the ion linewidths are indeed narrower than the neutral linewidths with an average ion-to-neutral linewidth ratio of 0.82$\pm$0.04. We use our results, along with the existing Zeeman and dust/CO polarization data on small scales, to derive the 3-D magnetic field structure. We obtain a field strength of $0.44\pm$0.12 mG with inclination of 36\deg to the line of sight, directed toward the observer, and a position angle of $-$75\deg in the plane of the sky. With the full magnetic field strength derived here, we are able to conclude that the MM1 core of DR21(OH) is magnetically supercritical; although turbulence provides the dominant support.

We present the polarization detections in DR21(OH) from both the thermal dust emission at 1.3 mm and the CO J=2-1 line obtained with the Berkeley-Illinois-Maryland Association (BIMA) array. Our results are consistent with the prediction of the Goldreich-Kylafis effect that the CO polarization is either parallel or perpendicular to the magnetic field direction. The detection of the polarized CO emission is over a more extended region than the dust polarization, while the dust polarization provides an aide in resolving the ambiguity of the CO polarization. The combined results suggest that the magnetic field direction in DR21(OH) is parallel to the CO polarization and therefore parallel to the major axis of DR21(OH). The strong correlation between the CO and dust polarization suggests that magnetic fields are remarkably uniform throughout the envelope and the cores. The dispersion in polarization position angles implies a magnetic field strength in the plane of the sky of about 1 mG, compared with about 0.5 mG inferred for the line-of-sight field from previous CN Zeeman observations. Our CO data also show that both MM1 and MM2 power high-velocity outflows with $v\gtrsim$25 km/s relative to the systematic velocity.

We have investigated the physical and chemical status of the pre-protostellar core B68. A previous extinction study suggested that the density profile of B68 is remarkably consistent with a Bonnor-Ebert sphere with 2.1 $M_{\odot}$ at 16 K. We mapped B68 in C3H2, CCS, and NH3 with the Deep Space Network (DSN) 70m telescope at Goldstone. Our results show that the B68 chemical structure is consistent with that seen in other pre-protostellar cores (L1498, L1544) and is explained by time dependent chemical models that include depletion. We measured a kinetic temperature of 11K with NH3 (1,1) and (2,2) spectra. We also show that B68 is thermally dominated with little contribution from turbulence support (< 10%). We consider a modified Bonnor-Ebert sphere to include effects of turbulence and magnetic fields and use it to constrain the uncertainties in its distance determination. We conclude that the distance to B68 is $\sim$ 95pc with a corresponding mass of $\sim$ 1.0 $M_{\odot}$. If some magnetic field is present it can be further away (beyond $\sim$ 100pc) and still satisfy the density structure of a Bonnor-Ebert sphere. Our observations do not preclude any instability such as the onset of collapse, or slow contraction, occurring in the center of the core, which cannot be resolved with our beam size (45").

NGC 5471B has been suggested to contain a hypernova remnant because of its extraordinarily bright X-ray emission. To assess its true nature, we have obtained high-resolution images in continuum bands and nebular lines with the Hubble Space Telescope, and high-dispersion long-slit spectra with the Kitt Peak National Observatory 4-m echelle spectrograph. The images reveal three supernova remnant (SNR) candidates in the giant HII region NGC 5471, with the brightest one being the 77x60 pc shell in NGC 5471B. The Ha velocity profile of NGC 5471B can be decomposed into a narrow component (FWHM = 41 km/s) from the background HII region and a broad component (FWHM = 148 km/s) from the SNR shell. Using the brightness ratio of the broad to narrow components and the Ha flux measured from the WFPC2 Ha image, we derive an Ha luminosity of (1.4+-0.1)x10^39 ergs/s for the SNR shell. The [SII]6716,6731 doublet ratio of the broad velocity component is used to derive an electron density of ~700 cm^-3 in the SNR shell. The mass of the SNR shell is thus 4600+-500 Mo. With a \~330 km/s expansion velocity implied by the extreme velocity extent of the broad component, the kinetic energy of the SNR shell is determined to be 5x10^51 ergs. This requires an explosion energy greater than 10^52 ergs, which can be provided by one hypernova or multiple supernovae. Comparing to SNRs in nearby active star formation regions, the SNR shell in NGC 5471B appears truly unique and energetic. We conclude that the optical observations support the existence of a hypernova remnant in NGC 5471B.

We present the first interferometric polarization maps of the NGC2024 FIR5 molecular core obtained with the BIMA array at approximately 2 arcsec resolution. We measure an average position angle of -60+-6 degrees in the main core of FIR5 and 54+-9 degrees in the eastern wing of FIR5. The morphology of the polarization angles in the main core of FIR5 suggests that the field lines are parabolic with a symmetry axis approximately parallel to the major axis of the putative disk in FIR5, which is consistent with the theoretical scenario that the gravitational collapse pulled the field lines into an hour-glass shape. The polarization percentage decreases toward regions with high intensity and close to the center of the core, suggesting that the dust alignment efficiency may decrease at high density. The plane-of-sky field strength can be estimated with the modified Chandrasekhar-Fermi formula, and the small dispersion of the polarization angles in FIR5 suggests that the magnetic field is strong ($\gtrsim$ 2mG) and perhaps dominates the turbulent motions in the core.

We present the first interferometric polarization map of the W51 e1/e2 molecular cores obtained with the BIMA array at 1.3 mm wavelength with approximately 3 arcsecond resolution. The polarization angle varies smoothly across the double cores with an average position angle of 23+-5 degrees for W51 e1 and 15+-7 degrees for W51 e2. The inferred magnetic field direction is parallel to the minor axis of the double cores, which is consistent with the theoretical picture that clouds collapse along the field lines. However, the magnetic field may not determine the axis of angular momentum of these two cores as the field directions of the two cores significantly differ with the previously measured directions of rotational axes. The polarization percentage decreases toward regions with high intensity, suggesting that the dust alignment efficiency decreases toward high density regions. The field directions are highly ordered, and the small dispersion of the polarization angles implies that magnetic fields are strong ($\gtrsim$ 1 mG) and perhaps dominate turbulence in W51 e1/e2.

The SNR candidate MF83 in M101 is coincident with a very luminous X-ray source. Based on the high X-ray luminosity, it has been suggested that MF83 is a "hypernova remnant" requiring an explosion energy about two orders of magnitude higher than normal supernovae. We have analyzed high-quality ground-based and HST observations of MF83, and find that MF83 is a star formation region, consisting of a large ionized gas shell and four HII regions along its periphery. Continuum images show OB associations in these HII regions and within the large shell. The shell has an expansion velocity of $\sim$50 km/s and a diameter of $\sim$270 pc. The optical properties of this shell in MF83 are similar to those of X-ray-bright superbubbles in the Large Magellanic Cloud. If the X-ray emission is indeed diffuse, the implied thermal energy in MF83 is high, a few $\times10^{52}$ ergs. This amount of thermal energy requires a large number of concentrated supernova explosions or one powerful explosion. Future X-ray observations with a high angular resolution are needed to resolve the diffuse emission and point sources in MF83, in order to determine more accurately the thermal energy in the shell interior and its required explosion energy.

We have conducted a survey of the CCS $J_N = 3_2-2_1$ line toward 11 dark clouds and star-forming regions at 30 arcsec spatial resolution and 0.054 km/s velocity resolution. CCS was only detected in quiescent clouds, not in active star-forming regions. The CCS distribution shows remarkable clumpy structure, and 25 clumps are identified in 7 clouds. Seven clumps with extremely narrow nonthermal linewidths < 0.1 km/s are among the most quiescent clumps ever found. The CCS clumps tend to exist around the higher density regions traced by NH_3 emission or submillimeter continuum sources, and the distribution is not spherically symmetric. Variation of the CCS abundance was suggested as an indicator of the evolutionary status of star formation. However, we can only find a weak correlation between N(CCS) and $n_{H_2,vir}$. The velocity distributions of CCS clouds reveal that a systematic velocity pattern generally exists. The most striking feature in our data is a ring structure in the position-velocity diagram of L1544 with an well-resolved inner hole of 0.04 pc x 0.13 km/s and an outer boundary of 0.16 pc x 0.55 km/s. This position-velocity structure clearly indicates an edge-on disk or ring geometry, and it can be interpreted as a collapsing disk with an infall velocity $\gtrsim$ 0.1 km/s and a rotational velocity less than our velocity resolution. Nonthermal linewidth distribution is generally coherent in CCS clouds, which could be evidence for the termination of Larson's Law at small scales, $\sim$ 0.1 pc.

Some extragalactic SNRs are more than two orders of magnitude more luminous than the young Galactic SNR Cas A. These SNRs are called super-luminous or ultra-luminous SNRs. Their high luminosities can be caused by chance superpositions of multiple objects, interactions with a very dense environment, or unusually powerful supernova explosions. Four super-luminous SNRs are known: one in NGC 4449, one in NGC 6946, and two in M101. The two remnants in M101, NGC 5471B and MF83, are recently suggested to be "hypernova remnants" possibly connected to the GRBs. We have obtained new or archival HST WFPC2 images and new high-dispersion echelle spectra of these super-luminous SNRs, in order to examine their stellar and interstellar environments and to analyze their energetics. We discuss the physical nature of these four SNRs, with a special emphasis on the two ``hypernova remnants'' in M101.