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

Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and size scales. During the last decade a variety of technological and computing advances have transformed our understanding of star formation through the use of multi-wavelength observations, large scale observational surveys, and multi-physics multi-dimensional numerical simulations. Additionally, the use of synthetic observations of simulations have provided a useful tool to interpret observational data and evaluate the importance of various physical processes on different scales in star formation. Here, we review these recent advancements in both high- ($M \gtrsim 8 \, M_{\rm \odot}$) and low-mass star formation.

Previous continuum observations from the MUSTANG camera on the Green Bank Telescope (GBT) of the nearby star-forming filament OMC 2/3 found elevated emission at 3.3 mm relative to shorter wavelength data. As a consequence, the inferred dust emissivity index obtained from modified black body dust spectra was considerably lower than what is typically measured on $\sim 0.1 \, {\rm pc}$ scales in nearby molecular clouds. Here we present new observations of OMC 2/3 collected with the MUSTANG-2 camera on the GBT which confirm this elevated emission. We also present for the first time sensitive 1 cm observations made with the Ka-band receiver on the GBT which also show higher than expected emission. We use these observations--- along with Herschel, JCMT, Mambo, and GISMO data--- to assemble spectral energy distributions (SEDs) of a variety of structures in OMC 2/3 spanning the range $160 \, {\rm \mu m}$ to $1 \, {\rm cm}$. The data at 2 mm and shorter are generally consistent with a modified black body spectrum and a single value of $\beta \sim 1.6$. The 3 mm and 1 cm data, however, lie well above such an SED. The spectrum of the long wavelength excess is inconsistent with both free-free emission and standard "Spinning Dust" models for Anomalous Microwave Emission (AME). The 3 mm and 1 cm data could be explained by a flatter dust emissivity at wavelengths shorter than 2 mm, potentially in concert with AME in some regions.

Magnetic fields are ubiquitous in our Universe, but remain poorly understood in many branches of astrophysics. A key tool for inferring astrophysical magnetic field properties is dust emission polarimetry. The James Clerk Maxwell Telescope (JCMT) is planning a new 850$\mu$m camera consisting of an array of 7272 paired Microwave Kinetic Inductance Detectors (MKIDs), which will inherently acquire linear polarization information. The camera will allow wide-area polarization mapping of dust emission at 14$^{\prime\prime}$-resolution, allowing magnetic field properties to be studied in a wide range of environments, including all stages of the star formation process, Asymptotic Giant Branch stellar envelopes and planetary nebula, external galaxies including starburst galaxies and analogues for the Milky Way, and the environments of active galactic nuclei (AGN). Time domain studies of AGN and protostellar polarization variability will also become practicable. Studies of the polarization properties of the interstellar medium will also allow detailed investigation of dust grain properties and physics. These investigations would benefit from a potential future upgrade adding 450$\mu$m capability to the camera, which would allow inference of spectral indices for polarized dust emission in a range of environments. The enhanced mapping speed and polarization capabilities of the new camera will transform the JCMT into a true submillimetre polarization survey instrument, offering the potential to revolutionize our understanding of magnetic fields in the cold Universe.

We report ALMA and VLA continuum observations that potentially identify the four youngest protostars in the Orion Molecular Clouds taken as part of the Orion VANDAM program. These are distinguished by bright, extended, irregular emission at 0.87 mm and 8 mm and are optically thick at 0.87 mm. These structures are distinct from the disk or point-like morphologies seen toward the other Orion protostars. The 0.87 mm emission implies temperatures of 41-170 K, requiring internal heating. The bright 8 mm emission implies masses of 0.5 to 1.2 M_sun assuming standard dust opacity models. One source has a Class 0 companion, while another exhibits substructure indicating a companion-candidate. Three compact outflows are detected, two of which may be driven by companions, with dynamical times of ~300 to ~400 years. The slowest outflow may be driven by a first hydrostatic core. These protostars appear to trace an early phase when the centers of collapsing fragments become optically thick to their own radiation and compression raises the gas temperature. This phase is thought to accompany the formation of hydrostatic cores. A key question is whether these structures are evolving on free fall times of ~100 years, or whether they are evolving on Kelvin-Helmholtz times of several thousand years. The number of these sources imply a lifetime of ~6000 years, in closer agreement with the Kelvin-Helmholtz time. In this case, rotational and/or magnetic support could be slowing the collapse.

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.

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $\mu$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/.

Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

The next generation Very Large Array (ngVLA) is a transformational radio observatory being designed by the U.S. National Radio Astronomy Observatory (NRAO). It will provide order of magnitude improvements in sensitivity, resolution, and uv coverage over the current Jansky Very Large Array (VLA) at ~1.2-50 GHz and extend the frequency range up to 70-115 GHz. This document is a white paper written by members of the Canadian community for the 2020 Long Range Plan panel, which will be making recommendations on Canada's future directions in astronomy. Since Canadians have been historically major users of the VLA and have been valued partners with NRAO for ALMA, Canada's participation in ngVLA is welcome. Canadians have been actually involved in ngVLA discussions for the past five years, and have played leadership roles in the ngVLA Science and Technical Advisory Councils. Canadian technologies are also very attractive for the ngVLA, in particular our designs for radio antennas, receivers, correlates, and data archives, and our industrial capacities to realize them. Indeed, the Canadian designs for the ngVLA antennas and correlator/beamformer are presently the baseline models for the project. Given the size of Canada's radio community and earlier use of the VLA (and ALMA), we recommend Canadian participation in the ngVLA at the 7% level. Such participation would be significant enough to allow Canadian leadership in gVLA's construction and usage. Canada's participation in ngVLA should not preclude its participation in SKA; access to both facilities is necessary to meet Canada's radio astronomy needs. Indeed, ngVLA will fill the gap between those radio frequencies observable with the SKA and ALMA at high sensitivities and resolutions. Canada's partnership in ngVLA will give it access to cutting-edge facilities together covering approximately three orders of magnitude in frequency.

(abridged) The Atacama Large Millimeter/submillimeter Array (ALMA) was the top-ranked priority for a new ground-based facility in the 2000 Canadian Long Range Plan. Ten years later, at the time of LRP2010, ALMA construction was well underway, with first science observations anticipated for 2011. In the past 8 years, ALMA has proved itself to be a high-impact, high-demand observatory, with record numbers of proposals submitted to the annual calls and large numbers of highly cited scientific papers across fields from protoplanetary disks to high-redshift galaxies and quasars. The LRP2010 ALMA white paper laid out 8 specific metrics that could be used to judge the success of Canada's participation in ALMA. Among these metrics were publications (number; impact), collaborations (international; multi-wavelength), and student training. To call out one particular metric, Canadians are making excellent use of ALMA in training graduate students and postdocs: as of June 2018, 12 of 23 Canadian first-author papers were led by a graduate student, and a further 4 papers were led by postdocs. All 8 metrics argue for Canada's involvement in ALMA over the past decade to be judged a success. The successful achievement of these wide-ranging goals argues strongly for Canada's continuing participation in ALMA over the next decade and beyond. Looking forward, our community needs to: (1) maintain Canadian access to ALMA and our competitiveness in using ALMA; (2) preserve full Canadian funding for our share of ALMA operations; (3) identify components of ALMA development in which Canada can play a significant role, including stimulating expertise in submillimetre instrumentation to capitalize on future opportunities; and (4) keep Canadians fully trained and engaged in ALMA, as new capabilities become available, reaching the widest possible community of potential users.

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

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

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

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.

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

It is usually assumed that the stellar initial mass function (IMF) takes a universal form and that there exists a direct mapping between this and the distribution of natal core masses (the core mass function, CMF). The IMF and CMF have been best characterized in the Solar neighborhood. Beyond 500~pc from the Sun, in diverse environments where metallicity varies and massive star feedback may dominate, the IMF has been measured only incompletely and imprecisely, while the CMF has hardly been measured at all. In order to establish if the IMF and CMF are indeed universal and related to each other, it is necessary to: 1) perform multi-wavelength large-scale imaging and spectroscopic surveys of different environments across the Galaxy; 2) require an angular resolution of < 0.1'' in the optical/near-IR for stars and < 5'' in the far-IR for cores; 3) achieve far-IR sensitivities to probe 0.1~Msol cores at 2--3 kpc.

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.

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

Background: low-mass stars are the dominant product of the star formation process, and they trace star formation over the full range of environments, from isolated globules to clusters in the central molecular zone. In the past two decades, our understanding of the spatial distribution and properties of young low-mass stars and protostars has been revolutionized by sensitive space-based observations at X-ray and IR wavelengths. By surveying spatial scales from clusters to molecular clouds, these data provide robust measurements of key star formation properties. Goal: with their large numbers and their presence in diverse environments, censuses of low mass stars and protostars can be used to measure the dependence of star formation on environmental properties, such as the density and temperature of the natal gas, strengths of the magnetic and radiation fields, and the density of stars. Here we summarize how such censuses can answer three basic questions: i.) how is the star formation rate influenced by environment, ii.) does the IMF vary with environment, and iii.) how does the environment shape the formation of bound clusters? Answering these questions is an important step toward understanding star and cluster formation across the extreme range of environments found in the Universe. Requirements: sensitivity and angular resolution improvements will allow us to study the full range of environments found in the Milky Way. High spatial dynamic range (< 1arcsec to > 1degree scales) imaging with space-based telescopes at X-ray, mid-IR, and far-IR and ground-based facilities at near-IR and sub-mm wavelengths are needed to identify and characterize young stars.

We present the severe, systematic uncertainties currently facing our understanding of dust emission, which stymie our ability to truly exploit dust as a tool for studying galaxy evolution. We propose a program of study to tackle these uncertainties, describe the necessary facilities, and discuss the potential science gains that will result. This white paper was submitted to the US National Academies' Astro2020 Decadal Survey on Astronomy and Astrophysics.

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.

Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variability due to accretion bursts on dynamical and chemical evolution of the source is largely unknown. The lasting effects on protostellar envelopes and disks are tracked through multi-wavelength and time domain observational campaigns, requiring deep X-ray, infrared, and radio imaging and spectroscopy, at a sufficient level of spatial detail to distinguish contributions from the various substructures (i.e., envelope from disk from star from outflow). Protostellar models derived from these campaigns will illuminate the initial chemical state of protoplanetary disks during the epoch of giant planet formation. Insight from individual star formation in the Milky Way is also necessary to understand star formation rates in extragalactic sources. This cannot be achieved with ground-based observatories and is not covered by currently approved instrumentation. Requirements: High (v < 10 km/s for survey; v < 1 km/s for followup) spectral resolution capabilities with relatively rapid response times in the IR (3-500 um), X-ray (0.1-10 keV), and radio (cm) are critical to follow the course of accretion and outflow during an outburst. Complementary, AU-scale radio observations are needed to probe the disk accretion zone, and 10 AU-scale to probe chemical and kinematic structures of the disk-forming regions, and track changes in the dust, ice, and gas within protostellar envelopes.

The past two decades have seen extensive surveys of the far-infrared to submillimeter continuum emission in the plane of our Galaxy. We line out prospects for the coming decade for corresponding molecular and atomic line surveys which are needed to fully understand the formation of the dense structures that give birth to clusters and stars out of the diffuse interstellar medium. We propose to work towards Galaxy wide surveys in mid-J CO lines to trace shocks from colliding clouds, Galaxy-wide surveys for atomic Carbon lines in order to get a detailed understanding of the relation of atomic and molecular gas in clouds, and to perform extensive surveys of the structure of the dense parts of molecular clouds to understand the importance of filaments/fibers over the full range of Galactic environments and to study how dense cloud cores are formed from the filaments. This work will require a large (50m) Single Dish submillimeter telescope equipped with massively multipixel spectrometer arrays, such as envisaged by the AtLAST project.

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

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

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.

The JCMT Gould Belt Survey was one of the first Legacy Surveys with the James Clerk Maxwell Telescope in Hawaii, mapping 47 square degrees of nearby (< 500 pc) molecular clouds in both dust continuum emission at 850 $\mu$m and 450 $\mu$m, as well as a more-limited area in lines of various CO isotopologues. While molecular clouds and the material that forms stars have structures on many size scales, their larger-scale structures are difficult to observe reliably in the submillimetre regime using ground-based facilities. In this paper, we quantify the extent to which three subsequent data-reduction methods employed by the JCMT GBS accurately recover emission structures of various size scales, in particular, dense cores which are the focus of many GBS science goals. With our current best data-reduction procedure, we expect to recover 100% of structures with Gaussian sigma sizes of $\le$30" and intensity peaks of at least five times the local noise for isolated peaks of emission. The measured sizes and peak fluxes of these compact structures are reliable (within 15% of the input values), but source recovery and reliability both decrease significantly for larger emission structures and for fainter peaks. Additional factors such as source crowding have not been tested in our analysis. The most recent JCMT GBS data release includes pointing corrections, and we demonstrate that these tend to decrease the sizes and increase the peak intensities of compact sources in our dataset, mostly at a low level (several percent), but occasionally with notable improvement.

We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES .

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

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

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

We present a study of hierarchical structure in the Perseus molecular cloud, from the scale of the entire cloud ($\gtrsim$10 pc) to smaller clumps ($\sim$1 pc), cores ($\sim$0.05-0.1 pc), envelopes ($\sim$300-3000 AU) and protostellar objects ($\sim$15 AU). We use new observations from the Submillimeter Array (SMA) large project "Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES)" to probe the envelopes, and recent single-dish and interferometric observations from the literature for the remaining scales. This is the first study to analyze hierarchical structure over five scales in the same cloud complex. We compare the number of fragments with the number of Jeans masses in each scale to calculate the Jeans efficiency, or the ratio of observed to expected number of fragments. The velocity dispersion is assumed to arise either from purely thermal motions, or from combined thermal and non-thermal motions inferred from observed spectral line widths. For each scale, thermal Jeans fragmentation predicts more fragments than observed, corresponding to inefficient thermal Jeans fragmentation. For the smallest scale, thermal plus non-thermal Jeans fragmentation also predicts too many protostellar objects. However at each of the larger scales thermal plus non-thermal Jeans fragmentation predicts fewer than one fragment, corresponding to no fragmentation into envelopes, cores, and clumps. Over all scales, the results are inconsistent with complete Jeans fragmentation based on either thermal or thermal plus non-thermal motions. They are more nearly consistent with inefficient thermal Jeans fragmentation, where the thermal Jeans efficiency increases from the largest to the smallest scale.

We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point corrections for small Herschel maps. We use radial profiles to identify offsets between the observed HSA intensities and the expected intensities from Planck. Most clouds have reliable offset measurements with this technique. In addition, we find that roughly half of the clouds have underestimated HSA-SPIRE intensities in their outer envelopes relative to Planck, even though the HSA-SPIRE maps were previously zero-point corrected. Using our technique, we produce corrected Herschel intensity maps for all 56 clouds and determine their line-of-sight average dust temperatures and optical depths from modified black body fits. The clouds have typical temperatures of ~ 14-20 K and optical depths of ~ 1e-5 - 1e-3. Across the whole sample, we find an anti-correlation between temperature and optical depth. We also find lower temperatures than what was measured in previous Herschel studies, which subtracted out a background level from their intensity maps to circumvent the zero-point correction. Accurate Herschel observations of clouds are key to obtain accurate density and temperature profiles. To make such future analyses possible, intensity-corrected maps for all 56 clouds are publicly available in the electronic version.

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.

Density profiles of isolated cores derived from thermal dust continuum emission rely on models of dust properties, such as mass opacity, which are poorly constrained. With complementary measures from near-infrared extinction maps, we can assess the reliability of commonly-used dust models. In this work, we compare Herschel-derived maps of the optical depth with equivalent maps derived from CFHT WIRCAM near-infrared observations for three isolated cores: CB68, L429, and L1552. We assess the dust opacities provided from four models: OH1a, OH5a, Orm1, and Orm4. Although the consistency of the models differs between the three sources, the results suggest that the optical properties of dust in the envelopes of the cores are best described by either silicate and bare graphite grains (e.g., Orm1) or carbonaceous grains with some coagulation and either thin or no ice mantles (e.g., OH5a). None of the models, however, individually produced the most consistent optical depth maps for every source. The results suggest that either the dust in the cores is not well described by any one dust property model, the application of the dust models cannot be extended beyond the very center of the cores, or more complex SED fitting functions are necessary.

We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from $Herschel$ observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of $\sim$3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation.

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

We explore the relationship between young, embedded binaries and their parent cores, using observations within the Perseus Molecular Cloud. We combine recently published VLA observations of young stars with core properties obtained from SCUBA-2 observations at 850 um. Most embedded binary systems are found toward the centres of their parent cores, although several systems have components closer to the core edge. Wide binaries, defined as those systems with physical separations greater than 500 au, show a tendency to be aligned with the long axes of their parent cores, whereas tight binaries show no preferred orientation. We test a number of simple, evolutionary models to account for the observed populations of Class 0 and I sources, both single and binary. In the model that best explains the observations, all stars form initially as wide binaries. These binaries either break up into separate stars or else shrink into tighter orbits. Under the assumption that both stars remain embedded following binary breakup, we find a total star formation rate of 168 Myr^-1. Alternatively, one star may be ejected from the dense core due to binary breakup. This latter assumption results in a star formation rate of 247 Myr^-1. Both production rates are in satisfactory agreement with current estimates from other studies of Perseus. Future observations should be able to distinguish between these two possibilities. If our model continues to provide a good fit to other star-forming regions, then the mass fraction of dense cores that becomes stars is double what is currently believed.

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

We present the results of our investigation of the star-forming potential in the Perseus star-forming complex. We build on previous starless core, protostellar core, and young stellar object (YSO) catalogs from Spitzer, Herschel, and SCUBA observations in the literature. We place the cores and YSOs within seven star-forming clumps based on column densities greater than 5x10^21 cm^-2. We calculate the mean density and free-fall time for 69 starless cores as 5.55x10^-19 gcm^-3 and 0.1 Myr,respectively, and we estimate the star formation rate for the near future as 150 Msun Myr^-1. According to Bonnor Ebert stability analysis, we find that majority of starless cores in Perseus are unstable. Broadly, these cores can collapse to form the next generation of stars. We found a relation between starless cores and YSOs, where the numbers of young protostars (Class 0 + Class I) are similar to the numbers of starless cores. This similarity, which shows a one-to-one relation, suggests that these starless cores may form the next generation of stars with approximately the same formation rate as the current generation, as identified by the Class 0 and Class I protostars. It follows that if such a relation between starless cores and any YSO stage exists, the SFR values of these two populations must be nearly constant. In brief, we propose that this one-to-one relation is an important factor in better understanding the star formation process within a cloud.

Compact substructure is expected to arise in a starless core as mass becomes concentrated in the central region likely to form a protostar. Additionally, multiple peaks may form if fragmentation occurs. We present ALMA Cycle 2 observations of 60 starless and protostellar cores in the Ophiuchus molecular cloud. We detect eight compact substructures which are >15 arcsec from the nearest Spitzer YSO. Only one of these has strong evidence for being truly starless after considering ancillary data, e.g., from Herschel and X-ray telescopes. An additional extended emission structure has tentative evidence for starlessness. The number of our detections is consistent with estimates from a combination of synthetic observations of numerical simulations and analytical arguments. This result suggests that a similar ALMA study in the Chamaeleon I cloud, which detected no compact substructure in starless cores, may be due to the peculiar evolutionary state of cores in that cloud.

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

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 observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

The Orion A molecular cloud is one of the most well-studied nearby star-forming regions, and includes regions of both highly clustered and more dispersed star formation across its full extent. Here, we analyze dense, star-forming cores identified in the 850 {\mu}m and 450 {\mu}m SCUBA-2 maps from the JCMT Gould Belt Legacy Survey. We identify dense cores in a uniform manner across the Orion A cloud and analyze their clustering properties. Using two independent lines of analysis, we find evidence that clusters of dense cores tend to be mass segregated, suggesting that stellar clusters may have some amount of primordial mass segregation already imprinted in them at an early stage. We also demonstrate that the dense core clusters have a tendency to be elongated, perhaps indicating a formation mechanism linked to the filamentary structure within molecular clouds.

We present a comparison of SCUBA-2 850-$\mu$m and Herschel 70--500-$\mu$m observations of the L1495 filament in the Taurus Molecular Cloud with the goal of characterising the SCUBA-2 Gould Belt Survey (GBS) data set. We identify and characterise starless cores in three data sets: SCUBA-2 850-$\mu$m, Herschel 250-$\mu$m, and Herschel 250-$\mu$m spatially filtered to mimic the SCUBA-2 data. SCUBA-2 detects only the highest-surface-brightness sources, principally detecting protostellar sources and starless cores embedded in filaments, while Herschel is sensitive to most of the cloud structure, including extended low-surface-brightness emission. Herschel detects considerably more sources than SCUBA-2 even after spatial filtering. We investigate which properties of a starless core detected by Herschel determine its detectability by SCUBA-2, and find that they are the core's temperature and column density (for given dust properties). For similar-temperature cores, such as those seen in L1495, the surface brightnesses of the cores are determined by their column densities, with the highest-column-density cores being detected by SCUBA-2. For roughly spherical geometries, column density corresponds to volume density, and so SCUBA-2 selects the densest cores from a population at a given temperature. This selection effect, which we quantify as a function of distance, makes SCUBA-2 ideal for identifying those cores in Herschel catalogues that are closest to forming stars. Our results can now be used by anyone wishing to use the SCUBA-2 GBS data set.

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 the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

Constraining the temperature and density structure of dense molecular cloud cores is fundamental for understanding the initial conditions of star formation. We use Herschel observations of the thermal FIR dust emission from nearby isolated molecular cloud cores and combine them with ground-based submillimeter continuum data to derive observational constraints on their temperature and density structure. The aim of this study is to verify the validity of a ray-tracing inversion technique developed to derive the dust temperature and density structure of isolated starless cores directly from the dust emission maps and to test if the resulting temperature and density profiles are consistent with physical models. Using this ray-tracing inversion technique, we derive the dust temperature and density structure of six isolated starless cloud cores. We employ self-consistent radiative transfer modeling to the derived density profiles, treating the ISRF as the only heating source. The best-fit values of local strength of the ISRF and the extinction by the outer envelope are derived by comparing the self-consistently calculated temperature profiles with those derived by the ray-tracing method. We find that all starless cores are significantly colder inside than outside, with the core temperatures showing a strong negative correlation with peak column density. This suggests that their thermal structure is dominated by external heating from the ISRF and shielding by dusty envelopes. The temperature profiles derived with the ray-tracing inversion method can be well-reproduced with self-consistent radiative transfer models.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

We present an Atacama Large Millimeter/submillimeter Array (ALMA) 106 GHz (Band 3) continuum survey of the complete population of dense cores in the Chamaeleon I molecular cloud. We detect a total of 24 continuum sources in 19 different target fields. All previously known Class 0 and Class I protostars in Chamaeleon I are detected, whereas all of the 56 starless cores in our sample are undetected. We show that the Spitzer+Herschel census of protostars in Chamaeleon I is complete, with the rate at which protostellar cores have been misclassified as starless cores calculated as <1/56, or < 2%. We use synthetic observations to show that starless cores collapsing following the turbulent fragmentation scenario are detectable by our ALMA observations when their central densities exceed ~10^8 cm^-3, with the exact density dependent on the viewing geometry. Bonnor-Ebert spheres, on the other hand, remain undetected to central densities at least as high as 10^10 cm^-3. Our starless core non-detections are used to infer that either the star formation rate is declining in Chamaeleon I and most of the starless cores are not collapsing, matching the findings of previous studies, or that the evolution of starless cores are more accurately described by models that develop less substructure than predicted by the turbulent fragmentation scenario, such as Bonnor-Ebert spheres. We outline future work necessary to distinguish between these two possibilities.

We present a catalogue of dense cores in a $\sim 4^\circ\times2^\circ$ field of the Taurus star-forming region, inclusive of the L1495 cloud, derived from Herschel SPIRE and PACS observations in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, and 500 $\mu$m continuum bands. Estimates of mean dust temperature and total mass are derived using modified blackbody fits to the spectral energy distributions. We detect 525 starless cores of which $\sim10$-20% are gravitationally bound and therefore presumably prestellar. Our census of unbound objects is $\sim85$% complete for $M>0.015\,M_\odot$ in low density regions ($A_V\stackrel{<}{_\sim}5$ mag), while the bound (prestellar) subset is $\sim85$% complete for $M>0.1\,M_\odot$ overall. The prestellar core mass function (CMF) is consistent with lognormal form, resembling the stellar system initial mass function, as has been reported previously. All of the inferred prestellar cores lie on filamentary structures whose column densities exceed the expected threshold for filamentary collapse, in agreement with previous reports. Unlike the prestellar CMF, the unbound starless CMF is not lognormal, but instead is consistent with a power-law form below $0.3\,M_\odot$ and shows no evidence for a low-mass turnover. It resembles previously reported mass distributions for CO clumps at low masses ($M\stackrel{<}{_\sim}0.3\,M_\odot$). The volume density PDF, however, is accurately lognormal except at high densities. It is consistent with the effects of self-gravity on magnetized supersonic turbulence. The only significant deviation from lognormality is a high-density tail which can be attributed unambiguously to prestellar cores.

We present new measurements of the dust emissivity index, beta, for the high-mass, star-forming OMC 2/3 filament. We combine 160-500 um data from Herschel with long-wavelength observations at 2 mm and fit the spectral energy distributions across a ~ 2 pc long, continuous section of OMC 2/3 at 15000 AU (0.08 pc) resolution. With these data, we measure beta and reconstruct simultaneously the filtered-out large-scale emission at 2 mm. We implement both variable and fixed values of beta, finding that beta = 1.7 - 1.8 provides the best fit across most of OMC 2/3. These beta values are consistent with a similar analysis carried out with filtered Herschel data. Thus, we show that beta values derived from spatial filtered emission maps agree well with those values from unfiltered data at the same resolution. Our results contradict the very low beta values (~ 0.9) previously measured in OMC 2/3 between 1.2 mm and 3.3 mm data, which we attribute to elevated fluxes in the 3.3 mm observations. Therefore, we find no evidence or rapid, extensive dust grain growth in OMC 2/3. Future studies with Herschel data and complementary ground-based long-wavelength data can apply our technique to obtain robust determinations of beta in nearby cold molecular clouds.

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 a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

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.

We present and discuss the results of the Herschel Gould Belt survey observations in a ~11 deg^2 area of the Aquila molecular cloud complex at d~260 pc, imaged with the SPIRE/PACS cameras from 70 to 500 micron. We identify a complete sample of starless dense cores and embedded protostars in this region, and analyze their global properties and spatial distributions. We find a total of 651 starless cores, ~60% of which are gravitationally bound prestellar cores, and they will likely form stars in the future. We also detect 58 protostellar cores. The core mass function (CMF) derived for the prestellar cores is very similar in shape to the stellar initial mass function (IMF), supporting the earlier view that there is a close physical link between the IMF and the CMF. The global shift in mass scale observed between the CMF and the IMF is consistent with a typical star formation efficiency of ~40%. By comparing the numbers of starless cores to the number of young stellar objects, we estimate that the lifetime of prestellar cores is ~1 Myr. We find a strong correlation between the spatial distribution of prestellar cores and the densest filaments. About 90% of the Herschel-identified prestellar cores are located above a background column density corresponding to A_V~7, and ~75% of them lie within filamentary structures with supercritical masses per unit length >~16 M_sun/pc. These findings support a picture wherein the cores making up the peak of the CMF (and probably responsible for the base of the IMF) result primarily from the gravitational fragmentation of marginally supercritical filaments. Given that filaments appear to dominate the mass budget of dense gas at A_V>7, our findings also suggest that the physics of prestellar core formation within filaments is responsible for a characteristic "efficiency" SFR/M_dense ~5+-2 x 10^-8 yr^-1 for the star formation process in dense gas.

We present 13CO and C18O (1-0), (2-1), and (3-2) maps towards the core-forming Perseus B1-E clump using observations from the James Clerk Maxwell Telescope (JCMT), Submillimeter Telescope (SMT) of the Arizona Radio Observatory, and IRAM 30 m telescope. We find that the 13CO and C18O line emission both have very complex velocity structures, indicative of multiple velocity components within the ambient gas. The (1-0) transitions reveal a radial velocity gradient across B1-E of 1 km/s/pc that increases from north-west to south-east, whereas the majority of the Perseus cloud has a radial velocity gradient increasing from south-west to north-east. In contrast, we see no evidence of a velocity gradient associated with the denser Herschel-identified substructures in B1-E. Additionally, the denser substructures have much lower systemic motions than the ambient clump material, which indicates that they are likely decoupled from the large-scale gas. Nevertheless, these substructures themselves have broad line widths (0.4 km/s) similar to that of the C18O gas in the clump, which suggests they inherited their kinematic properties from the larger-scale, moderately dense gas. Finally, we find evidence of C18O depletion only toward one substructure, B1-E2, which is also the only object with narrow (transonic) line widths. We suggest that as prestellar cores form, their chemical and kinematic properties are linked in evolution, such that these objects must first dissipate their turbulence before they deplete in CO.

We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

CO, $^{13}$CO and C$^{18}$O ${\it J}$ = 3--2 observations are presented of the Ophiuchus molecular cloud. The $^{13}$CO and C$^{18}$O emission is dominated by the Oph A clump, and the Oph B1, B2, C, E, F and J regions. The optically thin(ner) C$^{18}$O line is used as a column density tracer, from which the gravitational binding energy is estimated to be $4.5 \times 10^{39}$ J (2282 $M_\odot$ km$^2$ s$^{-2}$). The turbulent kinetic energy is $6.3 \times 10^{38}$ J (320 $M_\odot$ km$^2$ s$^{-2}$), or 7 times less than this, and therefore the Oph cloud as a whole is gravitationally bound. Thirty protostars were searched for high velocity gas, with eight showing outflows, and twenty more having evidence of high velocity gas along their lines-of-sight. The total outflow kinetic energy is $1.3 \times 10^{38}$ J (67 $M_\odot$ km$^2$ s$^{-2}$), corresponding to 21$\%$ of the cloud's turbulent kinetic energy. Although turbulent injection by outflows is significant, but does ${\it not}$ appear to be the dominant source of turbulence in the cloud. 105 dense molecular clumplets were identified, which had radii $\sim$ 0.01--0.05 pc, virial masses $\sim$ 0.1--12 $M_\odot$, luminosities $\sim$ 0.001--0.1 K~km s$^{-1}$ pc$^{-2}$, and excitation temperatures $\sim$ 10--50K. These are consistent with the standard GMC based size-line width relationships, showing that the scaling laws extend down to size scales of hundredths of a parsec, and to sub solar-mass condensations. There is however no compelling evidence that the majority of clumplets are undergoing free-fall collapse, nor that they are pressure confined.

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 a new 3.3 mm continuum map of the OMC-2/3 region. When paired with previously published maps of 1.2mm continuum and NH3-derived temperature, we derive the emissivity spectral index of dust emission in this region, tracking its changes across the filament and cores. We find that the median value of the emissivity spectral index is 0.9, much shallower than previous estimates in other nearby molecular clouds. We find no significant difference between the emissivity spectral index of dust in the OMC-2/3 filament and the starless or protostellar cores. Furthermore, the temperature and emissivity spectral index, beta, are anti-correlated at the 4 sigma level. The low values of the emissivity spectral index found in OMC-2/3 can be explained by the presence of millimeter-sized dust grains in the dense regions of the filaments to which these maps are most sensitive. Alternatively, a shallow dust emissivity spectral index may indicate non-powerlaw spectral energy distributions, significant free-free emission, or anomalous microwave emission. We discuss the possible implications of millimeter-sized dust grains compared to the alternatives.

We use PACS and SPIRE continuum data at 160 um, 250 um, 350 um, and 500 um from the Herschel Gould Belt Survey to sample seven clumps in Perseus: B1, B1-E, B5, IC348, L1448, L1455, and NGC1333. Additionally, we identify and characterize the embedded Class 0 protostars using detections of compact Herschel sources at 70 um as well as archival Spitzer catalogues and SCUBA 850 um photometric data. We identify 28 candidate Class 0 protostars, four of which are newly discovered sources not identified with Spitzer. We find that the star formation efficiency of clumps, as traced by Class 0 protostars, correlates strongly with the flatness of their respective column density distributions at high values. This correlation suggests that the fraction of high column density material in a clump reflects only its youngest protostellar population rather than its entire source population. We propose that feedback from either the formation or evolution of protostars changes the local density structure of clumps.

The whole of the Taurus region (a total area of 52 sq. deg.) has been observed by the Herschel SPIRE and PACS instruments at wavelengths of 70, 160, 250, 350 and 500 {\mu}m as part of the Herschel Gould Belt Survey. In this paper we present the first results from the part of the Taurus region that includes the Barnard 18 and L1536 clouds. A new source-finding routine, the Cardiff Source-finding AlgoRithm (CSAR), is introduced, which is loosely based on CLUMPFIND, but that also generates a structure tree, or dendrogram, which can be used to interpret hierarchical clump structure in a complex region. Sources were extracted from the data using the hierarchical version of CSAR and plotted on a mass-size diagram. We found a hierarchy of objects with sizes in the range 0.024-2.7 pc. Previous studies showed that gravitationally bound prestellar cores and unbound starless clumps appeared in different places on the mass-size diagram. However, it was unclear whether this was due to a lack of instrumental dynamic range or whether they were actually two distinct populations. The excellent sensitivity of Herschel shows that our sources fill the gap in the mass-size plane between starless and pre-stellar cores, and gives the first clear supporting observational evidence for the theory that unbound clumps and (gravitationally bound) prestellar cores are all part of the same population, and hence presumably part of the same evolutionary sequence (c.f. Simpson et al. 2011).

We present Herschel observations from the Herschel Gould Belt Survey and SCUBA-2 science verification observations from the JCMT Gould Belt Survey of the B1 clump in the Perseus molecular cloud. We determined the dust emissivity index using four different techniques to combine the Herschel PACS+SPIRE data at 160 - 500 microns with the SCUBA-2 data at 450 microns and 850 microns. Of our four techniques, we found the most robust method was to filter-out the large-scale emission in the Herschel bands to match the spatial scales recovered by the SCUBA-2 reduction pipeline. Using this method, we find beta ~ 2 towards the filament region and moderately dense material and lower beta values (beta > 1.6) towards the dense protostellar cores, possibly due to dust grain growth. We find that beta and temperature are more robust with the inclusion of the SCUBA-2 data, improving estimates from Herschel data alone by factors of ~ 2 for beta and by ~ 40% for temperature. Furthermore, we find core mass differences of < 30% compared to Herschel-only estimates with an adopted beta = 2, highlighting the necessity of long wavelength submillimeter data for deriving accurate masses of prestellar and protostellar cores.

The W3 GMC is a prime target for the study of the early stages of high-mass star formation. We have used Herschel data from the HOBYS key program to produce and analyze column density and temperature maps. Two preliminary catalogs were produced by extracting sources from the column density map and from Herschel maps convolved to the 500 micron resolution. Herschel reveals that among the compact sources (FWHM<0.45 pc), W3 East, W3 West, and W3 (OH) are the most massive and luminous and have the highest column density. Considering the unique properties of W3 East and W3 West, the only clumps with on-going high-mass star formation, we suggest a 'convergent constructive feedback' scenario to account for the formation of a cluster with decreasing age and increasing system/source mass toward the innermost regions. This process, which relies on feedback by high-mass stars to ensure the availability of material during cluster formation, could also lead to the creation of an environment suitable for the formation of Trapezium-like systems. In common with other scenarios proposed in other HOBYS studies, our results indicate that an active/dynamic process aiding in the accumulation, compression, and confinement of material is a critical feature of the high-mass star/cluster formation, distinguishing it from classical low-mass star formation. The environmental conditions and availability of triggers determine the form in which this process occurs, implying that high-mass star/cluster formation could arise from a range of scenarios: from large scale convergence of turbulent flows, to convergent constructive feedback or mergers of filaments.

We present a study of the star formation histories of the Lupus I, III, and IV clouds using the Herschel 70-500 micron maps obtained by the Herschel Gould Belt Survey Key Project. By combining the new Herschel data with the existing Spitzer catalog we obtained an unprecedented census of prestellar sources and young stellar objects in the Lupus clouds, which allowed us to study the overall star formation rate (SFR) and efficiency (SFE). The high SFE of Lupus III, its decreasing SFR, and its large number of pre-main sequence stars with respect to proto- and prestellar sources, suggest that Lupus III is the most evolved cloud, and after having experienced a major star formation event in the past, is now approaching the end of its current star-forming cycle. Lupus I is currently undergoing a large star formation event, apparent by the increasing SFR, the large number of prestellar objects with respect to more evolved objects, and the high percentage of material at high extinction (e.g., above A_V=8 mag). Also Lupus IV has an increasing SFR; however, the relative number of prestellar sources is much lower, suggesting that its star formation has not yet reached its peak.

We report far-IR Herschel observations obtained between 70 $\mu$m and 500 $\mu$m of two star-forming dusty condensations, B1-bS and B1-bN, in the B1 region of the Perseus star-forming cloud. In the Western part of the Perseus cloud, B1-bS is the only source detected in all of the 6 PACS and SPIRE photometric bands without being visible in the Spitzer map at 24 $\mu$m. B1-bN is clearly detected between 100 $\mu$m and 250 $\mu$m. We have fitted the spectral energy distributions of these sources to derive their physical properties, and find that a simple greybody model fails to reproduce the observed SEDs. At least a two-component model, consisting of a central source surrounded by a dusty envelope is required. The properties derived from the fit, however, suggest that the central source is not a Class 0 object. We then conclude that while B1-bS and B1-bN appear to be more evolved than a pre-stellar core, the best-fit models suggest that their central objects are younger than a Class 0 source. Hence, they may be good candidates to be examples of the first hydrostatic core phase. The projected distance between B1-bS and B1-bN is a few Jeans lengths. If their physical separation is close to this value, this pair would allow the mutual interactions between two forming stars at a very early stage of their evolution to be studied.

We present SMA and CARMA continuum and spectral line observations of five dense cores located in the Perseus and Ophiuchus molecular clouds whose masses exceed their thermal Jeans masses. Three of these cores have previously been identified as being starless and two have been classified as being possibly protostellar. We find that one core is certainly protostellar. The other four cores, however, are starless and undetected in both C18O (2-1) and 1.3mm continuum emission. These four starless cores have flat density profiles out to at least 0.006 pc, which is typical for starless cores in general. Density profiles predicted by some collapse models, especially in the early stages of infall, are consistent with our observations. Archival data reveal that these starless cores have significant non-thermal support against collapse, although they may still be unstable.

We present continuum observations of the Perseus B1-E region from the Herschel Gould Belt Survey. These Herschel data reveal a loose grouping of substructures at 160 - 500 micron not seen in previous submillimetre observations. We measure temperature and column density from these data and select the nine densest and coolest substructures for follow-up spectral line observations with the Green Bank Telescope. We find that the B1-E clump has a mass of ~ 100 solar masses and appears to be gravitationally bound. Furthermore, of the nine substructures examined here, one substructure (B1-E2) appears to be itself bound. The substructures are typically less than a Jeans length from their nearest neighbour and thus, may interact on a timescale of ~ 1 Myr. We propose that B1-E may be forming a first generation of dense cores, which could provide important constraints on the initial conditions of prestellar core formation. Our results suggest that B1-E may be influenced by a strong, localized magnetic field, but further observations are still required.

For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.

The Gould Belt Legacy Survey will map star-forming regions within 500 pc, using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common-User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes HARP observations of the J = 3-2 transitions of 13CO and C18O towards Orion A. The 1500-resolution observations cover 5 pc of the Orion filament, including OMC1 (inc. BN-KL and Orion Bar), OMC 2/3 and OMC 4, and allow a comparative study of the molecular gas properties throughout the star-forming cloud. The filament shows a velocity gradient of ~1 km/s /pc between OMC 1, 2 and 3, and high velocity emission is detected in both isotopologues. The Orion Nebula and Bar have the largest masses and line widths, and dominate the mass and energetics of the high velocity material. Compact, spatially resolved emission from CH3CN, 13CH3OH, SO, HCOOCH3, C2H5OH, CH3CHO and CH3OCHO is detected towards the Orion Hot Core. The cloud is warm, with a median excitation temperature of ~24 K; the Orion Bar has the highest excitation temperature gas, at >80 K. The C18O excitation temperature correlates well with the dust temperature (to within 40%). The C18O emission is optically thin, and the 13CO emission is marginally optically thick; despite its high mass, OMC 1 shows the lowest opacities. A virial analysis indicates that Orion A is too massive for thermal or turbulent support, but is consistent with a model of a filamentary cloud that is threaded by helical magnetic fields. The variation of physical conditions across the cloud is reflected in the physical characteristics of the dust cores....continued

In this paper, we present the results of CARMA continuum and spectral line observations of the dense core Per-Bolo 45. Although this core has previously been classified as starless, we find evidence for an outflow and conclude that Per-Bolo 45 is actually an embedded, low-luminosity protostar. We discuss the impact of newly discovered, low-luminosity, embedded objects in the Perseus molecular cloud on starless core and protostar lifetimes. We estimate that the starless core lifetime has been overestimated by 4-18% and the Class 0/I protostellar lifetime has been underestimated by 5-20%. Given the relatively large systematic uncertainties involved in these calculations, variations on the order of 10% do not significantly change either core lifetimes or the expected protostellar luminosity function. Finally, we suggest that high resolution (sub)millimeter surveys of known cores lacking near-infrared and mid-infrared emission are necessary to make an accurate census of starless cores.

From a survey of 729 cores based on JCMT/SCUBA data, we present an analysis of 17 candidate starless cores with masses that exceed their stable Jeans masses. We re-examine the classification of these super-Jeans cores using Spitzer maps and find that 3 are re-classified as protostellar, 11 have ambiguous emission near the core positions, and 3 appear to be genuinely starless. We suggest the 3 starless and 11 undetermined super-Jeans cores represent excellent targets for future observational and computational study to understand the evolution of dense cores and the process of star formation.

The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region

The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investigating possible gradients of the dust temperature will help to test the 'triggered star-formation' scenario in Rosette. We use continuum data obtained with the PACS (70 and 160 micron) and SPIRE instruments (250, 350, 500 micron) of the Herschel telescope during the Science Demonstration Phase of HOBYS. Three-color images of Rosette impressively show how the molecular gas is heated by the radiative impact of the NGC 2244 cluster. A clear negative temperature gradient and a positive density gradient (running from the HII-region/molecular cloud interface into the cloud) are detected. Studying the spatial distribution of the most massive dense cores (size scale 0.05 to 0.3 pc), we find an age-sequence (from more evolved to younger) with increasing distance to the cluster NGC 2244. No clear gradient is found for the clump (size-scale up to 1 pc) distribution. The existence of temperature and density gradients and the observed age-sequence imply that star formation in Rosette may indeed be influenced by the radiative impact of the central NGC 2244 cluster. A more complete overview of the prestellar and protostellar population in Rosette is required to obtain a firmer result.

We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses (M_C) over the range -0.75 <= log (M_C/M_sun) <= 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups.

Using data from the SCUBA Legacy Catalogue (850 um) and Spitzer Space Telescope (3.6 - 70 um), we explore dense cores in the Ophiuchus, Taurus, Perseus, Serpens, and Orion molecular clouds. We develop a new method to discriminate submillimeter cores found by SCUBA as starless or protostellar, using point source photometry from Spitzer wide field surveys. First, we identify infrared sources with red colors associated with embedded young stellar objects (YSOs). Second, we compare the positions of these YSO-candidates to our submillimeter cores. With these identifications, we construct new, self-consistent starless and protostellar core mass functions (CMFs) for the five clouds. We find best fit slopes to the high-mass end of the CMFs of -1.26 +/- 0.20, -1.22 +/- 0.06, -0.95 +/- 0.20, and -1.67 +/- 0.72 for Ophiuchus, Taurus, Perseus, and Orion, respectively. Broadly, these slopes are each consistent with the -1.35 power-law slope of the Salpeter IMF at higher masses, but suggest some differences. We examine a variety of trends between these CMF shapes and their parent cloud properties, potentially finding a correlation between the high-mass slope and core temperature. We also find a trend between core mass and effective size, but we are very limited by sensitivity. We make similar comparisons between core mass and size with visual extinction (for A_V >= 3) and find no obvious trends. We also predict the numbers and mass distributions of cores that future surveys with SCUBA-2 may detect in each of these clouds.

As part of a JCMT Legacy Survey of star formation in the Gould Belt, we present early science results for Taurus. CO J=3-2 maps have been secured along the north-west ridge and bowl, collectively known as L 1495, along with deep 13CO and C18O J=3-2 maps in two sub-regions. With these data we search for molecular outflows, and use the distribution of flows, HH objects and shocked H2 line emission features, together with the population of young stars, protostellar cores and starless condensations to map star formation across this extensive region. In total 21 outflows are identified. It is clear that the bowl is more evolved than the ridge, harbouring a greater population of T Tauri stars and a more diffuse, more turbulent ambient medium. By comparison, the ridge contains a much younger, less widely distributed population of protostars which, in turn, is associated with a greater number of molecular outflows. We estimate the ratio of the numbers of prestellar to protostellar cores in L 1495 to be ~ 1.3-2.3, and of gravitationally unbound starless cores to (gravitationally bound) prestellar cores to be ~ 1. If we take previous estimates of the protostellar lifetime of ~ 5 x 10^5 yrs, this indicates a prestellar lifetime of 9(+/-3) x 10^5 yrs. From the number of outflows we also crudely estimate the star formation efficiency in L 1495, finding it to be compatible with a canonical value of 10-15 %. We note that molecular outflow-driving sources have redder near-IR colours than their HH jet-driving counterparts. We also find that the smaller, denser cores are associated with the more massive outflows, as one might expect if mass build-up in the flow increases with the collapse and contraction of the protostellar envelope.

The Gould Belt Legacy Survey will survey nearby star-forming regions (within 500 pc), using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common- User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes the initial data obtained using HARP to observe 12CO, 13CO and C18O J = 3 - 2 towards two regions in Orion B, NGC 2024 and NGC 2071. We describe the physical characteristics of the two clouds, calculating temperatures and opacities utilizing all three isotopologues. We find good agreement between temperatures calculated from CO and from dust emission in the dense, energetic regions. We determine the mass and energetics of the clouds, and of the high-velocity material seen in 12CO emission, and compare the relative energetics of the high- and low-velocity material in the two clouds. We present a CLUMPFIND analysis of the 13CO condensations. The slope of the condensation mass functions, at the high-mass ends, is similar to the slope of the initial mass function.

We present observations of probable line-of-sight acceleration of a broad absorption trough of C IV in the quasar SDSS J024221.87+004912.6. We also discuss how the velocity overlap of two other outflowing systems in the same object constrains the properties of the outflows. The Si IV doublet in each system has one unblended transition and one transition which overlaps with absorption from the other system. The residual flux in the overlapping trough is well fit by the product of the residual fluxes in the unblended troughs. For these optically thick systems to yield such a result, at least one of them must consist of individual subunits rather than being a single structure with velocity-dependent coverage of the source. If these subunits are identical, opaque, spherical clouds, we estimate the cloud radius to be r = 3.9 10^15 cm. If they are identical, opaque, linear filaments, we estimate their width to be w = 6.5 10^14 cm. These subunits are observed to cover the Mg II broad emission line region of the quasar, at which distance from the black hole the above filament width is equal to the predicted scale height of the outer atmosphere of a thin accretion disk. Insofar as that scale height is a natural size scale for structures originating in an accretion disk, these observations are evidence that the accretion disk can be a source of quasar absorption systems. Based on data from ESO program 075.B-0190(A).