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

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

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

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.

The JCMT Gould Belt Legacy Survey obtained SCUBA-2 observations of dense cores within three sub-regions of Orion B: LDN 1622, NGC 2023/2024, and NGC 2068/2071, all of which contain clusters of cores. We present an analysis of the clustering properties of these cores, including the two-point correlation function and Cartwright's Q parameter. We identify individual clusters of dense cores across all three regions using a minimal spanning tree technique, and find that in each cluster, the most massive cores tend to be centrally located. We also apply the independent M-Sigma technique and find a strong correlation between core mass and the local surface density of cores. These two lines of evidence jointly suggest that some amount of mass segregation in clusters has happened already at the dense core stage.

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.

Performing ground-based submillimetre observations is a difficult task as the measurements are subject to absorption and emission from water vapour in the Earth's atmosphere and time variation in weather and instrument stability. Removing these features and other artifacts from the data is a vital process which affects the characteristics of the recovered astronomical structure we seek to study. In this paper, we explore two data reduction methods for data taken with the Submillimetre Common-User Bolometer Array-2 (SCUBA-2) at the James Clerk Maxwell Telescope (JCMT). The JCMT Legacy Reduction 1 (JCMT LR1) and The Gould Belt Legacy Survey Legacy Release 1 (GBS LR1) reduction both use the same software, Starlink, but differ in their choice of data reduction parameters. We find that the JCMT LR1 reduction is suitable for determining whether or not compact emission is present in a given region and the GBS LR1 reduction is tuned in a robust way to uncover more extended emission, which better serves more in-depth physical analyses of star-forming regions. Using the GBS LR1 method, we find that compact sources are recovered well, even at a peak brightness of only 3 times the noise, whereas the reconstruction of larger objects requires much care when drawing boundaries around the expected astronomical signal in the data reduction process. Incorrect boundaries can lead to false structure identification or it can cause structure to be missed. In the JCMT LR1 reduction, the extent of the true structure of objects larger than a point source is never fully recovered.

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.

An attitude determination system for balloon-borne experiments is presented. The system provides pointing information in azimuth and elevation for instruments flying on stratospheric balloons over Antarctica. In-flight attitude is given by the real-time combination of readings from star cameras, a magnetometer, sun sensors, GPS, gyroscopes, tilt sensors and an elevation encoder. Post-flight attitude reconstruction is determined from star camera solutions, interpolated by the gyroscopes using an extended Kalman Filter. The multi-sensor system was employed by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol), an experiment that measures polarized thermal emission from interstellar dust clouds. A similar system was designed for the upcoming flight of SPIDER, a Cosmic Microwave Background polarization experiment. The pointing requirements for these experiments are discussed, as well as the challenges in designing attitude reconstruction systems for high altitude balloon flights. In the 2010 and 2012 BLASTPol flights from McMurdo Station, Antarctica, the system demonstrated an accuracy of <5' rms in-flight, and <5" rms post-flight.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) is a suborbital mapping experiment, designed to study the role played by magnetic fields in the star formation process. BLASTPol observes polarized light using a total power instrument, photolithographic polarizing grids, and an achromatic half-wave plate to modulate the polarization signal. During its second flight from Antarctica in December 2012, BLASTPol made degree scale maps of linearly polarized dust emission from molecular clouds in three wavebands, centered at 250, 350, and 500 microns. The instrumental performance was an improvement over the 2010 BLASTPol flight, with decreased systematics resulting in a higher number of confirmed polarization vectors. The resultant dataset allows BLASTPol to trace magnetic fields in star-forming regions at scales ranging from cores to entire molecular cloud complexes.

We present the thermal model of the Balloon-borne Large-Aperture Submillimeter Telescope for Polarimetry (BLASTPol). This instrument was successfully flown in two circumpolar flights from McMurdo, Antarctica in 2010 and 2012. During these two flights, BLASTPol obtained unprecedented information about the magnetic field in molecular clouds through the measurement of the polarized thermal emission of interstellar dust grains. The thermal design of the experiment addresses the stability and control of the payload necessary for this kind of measurement. We describe the thermal modeling of the payload including the sun-shielding strategy. We present the in-flight thermal performance of the instrument and compare the predictions of the model with the temperatures registered during the flight. We describe the difficulties of modeling the thermal behavior of the balloon-borne platform and establish landmarks that can be used in the design of future balloon-borne instruments.

We present the second generation BLASTbus electronics. The primary purposes of this system are detector readout, attitude control, and cryogenic housekeeping, for balloon-borne telescopes. Readout of neutron transmutation doped germanium (NTD-Ge) bolometers requires low noise and parallel acquisition of hundreds of analog signals. Controlling a telescope's attitude requires the capability to interface to a wide variety of sensors and motors, and to use them together in a fast, closed loop. To achieve these different goals, the BLASTbus system employs a flexible motherboard-daughterboard architecture. The programmable motherboard features a digital signal processor (DSP) and field-programmable gate array (FPGA), as well as slots for three daughterboards. The daughterboards provide the interface to the outside world, with versions for analog to digital conversion, and optoisolated digital input/output. With the versatility afforded by this design, the BLASTbus also finds uses in cryogenic, thermometry, and power systems. For accurate timing control to tie everything together, the system operates in a fully synchronous manner. BLASTbus electronics have been successfully deployed to the South Pole, and flown on stratospheric balloons.

Turbulence and magnetic fields are expected to be important for regulating molecular cloud formation and evolution. However, their effects on subparsec to 100 parsec scales, leading to the formation of starless cores, is not well understood. We investigate the prestellar core structure morphologies obtained from analysis of the Herschel-SPIRE 350 $\mu$m maps of the Lupus I cloud. This distribution is first compared on a statistical basis to the large scale shape of the main filament. We find the distribution of the elongation position angle of the cores to be consistent with a random distribution, which means no specific orientation of the morphology of the cores is observed with respect to a large-scale filament shape model for Lupus I, or relative to a large-scale bent filament model. This distribution is also compared to the mean orientation of the large-scale magnetic fields probed at 350 $\mu$m with the Balloon-borne Large Aperture Telescope for Polarimetry (BLASTPol) during its 2010 campaign. Here again we do not find any correlation between the core morphology distribution and the average orientation of the magnetic fields on parsec scales. Our main conclusion is that the local filament dynamics - including secondary filaments that often run orthogonally to the primary filament - and possibly small-scale variations in the local magnetic field direction, could be the dominant factors for explaining the final orientation of each core.

We present observations of the Auriga-California Molecular Cloud (AMC) at 3.6, 4.5, 5.8, 8.0, 24, 70 and 160 micron observed with the IRAC and MIPS detectors as part of the Spitzer Gould Belt Legacy Survey. The total mapped areas are 2.5 sq-deg with IRAC and 10.47 sq-deg with MIPS. This giant molecular cloud is one of two in the nearby Gould Belt of star-forming regions, the other being the Orion A Molecular Cloud (OMC). We compare source counts, colors and magnitudes in our observed region to a subset of the SWIRE data that was processed through our pipeline. Using color-magnitude and color-color diagrams, we find evidence for a substantial population of 166 young stellar objects (YSOs) in the cloud, many of which were previously unknown. Most of this population is concentrated around the LkHalpha 101 cluster and the filament extending from it. We present a quantitative description of the degree of clustering and discuss the fraction of YSOs in the region with disks relative to an estimate of the diskless YSO population. Although the AMC is similar in mass, size and distance to the OMC, it is forming about 15 - 20 times fewer stars.

A cryogenic achromatic half-wave plate (HWP) for submillimetre astronomical polarimetry has been designed, manufactured, tested, and deployed in the Balloon-borne Large-Aperture Submillimeter Telescope for Polarimetry (BLASTPol). The design is based on the five-slab Pancharatnam recipe and it works in the wavelength range 200-600 micron, making it the broadest-band HWP built to date at (sub)millimetre wavelengths. The frequency behaviour of the HWP has been fully characterised at room and cryogenic temperatures with incoherent radiation from a polarising Fourier transform spectrometer. We develop a novel empirical model, complementary to the physical and analytical ones available in the literature, that allows us to recover the HWP Mueller matrix and phase shift as a function of frequency and extrapolated to 4K. We show that most of the HWP non-idealities can be modelled by quantifying one wavelength-dependent parameter, the position of the HWP equivalent axes, which is then readily implemented in a map-making algorithm. We derive this parameter for a range of spectral signatures of input astronomical sources relevant to BLASTPol, and provide a benchmark example of how our method can yield improved accuracy on measurements of the polarisation angle on the sky at submillimetre wavelengths.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was created by adding polarimetric capability to the BLAST experiment that was flown in 2003, 2005, and 2006. BLASTPol inherited BLAST's 1.8 m primary and its Herschel/SPIRE heritage focal plane that allows simultaneous observation at 250, 350, and 500 {\mu}m. We flew BLASTPol in 2010 and again in 2012. Both were long duration Antarctic flights. Here we present polarimetry of the nearby filamentary dark cloud Lupus I obtained during the 2010 flight. Despite limitations imposed by the effects of a damaged optical component, we were able to clearly detect submillimeter polarization on degree scales. We compare the resulting BLASTPol magnetic field map with a similar map made via optical polarimetry (The optical data were published in 1998 by J. Rizzo and collaborators.). The two maps partially overlap and are reasonably consistent with one another. We compare these magnetic field maps to the orientations of filaments in Lupus I, and we find that the dominant filament in the cloud is approximately perpendicular to the large-scale field, while secondary filaments appear to run parallel to the magnetic fields in their vicinities. This is similar to what is observed in Serpens South via near-IR polarimetry, and consistent with what is seen in MHD simulations by F. Nakamura and Z. Li.

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

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

We present first results from the Herschel Gould Belt survey for the B211/L1495 region in the Taurus molecular cloud. Thanks to their high sensitivity and dynamic range, the Herschel images reveal the structure of the dense, star-forming filament B211 with unprecedented detail, along with the presence of striations perpendicular to the filament and generally oriented along the magnetic field direction as traced by optical polarization vectors. Based on the column density and dust temperature maps derived from the Herschel data, we find that the radial density profile of the B211 filament approaches a power-law behavior {\rho} {\propto} r^(-2.0{\pm}0.4) at large radii and that the temperature profile exhibits a marked drop at small radii. The observed density and temperature profiles of the B211 filament are in good agreement with a theoretical model of a cylindrical filament undergoing gravitational contraction with a polytropic equation of state: P {\propto} {\rho}^{\gamma} and T {\propto} {\rho}^({\gamma}-1), with {\gamma}=0.97{\pm}0.01<1 (i.e. not strictly isothermal). The morphology of the column density map, where some of the perpendicular striations are apparently connected to the B211 filament, further suggests that the material may be accreting along the striations onto the main filament. The typical velocities expected for the infalling material in this picture are ~0.5-1 km/s, which are consistent with the existing kinematical constraints from previous CO observations.

We present observations of NGC1333 from SCUBA-2 on JCMT, observed as a JCMT Gould Belt Survey pilot project during the shared risk campaign when the first of four arrays was installed at each of 450 and 850 microns. Temperature maps are derived from 450 micron and 850 micron ratios under the assumption of constant dust opacity spectral index beta=1.8. Temperatures indicate that the dust in the northern (IRAS 6/8) region of NGC1333 is hot, 20-40 K, due to heating by the B star SVS3, other young stars in the IR/optically visible cluster, and embedded protostars. Other luminous protostars are also identified by temperature rises at the 17" resolution of the ratio maps (0.02 pc assuming a distance of 250 pc for Perseus). The extensive heating raises the possibility that the radiative feedback may lead to increased masses for the next generation of stars.

Observations of the dust emission using millimetre/submillimetre bolometer arrays can be contaminated by molecular line flux, such as flux from 12CO. As the brightest molecular line in the submillimetre, it is important to quantify the contribution of CO flux to the dust continuum bands. Conversion factors were used to convert molecular line integrated intensities to flux detected by bolometer arrays in mJy per beam. These factors were calculated for 12CO line integrated intensities to the SCUBA-2 850 {\mu}m and 450 {\mu}m bands. The conversion factors were then applied to HARP 12CO 3-2 maps of NGC 1333 in the Perseus complex and NGC 2071 and NGC 2024 in the Orion B molecular cloud complex to quantify the respective 12CO flux contribution to the 850 {\mu}m dust continuum emission. Sources with high molecular line contamination were analysed in further detail for molecular outflows and heating by nearby stars to determine the cause of the 12CO contribution. The majority of sources had a 12CO 3-2 flux contribution under 20 per cent. However, in regions of molecular outflows, the 12CO can dominate the source dust continuum (up to 79 per cent contamination) with 12CO fluxes reaching \sim 68 mJy per beam.

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

We propose an evolutionary path for prestellar cores on the radius-mass diagram, which is analogous to stellar evolutionary paths on the Hertzsprung-Russell Diagram. Using James Clerk Maxwell Telescope (JCMT) observations of L1688 in the Ophiuchus star-forming complex, we analyse the HCO+ (J=4\rightarrow3) spectral line profiles of prestellar cores. We find that of the 58 cores observed, 14 show signs of infall in the form of a blue-asymmetric double-peaked line profile. These 14 cores all lie beyond the Jeans mass line for the region on a radius-mass plot. Furthermore another 10 cores showing tentative signs of infall, in their spectral line profile shapes, appear on or just over the Jeans mass line. We therefore propose the manner in which a prestellar core evolves across this diagram. We hypothesise that a core is formed in the low-mass, low-radius region of the plot. It then accretes quasistatically, increasing in both mass and radius. When it crosses the limit of gravitational instability it begins to collapse, decreasing in radius, towards the region of the diagram where protostellar cores are seen.

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

We present the results of a comprehensive infrared, submillimetre, and millimetre continuum emission study of isolated low-mass star-forming cores in 32 Bok globules, with the aim to investigate the process of star formation in these regions. The submillimetre and millimetre dust continuum emission maps together with the spectral energy distributions are used to model and derive the physical properties of the star-forming cores, such as luminosities, sizes, masses, densities, etc. Comparisons with ground-based near-infrared and space-based mid and far-infrared images from Spitzer are used to reveal the stellar content of the Bok globules, association of embedded young stellar objects with the submm dust cores, and the evolutionary stages of the individual sources. Submm dust continuum emission was detected in 26 out of the 32 globule cores observed. For 18 globules with detected (sub)mm cores we derive evolutionary stages and physical parameters of the embedded sources. We identify nine starless cores, most of which are presumably prestellar, nine Class 0 protostars, and twelve Class I YSOs. Specific source properties like bolometric temperature, core size, and central densities are discussed as function of evolutionary stage. We find that at least two thirds (16 out of 24) of the star-forming globules studied here show evidence of forming multiple stars on scales between 1,000 and 50,000 AU. However, we also find that most of these small prototstar and star groups are comprised of sources with different evolutionary stages, suggesting a picture of slow and sequential star formation in isolated globules

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

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

We present optical and submillimetre polarimetry data of the Bok globule CB3 and optical polarimetry data of the Bok globule CB246. We use each set of polarimetry data to infer the B-field orientation in each of the clouds. The optical data can only be used in the low density, low extinction edge regions of clouds. The submillimetre data can only be used in the high column-density, central regions of the clouds. It has previously been found that near-infrared polarisation mapping of background stars does not accurately trace the magnetic field in dense cloud regions. This may be due to a lack of aligned grains in dense regions. We test this by comparing the field orientations measured by our two independent polarimetry methods. We find that the field orientation deduced from the optical data matches up well with the orientation estimated from the submillimetre data. We therefore claim that both methods are accurately tracing the same magnetic field in CB3. Hence, in this case, there must be significant numbers of aligned dust grains in the high density region, and they do indeed trace the magnetic field in the submillimetre. We find an offset of 40$\pm$14 degrees between the magnetic field orientation and the short axis of the globule. This is consistent with the mean value of 31$\pm$3 degrees found in our previous work on prestellar cores, even though CB3 is a protostellar core. Taken together, the six prestellar cores that we have now studied in this way show a mean offset between magnetic field orientation and core short axis of $\sim30\pm$3 degrees, in apparent contradiction with some models of magnetically dominated star formation.

We present observations of L1155 and L1148 in the Cepheus molecular cloud, taken using the FIS instrument on the Akari satellite. We compare these data to submillimetre data taken using the SCUBA camera on the JCMT, and far-infrared data taken with the ISOPHOT camera on board the ISO satellite. All of the data show a relation between the position of the peak of emission and the wavelength for the core of L1155. We interpret this as a temperature gradient. We fit modified blackbody curves to the spectral energy distributions at two positions in the core and see that the central core in L1155 (L1155C) is approximately 2 degrees warmer at one edge than it is in the centre. We consider a number of possible heating sources and conclude that the A6V star BD+67 1263 is the most likely candidate. This star is at a distance of 0.7 pc from the front of L1155C in the plane of the sky. We carry out radiative transfer modelling of the L1155C core including the effects from the nearby star. We find that we can generate a good fit to the observed data at all wavelengths, and demonstrate that the different morphologies of the core at different wavelengths can be explained by the observed 2 degree temperature gradient. The L1148 core exhibits a similar morphology to that of L1155C, and the data are also consistent with a temperature gradient across the core. In this case, the most likely heating source is the star BD197053. Our findings illustrate very clearly that the apparent observed morphology of a pre-stellar core can be highly dependent on the wavelength of the observation, and that temperature gradients must be taken into account before converting images into column density distributions. This is important to note when interpreting Akari and Spitzer data and will also be significant for Herschel data.

We re-analyse all of the archive observations of the Ophiuchus dark cloud L1688 that were carried out with the submillimetre common-user bolometer array (SCUBA) at the James Clerk Maxwell Telescope (JCMT). For the first time we put together all of the data that were taken of this cloud at different times to make a deeper map at 850 microns than has ever previously been published. Using this new, deeper map we extract the pre-stellar cores from the data. We use updated values for the distance to the cloud complex, and also for the internal temperatures of the pre-stellar cores to generate an updated core mass function (CMF). This updated CMF is consistent with previous results in so far as they went, but our deeper map gives an improved completeness limit of 0.1 Mo (0.16 Jy), which enables us to show that a turnover exists in the low-mass regime of the CMF. The L1688 CMF shows the same form as the stellar IMF and can be mapped onto the stellar IMF, showing that the IMF is determined at the prestellar core stage. We compare L1688 with the Orion star-forming region and find that the turnover in the L1688 CMF occurs at a mass roughly a factor of two lower than the CMF turnover in Orion. This suggests that the position of the CMF turnover may be a function of environment.

We present continuum data from the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), and the Mid-Infrared Photometer for Spitzer (MIPS) on the Spitzer Space Telescope, at submillimetre and infrared wavelengths respectively. We study the Taurus molecular cloud 1 (TMC1), and in particular the region of the Taurus Molecular Ring (TMR). In the continuum data we see no real evidence for a ring, but rather we see one side of it only, appearing as a filament. We name the filament `the bull's tail'. The filament is seen in emission at 850, 450 and 160um, and in absorption at 70um. We compare the data with archive data from the Infra-Red Astronomical Satellite (IRAS) at 12, 25, 60, 100um, in which the filament is also seen in absorption. We find that the emission from the filament consists of two components: a narrow, cold (~8K), central core; and a broader, slightly warmer (~12K), shoulder of emission. We use a radiative transfer code to model the filament's appearance, either in emission or absorption, simultaneously at each of the different wavelengths. Our best fit model uses a Plummer-like density profile and a homogeneous interstellar dust grain population. Unlike previous work on a similar, but different filament in Taurus, we require no grain coagulation to explain our data.

Stars form from dense molecular cores, and the mass function of these cores (the CMF) is often found to be similar to the form of the stellar initial mass function (IMF). This suggests that the form of the IMF is the result of the form of the CMF. However, most stars are thought to form in binary and multiple systems, therefore the relationship between the IMF and the CMF cannot be trivial. We test two star formation scenarios - one in which all stars form as binary or triple systems, and one in which low-mass stars form in a predominantly single mode. We show that from a log-normal CMF, similar to those observed, and expected on theoretical grounds, the model in which all stars form as multiples gives a better fit to the IMF.

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.

The sub-millimetre wavelength regime is perhaps the most poorly explored over large areas of the sky, despite the considerable effort that has been expended in making deep maps over small regions. As a consequence the properties of the sub-millimetre sky as a whole, and of rare bright objects in particular, remains largely unknown. Here we describe a forthcoming survey (the SCUBA-2 ``All-Sky'' Survey, or SASSy) designed to address this issue by making a large-area map of approximately one-fifth of the sky visible from the JCMT (4800 square degrees) down to a 1 sigma noise level of 30 mJy/beam. This map forms the pilot for a much larger survey, which will potentially map the remaining sky visible from the JCMT, with the region also visible to ALMA as a priority. SASSy has been awarded 500 hours for the 4800 square degree pilot phase and will commence after the commissioning of SCUBA-2, expected in early 2008.

We have re-analysed all of the SCUBA archive data of the Orion star-forming regions. We have put together all of the data taken at different times by different groups. Consequently we have constructed the deepest submillimetre maps of these regions ever made. There are four regions that have been mapped: Orion A North & South, and Orion B North & South. We find that two of the regions, Orion A North and Orion B North, have deeper sensitivity and completeness limits, and contain a larger number of sources, so we concentrate on these two. We compare the data with archive data from the Spitzer Space Telescope to determine whether or not a core detected in the submillimetre is pre-stellar in nature. We extract all of the pre-stellar cores from the data and make a histogram of the core masses. This can be compared to the stellar initial mass function (IMF). We find the high-mass core mass function follows a roughly Salpeter-like slope, just like the IMF, as seen in previous work. Our deeper maps allow us to see that the core mass function (CMF) turns over at ~ 1.3 Mo, about a factor of 4 higher than our completeness limit. This turnover has never previously been observed, and is only visible here due to our much deeper maps. It mimics the turnover seen in the stellar IMF at ~ 0.1 Mo. The low-mass side of the CMF is a power-law with an exponent of 0.35 +/- 0.2, which is consistent with the low-mass slope of the young cluster IMF of 0.3 +/- 0.1. This shows that the CMF continues to mimic the shape of the IMF all the way down to the lower completeness limit of these data at ~ 0.3 Mo.

We present observations taken with SCUBA on the JCMT of the Horsehead Nebula in Orion (B33), at wavelengths of 450 and 850 \mum. We see bright emission from that part of the cloud associated with the photon-dominated region (PDR) at the `top' of the horse's head, which we label B33-SMM1. We characterise the physical parameters of the extended dust responsible for this emission, and find that B33-SMM1 contains a more dense core than was previously suspected. We compare the SCUBA data with data from the Infrared Space Observatory (ISO) and find that the emission at 6.75-\mum is offset towards the west, indicating that the mid-infrared emission is tracing the PDR while the submillimetre emission comes from the molecular cloud core behind the PDR. We calculate the virial balance of this core and find that it is not gravitationally bound but is being confined by the external pressure from the HII region IC434, and that it will either be destroyed by the ionising radiation, or else may undergo triggered star formation. Furthermore we find evidence for a lozenge-shaped clump in the `throat' of the horse, which is not seen in emission at shorter wavelengths. We label this source B33-SMM2 and find that it is brighter at submillimetre wavelengths than B33-SMM1. SMM2 is seen in absorption in the 6.75-\mum ISO data, from which we obtain an independent estimate of the column density in excellent agreement with that calculated from the submillimetre emission. We calculate the stability of this core against collapse and find that it is in approximate gravitational virial equilibrium. This is consistent with it being a pre-existing core in B33, possibly pre-stellar in nature, but that it may also eventually undergo collapse under the effects of the HII region.

We present submillimetre data for the L1689 cloud in the rho-Ophiuchi molecular cloud complex. We detect a number of starless and prestellar cores and protostellar envelopes. We also detect a number of filaments for the first time in the submillimetre continuum that are parallel both to each other, and to filaments observed in the neighbouring L1688 cloud. These filaments are also seen in the 13CO observations of L1689. The filaments contain all of the star-formation activity in the cloud. L1689 lies next to the well studied L1688 cloud that contains the rho Oph-A core. L1688 has a much more active star-formation history than L1689 despite their apparent similarity in 13CO data. Hence we label L1689 as the dog that didn't bark. We endeavour to explain this apparent anomaly by comparing the total mass of each cloud that is currently in the form of dense material such as prestellar cores. We note firstly that L1688 is more massive than L1689, but we also find that when normalised to the total mass of each cloud, the L1689 cloud has a much lower percentage of mass in dense cores than L1688. We attribute this to the hypothesis of Loren (1989) that the star formation in the rho-Ophiuchi complex is being affected and probably dominated by the external influence of the nearby Upper Scorpius OB association and predominantly by sigma-Sco. L1689 is further from sigma-Sco and is therefore less active. The influence of sigma-Sco appears nonetheless to have created the filaments that we observe in L1689. Accepted by MNRAS.

We present 450 and 850 um maps of R Coronae Australis. We compare the maps to previous surveys of the region, and shed new light on the previously unknown nature of the protostellar sources at the centre of the cloud. We clarify the nature of two millimetre sources previously discovered in lower resolution data. We identify one new Class 0 protostar that we label SMM 1B, and we measure the envelope masses of a number of more evolved protostars. We identify two new prestellar cores that we call SMM 1A and SMM 6.

We estimate the levels of turbulence in the envelopes of Class 0 and I protostars using a model based on measurements of the peak separation of double-peaked asymmetric line profiles. We use observations of 20 protostars of both Class 0 and I taken in the HCO+ (J=3-2) line that show the classic double-peaked profile. We find that some Class 0 sources show high levels of turbulence whilst others demonstrate much lower levels. In Class I protostars we find predominantly low levels of turbulence. The observations are consistent with a scenario in which Class 0 protostars form in a variety of environments and subsequently evolve into Class I protostars. The data do not appear to be consistent with a recently proposed scenario in which Class 0 protostars can only form in extreme environments.

We have mapped linearly polarized dust emission from L183 with the JCMT SCUBA polarimeter and have analyzed these and our previously published data for the prestellar cores L183, L1544, and L43 in order to estimate magnetic field strengths in the plane of the sky, $B_{pos}$. The analysis used the Chandrasekhar-Fermi technique, which relates the dispersion in polarization position angles to $B_{pos}$. We have used these estimates of the field strengths (neglecting the unmeasured line-of-sight component) to find the mass-to-magnetic flux ratios $\lambda$ (in units of the critical ratio for magnetic support). Results are $B_{pos} \approx 80$ $\mu$G and $\lambda \approx 2.6$ for L183, $B_{pos} \approx 140$ $\mu$G and $\lambda \approx 2.3$ for L1544, and $B_{pos} \approx 160$ $\mu$G and $\lambda \approx 1.9$ for L43. Hence, without correction for geometrical biases, for all three cores the mass-to-flux ratios are supercritical by a factor of $\sim 2$, and magnetic support cannot prevent collapse. However, a statistical mean correction for geometrical bias may be up to a factor of three; this correction would reduce the individual $\lambda$'s to $\lambda_{cor} \approx 0.9$, 0.8, and 0.6, respectively; these values are approximately critical or slightly subcritical. These data are consistent with models of star formation driven by ambipolar diffusion in a weakly turbulent medium, but cannot rule out models of star formation driven by turbulence.

A brief summary is presented of our current knowledge of the structure of cold molecular cloud cores that do not contain protostars, sometimes known as starless cores. The most centrally condensed starless cores are known as pre-stellar cores. These cores probably represent observationally the initial conditions for protostellar collapse that must be input into all models of star formation. The current debate over the nature of core density profiles is summarised. A cautionary note is sounded over the use of such profiles to ascertain the equilibrium status of cores. The magnetic field structure of pre-stellar cores is also discussed.

C17O (J=2-1) observations have been carried out towards the pre-stellar core L1689B. By comparing the relative strengths of the hyperfine components of this line, the emission is shown to be optically thin. This allows accurate CO column densities to be determined and, for reference, this calculation is described in detail. The hydrogen column densities that these measurements imply are substantially smaller than those calculated from SCUBA dust emission data. Furthermore, the C17O column densities are approximately constant across L1689B whereas the SCUBA column densities are peaked towards the centre. The most likely explanation is that CO is depleted from the central regions of L1689B. Simple models of pre-stellar cores with an inner depleted region are compared with the results. This enables the magnitude of the CO depletion to be quantified and also allows the spatial extent of the freeze-out to be firmly established. We estimate that within about 5000 AU of the centre of L1689B, over 90% of the CO has frozen onto grains. This level of depletion can only be achieved after a duration that is at least comparable to the free-fall timescale.