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.

We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

We present 0.25 arcsec resolution (35 au) ALMA 1.3 mm dust polarization observations for 37 young stellar objects (YSOs) in the Ophiuchus cloud. These data encompass all the embedded protostars in the cloud and several Flat and Class II objects to produce the largest, homogeneous study of dust polarization on disk scales to date. The goal of this study is to study dust polarization down to disk scales. We find that 14/37 (38%) of the YSOs are detected in polarization. Nine of these sources have uniform polarization angles and four sources have azimuthal polarization. The sources with uniform polarization tend to have steeper inclinations (> 60 degree) than those with azimuthal polarization (< 60 degree). The majority (9/14) of the detected sources have polarization morphologies and disk properties consistent with dust self-scattering in optically thick disks. The remaining sources may be instead tracing magnetic fields. Their inferred field directions from rotating the polarization vectors by 90 degree are mainly poloidal or hourglass shaped. We find no evidence of a strong toroidal field component toward any of our disks. For the 23 YSOs that are undetected in polarization, roughly half of them have 3-sigma upper limits of < 2%. These sources also tend to have inclinations < 60 degree and they are generally compact. Since lower inclination sources tend to have azimuthal polarization, these YSOs may be undetected in polarization due to unresolved polarization structure within our beam. We propose that disks with inclinations > 60 degree are the best candidates for future polarization studies of dust self-scattering as these systems will generally show uniform polarization vectors that do not require very high resolution to resolve. We release the continuum and polarization images for all the sources with this publication. Data from the entire survey can be obtained from Dataverse.

L1521F is found to be forming multiple cores and it is cited as an example of the densest core with an embedded VeLLO in a highly dynamical environment. We present the core-scale magnetic fields (B-fields) in the near vicinity of the VeLLO L1521F-IRS using submm polarization measurements at 850$~\mu$m using JCMT POL-2. This is the first attempt to use high-sensitivity observations to map the sub-parsec scale B-fields in a core with a VeLLO. The B-fields are ordered and very well connected to the parsec-scale field geometry seen in our earlier optical polarization observations and the large-scale structure seen in Planck dust polarization. The core scale B-field strength estimated using Davis-Chandrasekhar-Fermi relation is $\rm 330\pm100~\mu$G which is more than ten times of the value we obtained in the envelope (envelope in this paper is "core envelope"). This indicates that B-fields are getting stronger on smaller scales. The magnetic energies are found to be 1 to 2 orders of magnitude higher than non-thermal kinetic energies in the envelope and core. This suggests that magnetic fields are more important than turbulence in the energy budget of L1521F. The mass-to-flux ratio of 2.3$\pm$0.7 suggests that the core is magnetically-supercritical. The degree of polarization is steadily decreasing towards the denser part of the core with a power law slope of -0.86.

This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Galactic Plane survey, a larger Gould Belt survey, and a follow-up survey of all Planck compact objects visible from the northern hemisphere. These surveys will provide a complete census of molecular clouds, filaments, and dense cores across the Galaxy, vital for studying star formation in various environments.

In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt.

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.

We present polarimetric observations of four Class II protoplanetary disks (DG Tau, Haro 6-13, RY Tau, and MWC 480) taken with the Atacama Large Millimeter/submillimeter Array (ALMA) at 3 mm. The polarization morphologies observed fall into two distinct categories: azimuthal polarization (DG Tau and Haro 6-13) and polarization parallel to the disk minor axis (RY Tau and MWC 480). The mechanisms responsible for disk polarization at millimeter and submillimeter wavelengths are currently under debate. In this Letter, we investigate two mechanisms capable of producing polarized dust emission in disks: self-scattering and grain alignment to the radiation anisotropy. The polarization morphologies of DG Tau and Haro 6-13 are broadly consistent with that expected from radiation alignment (though radiative alignment still does not account for all of the features seen in these disks), while RY Tau and MWC 480 are more consistent with self-scattering. Such a variation in the polarized morphology may provide evidence of dust grain size differences between the sources.

We have carried out polarimetric observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostellar system L1448 IRS 2, which is a proto-binary embedded within a flattened, rotating structure, and for which a hint of a central disk has been suggested, but whose magnetic fields are aligned with the bipolar outflow on the cloud core scale. Our high sensitivity and high resolution ($\sim 100$ au) observations show a clear hourglass magnetic field morphology centered on the protostellar system, but the central pattern is consistent with a toroidal field indicative of a circumstellar disk, although other interpretations are also possible, including field lines dragged by an equatorial accretion flow into a configuration parallel to the midplane. If a relatively large disk does exist, it would suggest that the magnetic braking catastrophe is averted in this system, not through a large misalignment between the magnetic and rotation axes, but rather through some other mechanisms, such as non-ideal magneto-hydrodynamic effects and/or turbulence. We have also found a relationship of decreasing polarization fractions with intensities and the various slopes of this relationship can be understood as multiple polarization mechanisms and/or depolarization from a changing field morphology. In addition, we found a prominent clumpy depolarization strip crossing the center perpendicular to the bipolar outflow. Moreover, a rough estimate of the magnetic field strength indicates that the field is strong enough to hinder formation of a rotationally supported disk, which is inconsistent with the feature of a central toroidal field.

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.

The combination of high resolution and sensitivity offered by ALMA is revolutionizing our understanding of protoplanetary discs, as their bulk gas and dust distributions can be studied independently. In this paper we present resolved ALMA observations of the continuum emission ($\lambda=1.3$ mm) and CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $J=2-1$) integrated intensity from the disc around the nearby ($d = 162$ pc), intermediate mass ($M_{\star}=1.67\,M_{\odot}$) pre-main-sequence star CQ Tau. The data show an inner depression in continuum, and in both $^{13}$CO and C$^{18}$O emission. We employ a thermo-chemical model of the disc reproducing both continuum and gas radial intensity profiles, together with the disc SED. The models show that a gas inner cavity with size between 15 and 25 au is needed to reproduce the data with a density depletion factor between $\sim 10^{-1}$ and $\sim 10^{-3}$. The radial profile of the distinct cavity in the dust continuum is described by a Gaussian ring centered at $R_{\rm dust}=53\,$au and with a width of $\sigma=13\,$au. Three dimensional gas and dust numerical simulations of a disc with an embedded planet at a separation from the central star of $\sim20\,$au and with a mass of $\sim 6\textrm{-} 9\,M_{\rm Jup}$ reproduce qualitatively the gas and dust profiles of the CQ Tau disc. However, a one planet model appears not to be able to reproduce the dust Gaussian density profile predicted using the thermo-chemical modeling.

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.

"TRAO FUNS" is a project to survey Gould Belt's clouds in molecular lines. This paper presents its first results on the central region of the California molecular cloud, L1478. We performed On-The-Fly mapping observations using the Taedeok Radio Astronomy Observatory (TRAO) 14m single dish telescope equipped with a 16 multi-beam array covering $\sim$1.0 square degree area of this region using C$^{18}$O (1-0) mainly tracing low density cloud and about 460 square arcminute area using N$_{2}$H$^{+}$ (1-0) mainly tracing dense cores. CS (2-1) and SO $(3_{2}-2_{1})$ were also used simultaneously to map $\sim$440 square arcminute area of this region. We identified 10 filaments by applying the dendrogram technique to the C$^{18}$O data-cube and 8 dense N$_{2}$H$^{+}$ cores by using {\sc FellWalker}. Basic physical properties of filaments such as mass, length, width, velocity field, and velocity dispersion are derived. It is found that L1478 consists of several filaments with slightly different velocities. Especially the filaments which are supercritical are found to contain dense cores detected in N$_{2}$H$^{+}$. Comparison of non-thermal velocity dispersions derived from C$^{18}$O and N$_{2}$H$^{+}$ for the filaments and dense cores indicates that some of dense cores share similar kinematics with those of the surrounding filaments while several dense cores have different kinematics with those of their filaments. This suggests that the formation mechanism of dense cores and filaments can be different in individual filaments depending on their morphologies and environments.

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.

Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process.

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

We present high-resolution (~ 35 au) ALMA Band 6 1.3 mm dust polarization observations of IRAS 16293. These observations spatially resolve the dust polarization across the two protostellar sources and toward the filamentary structures between them. The dust polarization and inferred magnetic field have complicated structures throughout the region. In particular, we find that the magnetic field is aligned parallel to three filamentary structures. We characterize the physical properties of the filamentary structure that bridges IRAS 16293A and IRAS 16293B and estimate a magnetic field strength of 23-78 mG using the Davis-Chandrasekhar-Fermi method. We construct a toy model for the bridge material assuming that the young stars dominate the mass and gravitational potential of the system. We find that the expected gas flow to each star is of comparable order to the Alfven speed, which suggests that the field may be regulating the gas flow. We also find that the bridging material should be depleted in ~ 1000 yr. If the bridge is part of the natal filament that formed the stars, then it must have accreted new material. Alternatively, the bridge could be a transient structure. Finally, we show that the 1.3 mm polarization morphology of the optically thick IRAS 16293B system is qualitatively similar to dust self-scattering. Based on similar polarization measurements at 6.9 mm, we propose that IRAS 16293B has produced a substantial population of large dust grains with sizes between 200 and 2000 um.

This paper presents our observational attempts to precisely measure the central mass of a proto-brown dwarf candidate, L328-IRS, in order to investigate whether L328-IRS is in the substellar mass regime. Observations were made for the central region of L328-IRS with the dust continuum and CO isotopologue line emission at ALMA band 6, discovering the detailed outflow activities and a deconvolved disk structure of a size of $\sim$87 AU $\times$ 37 AU. We investigated the rotational velocities as a function of the disk radius, finding that its motions between 130 AU and 60 AU are partially fitted with a Keplerian orbit by a stellar object of $\sim0.30 ~\rm\,M_\odot$, while the motions within 60 AU do not follow any Keplerian orbit at all. This makes it difficult to lead a reliable estimation of the mass of L328-IRS. Nonetheless, our ALMA observations were useful enough to well constrain the inclination angle of the outflow cavity of L328-IRS as $\sim 66^{\circ}$, enabling us to better determine the mass accretion rate of $\sim 8.9\times 10^{-7} \rm\,M_\odot~yr^{-1}$. From assumptions that the internal luminosity of L328-IRS is mostly due to this mass accretion process in the disk, or that L328-IRS has mostly accumulated the mass through this constant accretion rate during its outflow activity, its mass was estimated to be $\sim 0.012 - 0.023\rm\,M_\odot$, suggesting L328-IRS to be a substellar object. However, we leave our identification of L328-IRS as a proto-brown dwarf to be tentative because of various uncertainties especially regarding the mass accretion rate.

The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts.

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 high resolution (~ 30 au) ALMA Band 6 dust polarization observations of VLA 1623. The VLA 1623 data resolve compact ~ 40 au inner disks around the two protobinary sources, VLA 1623-A and VLA 1623-B, and also an extended ~ 180 au ring of dust around VLA 1623-A. This dust ring was previously identified as a large disk in lower-resolution observations. We detect highly structured dust polarization toward the inner disks and the extended ring with typical polarization fractions ~ 1.7% and ~ 2.4%, respectively. The two components also show distinct polarization morphologies. The inner disks have uniform polarization angles aligned with their minor axes. This morphology is consistent with expectations from dust scattering. By contrast, the extended dust ring has an azimuthal polarization morphology not previously seen in lower-resolution observations. We find that our observations are well fit by a static, oblate spheroid model with a flux-frozen, poloidal magnetic field. We propose that the polarization traces magnetic grain alignment likely from flux freezing on large scales and magnetic diffusion on small scales. Alternatively, the azimuthal polarization may be attributed to grain alignment by the anisotropic radiation field. If the grains are radiatively aligned, then our observations indicate that large (~ 100 um) dust grains grow quickly at large angular extents. Finally, we identify significant proper motion of VLA 1623 using our observations and those in the literature. This result indicates that the proper motion of nearby systems must be corrected for when combining ALMA data from different epochs.

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 present high-sensitivity ($\sigma_I\sim0.2-0.5$ mJy, $\sigma_{QU}\sim0.05$ mJy), high-resolution (0.12-0.2 arcsec) observations of polarized 872 $\mu$m dust emission from the young multiple system VLA 1623 in $\rho$ Ophiuchus and the protostar L1527 in Taurus. We detect the circumstellar material of VLA 1623A, the extended Keplerian disk surrounding VLA 1623A which we call VLA 1623CBdisk, VLA 1623B, VLA 1623W, and L1527 strongly in the polarized emission, at the $\sim1-3\%$ level. We spatially resolve VLA 1623A into two sources, VLA 1623Aa and VLA 1623Ab, separated by $\sim30$ au and located within a cavity of radius $\sim50$ au within the circumbinary Keplerian disk, as well as the edge-on disk of VLA 1623W. The polarization angle of the emission is uniform across each protostellar source and nearly coincides with each disk's minor axis. The offsets between the minor axis and the polarization angle are not uniformly distributed at the $P\lesssim2\times10^{-4}$ level. The circumbinary disk surrounding VLA 1623Aab is azimuthally symmetrically polarized. Each compact source's emission is partially optically thick ($\tau\gtrsim1$) at 872 $\mu$m, complicating interpretations of polarization involving aligned grains. We find evidence against alignment by radiative flux in each source, particularly in the edge-on VLA 1623W and L1527. We detect astrometric offsets between the polarized emission and the total intensity in VLA 1623Aa, VLA 1623Ab, and VLA 1623B, as predicted if self-scattering in the optically thick limit operates. We conclude that self-scattering is likely responsible for disk-scale polarization at 872 $\mu$m in these systems.

We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence that the spectrum steepens with radius in the disk, manifested as a positive correlation between the observing frequency and the radius that encircles a fixed fraction of the emission ($R_{eff} \propto \nu^{0.34 \pm 0.08}$). The origins of this size--frequency relation are explored in the context of a theoretical framework for the growth and migration of disk solids. While that framework can reproduce a similar size--frequency relation, it predicts a steeper spectrum than is observed. Moreover, it comes closest to matching the data only on timescales much shorter ($\le 1$ Myr) than the putative UZ Tau age (~2-3 Myr). These discrepancies are the direct consequences of the rapid radial drift rates predicted by models of dust evolution in a smooth gas disk. One way to mitigate that efficiency problem is to invoke small-scale gas pressure modulations that locally concentrate drifting solids. If such particle traps reach high continuum optical depths at 30-340 GHz with a ~30-60% filling fraction in the inner disk ($r \lesssim20$ au), they can also explain the observed spatial gradient in the UZ Tau E disk spectrum.

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.

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle 2 observations of the 1.3 mm dust continuum emission of the protoplanetary disc surrounding the T Tauri star Elias 24 with an angular resolution of $\sim 0.2"$ ($\sim 28$ au). The dust continuum emission map reveals a dark ring at a radial distance of $0.47"$ ($\sim 65$ au) from the central star, surrounded by a bright ring at $0.58"$ ($\sim 81$ au). In the outer disc, the radial intensity profile shows two inflection points at $0.71"$ and $0.87"$ ($\sim 99$ and $121$ au respectively). We perform global three-dimensional smoothed particle hydrodynamic gas/dust simulations of discs hosting a migrating and accreting planet. Combining the dust density maps of small and large grains with three dimensional radiative transfer calculations, we produce synthetic ALMA observations of a variety of disc models in order to reproduce the gap- and ring-like features observed in Elias 24. We find that the dust emission across the disc is consistent with the presence of an embedded planet with a mass of $\sim 0.7\, \mathrm{M_{\mathrm{J}}}$ at an orbital radius of $\sim$ 60 au. Our model suggests that the two inflection points in the radial intensity profile are due to the inward radial motion of large dust grains from the outer disc. The surface brightness map of our disc model provides a reasonable match to the gap- and ring-like structures observed in Elias 24, with an average discrepancy of $\sim$ 5% of the observed fluxes around the gap region.

The mechanism for producing polarized emission from protostellar disks at (sub)millimeter wavelengths is currently uncertain. Classically, polarization is expected from non-spherical grains aligned with the magnetic field. Recently, two alternatives have been suggested. One polarization mechanism is caused by self-scattering from dust grains of sizes comparable to the wavelength while the other mechanism is due to grains aligned with their short axes along the direction of radiation anisotropy. The latter has recently been shown as a likely mechanism for causing the dust polarization detected in HL Tau at 3.1 mm. In this paper, we present ALMA polarization observations of HL Tau for two more wavelengths: 870 $\mu$m and 1.3 mm. The morphology at 870 $\mu$m matches the expectation for self-scattering, while that at 1.3 mm shows a mix between self-scattering and grains aligned with the radiation anisotropy. The observations cast doubt on the ability of (sub)millimeter continuum polarization to probe disk magnetic fields for at least HL Tau. By showing two distinct polarization morphologies at 870 $\mu$m and 3.1 mm and a transition between the two at 1.3 mm, this paper provides definitive evidence that the dominant (sub)millimeter polarization mechanism transitions with wavelength. In addition, if the polarization at 870 $\mu$m is due to scattering, the lack of polarization asymmetry along the minor axis of the inclined disk implies that the large grains responsible for the scattering have already settled into a geometrically thin layer, and the presence of asymmetry along the major axis indicates that the HL Tau disk is not completely axisymmetric.

{Context}. The HIFI instrument on the Herschel Space Observatory performed over 9100 astronomical observations, almost 900 of which were calibration observations in the course of the nearly four-year Herschel mission. The data from each observation had to be converted from raw telemetry into calibrated products and were included in the Herschel Science Archive. {Aims}. The HIFI pipeline was designed to provide robust conversion from raw telemetry into calibrated data throughout all phases of the HIFI missions. Pre-launch laboratory testing was supported as were routine mission operations. {Methods}. A modular software design allowed components to be easily added, removed, amended and/or extended as the understanding of the HIFI data developed during and after mission operations. {Results}. The HIFI pipeline processed data from all HIFI observing modes within the Herschel automated processing environment as well as within an interactive environment. The same software can be used by the general astronomical community to reprocess any standard HIFI observation. The pipeline also recorded the consistency of processing results and provided automated quality reports. Many pipeline modules were in use since the HIFI pre-launch instrument level testing. {Conclusions}. Processing in steps facilitated data analysis to discover and address instrument artefacts and uncertainties. The availability of the same pipeline components from pre-launch throughout the mission made for well-understood, tested, and stable processing. A smooth transition from one phase to the next significantly enhanced processing reliability and robustness.

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

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk.

Protostellar evolution, following the formation of the protostar is becoming reasonably well characterized, but the evolution from a prestellar core to a protostar is not well known, although the first hydrostatic core (FHSC) must be a pivotal step. NGC 1333 IRAS 4C is a potentially very young object, that we directly compare with the nearby Class 0 IRAS 4A and IRAS 4B. Observational constraints are provided by spectral imaging from the JCMT Spectral Legacy Survey (330-373 GHz) and continuum and line observations from CARMA. We present integrated intensity and velocity maps of several species, including CO, H2CO and CH3OH. The velocity of an observed outflow, the degree of CO depletion, the deuteration of DCO+/HCO+ and gas kinetic temperatures are observational signatures that we present. We report differences between the three sources in four aspects: a) the kinetic temperature is much lower towards IRAS 4C, b) the line profiles of the detected species show strong outflow activity towards IRAS 4A and IRAS 4B but not towards IRAS 4C, c) the HCN/HNC is < 1 towards IRAS 4C, which confirms the cold nature of the source, d) the degree of CO depletion and the deuteration are the lowest towards the warmest of the sources, IRAS 4B. IRAS 4C seems to be in a different evolutionary state than the IRAS 4A and IRAS 4B sources. We can probably exclude the FHSC stage due to the relatively low Lsmm/Lbol (~6\%) and we investigate the earliest accretion phase of Class 0 stage and the transition between Class 0 to Class I. Our results do not show a consistent scenario for either case, with the major issue being the absence of outflow activity and the cold nature of IRAS 4C. The number of FHSC candidates in Perseus is ~10 times higher than current models predict, which suggests that the lifespan of these objects is greater than 1000 years, possibly due to an accretion rate lower than 4x10^-5 msun/yr.

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.

To characterize the mechanisms of planet formation it is crucial to investigate the properties and evolution of protoplanetary disks around young stars, where the initial conditions for the growth of planets are set. Our goal is to study grain growth in the disk of the young, intermediate mass star HD163296 where dust processing has already been observed, and to look for evidence of growth by ice condensation across the CO snowline, already identified in this disk with ALMA. Under the hypothesis of optically thin emission we compare images at different wavelengths from ALMA and VLA to measure the opacity spectral index across the disk and thus the maximum grain size. We also use a Bayesian tool based on a two-layer disk model to fit the observations and constrain the dust surface density. The measurements of the opacity spectral index indicate the presence of large grains and pebbles ($\geq$1 cm) in the inner regions of the disk (inside $\sim$50 AU) and smaller grains, consistent with ISM sizes, in the outer disk (beyond 150 AU). Re-analysing ALMA Band 7 Science Verification data we find (radially) unresolved excess continuum emission centered near the location of the CO snowline at $\sim$90 AU. Our analysis suggests a grain size distribution consistent with an enhanced production of large grains at the CO snowline and consequent transport to the inner regions. Our results combined with the excess in infrared scattered light found by Garufi et al. (2014) suggests the presence of a structure at 90~AU involving the whole vertical extent of the disk. This could be evidence for small scale processing of dust at the CO snowline.

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.

Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.

A massive envelope and a strong bipolar outflow are the two main structures characterizing the youngest protostellar systems. In order to understand the physical properties of a bipolar outflow and the relationship with those of the envelope, we obtained a mosaic map covering the whole bipolar outflow of the youngest protostellar system L1157 with about $5"$ angular resolution in CO J=2-1 using the Combined Array for Research in Millimeter-wave Astronomy. By utilizing these observations of the whole bipolar outflow, we estimate its physical properties and show that they are consistent with multiple jets. We also constrain a preferred precession direction. In addition, we observed the central envelope structure with $2"$ resolution in the $\lambda=1.3$ and 3 mm continua and various molecular lines: C$^{17}$O, C$^{18}$O, $^{13}$CO, CS, CN, N$_2$H$^+$, CH$_3$OH, H$_2$O, SO, and SO$_2$. All the CO isotopes and CS, CN, and N$_2$H$^+$ have been detected and imaged. We marginally detected the features that can be interpreted as a rotating inner envelope in C$^{17}$O and C$^{18}$O and as an infalling outer envelope in N$_2$H$^+$. We also estimated the envelope and central protostellar masses and found that the dust opacity spectral index changes with radius.

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at {\lambda} = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at {\lambda} = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, {\beta}, is consistent with the observed dust emission in both disks, with low-{\beta} in the inner disk and high-{\beta} in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

We present continuum observations at 1.3 and 2.7 mm using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) toward six protoplanetary disks in the Taurus molecular cloud: CI Tau, DL Tau, DO Tau, FT Tau, Haro 6-13, and HL Tau. We constrain physical properties of the disks with Bayesian inference using two disk models; flared power-law disk model and flared accretion disk model. Comparing the physical properties, we find that the more extended disks are less flared and that the dust opacity spectral index (beta) is smaller in the less massive disks. In addition, disks with a steeper mid-plane density gradient have a smaller beta, which suggests that grains grow and radially move. Furthermore, we compare the two disk models quantitatively and find that the accretion disk model provides a better fit overall. We also discuss the possibilities of substructures on three extended protoplanetary disks.

We present a CARMA 1.3 mm continuum survey toward 9 Class 0 protostars in the Perseus molecular cloud at $\sim$0.3$^{\prime\prime}$ (70 AU) resolution. This study approximately doubles the number of Class 0 protostars observed with spatial resolutions $<$ 100 AU at millimeter wavelengths, enabling the presence of protostellar disks and proto-binary systems to be probed. We detect flattened structures with radii $>$ 100 AU around 2 sources (L1448 IRS2 and Per-emb-14) and these sources may be strong disk candidates. Marginally-resolved structures with position angles within 30$^{\circ}$ of perpendicular to the outflow are found toward 3 protostars (L1448 IRS3C, IRAS 03282+3035, and L1448C) and are considered disk candidates. Two others (L1448 IRS3B and IRAS 03292+3039) have resolved structure, possibly indicative of massive inner envelopes or disks; L1448 IRS3B also has a companion separated by 0.9$^{\prime\prime}$ ($\sim$210 AU). IC348-MMS does not have well-resolved structure and the candidate first hydrostatic core L1451-MMS is marginally resolved on 1$^{\prime\prime}$ scales. The strong disk candidate sources were followed-up with C$^{18}$O ($J=2\rightarrow1$) observations, detecting velocity gradients consistent with rotation, but it is unclear if the rotation is Keplerian. We compare the observed visibility amplitudes to radiative transfer models, finding that visibility amplitude ratios suggest a compact component (possibly a disk) is necessary for 5 of 9 Class 0 sources; envelopes alone may explain the other 4 systems. We conclude that there is evidence for the formation of large disks in the Class 0 phase with a range of radii and masses dependent upon their initial formation conditions.

We present subarcsecond (~0.35") resolved observations of the 1.3 mm dust polarization from the edge-on circumstellar disk around the Class 0 protostar L1527. The inferred magnetic field is consistent with a dominantly toroidal morphology; there is no significantly detected vertical poloidal component to which observations of an edge-on disk are most sensitive. This suggests that angular momentum transport in Class 0 protostars (when large amounts of material are fed down to the disk from the envelope and accreted onto the protostar) is driven mainly by magnetorotational instability rather than magnetocentrifugal winds at 50 AU scales. In addition, with the data to date there is an early, tentative trend that R>30 AU disks have so far been found in Class 0 systems with average magnetic fields on the 1000 AU scale strongly misaligned with the rotation axis. The absence of such a disk in the aligned case could be due to efficient magnetic braking that disrupts disk formation. If this is the case, this implies that candidate Class 0 disk systems could be identified by the average magnetic field direction at ~1000 AU spatial scales.

For the Class 0 protostar, L1527, we compare 131 polarization vectors from SCUPOL/JCMT, SHARP/CSO and TADPOL/CARMA observations with the corresponding model polarization vectors of four ideal-MHD, non-turbulent, cloud core collapse models. These four models differ by their initial magnetic fields before collapse; two initially have aligned fields (strong and weak) and two initially have orthogonal fields (strong and weak) with respect to the rotation axis of the L1527 core. Only the initial weak orthogonal field model produces the observed circumstellar disk within L1527. This is a characteristic of nearly all ideal-MHD, non-turbulent, core collapse models. In this paper we test whether this weak orthogonal model also has the best agreement between its magnetic field structure and that inferred from the polarimetry observations of L1527. We found that this is not the case; based on the polarimetry observations the most favored model of the four is the weak aligned model. However, this model does not produce a circumstellar disk, so our result implies that a non-turbulent, ideal-MHD global collapse model probably does not represent the core collapse that has occurred in L1527. Our study also illustrates the importance of using polarization vectors covering a large area of a cloud core to determine the initial magnetic field orientation before collapse; the inner core magnetic field structure can be highly altered by a collapse and so measurements from this region alone can give unreliable estimates of the initial field configuration before collapse.

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

Magnetic fields in accretion disks play a dominant role during the star formation process but have hitherto been observationally poorly constrained. Field strengths have been inferred on T Tauri stars themselves and possibly in the innermost part of the accretion disk, but the strength and morphology of the field in the bulk of the disk have not been observed. Unresolved measurements of polarized emission (arising from elongated dust grains aligned perpendicular to the field) imply average fields aligned with the disks. Theoretically, the fields are expected to be largely toroidal, poloidal, or a mixture of the two, which imply different mechanisms for transporting angular momentum in the disks of actively accreting young stars such as HL Tau. Here we report resolved measurements of the polarized 1.25 mm continuum emission from HL Tau's disk. The magnetic field on a scale of 80 AU is coincident with the major axis (~210 AU diameter) of the disk. From this we conclude that the magnetic field inside the disk at this scale cannot be dominated by a vertical component, though a purely toroidal field does not fit the data well either. The unexpected morphology suggests that the magnetic field's role for the accretion of a T Tauri star is more complex than the current theoretical understanding.

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

We present the N2H+(J=1-0) map of the Serpens South molecular cloud obtained as part of the CARMA Large Area Star Formation Survey (CLASSy). The observations cover 250 square arcminutes and fully sample structures from 3000 AU to 3 pc with a velocity resolution of 0.16 km/s, and they can be used to constrain the origin and evolution of molecular cloud filaments. The spatial distribution of the N2H+ emission is characterized by long filaments that resemble those observed in the dust continuum emission by Herschel. However, the gas filaments are typically narrower such that, in some cases, two or three quasi-parallel N2H+ filaments comprise a single observed dust continuum filament. The difference between the dust and gas filament widths casts doubt on Herschel ability to resolve the Serpens South filaments. Some molecular filaments show velocity gradients along their major axis, and two are characterized by a steep velocity gradient in the direction perpendicular to the filament axis. The observed velocity gradient along one of these filaments was previously postulated as evidence for mass infall toward the central cluster, but these kind of gradients can be interpreted as projection of large-scale turbulence.

We present {\lambda}1.3 mm CARMA observations of dust polarization toward 30 star-forming cores and 8 star-forming regions from the TADPOL survey. We show maps of all sources, and compare the ~2.5" resolution TADPOL maps with ~20" resolution polarization maps from single-dish submillimeter telescopes. Here we do not attempt to interpret the detailed B-field morphology of each object. Rather, we use average B-field orientations to derive conclusions in a statistical sense from the ensemble of sources, bearing in mind that these average orientations can be quite uncertain. We discuss three main findings: (1) A subset of the sources have consistent magnetic field (B-field) orientations between large (~20") and small (~2.5") scales. Those same sources also tend to have higher fractional polarizations than the sources with inconsistent large-to-small-scale fields. We interpret this to mean that in at least some cases B-fields play a role in regulating the infall of material all the way down to the ~1000 AU scales of protostellar envelopes. (2) Outflows appear to be randomly aligned with B-fields; although, in sources with low polarization fractions there is a hint that outflows are preferentially perpendicular to small-scale B-fields, which suggests that in these sources the fields have been wrapped up by envelope rotation. (3) Finally, even at ~2.5" resolution we see the so-called "polarization hole" effect, where the fractional polarization drops significantly near the total intensity peak. All data are publicly available in the electronic edition of this article.

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

We present observations of three Class 0/I protostars (L1157-mm, CB230 IRS1, and L1165-SMM1) using the Karl G. Jansky Very Large Array (VLA) and observations of two (L1165-SMM1 and CB230 IRS1) with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The VLA observations were taken at wavelengths of $\lambda = 7.3$ mm, 1.4 cm, 3.3 cm, 4.0 cm, and 6.5 cm with a best resolution of $\sim$0\farcs06 (18 AU) at 7.3 mm. The L1165-SMM1 CARMA observations were taken at $\lambda = 1.3$ mm with a best resolution of $\sim0\farcs3$ (100 AU), and the CB230 IRS1 observations were taken at $\lambda = 3.4$ mm with a best resolution of $\sim$3\arcsec\ (900 AU). We find that L1165-SMM1 and CB230 IRS1 have probable binary companions at separations of $\sim$0\farcs3 (100 AU) from detections of secondary peaks at multiple wavelengths. The position angles of these companions are nearly orthogonal to the direction of the observed bipolar outflows, consistent with the expected protostellar disk orientations. We suggest that these companions may have formed from disk fragmentation; turbulent fragmentation would not preferentially arrange the binary companions to be orthogonal to the outflow direction. For L1165-SMM1, both the 7.3 mm and 1.3 mm emission show evidence of a large (R $>$ 100 AU) disk. For the L1165-SMM1 primary protostar and the CB230 IRS1 secondary protostar, the 7.3 mm emission is resolved into structures consistent with $\sim20$ AU radius disks. For the other protostars, including L1157-mm, the emission is unresolved, suggesting disks with radii $< 20$ AU.

We present 350 micron polarization observations of four low-mass cores containing Class 0 protostars: L483, L1157, L1448-IRS2, and Serp-FIR1. This is the second paper in a larger survey aimed at testing magnetically regulated models for core-collapse. One key prediction of these models is that the mean magnetic field in a core should be aligned with the symmetry axis (minor axis) of the flattened YSO inner envelope (aka pseudodisk). Furthermore, the field should exhibit a pinched or hour-glass shaped morphology as gravity drags the field inward towards the central protostar. We combine our results for the four cores with results for three similar cores that were published in the first paper from our survey. An analysis of the 350 micron polarization data for the seven cores yields evidence of a positive correlation between mean field direction and pseudodisk symmetry axis. Our rough estimate for the probability of obtaining by pure chance a correlation as strong as the one we found is about 5%. In addition, we combine together data for multiple cores to create a source-averaged magnetic field map having improved signal-to-noise ratio, and this map shows good agreement between mean field direction and pseudodisk axis (they are within 15 degrees). We also see hints of a magnetic pinch in the source-averaged map. We conclude that core-scale magnetic fields appear to be strong enough to guide gas infall, as predicted by the magnetically regulated models. Finally, we find evidence of a positive correlation between core magnetic field direction and bipolar outflow axis.

We present results of 1.3 mm dust polarization observations toward 16 nearby, low-mass protostars, mapped with ~2.5" resolution at CARMA. The results show that magnetic fields in protostellar cores on scales of ~1000 AU are not tightly aligned with outflows from the protostars. Rather, the data are consistent with scenarios where outflows and magnetic fields are preferentially misaligned (perpendicular), or where they are randomly aligned. If one assumes that outflows emerge along the rotation axes of circumstellar disks, and that the outflows have not disrupted the fields in the surrounding material, then our results imply that the disks are not aligned with the fields in the cores from which they formed.

We present the first detection of polarization around the Class 0 low-mass protostar L1157-mm at two different wavelengths. We show polarimetric maps at large scales (10" resolution at 350 um) from the SHARC-II Polarimeter and at smaller scales (1.2"-4.5" at 1.3 mm) from the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The observations are consistent with each other and show inferred magnetic field lines aligned with the outflow. The CARMA observations suggest a full hourglass magnetic field morphology centered about the core; this is only the second well-defined hourglass detected around a low-mass protostar to date. We apply two different methods to CARMA polarimetric observations to estimate the plane-of-sky magnetic field magnitude, finding values of 1.4 and 3.4 mG.

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity \kappa_{\nu}. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (\beta), we find that \beta(R) increases from \beta<0.5 at \sim20 AU to \beta>1.5 for R>80 AU, inconsistent with a constant value of \beta\ across the disk (at the 10\sigma\ level). Furthermore, if radial variations of \kappa_{\nu} are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.

We present high-resolution CARMA 230 GHz continuum imaging of nine deeply embedded protostars in the Serpens Molecular Cloud, including six of the nine known Class 0 protostars in Serpens. This work is part of a program to characterize disk and envelope properties for a complete sample of Class 0 protostars in nearby low-mass star forming regions. Here we present CARMA maps and visibility amplitudes as a function of uv-distance for the Serpens sample. Observations are made in the B, C, D, and E antenna configurations, with B configuration observations utilizing the CARMA Paired Antenna Calibration System. Combining data from multiple configurations provides excellent uv-coverage (4-500 klam), allowing us to trace spatial scales from 1e2 to 1e4 AU. We find evidence for compact disk components in all of the observed Class 0 protostars, suggesting that disks form at very early times (t<0.2 Myr) in Serpens. We make a first estimate of disk masses using the flux at 50 klam, where the contribution from the envelope should be negligible, assuming an unresolved disk. The resulting disk masses range from 0.04 Msun to 1.7 Msun, with a mean of approximately 0.2 Msun. Our high resolution maps are also sensitive to binary or multiple sources with separations > 250 AU, but significant evidence of multiplicity on scales <2000 AU is seen in only one source.

We present results of high-resolution imaging toward HL Tau by the Combined Array for Research in Millimeter-wave Astronomy (CARMA). We have obtained 1.3 and 2.7 mm dust continua with an angular resolution down to 0.13 arc second. Through model fitting to the two wavelength data simultaneously in Bayesian inference using a flared viscous accretion disk model, we estimate the physical properties of HL Tau, such as density distribution, dust opacity spectral index, disk mass, disk size, inclination angle, position angle, and disk thickness. HL Tau has a circumstellar disk mass of 0.13 solar mass, a characteristic radius of 79 AU, an inclination of 40 degree, and a position angle of 136 degree. Although a thin disk model is preferred by our two wavelength data, a thick disk model is needed to explain the high mid- and far-infrared emission of the HL Tau spectral energy distribution. This could imply large dust grains settled down on the mid plane with fine dust grains mixed with gas. The HL Tau disk is likely gravitationally unstable and can be fragmented between 50 and 100 AU of radius. However, we did not detect dust thermal continuum supporting the protoplanet candidate claimed by a previous study using observations of the Very Large Array at 1.3 cm.

We present 0.15" resolution observations of the 227 GHz continuum emission from the circumstellar disk around the FU-Orionis star PP 13S*. The data were obtained with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) Paired Antenna Calibration System (C-PACS), which measures and corrects the atmospheric delay fluctuations on the longest baselines of the array in order to improve the sensitivity and angular resolution of the observations. A description of the C-PACS technique and the data reduction procedures are presented. C-PACS was applied to CARMA observations of PP 13S*, which led to a factor of 1.6 increase in the observed peak flux of the source, a 36% reduction in the noise of the image, and a 52% decrease in the measured size of the source major axis. The calibrated complex visibilities were fitted with a theoretical disk model to constrain the disk surface density. The total disk mass from the best fit model corresponds to 0.06 \msun, which is larger than the median mass of a disk around a classical T Tauri star. The disk is optically thick at a wavelength of 1.3 mm for orbital radii less than 48 AU. At larger radii, the inferred surface density of the PP 13S* disk is an order of magnitude lower than that needed to develop a gravitational instability.

We present dust opacity spectral indexes (beta) of the youngest protostellar systems (so-called Class 0 sources), L1448 IRS 2, L1448 IRS 3, and L1157, obtained between 1.3 mm and 2.7 mm continua, using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The unprecedented compact configuration and image fidelity of CARMA allow a better detection of the dust continuum emission from Class 0 sources, with a less serious missing flux problem normally associated with interferometry. Through visibility-modeling at both 1.3 mm and 2.7 mm simultaneously, as well as image- and visibility-comparison, we show that beta of the three Class 0 sources are around or smaller than 1, indicating that dust grains have already significantly grown at the Class 0 stage. In addition, we find a radial dependence of beta, which implies faster grain growth in the denser central regions and/or dust segregation. Density distributions of the Class 0 sources are also addressed by visibility-modeling.

Deep Spitzer IRAC images of L1157 reveal many of the details of the outflow and the circumstellar environment of this Class 0 protostar. In IRAC band 4, 8 microns, there is a flattened structure seen in absorption against the background emission. The structure is perpendicular to the outflow and is extended to a diameter of 2 arcminutes. This structure is the first clear detection of a flattened circumstellar envelope or pseudo-disk around a Class 0 protostar. Such a flattened morphology is an expected outcome for many collapse theories that include magnetic fields or rotation. We construct an extinction model for a power-law density profile, but we do not constrain the density power-law index.

We present deep IRAC images that highlight the scattered light emission around many of the youngest protostars, the so-called Class 0 sources, in L1448. By comparison of the data with a Monte Carlo radiative transfer code \citep{whitney2003a}, we demonstrate for the first time that the observed infrared light from these objects is consistent with scattered light from the central protostar. The scattered light escapes out the cavity, carved by molecular outflows, in the circumstellar envelope. In particular, we observe prominent scattered light nebulae associated with the Class 0 sources: L1448-mm, L1448 IRS 2, and 3B, as well as a Class I source: IRS 3A. We use a grid of models with probable protostellar properties to generate model spectral energy distributions (SEDs) and images for bands sensitive to this scattered light: J, H, Ks, and \textit{Spitzer} IRAC bands. By simultaneously fitting SEDs and images of the outflow cavities, we are able to model geometric parameters, i.e. inclination angle and opening angle, and loosely constrain physical parameters. The opening angle may be an important indicator of the evolutionary state of a source. We compare our results for Class 0 sources to similar studies of Class I sources. There may be a transition phase from Class 0 to Class I when a source has an opening angle between 20$^o$ to 30$^o$. It is important to note that while the best fit model image and SED do not fully describe the sources, the fits generally describe the circumstellar structure of Class 0 sources in L1448.

We performed spectral line observations of CO J=2-1, 13CO J=1-0, and C18O J=1-0 and polarimetric observations in the 1.3 mm continuum and CO J=2-1 toward a multiple protostar system, L1448 IRS 3, in the Perseus molecular complex at a distance of ~250 pc, using the BIMA array. In the 1.3 mm continuum, two sources (IRS 3A and 3B) were clearly detected with estimated envelope masses of 0.21 and 1.15 solar masses, and one source (IRS 3C) was marginally detected with an upper mass limit of 0.03 solar masses. In CO J=2-1, we revealed two outflows originating from IRS 3A and 3B. The masses, mean number densities, momentums, and kinetic energies of outflow lobes were estimated. Based on those estimates and outflow features, we concluded that the two outflows are interacting and that the IRS 3A outflow is nearly perpendicular to the line of sight. In addition, we estimated the velocity, inclination, and opening of the IRS 3B outflow using Bayesian statistics. When the opening angle is ~20 arcdeg, we constrain the velocity to ~45 km/s and the inclination angle to ~57 arcdeg. Linear polarization was detected in both the 1.3 mm continuum and CO J=2-1. The linear polarization in the continuum shows a magnetic field at the central source (IRS 3B) perpendicular to the outflow direction, and the linear polarization in the CO J=2-1 was detected in the outflow regions, parallel or perpendicular to the outflow direction. Moreover, we comprehensively discuss whether the binary system of IRS 3A and 3B is gravitationally bound, based on the velocity differences detected in 13CO J=1-0 and C18O J=1-0 observations and on the outflow features. The specific angular momentum of the system was estimated as ~3e20 cm^2/s, comparable to the values obtained from previous studies on binaries and molecular clouds in Taurus.