We report ALMA Cycle 4 observations of the Class I binary protostellar system L1551 IRS 5 in the 0.9-mm continuum emission, C18O (J=3-2), OCS (J=28-27), and four other Band 7 lines. At ~0.07" (= 10 au) resolution in the 0.9 mm emission, two circumstellar disks (CSDs) associated with the binary protostars are separated from the circumbinary disk (CBD). The CBD is resolved into two spiral arms, one connecting to the CSD around the northern binary source, Source N, and the other to Source S. As compared to the CBD in the neighboring protobinary system L1551 NE, the CBD in L1551 IRS 5 is more compact (r ~150 au) and the m=1 mode of the spirals found in L1551 NE is less obvious in L1551 IRS 5. Furthermore, the dust and molecular-line brightness temperatures of CSDs and CBD reach >260 K and >100 K, respectively, in L1551 IRS 5, much hotter than those in L1551 NE. The gas motions in the spiral arms are characterized by rotation and expansion. Furthermore, the transitions from the CBD to the CSD rotations at around the L2 and L3 Lagrangian points and gas motions around the L1 point are identified. Our numerical simulations reproduce the observed two spiral arms and expanding gas motion as a result of gravitational torques from the binary, transitions from the CBD to the CSD rotations, and the gas motion around the L1 point. The higher temperature in L1551 IRS 5 likely reflects the inferred FU-Ori event.

An increasing number of young circumstellar disks show strikingly ordered (sub)millimeter polarization orientations along the minor axis, which is strong evidence for polarization due to scattering by ~0.1 mm sized grains. To test this mechanism further, we model the dust continuum and polarization data of HD 163296, one of the best observed disks with prominent rings and gaps, using the RADMC-3D radiative transfer code. We find that scattering by grains with a maximum size of 90$\mu$m can simultaneously reproduce the polarization observed at 0.87 mm (ALMA Band 7) and the unusually low spectral index of $\alpha$ ~ 1.5 between 0.87 and 1.25 mm (ALMA Band 6) in the optically thick inner disk as a result of more efficient scattering at a shorter wavelength. The relatively low spectral index of $\alpha$ ~ 2.5 inferred for the optically thin gaps is also reproduced by the same (relatively small) grains, as a result of telescope beam averaging of the gaps (with an intrinsic $\alpha$ ~ 4) and their adjacent optically thick rings (where $\alpha$ << 2). In this case, the long-standing tension between the grain sizes inferred from polarization and spectral index disappears because the relatively low $\alpha$ values are illusory and do not require large mm-sized grains. In addition, the polarization fraction has a unique pattern of azimuthal variation: higher along the major axis than the minor axis in the gaps but higher along the minor axis than the major axis in the rings. We find a rapidly declining polarization spectrum (with the fraction $p \propto \lambda^{-3}$ approximately) in the gaps, which becomes flattened or even inverted towards short wavelengths in the optically thick rings. These contrasting behaviors in the rings and gaps provide further tests of scattering-induced polarization that can be tested via multi-wavelength observations that resolve the disk substructure.

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

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

It has been asserted that the primary star in the HD 189733 system steadily accretes evaporated exospheric gases from its ``hot Jupiter'' companion, rather like a T Tauri star accreting from a disk. We conduct statistical and periodogram analyses of the photometric time series of the primary, as acquired by the automated photoelectric telescope (APT), Microvariability and Oscillations of Stars (MOST), and Wise Observatory, to investigate this claim with the goal of revealing the presence of accretion shocks or photospheric accretion hotspots as are found in T Tauri systems such as AA Tau. None of the anticipated features were found. We re-analyze existing radio, optical, ultraviolet, and X-ray data within the framework of accreting T Tauri systems to determine physical quantities such as plasma density and temperature, accretion rate, and flare lengths. We find that with an accretion rate of $\dot{M}\sim10^{9}$ to 10$^{11}$ g s$^{-1}$, the star is more similar to a system that intermittently absorbs gas from sungrazing comets in outburst than classical T Tauri systems, which have accretion rates at least two orders of magnitude larger. If such accretion exists, it would result in undetectably low activity at all wavelengths. Alternatively, all of the emission properties observed thus far are in agreement with stellar activity from a magnetically active star.

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $\mu$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $\mu$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

We present the first results of high-spectral resolution (0.023 km/s) N$_2$H$^+$ observations of dense gas dynamics at core scales (~0.01 pc) using the recently commissioned Argus instrument on the Green Bank Telescope (GBT). While the fitted linear velocity gradients across the cores measured in our targets nicely agree with the well-known power-law correlation between the specific angular momentum and core size, it is unclear if the observed gradients represent core-scale rotation. In addition, our Argus data reveal detailed and intriguing gas structures in position-velocity (PV) space for all 5 targets studied in this project, which could suggest that the velocity gradients previously observed in many dense cores actually originate from large-scale turbulence or convergent flow compression instead of rigid-body rotation. We also note that there are targets in this study with their star-forming disks nearly perpendicular to the local velocity gradients, which, assuming the velocity gradient represents the direction of rotation, is opposite to what is described by the classical theory of star formation. This provides important insight on the transport of angular momentum within star-forming cores, which is a critical topic on studying protostellar disk formation.

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.

We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $\mu$m) and D (154 $\mu$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_DC and the molecular hydrogen column density across the cloud. A positive slope (R_DC > 1) dominates the lower density and well illuminated portions of the cloud, that are heated by the high mass star Oph S1, whereas a transition to a negative slope (R_DC < 1) is observed toward the denser and less evenly illuminated cloud core. We interpret the trends as due to a combination of: (1) Warm grains at the cloud outskirts, which are efficiently aligned by the abundant exposure to radiation from Oph S1, as proposed in the radiative torques theory; and (2) Cold grains deep in the cloud core, which are poorly aligned due to shielding from external radiation. To assess this interpretation, we developed a very simple toy model using a spherically symmetric cloud core based on Herschel data, and verified that the predicted variation of R_DC is consistent with the observations. This result introduces a new method that can be used to probe the grain alignment efficiency in molecular clouds, based on the analysis of trends in the far-infrared polarization spectrum.

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.

Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the building blocks for planets, and changes in dust grain properties will affect all of these phenomena. A full understanding of dust in all of its forms and stages requires a multi-disciplinary investigation of the dust life cycle. Such an investigation can be achieved with a statistical study of dust properties across stellar evolution, star and planet formation, and redshift. Current and future instrumentation will enable this investigation through fast and sensitive observations in dust continuum, polarization, and spectroscopy from near-infrared to millimeter wavelengths.

Significant advances have been made over the past decade in the characterization of multiple protostar systems, enabled by the Karl G. Jansky Very Large Array (VLA), high-resolution infrared observations with the Hubble Space Telescope, and ground-based facilities. To further understand the mechanism(s) of multiple star formation, a combination of statistics, high-angular resolution radio/millimeter continuum imaging, characterization of kinematic structure, magnetic fields via polarimetry, and comparison with numerical simulations are needed. Thus, understanding the origin of stellar multiplicity in different regimes of companion separation will soon be within reach. However, to overcome challenges that studies in this field are now confronted with, a range of new capabilities are required: a new millimeter/centimeter wave facility with 10 mas resolution at {\lambda}=1 cm, space-based near to far-infrared observatories, continued development of low to high resolution spectroscopy on 3m to 10m class telescopes, and an ELT-class telescope with near to mid-infrared imaging/spectroscopic capability.

Knowledge of protostellar evolution has been revolutionized with the advent of surveys at near-infrared to submillimeter wavelengths. This has enabled the bolometric luminosities and bolometric temperatures (traditional protostellar evolution diagnostics) to be measured for large numbers of protostars. However, further progress is difficult without knowing the masses of the central protostars. Protostar masses can be most accurately determined via molecular line kinematics from millimeter interferometers (i.e., ALMA). Theoretical investigations have predicted the protostellar mass function (PMF) for various protostellar mass accretion models, and it is now imperative to observationally constrain its functional form. While ALMA has enabled protostellar mass measurements, samples approaching 100 sources are necessary to constrain the functional form of the PMF, and upgrades to ALMA and/or a new mm/cm facility will increase the feasibility of measuring such a large number of protostar masses. The masses of protostars will enable their stellar structure (radius and intrinsic luminosity), evolution, and accretion histories to be better understood. This is made more robust when effective temperatures and accretion rates can be measured via ground/space-based near to mid-infrared spectroscopy. Furthermore, access to supercomputing facilities is essential to fit the protostar masses via radiative transfer modeling and updated theoretical/numerical modeling of stellar structure may also be required.

The majority of the ultimate main-sequence mass of a star is assembled in the protostellar phase, where a forming star is embedded in an infalling envelope and encircled by a protoplanetary disk. Studying mass accretion in protostars is thus a key to understanding how stars gain their mass and ultimately how their disks and planets form and evolve. At this early stage, the dense envelope reprocesses most of the luminosity generated by accretion to far-infrared and submillimeter wavelengths. Time-domain photometry at these wavelengths is needed to probe the physics of accretion onto protostars, but variability studies have so far been limited, in large part because of the difficulty in accessing these wavelengths from the ground. We discuss the scientific progress that would be enabled with far-infrared and submillimeter programs to probe protostellar variability in the nearest kiloparsec.

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 report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

We present Far-Infrared polarimetry observations of M82 at 53 and $154~\mu \rm{m}$ and NGC 253 at $89~\mu \rm{m}$, which were taken with HAWC+ in polarimetry mode on the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization of M82 at $53~\mu \rm{m}$ clearly shows a magnetic field geometry perpendicular to the disk in the hot dust emission. For M82 the polarization at $154~\mu \rm{m}$ shows a combination of field geometry perpendicular to the disk in the nuclear region, but closer to parallel to the disk away from the nucleus. The fractional polarization at $53~\mu \rm{m}$ $(154~\mu \rm{m})$ ranges from 7% (3%) off nucleus to 0.5% (0.3%) near the nucleus. A simple interpretation of the observations of M82 invokes a massive polar outflow, dragging the field along, from a region $\sim 700$~pc in diameter that has entrained some of the gas and dust, creating a vertical field geometry seen mostly in the hotter $(53~\mu \rm{m})$ dust emission. This outflow sits within a larger disk with a more typical planar geometry that more strongly contributes to the cooler $(154~\mu \rm{m})$ dust emission. For NGC 253, the polarization at $89~\mu \rm{m}$ is dominated by a planar geometry in the tilted disk, with weak indication of a vertical geometry above and below the plane from the nucleus. The polarization observations of NGC 253 at $53~\mu \rm{m}$ were of insufficient S/N for detailed analysis.

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.

Disk polarization in (sub)millimeter dust continuum is a rapidly growing field in the ALMA era. It opens up the exciting possibility of detecting and characterizing magnetic fields and grain growth in disks around young stellar objects. However, to use polarization for probing the disk properties, its production mechanism must be ascertained first. To date, the conventional mechanism involving magnetically aligned grains fails to explain the polarization patterns detected in most disks. This is especially true for the inclined disk of HL Tau in ALMA Band 3 (wavelength $\sim 3$ mm), which has an elliptical polarization pattern. The elliptical pattern was taken as evidence for polarized emission by dust grains aligned with their long axes perpendicular the direction of the radiative flux. We show that the radiatively aligned grains produce a circular, rather than elliptical, polarization pattern even in inclined disks such as HL Tau. An elliptical polarization pattern can be produced if the grains are aligned aerodynamically by the difference in rotation speed between the dust and gas through the Gold mechanism. However, a strong azimuthal variation in polarized intensity is expected for both the radiative and aerodynamic alignment, but not observed in the HL Tau disk in ALMA Band 3. We conclude that neither of these two mechanisms alone can explain the data and the origin of the 3 mm polarization remains a mystery. We speculate that this mystery may be resolved by a combination of both direct emission and scattering by aerodynamically aligned grains.

Spatially resolved polarized millimeter/submillimeter emission has been observed in the disk of HL Tau and two other young stellar objects. It is usually interpreted as coming from magnetically aligned grains, but can also be produced by dust scattering, as demonstrated explicitly by Kataoka et al. for face-on disks. We extend their work by including the polarization induced by disk inclination with respect to the line of sight. Using a physically motivated, semi-analytic model, we show that the polarization fraction of the scattered light increases with the inclination angle $i$, reaching $1/3$ for edge-on disks. The inclination-induced polarization can easily dominate that intrinsic to the disk in the face-on view. It provides a natural explanation for the two main features of the polarization pattern observed in the tilted disk of HL Tau ($i \sim 45^\circ$): the polarized intensity concentrating in a region elongated more or less along the major axis, and polarization in this region roughly parallel to the minor axis. This broad agreement provides support to dust scattering as a viable mechanism for producing, at least in part, polarized millimeter radiation. In order to produce polarization at the observed level ($\sim 1\%$), the scattering grains must have grown to a maximum size of tens of microns. However, such grains may be too small to produce the opacity spectral index of $\beta \lesssim 1$ observed in HL Tau and other sources; another population of larger, millimeter/centimeter-sized, grains may be needed to explain the bulk of the unpolarized continuum emission.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of multiple protostar systems in the Perseus molecular cloud previously detected by the Karl G. Jansky Very Large Array (VLA). We observed 17 close ($<$600~AU separation) multiple systems at 1.3~mm in continuum and five molecular lines (i.e., \twco, \cateo, \thco, H$_2$CO, SO) to characterize the circum-multiple environments in which these systems are forming. We detect at least one component in the continuum for the 17 multiple systems. In three systems, one companion is not detected, and for two systems the companions are unresolved at our observed resolution. We also detect circum-multiple dust emission toward 8 out of 9 Class 0 multiples. Circum-multiple dust emission is not detected toward any of the 8 Class I multiples. Twelve systems are detected in the dense gas tracers toward their disks/inner envelopes. For these 12 systems, we use the dense gas observations to characterize their formation mechanism. The velocity gradients in the circum-multiple gas are clearly orthogonal to the outflow directions in 8 out of the 12 systems, consistent with disk fragmentation. Moreover, only two systems with separations $<$200~AU are \textit{inconsistent} with disk fragmentation, in addition to the two widest systems ($>$500~AU). Our results suggest that disk fragmentation via gravitational instability is an important formation mechanism for close multiple systems, but further statistics are needed to better determine the relative fraction formed via this method.

We present the full disk-fit results VANDAM survey of all Class 0 and I protostars in the Perseus molecular cloud. We have 18 new protostellar disk candidates around Class 0 and I sources, which are well described by a simple, parametrized disk model fit to the 8 mm VLA dust-continuum observations. 33% of Class 0 protostars and just 11% of Class I protostars have candidate disks, while 78% of Class 0 and I protostars do not have signs of disks within our 12 AU disk diameter resolution limit, indicating that at 8 mm most disks in the Class 0 and I phases are <10 AU in radius. These small radii may be a result of surface brightness sensitivity limits. Modeled 8 mm radii are similar to the radii of known Class 0 disks with detected Keplerian rotation. Since our 8 mm data trace a population of larger dust grains which radially drift towards the protostar and are lower limits on true disk sizes, large disks at early times do not seem to be particularly rare. We find statistical evidence that Class 0 and I disks are likely drawn from the same distribution, meaning disk properties may be defined early in the Class 0 phase and do not undergo large changes through the Class I phase. By combining our candidate disk properties with previous polarization observations, we find a qualitative indication that misalignment between inferred envelope-scale magnetic fields and outflows may indicate disks on smaller scales in Class 0 sources.

Emission from protostars at centimeter radio wavelengths has been shown to trace the free-free emission arising from ionizing shocks as a result of jets and outflows driven by protostars. Therefore, measuring properties of protostars at radio frequencies can provide valuable insights into the nature of their outflows and jets. We present a C-band (4.1 cm and 6.4 cm) survey of all known protostars (Class 0 and Class I) in Perseus as part of the VLA Nascent Disk and Multiplicity (VANDAM) Survey. We examine the known correlations between radio flux density and protostellar parameters such as bolometric luminosity and outflow force, for our sample. We also investigate the relationship between radio flux density and far-infrared line luminosities from Herschel. We show that free-free emission originates most likely from J-type shocks; however, the large scatter indicates that those two types of emission probe different time and spatial scales. Using C-band fluxes, we removed an estimation of free-free contamination from the corresponding Ka-band (9 mm) flux densities that primarily probe dust emission from embedded disks. We find that the compact ($<$~1") dust emission is lower for Class I sources (median dust mass 96 M$_{\oplus}$) relative to Class 0 (248 M$_{\oplus}$), but several times higher than in Class II (5-15 M$_{\oplus}$). If this compact dust emission is tracing primarily the embedded disk, as is likely for many sources, this result provides evidence for decreasing disk masses with protostellar evolution, with sufficient mass for forming giant planet cores primarily at early times.

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 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 870 $\mu$m ALMA observations of polarized dust emission toward the Class II protoplanetary disk IM Lup. We find that the orientation of the polarized emission is along the minor axis of the disk, and that the value of the polarization fraction increases steadily toward the center of the disk, reaching a peak value of ~1.1%. All of these characteristics are consistent with models of self-scattering of submillimeter-wave emission from an optically thin inclined disk. The distribution of the polarization position angles across the disk reveals that while the average orientation is along the minor axis, the polarization orientations show a significant spread in angles; this can also be explained by models of pure scattering. We compare the polarization with that of the Class I/II source HL Tau. A comparison of cuts of the polarization fraction across the major and minor axes of both sources reveals that IM Lup has a substantially higher polarization fraction than HL Tau toward the center of the disk. This enhanced polarization fraction could be due a number of factors, including higher optical depth in HL Tau, or scattering by larger dust grains in the more evolved IM Lup disk. However, models yield similar maximum grain sizes for both HL Tau (72 $\mu$m) and IM Lup (61 $\mu$m, this work). This reveals continued tension between grain-size estimates from scattering models and from models of the dust emission spectrum, which find that the bulk of the (unpolarized) emission in disks is most likely due to millimeter (or even centimeter) sized grains.

We present 870 $\mu$m ALMA dust polarization observations of 10 young Class 0/I protostars in the Perseus Molecular Cloud. At $\sim$ 0.35$"$ (80 au) resolution, all of our sources show some degree of polarization, with most (9/10) showing significantly extended emission in the polarized continuum. Each source has incredibly intricate polarization signatures. In particular, all three disk-candidates have polarization vectors roughly along the minor axis, which is indicative of polarization produced by dust scattering. On $\sim$ 100 au scales, the polarization is at a relatively low level ($\lesssim 1\%$) and is quite ordered. In sources with significant envelope emission, the envelope is typically polarized at a much higher ($\gtrsim 5\%$) level and has a far more disordered morphology. We compute the cumulative probability distributions for both the small (disk-scale) and large (envelope-scale) polarization percentage. We find that the two are intrinsically different, even after accounting for the different detection thresholds in the high/low surface brightness regions. We perform Kolmogorov-Smirnov and Anderson-Darling tests on the distributions of angle offsets of the polarization from the outflow axis. We find disk-candidate sources are different from the non-disk-candidate sources. We conclude that the polarization on the 100 au scale is consistent with the signature of dust scattering for disk-candidates and that the polarization on the envelope-scale in all sources may come from another mechanism, most likely magnetically aligned grains.

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.

Centimeter continuum emission from protostars offers insight into the innermost part of the outflows, as shock-ionized gas produces free-free emission. We observed a complete population of Class 0 and I protostars in the Perseus molecular cloud at 4.1 cm and 6.4 cm with resolution and sensitivity superior to previous surveys. From a total of 71 detections, 8 sources exhibit resolved emission at 4.1 cm and/or 6.4 cm. In this paper we focus on this sub-sample, analyzing their spectral indices along the jet, and their alignment with respect to the large-scale molecular outflow. Spectral indices for fluxes integrated toward the position of the protostar are consistent with free-free thermal emission. The value of the spectral index along a radio jet decreases with distance from the protostar. For six sources, emission is well aligned with the outflow central axis, showing that we observe the ionized base of the jet. This is not the case for two sources, where we note misalignment of the emission with respect to the large-scale outflow. This might indicate that the emission does not originate in the radio jet, but rather in an ionized outflow cavity wall or disk surface. For five of the sources, the spectral indices along the jet decrease well below the thermal free-free limit of -0.1 with $>2\sigma$ significance. This is indicative of synchrotron emission, meaning that high energy electrons are being produced in the outflows close to the disk. This result can have far-reaching implications for the chemical composition of the embedded disks.

We present a high angular resolution ($\sim 0.2^{\prime\prime}$), high sensitivity ($\sigma \sim 0.2$ mJy) survey of the 870 $\mu$m continuum emission from the circumstellar material around 49 pre-main sequence stars in the $\rho$ Ophiuchus molecular cloud. Because most millimeter instruments have resided in the northern hemisphere, this represents the largest high-resolution, millimeter-wave survey of the circumstellar disk content of this cloud. Our survey of 49 systems comprises 63 stars; we detect disks associated with 29 single sources, 11 binaries, 3 triple systems and 4 transition disks. We present flux and radius distributions for these systems; in particular, this is the first presentation of a reasonably complete probability distribution of disk radii at millimeter-wavelengths. We also compare the flux distribution of these protoplanetary disks with that of the disk population of the Taurus-Auriga molecular cloud. We find that disks in binaries are both significantly smaller and have much less flux than their counterparts around isolated stars. We compute truncation calculations on our binary sources and find that these disks are too small to have been affected by tidal truncation and posit some explanations for this. Lastly, our survey found 3 candidate gapped disks, one of which is a newly identified transition disk with no signature of a dip in infrared excess in extant observations.

Orion-KL is a well known high mass star forming region that has long been the target of spectral line surveys and searches for complex molecules. One spectral window where the region had never been surveyed is around wavelengths of $\lambda$=1 cm. This is an important window to observe due to the fundamental and low energy transitions of numerous complex molecules that indicate the maximum spatial extent of the molecular species; knowing the spatial distribution of a molecule aids in determining the formation mechanism(s) of that molecule. Additionally, there are fewer transitions in this window, reducing confusion caused by blended lines that can be very problematic at shorter wavelengths ($\lambda<$3 mm). In this work, we present the first spectral line survey at $\lambda$=1 cm of the Orion-KL region. A total of 89 transitions were detected from 14 molecular species and isotopologues and two atomic species. The observations were conducted with the Combined Array for Research in Millimeter-wave Astronomy in both interferometric and single dish modes.

Disk polarization at (sub)millimeter wavelengths is being revolutionized by ALMA observationally, but its origin remains uncertain. Dust scattering was recently recognized as a potential contributor to polarization, although its basic properties have yet to be thoroughly explored. Here, we quantify the effects of optical depth on the scattering-induced polarization in inclined disks through a combination of analytical illustration, approximate semi-analytical modeling using formal solution to the radiative transfer equation, and Monte Carlo simulations. We find that the near-side of the disk is significantly brighter in polarized intensity than the far-side, provided that the disk is optically thick and that the scattering grains have yet to settle to the midplane. This asymmetry is the consequence of a simple geometric effect: the near-side of the disk surface is viewed more edge-on than the far-side. It is a robust signature that may be used to distinguish the scattering-induced polarization from that by other mechanisms, such as aligned grains. The asymmetry is weaker for a geometrically thinner dust disk. As such, it opens an exciting new window on dust settling. We find anecdotal evidence from dust continuum imaging of edge-on disks that large grains are not yet settled in the youngest (Class 0) disks, but become more so in older disks. This trend is corroborated by the polarization data in inclined disks showing that younger disks have more pronounced near-far side asymmetry and thus less grain settling. If confirmed, the trend would have far-reaching implications for grain evolution and, ultimately, the formation of planetesimals and planets.

Binary and multiple star systems are a frequent outcome of the star formation process, and as a result, almost half of all sun-like stars have at least one companion star. Theoretical studies indicate that there are two main pathways that can operate concurrently to form binary/multiple star systems: large scale fragmentation of turbulent gas cores and filaments or smaller scale fragmentation of a massive protostellar disk due to gravitational instability. Observational evidence for turbulent fragmentation on scales of $>$1000~AU has recently emerged. Previous evidence for disk fragmentation was limited to inferences based on the separations of more-evolved pre-main sequence and protostellar multiple systems. The triple protostar system L1448 IRS3B is an ideal candidate to search for evidence of disk fragmentation. L1448 IRS3B is in an early phase of the star formation process, likely less than 150,000 years in age, and all protostars in the system are separated by $<$200~AU. Here we report observations of dust and molecular gas emission that reveal a disk with spiral structure surrounding the three protostars. Two protostars near the center of the disk are separated by 61 AU, and a tertiary protostar is coincident with a spiral arm in the outer disk at a 183 AU separation. The inferred mass of the central pair of protostellar objects is $\sim$1 M$_{sun}$, while the disk surrounding the three protostars has a total mass of $\sim$0.30 M$_{\sun}$. The tertiary protostar itself has a minimum mass of $\sim$0.085 M$_{sun}$. We demonstrate that the disk around L1448 IRS3B appears susceptible to disk fragmentation at radii between 150~AU and 320~AU, overlapping with the location of the tertiary protostar. This is consistent with models for a protostellar disk that has recently undergone gravitational instability, spawning one or two companion stars.

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

Observations suggest that there is a significant fraction of O-stars in the field of the Milky Way that appear to have formed in isolation or in low mass clusters ($<$100 $M_\odot$). The existence of these high-mass stars that apparently formed in the field challenges the generally accepted paradigm, which requires star formation to occur in clustered environments. In order to understand the physical conditions for the formation of these stars, it is necessary to observe isolated high-mass stars while they are still forming. With the $Hubble$ $Space$ $Telescope$, we observe the seven most isolated massive ($>$8 $M_\odot$) young stellar objects (MYSOs) in the Large Magellanic Cloud (LMC). The observations show that while these MYSOs are remote from other MYSOs, OB associations, and even from known giant molecular clouds, they are actually not isolated at all. Imaging reveals $\sim$100 to several hundred pre--main-sequence (PMS) stars in the vicinity of each MYSO. These previously undetected PMS stars form prominent compact clusters around the MYSOs, and in most cases they are also distributed sparsely across the observed regions. Contrary to what previous high-mass field star studies show, these observations suggest that high-mass stars may not be able to form in clusters with masses less than 100 $M_\odot$. If these MYSOs are indeed the best candidates for isolated high-mass star formation, then the lack of isolation is at odds with random sampling of the IMF. Moreover, while isolated MYSOs may not exist, we find evidence that isolated clusters containing O-stars can exist, which in itself is rare.

There is substantial evidence for disk formation taking place during the early stages of star formation and for most stars being born in multiple systems; however, protostellar multiplicity and disk searches have been hampered by low resolution, sample bias, and variable sensitivity. We have conducted an unbiased, high-sensitivity Karl G. Jansky Very Large Array (VLA) survey toward all known protostars (n = 94) in the Perseus molecular cloud (d~230 pc), with a resolution of ~15 AU (0.06") at 8 mm. We have detected candidate protostellar disks toward 17 sources (with 12 of those in the Class 0 stage) and we have found substructure on < 50AU scales for three Class 0 disk candidates, possibly evidence for disk fragmentation. We have discovered 16 new multiple systems (or new components) in this survey; the new systems have separations < 500 AU and 3 by < 30 AU. We also found a bi-modal distribution of separations, with peaks at ~75 AU and ~3000 AU, suggestive of formation through two distinct mechanisms: disk and turbulent fragmentation. The results from this survey demonstrate the necessity and utility of uniform, unbiased surveys of protostellar systems at millimeter and centimeter wavelengths.

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.

Dust polarization in millimeter (and centimeter) has been mapped in disks around an increasing number of young stellar objects. It is usually thought to come from emission by magnetically aligned (non-spherical) grains, but can also be produced by dust scattering. We present a semi-analytic theory of disk polarization that includes both the direction emission and scattering, with an emphasis on their relative importance and how they are affected by the disk inclination. For face-on disks, both emission and scattering tend to produce polarization in the radial direction, making them difficult to distinguish, although the scattering-induced polarization can switch to the azimuthal direction if the incident radiation is beamed strongly enough in the radial direction in the disk plane. Disk inclination affects the polarizations from emission and scattering differently, especially on the major axis where, in the edge-on limit, the former vanishes while the latter reaches a polarization fraction as large as 1/3. The polarizations from the two competing mechanisms tend to cancel each other on the major axis, producing two low polarization "holes" (one on each side of the center) under certain conditions. We find tantalizing evidence for at least one such "hole" in NGC1333 IRAS4A1, whose polarization observed at 8 mm on the 100 AU scale is indicative of a pattern dominated by scattering close to the center and by direction emission in the outer region. If true, it would imply not only that a magnetic field exists on the disk scale, but that it is strong enough to align large, possibly mm-sized, grains.

We present a multiplicity study of all known protostars (94) in the Perseus molecular cloud from a Karl G. Jansky Very Large Array (VLA) survey at Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.6 cm). The observed sample has a bolometric luminosity range between 0.1 L$_{\odot}$ and $\sim$33 L$_{\odot}$, with a median of 0.7 L$_{\odot}$. This multiplicity study is based on the Ka-band data, having a best resolution of $\sim$0.065" (15 AU) and separations out to $\sim$43" (10000 AU) can be probed. The overall multiplicity fraction (MF) is found to be of 0.40$\pm$0.06 and the companion star fraction (CSF) is 0.71$\pm$0.06. The MF and CSF of the Class 0 protostars are 0.57$\pm$0.09 and 1.2$\pm$0.2, and the MF and CSF of Class I protostars are both 0.23$\pm$0.08. The distribution of companion separations appears bi-modal, with a peak at $\sim$75 AU and another peak at $\sim$3000 AU. Turbulent fragmentation is likely the dominant mechanism on $>$1000 AU scales and disk fragmentation is likely to be the dominant mechanism on $<$200 AU scales. Toward three Class 0 sources we find companions separated by $<$30 AU. These systems have the smallest separations of currently known Class 0 protostellar binary systems. Moreover, these close systems are embedded within larger (50 AU to 400 AU) structures and may be candidates for ongoing disk fragmentation.

We present the first dust emission results toward a sample of seven protostellar disk candidates around Class 0 and I sources in the Perseus molecular cloud from the VLA Nascent Disk and Multiplicity (VANDAM) survey with ~0.05'' or 12 AU resolution. To examine the surface brightness profiles of these sources, we fit the Ka-band 8 mm dust-continuum data in the u,v-plane to a simple, parametrized model based on the Shakura-Sunyaev disk model. The candidate disks are well-fit by a model with a disk-shaped profile and have masses consistent with known Class 0 and I disks. The inner-disk surface densities of the VANDAM candidate disks have shallower density profiles compared to disks around more evolved Class II systems. The best-fit model radii of the seven early-result candidate disks are R_c > 10 AU; at 8 mm, the radii reflect lower limits on the disk size since dust continuum emission is tied to grain size and large grains radially drift inwards. These relatively large disks, if confirmed kinematically, are inconsistent with theoretical models where the disk size is limited by strong magnetic braking to < 10 AU at early times.

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

Magnetic fields can regulate disk formation, accretion and jet launching. Until recently, it has been difficult to obtain high resolution observations of the magnetic fields of the youngest protostars in the critical region near the protostar. The VANDAM survey is observing all known protostars in the Perseus Molecular Cloud. Here we present the polarization data of IRAS 4A. We find that with ~ 0.2'' (50 AU) resolution at {\lambda} = 8.1 and 10.3 mm, the inferred magnetic field is consistent with a circular morphology, in marked contrast with the hourglass morphology seen on larger scales. This morphology is consistent with frozen-in field lines that were dragged in by rotating material entering the infall region. The field morphology is reminiscent of rotating circumstellar material near the protostar. This is the first polarization detection of a protostar at these wavelengths. We conclude from our observations that the dust emission is optically thin with {\beta} ~ 1.3, suggesting that mm/cm-sized grains have grown and survived in the short lifetime of the protostar.

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 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 are conducting a Jansky VLA Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.4 cm) survey of all known protostars in the Perseus Molecular Cloud, providing resolution down to $\sim$0.06'' and $\sim$0.35" in Ka-band and C-band, respectively. Here we present first results from this survey that enable us to examine the source NGC 1333 IRAS2A in unprecedented detail and resolve it into a proto-binary system separated by 0.621"$\pm$0.006" ($\sim$143 AU) at 8 mm, 1 cm, and 4 cm. These 2 sources (IRAS2A VLA1 and VLA2) are likely driving the two orthogonal outflows known to originate from IRAS2A. The brighter source IRAS2A VLA1 is extended perpendicular to its outflow in the VLA data, with a deconvolved size of 0.055" ($\sim$13 AU), possibly tracing a protostellar disk. The recently reported candidate companions (IRAS2A MM2 and MM3) are not detected in either our VLA data, CARMA 1.3 mm data, or SMA 850 $\mu$m data. SMA CO ($J=3\rightarrow2$), CARMA CO ($J=2\rightarrow1$), and lower resolution CARMA CO ($J=1\rightarrow0$) observations are used to examine the outflow origins and the nature of the candidate companions to IRAS2A VLA1. The CO ($J=3\rightarrow2$) and ($J=2\rightarrow1$) data show that IRAS2A MM2 is coincident with a bright CO emission spot in the east-west outflow, and IRAS2A MM3 is within the north-south outflow. In contrast, IRAS2A VLA2 lies at the east-west outflow symmetry point. We propose that IRAS2A VLA2 is the driving source of the East-West outflow and a true companion to IRAS2A VLA1, whereas IRAS2A MM2 and MM3 may not be protostellar.

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 parsec scale interferometric maps of HCN(1-0) and HCO$^{+}$(1-0) emission from dense gas in the star-forming region 30 Doradus, obtained using the Australia Telescope Compact Array. This extreme star-forming region, located in the Large Magellanic Cloud (LMC), is characterized by a very intense ultraviolet ionizing radiation field and sub-solar metallicity, both of which are expected to impact molecular cloud structure. We detect 13 bright, dense clumps within the 30 Doradus-10 giant molecular cloud. Some of the clumps are aligned along a filamentary structure with a characteristic spacing that is consistent with formation via the varicose fluid instability. Our analysis shows that the filament is gravitationally unstable and collapsing to form stars. There is a good correlation between HCO$^{+}$ emission in the filament and signatures of recent star formation activity including H$_{2}$O masers and young stellar objects (YSOs). YSOs seem to continue along the same direction of the filament toward the massive compact star cluster R136 in the southwest. We present detailed comparisons of clump properities (masses, linewidths, sizes) in 30Dor-10 to those in other star forming regions of the LMC (N159, N113, N105, N44). Our analysis shows that the 30Dor-10 clumps have similar mass but wider linewidths and similar HCN/HCO$^{+}$ (1-0) line ratios as clumps detected in other LMC star-forming regions. Our results suggest that the dense molecular gas clumps in the interior of 30Dor-10 are well-shielded against the intense ionizing field that is present in the 30Doradus region.

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.

In this paper, we present the results of an observational search for gas phase urea [(NH$_2$)$_2$CO] observed towards the Sgr B2(N-LMH) region. We show data covering urea transitions from $\sim$100 GHz to 250 GHz from five different observational facilities: BIMA, CARMA, the NRAO 12 m telescope, the IRAM 30 m telescope, and SEST. The results show that the features ascribed to urea can be reproduced across the entire observed bandwidth and all facilities by best fit column density, temperature, and source size parameters which vary by less than a factor of 2 between observations merely by adjusting for telescope-specific parameters. Interferometric observations show that the emission arising from these transitions is cospatial and compact, consistent with the derived source sizes and emission from a single species. Despite this evidence, the spectral complexity, both of (NH$_2$)$_2$CO and of Sgr B2(N), makes the definitive identification of this molecule challenging. We present observational spectra, laboratory data, and models, and discuss our results in the context of a possible molecular detection of urea.

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 study the structure and kinematics of nine 0.1 pc-scale cores in Orion with the IRAM 30-m telescope and at higher resolution eight of the cores with CARMA, using CS(2-1) as the main tracer. The single-dish moment zero maps of the starless cores show single structures with central column densities ranging from 7 to 42 times 10^23 cm^-2 and LTE masses from 20 solar masses to 154 solar masses. However, at the higher CARMA resolution (5 arcsec), all of the cores except one fragment into 3 - 5 components. The number of fragments is small compared to that found in some turbulent fragmentation models, although inclusion of magnetic fields may reduce the predicted fragment number and improve the model agreement. This result demonstrates that fragmentation from parsec-scale molecular clouds to sub-parsec cores continues to take place inside the starless cores. The starless cores and their fragments are embedded in larger filamentary structures, which likely played a role in the core formation and fragmentation. Most cores show clear velocity gradients, with magnitudes ranging from 1.7 to 14.3 km/s/pc. We modeled one of them in detail, and found that its spectra are best explained by a converging flow along a filament toward the core center; the gradients in other cores may be modeled similarly. We infer a mass inflow rate of ~ 2 x 10^{-3} Msolar/yr, which is in principle high enough to overcome radiation pressure and allow for massive star formation. However, the core contains multiple fragments, and it is unclear whether the rapid inflow would feed the growth of primarily a single massive star or a cluster of lower mass objects. We conclude that fast, supersonic converging flow along filaments play an important role in massive star and cluster formation.

We present high-resolution sub/millimeter interferometric imaging of the Class 0 protostar L1527 IRS (IRAS 04368+2557) at 870 micron and 3.4 mm from the Submillimeter Array (SMA) and Combined Array for Research in Millimeter Astronomy (CARMA). We detect the signature of an edge-on disk surrounding the protostar with an observed diameter of 180 AU in the sub/millimeter images. The mass of the disk is estimated to be 0.007 M_sun, assuming optically thin, isothermal dust emission. The millimeter spectral index is observed to be quite shallow at all the spatial scales probed; alpha ~ 2, implying a dust opacity spectral index beta ~ 0. We model the emission from the disk and surrounding envelope using Monte Carlo radiative transfer codes, simultaneously fitting the sub/millimeter visibility amplitudes, sub/millimeter images, resolved L\arcmin\ image, spectral energy distribution, and mid-infrared spectrum. The best fitting model has a disk radius of R = 125 AU, is highly flared (H ~ R^1.3), has a radial density profile rho ~ R^-2.5, and has a mass of 0.0075 M_sun. The scale height at 100 AU is 48 AU, about a factor of two greater than vertical hydrostatic equilibrium. The resolved millimeter observations indicate that disks may grow rapidly throughout the Class 0 phase. The mass and radius of the young disk around L1527 is comparable to disks around pre-main sequence stars; however, the disk is considerably more vertically extended, possibly due to a combination of lower protostellar mass, infall onto the disk upper layers, and little settling of ~1 micron-sized dust grains.

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 a study on the spatial distribution of N2D+ and N2H+ in thirteen protostellar systems. Eight of thirteen objects observed with the IRAM 30m telescope show relative offsets between the peak N2D+ (J=2-1) and N2H+ (J=1-0) emission. We highlight the case of L1157 using interferometric observations from the Submillimeter Array and Plateau de Bure Interferometer of the N2D+ (J=3-2) and N2H+ (J=1-0) transitions respectively. Depletion of N2D+ in L1157 is clearly observed inside a radius of ~2000 AU (7") and the N2H+ emission is resolved into two peaks at radii of ~1000 AU (3.5"), inside the depletion region of N2D+. Chemical models predict a depletion zone in N2H+ and N2D+ due to destruction of H2D+ at T ~ 20 K and the evaporation of CO off dust grains at the same temperature. However, the abundance offsets of 1000 AU between the two species are not reproduced by chemical models, including a model that follows the infall of the protostellar envelope. The average abundance ratios of N2D+ to N2H+ have been shown to decrease as protostars evolve by Emprechtinger et al., but this is the first time depletion zones of N2D+ have been spatially resolved. We suggest that the difference in depletion zone radii for N2H+ and N2D+ is caused by either the CO evaporation temperature being above 20 K or an H2 ortho-to-para ratio gradient in the inner envelope.

In their earliest stages, protostars accrete mass from their surrounding envelopes through circumstellar disks. Until now, the smallest observed protostar/envelope mass ratio was ~2.1. The protostar L1527 IRS is thought to be in the earliest stages of star formation. Its envelope contains ~1 solar mass of material within a ~0.05 pc radius, and earlier observations suggested the presence of an edge-on disk. Here we report observations of dust continuum emission and 13CO (J=2-1) line emission from the disk around L1527, from which we determine a protostellar mass of M = 0.19 +/- 0.04 solar masses and a protostar/envelope mass ratio of ~0.2. We conclude that most of the luminosity is generated through the accretion process, with an accretion rate of ~6.6 x 10^-7 solar masses per year. If it has been accreting at that rate through much of its life, its age is ~300,000 yr, though theory suggests larger accretion rates earlier, so it may be younger. The presence of a rotationally--supported disk is confirmed and significantly more mass may be added to its planet-forming region as well as the protostar itself.

Filamentary structures are ubiquitous from large-scale molecular clouds (few parsecs) to small-scale circumstellar envelopes around Class 0 sources (~1000 AU to ~0.1 pc). In particular, recent observations with the Herschel Space Observatory emphasize the importance of large-scale filaments (few parsecs) and star formation. The small-scale flattened envelopes around Class 0 sources are reminiscent of the large-scale filaments. We propose an observationally derived scenario for filamentary star formation that describes the evolution of filaments as part of the process for formation of cores and circumstellar envelopes. If such a scenario is correct, small-scale filamentary structures (0.1 pc in length) with higher densities embedded in starless cores should exist, although to date almost all the interferometers have failed to observe such structures. We perform synthetic observations of filaments at the prestellar stage by modeling the known Class 0 flattened envelope in L1157 using both the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Atacama Large Millimeter/Submillimeter Array (ALMA). We show that with reasonable estimates for the column density through the flattened envelope, the CARMA D-array at 3mm wavelengths is not able to detect such filamentary structure, so previous studies would not have detected them. However, the substructures may be detected with CARMA D+E array at 3 mm and CARMA E array at 1 mm as a result of more appropriate resolution and sensitivity. ALMA is also capable of detecting the substructures and showing the structures in detail compared to the CARMA results with its unprecedented sensitivity. Such detection will confirm the new proposed paradigm of non-spherical star formation.

We present dual-wavelength observations and modeling of the nearly edge-on Class 0 young stellar object L1157-mm. Using the Combined Array for Research in Millimeter-wave Astronomy, a nearly spherical structure is seen from the circumstellar envelope at the size scale of 10^2 to 10^3 AU in both 1 mm and 3 mm dust emission. Radiative transfer modeling is performed to compare data with theoretical envelope models, including a power-law envelope model and the Terebey-Shu-Cassen model. Bayesian inference is applied for parameter estimation and information criteria is used for model selection. The results prefer the power-law envelope model against the Terebey-Shu-Cassen model. In particular, for the power-law envelope model, a steep density profile with an index of ~2 is inferred. Moreover, the dust opacity spectral index (beta) is estimated to be ~0.9, implying that grain growth has started at L1157-mm. Also, the unresolved disk component is constrained to be < 40 AU in radius and < 4-25 M_Jup in mass. However, the estimate of the embedded disk component relies on the assumed envelope model.

We report the results of a high spatial (parsec) resolution HCO+ (J = 1-0) and HCN (J = 1-0) emission survey toward the giant molecular clouds of the star formation regions N105, N113, N159, and N44 in the Large Magellanic Cloud. The HCO+ and HCN observations at 89.2 and 88.6 GHz, respectively, were conducted in the compact configuration of the Australia Telescope Compact Array. The emission is imaged into individual clumps with masses between 10^2 and 10^4 solar masses and radii of <1 pc to ~2 pc. Many of the clumps are coincident with indicators of current massive star formation, indicating that many of the clumps are associated with deeply-embedded forming stars and star clusters. We find that massive YSO-bearing clumps tend to be larger (>1 pc), more massive (M > 10^3 solar masses), and have higher surface densities (~1 g cm^-2), while clumps without signs of star formation are smaller (<1 pc), less massive (M < 10^3 solar masses), and have lower surface densities (~0.1 g cm^-2). The dearth of massive (M >10^3 solar masses) clumps not bearing massive YSOs suggests the onset of star formation occurs rapidly once the clump has attained physical properties favorable to massive star formation. Using a large sample of LMC massive YSO mid-IR spectra, we estimate that ~2/3 of the massive YSOs for which there are Spitzer mid-IR spectra are no longer located in molecular clumps; we estimate that these young stars/clusters have destroyed their natal clumps on a time scale of at least 3 x 10^{5}$ yrs.

We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources with well-ordered velocity fields and envelope structures are now analyzed in more detail. We compare the interferometric data to position-velocity diagrams of kinematic models for spherical rotating collapse and filamentary rotating collapse. For this purpose, we developed a filamentary parametrization of the rotating collapse model to explore the effects of geometric projection on the observed velocity structures. We find that most envelopes in our sample have PV structures that can be reproduced by an infalling filamentary envelope projected at different angles within the plane of the sky. The infalling filament produces velocity shifts across the envelope that can mimic rotation, especially when viewed at single-dish resolutions and the axisymmetric rotating collapse model does not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect rotation, then the inferred centrifugal radii are quite large in most cases, indicating significant fragmentation potential or more likely another component to the line-center velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation alone when envelopes are non-axisymmetric and that projected infall velocities likely dominate the velocity field on scales larger than 1000 AU.

We present the properties of an extensive sample of molecular clouds in the Large Magellanic Cloud (LMC) mapped at 11 pc resolution in the CO(1-0) line. We identify clouds as regions of connected CO emission, and find that the distributions of cloud sizes, fluxes and masses are sensitive to the choice of decomposition parameters. In all cases, however, the luminosity function of CO clouds is steeper than dN/dL \propto L^{-2}, suggesting that a substantial fraction of mass is in low-mass clouds. A correlation between size and linewidth, while apparent for the largest emission structures, breaks down when those structures are decomposed into smaller structures. We argue that the correlation between virial mass and CO luminosity is the result of comparing two covariant quantities, with the correlation appearing tighter on larger scales where a size-linewidth relation holds. The virial parameter (the ratio of a cloud's kinetic to self-gravitational energy) shows a wide range of values and exhibits no clear trends with the CO luminosity or the likelihood of hosting young stellar object (YSO) candidates, casting further doubt on the assumption of virialization for molecular clouds in the LMC. Higher CO luminosity increases the likelihood of a cloud harboring a YSO candidate, and more luminous YSOs are more likely to be coincident with detectable CO emission, confirming the close link between giant molecular clouds and massive star formation.

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 a study of dense molecular gas kinematics in seventeen nearby protostellar systems using single-dish and interferometric molecular line observations. The non-axisymmetric envelopes around a sample of Class 0/I protostars were mapped in the N2H+ (J=1-0) tracer with the IRAM 30m, CARMA and PdBI as well as NH3 (1,1) with the VLA. The molecular line emission is used to construct line-center velocity and linewidth maps for all sources to examine the kinematic structure in the envelopes on spatial scales from 0.1 pc to ~1000 AU. The direction of the large-scale velocity gradients from single-dish mapping is within 45 degrees of normal to the outflow axis in more than half the sample. Furthermore, the velocity gradients are often quite substantial, the average being ~2.3 km\s\pc. The interferometric data often reveal small-scale velocity structure, departing from the more gradual large-scale velocity gradients. In some cases, this likely indicates accelerating infall and/or rotational spin-up in the inner envelope; the median velocity gradient from the interferometric data is ~10.7 km/s/pc. In two systems, we detect high-velocity HCO+ (J=1-0) emission inside the highest-velocity \nthp\ emission. This enables us to study the infall and rotation close to the disk and estimate the central object masses. The velocity fields observed on large and small-scales are more complex than would be expected from rotation alone, suggesting that complex envelope structure enables other dynamical processes (i.e. infall) to affect the velocity field.

We present observations of the intermediate to massive star-forming region I05345+3157 using the molecular line tracer CS(2-1) with CARMA to reveal the properties of the dense gas cores. Seven gas cores are identified in the integrated intensity map of CS(2-1). Among these, core 1 and core 3 have counterparts in the 2.7 millimeter continuum data. We suggest that core 1 and core 3 are star-forming cores that may already or will very soon harbor young massive protostars. The total masses of core 1 estimated from the LTE method and dust emission by assuming a gas-to-dust ratio are 5 +- 1 solar masses and 18 +- 6 solar masses, and that of core 3 are 15 +- 7 solar masses and 11 +- 3 solar masses. The spectrum of core 3 shows blue-skewed self-absorption, which suggests gas infall -- a collapsing core. The observed broad linewidths of the seven gas cores indicate non-thermal motions. These non-thermal motions can be interactions with nearby outflows or due to the initial turbulence; the former is observed, while the role of initial turbulence is less certain. Finally, the virial masses of the gas cores are larger than the LTE masses, which for a bound core implies a requirement on the external pressure of ~ 10^8 K/cm^3. The cores have the potential to further form massive stars.

We investigate the effect of the Milky Way's magnetic field in star forming regions using archived 350 micron polarization data on 52 Galactic star formation regions from the Hertz polarimeter module. The polarization angles and percentages for individual telescope beams were combined in order to produce a large-scale average for each source and for complexes of sources. In more than 80% of the sources, we find a meaningful mean magnetic field direction, implying the existence of an ordered magnetic field component at the scale of these sources. The average polarization angles were analyzed with respect to the Galactic coordinates in order to test for correlations between polarization percentage, polarization angle, intensity, and Galactic location. No correlation was found, which suggests that the magnetic field in dense molecular clouds is decoupled from the large-scale Galactic magnetic field. Finally, we show that the magnetic field directions in the complexes are consistent with a random distribution on the sky.

We present an analytical study of a large sample of ~109 young stellar objects in the X-ray. Our objects were detected in X-ray independent of age. Unexpectedly, the X-ray energy is somewhat correlated with the ages. It decreases with time and with column density, while it should increase. We conclude that the youngest protostars, Class 0/I, emit X-rays in the 1-8 keV band. The deeply embedded sources with the strongest accretion activity are detected in the hard-band (> 2keV) only. Due to extinction, their soft X-rays are not detected. To explain the decline in energy, we suggest that within a timescale of few Myrs the corona cools down via the accretion material.

Extinction maps at 8 micron from the Spitzer Space Telescope show that many Class 0 protostars exhibit complex, irregular, and on-axisymmetric structure within the densest regions of their dusty envelopes. Many of the systems have highly irregular and on-axisymmetric morphologies on scales $\sim$1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. We suggest that the observed envelope complexity is the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We then show that the kinematics around L1165 as probed with N2H+ are indicative of asymmetric infall; the velocity gradient is not perpendicular to the outflow.

It has often been suggested that supernova remnants (SNRs) can trigger star formation. To investigate the relationship between SNRs and star formation, we have examined the known sample of 45 SNRs in the Large Magellanic Cloud to search for associated young stellar objects (YSOs) and molecular clouds. We find seven SNRs associated with both YSOs and molecular clouds, three SNRs associated with YSOs but not molecular clouds, and eight SNRs near molecular clouds but not associated with YSOs. Among the 10 SNRs associated with YSOs, the association between the YSOs and SNRs can be either rejected or cannot be convincingly established for eight cases. Only two SNRs have YSOs closely aligned along their rims; however, the time elapsed since the SNR began to interact with the YSOs' natal clouds is much shorter than the contraction timescales of the YSOs, and thus we do not see any evidence of SNR-triggered star formation in the LMC. The 15 SNRs that are near molecular clouds may trigger star formation in the future when the SNR shocks have slowed down to <45 km/s. We discuss how SNRs can alter the physical properties and abundances of YSOs.

Acetic acid (CH$_3$COOH) has been detected mainly in hot molecular cores where the distribution between oxygen (O) and nitrogen (N) containing molecular species is co-spatial within the telescope beam. Previous work has presumed that similar cores with co-spatial O and N species may be an indicator for detecting acetic acid. However, does this presumption hold as higher spatial resolution observations become available of large O and N-containing molecules? As the number of detected acetic acid sources is still low, more observations are needed to support this postulate. In this paper, we report the first acetic acid survey conducted with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at 3 mm wavelengths towards G19.61-0.23, G29.96-0.02 and IRAS 16293-2422. We have successfully detected CH$_3$COOH via two transitions toward G19.61-0.23 and tentatively confirmed the detection toward IRAS 16293-2422 A. The determined column density of CH$_3$COOH is 2.0(1.0)$\times 10^{16}$ cm$^{-2}$ and the abundance ratio of CH$_3$COOH to methyl formate (HCOOCH$_3$) is 2.2(0.1)$\times 10^{-1}$ toward G19.61-0.23. Toward IRAS 16293 A, the determined column density of CH$_3$COOH is $\sim$ 1.6 $\times 10^{15}$ cm$^{-2}$ and the abundance ratio of CH$_3$COOH to methyl formate (HCOOCH$_3$) is $\sim$ 1.0 $\times 10^{-1}$ both of which are consistent with abundance ratios determined toward other hot cores. Finally, we model all known line emission in our passband to determine physical conditions in the regions and introduce a new metric to better reveal weak spectral features that are blended with stronger lines or that may be near the 1-2$\sigma$ detection limit.

We use archived IRAC images from the Spitzer Space Telescope to show that many Class 0 protostars exhibit complex, irregular, and non-axisymmetric structure within their dusty envelopes. Our 8 $\mu$m extinction maps probe some of the densest regions in these protostellar envelopes. Many of the systems are observed to have highly irregular and non-axisymmetric morphologies on scales >= 1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. These results indicate that mass ejection is not responsible for much of the irregular morphologies we detect; rather, we suggest that the observed envelope complexity is mostly the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We also note that protostars in the sample appear to be formed preferentially near the edges of clouds or bends in filaments, suggesting formation by gravitational focusing.

We present observations of the Class 0 protostar L1157-mm using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) in 3 mm dust continuum and N2H+ line emission. In the N2H+ line, we detect a large-scale envelope extended over a linear size of ~20,000AU flattened in the direction perpendicular to the outflow. This N2H+ feature coincides with the outer envelope seen in the 8 micron extinction by Looney et al. Meanwhile, the dust continuum traces the compact, nearly spherical structure of the inner envelope, where N2H+ becomes depleted. This highly flattened N2H+ envelope also shows dynamical signatures consistent with gravitational infall in the inner region, but a slow, solid-body rotation at large scales. This flattened structure is not a rotationally supported circumstellar disk; instead, it resembles a prestellar core both morphologically and kinematically, representing the early phase of a Class 0 system. In this paper, we construct a simple model to interpret both the dust continuum and N2H+ emission and suggest a possible dynamical scenario for the overall properties of the envelope.

We have used archival HST H$\alpha$ images to study the immediate environments of massive and intermediate-mass young stellar object (YSO) candidates in the Large Magellanic Cloud (LMC). The sample of YSO candidates, taken from Gruendl & Chu (2009), was selected based on Spitzer IRAC and MIPS observations of the entire LMC and complementary ground-based optical and near-infrared observations. We found HST H$\alpha$ images for 99 YSO candidates in the LMC, of which 82 appear to be genuine YSOs. More than 95% of the YSOs are found to be associated with molecular clouds. YSOs are seen in three different kinds of environments in the H$\alpha$ images: in dark clouds, inside or on the tip of bright-rimmed dust pillars, and in small H II regions. Comparisons of spectral energy distributions for YSOs in these three different kinds of environments suggest that YSOs in dark clouds are the youngest, YSOs with small H II regions are the most evolved, and YSOs in bright-rimmed dust pillars span a range of intermediate evolutionary stages. This rough evolutionary sequence is substantiated by the presence of silicate absorption features in the Spitzer IRS spectra of some YSOs in dark clouds and in bright-rimmed dust pillars, but not those of YSOs in small H II regions. We present a discussion on triggered star formation for YSOs in bright-rimmed dust pillars or in dark clouds adjacent to H II regions. As many as 50% of the YSOs are resolved into multiple sources in high-resolution HST images. This illustrates the importance of using high-resolution images to probe the true nature and physical properties of YSOs in the LMC.

We present and categorize Spitzer IRS spectra of 294 objects in the Large Magellanic Cloud (LMC) to create the largest and most complete catalog of massive young stellar object (YSO) spectra in the LMC. Target sources were identified from infrared photometry and multi-wavelength images indicative of young, massive stars highly enshrouded in their natal gas and dust clouds. Our sample primarily consists of 277 objects we identify as having spectral features indicative of embedded YSOs. The remaining sources are comprised of 7 C-rich evolved sources, 8 sources dominated by broad silicate emission, and 1 source with multiple broad emission features. Those with YSO-like spectra show a range of spectral features including polycyclic aromatic hydrocarbon emission, deep silicate absorption, fine-structure lines, and ice absorption features. Based upon the relative strengths of these features, we have classified the YSO candidates into several distinct categories using the widely-used statistical procedure known as principal component analysis. We propose that these categories represent a spectrum of evolutionary stages during massive YSO formation. We conclude that massive pre-main sequence stars spend a majority of their massive, embedded lives emitting in the UV. Half of the sources in our study have features typical of compact HII regions, suggesting that massive YSOs can create a detectable compact HII region half-way through the formation time present in our sample. This study also provides a check on commonly used source-selection procedures including the use of photometry to identify YSOs. We determine a high success rate (>95%) of identifying objects with YSO-like spectra can be achieved through careful use of infrared CMDs, SEDs, and image inspections.

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.

Using Spitzer IRAC and MIPS observations of the Large Magellanic Cloud, we have identified 13 objects that have extremely red mid-IR colors. Follow-up Spitzer IRS observations of seven of these sources reveal varying amounts of SiC and C2H2 absorption as well as the presence of a broad MgS feature in at least two cases, indicating that these are extreme carbon stars. Preliminary estimates find these objects have luminosities of 4-11x10^3 Lsol and preliminary model fitting gives mass-loss rates between 4x10^-5 and 2x10^-4 Msol/yr, higher than any known carbon-rich AGB star in the LMC. These spectral and physical properties require careful reconsideration of dust condensation and mass-loss processes for carbon stars in low metallicity environments.

We have observed a small sample of powerful double radio sources (radio galaxies and quasars) at frequencies around 90 GHz with the BIMA millimetre array, with the intention of constraining the resolved high-frequency spectra of radio galaxies. When combined with other sources we have previously observed and with data from the BIMA archive, these observations allow us for the first time to make general statements about the high-frequency behaviour of compact components of radio galaxies -- cores, jets and hotspots. We find that cores in our sample remain flat-spectrum up to 90 GHz; jets in some of our targets are detected at 90 GHz for the first time in our new observations; and hotspots are found to be almost universal, but show a wide range of spectral properties. Emission from the extended lobes of radio galaxies is detected in a few cases and shows rough consistency with the expectations from standard spectral ageing models, though our ability to probe this in detail is limited by the sensitivity of BIMA. We briefly discuss the prospects for radio-galaxy astrophysics with ALMA.

Using interferometric data from BIMA observations, combined with detailed modeling in Fourier space of the physical structures predicted by models, we constrain the circumstellar envelope parameters for four Class 0 young stellar objects, as well as their embedded circumstellar disks. The envelopes of these objects are still undergoing collapse, and theoretical collapse models can be compared to the observations. Since it has been suggested in a previous study that both the Larson-Penston and Shu similarity solutions underestimate the age of the system, we adopt Tassis & Mouschovias' model of the collapse process, which includes all relevant magnetic fields effects. The results of the model fitting show a good consistency between theory and data; furthermore, no age problem exists since the Tassis & Mouschovias' model is age independent for the first 255 kyr. Although the majority of the continuum dust emission arises from the circumstellar envelopes, these objects have well known outflows, which suggest the presence of circumstellar disks. At the highest resolution, most of the large-scale envelope emission is resolved out by interferometry, but the small-scale residual emission remains, making it difficult to observe only the compact disk component. By modeling the emission of the envelope and subtracting it from the total emission, we constrain the disk masses in our four systems to be comparable to or smaller than the typical disk masses for T Tauri systems.

We present a statistical study of the X-Ray emission toward 22 Herbig AeBe stars using the Chandra archive. We probe the origin of the X-Rays toward Herbig stars: are they intrinsic? This question is addressed by correlations between the physical stellar properties and the X-Ray emission. There is a weak correlation between the continuum radio emission at lambda = 3.6cm and Lx, which suggests that the X-Ray emission depends upon the source. On the other hand, no correlation was found with the stellar rotational period, but that only excludes solar-like magnetic activity as the origin of the X-Rays. Most importantly, the X-Ray luminosity of Herbig AeBe stars have a different distribution than T Tauri stars, suggesting X-Ray emission from an unseen late type star companion can be ruled out with an 80% confidence level. This implies that the Herbig AeBe stars must have magnetic activity. In addition, we report the observation of five sources for the first time, three detections.

We present Spitzer IRAC images that indicate the presence of cavities cut into the dense outer envelope surrounding very young pre-main sequence stars. These young stellar objects (YSOs) characterized by an outflow represent the earliest stages of star formation. Mid-infrared photons thermally created by the central protostar/disk are scattered by dust particles within the outflow cavity itself into the line of sight. We observed this scattered light from 27 nearby, cavity-resolved YSOs, and quantified the shape of the outflow cavities. Using the grid models of Robitaille et al. (2006), we matched model spectral energy distributions (SEDs) to the observed SEDs of the 27 cataloged YSOs using photometry from IRAC, MIPS, and IRAS. This allows for the estimation of geometric and physical properties such as inclination angle, cavity density, and accretion rate. By using the relative parameter estimates determined by the models, we are able to deduce an evolutionary picture for outflows. Our work supports the concept that cavities widen with time, beginning as a thin jet-like outflow that widens to reveal the central protostar and disk until the protostellar envelope is completely dispersed by outflow and accretion.

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 CS $J = 2 \to 1$, $^{13}$CO $J = 1 \to 0$, and C$^{18}$O $J = 1 \to 0$, observations with the 10-element Berkeley Illinois Maryland Association (BIMA) Array toward the young cluster around the Be star MWC 1080. These observations reveal a biconical outflow cavity with size $\sim$ 0.3 and 0.05 pc for the semimajor and semiminor axis and $\sim$ 45$\arcdeg$ position angle. These transitions trace the dense gas, which is likely the swept-up gas of the outflow cavity, rather than the remaining natal gas or the outflow gas. The gas is clumpy; thirty-two clumps are identified. The identified clumps are approximately gravitationally bound and consistent with a standard isothermal sphere density, which suggests that they are likely collapsing protostellar cores. The gas kinematics suggests that there exists velocity gradients implying effects from the inclination of the cavity and MWC 1080. The kinematics of dense gas has also been affected by either outflows or stellar winds from MWC 1080, and lower-mass clumps are possibly under stronger effects from MWC 1080 than higher-mass clumps. In addition, low-mass cluster members tend to be formed in the denser and more turbulent cores, compared to isolated low-mass star-forming cores. This results from contributions of nearby forming massive stars, such as outflows or stellar winds. Therefore, we conclude that in clusters like the MWC 1080 system, effects from massive stars dominate the star-forming environment in both the kinematics and dynamics of the natal cloud and the formation of low-mass cluster members. This study provides insights into the effects of MWC 1080 on its natal cloud, and suggests a different low-mass star forming environment in clusters compared to isolated star formation.

(abridged) We present high angular resolution (~ 1") multi-tracer spectral line observations toward the hot core associated with G34.26+0.15 between 87--109 GHz. We have mapped emission from (i) complex nitrogen- and oxygen-rich molecules like CH3OH, HC3N, C2H5CN, NH2CHO, CH3OCH3, HCOOCH3; (ii) sulfur-bearing molecules like OCS, SO and SO2; and (iii) the recombination line H53 beta. The high angular resolution enables us to directly probe the hot molecular core associated with G34.26+0.15 at spatial scales of 0.018 pc. At this resolution we find no evidence for the hot core being internally heated. The continuum peak detected at lambda=2.8 mm is consistent with the free-free emission from component C of the ultracompact H II region. Velocity structure and morphology outlined by the different tracers suggest that the hot core is primarily energized by component C. Emission from the N- and O-bearing molecules peak at different positions within the innermost regions of the core; none is coincident with the continuum peak. Based on the brightness temperatures of optically thick lines in our sample, we estimate the kinetic temperature of the inner regions of the HMC to be 160+-30 K. Comparison of the observed abundances of the different species in G34.26+0.15 with existing models does not produce a consistent picture.

We present NIR and MIR observations of eight embedded young stellar groups around Herbig Ae/Be stars (HAEBEs) using archived Spitzer IRAC data and 2MASS data. These young stellar groups are nearby ($\leq$ 1 kpc) and still embedded within their molecular clouds. In order to identify the young stellar objects in our sample, we use the color-color diagram of J - [3.6] vs. Ks - [4.5]. The Spitzer images of our sample show that the groups around HAEBEs, spectral types earlier than B8, are usually associated with bright infrared nebulosity. Within this, there are normally 10 - 50 young stars distributed close to the HAEBEs ($<$ 1 pc). Not only are there young stars around the HAEBEs, there are also young stellar populations throughout the whole cloud, some are distributed and some are clumped. The groups around the HAEBEs are sub-structures of the large young population within the molecular cloud. The sizes of groups are also comparable with those sub-structures seen in massive clusters. Young stars in groups around HAEBEs have generally larger SED slopes compared to those outside, which suggests that the young stars in groups are probably younger than the distributed systems. This might imply that there is usually a higher and more continuous star forming rate in groups, that the formation of groups initiates later, or that low mass stars in groups form slower than those outside. Finally, there is no obvious trend between the SED slopes and the distance to the HAEBEs for those young stars within the groups. This suggests that the clustering of young stars dominates over the effect of massive stars on the low-mass young stars at the scale of our study.

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.

We construct a simple model for radioisotopic enrichment of the protosolar nebula by injection from a nearby supernova, based on the inverse square law for ejecta dispersion. We find that the presolar radioisotopes abundances (i.e., in solar masses) demand a nearby supernova: its distance can be no larger than 66 times the size of the protosolar nebula, at a 90% confidence level, assuming 1 solar mass of protosolar material. The relevant size of the nebula depends on its state of evolution at the time of radioactivity injection. In one scenario, a collection of low-mass stars, including our sun, formed in a group or cluster with an intermediate- to high-mass star that ended its life as a supernova while our sun was still a protostar, a starless core, or perhaps a diffuse cloud. Using recent observations of protostars to estimate the size of the protosolar nebula constrains the distance of the supernova at 0.02 to 1.6 pc. The supernova distance limit is consistent with the scales of low-mass stars formation around one or more massive stars, but it is closer than expected were the sun formed in an isolated, solitary state. Consequently, if any presolar radioactivities originated via supernova injection, we must conclude that our sun was a member of such a group or cluster that has since dispersed, and thus that solar system formation should be understood in this context. In addition, we show that the timescale from explosion to the creation of small bodies was on the order of 1.8 Myr (formal 90% confidence range of 0 to 2.2 Myr), and thus the temporal choreography from supernova ejecta to meteorites is important. Finally, we can not distinguish between progenitor masses from 15 to 25 solar masses in the nucleosynthesis models; however, the 20 solar mass model is somewhat preferred.

We present high resolution 0.45" x 0.32" observations from the BIMA array toward the Herbig Ae system MWC 480 in the lambda = 1.4 mm dust continuum. We resolve a circumstellar disk of radius ~170 AU and constrain the disk parameters by comparing the observations to flat disk models. These results show that the typical fit parameters of the disk, such as the mass, Md ~ 0.04-0.18 Mo, and the surface density power law index, p=0.5 or 1, are comparable to those of the lower mass T Tauri stars. The dust in the MWC 480 disk can be modeled as processed dust material (beta ~ 0.8), similar to the Herbig Ae star CQ Tau disk; the fitted disk parameters are also consistent with less-evolved dust (beta ~ 1.2). The possibility of grain growth in the MWC 480 circumstellar disk is supported by the acceptable fits with beta ~ 0.8. The surface density power-law profiles of p=0.5 and p=1 can be easily fit to the MWC 480 disk; however, a surface density power-law profile similiar to the minimum mass solar nebula model p=1.5 is ruled out at an 80% confidence level.

We present an interferometric search for large molecules, including methanol, methyl cyanide, ethyl cyanide, ethanol, and methyl formate in comets LINEAR (C/2002 T7) and NEAT (C/2001 Q4) with the Berkeley-Illinois-Maryland Association (BIMA) array. In addition, we also searched for transitions of the simpler molecules CS, SiO, HNC, HN13C and 13CO . We detected transitions of methanol and CS around Comet LINEAR and one transition of methanol around Comet NEAT within a synthesized beam of ~20''. We calculated the total column density and production rate of each molecular species using the variable temperature and outflow velocity (VTOV) model described by Friedel et al.(2005).Considering the molecular production rate ratios with respect to water, Comet T7 LINEAR is more similar to Comet Hale-Bopp while Comet Q4 NEAT is more similar to Comet Hyakutake. It is unclear, however, due to such a small sample size, whether there is a clear distinction between a Hale-Bopp and Hyakutake class of comet or whether comets have a continuous range of molecular production rate ratios.

We present BIMA and SCUBA observations of the young cluster associated with BD +40 4124 in the dense molecular gas tracer CS J=2-1 and the continuum dust continuum at 3.1 mm and 850 micron. The dense gas and dust in the system is aligned into a long ridge morphology extending ~0.4 pc with 16 gas clumps of estimated masses ranging from 0.14 - 1.8 Solar masses. A north-south variation in the CS center line velocity can be explained with a two cloud model. We posit that the BD +40 4124 stellar cluster formed from a cloud-cloud collision. The largest linewidths occur near V1318 Cyg-S, a massive star affecting its natal environment. In contrast, the dense gas near the other, more evolved, massive stars displays no evidence for disruption; the material must either be processed into the star, dissipate, or relax, fairly quickly. The more evolved low-mass protostars are more likely to be found near the massive stars. If the majority of low-mass stars are coeval, the seemingly evolved low-mass protostars are not older: the massive stars have eroded their structures. Finally, at the highest resolution, the 3.1 mm dust emission is resolved into a flattened structure 3100 by 1500 AU with an estimated mass of 3.4 Solar masses. The continuum and CS emission are offset by 1.1 arcsec from the southern binary source. A simple estimate of the extinction due to the continuum emission structure is AV = 700 mag. From the offset and as the southern source is detected in the optical, the continuum emission is from a previously unknown very young, intermediate-mass, embedded stellar object. (see high resolution figures at http://eeyore.astro.uiuc.edu/~lwl/recent/bd40.pdf).

Using Spitzer Space Telescope and Hubble Space Telescope observations of the superbubble N51D, we have identified three young stellar objects (YSOs) in dust globules, and made the first detection of a Herbig-Haro object outside the Galaxy. The spectral energy distributions of these YSOs suggest young massive stars with disk, envelope, and outflow cavities. The interstellar conditions are used to assess whether the star formation was spontaneous or induced by external pressure.

We present interferometric detections of HCN in comets LINEAR (C/2002 T7) and NEAT (C/2001 Q4) with the Berkeley-Illinois-Maryland Association (BIMA) Array. With a 25.4"x20.3" synthesized beam around Comet LINEAR and using a variable temperature and outflow velocity (VTOV) model, we found an HCN column density of N_T=6.4+-2.1x10^12 cm^-2, and a production rate of Q(HCN)=6.5+-2.2x10^26 s^-1, giving a production rate ratio of HCN relative to H_2O of ~3.3+-1.1x10^-3 and relative to CN of ~4.6+-1.5. With a 21.3"x17.5" synthesized beam around Comet NEAT and using a VTOV model, we found an HCN column density of N_T=8.5+-4.5x10^11 cm^-2, and a production rate of Q(HCN)=8.9+-4.7x10^25 s^-1, giving a production rate ratio of HCN relative to H_2O of ~7.4+-3.9x10^-4 and relative to CN of ~0.3+-0.2. For both comets, the production rates relative to H_2O are similar to those found in previous comet observations. For Comet LINEAR the production rate relative to CN is consistent with HCN being the primary parent species of CN, while for Comet NEAT it is too low for this to be the case.

We present modeling results for six of the eleven deeply embedded systems from our sub-arcsecond 2.7 mm wavelength continuum interferometric survey. The modeling, performed in the uv plane, assumes dust properties, allows for a power-law density profile, uses a self-consistent, luminosity conserving temperature profile, and has an embedded point source to represent a circumstellar disk. Even though we have the highest spatial resolution to date at these wavelengths, only the highest signal-to-noise systems can adequately constrain the simple self-similar collapse models. Of the six sources modeled, all six were fit with a density power-law index of 2.0; however, in half of the systems, those with the highest signal-to-noise, a density power-law index of 1.5 can be rejected at the 95% confidence level. Further, we modeled the systems using the pure Larson-Penston (LP) and Shu solutions with only age and sound speed as parameters. Overall, the LP solution provides a better fit to the data, both in likelihood and providing the observed luminosity, but the age of the systems required by the fits are surprising low (1000-2000 yrs). We suggest that either there is some overall time scaling of the self-similar solutions that invalidate the age estimates, or more likely we are at the limit of the usefulness of these models. With our observations we have begun to reach the stage where models need to incorporate more of the fundamental physics of the collapse process, probably including magnetic fields and/or turbulence. In addition to constraining collapse solutions, our modeling allows the separation of large-scale emission from compact emission, enabling the probing of the circumstellar disk component embedded within the protostellar envelope.

We have made deep Chandra observations of three powerful FRII radio sources: two quasars (3C263 and 3C351) and one radio galaxy (3C330). X-ray emission from hotspots and lobes, as well as from the active nucleus, is detected in each source. We model the hotspots' synchrotron spectra using VLA, BIMA and HST data. In 3C263 and 3C330, the hotspots' X-ray emission is at a level consistent with being synchrotron self-Compton (SSC) emission, with a hotspot magnetic field close to the equipartition value. In the two hotspots of 3C351, however, an SSC origin for the X-rays would require the magnetic field strength to be an order of magnitude below the equipartition value in our models: in addition, there are offsets between the radio, optical and X-ray emission from the secondary hotspot which are hard to explain in a simple SSC model. We discuss the emission mechanisms that may be responsible for these observations. On our preferred model, the X-ray emission from the radio lobes of the three sources is due to inverse-Compton scattering of the microwave background radiation. If this is the case, the magnetic field strengths in the lobes are typically about a factor 2 below the equipartition values, assuming uniform lobe electron and magnetic field distributions. We detect extended X-ray emission, which we attribute to a cluster/group environment, around 3C263 and 3C330. This detection allows us to show that the lobes are close to pressure balance with their surroundings, as long as no non-radiating particles contribute to the internal pressure of the lobes.

Progress in understanding the embedded stars in LkHa225 has been hampered by their variability, making it hard to compare data taken at different times, and by the limited resolution of the available data, which cannot probe the small scales between the two stars. In an attempt to overcome these difficulties, we present new near-infrared data on this object taken using the ALFA adaptive optics system with the MPE 3D integral field spectrometer and the near-infrared camera Omega-Cass. The stars themselves have K-band spectra which are dominated by warm dust emission, analagous to class I-II for low mass YSOs, suggesting that the stars are in a phase where they are still accreting matter. On the other hand, the ridge of continuum emission between them is rather bluer, suggestive of extincted and/or scattered stellar light rather than direct dust emission. The compactness of the CO emission seen toward each star argues for accretion disks (which can also account for much of the K-band veiling) rather than a neutral wind. In contrast to other YSOs with CO emission, LkHa225 has no detectable Br_gamma emission. Additionally there is no H_2 detected on the northern star, although we do confirm that the strongest H_2 emission is on the southern star, where we find it is excited primarily by thermal mechanisms. A second knot of H_2 is observed to its northeast, with a velocity shift of -75kms and a higher fraction of non-thermal emission. This is discussed with reference to the H2O maser, the molecular outflow, and [S II] emission observed between the stars.

We present 85-GHz observations of the archetypal double-hotspot radio source 3C20 made with the BIMA millimetre array. The resolution of BIMA allows us to separate the two components of the eastern hotspot. By comparing the BIMA observations with existing VLA data, we show that the spectra of the two hotspot components are very similar, despite the clear differences in their radio structure and their wide separation. We discuss the implications for models of double hotspot formation. Weak emission from the lobes of 3C20 is detected at 85 GHz, at a level consistent with the predictions of standard spectral ageing models.

We are developing two-dimensional 16 x 25 pixel detector arrays of both unstressed and stressed Ge:Ga photoconductive detectors for far-infrared astronomy from SOFIA. The arrays, based on earlier 5 x 5 detector arrays used on the KAO, will be for our new instrument, the Far Infrared Field Imaging Line Spectrometer (FIFI LS). The unstressed Ge:Ga detector array will cover the wavelength range from 40 to 120 microns, and the stressed Ge:Ga detector array from 120 to 210 microns. The detector arrays will be operated with multiplexed integrating amplifiers with cryogenic readout electronics located close to the detector arrays. The design of the stressed detector array and results of current measurements on several prototype 16 pixel linear arrays are reported. They demonstrate the feasibility of the current concept. ***This paper does not include Figures due to astro-ph size limitations. Please download entire file at http://fifi-ls.mpe-garching.mpg.de/spie.det.ps.gz ***

We present a progress report on the design and construction of the Field-Imaging Far-Infrared Line Spectrometer (FIFI LS) for the SOFIA airborne observatory. The design of the instrument is driven by the goal of maximizing observing efficiency, especially for observations of faint, extragalactic objects. Thus, FIFI LS utilizes an integral field technique that uses slicer mirrors to optically re-arrange the two-dimensional field into a single slit for a long slit spectrometer. Effectively, a 5x5 pixel spatial field of view is imaged to a 25x1 pixel slit and dispersed to a 25x16 pixel, two-dimensional detector array, providing diffraction-limited spatial and spectral multiplexing. In this manner, the instrument employs two parallel, medium resolution (R~2000) grating spectrometers for simultaneous observations in two bands: a short wavelength band (42 to 110 micron) and a long wavelength band (110 to 210 micron). Overall, for each of the 25 spatial pixels, the instrument can cover a velocity range of ~1500 km/s around selected far-infrared lines with an estimated sensitivity of 2x10^-15 W Hz^1/2 per pixel. This arrangement provides good spectral coverage with high responsivity. ***This paper does not include Figures due to astro-ph size limitations. Please download entire file at http://fifi-ls.mpe-garching.mpg.de/fifils.ps.gz ****

We present the results of sub-arcsecond 108 GHz continuum interferometric observations toward the radio luminous galaxy 3C123. Using multi-array observations, we utilize the high u,v dynamic range of the BIMA millimeter array to sample fully spatial scales ranging from 0.5" to 50". This allows us to make one-to-one comparisons of millimeter-wavelength emission in the radio lobes and hotspots to VLA centimeter observations at 1.4, 4.9, 8.4, and 15 GHz. At 108 GHz, the bright, eastern double hotspot in the southern lobe is resolved. This is only the second time that a multiple hotspot region has been resolved in the millimeter regime. We model the synchrotron spectra of the hotspots and radio lobes using simple broken power-law models with high energy cutoffs, and discuss the hotspot spectra and their implications for models of multiple hotspot formation.

We present the results of a 2.7 mm continuum interferometric survey of 24 young stellar objects in 11 fields. The target objects range from deeply embedded Class 0 sources to optical T Tauri sources. This is the first sub-arcsecond survey of the 2.7 mm dust continuum emission from young, embedded stellar systems. The images show a diversity of structure and complexity. The optically visible T Tauri stars (DG Tauri, HL Tauri, GG Tauri,and GM Aurigae) have continuum emission dominated by compact, less than 1", circumstellar disks. The more embedded near-infrared sources (SVS13 and L1551 IRS5) have continuum emission that is extended and compact. The embedded sources (L1448 IRS3, NGC1333 IRAS2, NGC1333 IRAS4, VLA1623, and IRAS 16293-2422) have continuum emission dominated by the extended envelope, typically more than 85%. In fact, in many of the deeply embedded systems it is difficult to uniquely isolate the disk emission component from the envelope extending inward to AU size scales. All of the target embedded objects are in multiple systems with separations on scales of 30" or less. Based on the system separation, we place the objects into three categories: separate envelope (separation > 6500 AU), common envelope (separation 150-3000 AU), and common disk (separation < 100 AU). These three groups can be linked with fragmentation events during the star formation process: separate envelopes from prompt initial fragmentation and the separate collapse of a loosely condensed cloud, common envelopes from fragmentation of a moderately centrally condensed spherical system, and common disk from fragmentation of a high angular momentum circumstellar disk.

We present sub-arcsecond resolution imaging of the $\lambda$ = 2.7 mm continuum emission from the young, embedded system L1551 IRS5 using the nine-element, high-resolution configuration of the BIMA array. The observed emission arises from two compact sources separated by $0\farcs$35, coinciding with the two sources seen at $\lambda$ = 2 cm and $\lambda$ = 1.3 cm. When the high resolution data is combined with data from two compact configurations, L1551 IRS5 is argued to consist of a protobinary system separated by $\sim$50 AU with individual circumstellar disks, a circumbinary structure, and a large-scale envelope. The characteristic masses of the components are: 0.024 M$_{\sun}$ for the northern circumstellar disk, 0.009 M$_{\sun}$ for the southern circumstellar disk, 0.04 M$_{\sun}$ for the circumbinary material, and 0.28 M$_{\sun}$ for the envelope.