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

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments.

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.

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.

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.

We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that the droplets follow the same power-law relation between the velocity gradient and size found for larger-scale dense cores. Assuming that rotational motion giving rise to the observed velocity gradient in each core is a solid-body rotation of a rotating body with a uniform density, we derive the "net rotational motions" of the droplets. We find a ratio between rotational and gravitational energies, $\beta$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $\beta \sim 0.039$. We then examine the alignment between the velocity gradient and the major axis of each droplet, using methods adapted from the histogram of relative orientations (HRO) introduced by Soler et al. (2013). We find no definitive correlation between the directions of velocity gradients and the elongations of the cores. Lastly, we discuss physical processes other than rotation that may give rise to the observed velocity field.

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

We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$\mu$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars.

We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.

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.

We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$\mu$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud.

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.

Massive young stellar objects are known to undergo an evolutionary phase in which high mass accretion rates drive strong outflows. A class of objects believed to trace this phase accurately is the GLIMPSE Extended Green Object (EGO) sample, so named for the presence of extended 4.5 $\mu$m emission on sizescales of $\sim$0.1 pc in \textit{Spitzer} images. We have been conducting a multi-wavelength examination of a sample of 12 EGOs with distances of 1 to 5 kpc. In this paper, we present mid-infrared images and photometry of these EGOs obtained with the SOFIA telescope, and subsequently construct SEDs for these sources from the near-IR to sub-millimeter regimes using additional archival data. We compare the results from greybody models and several publicly-available software packages which produce model SEDs in the context of a single massive protostar. The models yield typical \rstar\/$\sim$10 \rsun, \tstar\/$\sim$10$^3$ to 10$^4$ K, and \lstar\/$\sim$1~$-$~40~$\times$~10$^3$ \lsun; the median $L/M$ for our sample is 24.7 \lsun/\msun. Model results rarely converge for \rstar\/ and \tstar, but do for \lstar, which we take to be an indication of the multiplicity and inherently clustered nature of these sources even though, typically, only a single source dominates in the mid-infrared. The median $L/M$ value for the sample suggests that these objects may be in a transitional stage between the commonly described ``IR-quiet'' and ``IR-bright'' stages of MYSO evolution. The median $T_{dust}$ for the sample is less conclusive, but suggests that these objects are either in this transitional stage or occupy the cooler (and presumably younger) part of the IR-bright stage.

We perform a virial analysis of starless dense cores in three nearby star-forming regions : L1688 in Ophiuchus, NGC 1333 in Perseus, and B18 in Taurus. Our analysis takes advantage of comprehensive kinematic information for the dense gas in all of these regions made publicly available through the Green Bank Ammonia Survey Data Release 1, which used to estimate internal support against collapse. We combine this information with ancillary data used to estimate other important properties of the cores, including continuum data from the James Clerk Maxwell Telescope Gould Belt Survey for core identification, core masses, and core sizes. Additionally, we used \textit{Planck} and \textit{Herschel}-based column density maps for external cloud weight pressure, and Five College Radio Astronomy Observatory $^{13}$CO observations for external turbulent pressure. Our self-consistent analysis suggests that many dense cores in all three star-forming regions are not bound by gravity alone, but rather require additional pressure confinement to remain bound. Unlike a recent, similar study in Orion~A, we find that turbulent pressure represents a significant portion of the external pressure budget. Our broad conclusion emphasizing the importance of pressure confinement in dense core evolution, however, agrees with earlier work.

GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from Av~3 to ~20mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of (1-5)x10$^3$ cm$^{-3}$; and ii) the dense phase, located at Av>10mag, with molecular hydrogen densities >10$^4$ cm$^{-3}$. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone where C/H~8x10$^{-5}$ and C/O~0.8-1, until the beginning of the dense phase at Av~10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H~(0.4 - 2.2) x10$^{-6}$, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H~8x10$^{-8}$). Based on our chemical modeling, we constrain the value of $\zeta_{\rm H_2}$ to ~(0.5 - 1.8) x10$^{-16}$ s$^{-1}$ in the translucent cloud.

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

Stars like our Sun form in self-gravitating dense and cold structures within interstellar clouds, called pre-stellar cores. Although much is known about the physical structure of dense clouds just before and soon after the switch-on of a protostar, the central few thousand astronomical units (au) of pre-stellar cores are unexplored. It is within these central regions that stellar systems assemble and fragmentation may take place, with the consequent formation of binaries and multiple systems. We present ALMA Band 6 observations (ACA and 12m array) of the dust continuum emission of the 8 Msun pre-stellar core L1544, with angular resolution of 2'' x 1.6'' (linear resolution 270 au x 216 au). Within the primary beam, a compact region of 0.1 Msun, which we call a "kernel", has been unveiled. The kernel is elongated, with a central flat zone with radius Rker ~ 10'' (~ 1400 au). The average number density within Rker is ~1 x 10^6 cm^{-3}, with possible local density enhancements. The region within Rker appears to have fragmented, but detailed analysis shows that similar substructure can be reproduced by synthetic interferometric observations of a smooth centrally concentrated dense core with a similar central flat zone. The presence of a smooth kernel within a dense core is in agreement with non-ideal magneto-hydro-dynamical simulations of a contracting cloud core with a peak number density of 1 x 10^7 cm^{-3}. Dense cores with lower central densities are completely filtered out when simulated 12m-array observations are carried out. These observations demonstrate that the kernel of dynamically evolved dense cores can be investigated at high angular resolution with ALMA.

We develop a new ''core field structure'' (CFS) model to predict the magnetic field strength and magnetic field fluctuation profile of dense cores using gas kinematics. We use spatially resolved observations of the nonthermal velocity dispersion from the Green Bank Ammonia survey along with column density maps from SCUBA-2 to estimate the magnetic field strength across seven dense cores located in the L1688 region of Ophiuchus. The CFS model predicts the profile of the relative field fluctuation, which is related to the observable dispersion in direction of the polarization vectors. Within the context of our model we find that all the cores have a transcritical mass-to-flux ratio.

We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL region in a field covering 6' x 7'. We find a 0.012 pc wide filament within OMC1, with an aspect ratio of ~37:1, that was missed in previous studies. Its orientation is directly compared to the relative orientation of the magnetic field from the James Clerk Maxwell Telescope BISTRO survey in Orion A. We find a small deviation of ~11 deg between the mean orientation of the filament and the magnetic field, suggesting that they are almost parallel to one another. The filament's column density is estimated to be 2-3 orders of magnitude larger than the filaments studied with Herschel and is possibly self-gravitating given the low values of turbulence found. We further produce maps of the gas kinematics by forward modeling the hyperfine structure of the NH3 (J,K) = (1,1) and (2,2) lines. The resulting distribution of velocity dispersions peaks at ~0.5 km/s, close to the subsonic regime of the gas. This value is about 0.2 km/s smaller than previously measured in single-dish observations of the same region, suggesting that higher angular and spectral resolution observations will identify even lower velocity dispersions that might reach the subsonic turbulence regime in dense gas filaments.

We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both evolved and young regions and weak emission in two evolved regions. HC$_7$N emission is observed only in L1495A-N and L1521D. We find that CCS and HC$_7$N intensity peaks do not coincide with NH$_3$ or dust continuum intensity peaks. We also find that the fractional abundance of CCS does not show a clear correlation with the dynamical evolutionary stage of dense cores. Our findings and chemical modeling indicate that the fractional abundances of CCS and HC$_7$N are sensitive to the initial gas-phase C/O ratio, and they are good tracers of young condensed gas only when the initial C/O is close to solar value. Kinematic analysis using multiple lines including NH$_3$, HC$_7$N, CCS, CO, HCN, \& HCO$^+$ suggests that there may be three different star formation modes in the L1495-B218 filaments. At the hub of the filaments, L1495A/B7N has formed a stellar cluster with large-scale inward flows (fast mode), while L1521D, a core embedded in a filament, is slowly contracting due to its self-gravity (slow mode). There is also one isolated core that appears to be marginally stable and may undergo quasi-static evolution (isolated mode).

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.1 mm dust continuum and CO 2-1 emission toward six dense cores within the Ophiuchus molecular cloud. We detect compact, sub-arcsecond continuum structures toward three targets, two of which (Oph A N6 and SM1) are located in the Ophiuchus A ridge. Two targets, SM1 and GSS 30, contain two compact sources within the ALMA primary beam. We argue that several of the compact structures are small ($R \lesssim 80$ au) accretion disks around young protostars, due to their resolved, elongated structures, coincident radio and x-ray detections, or bipolar outflow detections. While CO line wings extend to $\pm 10-20$ km s$^{-1}$ for the more evolved sources GSS 30 IRS3 and IRS1, CO emission toward other sources, where detected, only extends a few km s$^{-1}$ from the cloud $v_\mathrm{LSR}$. The dust spectral index toward the compact objects suggests that the disks are either optically thick at 1.1 mm, or that significant grain growth has already occurred. We identify, for the first time, a single compact continuum source ($R \sim 100$ au) toward N6 embedded within a larger continuum structure. SM1N is extended in the continuum but is highly centrally concentrated, with a density profile that follows a $r^{-1.3}$ power law within 200 au, and additional structure suggested by the uv-data. Both N6 and SM1N show no clear bipolar outflows with velocities greater than a few km s$^{-1}$ from the cloud velocity. These sources are candidates to be the youngest protostars or first hydrostatic cores in the Ophiuchus molecular cloud.

In dense starless and protostellar cores, the relative abundance of deuterated species to their non-deuterated counterparts can become orders of magnitude greater than in the local interstellar medium. This enhancement proceeds through multiple pathways in the gas phase and on dust grains, where the chemistry is strongly dependent on the physical conditions. In this Chapter, we discuss how sensitive, high resolution observations with the ngVLA of emission from deuterated molecules will trace both the dense gas structure and kinematics on the compact physical scales required to track the gravitational collapse of star-forming cores and the subsequent formation of young protostars and circumstellar accretion regions. Simultaneously, such observations will play a critical role in tracing the chemical history throughout the various phases of star and planet formation. Many low-J transitions of key deuterated species, along with their undeuterated counterparts, lie within the 60-110 GHz frequency window, the lower end of which is largely unavailable with current facilities and instrumentation. The combination of sensitivity and angular resolution provided only by the ngVLA will enable unparalleled detailed studies of the physics and chemistry of the earliest stages of star formation.

Continuum observations of molecular clouds have revealed a surprising amount of substructure in the form of filaments of a few pc length and cores of ~0.1 pc diameter. Understanding the evolution of these substructures towards star formation requires the kinematic and dynamical insights provided uniquely by sensitive line observations at high angular and spectral resolution. In this short paper, we describe how an ngVLA can probe effectively the dynamics of filaments and cores in nearby star-forming molecular clouds using the NH3 rotation-inversion transitions at 24 GHz. Such emission has been proven to trace well the high column density environments of star-forming cores and filaments but higher-resolution observations are needed to reveal important details of how dense gas is flowing within and onto these substructures. In particular, we describe how 150 x 18-m antennas with a maximum baseline of 1 km can be used to map sensitively NH3 emission across high column density locations in clouds in roughly an order of magnitude less time than with the current Jansky VLA.

On the basis of various data such as ALMA, JVLA, Chandra, {\it Herschel}, and {\it Spitzer}, we confirmed that two protostellar candidates in Oph-A are bona fide protostars or proto-brown dwarfs (proto-BDs) in extremely early evolutionary stages. Both objects are barely visible across infrared (IR, i.e., near-IR to far-IR) bands. The physical nature of the cores is very similar to that expected in first hydrostatic cores (FHSCs), objects theoretically predicted in the evolutionary phase prior to stellar core formation with gas densities of $\sim$ 10$^{11-12}$ cm$^{-3}$. This suggests that the evolutionary stage is close to the FHSC formation phase. The two objects are associated with faint X-ray sources, suggesting that they are in very early phase of stellar core formation with magnetic activity. In addition, we found the CO outflow components around both sources which may originate from the young outflows driven by these sources. The masses of these objects are calculated to be $\sim 0.01-0.03$ $M_\odot$ from the dust continuum emission. Their physical properties are consistent with that expected from the numerical model of forming brown dwarfs. These facts (the X-ray detection, CO outflow association, and FHSC-like spectral energy distributions) strongly indicate that the two objects are proto-BDs or will be in the very early phase of protostars which will evolve more massive protostars if they gain enough mass from the surroundings. The ages of these two objects are likely to be within $\sim 10^3$ years after the protostellar core (or second core) formation, taking into account the outflow dynamical times ($\lesssim$ 500 yrs).

We compare multi-epoch sub-arcsecond VLA imaging of the 22 GHz water masers toward the massive protocluster NGC6334I observed before and after the recent outburst of MM1B in (sub)millimeter continuum. Since the outburst, the water maser emission toward MM1 has substantially weakened. Simultaneously, the strong water masers associated with the synchrotron continuum point source CM2 have flared by a mean factor of 6.5 (to 4.2 kJy) with highly-blueshifted features (up to 70 km/s from LSR) becoming more prominent. The strongest flaring water masers reside 3000 au north of MM1B and form a remarkable bow shock pattern whose vertex coincides with CM2 and tail points back to MM1B. Excited OH masers trace a secondary bow shock located ~120 au downstream. ALMA images of CS (6-5) reveal a highly-collimated north-south structure encompassing the flaring masers to the north and the non-flaring masers to the south seen in projection toward the MM3-UCHII region. Proper motions of the southern water masers over 5.3 years indicate a bulk projected motion of 117 km/s southward from MM1B with a dynamical time of 170 yr. We conclude that CM2, the water masers, and many of the excited OH masers trace the interaction of the high velocity bipolar outflow from MM1B with ambient molecular gas. The previously-excavated outflow cavity has apparently allowed the radiative energy of the current outburst to propagate freely until terminating at the northern bow shock where it strengthened the masers. Additionally, water masers have been detected toward MM7 for the first time, and a highly-collimated CS (6-5) outflow has been detected toward MM4.

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

We present 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 report the first sub-arcsecond VLA imaging of 6 GHz continuum, methanol maser, and excited-state hydroxyl maser emission toward the massive protostellar cluster NGC6334I following the recent 2015 outburst in (sub)millimeter continuum toward MM1, the strongest (sub)millimeter source in the protocluster. In addition to detections toward the previously known 6.7 GHz Class II methanol maser sites in the hot core MM2 and the UCHII region MM3 (NGC6334F), we find new maser features toward several components of MM1, along with weaker features $\sim1''$ north, west, and southwest of MM1, and toward the non-thermal radio continuum source CM2. None of these areas have heretofore exhibited Class II methanol maser emission in three decades of observations. The strongest MM1 masers trace a dust cavity, while no masers are seen toward the strongest dust sources MM1A, 1B and 1D. The locations of the masers are consistent with a combination of increased radiative pumping due to elevated dust grain temperature following the outburst, the presence of infrared photon propagation cavities, and the presence of high methanol column densities as indicated by ALMA images of thermal transitions. The non-thermal radio emission source CM2 ($2''$ north of MM1) also exhibits new maser emission from the excited 6.035 and 6.030 GHz OH lines. Using the Zeeman effect, we measure a line-of-sight magnetic field of +0.5 to +3.7 mG toward CM2. In agreement with previous studies, we also detect numerous methanol and excited OH maser spots toward the UCHII region MM3, with predominantly negative line-of-sight magnetic field strengths of -2 to -5 mG and an intriguing south-north field reversal.

Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency.

We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.

We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.

We present $\sim1-4"$ resolution Very Large Array (VLA) observations of four CH$_3$OH $J_2-J_1$-$E$ 25~GHz transitions ($J$=3, 5, 8, 10) along with 1.3~cm continuum toward 20 regions of active massive star formation containing Extended Green Objects (EGOs), 14 of which we have previously studied with the VLA in the Class~I 44~GHz and Class~II 6.7~GHz maser lines (Cyganowski et al. 2009). Sixteen regions are detected in at least one 25~GHz line ($J$=5), with 13 of 16 exhibiting maser emission. In total, we report 34 new sites of CH$_3$OH maser emission and ten new sites of thermal CH$_3$OH emission, significantly increasing the number of 25~GHz Class I CH$_3$OH masers observed at high angular resolution. We identify probable or likely maser counterparts at 44~GHz for all 15 of the 25~GHz masers for which we have complementary data, providing further evidence that these masers trace similar physical conditions despite uncorrelated flux densities. The sites of thermal and maser emission of CH$_3$OH are both predominantly associated with the 4.5 $\mu$m emission from the EGO, and the presence of thermal CH$_3$OH emission is accompanied by 1.3~cm continuum emission in 9 out of 10 cases. Of the 19 regions that exhibit 1.3~cm continuum emission, it is associated with the EGO in 16 cases (out of a total of 20 sites), 13 of which are new detections at 1.3~cm. Twelve of the 1.3~cm continuum sources are associated with 6.7~GHz maser emission and likely trace deeply-embedded massive protostars.

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 an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.

Based on sub-arcsecond Atacama Large Millimeter/submillimeter Array (ALMA) and Submillimeter Array (SMA) 1.3 mm continuum images of the massive protocluster NGC 6334I obtained in 2015 and 2008, we find that the dust emission from MM1 has increased by a factor of 4.0$\pm$0.3 during the intervening years, and undergone a significant change in morphology. The continuum emission from the other cluster members (MM2, MM4 and the UCHII region MM3 = NGC 6334F) has remained constant. Long term single-dish maser monitoring at HartRAO finds that multiple maser species toward NGC 6334I flared beginning in early 2015, a few months before our ALMA observation, and some persist in that state. New ALMA images obtained in 2016 July-August at 1.1 and 0.87 mm confirm the changes with respect to SMA 0.87 mm images from 2008, and indicate that the (sub)millimeter flaring has continued for at least a year. The excess continuum emission, centered on the hypercompact HII region MM1B, is extended and elongated ($1.6" \times 1.0" \approx 2100 \times 1300$ au) with multiple peaks, suggestive of general heating of the surrounding subcomponents of MM1, some of which may trace clumps in a fragmented disk rather than separate protostars. In either case, these remarkable increases in maser and dust emission provide direct observational evidence of a sudden accretion event in the growth of a massive protostar yielding a sustained luminosity surge by a factor of $70\pm20$, analogous to the largest events in simulations by Meyer et al. (2017). This target provides an excellent opportunity to assess the impact of such a rare event on a protocluster over many years.

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.

Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.

Filamentary structures are ubiquitous in molecular clouds, and have been recently argued to play an important role in regulating the size and mass of embedded clumps through fragmentation and mass accretion. Here, we reveal the dynamical state and fragmentation of filamentary molecular gas associated with the Serpens South protocluster through analysis of wide (~4 x 4 pc) observations of NH3 (1,1) and (2,2) inversion transitions with the Green Bank Telescope. Detailed modeling of the NH3 lines reveals that the kinematics of the cluster and surrounding filaments are complex. We identify hierarchical structure using a dendrogram analysis of the NH3 emission. The distance between neighbour structures that are embedded within the same parent structure is generally greater than expected from a spherical Jeans analysis, and is in better agreement with cylindrical fragmentation models. The NH3 line width-size relation is flat, and average gas motions are sub- or trans-sonic over all physical scales observed. Subsonic regions extend far beyond the typical 0.1 pc scale previously identified in star-forming cores. As a result, we find a strong trend of decreasing virial parameter with increasing structure mass in Serpens South. Extremely low virial parameters on the largest scales probed by our data suggest that the previously observed, ordered magnetic field is insufficient to support the region against collapse, in agreement with large radial infall motions previously measured toward some of the filaments. A more complex magnetic field configuration in the dense gas, however, may be able to support the filaments.

Although surveys of infall motions in dense cores have been carried out for years, few surveys have focused on mapping infall across cores using multiple spectral line observations. To fill this gap, we present IRAM 30-m Telescope maps of N2H+(1-0), DCO+(2-1), DCO+(3-2), and HCO+(3-2) emission towards two prestellar cores (L492 and L694-2) and one protostellar core (L1521F). We find that the measured infall velocity varies with position across each core and choice of molecular line, likely as a result of radial variations in core chemistry and dynamics. Line-of-sight infall speeds estimated from DCO+(2-1) line profiles can decrease by 40-50 m/s when observing at a radial offset >= 0.04 pc from the core's dust continuum emission peak. Median infall speeds calculated from all observed positions across a core can also vary by as much as 65 m/s depending on the transition. These results show that while single-pointing, single-transition surveys of core infall velocities may be good indicators of whether a core is either contracting or expanding, the magnitude of the velocities they measure are significantly impacted by the choice of molecular line, proximity to the core center, and core evolutionary state.

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.

We present VLA and ALMA imaging of the deeply-embedded protostellar cluster NGC6334I from 5 cm to 1.3 mm at angular resolutions as fine as 0.17 (220 AU). The dominant hot core MM1 is resolved into seven components at 1.3 mm, clustered within a radius of 1000 AU. Four of the components have brightness temperatures >200 K, radii ~300 AU, minimum luminosities ~10000 Lsun, and must be centrally heated. We term this new phenomenon a "hot multi-core". Two of these objects also exhibit compact free-free emission at longer wavelengths, consistent with a hypercompact HII region (MM1B) and a jet (MM1D). The spatial kinematics of the water maser emission centered on MM1D are consistent with it being the origin of the high-velocity bipolar molecular outflow seen in CO. The close proximity of MM1B and MM1D (440~AU) suggests a proto-binary or a transient bound system. Several components of MM1 exhibit steep millimeter SEDs indicative of either unusual dust spectral properties or time variability. In addition to resolving MM1 and the other hot core (MM2) into multiple components, we detect five new millimeter and two new centimeter sources. Water masers are detected for the first time toward MM4A, confirming its membership in the protocluster. With a 1.3 mm brightness temperature of 97 K coupled with a lack of thermal molecular line emission, MM4A appears to be a highly optically-thick 240 Lsun dust core, possibly tracing a transient stage of massive protostellar evolution. The nature of the strongest water maser source CM2 remains unclear due to its combination of non-thermal radio continuum and lack of dust emission.

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.

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 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 report the results of an HCO+ (3-2) and N2D+ (3-2) molecular line survey performed toward 91 dense cores in the Perseus molecular cloud using the James Clerk Maxwell Telescope, to identify the fraction of starless and protostellar cores with systematic radial motions. We quantify the HCO+ asymmetry using a dimensionless asymmetry parameter $\delta_v$, and identify 20 cores with significant blue or red line asymmetries in optically-thick emission indicative of collapsing or expanding motions, respectively. We separately fit the HCO+ profiles with an analytic collapse model and determine contraction (expansion) speeds toward 22 cores. Comparing the $\delta_v$ and collapse model results, we find that $\delta_v$ is a good tracer of core contraction if the optically-thin emission is aligned with the model-derived systemic velocity. The contraction speeds range from subsonic (0.03 km/s) to supersonic (0.4 km/s), where the supersonic contraction speeds may trace global rather than local core contraction. Most cores have contraction speeds significantly less than their free-fall speeds. Only 7 of 28 starless cores have spectra well-fit by the collapse model, which more than doubles (15 of 28) for protostellar cores. Starless cores with masses greater than the Jeans mass (M/M$_J$ > 1) are somewhat more likely to show contraction motions. We find no trend of optically-thin non-thermal line width with M/M$_J$, suggesting that any undetected contraction motions are small and subsonic. Most starless cores in Perseus are either not in a state of collapse or expansion, or are in a very early stage of collapse.

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.

Infrared Dark Clouds (IRDCs) harbor the earliest phases of massive star formation, and many of the compact cores in IRDCs, traced by millimeter continuum or by molecular emission in high critical density lines, host massive young stellar objects (YSOs). We used the Robert C. Byrd Green Bank Telescope (GBT) and the Karl G. Jansky Very Large Array (VLA) to map NH$_{3}$ and CCS in nine IRDCs to reveal the temperature, density, and velocity structures and explore chemical evolution in the dense ($>10^{22}$ cm$^{-2}$) gas. Ammonia is an excellent molecular tracer for these cold, dense environments. The internal structure and kinematics of the IRDCs include velocity gradients, filaments, and possibly colliding clumps that elucidate the formation process of these structures and their YSOs. We find a wide variety of substructure including filaments and globules at distinct velocities, sometimes overlapping at sites of ongoing star formation. It appears that these IRDCs are still being assembled from molecular gas clumps even as star formation has already begun, and at least three of them appear consistent with the morphology of ``hub-filament structures'' discussed in the literature. Furthermore, we find that these clumps are typically near equipartition between gravitational and kinetic energies, so these structures may survive for multiple free-fall times.Keywords: molecular data -- ISM: clouds -- (ISM:) dust, extinction -- ISM: molecules -- Stars: formation -- radio lines: ISM

We present deep NH$_3$ observations of the L1495-B218 filaments in the Taurus molecular cloud covering over a 3 degree angular range using the K-band focal plane array on the 100m Green Bank Telescope. The L1495-B218 filaments form an interconnected, nearby, large complex extending over 8 pc. We observed NH$_3$ (1,1) and (2,2) with a spectral resolution of 0.038 km/s and a spatial resolution of 31$"$. Most of the ammonia peaks coincide with intensity peaks in dust continuum maps at 350 $\mu$m and 500 $\mu$m. We deduced physical properties by fitting a model to the observed spectra. We find gas kinetic temperatures of 8 $-$ 15 K, velocity dispersions of 0.05 $-$ 0.25 km/s, and NH$_3$ column densities of 5$\times$10$^{12}$ $-$ 1$\times$10$^{14}$ cm$^{-2}$. The CSAR algorithm, which is a hybrid of seeded-watershed and binary dendrogram algorithms, identifies a total of 55 NH$_3$ structures including 39 leaves and 16 branches. The masses of the NH$_3$ sources range from 0.05 M$_\odot$ to 9.5 M$_\odot$. The masses of NH$_3$ leaves are mostly smaller than their corresponding virial mass estimated from their internal and gravitational energies, which suggests these leaves are gravitationally unbound structures. 9 out of 39 NH$_3$ leaves are gravitationally bound and 7 out of 9 gravitationally bound NH$_3$ leaves are associated with star formation. We also found that 12 out of 30 gravitationally unbound leaves are pressure-confined. Our data suggest that a dense core may form as a pressure-confined structure, evolve to a gravitationally bound core, and undergo collapse to form a protostar.

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 present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the submillimeter dust continuum and H2D+ 1_{10}-1_{11} emission toward two evolved, potentially protostellar cores within the Ophiuchus molecular cloud, Oph A SM1 and SM1N. The data reveal small-scale condensations within both cores, with mass upper limits of M <~ 0.02M_Sun (~ 20 M_Jup). The SM1 condensation is consistent with a nearly-symmetric Gaussian source with a width of only 37 AU. The SM1N condensation is elongated, and extends 500 AU along its major axis. No evidence for substructure is seen in either source. A Jeans analysis indicates these sources are unlikely to fragment, suggesting that both will form single stars. H2D+ is only detected toward SM1N, offset from the continuum peak by ~150-200 AU. This offset may be due to either heating from an undetected, young, low luminosity protostellar source or first hydrostatic core, or HD (and consequently H2D+) depletion in the cold centre of the condensation. We propose that SM1 is protostellar, and that the condensation detected by ALMA is a warm (T ~ 30-50 K) accretion disk. The less concentrated emission of the SM1N condensation suggests that it is still starless, but we cannot rule out the presence of a low-luminosity source, perhaps surrounded by a pseudodisk. These data reveal observationally the earliest stages of the formation of circumstellar accretion regions, and agree with theoretical predictions that disk formation can occur very early in the star formation process, coeval with or just after the formation of a first hydrostatic core or protostar.

Core accretion models of massive star formation require the existence of stable massive starless cores, but robust observational examples of such objects have proven elusive. We report subarcsecond-resolution SMA 1.3 mm, 1.1 mm, and 0.88 mm and VLA 1.3 cm observations of an excellent massive starless core candidate, G11.92-0.61-MM2, initially identified in the course of studies of GLIMPSE Extended Green Objects (EGOs). Separated by ~7.2" from the nearby MM1 protostellar hot core, MM2 is a strong, compact dust continuum source (submillimeter spectral index alpha=2.6+/-0.1), but is devoid of star formation indicators. In contrast to MM1, MM2 has no masers, no centimeter continuum, and no (sub)millimeter wavelength line emission in ~24 GHz of bandwidth observed with the SMA, including N2H+(3-2), HCO+(3-2), and HCN(3-2). Additionally, there is no evidence for an outflow driven by MM2. The (sub)millimeter spectral energy distribution (SED) of MM2 is best fit with a dust temperature of ~17-19 K and luminosity of ~5-7 L_sun. The combined physical properties of MM2, as inferred from its dust continuum emission, are extreme: M>30 M_sun within a radius<1000 AU, N(H2)>10^25 cm^-2 and n(H2)>10^9 cm^-3. Comparison of the molecular abundance limits derived from our SMA observations with gas-grain chemical models indicates that extremely dense (n(H)>>10^8 cm^-3), cold (<20 K) conditions are required to explain the lack of observed (sub)millimeter line emission, consistent with the dust continuum results. Our data suggest that G11.92-0.61-MM2 is the best candidate for a bonafide massive prestellar core found to date, and a promising target for future, higher-sensitivity observations.

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 present a survey of HCO+ (3-2) observations pointed towards dense cores with previous measurements of N(N2D+)/N(N2H+). Of the 26 cores in this survey, five show the spectroscopic signature of outward motion, nine exhibit neither inward nor outward motion, eleven appear to be infalling, and one is not detected. We compare the degree of deuterium fractionation with infall velocities calculated from the HCO+ spectra and find that those cores with [D]/[H] > 0.1 are more likely to have the signature of inward motions than cores with smaller [D]/[H] ratios. Infall motions are also much more common in cores with masses exceeding their thermal Jeans masses. The fastest infall velocity measured belongs to one of the two protostellar cores in our survey, L1521F, and the observed motions are typically on the order of the sound speed.

We have detected bright HC7N J = 21-20 emission toward multiple locations in the Serpens South cluster-forming region using the K-Band Focal Plane Array at the Robert C. Byrd Green Bank Telescope. HC7N is seen primarily toward cold filamentary structures that have yet to form stars, largely avoiding the dense gas associated with small protostellar groups and the main central cluster of Serpens South. Where detected, the HC7N abundances are similar to those found in other nearby star forming regions. Toward some HC7N `clumps', we find consistent variations in the line centroids relative to NH3 (1,1) emission, as well as systematic increases in the HC7N non-thermal line widths, which we argue reveal infall motions onto dense filaments within Serpens South with minimum mass accretion rates of M ~ 2-5 M_sun Myr^-1. The relative abundance of NH3 to HC7N suggests that the HC7N is tracing gas that has been at densities n ~ 10^4 cm^-3, for timescales t < 1-2 x 10^5 yr. Since HC7N emission peaks are rarely co-located with those of either NH3 or continuum, it is likely that Serpens South is not particularly remarkable in its abundance of HC7N, but instead the serendipitous mapping of HC7N simultaneously with NH3 has allowed us to detect HC7N at low abundances in regions where it otherwise may not have been looked for. This result extends the known star-forming regions containing significant HC7N emission from typically quiescent regions, like the Taurus molecular cloud, to more complex, active environments.

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.

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

We have performed a pointed survey of N2D+ 2-1 and N2D+ 3-2 emission toward 64 N2H+-bright starless and protostellar cores in the Perseus molecular cloud using the Arizona Radio Observatory Submillimeter Telescope and Kitt Peak 12 m telescope. We find a mean deuterium fractionation in N2H+, R_D = N(N2D+)/N(N2H+), of 0.08, with a maximum R_D = 0.2. In detected sources, we find no significant difference in the deuterium fractionation between starless and protostellar cores, nor between cores in clustered or isolated environments. We compare the deuterium fraction in N2H+ with parameters linked to advanced core evolution. We only find significant correlations between the deuterium fraction and increased H_2 column density, as well as with increased central core density, for all cores. Towards protostellar sources, we additionally find a significant anti-correlation between R_D and bolometric temperature. We show that the Perseus cores are characterized by low CO depletion values relative to previous studies of star forming cores, similar to recent results in the Ophiuchus molecular cloud. We suggest that the low average CO depletion is the dominant mechanism that constrains the average deuterium fractionation in the Perseus cores to small values. While current equilibrium and dynamic chemical models are able to reproduce the range of deuterium fractionation values we find in Perseus, reproducing the scatter across the cores requires variation in parameters such as the ionization fraction or the ortho- to para-H_2 ratio across the cloud, or a range in core evolution timescales.

We investigated six H II regions with infrared, bright rimmed bubble or cometary morphology, in search of quantitative evidence for triggered star formation, both collect and collapse and radiatively driven implosion. We identified and classified 458 Young Stellar Objects (YSOs) in and around the H II regions. YSOs were determined by fitting a collection of radiative transfer model spectral energy distributions (SEDs) to infrared photometry for a large sample of point sources. We determined areas where there exist enhanced populations of relatively unevolved YSOs on the bright rims of these regions, suggesting that star formation has been triggered there. We further investigated the physical properties of the regions by using radio continuum emission as a proxy for ionizing flux powering the H II regions, and 13CO (1-0) observations to measure masses and gravitational stability of molecular clumps. We used an analytical model of collect and collapse triggered star formation, as well as a simulation of radiatively driven implosion, and thus we compare the observed properties of the molecular gas with those predicted in the triggering scenarios. Notably, those regions in our sample that show evidence of cometary, or "blister," morphology are more likely to show evidence of triggering. Keywords: HII regions -- ISM: bubbles -- Stars: formation -- Stars: protostars

We present high angular resolution Submillimeter Array (SMA) and Karl G. Jansky Very Large Array (VLA) observations of the massive protocluster G18.67+0.03. Previously targeted in maser surveys of GLIMPSE Extended Green Objects (EGOs), this cluster contains three Class I methanol maser sources, providing a unique opportunity to test the proposed role of Class I masers as evolutionary indicators for massive star formation. The millimeter observations reveal bipolar molecular outflows, traced by 13CO(2-1) emission, associated with all three Class I maser sources. Two of these sources (including the EGO) are also associated with 6.7 GHz Class II methanol masers; the Class II masers are coincident with millimeter continuum cores that exhibit hot core line emission and drive active outflows, as indicated by the detection of SiO(5-4). In these cases, the Class I masers are coincident with outflow lobes, and appear as clear cases of excitation by active outflows. In contrast, the third Class I source is associated with an ultracompact HII region, and not with Class II masers. The lack of SiO emission suggests the 13CO outflow is a relic, consistent with its longer dynamical timescale. Our data show that massive young stellar objects associated only with Class I masers are not necessarily young, and provide the first unambiguous evidence that Class I masers may be excited by both young (hot core) and older (UC HII) MYSOs within the same protocluster.

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.

We have performed a 1.3 centimeter survey of 24 massive young stellar objects (MYSOs) using the Expanded Very Large Array (EVLA). The sources in the sample exhibit a broad range of massive star formation signposts including Infrared Dark Clouds (IRDCs), UCHII regions, and extended 4.5 micron emission in the form of Extended Green Objects (EGOs). In this work, we present results for G35.03+0.35 which exhibits all of these phenomena. We simultaneously image the 1.3 cm ammonia (1,1) through (6,6) inversion lines, four methanol transitions, two H recombination lines, plus continuum at 0.05 pc resolution. We find three areas of thermal ammonia emission, two within the EGO (designated the NE and SW cores) and one toward an adjacent IRDC. The NE core contains an UCHII region (CM1) and a candidate HCHII region (CM2). A region of non-thermal, likely masing ammonia (3,3) and (6,6) emission is coincident with an arc of 44 GHz methanol masers. We also detect two new 25 GHz Class I methanol masers. A complementary Submillimeter Array 1.3 mm continuum image shows that the distribution of dust emission is similar to the lower-lying ammonia lines, all peaking to the NW of CM2, indicating the likely presence of an additional MYSO in this protocluster. By modeling the ammonia and 1.3 mm continuum data, we obtain gas temperatures of 20-220 K and masses of 20-130 solar. The diversity of continuum emission properties and gas temperatures suggest that objects in a range of evolutionary states exist concurrently in this protocluster.

We present N2D+ 3-2 (IRAM) and H2D+ 1_11 - 1_10 and N2H+ 4-3 (JCMT) maps of the small cluster-forming Ophiuchus B2 core in the nearby Ophiuchus molecular cloud. In conjunction with previously published N2H+ 1-0 observations, the N2D+ data reveal the deuterium fractionation in the high density gas across Oph B2. The average deuterium fractionation R_D = N(N2D+)/N(N2H+) ~ 0.03 over Oph B2, with several small scale R_D peaks and a maximum R_D = 0.1. The mean R_D is consistent with previous results in isolated starless and protostellar cores. The column density distributions of both H2D+ and N2D+ show no correlation with total H2 column density. We find, however, an anticorrelation in deuterium fractionation with proximity to the embedded protostars in Oph B2 to distances >= 0.04 pc. Destruction mechanisms for deuterated molecules require gas temperatures greater than those previously determined through NH3 observations of Oph B2 to proceed. We present temperatures calculated for the dense core gas through the equating of non-thermal line widths for molecules (i.e., N2D+ and H2D+) expected to trace the same core regions, but the observed complex line structures in B2 preclude finding a reasonable result in many locations. This method may, however, work well in isolated cores with less complicated velocity structures. Finally, we use R_D and the H2D+ column density across Oph B2 to set a lower limit on the ionization fraction across the core, finding a mean x_e, lim >= few x 10^{-8}. Our results show that care must be taken when using deuterated species as a probe of the physical conditions of dense gas in star-forming regions.

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

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.

We present a Nobeyama 45 m Radio Telescope map and Australia Telescope Compact Array pointed observations of N2H+ 1-0 emission towards the clustered, low mass star forming Oph B Core within the Ophiuchus molecular cloud. We compare these data with previously published results of high resolution NH3 (1,1) and (2,2) observations in Oph B. We use 3D Clumpfind to identify emission features in the single-dish N2H+ map, and find that the N2H+ `clumps' match well similar features previously identified in NH3 (1,1) emission, but are frequently offset to clumps identified at similar resolution in 850 micron continuum emission. Wide line widths in the Oph B2 sub-Core indicate non-thermal motions dominate the Core kinematics, and remain transonic at densities n ~ 3 x 10^5 cm^-3 with large scatter and no trend with N(H2). Non-thermal motions in Oph B1 and B3 are subsonic with little variation, but also show no trend with H2 column density. Over all Oph B, non-thermal N2H+ line widths are substantially narrower than those traced by NH3, making it unlikely NH3 and N2H+ trace the same material, but the v_LSR of both species agree well. We find evidence for accretion in Oph B1 from the surrounding ambient gas. The NH3/N2H+ abundance ratio is larger towards starless Oph B1 than towards protostellar Oph B2, similar to recent observational results in other star-forming regions. Small-scale structure is found in the ATCA N2H+ 1-0 emission, where emission peaks are again offset from continuum emission. In particular, the ~1 M_Sun B2-MM8 clump is associated with a N2H+ emission minimum and surrounded by a broken ring-like N2H+ emission structure, suggestive of N2H+ depletion. We find a strong general trend of decreasing N2H+ abundance with increasing N(H2) in Oph B which matches that found for NH3.

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 have observed the Oph A-N6 prestellar core in the following transitions: N2D+ J=3 to 2, DCO+ J=3 to 2 and J=5 to 4, HCO+ J=3 to 2, CS J=5 to 4 and J=7 to 6, and H13CO+ J=3 to 2 and J=4 to 3, using the James Clerk Maxwell Telescope. We also observed the NH3 (1,1) and (2,2) inversion transitions towards the Oph A-N6 peak with the Green Bank Telescope. We have found that the N6 core is composed of shells of different chemical composition due to the freezing out of chemical species at different densities. The undepleted species N2D+ appears to trace the high-density interior of the core, DCO+ and H13CO+ trace an intermediate region, and CS traces the outermost edges of the core. A distinct blue-red spectral asymmetry, indicative of infall motion, is clearly detected in the HCO+ spectra, suggesting that N6 is undergoing gravitational collapse. This collapse was possibly initiated by a decrease in turbulent support suggested by the fact that the non-thermal contribution to the line widths is smaller for the molecular species closer to the center of the core. We also present a temperature profile for the core. These observations support the claim that the Oph A-N6 core is an extremely young prestellar core, which may have been recently cut off from MHD support and begun to collapse.

We present combined interferometer and single dish telescope data of NH3 (J,K) = (1,1) and (2,2) emission towards the clustered star forming Ophiuchus B, C and F Cores at high spatial resolution (~1200 AU) using the Australia Telescope Compact Array, the Very Large Array, and the Green Bank Telescope. While the large scale features of the NH3 (1,1) integrated intensity appear similar to 850 micron continuum emission maps of the Cores, on 15" (1800 AU) scales we find significant discrepancies between the dense gas tracers in Oph B, but good correspondence in Oph C and F. Using the Clumpfind structure identifying algorithm, we identify 15 NH3 clumps in Oph B, and 3 each in Oph C and F. Only five of the Oph B NH3 clumps are coincident within 30" (3600 AU) of a submillimeter clump. We find v_LSR varies little across any of the Cores, and additionally varies by only ~1.5 km/s between them. The observed NH3 line widths within the Oph B and F Cores are generally large and often mildly supersonic, while Oph C is characterized by narrow line widths which decrease to nearly thermal values. We find several regions of localized narrow line emission (\Delta v < 0.4 km/s), some of which are associated with NH3 clumps. We derive the kinetic temperatures of the gas, and find they are remarkably constant across Oph B and F, with a warmer mean value (T_K = 15 K) than typically found in isolated regions and consistent with previous results in clustered regions. Oph C, however, has a mean T_K = 12 K, decreasing to a minimum T_K = 9.4 K towards the submillimeter continuum peak, similar to previous studies of isolated starless cores. There is no significant difference in temperature towards protostars embedded in the Cores. [Abridged]

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

We have used a Fourier Transform Spectrometer on the James Clerk Maxwell Telescope to study the submillimetre continuum emission from dust in three hot molecular cores (HMC). The spectral index beta of the dust emission for these sources has been determined solely within the 30 GHz wide 350 GHz (850 micron) passband to an accuracy comparable to those determined through multi-wavelength observations. We find an average beta = 1.6, in agreement with spectral indices determined from previous submillimetre observations of these sources and with those determined for HMC in general. The largest single source of uncertainty in these results is the variability of the atmosphere at 350 GHz, and with better sky subtraction techniques we show that the dust spectral index can clearly be determined within one passband to high accuracy with a submillimetre FTS. Using an imaging FTS on SCUBA-2, the next generation wide-field submillimetre camera currently under development to replace SCUBA at the JCMT in 2006, we calculate that at 350 GHz it will be possible to determine beta to +/- 0.1 for sources as faint as 400 mJy/beam and to +/- 0.3 for sources as faint as 140 mJy/beam.