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.

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.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We present the properties of magnetic fields (B-fields) in two clumps (clump 1 and clump 2), located at the waist of the bipolar H II region Sh2-201, based on JCMT SCUBA-2/POL-2 observations of 850 $\mu$m polarized dust emission. We find that B-fields in the direction of the clumps are bent and compressed, showing bow-like morphologies, which we attribute to the feedback effect of the H II region on the surface of the clumps. Using the modified Davis-Chandrasekhar-Fermi method we estimate B-fields strengths of 266 $\mu$G and 65 $\mu$G for clump 1 and clump 2, respectively. From virial analyses and critical mass ratio estimates, we argue that clump 1 is gravitationally bound and could be undergoing collapse, whereas clump 2 is unbound and stable. We hypothesize that the interplay between thermal pressure imparted by the H II region, B-field morphologies, and the various internal pressures of the clumps (such as magnetic, turbulent, and gas thermal pressure), has the following consequences: (a) formation of clumps at the waist of the H II region; (b) progressive compression and enhancement of the B-fields in the clumps; (c) stronger B-fields will shield the clumps from erosion by the H II region and cause pressure equilibrium between the clumps and the H II region, thereby allowing expanding I-fronts to blow away from the filament ridge, forming bipolar H II regions; and (d) stronger B-fields and turbulence will be able to stabilize the clumps. A study of a larger sample of bipolar H II regions would help to determine whether our hypotheses are widely applicable.

We investigate the spatial properties of sources from the GAIA catalogue previously identified as being members of the Taurus star forming region and which appear in the Spitzer catalogue. We study an area of sky of 10x15 deg, centred on Right Ascension (2000.0) = 68.5 deg and Declination (2000.0)= 27.0 deg, this being an area surrounding the Taurus molecular cloud. We use data obtained from the GAIA DR2 release. By using an inversion of GAIA parallax measurements to obtain distance values and by defining limits to the proper motions of the Taurus moving group, we are able to show that there are substantial differences in depth within the Taurus complex. Our results suggest that the Taurus cloud has significant depth and that there are two main associations centred at ~130.6+/-0.7 pc and ~160.2+/-0.9 pc at 1{\sigma}. These two associations also have different proper motions, of 24.5+/-2.8 and 20.1+/-2.4 mas/yr respectively. We here label them the 'Two Horns' of Taurus.

Herschel observations of nearby clouds in the Gould Belt support a paradigm for low-mass star formation, starting with the generation of molecular filaments, followed by filament fragmentation, and the concentration of mass into self-gravitating prestellar cores. With the unique far-infrared and submillimeter continuum imaging capabilities of the Herschel Space observatory, the closeby (d = 139 pc) Ophiuchus cloud was mapped at five wavelengths from 70 ${\mu}m$ to 500 ${\mu}m$ with the aim of providing a complete census of dense cores in this region, including unbound starless cores, bound prestellar cores, and protostellar cores. Taking advantage of the high dynamic range and multi-wavelength nature of the Herschel data, we used the multi-scale decomposition algorithms getsources and getfilaments to identify a complete sample of dense cores and filaments in the cloud and study their properties. The densest clouds of the Ophiuchus complex, L1688 and L1689, which thus far are only indirectly described as filamentary regions owing to the spatial distribution of their young stellar objects (YSOs), are confirmed to be dominated by filamentary structures. The tight correlation observed between prestellar cores and filamentary structures in L1688 and L1689 supports the view that solar-type star formation occurs primarily in dense filaments. While the sub clouds of the complex show disparities, L1689 being less efficient than L1688 at forming stars when considering their total mass budgets, both sub clouds share almost the same prestellar core formation efficiency in dense molecular gas. We also find evidence in the Herschel data for a remarkable concentric geometrical configuration in L1688 which is dominated by up to three arc-like compression fronts and presumably created by shockwave events emanating from the Sco OB2 association, including the neighboring massive (O9V) star ${\sigma}$ Sco.

We present a detailed study of the Orion B clouds (d~400 pc), imaged with the PACS/SPIRE cameras at 70-500 $\mu$m by the Herschel Gould Belt survey (HGBS). We release new high-res. maps of column density and dust temperature. In the filamentary sub-regions NGC2023/2024, NGC2068/2071, and L1622, 1768 starless dense cores were identified, ~28-45% of which are self-gravitating prestellar cores. A total of 76 protostellar dense cores were also found. The typical lifetime of the prestellar cores was found to be $t_{\rm pre}=1.7_{-0.6}^{+0.8}$ Myr. The prestellar core mass function (CMF) peaks at ~0.5 $M_\odot$ and is consistent with a power law with log slope -1.27$\pm$0.24 at the high-mass end, compared to the Salpeter slope of -1.35. In this region, we confirm the existence of a transition in prestellar core formation efficiency (CFE) around a fiducial value A_V_bg~7 mag in background visual extinction, similar to the trend observed with Herschel in other clouds. This is not a sharp threshold, but a smooth transition between a regime with very low prestellar CFE at A_V_bg<5 and a regime with higher, roughly constant CFE at A_V_bg$\gtrsim$10. The total mass in the form of prestellar cores represents only ~20% of the dense molecular cloud gas at A_V_bg$\gtrsim$7 mag. About 60-80% of the prestellar cores are closely associated with filaments, and this fraction increases up to >90% when a more complete sample of filamentary structures is considered. Interestingly, the median separation between nearest core neighbors corresponds to the typical inner filament width of ~0.1 pc commonly observed in nearby molecular clouds. Analysis of the CMF observed as a function of background cloud column density shows that the most massive prestellar cores are spatially segregated in the highest column density areas, and suggests that both higher- and lower-mass prestellar cores may form in denser filaments.

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

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

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

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

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

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.

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 present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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

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

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

Linear polarization maps of the Carina Nebula were obtained at 250, 350, and 500 $\mu$m during the 2012 flight of the BLASTPol balloon-borne telescope. These measurements are combined with Planck 850 $\mu$m data in order to produce a submillimeter spectrum of the polarization fraction of the dust emission, averaged over the cloud. This spectrum is flat to within $\pm$15% (relative to the 350 $\mu$m polarization fraction). In particular, there is no evidence for a pronounced minimum of the spectrum near 350 $\mu$m, as suggested by previous ground-based measurements of other molecular clouds. This result of a flat polarization spectrum in Carina is consistent with recently-published BLASTPol measurements of the Vela C molecular cloud, and also agrees with a published model for an externally-illuminated, dense molecular cloud by Bethell and collaborators. The shape of the spectrum in Carina does not show any dependence on the radiative environment of the dust, as quantified by the Planck-derived dust temperature or dust optical depth at 353 GHz.

The Herschel Gould Belt Survey key project mapped the bulk of nearby star-forming molecular clouds in five far-infrared bands with the aim of compiling complete census of prestellar cores and young, embedded protostars. In this paper, we present the catalogue of the dense cores and YSOs/protostars extracted from the Herschel maps of the Lupus I, III, and IV molecular clouds. The physical properties of the detected objects were derived by fitting their spectral energy distributions. A total of 532 dense cores, out of which 103 are presumably prestellar in nature, and 38 YSOs/protostars have been detected in the three clouds. Almost all the prestellar cores are associated with filaments against only about one third of the unbound cores and YSOs/protostars. Prestellar core candidates are found even in filaments that are on average thermally sub-critical and over a background column density lower than that measured in other star forming regions so far. The core mass function of the prestellar cores peaks between 0.2 and 0.3 solar masses and it is compatible with the log-normal shape found in other regions. Herschel data reveal several, previously undetected, protostars and new candidates of Class 0 and Class II with transitional disks. We estimate the evolutionary status of the YSOs/protostars using two independent indicators: the $\alpha$ index and the fitting of the spectral energy distribution from near- to far-infrared wavelengths. For 70% of the objects, the evolutionary stages derived with the two methods are in agreement. Lupus is confirmed to be a very low-mass star-forming region, both in terms of the prestellar condensations and of the diffuse medium. Noticeably, in the Lupus clouds we have found star formation activity associated with interstellar medium at low column density, usually quiescent in other (more massive) star forming regions.

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

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

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

We present the first high-resolution, submillimeter-wavelength polarimetric observations of -- and thus direct observations of the magnetic field morphology within -- the dense gas of the Pillars of Creation in M16. These 850$\,\mu$m observations, taken as part of the BISTRO (B-Fields in Star-forming Region Observations) Survey using the POL-2 polarimeter on the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT), show that the magnetic field runs along the length of the pillars, perpendicular to, and decoupled from, the field in the surrounding photoionized cloud. Using the Chandrasekhar-Fermi method we estimate a plane-of-sky magnetic field strength of $170-320\,\mu$G in the Pillars, consistent with their having been formed through compression of gas with initially weak magnetization. The observed magnetic field strength and morphology suggests that the magnetic field may be slowing the pillars' evolution into cometary globules. We thus hypothesize that the evolution and lifetime of the Pillars may be strongly influenced by the strength of the coupling of their magnetic field to that of their parent photoionized cloud -- i.e. that the Pillars' longevity results from magnetic support

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

Planck Galactic Cold Clumps (PGCCs) possibly represent the early stages of star formation. To understand better the properties of PGCCs, we studied 16 PGCCs in the L1495 cloud with molecular lines and continuum data from Herschel, JCMT/SCUBA-2 and the PMO 13.7 m telescope. Thirty dense cores were identified in 16 PGCCs from 2-D Gaussian fitting. The dense cores have dust temperatures of $T_{\rm d}$ = 11-14 K, and H$_{2}$ column densities of $N_{\rm H_{2}}$ = 0.36-2.5$\times10^{22}$ cm$^{-2}$. We found that not all PGCCs contain prestellar objects. In general, the dense cores in PGCCs are usually at their earliest evolutionary stages. All the dense cores have non-thermal velocity dispersions larger than the thermal velocity dispersions from molecular line data, suggesting that the dense cores may be turbulence-dominated. We have calculated the virial parameter $\alpha$ and found that 14 of the dense cores have $\alpha$ $<$ 2, while 16 of the dense cores have $\alpha$ $>$ 2. This suggests that some of the dense cores are not bound in the absence of external pressure and magnetic fields. The column density profiles of dense cores were fitted. The sizes of the flat regions and core radii decrease with the evolution of dense cores. CO depletion was found to occur in all the dense cores, but is more significant in prestellar core candidates than in protostellar or starless cores. The protostellar cores inside the PGCCs are still at a very early evolutionary stage, sharing similar physical and chemical properties with the prestellar core candidates.

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

We present a catalogue of prestellar and starless cores within the Corona Australis molecular cloud using photometric data from the Herschel Space Observatory. At a distance of d~130 pc, Corona Australis is one of the closest star-forming regions. Herschel has taken multi-wavelength data of Corona Australis with both the SPIRE and PACS photometric cameras in a parallel mode with wavelengths in the range 70 {\mu}m to 500 {\mu}m. A complete sample of starless and prestellar cores and embedded protostars is identified. Other results from the Herschel Gould Belt Survey have shown spatial correlation between the distribution of dense cores and the filamentary structure within the molecular clouds. We go further and show correlations between the properties of these cores and their spatial distribution within the clouds, with a particular focus on the mass distribution of the dense cores with respect to their filamentary proximity. We find that only lower-mass starless cores form away from filaments, while all of the higher-mass prestellar cores form in close proximity to, or directly on the filamentary structure. This result supports the paradigm that prestellar cores mostly form on filaments. We analyse the mass distribution across the molecular cloud, finding evidence that the region around the Coronet appears to be at a more dynamically advanced evolutionary stage to the rest of the clumps within the cloud.

Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths.

We present far-infrared observations of Monoceros R2 (a giant molecular cloud at approximately 830 pc distance, containing several sites of active star formation), as observed at 70 {\mu}m, 160 {\mu}m, 250 {\mu}m, 350 {\mu}m, and 500 {\mu}m by the Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver (SPIRE) instruments on the Herschel Space Observatory as part of the Herschel imaging survey of OB young stellar objects (HOBYS) Key programme. The Herschel data are complemented by SCUBA-2 data in the submillimetre range, and WISE and Spitzer data in the mid-infrared. In addition, C18O data from the IRAM 30-m Telescope are presented, and used for kinematic information. Sources were extracted from the maps with getsources, and from the fluxes measured, spectral energy distributions were constructed, allowing measurements of source mass and dust temperature. Of 177 Herschel sources robustly detected in the region (a detection with high signal-to-noise and low axis ratio at multiple wavelengths), including protostars and starless cores, 29 are found in a filamentary hub at the centre of the region (a little over 1% of the observed area). These objects are on average smaller, more massive, and more luminous than those in the surrounding regions (which together suggest that they are at a later stage of evolution), a result that cannot be explained entirely by selection effects. These results suggest a picture in which the hub may have begun star formation at a point significantly earlier than the outer regions, possibly forming as a result of feedback from earlier star formation. Furthermore, the hub may be sustaining its star formation by accreting material from the surrounding filaments.

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

Polarized emission from aligned dust is a crucial tool for studies of magnetism in the ISM and a troublesome contaminant for studies of CMB polarization. In each case, an understanding of the significance of the polarization signal requires well-calibrated physical models of dust grains. Despite decades of progress in theory and observation, polarized dust models remain largely underconstrained. During its 2012 flight, the balloon-borne telescope BLASTPol obtained simultaneous broad-band polarimetric maps of a translucent molecular cloud at 250, 350, and 500 microns. Combining these data with polarimetry from the Planck 850 micron band, we have produced a submillimeter polarization spectrum for a cloud of this type for the first time. We find the polarization degree to be largely constant across the four bands. This result introduces a new observable with the potential to place strong empirical constraints on ISM dust polarization models in a previously inaccessible density regime. Comparing with models by Draine and Fraisse (2009), our result disfavors two of their models for which all polarization arises due only to aligned silicate grains. By creating simple models for polarized emission in a translucent cloud, we verify that extinction within the cloud should have only a small effect on the polarization spectrum shape compared to the diffuse ISM. Thus we expect the measured polarization spectrum to be a valid check on diffuse ISM dust models. The general flatness of the observed polarization spectrum suggests a challenge to models where temperature and alignment degree are strongly correlated across major dust components.

Debris discs are evidence of the ongoing destructive collisions between planetesimals, and their presence around stars also suggests that planets exist in these systems. In this paper, we present submillimetre images of the thermal emission from debris discs that formed the SCUBA-2 Observations of Nearby Stars (SONS) survey, one of seven legacy surveys undertaken on the James Clerk Maxwell telescope between 2012 and 2015. The overall results of the survey are presented in the form of 850 microns (and 450 microns, where possible) images and fluxes for the observed fields. Excess thermal emission, over that expected from the stellar photosphere, is detected around 49 stars out of the 100 observed fields. The discs are characterised in terms of their flux density, size (radial distribution of the dust) and derived dust properties from their spectral energy distributions. The results show discs over a range of sizes, typically 1-10 times the diameter of the Edgeworth-Kuiper Belt in our Solar System. The mass of a disc, for particles up to a few millimetres in size, is uniquely obtainable with submillimetre observations and this quantity is presented as a function of the host stars' age, showing a tentative decline in mass with age. Having doubled the number of imaged discs at submillimetre wavelengths from ground-based, single dish telescope observations, one of the key legacy products from the SONS survey is to provide a comprehensive target list to observe at high angular resolution using submillimetre/millimetre interferometers (e.g., ALMA, SMA).

Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

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 a large-scale combination of near-infrared (near-IR) interstellar polarization data from background starlight with polarized emission data at submillimeter (sub-mm) wavelengths for the Vela C molecular cloud. The near-IR data consist of more than 6700 detections probing a range of visual extinctions between $2$ and $20\,$mag in and around the cloud. The sub-mm data was collected in Antartica by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). This is the first direct combination of near-IR and sub-mm polarization data for a molecular cloud aimed at measuring the "polarization efficiency ratio" ($R_{\mathrm{eff}}$), a quantity that is expected to depend only on grain intrinsic physical properties. It is defined as $p_{500}/(p_{I}/\tau_{V})$, where $p_{500}$ and $p_{I}$ are polarization fractions at $500\,\mu$m and $I$-band, respectively, and $\tau_{V}$ is the optical depth. To ensure that the same column density of material is producing both polarization from emission and from extinction, we conducted a careful selection of near-background stars using 2MASS, $Herschel$ and $Planck$ data. This selection excludes objects contaminated by the Galactic diffuse background material as well as objects located in the foreground. Accounting for statistical and systematic uncertainties, we estimate an average $R_{\mathrm{eff}}$ value of $2.4\pm0.8$, which can be used to test the predictions of dust grain models designed for molecular clouds when such predictions become available. $R_{\mathrm{eff}}$ appears to be relatively flat as a function of the cloud depth for the range of visual extinctions probed.

Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two main components: the central (Gaussian) region (core component), and the power-law like region dominating the filament column density profile at larger radii (wing component). The properties and behaviour of these components relative to the total linear mass density of the filament and its environmental column density were compared with theoretical models describing the evolution of filaments under gravity-dominated conditions. The feasibility of a transition to supercritical state by accretion is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high Mline-core) filaments in dense environments (av\sim3mag) can become supercritical in timescales of t\sim1Myr by accreting mass at constant or decreasing width. The trend of increasing Mline-tot (Mline-core and Mline-wing), and ridge Av with background also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the association between dense environments and high Mline-core values, and argues against filaments remaining in constant single-pressure equilibrium states. The simultaneous growth of filament and background in locations with efficient mass assembly, predicted in numerical models of collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations.

We statistically evaluate the relative orientation between gas column density structures, inferred from Herschel submillimetre observations, and the magnetic field projected on the plane of sky, inferred from polarized thermal emission of Galactic dust observed by BLASTPol at 250, 350, and 500 micron, towards the Vela C molecular complex. First, we find very good agreement between the polarization orientations in the three wavelength-bands, suggesting that, at the considered common angular resolution of 3.0 arcminutes that corresponds to a physical scale of approximately 0.61 pc, the inferred magnetic field orientation is not significantly affected by temperature or dust grain alignment effects. Second, we find that the relative orientation between gas column density structures and the magnetic field changes progressively with increasing gas column density, from mostly parallel or having no preferred orientation at low column densities to mostly perpendicular at the highest column densities. This observation is in agreement with previous studies by the Planck collaboration towards more nearby molecular clouds. Finally, we find a correspondence between the trends in relative orientation and the shape of the column density probability distribution functions. In the sub-regions of Vela C dominated by one clear filamentary structure, or "ridges", we find a sharp transition from preferentially parallel or having no preferred relative orientation at low column densities to preferentially perpendicular at highest column densities. In the sub-regions of Vela C dominated by several filamentary structures with multiple orientations, or "nests", such a transition is also present, but it is clearly less sharp than in the ridge-like sub-regions. Both of these results suggest that the magnetic field is dynamically important for the formation of density structures in this region.

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.

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

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

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.

Polarization maps of the Vela C molecular cloud were obtained at 250, 350, and 500um during the 2012 flight of the balloon-borne telescope BLASTPol. These measurements are used in conjunction with 850um data from Planck to study the submillimeter spectrum of the polarization fraction for this cloud. The spectrum is relatively flat and does not exhibit a pronounced minimum at \lambda ~350um as suggested by previous measurements of other molecular clouds. The shape of the spectrum does not depend strongly on the radiative environment of the dust, as quantified by the column density or the dust temperature obtained from Herschel data. The polarization ratios observed in Vela C are consistent with a model of a porous clumpy molecular cloud being uniformly heated by the interstellar radiation field.

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 results for Vela C obtained during the 2012 flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). We mapped polarized intensity across almost the entire extent of this giant molecular cloud, in bands centered at 250, 350, and 500 {\mu}m. In this initial paper, we show our 500 {\mu}m data smoothed to a resolution of 2.5 arcminutes (approximately 0.5 pc). We show that the mean level of the fractional polarization p and most of its spatial variations can be accounted for using an empirical three-parameter power-law fit, p = p_0 N^(-0.4) S^(-0.6), where N is the hydrogen column density and S is the polarization-angle dispersion on 0.5 pc scales. The decrease of p with increasing S is expected because changes in the magnetic field direction within the cloud volume sampled by each measurement will lead to cancellation of polarization signals. The decrease of p with increasing N might be caused by the same effect, if magnetic field disorder increases for high column density sightlines. Alternatively, the intrinsic polarization efficiency of the dust grain population might be lower for material along higher density sightlines. We find no significant correlation between N and S. Comparison of observed submillimeter polarization maps with synthetic polarization maps derived from numerical simulations provides a promising method for testing star formation theories. Realistic simulations should allow for the possibility of variable intrinsic polarization efficiency. The measured levels of correlation among p, N, and S provide points of comparison between observations and simulations.

(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68{\deg} > l > -70{\deg} in a |b|<1{\deg} latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.

The association of filaments with protostellar objects has made these structures a priority target in star formation studies. The datasets of the Herschel Galactic Cold Cores Key Programme allow for a statistical study of filaments with a wide range of intrinsic and environmental characteristics. Characterisation of this sample can be used to identify key physical parameters and quantify the role of environment in the formation of supercritical filaments. Filaments were extracted from fields at D<500pc with the getfilaments algorithm and characterised according to their column density profiles and intrinsic properties. Each profile was fitted with a beam-convolved Plummer-like function and quantified based on the relative contributions from the filament 'core', represented by a Gaussian, and 'wing' component, dominated by the power-law of the Plummer-like function. These parameters were examined for populations associated with different background levels. We find that filaments increase their core (Mcore) and wing (Mwing) contributions while increasing their total linear mass density (Mtot). Both components appear to be linked to the local environment, with filaments in higher backgrounds having systematically more massive Mcore and Mwing. This dependence on the environment supports an accretion-based model for filament evolution in the local neighbourhood (D<500pc). Structures located in the highest backgrounds develop the highest central Av, Mcore, and Mwing as Mtot increases with time, favoured by the local availability of material and the enhanced gravitational potential. Our results indicate that filaments acquiring a significantly massive central region with Mcore>Mcrit/2 may become supercritical and form stars. This translates into a need for filaments to become at least moderately self-gravitating in order to undergo localised star formation or become star-forming filaments.

We observed high S/N, high velocity resolution NH$_3$(1,1) and (2,2) emission on an extended map in TMC-1, a filamentary cloud in a nearby quiescent star forming area. By fitting multiple hyperfine-split line profiles to the NH$_3$(1,1) spectra we derived the velocity distribution of the line components and calculated gas parameters on several positions. Herschel SPIRE continuum observations were reduced and used to calculate the physical parameters of the Planck Galactic Cold Clumps in the region. The Herschel-based column density map of TMC-1 shows a main ridge with two local maxima and a separated peak to the south-west. H$_2$-column densities and dust temperatures are in the range of 0.5-3.3 $\times$ 10$^{22}$ cm$^{-2}$ and 10.5-12 K, respectively. NH$_3$-column densities are 2.8-14.2 $\times$ 10$^{14}$ cm$^{-2}$ and and H$_2$-volume densities are 0.4-2.8 $\times$ 10$^4$ cm$^{-3}$. Kinetic temperatures are typically very low with a minimum of 9 K, and a maximum of 13.7 K was found at the Class I protostar IRAS 04381+2540. The kinetic temperatures vary similarly as the dust temperatures in spite of the fact that densities are lower than the critical density for coupling between the gas and dust phase. The k-means clustering method separated four sub-filaments in TMC-1 in the position-velocity-column density parameter space. They have masses of 32.5, 19.6, 28.9 and 45.9 M$_{\odot}$, low turbulent velocity dispersion (0.13-0.2 kms$^{-1}$) and they are close to gravitational equilibrium. We label them TMC-1F1 through F4. TMC-1F1, TMC-1F2 and TMC-1F4 are very elongated, dense and cold. TMC-1F3 is a little less elongated and somewhat warmer, probably heated by IRAS 04381+2540 that is embedded in it. TMC-1F3 is $\approx$ 0.1 pc behind TMC1-F1. Because of its structure, TMC-1 is a good target to test filament evolution scenarios.

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

Dust grains play a key role in the physics of star-forming regions, even though they constitute only $\sim$1 % of the mass of the interstellar medium. The derivation of accurate dust parameters such as temperature ($T_{dust}$), emissivity spectral index ($\beta$) and column density requires broadband continuum observations at far-infrared wavelengths. We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) measurements of three starless cores: L1521E, L1521F and L1689B, covering wavelengths between 194 and 671 $\mu$m. This paper is the first to use our recently updated SPIRE-FTS intensity calibration, yielding a direct match with SPIRE photometer measurements of extended sources. In addition, we carefully assess the validity of calibration schemes depending on source extent and on the strength of background emission. The broadband far-infrared spectra for all three sources peak near 250 $\mu$m. Our observations therefore provide much tighter constraints on the spectral energy distribution (SED) shape than measurements that do not probe the SED peak. The spectra are fitted using modified blackbody functions, allowing both $T_{dust}$ and $\beta$ to vary as free parameters. This yields $T_{dust}$ of 9.8$\pm$0.2 K, 15.6$\pm$0.5 K and 10.9$\pm$0.2 K and corresponding $\beta$ of 2.6$\mp$0.9, 0.8$\mp$0.1 and 2.4$\mp$0.8 for L1521E, L1521F and L1689B respectively. The derived core masses are 1.0$\pm$0.1, 0.10$\pm$0.01 and 0.49$\pm$0.05 $M_{\odot}$, respectively. The core mass/Jeans mass ratios for L1521E and L1689B exceed unity indicating that they are unstable to gravitational collapse, and thus pre-stellar cores. By comparison, the elevated temperature and gravitational stability of L1521F support previous arguments that this source is more evolved and likely a protostar.

Galactic plane surveys of pristine molecular clouds are key for establishing a Galactic-scale view of the earliest stages of star formation. For this reason Peretto & Fuller (2009) built an unbiased sample of IRDCs in the 10 deg < |l| < 65 deg, |b|<1 deg region of the Galactic plane using Spitzer 8micron extinction. However, in absorption studies, intrinsic fluctuations in the mid-infrared background can be mis-interpreted as foreground clouds. The main goal of the study presented here is to disentangle real clouds in the Spitzer Dark Cloud (SDC) catalogue from artefacts due to fluctuations in the mid-infrared background. We constructed H_2 column density maps at ~18" resolution using the 160micron and 250micron data from the Herschel Galactic plane survey Hi-GAL. We also developed an automated detection scheme that confirms the existence of a SDC through its association with a peak on these Herschel column density maps. Detection simulations, along with visual inspection of a small sub-sample of SDCs, have been performed to get better insight into the limitations of our automated identification scheme. Our analysis shows that 76(+/-19)% of the catalogued SDCs are real. This fraction drops to 55(+/-12)% for clouds with angular diameters larger than ~1 arcminute. The contamination of the PF09 catalogue by large spurious sources reflect the large uncertainties associated to the construction of the 8micron background emission, a key stage towards the identification of SDCs. A comparison of the Herschel confirmed SDC sample with the BGPS and ATLASGAL samples shows that SDCs probe a unique range of cloud properties, reaching down to more compact and lower column density clouds than any of these two (sub-)millimetre Galactic plane surveys.

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.

The surroundings of HII regions can have a profound influence on their development, morphology, and evolution. This paper explores the effect of the environment on H II regions in the MonR2 molecular cloud. We aim to investigate the density structure of envelopes surrounding HII regions and to determine their collapse and ionisation expansion ages. The Mon R2 molecular cloud is an ideal target since it hosts an H II region association. Column density and temperature images derived from Herschel data were used together to model the structure of HII bubbles and their surrounding envelopes. The resulting observational constraints were used to follow the development of the Mon R2 ionised regions with analytical calculations and numerical simulations. The four hot bubbles associated with H II regions are surrounded by dense, cold, and neutral gas envelopes. The radial density profiles are reminiscent of those of low-mass protostellar envelopes. The inner parts of envelopes of all four HII regions could be free-falling because they display shallow density profiles. As for their outer parts, the two compact HII regions show a density profile, which is typical of the equilibrium structure of an isothermal sphere. In contrast, the central UCHii region shows a steeper outer profile, that could be interpreted as material being forced to collapse. The size of the heated bubbles, the spectral type of the irradiating stars, and the mean initial neutral gas density are used to estimate the ionisation expansion time, texp, 0.1Myr,for the dense UCHII and compact HII regions and 0.35 Myr for the extended HII region. The envelope transition radii between the shallow and steeper density profiles are used to estimate the time elapsed since the formation of the first proto stellar embryo, Tinf : 1Myr, for the ultra-compact, 1.5 / 3Myr for the compact, and greater than 6Myr for the extended HII regions.

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.

Two major features of the prestellar CMF are: 1) a broad peak below 1 Msun, presumably corresponding to a mean gravitational fragmentation scale, and 2) a characteristic power-law slope, very similar to the Salpeter slope of the stellar initial mass function (IMF) at the high-mass end. While recent Herschel observations have shown that the peak of the prestellar CMF is close to the thermal Jeans mass in marginally supercritical filaments, the origin of the power-law tail of the CMF/IMF at the high-mass end is less clear. Inutsuka (2001) proposed a theoretical scenario in which the origin of the power-law tail can be understood as resulting from the growth of an initial spectrum of density perturbations seeded along the long axis of filaments by interstellar turbulence. Here, we report the statistical properties of the line-mass fluctuations of filaments in nearby molecular clouds observed with Herschel using a 1-D power spectrum analysis. The observed filament power spectra were fitted by a power-law function $(P_{true}(s) \propto s^{\alpha})$ after removing the effect of beam convolution at small scales. A Gaussian-like distribution of power-spectrum slopes was found centered at -1.6, close to that of the one-dimensional velocity power spectrum generated by subsonic Kolomogorov turbulence (-1.67). An empirical correlation, $P^{0.5}(s_0) \propto <N_{\rm H_2}>^{1.4 \pm 0.1} $, was also found between the amplitude of each filament power spectrum $P(s_0)$ and the mean column density along the filament $<N_{\rm H_2}>$. Finally, the dispersion of line-mass fluctuations along each filament $\sigma_{\rm M_{line}}$ was found to scale with the physical length $L$ of the filament, roughly as $\sigma_{M_{line}} \propto L^{0.7}$. Overall, our results are consistent with the suggestion that the bulk of the CMF/IMF results from the gravitational fragmentation of filaments.

We have studied the filaments extracted from the column density maps of the nearby Lupus 1, 3, and 4 molecular clouds, derived from photometric maps observed with the Herschel satellite. Filaments in the Lupus clouds have quite low column densities, with a median value of $\sim$1.5$\times$10$^{21}$ cm$^{-2}$ and most have masses per unit length lower than the maximum critical value for radial gravitational collapse. Indeed, no evidence of filament contraction has been seen in the gas kinematics. We find that some filaments, that on average are thermally subcritical, contain dense cores that may eventually form stars. This is an indication that in the low column density regime, the critical condition for the formation of stars may be reached only locally and this condition is not a global property of the filament. Finally, in Lupus we find multiple observational evidences of the key role that the magnetic field plays in forming filaments, and determining their confinement and dynamical evolution.

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

We show Akari data, Herschel data and data from the SCUBA2 camera on JCMT, of molecular clouds. We focus on pre-stellar cores within the clouds. We present Akari data of the L1147-1157 ring in Cepheus and show how the data indicate that the cores are being externally heated. We present SCUBA2 and Herschel data of the Ophiuchus region and show how the environment is also affecting core evolution in this region. We discuss the effects of the magnetic field in the Lupus I region, and how this lends support to a model for the formation and evolution of cores in filamentary molecular clouds.

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.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

We present a continuing study of a sample 44 molecular outflows, observed in 13CO lines, closely associated with 6.7GHz methanol masers, hence called Methanol Maser Associated Outflows (MMAOs). We compare MMAO properties with those of outflows from other surveys in the literature. In general, MMAOs follow similar trends, but show a deficit in number at low masses and momenta, with a corresponding higher fraction at the high end of the distributions. A similar trend is seen for the dynamical timescales of MMAOs. We argue that the lack of relatively low mass and young flows in MMAOs is due to the inherent selection-bias in the sample, i.e. its direct association with 6.7GHz methanol masers. This implies that methanol masers must switch on after the onset of outflows (hence accretion), and not before a sufficient abundance of methanol is liberated from icy dust mantles. Consequently the average dynamical age of MMAOs is older than for the general population of molecular outflows. We propose an adjusted evolutionary sequence of outflow and maser occurrence in the hot core phase, where methanol masers turn on after the onset of the outflow phase.

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.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was a suborbital experiment designed to map magnetic fields in order to study their role in star formation processes. BLASTPol made detailed polarization maps of a number of molecular clouds during its successful flights from Antarctica in 2010 and 2012. We present the next-generation BLASTPol instrument (BLAST-TNG) that will build off the success of the previous experiment and continue its role as a unique instrument and a test bed for new technologies. With a 16-fold increase in mapping speed, BLAST-TNG will make larger and deeper maps. Major improvements include a 2.5 m carbon fiber mirror that is 40% wider than the BLASTPol mirror and ~3000 polarization sensitive detectors. BLAST-TNG will observe in three bands at 250, 350, and 500 microns. The telescope will serve as a pathfinder project for microwave kinetic inductance detector (MKID) technology, as applied to feedhorn coupled submillimeter detector arrays. The liquid helium cooled cryostat will have a 28-day hold time and will utilize a closed-cycle $^3$He refrigerator to cool the detector arrays to 270 mK. This will enable a detailed mapping of more targets with higher polarization resolution than any other submillimeter experiment to date. BLAST-TNG will also be the first balloon-borne telescope to offer shared risk observing time to the community. This paper outlines the motivation for the project and the instrumental design.

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 have selected the positions of 54 6.7GHz methanol masers from the Methanol Multibeam Survey catalogue, covering a range of longitudes between $20^{\circ}$ and $34^{\circ}$ of the Galactic Plane. These positions were mapped in the J=3-2 transition of both the $\rm{^{13}CO}$ and $\rm{C^{18}O}$ lines. A total of 58 $\rm{^{13}CO}$ emission peaks are found in the vicinity of these maser positions. We search for outflows around all $\rm{^{13}CO}$ peaks, and find evidence for high-velocity gas in all cases, spatially resolving the red and blue outflow lobes in 55 cases. Of these sources, 44 have resolved kinematic distances, and are closely associated with the 6.7GHz masers, a sub-set referred to as Methanol Maser Associated Outflows (MMAOs). We calculate the masses of the clumps associated with each peak using 870 $\rm{\mu m}$ continuum emission from the ATLASGAL survey. A strong correlation is seen between the clump mass and both outflow mass and mechanical force, lending support to models in which accretion is strongly linked to outflow. We find that the scaling law between outflow activity and clump masses observed for low-mass objects, is also followed by the MMAOs in this study, indicating a commonality in the formation processes of low-mass and high-mass stars.

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 magnetic fields play in star formation. BLASTPol has had two science flights from McMurdo Station, Antarctica in 2010 and 2012. These flights have produced thousands of polarization vectors at 250, 350 and 500 microns in several molecular cloud targets. We present the design, specifications, and progress towards the next-generation BLASTPol experiment (BLAST-TNG). BLAST-TNG will fly a 40% larger diameter primary mirror, with almost 8 times the number of polarization-sensitive detectors resulting in a factor of 16 increase in mapping speed. With a spatial resolution of 22 arcseconds and four times the field of view of BLASTPol, BLAST-TNG will bridge the angular scales between Planck's low resolution all-sky maps and ALMA's ultra-high resolution narrow fields. The new receiver has a larger cryogenics volume, allowing for a 28 day hold time. BLAST-TNG employs three arrays of Microwave Kinetic Inductance Detectors (MKIDs) with 30% fractional bandwidth at 250, 350 and 500 microns. In this paper, we will present the new BLAST-TNG instrument and science objectives.

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 the results of Herschel HOBYS photometric mapping combined with BIMA observations and additional archival data, and perform an in-depth study of the evolutionary phases of the star-forming clumps in W 48A and their surroundings. Age estimates for the compact sources were derived from bolometric luminosities and envelope masses, which were obtained from the dust continuum emission, and agree within an order of magnitude with age estimates from molecular line and radio data. The clumps in W 48A are linearly aligned by age (east-old to west-young): we find a ultra compact (UC) HII region, a young stellar object (YSO) with class II methanol maser emission, a YSO with a massive outflow, and finally the NH_2D prestellar cores from Pillai et al. This remarkable positioning reflects the (star) formation history of the region. We find that it is unlikely that the star formation in the W 48A molecular cloud was triggered by the UCHII region and discuss the Aquila supershell expansion as a mayor influence on the evolution of W 48A. We conclude that the combination of Herschel continuum data with interferometric molecular line and radio continuum data is important to derive trustworthy age estimates and interpret the origin of large scale structures through kinematic information.

The density and temperature structures of dense cores in the L1495 cloud of the Taurus star-forming region are investigated using Herschel SPIRE and PACS images in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m and 500 $\mu$m continuum bands. A sample consisting of 20 cores, selected using spectral and spatial criteria, is analysed using a new maximum likelihood technique, COREFIT, which takes full account of the instrumental point spread functions. We obtain central dust temperatures, $T_0$, in the range 6-12 K and find that, in the majority of cases, the radial density falloff at large radial distances is consistent with the $r^{-2}$ variation expected for Bonnor-Ebert spheres. Two of our cores exhibit a significantly steeper falloff, however, and since both appear to be gravitationally unstable, such behaviour may have implications for collapse models. We find a strong negative correlation between $T_0$ and peak column density, as expected if the dust is heated predominantly by the interstellar radiation field. At the temperatures we estimate for the core centres, carbon-bearing molecules freeze out as ice mantles on dust grains, and this behaviour is supported here by the lack of correspondence between our estimated core locations and the previously-published positions of H$^{13}$CO$^+$ peaks. On this basis, our observations suggest a sublimation-zone radius typically $\sim 10^4$ AU. Comparison with previously-published N$_2$H$^+$ data at 8400 AU resolution, however, shows no evidence for N$_2$H$^+$ depletion at that resolution.

Recent studies of the nearest star-forming clouds of the Galaxy at submillimeter wavelengths with the Herschel Space Observatory have provided us with unprecedented images of the initial and boundary conditions of the star formation process. The Herschel results emphasize the role of interstellar filaments in the star formation process and connect remarkably well with nearly a decade's worth of numerical simulations and theory that have consistently shown that the ISM should be highly filamentary on all scales and star formation is intimately related to self-gravitating filaments. In this review, we trace how the apparent complexity of cloud structure and star formation is governed by relatively simple universal processes - from filamentary clumps to galactic scales. We emphasize two crucial and complementary aspects: (i) the key observational results obtained with Herschel over the past three years, along with relevant new results obtained from the ground on the kinematics of interstellar structures, and (ii) the key existing theoretical models and the many numerical simulations of interstellar cloud structure and star formation. We then synthesize a comprehensive physical picture that arises from the confrontation of these observations and simulations.

Utilizing multi-wavelength dust emission maps acquired with $Herschel$, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using inverse-Abel transform based technique. We present intrinsic radial dust temperature profiles of starless cores directly from dust continuum emission maps disentangling the effect of temperature variations along the line of sight which was previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in core center, a flat inner part, and a rising outward trend until the background cloud temperature is reached. The central beam-averaged dust temperatures obtained for B68 and L1689B are 9.3 $\pm$ 0.5 K and 9.8 $\pm$0.5 K, respectively, which are lower than the temperatures of 11.3 K and 11.6 K obtained from direct SED fitting. The best mass estimates derived by integrating the volume density profiles of B68 and L1689B are 1.6 M_sol and 11 M_sol, respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, $\kappa_{\lambda} =0.1(\lambda/300 \mu m)^{-2}$, agrees to within 50% with the dust extinction constraints

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.

We present the first overview of the Herschel observations of the nearby high-mass star-forming region NGC 7538, taken as part of the Herschel imaging study of OB Young Stellar objects (HOBYS) Key Programme. These PACS and SPIRE maps cover an approximate area of one square degree at five submillimeter and far-infrared wavebands. We have identified 780 dense sources and classified 224 of those. With the intention of investigating the existence of cold massive starless or class 0-like clumps that would have the potential to form intermediate- to high-mass stars, we further isolate 13 clumps as the most likely candidates for followup studies. These 13 clumps have masses in excess of 40 M_sun and temperatures below 15 K. They range in size from 0.4 pc to 2.5 pc and have densities between 3x10^3 cm^-3 to 4x10^4 cm^-3. Spectral energy distributions are then used to characterize their energetics and evolutionary state through a luminosity-mass diagram. NGC 7538 has a highly filamentary structure, previously unseen in the dust continuum of existing submillimeter surveys. We report the most complete imaging to date of a large, evacuated ring of material in NGC 7538 which is bordered by many cool sources.

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

A key parameter to the description of all star formation processes is the density structure of the gas. In this letter, we make use of probability distribution functions (PDFs) of Herschel column density maps of Orion B, Aquila, and Polaris, obtained with the Herschel Gould Belt survey (HGBS). We aim to understand which physical processes influence the PDF shape, and with which signatures. The PDFs of Orion B (Aquila) show a lognormal distribution for low column densities until Av 3 (6), and a power-law tail for high column densities, consistent with a rho r^-2 profile for the equivalent spherical density distribution. The PDF of Orion B is broadened by external compression due to the nearby OB stellar aggregates. The PDF of a quiescent subregion of the non-star-forming Polaris cloud is nearly lognormal, indicating that supersonic turbulence governs the density distribution. But we also observe a deviation from the lognormal shape at Av>1 for a subregion in Polaris that includes a prominent filament. We conclude that (i) the point where the PDF deviates from the lognormal form does not trace a universal Av-threshold for star formation, (ii) statistical density fluctuations, intermittency and magnetic fields can cause excess from the lognormal PDF at an early cloud formation stage, (iii) core formation and/or global collapse of filaments and a non-isothermal gas distribution lead to a power-law tail, and (iv) external compression broadens the column density PDF, consistent with numerical simulations.

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 of the mid-infrared environments and association with star formation tracers of 6.7 GHz methanol masers taken from the Methanol Multi-Beam (MMB) Survey. Our ultimate goal is to establish the mass of the host star and its evolutionary stage for each maser site. As a first step, the GLIMPSE survey of the Galactic Plane is utilised to investigate the environment of 776 methanol masers and we find that while the majority of the masers are associated with mid-infrared counterparts, a significant fraction (17%) are not associated with any detectable mid-infrared emission. A number of the maser counterparts are clearly extended with respect to the GLIMPSE point spread function and we implement an adaptive non-circular aperture photometry (ANCAP) technique to determine the fluxes of the maser counterparts. The ANCAP technique doubles the number of masers with flux information at all four wavelengths compared to the number of the corresponding counterparts obtained from the GLIMPSE Point Source Catalogue. The colours of the maser counterparts are found to be very similar to the smaller study carried out by Ellingsen (2006). The MMB masers are weakly associated with Extended Green Objects (EGOs) and Red MSX Survey (RMS) embedded sources (YSO and HII classifications) with 18% and 12% of masers associated with these objects respectively. The majority of MMB masers (60%) have detectable GLIMPSE infrared counterparts but have not been identified with previously recognised tracers of massive star formation; this confirms that the MMB survey has the potential to identify massive star forming regions independent of infrared selection.

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

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

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

It is not known whether brown dwarfs (stellar-like objects with masses less than the hydrogen-burning limit, 0.075 Msun) are formed in the same way as solar-type stars or by some other process. Here we report the clear-cut identification of a self-gravitating condensation of gas and dust with a mass in the brown-dwarf regime, made through millimeter interferometric observations. The level of thermal millimeter continuum emission detected from this object indicates a mass ~ 0.02-0.03 Msun, while the small radius < 460 AU and narrow spectral lines imply a dynamical mass of 0.015-0.02 Msun. The identification of such a pre-brown dwarf core supports models according to which brown dwarfs are formed in the same manner as hydrogen-burning stars.

Interstellar filaments are an important part of star formation. To understand the structure of filaments, cross-section profiles are often fitted with Plummer profiles. This profiling is often done with submm studies, such as Herschel. It would be convenient if filament properties could also be studied using groundbased NIR data. We compare the filament profiles obtained by NIR extinction and submm observations to find out if reliable profiles can be derived using NIR data. We use J-, H-, and K-band data of a filament north of TMC-1 to derive an extinction map from colour excesses of background stars. We compare the Plummer profiles obtained from extinction maps with Herschel dust emission maps. We present 2 methods to estimate profiles from NIR: Plummer profile fits to median Av of stars or directly to the Av of individual stars. We compare the methods by simulations. In simulations extinction maps and the new methods give correct results to within ~10-20 for modest densities. Direct fit to data on individual stars gives more accurate results than extinction map, and can work in higher density. In profile fits to real observations, values of Plummer parameters are generally similar to within a factor of ~2. Although parameter values can vary significantly, estimates of filament mass usually remain accurate to within some tens of per cent. Our results for TMC-1 are in agreement with earlier results. High resolution NIR data give more details, but 2MASS data can be used to estimate profiles. NIR extinction can be used as an alternative to submm observations to profile filaments. Direct fits of stars can also be a valuable tool. Plummer profile parameters are not always well constrained, and caution should be taken when making fits. In the evaluation of Plummer parameters, one can use the independence of dust emission and NIR data and the difference in the shapes of the confidence regions.

In order to characterise the cloud structures responsible for the formation of high-mass stars, we present Herschel observations of the DR21 environment. Maps of the column density and dust temperature unveil the structure of the DR21 ridge and several connected filaments. The ridge has column densities larger than 1e23/cm^2 over a region of 2.3 pc^2. It shows substructured column density profiles and branching into two major filaments in the north. The masses in the studied filaments range between 130 and 1400 Msun whereas the mass in the ridge is 15000 Msun. The accretion of these filaments onto the DR21 ridge, suggested by a previous molecular line study, could provide a continuous mass inflow to the ridge. In contrast to the striations seen in e.g., the Taurus region, these filaments are gravitationally unstable and form cores and protostars. These cores formed in the filaments potentially fall into the ridge. Both inflow and collisions of cores could be important to drive the observed high-mass star formation. The evolutionary gradient of star formation running from DR21 in the south to the northern branching is traced by decreasing dust temperature. This evolution and the ridge structure can be explained by two main filamentary components of the ridge that merged first in the south.

In a previous paper we identified cores within infrared dark clouds (IRDCs). We regarded those without embedded sources as the least evolved, and labelled them starless. Here we identify the most isolated starless cores and model them using a three-dimensional, multi-wavelength, Monte Carlo, radiative transfer code. We derive the cores' physical parameters and discuss the relation between the mass, temperature, density, size and the surrounding interstellar radiation field (ISRF) for the cores. The masses of the cores were found not to correlate with their radial size or central density. The temperature at the surface of a core was seen to depend almost entirely on the level of the ISRF surrounding the core. No correlation was found between the temperature at the centre of a core and its local ISRF. This was seen to depend, instead, on the density and mass of the core.

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

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

A growing body of evidence indicates that the formation of filaments in interstellar clouds is a key component of the star formation process. In this paper, we present new Herschel PACS and SPIRE observations of the B59 and Stem regions in the Pipe Nebula complex, revealing a rich, organized network of filaments. The asymmetric column density profiles observed for several filaments, along with the bow-like edge of B59, indicates that the Pipe Nebula is being compressed from its western side, most likely by the winds from the nearby Sco OB2 association. We suggest that this compressive flow has contributed to the formation of some of the observed filamentary structures. In B59, the only region of the entire Pipe complex showing star formation activity, the same compressive flow has likely enhanced the initial column density of the clump, allowing it to become globally gravitationally unstable. Although more speculative, we propose that gravity has also been responsible for shaping the converging filamentary pattern observed in B59. While the question of the relative impact of large-scale compression and gravity remains open in B59, large-scale compression appears to be a plausible mechanism for the initial formation of filamentary structures in the rest of the complex

As part of the Herschel guaranteed time key program 'HOBYS', we present the photometric survey of the star forming region Vela-C, one of the nearest sites of low-to-high-mass star formation in the Galactic plane. Vela-C has been observed with PACS and SPIRE in parallel mode between 70 um and 500 um over an area of about 3 square degrees. A photometric catalogue has been extracted from the detections in each band, using a threshold of 5 sigma over the local background. Out of this catalogue we have selected a robust sub-sample of 268 sources, of which 75% are cloud clumps and 25% are cores. Their Spectral Energy Distributions (SEDs) have been fitted with a modified black body function. We classify 48 sources as protostellar and 218 as starless. For two further sources, we do not provide a secure classification, but suggest they are Class 0 protostars. From SED fitting we have derived key physical parameters. Protostellar sources are in general warmer and more compact than starless sources. Both these evidences can be ascribed to the presence of an internal source(s) of moderate heating, which also causes a temperature gradient and hence a more peaked intensity distribution. Moreover, the reduced dimensions of protostellar sources may indicate that they will not fragment further. A virial analysis of the starless sources gives an upper limit of 90% for the sources gravitationally bound and therefore prestellar. We fit a power law N(logM) prop M^-1.1 to the linear portion of the mass distribution of prestellar sources. This is in between that typical of CO clumps and those of cores in nearby star-forming regions. We interpret this as a result of the inhomogeneity of our sample, which is composed of comparable fractions of clumps and cores.

We have used data taken as part of the Herschel infrared Galactic Plane survey (Hi-GAL) to study 3171 infrared-dark cloud (IRDC) candidates that were identified in the mid-infrared (8 {\mu}m) by Spitzer (we refer to these as 'Spitzer-dark' regions). They all lie in the range l=300 - 330 \circ and |b| 6 1 \circ. Of these, only 1205 were seen in emission in the far-infrared (250-500 {\mu}m) by Herschel (we call these 'Herschel-bright' clouds). It is predicted that a dense cloud will not only be seen in absorption in the mid-infrared, but will also be seen in emission in the far-infrared at the longest Herschel wavebands (250-500 {\mu}m). If a region is dark at all wavelengths throughout the mid-infrared and far-infrared, then it is most likely to be simply a region of lower background infrared emission (a 'hole in the sky'). Hence, it appears that previous surveys, based on Spitzer and other mid-infrared data alone, may have over-estimated the total IRDC population by a factor of 2. This has implications for estimates of the star formation rate in IRDCs in the Galaxy.We studied the 1205 Herschel-bright IRDCs at 250 {\mu}m, and found that 972 of them had at least one clearly defined 250-{\mu}m peak, indicating that they contained one or more dense cores. Of these, 653 (67 per cent) contained an 8-{\mu}m point source somewhere within the cloud, 149 (15 per cent) contained a 24-{\mu}m point source but no 8-{\mu}m source, and 170 (18 per cent) contained no 24-{\mu}m or 8-{\mu}m point sources. We use these statistics to make inferences about the lifetimes of the various evolutionary stages of IRDCs.

The filament IRDC G035.39--00.33 in the W48 molecular complex is one of the darkest infrared clouds observed by \textit{Spitzer}. It has been observed by the PACS (70 and 160\,$\micron$) and SPIRE (250, 350, and 500\,$\micron$) cameras of the \textit{Herschel} Space Observatory as part of the W48 molecular cloud complex in the framework of the HOBYS key programme. The observations reveal a sample of 28 compact sources (deconvolved FWHM sizes $<$0.3 pc) complete down to $\sim$$5 \msun$ in G035.39--00.33 and its surroundings. Among them, 13 compact sources are massive dense cores with masses $>$$20 \msun$. The cloud characteristics we derive from the analysis of their spectral energy distributions are masses of $20-50 \msun$, sizes of 0.1--0.2 pc, and average densities of $2-20 \times 10^{5} \cmc$, which make these massive dense cores excellent candidates to form intermediate- to high-mass stars. Most of the massive dense cores are located inside the G035.39--00.33 ridge and host IR-quiet high-mass protostars. The large number of protostars found in this filament suggests that we are witnessing a mini-burst of star formation with an efficiency of $\sim$15% and a rate density of $\sim$$40 \msun\,$yr$^{-1}\,$kpc$^{-2}$ within $\sim$8 pc$^2$, a large area covering the full ridge. Part of the extended SiO emission observed towards G035.39--00.33 is not associated with obvious protostars and may originate from low-velocity shocks within converging flows, as advocated by previous studies.

We present the first Herschel PACS and SPIRE results of the Vela C molecular complex in the far-infrared and submillimetre regimes at 70, 160, 250, 350, and 500 um, spanning the peak of emission of cold prestellar or protostellar cores. Column density and multi-resolution analysis (MRA) differentiates the Vela C complex into five distinct sub-regions. Each sub-region displays differences in their column density and temperature probability distribution functions (PDFs), in particular, the PDFs of the `Centre-Ridge' and `South-Nest' sub-regions appear in stark contrast to each other. The Centre-Ridge displays a bimodal temperature PDF representative of hot gas surrounding the HII region RCW 36 and the cold neighbouring filaments, whilst the South-Nest is dominated by cold filamentary structure. The column density PDF of the Centre-Ridge is flatter than the South-Nest, with a high column density tail, consistent with formation through large-scale flows, and regulation by self-gravity. At small to intermediate scales MRA indicates the Centre-Ridge to be twice as concentrated as the South-Nest, whilst on larger scales, a greater portion of the gas in the South-Nest is dominated by turbulence than in the Centre-Ridge. In Vela C, high-mass stars appear to be preferentially forming in ridges, i.e., dominant high column density filaments.

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.

Thermal images of cold dust in the Central Molecular Zone of the Milky Way, obtained with the far-infrared cameras on-board the Herschel satellite, reveal a 3x10^7 solar masses ring of dense and cold clouds orbiting the Galactic Center. Using a simple toy-model, an elliptical shape having semi-major axes of 100 and 60 parsecs is deduced. The major axis of this 100-pc ring is inclined by about 40 degrees with respect to the plane-of-the-sky and is oriented perpendicular to the major axes of the Galactic Bar. The 100-pc ring appears to trace the system of stable x_2 orbits predicted for the barred Galactic potential. Sgr A* is displaced with respect to the geometrical center of symmetry of the ring. The ring is twisted and its morphology suggests a flattening-ratio of 2 for the Galactic potential, which is in good agreement with the bulge flattening ratio derived from the 2MASS data.

We present high angular resolution observations, taken with the Very Large Array (VLA) and Multiple Element Radio Linked Interferometer Network (MERLIN) radio telescopes, at 7mm and 4.4cm respectively, of the prototype Class 0 protostar VLA1623. At 7mm we detect two sources (VLA1623A & B) coincident with the two previously detected components at the centre of this system. The separation between the two is 1.2arcsec, or ~170AU at an assumed distance of 139pc. The upper limit to the size of the source coincident with each component of VLA1623 is ~0.7arcsec, in agreement with previous findings. This corresponds to a diameter of ~100AU at an assumed distance of 139pc. Both components show the same general trend in their broadband continuum spectra, of a steeper dust continuum spectrum shortward of 7mm and a flatter spectrum longward of this. We estimate an upper limit to the VLA1623A disc mass of <0.13Msol and an upper limit to its radius of ~50AU. The longer wavelength data have a spectral index of \alpha~0.6+/-0.3. This is too steep to be explained by optically thin free-free emission. It is most likely due to optically thick free-free emission. Alternatively, we speculate that it might be due to the formation of larger grains or planetesimals in the circumstellar disc. We estimate the mass of VLA1623B to be <0.15M$sol. We can place a lower limit to its size of ~30x7 AU, and an upper limit to its diameter of ~100AU. The longer wavelength data of VLA1623B also have a spectral index of \alpha~0.6+/-0.3. The nature of VLA1623B remains a matter of debate. It could be a binary companion to the protostar, or a knot in the radio jet from VLA1623A.

We provide a first look at the results of the Herschel Gould Belt survey toward the IC5146 molecular cloud and present a preliminary analysis of the filamentary structure in this region. The column density map, derived from our 70-500 micron Herschel data, reveals a complex network of filaments, and confirms that these filaments are the main birth sites of prestellar cores. We analyze the column density profiles of 27 filaments and show that the underlying radial density profiles fall off as r^{-1.5} to r^{-2.5} at large radii. Our main result is that the filaments seem to be characterized by a narrow distribution of widths having a median value of 0.10 +- 0.03 pc, which is in stark contrast to a much broader distribution of central Jeans lengths. This characteristic width of ~0.1 pc corresponds to within a factor of ~2 to the sonic scale below which interstellar turbulence becomes subsonic in diffuse gas, supporting the argument that the filaments may form as a result of the dissipation of large-scale turbulence.

The densest infrared dark clouds (IRDCs) may represent the earliest observable stage of high-mass star formation. These clouds are very cold, hence they emit mainly at far-infrared and sub-mm wavelengths. For the first time, Herschel has provided multi-wavelength, spatially resolved observations of cores within IRDCs, which, when combined with radiative transfer modelling, can constrain their properties, such as mass, density profile and dust temperature. We use a 3D, multi-wavelength Monte Carlo radiative transfer code to model in detail the emission from six cores in three typical IRDCs seen in the Hi-GAL survey (G030.50+00.95, G031.03+00.26 and G031.03+00.76), and thereby to determine the properties of these cores and compare them with their low-mass equivalents. We found masses ranging from 90 to 290 solar masses with temperatures from 8 to 11K at the centre of each core and 18 to 28K at the surface. The maximum luminosity of an embedded star within each core was calculated, and we rule out the possibility of significant high mass star formation having yet occurred in three of our cores.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLAST-Pol) is a suborbital mapping experiment designed to study the role played by magnetic fields in the star formation process. BLAST-Pol is the reconstructed BLAST telescope, with the addition of linear polarization capability. Using a 1.8 m Cassegrain telescope, BLAST-Pol images the sky onto a focal plane that consists of 280 bolometric detectors in three arrays, observing simultaneously at 250, 350, and 500 um. The diffraction-limited optical system provides a resolution of 30'' at 250 um. The polarimeter consists of photolithographic polarizing grids mounted in front of each bolometer/detector array. A rotating 4 K achromatic half-wave plate provides additional polarization modulation. With its unprecedented mapping speed and resolution, BLAST-Pol will produce three-color polarization maps for a large number of molecular clouds. The instrument provides a much needed bridge in spatial coverage between larger-scale, coarse resolution surveys and narrow field of view, and high resolution observations of substructure within molecular cloud cores. The first science flight will be from McMurdo Station, Antarctica in December 2010.

We demonstrate the use of the 3D Monte Carlo radiative transfer code PHAETHON to model infrared-dark clouds (IRDCs) that are externally illuminated by the interstellar radiation field (ISRF). These clouds are believed to be the earliest observed phase of high-mass star formation, and may be the high-mass equivalent of lower-mass prestellar cores. We model three different cases as examples of the use of the code, in which we vary the mass, density, radius, morphology and internal velocity field of the IRDC. We show the predicted output of the models at different wavelengths chosen to match the observing wavebands of Herschel and Spitzer. For the wavebands of the long- wavelength SPIRE photometer on Herschel, we also pass the model output through the SPIRE simulator to generate output images that are as close as possible to the ones that would be seen using SPIRE. We then analyse the images as if they were real observations, and compare the results of this analysis with the results of the radiative transfer models. We find that detailed radiative transfer modelling is necessary to accurately determine the physical parameters of IRDCs (e.g. dust temperature, density profile). This method is applied to study G29.55+00.18, an IRDC observed by the Herschel Infrared Galactic Plane survey (Hi-GAL), and in the future it will be used to model a larger sample of IRDCs from the same survey.

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 summarize the first results from the Gould Belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demonstration phase' of Herschel. Our 70-500 micron images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-60 Class 0 protostars can be identified in the Aquila field, while ~ unbound starless cores and no protostars are observed in the Polaris field. The prestellar core mass function (CMF) derived for the Aquila region bears a strong resemblance to the stellar initial mass function (IMF), already confirming the close connection between the CMF and the IMF with much better statistics than earlier studies. Comparing and contrasting our Herschel results in Aquila and Polaris, we propose an observationally-driven scenario for core formation according to which complex networks of long, thin filaments form first within molecular clouds, and then the densest filaments fragment into a number of prestellar cores via gravitational instability.

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

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

We present a power spectrum analysis of the Herschel-SPIRE observations of the Polaris flare, a high Galactic latitude cirrus cloud midway between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the structure of the diffuse interstellar medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium is well described by a single power law with an exponent of -2.7 +- 0.1 at all scales from 30" to 8 degrees. That the power spectrum slope of the dust emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100 micron data of the Polaris flare clearly assesses the capability of SPIRE in maping diffuse emission over large areas.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

Our PACS and SPIRE images of the Aquila Rift and part of the Polaris Flare regions, taken during the science demonstration phase of Herschel discovered fascinating, omnipresent filamentary structures that appear to be physically related to compact cores. We briefly describe a new multi-scale, multi-wavelength source extraction method used to detect objects and measure their parameters in our Herschel images. All of the extracted starless cores (541 in Aquila and 302 in Polaris) appear to form in the long and very narrow filaments. With its combination of the far-IR resolution and sensitivity, Herschel directly reveals the filaments in which the dense cores are embedded; the filaments are resolved and have deconvolved widths of 35 arcsec in Aquila and 59 arcsec in Polaris (9000 AU in both regions). Our first results of observations with Herschel enable us to suggest that in general dense cores may originate in a process of fragmentation of complex networks of long, thin filaments, likely formed as a result of an interplay between gravity, interstellar turbulence, and magnetic fields. To unravel the roles of the processes, one has to obtain additional kinematic and polarization information; these follow-up observations are planned.

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star formation activity. These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to characterise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 Msun with luminosities between 1 and 150 Lsun, which extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared data are actually candidate Class 0 objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a single dense core, and they represent different evolutionary stages from Class 0 to Class I. Their mid- to far-infrared spectral slopes flatten towards the Class I stage, and the 160 to 70um flux ratio is greatest for the presumed Class 0 source. This shows that the Herschel observations characterise the earliest stages of protostellar evolution in detail.

The origin and possible universality of the stellar initial mass function (IMF) is a major issue in astrophysics. One of the main objectives of the Herschel Gould Belt Survey is to clarify the link between the prestellar core mass function (CMF) and the IMF. We present and discuss the core mass function derived from Herschel data for the large population of prestellar cores discovered with SPIRE and PACS in the Aquila Rift cloud complex at d ~ 260 pc. We detect a total of 541 starless cores in the entire ~11 deg^2 area of the field imaged at 70-500 micron with SPIRE/PACS. Most of these cores appear to be gravitationally bound, and thus prestellar in nature. Our Herschel results confirm that the shape of the prestellar CMF resembles the stellar IMF, with much higher quality statistics than earlier submillimeter continuum ground-based surveys.

We present far-infrared spectra and maps of the DR21 molecular cloud core between 196 and 671 microns, using the Herschel-SPIRE spectrometer. Nineteen molecular lines originating from CO, 13CO, HCO+ and H2O, plus lines of [N II] and [CI] were recorded, including several transitions not previously detected. The CO lines are excited in warm gas with Tkin ~ 125 K and nH2 ~ 7 x 10^4 cm-3, CO column density N(CO) ~ 3.5 x 10^18 cm^-2 and a filling factor of ~ 12%, and appear to trace gas associated with an outflow. The rotational temperature analysis incorporating observations from ground-based telescopes reveals an additional lower excitation CO compoment which has a temperature ~ 78 K and N(CO) ~ 4.5 x 10^21 cm^-2. Astronomy & Astrophysics HERSCHEL special Issue, in press.

Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index {\beta} and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 {\mu}m we determine T_dust and {\beta} throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of {\beta} ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between {\beta} and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.

The Polaris Flare cloud region contains a great deal of extended emission. It is at high declination and high Galactic latitude. It was previously seen strongly in IRAS Cirrus emission at 100 microns. We have detected it with both PACS and SPIRE on Herschel. We see filamentary and low-level structure. We identify the five densest cores within this structure. We present the results of a temperature, mass and density analysis of these cores. We compare their observed masses to their virial masses, and see that in all cases the observed masses lie close to the lower end of the range of estimated virial masses. Therefore, we cannot say whether they are gravitationally bound prestellar cores. Nevertheless, these are the best candidates to be potentialprestellar cores in the Polaris cloud region.

As part of the science demonstration phase of the Herschel mission of the Gould Belt Key Program, the Aquila Rift molecular complex has been observed. The complete ~ 3.3deg x 3.3deg imaging with SPIRE 250/350/500 micron and PACS 70/160 micron allows a deep investigation of embedded protostellar phases, probing of the dust emission from warm inner regions at 70 and 160 micron to the bulk of the cold envelopes between 250 and 500 micron. We used a systematic detection technique operating simultaneously on all Herschel bands to build a sample of protostars. Spectral energy distributions are derived to measure luminosities and envelope masses, and to place the protostars in an M_env - L_bol evolutionary diagram. The spatial distribution of protostars indicates three star-forming sites in Aquila, with W40/Sh2-64 HII region by far the richest. Most of the detected protostars are newly discovered. For a reduced area around the Serpens South cluster, we could compare the Herschel census of protostars with Spitzer results. The Herschel protostars are younger than in Spitzer with 7 Class 0 YSOs newly revealed by Herschel. For the entire Aquila field, we find a total of ~ 45-60 Class 0 YSOs discovered by Herschel. This confirms the global statistics of several hundred Class~0 YSOs that should be found in the whole Gould Belt survey.

Aims. To follow the species chemistry arising in diverse sources of the Galaxy with Herschel. Methods. SPIRE FTS sparse sampled maps of the Orion bar & compact HII regions G29.96-0.02 and G32.80+0.19 have been analyzed. Results. Beyond the wealth of atomic and molecular lines detected in the high-resolution spectra obtained with the FTS of SPIRE in the Orion Bar, one emission line is found to lie at the position of the fundamental rotational transition of CH+ as measured precisely in the laboratory (Pearson & Drouion 2006). This coincidence suggests that it is the first detection of the fundamental rotational transition of CH+. This claim is strengthened by the observation of the lambda doublet transitions arising from its relative, CH, which are also observed in the same spectrum. The broad spectral coverage of the SPIRE FTS allows for the simultaneous measurement of these closely related chemically species, under the same observing conditions. The importance of these lines are discussed and a comparison with results obtained from models of the Photon Dominated Region (PDR) of Orion are presented. The CH+ line also appears in absorption in the spectra of the two galactic compact HII regions G29.96-0.02 and G32.80+0.19, which is likely due to the presence of CH+ in the the Cold Neutral Medium of the galactic plane. These detections will shed light on the formation processes and on the existence of CH+, which are still outstanding questions in astrophysics.

We use the SPIRE Fourier-Transform Spectrometer (FTS) on-board the ESA Herschel Space Telescope to analyse the submillimetre spectrum of the Ultra-compact HII region G29.96-0.02. Spectral lines from species including 13CO, CO, [CI], and [NII] are detected. A sparse map of the [NII] emission shows at least one other HII region neighbouring the clump containing the UCHII. The FTS spectra are combined with ISO SWS and LWS spectra and fluxes from the literature to present a detailed spectrum of the source spanning three orders of magnitude in wavelength. The quality of the spectrum longwards of 100 {\mu}m allows us to fit a single temperature greybody with temperature 80.3\pm0.6K and dust emissivity index 1.73\pm0.02, an accuracy rarely obtained with previous instruments. We estimate a mass of 1500 Msol for the clump containing the HII region. The clump's bolometeric luminosity of 4 x 10^6 Lsol is comparable to, or slightly greater than, the known O-star powering the UCHII region.

Context. RCW 120 is a well-studied, nearby Galactic HII region with ongoing star formation in its surroundings. Previous work has shown that it displays a bubble morphology at mid-infrared wavelengths and has a massive layer of collected neutral material seen at sub-mm wavelengths. Given the well-defined photo-dissociation region (PDR) boundary and collected layer, it is an excellent laboratory to study the "collect and collapse" process of triggered star formation. Using Herschel Space Observatory data at 100, 160, 250, 350, and 500 micron, in combination with Spitzer and APEX-LABOCA data, we can for the first time map the entire spectral energy distribution of an HII region at high angular resolution. Aims. We seek a better understanding of RCW120 and its local environment by analysing its dust temperature distribution. Additionally, we wish to understand how the dust emissivity index, beta, is related to the dust temperature. Methods. We determine dust temperatures in selected regions of the RCW 120 field by fitting their spectral energy distribution (SED), derived using aperture photometry. Additionally, we fit the SED extracted from a grid of positions to create a temperature map. Results. We find a gradient in dust temperature, ranging from >30 K in the interior of RCW 120, to ~20K for the material collected in the PDR, to ~10K toward local infrared dark clouds and cold filaments. Our results suggest that RCW 120 is in the process of destroying the PDR delineating its bubble morphology. The leaked radiation from its interior may influence the creation of the next generation of stars. We find support for an anti-correlation between the fitted temperature and beta, in rough agreement with what has been found previously. The extended wavelength coverage of the Herschel data greatly increases the reliability of this result.

Infrared dark clouds (IRDCs) are cold and dense reservoirs of gas potentially available to form stars. Many of these clouds are likely to be pristine structures representing the initial conditions for star formation. The study presented here aims to construct and analyze accurate column density and dust temperature maps of IRDCs by using the first Herschel data from the Hi-GAL galactic plane survey. These fundamental quantities, are essential for understanding processes such as fragmentation in the early stages of the formation of stars in molecular clouds. We have developed a simple pixel-by-pixel SED fitting method, which accounts for the background emission. By fitting a grey-body function at each position, we recover the spatial variations in both the dust column density and temperature within the IRDCs. This method is applied to a sample of 22 IRDCs exhibiting a range of angular sizes and peak column densities. Our analysis shows that the dust temperature decreases significantly within IRDCs, from background temperatures of 20-30 K to minimum temperatures of 8-15 K within the clouds, showing that dense molecular clouds are not isothermal. Temperature gradients have most likely an important impact on the fragmentation of IRDCs. Local temperature minima are strongly correlated with column density peaks, which in a few cases reach NH2 = 1 x 10^{23} cm^{-2}, identifying these clouds as candidate massive prestellar cores. Applying this technique to the full Hi-GAL data set will provide important constraints on the fragmentation and thermal properties of IRDCs, and help identify hundreds of massive prestellar core candidates.

It has been shown that by means of different physical mechanisms the expansion of HII regions can trigger the formation of new stars of all masses. This process may be important to the formation of massive stars but has never been quantified in the Galaxy. We use Herschel-PACS and -SPIRE images from the Herschel Infrared survey of the Galactic plane, Hi-GAL, to perform this study. We combine the Spitzer-GLIMPSE and -MIPSGAL, radio-continuum and sub-millimeter surveys such as ATLASGAL with Hi-GAL to study Young Stellar Objects (YSOs) observed towards Galactic HII regions. We select a representative HII region, N49, located in the field centered on l=30 degr observed as part of the Hi-GAL Science Demonstration Phase, to demonstrate the importance Hi-GAL will have to this field of research. Hi-GAL PACS and SPIRE images reveal a new population of embedded young stars, coincident with bright ATLASGAL condensations. The Hi-GAL images also allow us, for the first time, to constrain the physical properties of the newly formed stars by means of fits to their spectral energy distribution. Massive young stellar objects are observed at the borders of the N49 region and represent second generation massive stars whose formation has been triggered by the expansion of the ionized region. Hi-GAL enables us to detect a population of young stars at different evolutionary stages, cold condensations only being detected in the SPIRE wavelength range. The far IR coverage of Hi-GAL strongly constrains the physical properties of the YSOs. The large and unbiased spatial coverage of this survey offers us a unique opportunity to lead, for the first time, a global study of star formation triggered by HII regions in our Galaxy.

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

We have conducted a Galactic plane survey of methanol masers at 6668 MHz using a 7-beam receiver on the Parkes telescope. Here we present results from the first part, which provides sensitive unbiased coverage of a large region around the Galactic Centre. Details are given for 183 methanol maser sites in the longitude range 345$^{\circ}$ through the Galactic Centre to 6$^{\circ}$. Within 6$^{\circ}$ of the Centre, we found 88 maser sites, of which more than half (48) are new discoveries. The masers are confined to a narrow Galactic latitude range, indicative of many sources at the Galactic Centre distance and beyond, and confined to a thin disk population; there is no high latitude population that might be ascribed to the Galactic Bulge. Within 2$^{\circ}$ of the Galactic Centre the maser velocities all lie between -60 and +77 \kms, a range much smaller than the 540 \kms range observed in CO. Elsewhere, the maser with highest positive velocity (+107 \kms) occurs, surprisingly, near longitude 355$^{\circ}$ and is probably attributable to the Galactic bar. The maser with the most negative velocity (-127 \kms) is near longitude 346$^{\circ}$, within the longitude-velocity locus of the near side of the `3-kpc arm'. It has the most extreme velocity of a clear population of masers associated with the near and far sides of the 3-kpc arm. Closer to the Galactic Centre the maser space density is generally low, except within 0.25 kpc of the Centre itself, the `Galactic Centre Zone', where it is 50 times higher, which is hinted at by the longitude distribution, and confirmed by the unusual velocities.

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 have used the Berkeley-Illinois-Maryland Array (BIMA) to make deep N2H+ 1-0 maps of the pre-stellar core L183, in order to study the spatial and kinematic substructure within the densest part of the core. Three spatially and kinematically distinct clumps are detected, which we label L183-N1, L183-N2 and L183-N3. L183-N2 is approximately coincident with the submillimetre dust peak and lies at the systemic velocity of L183. Thus we conclude that L183-N2 is the central dense core of L183. L183-N1 and 3 are newly-discovered fragments of L183, which are marked by velocity gradients that are parallel to, but far stronger than, the velocity gradient of L183 as a whole, as detected in previous single-dish data. Furthermore, the ratio of the large-scale and small-scale velocity gradients, and the ratio of their respective size-scales, are consistent with the conservation of angular momentum for a rotating, collapsing core undergoing spin-up. The inferred axis of rotation is parallel to the magnetic field direction, which is offset from its long axis, as we have seen in other pre-stellar cores. Therefore, we propose that we have detected a fragmenting, collapsing, filamentary, pre-stellar core, rotating about its B-field, which is spinning up as it collapses. It will presumably go on to form a multiple protostellar system.

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 Spitzer IRAC (~2 deg^2) and MIPS (~8 deg^2) observations of the Cepheus Flare which is associated with the Gould Belt, at an approximate distance of ~300 pc. Around 6500 sources are detected in all four IRAC bands, of which ~900 have MIPS 24 micron detections. We identify 133 YSO candidates using color-magnitude diagram techniques, a large number of the YSO candidates are associated with the NGC 7023 reflection nebula. Cross identifications were made with the Guide Star Catalog II and the IRAS Faint Source Catalog, and spectral energy distributions (SED) were constructed. SED modeling was conducted to estimate the degree of infrared excess. It was found that a large majority of disks were optically thick accreting disks, suggesting that there has been little disk evolution in these sources. Nearest-neighbor clustering analysis identified four small protostellar groups (L1228, L1228N, L1251A, and L1251B) with 5-15 members each and the larger NGC 7023 association with 32 YSO members. The star formation efficiency for cores with clusters of protostars and for those without clusters was found to be ~8% and ~1% respectively. The cores L1155, L1241, and L1247 are confirmed to be starless down to our luminosity limit of L_bol=0.06 L_sol.

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.

A new 7-beam 6-7 GHz receiver has been built to survey the Galaxy and the Magellanic Clouds for newly forming high-mass stars that are pinpointed by strong methanol maser emission at 6668 MHz. The receiver was jointly constructed by Jodrell Bank Observatory (JBO) and the Australia Telescope National Facility (ATNF) and allows simultaneous coverage at 6668 and 6035 MHz. It was successfully commissioned at Parkes in January 2006 and is now being used to conduct the Parkes-Jodrell multibeam maser survey of the Milky Way. This will be the first systematic survey of the entire Galactic plane for masers of not only 6668-MHz methanol, but also 6035-MHz excited-state hydroxyl. The survey is two orders of magnitude faster than most previous systematic surveys and has an rms noise level of ~0.17 Jy.This paper describes the observational strategy, techniques and reduction procedures of the Galactic and Magellanic Cloud surveys, together with deeper, pointed, follow-up observations and complementary observations with other instruments. It also includes an estimate of the survey detection efficiency. The 111 days of observations with the Parkes telescope have so far yielded >800 methanol sources, of which ~350 are new discoveries. The whole project will provide the first comprehensive Galaxy-wide catalogue of 6668-MHz and 6035-MHz masers.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a sub-orbital experiment designed to study the process of star formation in local galaxies (including the Milky Way) and in galaxies at cosmological distances. Using a 2-m Cassegrain telescope, BLAST images the sky onto a focal plane, which consists of 270 bolometric detectors split between three arrays, observing simultaneously in 30% wide bands, centered at 250, 350, and 500 microns. The diffraction-limited optical system provides a resolution of 30" at 250 microns. The pointing system enables raster-like scans with a positional accuracy of ~30", reconstructed to better than 5" rms in post-flight analysis. BLAST had two successful flights, from the Arctic in 2005, and from Antarctica in 2006, which provided the first high-resolution and large-area (~0.8-200 deg^2) submillimeter surveys at these wavelengths. As a pathfinder for the SPIRE instrument on Herschel, BLAST shares with the ESA satellite similar focal plane technology and scientific motivation. A third flight in 2009 will see the instrument modified to be polarization-sensitive (BLAST-Pol). With its unprecedented mapping speed and resolution, BLAST-Pol will provide insights into Galactic star-forming nurseries, and give the necessary link between the larger, coarse resolution surveys and the narrow, resolved observations of star-forming structures from space and ground based instruments being commissioned in the next 5 years.

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 the results of the first complete survey of the Large and Small Magellanic Clouds for 6668-MHz methanol and 6035-MHz excited-state hydroxyl masers. In addition to the survey, higher-sensitivity targeted searches towards known star-formation regions were conducted. The observations yielded the discovery of a fourth 6668-MHz methanol maser in the Large Magellanic Cloud (LMC), found towards the star-forming region N160a, and a second 6035-MHz excited-state hydroxyl maser, found towards N157a. We have also re-observed the three previously known 6668-MHz methanol masers and the single 6035-MHz hydroxyl maser. We failed to detect emission from either transition in the Small Magellanic Cloud. All observations were initially made using the Methanol Multibeam (MMB) survey receiver on the 64-m Parkes telescope as part of the MMB project and accurate positions have been measured with the Australia Telescope Compact Array (ATCA). We compare the maser populations in the Magellanic Clouds with those of our Galaxy and discuss their implications for the relative rates of massive star-formation, heavy metal abundance, and the abundance of complex molecules. The LMC maser populations are demonstrated to be smaller than their Milky Way counterparts. Methanol masers are under-abundant by a factor of ~45, whilst hydroxyl and water masers are a factor of ~10 less abundant than our Galaxy.

We report the discovery of a nearby, embedded cluster of young stellar objects, associated filamentary infrared dark cloud, and 4.5 micron shock emission knots from outflows detected in Spitzer/IRAC mid-infrared imaging of the Serpens-Aquila Rift obtained as part of the Spitzer Gould Belt Legacy Survey. We also present radial velocity measurements of the region from molecular line observations obtained with the Submillimeter Array (SMA) that suggest the cluster is co-moving with the Serpens Main embedded cluster 3 degrees to the north. We therefore assign it the same distance, 260 pc. The core of the new cluster, which we call Serpens South, is composed of an unusually large fraction of protostars (77%) at high mean surface density (>430 pc^-2) and short median nearest neighbor spacing (3700 AU). We perform basic cluster structure characterization using nearest neighbor surface density mapping of the YSOs and compare our findings to other known clusters with equivalent analyses available in the literature.

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.

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc).

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 estimate the dust temperatures of the clumps in the rho Oph main cloud taking into account the 3D geometry of the region, and external heating from the interstellar radiation field and from HD147879,a nearby luminous B2V star, which is believed to dominate the radiation field in the region. We find that the regions where prestellar cores are observed (i.e. at optical visual extinctions >7 mag) are colder than ~10-11K. These dust temperatures are smaller than those which previous studies of the same region have assumed. We use the new dust temperatures to estimate the masses of the prestellar cores in the rho Oph main cloud from mm observations, and we find core masses that are larger than previous estimates by a factor of ~2-3. This affects the core mass function (CMF) of the region; we find that the mass at which the core mass spectrum steepens from a slope alpha~1.5 to a slope alpha~2.5 has moved from ~0.5 Msun to ~1Msun. In contrast with the CMF in other star forming regions (e.g. Orion), there is no indication for a turnover down to the completeness limit (~0.2Msun), but the CMF may flatten at around ~0.4Msun. We generalize our results to the prestellar cores in Taurus and in Orion. In Taurus the ambient radiation field heating the prestellar cores is believed to be weaker than than that in rho Oph. Hence, the dust temperatures of the cores in Taurus are expected to be below ~ 10-11 K.In Orion the radiation field is believed to be 10^3 times stronger than the standard interstellar radiation field. Based on this assumption we estimate that the dust temperatures of the prestellar cores in Orion are around ~20-30K.

Four dense cores, L1582A, L1689A, B133 and B68, classified as prestellar in terms of the absence of detectable NIR emission, are observed at radio wavelengths to investigate whether they nurture very young protostars. No definite young protostars were discovered in any of the four cores observed. A few radio sources were discovered close to the observed cores, but these are most likely extragalactic sources or YSOs unrelated to the cores observed. In L1582A we discovered a weak radio source near the centre of the core with radio characteristics and offset from the peak of the submillimeter emission similar to that of the newly discovered protostar in the core L1014, indicating a possible protostellar nature for this source. This needs to be confirmed with near- and/or mid-infrared observations (e.g. with Spitzer). Hence based on the current observations we are unable to confirm unequivocally that L1582A is starless. In L1689A a possible 4.5-sigma radio source was discovered at the centre of the core, but needs to be confirmed with future observations. In B133 a weak radio source, possibly a protostar, was discovered at the edge of the core on a local peak of the core submm emission, but no source was detected at the centre of the core. Thus, B133 is probably starless, but may have a protostar at its edge. In B68 no radio sources were discovered inside or at the edge of the core, and thus B68 is indeed starless. Four more radio sources with spectral indices characteristic of young protostars were discovered outside the cores but within the extended clouds in which these cores reside. Conclusions: We conclude that the number of cores misclassified as prestellar is probably very small and does not significantly alter the estimated lifetime of the prestellar phase.

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 have retrieved Spitzer archive data of pre-stellar cores taken with the Multiband Imaging Photometer for Spitzer (MIPS) at a wavelength of 160 um. Seventeen images, containing eighteen cores, were constructed. Flux densities were measured for each core, and background estimates were made. Mean off-source backgrounds were found to be 48\pm10 MJy/sr in Taurus and 140\pm55 MJy/sr in Ophiuchus. Consistency was found between the MIPS 170um and ISOPHOT 160um calibration. Fourteen cores were detected both by MIPS and our previous submillimetre surveys. Spectral Energy Distributions (SEDs) were made for each core, using additional 24- and 70-um data from the Spitzer data archive, as well as previous infra-red and submillimetre data. Previous temperature estimates were refined, and new temperature estimates were made where no Infrared Space Observatory (ISO) data exist. A temperature range of 8-18K was found for the cores, with most lying in the range 10-13K. We discount recent claims that a large number of pre-stellar cores may have been misclassified and in fact contain low luminosity protostars detectable only by Spitzer. We find no new protostars in our sample other than that previously reported in L1521F. It is shown that this has a negligible effect on pre-stellar lifetime estimates.

We have mapped linearly polarized dust emission from the prestellar cores L1498 and L1517B with the James Clerk Maxwell Telescope (JCMT) using the Submillimetre Common User Bolometer Array (SCUBA) and its polarimeter SCUBAPOL at a wavelength of 850um. We use these measurements to determine the plane-of-sky magnetic field orientation in the cores. In L1498 we see a magnetic field across the peak of the core that lies at an offset of 19 degrees to the short axis of the core. This is similar to the offsets seen in previous observations of prestellar cores. To the southeast of the peak, in the filamentary tail of the core, we see that the magnetic field has rotated to lie almost parallel to the long axis of the filament. We hypothesise that the field in the core may have decoupled from the field in the filament that connects the core to the rest of the cloud. We use the Chandrasekhar-Fermi (CF) method to measure the plane-of-sky field strength in the core of L1498 to be 10 +/- 7 uG. In L1517B we see a more gradual turn in the field direction from the northern part of the core to the south. This appears to follow a twist in the filament in which the core is buried, with the field staying at a roughly constant 25 degree offset to the short axis of the filament, also consistent with previous observations of prestellar cores. We again use the CF method and calculate the magnetic field strength in L1517B also to be 30 +/- 10 uG. Both cores appear to be roughly virialised. Comparison with our previous work on somewhat denser cores shows that, for the denser cores, thermal and non-thermal (including magnetic) support are approximately equal, while for the lower density cores studied here, thermal support dominates.

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 review the properties of low mass dense molecular cloud cores, including starless, prestellar, and Class 0 protostellar cores, as derived from observations. In particular we discuss them in the context of the current debate surrounding the formation and evolution of cores. There exist several families of model scenarios to explain this evolution (with many variations of each) that can be thought of as a continuum of models lying between two extreme paradigms for the star and core formation process. At one extreme there is the dynamic, turbulent picture, while at the other extreme there is a slow, quasi-static vision of core evolution. In the latter view the magnetic field plays a dominant role, and it may also play some role in the former picture. Polarization and Zeeman measurements indicate that some, if not all, cores contain a significant magnetic field. Wide-field surveys constrain the timescales of the core formation and evolution processes, as well as the statistical distribution of core masses. The former indicates that prestellar cores typically live for 2--5 free-fall times, while the latter seems to determine the stellar initial mass function. In addition, multiple surveys allow one to compare core properties in different regions. From this it appears that aspects of different models may be relevant to different star-forming regions, depending on the environment. Prestellar cores in cluster-forming regions are smaller in radius and have higher column densities, by up to an order of magnitude, than isolated prestellar cores. This is probably due to the fact that in cluster-forming regions the prestellar cores are formed by fragmentation of larger, more turbulent cluster-forming cores, which in turn form as a result of strong external compression.

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 investigate the effect of the turbulent power spectrum (P(k) \propto k^{-n}, with n=3, 4 or 5) on the fragmentation of low-mass cores, by means of SPH simulations. We adopt initial density profiles and low levels of turbulence based on observation, and for each n-value we conduct an ensemble of simulations with different initial seeds for the turbulent velocity field, so as to obtain reasonable statistics. We find that when power is concentrated at larger scales (i.e. for larger n), more protostellar objects form and there is a higher proportion of low-mass stars and brown dwarfs. This is in direct contrast with the recent results of Delgado Donate et al., presumably because they adopted much higher levels of turbulence.

We present observations that were carried out with the Two HUndred Micron PhotometER (THUMPER) mounted on the James Clerk Maxwell Telescope (JCMT) in Hawaii, at a wavelength of 200 um (frequency 1.5 THz). The observations utilise a small atmospheric window that opens up at this wavelength under very dry conditions at high-altitude observing sites. The atmosphere was calibrated using the sky-dipping method and a relation was established between the optical depth, tau, at 1.5 THz and that at 225 GHz: tau_1.5THz = (95 +/- 10)*tau_225GHz. Mars and Jupiter were mapped from the ground at this wavelength for the first time, and the system characteristics measured. A noise equivalent flux density (NEFD) of ~65 +/- 10 Jy (1 sigma 1 second) was measured for the THUMPER-JCMT combination, consistent with predictions based upon our laboratory measurements. The main-beam resolution of 14 arcsec was confirmed and an extended error-beam detected at roughly two-thirds of the magnitude of the main beam. Measurements of the Sun allow us to estimate that the fraction of the power in the main beam is ~15%, consistent with predictions based on modelling the dish surface accuracy. It is therefore shown that the sky over Mauna Kea is suitable for astronomy at this wavelength under the best conditions. However, higher or drier sites should have a larger number of useable nights per year.

Observations have been carried out with SCUBA at the JCMT of 52 molecular cloud cores that do not contain any sign of protostellar activity. These are all therefore candidate prestellar cores, which are believed to represent the stage of star formation that precedes the formation of a protostar. 29 of the 52 cores were detected at 850 microns at varying levels of signal-to-noise ratio greater than 3 sigma at peak. The detected cores were split into 'bright' cores and `intermediate' cores, depending on their peak flux density at 850 microns. Cores with peak 850 microns flux densities greater than 170 mJy/beam were designated 'bright' (13 cores), while those flux densities below this value were designated 'intermediate' (16 cores). This dividing line corresponds to A_v~50 under typical assumptions. The data are combined with our previously published ISO data, and the physical parameters of the cores, such as density and temperature, are calculated. Detailed fitting of the bright core radial profiles shows that they are not critical Bonnor-Ebert spheres, in agreement with previous findings. However, we find that intermediate cores, such as B68 (which has previously been claimed to be a Bonnor-Ebert sphere), may in fact be consistent with the Bonnor-Ebert criterion, suggesting perhaps that cores pass through such a phase during their evolution. We make rough estimates of core lifetimes based on the statistics of detections and find that the lifetime of a prestellar core is roughly ~3x10^5 years, while that of a bright core is \~1.5x10^5 years. Comparisons with some magnetic and turbulence regulated collapse models show that no model can match all of the data. Models that are tuned to fit the total prestellar core lifetime, do not predict the relative numbers of cores seen at each stage.

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 examine the initial mass functions (IMFs) of stars produced by different molecular core mass functions. Simulations suggest that more massive cores produce more stars, so we propose a model in which the average number of stars formed in a core is equal to the initial number of Jeans masses in that core. Small-$N$ systems decay through dynamical interactions, ejecting low-mass stars and brown dwarfs which populate the low-mass tail of the IMF. Stars which remain in cores are able to competitively accrete more gas and become more massive. We deduce the forms of the core mass functions required to explain the IMFs of Taurus, Orion, IC 348 and NGC 2547. These core mass functions fall into two categories - one which peaks at a few $M_\odot$ to explain Taurus and NGC 2547, and one that peaks at around $0.2 M_\odot$ to explain Orion and IC 348.

(Abridged) We explore, by means of a large ensemble of SPH simulations, how the level of turbulence affects the collapse and fragmentation of a star-forming core. All our simulated cores have the same, except that we vary (a) the initial level of turbulence (as measured by the ratio of turbulent to gravitational energy, $\alpha_{\rm turb} \equiv U_{\rm turb}/|\Omega| = 0, 0.01, 0.025, 0.05, 0.10 {\rm and} 0.25$) and (b), for fixed $\alpha_{\rm turb}$, the details of the initial turbulent velocity field (so as to obtain good statistics). A low level of turbulence ($\alpha_{\rm turb} \sim 0.05$) suffices to produce multiple systems. As $\alpha_{\rm turb}$ is increased, the number of objects formed and the companion frequency both increase. The mass function is bimodal, with a flat low-mass segment representing single objects ejected from the core before they can accrete much, and a Gaussian high-mass segment representing objects which because they remain in the core grow by accretion and tend to pair up in multiple systems.

We present 2D Monte Carlo radiative transfer simulations of prestellar cores. We consider two types of asymmetry: disk-like asymmetry, in which the core is denser towards the equatorial plane than towards the poles; and axial asymmetry, in which the core is denser towards the south pole than the north pole. We limit our treatment to cores with mild asymmetries, which are exposed directly to the interstellar radiation field or are embedded inside molecular clouds. The isophotal maps of a core depend strongly on the viewing angle. Maps at wavelengths longer than the peak of the SED (e.g. 850 micron) essentially trace the column-density. Thus, for instance, cores with disk-like asymmetry appear elongated when mapped at 850 micron from close to the equatorial plane. However, at wavelengths near the peak of the SED (e.g. 200 micron), the emissivity is more strongly dependent on the temperature, and therefore, at particular viewing angles, there are characteristic features which reflect a more complicated convolution of the density and temperature fields within the core. These characteristic features are on scales 1/5 to 1/3 of the overall core size, and so high resolution observations are needed to observe them. They are also weaker if the core is embedded in a molecular cloud (because the range of temperature within the core is then smaller), and so high sensitivity is needed to detect them. Herschel, to be launched in 2007, will in principle provide the necessary resolution and sensitivity at 170 to 250 micron.

In comparison with other well studied star formation regions, Taurus is unusual in several respects. (i) Its stellar initial mass function (IMF) peaks at relatively high mass ($\sim 0.8 M_\odot$), but contains very few stars much more massive than $1 M_\odot$, and is relatively deficient in brown dwarfs. (ii) It has a higher binary fraction, particularly at large separations. (iii) Its core mass function is strongly peaked at a few $M_\odot$, and the cores have extended envelopes and relatively low levels of turbulence. We present here the results of an ensemble of hydrodynamic simulations which suggest that the unusual stellar IMF in Taurus is a direct consequence of the unusual properties of the cores there. By following the collapse and fragmentation of cores having properties typical of Taurus, we find that roughly 50% of the objects formed in a core, predominantly the low-mass ones, are ejected from the core to form a population of low-mass stars and brown dwarfs with a flat mass function. The remaining objects form multiple systems within the core, accreting until their masses approach $1 M_\odot$; this produces a population of intermediate-mass stars whose mass function peaks at $\sim 0.8 M_\odot$. Together these two populations reproduce the IMF in Taurus very well. This demonstrates, for the first time, a direct causal link between the core mass function and the stellar IMF in a star formation region.

Improving our understanding of the earliest stages of star formation is crucial to gain insight into the origin of stellar masses, multiple systems, and protoplanetary disks. We discuss recent advances made in this area thanks to detailed mapping observations at infrared and (sub)millimeter wavelengths. Although ambipolar diffusion appears to be too slow to play a direct role in the formation of dense cores, there is nevertheless good evidence that the gravitational collapse of isolated protostellar cores is strongly magnetically controlled. We also argue that the beginning of protostellar collapse is much more violent in cluster-forming clouds than in regions of distributed star formation.

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.

We present the results of an ensemble of simulations of the collapse and fragmentation of dense star-forming cores. We show that even with very low levels of turbulence the outcome is usually a binary, or higher-order multiple, system. We take as the initial conditions for these simulations a typical low-mass core, based on the average properties of a large sample of observed cores. All the simulated cores start with a mass of $M = 5.4 M_{\odot}$, a flattened central density profile, a ratio of thermal to gravitational energy $\alpha_{\rm therm} = 0.45$ and a ratio of turbulent to gravitational energy $\alpha_{\rm turb} = 0.05 $. Even this low level of turbulence is sufficient to produce multiple star formation in 80% of the cores; the mean number of stars and brown dwarfs formed from a single core is 4.55, and the maximum is 10. At the outset, the cores have no large-scale rotation. The only difference between each individual simulation is the detailed structure of the turbulent velocity field. The multiple systems formed in the simulations have properties consistent with observed multiple systems. Dynamical evolution tends preferentially to eject lower mass stars and brown dwarves whilst hardening the remaining binaries so that the median semi-major axis of binaries formed is $\sim 30$ au. Ejected objects are usually single low-mass stars and brown dwarfs, yielding a strong correlation between mass and multiplicity. Our simulations suggest a natural mechanism for forming binary stars that does not require large-scale rotation, capture, or large amounts of turbulence.

We evaluate the intrinsic three dimensional shapes of molecular cores, by analysing their projected shapes. We use the recent catalogue of molecular line observations of Jijina et al. and model the data by the method originally devised for elliptical galaxies. Our analysis broadly supports the conclusion of Jones et al. that molecular cores are better represented by triaxial intrinsic shapes (ellipsoids) than biaxial intrinsic shapes (spheroids). However, we find that the best fit to all of the data is obtained with more extreme axial ratios ($1:0.8:0.4$) than those derived by Jones et al. More surprisingly, we find that starless cores have more extreme axial ratios than protostellar cores -- starless cores appear more `flattened'. This is the opposite of what would be expected from modeling the freefall collapse of triaxial ellipsoids. The collapse of starless cores would be expected to proceed most swiftly along the shortest axis - as has been predicted by every modeller since Zel'dovich - which should produce more flattened cores around protostars, the opposite of what is seen.

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.

ISO data taken with the long-wavelength imaging photo-polarimeter ISOPHOT are presented of 18 pre-stellar cores at three far-infrared wavelengths - 90, 170 and 200 microns. Most of the cores are detected clearly at 170 and 200 but only one is detected strongly at 90 microns, indicating that mostly they are very cold, with typical temperatures of only 10-20K. Colour temperature images are constructed for each of the cores. Most of the cores are seen to be either isothermal, or to have associated temperature gradients from the core centres to their edges, with all except one being cooler at the centre. We compare the data with previous ISOCAM absorption data and calculate the energy balance for those cores in common between the two samples. We find that the energy radiated by each core in the far-IR is similar to that absorbed at shorter wavelengths. Hence there is no evidence for a central heating source in any of the cores - even those for which previous evidence for core contraction exists. This is all consistent with external heating of the cores by the local interstellar radiation field, confirming their pre-stellar nature.

We present hydrodynamical simulations of the formation, structure and evolution of photoionized columns, with parameters based on those observed in the Eagle Nebula. On the basis of these simulations we argue that there is no unequivocal evidence that the dense neutral clumps at heads of the columns were cores in the pre-existing molecular cloud. In our simulations, a variety of initial conditions leads to the formation and maintenance of near-equilibrium columns. Therefore, it is likely that narrow columns will often occur in regions with large-scale inhomogeneities, but that observations of such columns can tell us little about the processes by which they formed. The manner in which the columns in our simulations develop suggests that their evolution may result in extended sequences of radiation-induced star formation.

We present a radiative transfer model, which is applicable to the study of submillimetre spectral line observations of protostellar envelopes. The model uses an exact, non-LTE, spherically symmetric radiative transfer `Stenholm' method, which numerically solves the radiative transfer problem by the process of `Lambda-iteration'. We also present submillimetre spectral line data of the Class 0 protostars NGC1333-IRAS2 and Serpens SMM4. We examine the physical constraints which can be used to limit the number and range of parameters used in protostellar envelope models, and identify the turbulent velocity and tracer molecule abundance as the principle sources of uncertainty in the radiative transfer modelling. We explore the trends in the appearance of the predicted line profiles as key parameters in the models are varied. We find that the separation of the two peaks of a typical infall profile is dependent not on the evolutionary status of the collapsing protostar, but on the turbulent velocity dispersion in the envelope. We also find that the line shapes can be significantly altered by rotation. Fits are found for the observed line profiles of IRAS2 and SMM4 using plausible infall model parameters. The density and velocity profiles in our best fit models are inconsistent with a singular isothermal sphere model. We find better agreement with a form of collapse which assumes non-static initial conditions. We also find some evidence that the infall velocities are retarded from free-fall towards the centre of the cloud, probably by rotation, and that the envelope of SMM4 is rotationally flattened.

We describe the results of a sequence of simulations of gravitational collapse in a turbulent magnetized region. The parameters are chosen to be representative of molecular cloud material. We find that several protostellar cores and filamentary structures of higher than average density form. The filaments inter-connect the high density cores. Furthermore, the magnetic field strengths are found to correlate positively with the density, in agreement with recent observations. We make synthetic channel maps of the simulations and show that material accreting onto the cores is channelled along the magnetized filamentary structures. This is compared with recent observations of S106, and shown to be consistent with these data. We postulate that this mechanism of accretion along filaments may provide a means for molecular cloud cores to grow to the point where they become gravitationally unstable and collapse to form stars.

We report the results of a submillimeter dust continuum survey of the protoclusters NGC 2068 and NGC 2071 in Orion B carried out at 850 microns and 450 microns with SCUBA on JCMT. The mapped region is ~ 32' x 18' in size (~ 4 pc x 2 pc) and consists of filamentary dense cores which break up into small-scale (~ 5000 AU) fragments, including 70 starless condensations and 5 circumstellar envelopes/disks. The starless condensations, seen on the same spatial scales as protostellar envelopes, are likely to be gravitationally bound and pre-stellar in nature. Their mass spectrum, ranging from ~ 0.3 Msun to ~ 5 Msun, is reminiscent of the stellar initial mass function (IMF). Their mass-size relation suggests that they originate from gravitationally-driven fragmentation. We thus argue that pre-collapse cloud fragmentation plays a major role in shaping the IMF.

Observations of dark cloud cores have been carried out in the mid-infrared using ISOCAM and in the far-infrared using ISOPHOT, both aboard the Infrared Space Observatory. The cores are in most cases detected in emission at 200 and 170 micron, remain undetected at 90 micron and are seen in absorption against the diffuse mid-infrared background at 7 micron. The observations are consistent with the cores being pre-stellar and not having a central heating source, and yield core temperatures of ~ 11-13 K. We were able to determine the density structure of the cores up to radii that extend beyond the sensitivity limit of previous submillimetre studies. The column density profiles of the cores studied here flatten out in the centre, as shown by previous submillimetre continuum results, and interestingly, in a few cases, the profiles steepen with radius beyond ~ 5000-10000 AU and present sharp edges. This could indicate that these cores are decoupled from their parent molecular cloud and represent finite reservoirs of mass for subsequent accretion.

We present the results of short baseline interferometry observations at submillimetre wavelengths, using a two-element interferometer comprising the JCMT and CSO, of circumstellar discs around young YSOs. We model data for the Class 0 protostar IRAS03282 and the Class I protostar L1709B with discs and compare them with previously published results. We find evidence for the ratio of the disc mass to envelope mass to increase as the objects evolve.

We present C18O observations of the pre-stellar core L1689B, in the (J=3-2) and (J=2-1) rotational transitions, taken at the James Clerk Maxwell Telescope in Hawaii. We use a lambda-iteration radiative transfer code to model the data. We adopt a similar form of radial density profile to that which we have found in all pre-stellar cores, with a `flat' inner profile, steepening towards the edge, but we make the gradient of the `flat' region a free parameter. We find that the core is close to virial equilibrium, but there is tentative evidence for core contraction. We allow the temperature to vary with a power-law form and find we can consistently fit all of the CO data with an inverse temperature gradient that is warmer at the edge than the centre. However, when we combine the CO data with the previously published millimetre data we fail to find a simultaneous fit to both data-sets without additionally allowing the CO abundance to decrease towards the centre. This effect has been observed qualitatively many times before, as the CO freezes out onto the dust grains at high densities, but we quantify the effect. Hence we show that the combination of mm/submm continuum and spectral line data is a very powerful method of constraining the physical parameters of cores on the verge of forming stars.

Molecular outflows from young protostars are widely believed to be collimated by magnetic fields, but there has been little observational evidence to support this hypothesis. Using the new technique of millimetre-wavelength spectro-polarimetry, we demonstrate the existence of a magnetic field in the NGC2024-FIR5 outflow lobe. The 1.3mm J=2-1 transition of carbon monoxide (CO) is polarized at a level of approximately 1%, in a direction within 10-15 degrees of the outflow axis. This agrees with theoretical models where the magnetic field channels the outflowing gas, and shows that the process can be effective as far as 0.1pc from the protostar.

We present first results from a submillimetre continuum survey of the Galactic Centre `Central Molecular Zone' (CMZ), made with SCUBA on the James Clerk Maxwell Telescope. SCUBA's scan-map mode has allowed us to make extremely wide-field maps of thermal dust emission with unprecedented speed and sensitivity. We also discuss some issues related to the elimination of artefacts in scan-map data. Our simultaneous 850/450 micron maps have a total size of approximately 2.8 x 0.5 degrees (400 x 75 pc) elongated along the galactic plane. They cover the Sgr A region-including Sgr A*, the circumnuclear disc, and the +20 km/s and +50 km/s clouds; the area around the Pistol; Sgr B2-the brightest feature on the map; and at their Galactic Western and Eastern edges the Sgr C and Sgr D regions. There are many striking features such as filaments and shell-like structures, as well as point sources such as Sgr A* itself. The total mass in the Central Molecular Zone is greater than that revealed in previous optically-thin molecular line maps by a factor of ~3, and new details are revealed on scales down to 0.33 pc across this 400 pc wide region.

We propose a new analytic model for the initial conditions of protostellar collapse in relatively isolated regions of star formation. The model is non-magnetic, and is based on a Plummer-like radial density profile as its initial condition. It fits: the observed density profiles of pre-stellar cores and Class 0 protostars; recent observations in pre-stellar cores of roughly constant contraction velocities over a wide range of radii; and the lifetimes and accretion rates derived for Class 0 and Class I protostars. However, the model is very simple, having in effect only 2 free parameters, and so should provide a useful framework for interpreting observations of pre-stellar cores and protostars, and for calculations of radiation transport and time-dependent chemistry. As an example, we model the pre-stellar core L1544.

We present the results of a mid-infrared (7 micron) imaging survey of a sample of 24 starless dense cores carried out at an angular resolution of 6 arcsec with the ISOCAM camera aboard the Infrared Space Observatory (ISO). The targeted cores are believed to be pre-stellar in nature and to represent the initial conditions of low-mass, isolated star formation. In previous submillimeter dust continuum studies of such pre-stellar cores, it was found that the derived column density profiles did not follow a single power-law such as N[H2] \propto r^(-1) throughout their full extent but flattened out near their center. These submillimeter observations however could not constrain the density profiles at radii greater than ~ 10000 AU. The present absorption study uses ISOCAM's sensitivity to map these pre-stellar cores in absorption against the diffuse mid-infrared background. The goal was to determine their structure at radii that extend beyond the limits of sensitivity of the submillimeter continuum maps and at twiceas good an angular resolution. Among the 24 cores observed in our survey, a majority of them show deep absorption features. The starless cores studied here all show a column density profile that flattens in the center, which confirms the submillimeter emission results. Moreover, beyond a radius of ~ 5000-10000 AU, the typical column density profile steepens with distance from core center and gets steeper than N[H2] \propto r^(-1), until it eventually merges with the low-density ambient molecular cloud. At least three of the cores present sharp edges at R ~ 15000-30000 AU and appear to be decoupled from their parent clouds, providing finite reservoirs of mass for subsequent star formation.

We present the first published maps of magnetic fields in pre-stellar cores, to test theoretical ideas about the way in which the magnetic field geometry affects the star formation process. The observations are JCMT-SCUBA maps of 850 micron thermal emission from dust. Linear polarizations at typically ten or more independent positions in each of three objects, L1544, L183 and L43 were measured, and the geometries of the magnetic fields in the plane of the sky were mapped from the polarization directions. The observed polarizations in all three objects appear smooth and fairly uniform. In L1544 and L183 the mean magnetic fields are at an angle of around 30 degrees to the minor axes of the cores. The L43 B-field appears to have been influenced in its southern half, such that it is parallel to the wall of a cavity produced by a CO outflow from a nearby T Tauri star, whilst in the northern half the field appears less disturbed and has an angle of 44 degrees to the core minor axis. We briefly compare our results with published models of magnetized cloud cores and conclude that no current model can explain these observations simultaneously with previous ISOCAM data.

We present a catalogue of molecular cloud cores drawn from high latitude, medium opacity clouds, using the all-sky IRAS Sky Survey Atlas (ISSA) images at 60 and 100~$\mu$m. The typical column densities of the cores are $ N(H_2)\sim 3.8 \times 10^{21} $cm$^{-2}$ and the typical volume densities are $ n(H_2) \sim 2 \times 10^3$cm$^{-3}$. They are therefore significantly less dense than many other samples obtained in other ways. Those cloud cores with IRAS point sources are seen to be already forming stars, but this is found to be only a small fraction of the total number of cores. The fraction of the cores in the protostellar stage is used to estimate the prestellar timescale - the time until the formation of a hydrostatically supported protostellar object. We argue, on the basis of a comparison with other samples, that a trend exists for the prestellar lifetime of a cloud core to decrease with the mean column density and number density of the core. We compare this with model predictions and show that the data are consistent with star formation regulated by the ionisation fraction.

The last decade has witnessed significant advances in our observational understanding of the earliest stages of low-mass star formation. The advent of sensitive receivers on large radio telescopes such as the JCMT and IRAM 30m MRT has led to the identification of young protostars at the beginning of the main accretion phase (`Class 0' objects), and has made it possible to probe, for the first time, the inner density structure of pre-collapse cores. Class 0 objects are characterized by strong, centrally-condensed dust continuum emission at submillimeter wavelengths, very little emission shortward of 10 microns, and powerful jet-like outflows. Direct evidence for gravitational infall has been observed toward several of them. They are interpreted as accreting protostars which have not yet accumulated the majority of their final stellar mass. In contrast to protostars, pre-stellar cores have flat inner density profiles, suggesting the initial conditions for fast protostellar collapse depart sometimes significantly from a singular isothermal sphere. In the case of non-singular initial conditions, the beginning of protostellar evolution is expected to feature a brief phase of vigorous accretion/ejection which may coincide with Class 0 objects. In addition, submillimeter continuum imaging surveys of regions of multiple star formation such as Ophiuchus and Serpens suggest a picture according to which each star in an embedded cluster is built from a finite reservoir of mass and the associated IMF is primarily determined at the pre-stellar stage of evolution.

Continuum emission at 220 and 355 GHz from the prototype Class 0 source VLA 1623 has been detected using the JCMT-CSO interferometer. Gaussian fits to the data place an upper limit of 70 AU on the half-width half-maximum radius of the emission, which implies an upper limit of ~175 AU for the cutoff radius of the circumstellar disk in the system. In the context of existing collapse models, this disk could be magnetically supported on the largest scales and have an age of ~6x10^4 yr, consistent with previous suggestions that Class 0 sources are quite young. The innermost region of the disk within ~6 AU is likely to be in centrifugal support, which is likely large enough to provide a drive for the outflow according to current theoretical models. Alternatively, if 175 AU corresponds to the centrifugal radius of the disk, the age of the system is ~2x10^5 yr, closer to age estimates for Class I sources.

We present submm images of Sgr A* and its surroundings obtained at 800, 600 and 450 $\mu$m with the JCMT and derive flux densities of Sgr A* at all three wave- lengths. Combined with upper limits by Gezari and associates at MIR wavelengths a time averaged radio spectrum is obtained which increases $\propto \nu^{1/3}$, attains a maximum at $\nu_max \ga 600 GHz$ and must decrease rapidly at fre- quencies $\ga 10^4$GHz. This spectrum allows for about $3 M}_\odot$ of 50K dust and associated hydrogen in the telescope beam. While variations on time-scales of a few months are now well established for the frequency range $\nu \la 100 GHz$, our investigation of variability at higher frequencies still yields only marginal results. The Circum-Nuclear Disk (CND) extends over the central 12 pc. At the galactocentric radius $R \sim 1 pc$ dust and hydrogen column densities drop to low values and form the Central Cavity. Our submm images show that the bottom of this cavity is rather flat. Variations in the dust emission are just consistent with the detection of a `Tongue' of $\sim 200 M_\odot$ of atomic hydrogen reported by Jackson et al. (1993) to be located between the Northern and the Eastern Arm of the Minispiral. We present a revised submm/IR spectrum of the central $30\arcsec$ (R$\la 0.6 pc$) with flux densities corrected for an interstellar extinction of $A_v \sim 31 mag$. This spectrum attains its maximum at $\lambda \sim 20 \mu m$ and comes from dust with temperatures $\sim 170 - 400 K$ which is associated with the Eastern Arm and the East-West Bar. The in- tegrated luminosity is $\sim 5 10^6 L_\odot$ to which emission at $\lambda \la 30 \mu m$ contributes $\sim 80\%$. Heating of this dust is not provided by a central source but rather by a cluster of hot ($T eff \sim 3 - 3.5 10^4 K$) and luminous stars which could include the He{\sc i}/H{\sc i}-stars detected by