The formation and evolution of cold diffuse clouds (CDCs), the parent objects of dense molecular clouds, affects both the star formation process and that of larger-scale galactic evolution. We have begun a pilot study of one CDC's dust content, with the goal of quantifying the abundances of different types of dust and relating these to the relative abundance of molecular gas, the cloud's physical properties, and its general stage of development. Using photometry from AKARI and other surveys, we have extracted a sample spectral energy distribution (SED) of the CDC dust thermal emission over the near-infrared to submillimeter range. The extracted SED closely resembles others in the literature, confirming our isolation of the cloud emission from other sources along the sight line. We plan to fit this SED with dust models at each position in the cloud, automating our procedure to map out the structure of this CDC and others.

The interstellar medium is a complex 'ecosystem' with gas constituents in the atomic, molecular, and ionized states, dust, magnetic fields, and relativistic particles. The Canadian Galactic Plane Survey has imaged these constituents with angular resolution of the order of arcminutes. This paper presents radio continuum data at 408 MHz over the area 52 degrees < longitude < 193 degrees, -6.5 degrees < latitude < 8.5 degrees, with an extension to latitude = 21 degrees in the range 97 degrees < longitude < 120 degrees, with angular resolution 2.8' x 2.8' cosec(declination). Observations were made with the Synthesis Telescope at the Dominion Radio Astrophysical Observatory as part of the Canadian Galactic Plane Survey. The calibration of the survey using existing radio source catalogs is described. The accuracy of 408-MHz flux densities from the data is 6%. Information on large structures has been incorporated into the data using the single-antenna survey of Haslam (1982). The paper presents the data, describes how it can be accessed electronically, and gives examples of applications of the data to ISM research.

We present the first data release of the James Clerk Maxwell Telescope (JCMT) Plane Survey (JPS), the JPS Public Release 1 (JPSPR1). JPS is an 850-um continuum survey of six fields in the northern inner Galactic Plane in a longitude range of l=7-63, made with the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2). This first data release consists of emission maps of the six JPS regions with an average pixel-to-pixel noise of 7.19 mJy beam^-1, when smoothed over the beam, and a compact-source catalogue containing 7,813 sources. The 95 per cent completeness limits of the catalogue are estimated at 0.04 Jy beam^-1 and 0.3 Jy for the peak and integrated flux densities, respectively. The emission contained in the compact-source catalogue is 42 +- 5 per cent of the total and, apart from the large-scale (greater than 8') emission, there is excellent correspondence with features in the 500-um Herschel maps. We find that, with two-dimensional matching, 98 +- 2 per cent of sources within the fields centred at l=20, 30, 40 and 50 are associated with molecular clouds, with 91 +- 3 per cent of the l=30 and 40 sources associated with dense molecular clumps. Matching the JPS catalogue to Herschel 70-um sources, we find that 38 +- 1 per cent of sources show evidence of ongoing star formation. The images and catalogue will be a valuable resource for studies of star formation in the Galaxy and the role of environment and spiral arms in the star formation process.

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

Emission from carbon monoxide (CO) is ubiquitously used as a tracer of dense star forming molecular clouds. There is, however, growing evidence that a significant fraction of CO emission originates from diffuse molecular gas. Quantifying the contribution of diffuse CO-emitting gas is vital for understanding the relation between molecular gas and star formation. We examine the Galactic distribution of two CO-emitting gas components, a high column density component detected in 13CO and 12CO, and a low column density component detected in 12CO, but not in 13CO. The "diffuse" and "dense" components are identified using a combination of smoothing, masking, and erosion/dilation procedures, making use of three large-scale 12CO and 13CO surveys of the Inner and Outer Milky Way. The diffuse component, which globally represents 25 (1.5x1e8 Mo) of the total molecular gas mass (6.5x1e8 Mo), is more extended perpendicular to the Galactic plane. The fraction of diffuse gas increases from 15% at a galactocentric radius of 3 kpc to 50% at 15 kpc, and increases with decreasing surface density. In the Inner Galaxy, a yet denser component traced by CS emission represents 14% of the total molecular gas mass traced by 12CO emission. Only 14% of the molecular gas mass traced by 12CO emission is identified as part of molecular clouds in 12CO surveys by cloud identification algorithms. This study indicates that CO emission not only traces star forming clouds, but also a significant diffuse molecular ISM component.

We present the $^{13}$CO/C$^{18}$O (J=3-2) Heterodyne Inner Milky Way Plane Survey (CHIMPS) which has been carried out using the Heterodyne Array Receiver Program on the 15 m James Clerk Maxwell Telescope (JCMT) in Hawaii. The high-resolution spectral survey currently covers |b| < 0.5 deg and 28 < l < 46 deg, with an angular resolution of 15 arcsec in 0.5 km/s velocity channels. The spectra have a median rms of $\sim$ 0.6 K at this resolution, and for optically thin gas at an excitation temperature of 10 K, this sensitivity corresponds to column densities of $N_{\mathrm{H}_{2}} \sim 3 \times 10^{20}\,$cm$^{-2}$ and $N_{\mathrm{H}_{2}} \sim 4 \times 10^{21}\,$cm$^{-2}$ for $^{13}$CO and C$^{18}$O, respectively. The molecular gas that CHIMPS traces is at higher column densities and is also more optically thin than in other publicly available CO surveys due to its rarer isotopologues, and thus more representative of the three-dimensional structure of the clouds. The critical density of the J=3-2 transition of CO is $\gtrsim 10^{4}$ cm$^{-3}$ at temperatures of $\leq 20$ K, and so the higher density gas associated with star formation is well traced. These data complement other existing Galactic plane surveys, especially the JCMT Galactic Plane Survey which has similar spatial resolution and column density sensitivity, and the Herschel infrared Galactic Plane Survey. In this paper, we discuss the observations, data reduction and characteristics of the survey, presenting integrated emission maps for the region covered. Position-velocity diagrams allow comparison with Galactic structure models of the Milky Way, and while we find good agreement with a particular four arm model, there are some significant deviations.

We have mapped cold atomic gas in 21cm line HI self-absorption (HISA) at arcminute resolution over more than 90% of the Milky Way's disk. To probe the formation of H2 clouds, we have compared our HISA distribution with CO J=1-0 line emission. Few HISA features in the outer Galaxy have CO at the same position and velocity, while most inner-Galaxy HISA has overlapping CO. But many apparent inner-Galaxy HISA-CO associations can be explained as chance superpositions, so most inner-Galaxy HISA may also be CO-free. Since standard equilibrium cloud models cannot explain the very cold HI in many HISA features without molecules being present, these clouds may instead have significant CO-dark H2.

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-{\mu}m continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-{\sigma} surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average.

Power-law tails are often seen in probability density functions (PDFs) of molecular cloud column densities, and have been attributed to the effect of gravity. We show that extinction PDFs of a sample of five molecular clouds obtained at a few tenths of a parsec resolution, probing extinctions up to A$_{{\mathrm{V}}}$ $\sim$ 10 magnitudes, are very well described by lognormal functions provided that the field selection is tightly constrained to the cold, molecular zone and that noise and foreground contamination are appropriately accounted for. In general, field selections that incorporate warm, diffuse material in addition to the cold, molecular material will display apparent core+tail PDFs. The apparent tail, however, is best understood as the high extinction part of a lognormal PDF arising from the cold, molecular part of the cloud. We also describe the effects of noise and foreground/background contamination on the PDF structure, and show that these can, if not appropriately accounted for, induce spurious tails or amplify any that are truly present.

We present CO, H2, HI and HISA distributions from a set of simulations of grand design spirals including stellar feedback, self-gravity, heating and cooling. We replicate the emission of the 2nd Galactic Quadrant by placing the observer inside the modelled galaxies and post process the simulations using a radiative transfer code, so as to create synthetic observations. We compare the synthetic datacubes to observations of the 2nd Quadrant of the Milky Way to test the ability of the current models to reproduce the basic chemistry of the Galactic ISM, as well as to test how sensitive such galaxy models are to different recipes of chemistry and/or feedback. We find that models which include feedback and self-gravity can reproduce the production of CO with respect to H2 as observed in our Galaxy, as well as the distribution of the material perpendicular to the Galactic plane. While changes in the chemistry/feedback recipes do not have a huge impact on the statistical properties of the chemistry in the simulated galaxies, we find that the inclusion of both feedback and self-gravity are crucial ingredients, as our test without feedback failed to reproduce all of the observables. Finally, even though the transition from H2 to CO seems to be robust, we find that all models seem to underproduce molecular gas, and have a lower molecular to atomic gas fraction than is observed. Nevertheless, our fiducial model with feedback and self-gravity has shown to be robust in reproducing the statistical properties of the basic molecular gas components of the ISM in our Galaxy.

We present a new catalogue of spectrophotometric distances and line-of-sight systemic velocities to 103 HII regions between 90$^{\circ}\leq \ell \leq~$195$^{\circ}$ (longitude quadrants II and part of III). Two new velocities for each region are independently measured using 1-arcminute resolution 21~cm HI and 2.6~mm $^{12}$CO line maps (from the Canadian Galactic Plane Survey and FCRAO Outer Galaxy Surveys) that show where gaseous shells are observed around the periphery of the ionized gas. Known and neighbouring OB-type stars with published UBV photometry and MK classifications are overlaid onto 21~cm continuum maps, and those stars observed within the boundary of the HII emission (and whose distance is not more than 3 times the standard deviation of the others) are used to calculate new mean stellar distances to each of the 103 nebulae. Using this approach of excluding distance outliers from the mean distance to a group of many stars in each HII region lessens the impact of anomalous reddening for certain individuals. Final mean distances of 9 common objects with VLBI parallax distances show a 1:1 correspondence. Further, comparison with previous catalogues of HII regions in these quadrants shows a 50% reduction in scatter for the distance to Perseus spiral arm objects in the same region, and a reduction by $\sim$1/$\sqrt{2}$ in scatter around a common angular velocity relative to the Sun $\Omega-\Omega_0$(km~s$^{-1}$~kpc$^{-1}$). The purpose of the catalogue is to provide a foundation for more detailed large-scale Galactic spiral structure and dynamics (rotation curve, density wave streaming) studies in the 2$^{\textrm{nd}}$ and 3$^{\textrm{rd}}$ quadrants, which from the Sun's location is the most favourably viewed section of the Galaxy.

We introduce a new method for observationally estimating the fraction of momentum density (${\rho}{\mathbf{v}}$) power contained in solenoidal modes (for which $\nabla \cdot {\rho}{\mathbf{v}} = 0$) in molecular clouds. The method is successfully tested with numerical simulations of supersonic turbulence that produce the full range of possible solenoidal/compressible fractions. At present the method assumes statistical isotropy, and does not account for anisotropies caused by (e.g.) magnetic fields. We also introduce a framework for statistically describing density--velocity correlations in turbulent clouds.

We present Spitzer, NIR and millimeter observations of the massive star forming regions W5-east, S235, S252, S254-S258 and NGC7538. Spitzer data is combined with near-IR observations to identify and classify the young population while 12CO and 13CO observations are used to examine the parental molecular cloud. We detect in total 3021 young stellar objects (YSOs). Of those, 539 are classified as Class I, and 1186 as Class II sources. YSOs are distributed in groups surrounded by a more scattered population. Class I sources are more hierarchically organized than Class II and associated with the most dense molecular material. We identify in total 41 embedded clusters containing between 52 and 73% of the YSOs. Clusters are in general non-virialized, turbulent and have star formation efficiencies between 5 and 50%. We compare the physical properties of embedded clusters harboring massive stars (MEC) and low-mass embedded clusters (LEC) and find that both groups follow similar correlations where the MEC are an extrapolation of the LEC. The mean separation between MEC members is smaller compared to the cluster Jeans length than for LEC members. These results are in agreement with a scenario where stars are formed in hierarchically distributed dusty filaments where fragmentation is mainly driven by turbulence for the more massive clusters. We find several young OB-type stars having IR-excess emission which may be due to the presence of an accretion disk.

Principal component analysis is a powerful statistical system to investigate the structure and dynamics of the molecular interstellar medium, with particular emphasis on the study of turbulence, as revealed by spectroscopic imaging of molecular line emission. To-date, the method to retrieve the power law index of the velocity structure function or power spectrum has relied on an empirical calibration and testing with model turbulent velocity fields, while lacking a firm theoretical basis. In this paper, we present an analytic formulation that reveals the detailed mechanics of the method and confirms previous empirical calibrations of its recovery of the scale dependence of turbulent velocity fluctuations.

Infrared stellar photometry from 2MASS and spectral line imaging observations of 12CO and 13CO J = 1-0 line emission from the FCRAO 14m telescope are analysed to assess the variation of the CO abundance with physical conditions throughout the Orion A and Orion B molecular clouds. Three distinct Av regimes are identified in which the ratio between the 13CO column density and visual extinction changes corresponding to the photon dominated envelope, the strongly self-shielded interior, and the cold, dense volumes of the clouds. Within the strongly self-shielded interior of the Orion A cloud, the 13CO abundance varies by 100% with a peak value located near regions of enhanced star formation activity. The effect of CO depletion onto the ice mantles of dust grains is limited to regions with AV > 10 mag and gas temperatures less than 20 K as predicted by chemical models that consider thermal-evaporation to desorb molecules from grain surfaces. Values of the molecular mass of each cloud are independently derived from the distributions of Av and 13CO column densities with a constant 13CO-to-H2 abundance over various extinction ranges. Within the strongly self-shielded interior of the cloud (Av > 3 mag), 13CO provides a reliable tracer of H2 mass with the exception of the cold, dense volumes where depletion is important. However, owing to its reduced abundance, 13CO does not trace the H2 mass that resides in the extended cloud envelope, which comprises 40-50% of the molecular mass of each cloud. The implied CO luminosity to mass ratios, M/L_{CO}, are 3.2 and 2.9 for Orion A and Orion B respectively, which are comparable to the value (2.9), derived from gamma-ray observations of the Orion region. Our results emphasize the need to consider local conditions when applying CO observations to derive H2 column densities.

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

We generate synthetic HI Galactic plane surveys from spiral galaxy simulations which include stellar feedback processes. Compared to a model without feedback we find an increased scale height of HI emission (in better agreement with observations) and more realistic spatial structure (including supernova blown bubbles). The synthetic data show HI self-absorption with a morphology similar to that seen in observations. The density and temperature of the material responsible for HI self-absorption is consistent with observationally determined values, and is found to be only weakly dependent on absorption strength and star formation efficiency.

We present $^{12}$CO, $^{13}$CO and C$^{18}$O (J=3$-$2) observations of a new cluster of outflows in the Vulpecula Rift with HARP-B on the JCMT. The mass associated with the outflows, measured using the $^{12}$CO HARP-B observations and assuming a distance to the region of 2.3 kpc, is 129 \msol{}, while the mass associated with the dense gas from C$^{18}$O observations is 458 \msol{} and the associated sub-millimeter core has a mass of 327 $\pm$ 112 \msol{} independently determined from Bolocam 1.1mm data. The outflow-to-core mass ratio is therefore $\sim$0.4, making this region one of the most efficient observed thus far with more than an order of magnitude more mass in the outflow than would be expected based on previous results. The kinetic energy associated with the flows, 94$\times10^{45}$ ergs, is enough to drive the turbulence in the local clump, and potentially unbind the local region altogether. The detection of SiO (J=8$-$7) emission toward the outflows indicates that the flow is still active, and not simply a fossil flow. We also model the SEDs of the four YSOs associated with the molecular material, finding them all to be of mid to early B spectral type. The energetic nature of the outflows and significant reservoir of cold dust detected in the sub-mm suggest that these intermediate mass YSOs will continue to accrete and become massive, rather than reach the main sequence at their current mass.

Magnetically aligned velocity anisotropy over varying physical conditions and environments within the Taurus Molecular Cloud is evaluated from analysis of wide field spectroscopic imaging of 12CO and 13CO J=1-0 emission. Such anisotropy is a result of MHD turbulence in the strong magnetic field regime and provides an indirect measure of the role of magnetic fields upon the gas. Velocity anisotropy aligned with the local, projected mean magnetic field direction is limited to fields with low surface brightness 12CO emission corresponding to regions of low visual extinction and presumably, low gas volume density. The more optically thin 13CO J=1-0 emission shows little evidence for velocity anisotropy. We compare our results with computational simulations with varying degrees of magnetic field strength and Alfvenic Mach numbers. In the diffuse, molecular envelope of the cloud, a strong magnetic field and sub-Alfvenic turbulent motions are inferred. Super-Alfvenic motions are present within the high column density filaments of the Taurus cloud. From this trans-Alfvenic flow, we constrain the scaling exponent, kappa, of the magnetic field density relation (B ~ n^kappa) to be near zero as expected for ambipolar diffusion or material loading of magnetic flux tubes.

Turbulence plays a major role in the formation and evolution of molecular clouds. The problem is that turbulent velocities are convolved with the density of an observed region. To correct for this convolution, we investigate the relation between the turbulence spectrum of model clouds, and the statistics of their synthetic observations obtained from Principal Component Analysis (PCA). We apply PCA to spectral maps generated from simulated density and velocity fields, obtained from hydrodynamic simulations of supersonic turbulence, and from fractional Brownian motion fields with varying velocity, density spectra, and density dispersion. We examine the dependence of the slope of the PCA structure function, alpha_PCA, on intermittency, on the turbulence velocity (beta_v) and density (beta_n) spectral indexes, and on density dispersion. We find that PCA is insensitive to beta_n and to the log-density dispersion sigma_s, provided sigma_s < 2. For sigma_s>2, alpha_PCA increases with sigma_s due to the intermittent sampling of the velocity field by the density field. The PCA calibration also depends on intermittency. We derive a PCA calibration based on fBms with sigma_s<2 and apply it to 367 CO spectral maps of molecular clouds in the Galactic Ring Survey. The average slope of the PCA structure function, <alpha_PCA>=0.62\pm0.2, is consistent with the hydrodynamic simulations and leads to a turbulence velocity exponent <beta_v>=2.06\pm0.6 for a non-intermittent, low density dispersion flow. Accounting for intermittency and density dispersion, the coincidence between the PCA slope of the GRS clouds and the hydrodynamic simulations suggests beta_v~1.9, consistent with both Burgers and compressible intermittent turbulence.

Context: Despite their profound effect on the universe, the formation of massive stars and stellar clusters remains elusive. In the past decade, compelling evidence has emerged that suggests Infrared Dark Clouds (IRDCs) may be precursors to stellar clusters. However, the usual method for identifying IRDCs is biased by the requirement that they are seen in absorption against bright mid-IR emission, whereas dust continuum observations allow cold, dense pre-stellar-clusters to be identified anywhere. Aims: We aim to understand what physical properties characterize IRDCs, to explore the population of dust continuum sources that are not IRDCs, and to roughly characterize the star formation activity in dust continuum sources. Results: We present temperature and column density maps in the Hi-GAL l=30 and l=59 SDP fields, as well as a robust algorithm for cirrus subtraction and source identification using Hi-GAL. We report on the fraction of Hi-GAL sources which are mid-IR-dark, -neutral, or -bright in both fields. We find significant trends in column density and temperature between mid-IR-dark and -bright pixels; mid-IR-dark pixels are about 10 K colder and have a factor of 2 higher column density on average than mid-IR-bright pixels. We find that Hi-GAL dust continuum sources span a range of evolutionary states from pre- to star-forming, and that warmer sources are associated with more star formation tracers. There is a trend of increasing temperature with tracer type from mid-IR-dark at the coldest, to outflow/maser sources in the middle, and finally to 8 and 24 micron bright sources at the warmest. Finally, we identify five candidate IRDC-like sources on the far-side of the Galaxy. These are cold (~ 20 K), high column density (N(H2) > 10^22 cm^-2) clouds identified with Hi-GAL which, despite bright surrounding mid-IR emission, show little to no absorption at 8 micron. (abridged)

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.

We aim to better understand how the spatial structure of molecular clouds is governed by turbulence. For that, we study the large-scale spatial distribution of low density molecular gas and search for characteristic length scales. We employ a 35 square degrees 13CO 1-0 molecular line survey of Cygnus X and visual extinction (A_V) maps of 17 Galactic clouds to analyse the spatial structure using the Delta-variance method. This sample contains a large variety of different molecular cloud types with different star forming activity. The Delta-variance spectra obtained from the A_V maps show differences between low-mass star-forming (SF) clouds and massive giant molecular clouds (GMC) in terms of shape of the spectrum and its power-law exponent beta. Low-mass SF clouds have a double-peak structure with characteristic size scales around 1 pc (though with a large scatter around this value) and 4 pc. GMCs show no characteristic scale in the A_V-maps, which can partly be ascribed to a distance effect due to a larger line-of-sight (LOS) confusion. The Delta-variance for Cygnus, determined from the 13CO survey, shows characteristic scales at 4 pc and 40 pc, either reflecting the filament structure and large-scale turbulence forcing or - for the 4 pc scale - the scale below which the 13CO 1-0 line becomes optically thick. Though there are different processes that can introduce characteristic scales, i.e. geometry, decaying turbulence the transition scale from supersonic to subsonic turbulence (the sonic scale), line-of-sight effects and energy injection due to expanding supernova shells, outflows, HII-regions, and the relative contribution of these effects strongly varies from cloud to cloud, it is remarkable that the resulting turbulent structure of molecular clouds shows similar characteristics.

We examine the relation between the density variance and the mean-square Mach number in supersonic, isothermal turbulence, assumed in several recent analytic models of the star formation process. From a series of calculations of supersonic, hydrodynamic turbulence driven using purely solenoidal Fourier modes, we find that the `standard' relationship between the variance in the log of density and the Mach number squared, i.e., sigma^2_(ln rho/rhobar)=ln (1+b^2 M^2), with b = 1/3 is a good fit to the numerical results in the supersonic regime up to at least Mach 20, similar to previous determinations at lower Mach numbers. While direct measurements of the variance in linear density are found to be severely underestimated by finite resolution effects, it is possible to infer the linear density variance via the assumption of log-normality in the Probability Distribution Function. The inferred relationship with Mach number, consistent with sigma_(rho/rhobar) ~ b M with b=1/3, is, however, significantly shallower than observational determinations of the relationship in the Taurus Molecular Cloud and IC5146 (both consistent with b~ 0.5), implying that additional physics such as gravity is important in these clouds and/or that turbulent driving in the ISM contains a significant compressive component. Magnetic fields are not found to change this picture significantly, in general reducing the measured variances and thus worsening the discrepancy with observations.

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

We present an implementation of the iterative flux-conserving Lucy-Richardson (L-R) deconvolution method of image restoration for maps produced by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). We have analyzed its performance and convergence extensively through simulations and cross-correlations of the deconvolved images with available highresolution maps. We present new science results from two BLAST surveys, in the Galactic regions K3-50 and IC 5146, further demonstrating the benefits of performing this deconvolution. We have resolved three clumps within a radius of 4.'5 inside the star-forming molecular cloud containing K3-50. Combining the well-resolved dust emission map with available multi-wavelength data, we have constrained the Spectral Energy Distributions (SEDs) of five clumps to obtain masses (M), bolometric luminosities (L), and dust temperatures (T). The L-M diagram has been used as a diagnostic tool to estimate the evolutionary stages of the clumps. There are close relationships between dust continuum emission and both 21-cm radio continuum and 12CO molecular line emission. The restored extended large scale structures in the Northern Streamer of IC 5146 have a strong spatial correlation with both SCUBA and high resolution extinction images. A dust temperature of 12 K has been obtained for the central filament. We report physical properties of ten compact sources, including six associated protostars, by fitting SEDs to multi-wavelength data. All of these compact sources are still quite cold (typical temperature below ~ 16 K) and are above the critical Bonner-Ebert mass. They have associated low-power Young Stellar Objects (YSOs). Further evidence for starless clumps has also been found in the IC 5146 region.

(abridged) We report a study of the relation between dust and gas over a 100deg^2 area in the Taurus molecular cloud. We compare the H2 column density derived from dust extinction with the CO column density derived from the 12CO and 13CO J= 1-0 lines. We derive the visual extinction from reddening determined from 2MASS data. The comparison is done at an angular size of 200", corresponding to 0.14pc at a distance of 140pc. We find that the relation between visual extinction Av and N(CO) is linear between Av~3 and 10 mag in the region associated with the B213--L1495 filament. In other regions the linear relation is flattened for Av > 4 mag. We find that the presence of temperature gradients in the molecular gas affects the determination of N(CO) by ~30--70% with the largest difference occurring at large column densities. Adding a correction for this effect and accounting for the observed relation between the column density of CO and CO2 ices and Av, we find a linear relationship between the column of carbon monoxide and dust for observed visual extinctions up to the maximum value in our data 23mag. We have used these data to study a sample of dense cores in Taurus. Fitting an analytical column density profile to these cores we derive an average volume density of about 1.4e4 cm^-3 and a CO depletion age of about 4.2e5 years. We estimate the H2 mass of Taurus to be about 1.5e4 M_sun, independently derived from the Av and N(CO) maps. We derive a CO integrated intensity to H2 conversion factor of about 2.1e20 cm^-2 (K km/s)^-1, which applies even in the region where the [CO]/[H_2] ratio is reduced by up to two orders of magnitude. The distribution of column densities in our Taurus maps resembles a log--normal function but shows tails at large and low column densities.

The Exeter FCRAO CO Galactic Plane Survey consists of 12CO and 13CO (J=1-0) observations over the galactic plane covering 55 degrees <= l <= 102 degrees, |b| >= 1 degree and 141 degrees <= l <= 195 degrees, -3.5 degrees <= b <= 5.5 degrees with a spatial resolution of ~45" and a spectral resolution of ~0.15km/s. We will present the methodology of a threshold-based cloud and clump determination method which retains hierarchical information, then discuss associating sources with clouds in the catalogue. Once complete, this catalogue of clouds and clumps will encompass the majority of the Northern Galactic Plane, providing knowledge of the molecular structure of the galaxy and the starting point for studies of the variation in star formation efficiency. In addition, it will allow us to identify clouds that have no or little star formation taking place inside them, which are often overlooked in the study of the conditions required for star formation to take place.

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

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

In Herschel images of the Galactic plane and many star forming regions, a major factor limiting our ability to extract faint compact sources is cirrus confusion noise, operationally defined as the "statistical error to be expected in photometric measurements due to confusion in a background of fluctuating surface brightness". The histogram of the flux densities of extracted sources shows a distinctive faint-end cutoff below which the catalog suffers from incompleteness and the flux densities become unreliable. This empirical cutoff should be closely related to the estimated cirrus noise and we show that this is the case. We compute the cirrus noise directly, both on Herschel images from which the bright sources have been removed and on simulated images of cirrus with statistically similar fluctuations. We connect these direct estimates with those from power spectrum analysis, which has been used extensively to predict the cirrus noise and provides insight into how it depends on various statistical properties and photometric operational parameters. We report multi-wavelength power spectra of diffuse Galactic dust emission from Hi-GAL observations at 70 to 500 microns within Galactic plane fields at l= 30 degrees and l= 59 degrees. We find that the exponent of the power spectrum is about -3. At 250 microns, the amplitude of the power spectrum increases roughly as the square of the median brightness of the map and so the expected cirrus noise scales linearly with the median brightness. Generally, the confusion noise will be a worse problem at longer wavelengths, because of the combination of lower angular resolution and the rising power spectrum of cirrus toward lower spatial frequencies, but the photometric signal to noise will also depend on the relative spectral energy distribution of the source compared to the cirrus.

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 have created synthetic neutral hydrogen (HI) Galactic Plane Survey data cubes covering 90 degrees < l < 180 degrees, using a model spiral galaxy from SPH simulations and the radiative transfer code TORUS. The density, temperature and other physical parameters are fed from the SPH simulation into TORUS, where the HI emissivity and opacity are calculated before the 21-cm line emission profile is determined. Our main focus is the observation of Outer Galaxy `Perseus Arm' HI, with a view to tracing atomic gas as it encounters shock motions as it enters a spiral arm interface, an early step in the formation of molecular clouds. The observation of HI self-absorption features at these shock sites (in both real observations and our synthetic data) allows us to investigate further the connection between cold atomic gas and the onset of molecular cloud formation.

Using the Torus radiative transfer code we produce synthetic observations of the 21 cm neutral hydrogen line from an SPH simulation of a spiral galaxy. The SPH representation of the galaxy is mapped onto an AMR grid, and a ray tracing method is used to calculate 21 cm line emission for lines of sight through the AMR grid in different velocity channels and spatial pixels. The result is a synthetic spectral cube which can be directly compared to real observations. We compare our synthetic spectral cubes to observations of M31 and M33 and find good agreement, whereby increasing velocity channels trace the main disc of the galaxy. The synthetic data also show kinks in the velocity across the spiral arms, evidence of non-circular velocities. These are still present even when we blur our data to a similar resolution as the observations, but largely absent in M31 and M33, indicating those galaxies do not contain significant spiral shocks. Thus the detailed velocity structure of our maps better represent previous observations of the grand design spiral M81.

We introduce a method for calculating the probability density function (PDF) of a turbulent density field in three dimensions using only information contained in the projected two-dimensional column density field. We test the method by applying it to numerical simulations of hydrodynamic and magnetohydrodynamic turbulence in molecular clouds. To a good approximation, the PDF of log(normalised column density) is a compressed, shifted version of the PDF of log(normalised density). The degree of compression can be determined observationally from the column density power spectrum, under the assumption of statistical isotropy of the turbulence.

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.

Supersonic turbulence in molecular clouds is a key agent in generating density enhancements that may subsequently go on to form stars. The stronger the turbulence - the higher the Mach number - the more extreme the density fluctuations are expected to be. Numerical models predict an increase in density variance with rms Mach number of the form: sigma^{2}_{rho/rho_{0}} = b^{2}M^{2}, where b is a numerically-estimated parameter, and this prediction forms the basis of a large number of analytic models of star formation. We provide an estimate of the parameter b from 13CO J=1-0 spectral line imaging observations and extinction mapping of the Taurus molecular cloud, using a recently developed technique that needs information contained solely in the projected column density field to calculate sigma^{2}_{rho/rho_{0}}. We find b ~ 0.48, which is consistent with typical numerical estimates, and is characteristic of turbulent driving that includes a mixture of solenoidal and compressive modes. More conservatively, we constrain b to lie in the range 0.3-0.8, depending on the influence of sub-resolution structure and the role of diffuse atomic material in the column density budget. We also report a break in the Taurus column density power spectrum at a scale of ~1pc, and find that the break is associated with anisotropy in the power spectrum. The break is observed in both 13CO and dust extinction power spectra, which, remarkably, are effectively identical despite detailed spatial differences between the 13CO and dust extinction maps. [ abridged ]

We introduce and test an expression for calculating the variance of a physical field in three dimensions using only information contained in the two-dimensional projection of the field. The method is general but assumes statistical isotropy. To test the method we apply it to numerical simulations of hydrodynamic and magnetohydrodynamic turbulence in molecular clouds, and demonstrate that it can recover the 3D normalised density variance with ~10% accuracy if the assumption of isotropy is valid. We show that the assumption of isotropy breaks down at low sonic Mach number if the turbulence is sub-Alfvenic. Theoretical predictions suggest that the 3D density variance should increase proportionally to the square of the Mach number of the turbulence. Application of our method will allow this prediction to be tested observationally and therefore constrain a large body of analytic models of star formation that rely on it.

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.

Supersonic turbulence in molecular clouds is a dominant agent that strongly affects the clouds' evolution and star formation activity. Turbulence may be initiated and maintained by a number of processes, acting at a wide range of physical scales. By examining the dynamical state of molecular clouds, it is possible to assess the primary candidates for how the turbulent energy is injected. The aim of this paper is to constrain the scales at which turbulence is driven in the molecular interstellar medium, by comparing simulated molecular spectral line observations of numerical magnetohydrodynamic (MHD) models and molecular spectral line observations of real molecular clouds. We use principal component analysis, applied to both models and observational data, to extract a quantitative measure of the driving scale of turbulence. We find that only models driven at large scales (comparable to, or exceeding, the size of the cloud) are consistent with observations. This result applies also to clouds with little or no internal star formation activity. Astrophysical processes acting on large scales, including supernova-driven turbulence, magnetorotational instability, or spiral shock forcing, are viable candidates for the generation and maintenance of molecular cloud turbulence. Small scale driving by sources internal to molecular clouds, such as outflows, can be important on small scales, but cannot replicate the observed large-scale velocity fluctuations in the molecular interstellar medium.

Using the radiative transfer code Torus, we produce spectral-line cubes of the predicted HI profile from global SPH simulations of spiral galaxies. Torus grids the SPH galaxy using Adaptive Mesh Refinement, then applies a ray-tracing method to infer the HI profile along the line(s) of sight. The gridded galaxy can be observed from any direction, which enables us to model the observed HI profile for galaxies of any orientation. We can also place the observer inside the galaxy, to simulate HI observations taken from the Earth's position in the Milky Way.

We present images and initial results from our extensive Spitzer Space Telescope imaging survey of the W5 H II region with the Infrared Array Camera (IRAC) and Multiband Imaging Photometer for Spitzer (MIPS). We detect dense clusters of stars, centered on the O stars: HD 18326, BD +60 586, HD 17505 and HD 17520. At 24 microns substantial extended emission is visible, presumably from heated dust grains that survive in the strongly ionizing environment of the H {\sc ii} region. With photometry of more than 18000 point sources, we analyze the clustering properties of objects classified as young stars by their IR spectral energy distributions (a total of 2064 sources) across the region using a minimal-spanning-tree algorithm. We find ~40--70% of infrared excess sources belong to clusters with >10 members. We find that within the evacuated cavities of the H II regions that make up W5, the ratio of Class II to Class I sources is ~7 times higher than for objects coincident with molecular gas as traced by 12CO emission and near-IR extinction maps. We attribute this contrast to an age difference between the two locations, and postulate that at least two distinct generations of star formation are visible across W5. Our preliminary analysis shows that triggering is a plausible mechanism to explain the multiple generations of star formation in W5, and merits further investigation.

We present Spitzer-IRAC, NOAO 2.1meter-Flamingos, Keck-NIRC, and FCRAO-SEQUOIA observations of the massive star forming complex S254-S258, covering an area of 25x20 arc-minutes. Using a combination of the IRAC and NIR data, we identify and classify the young stellar objects (YSO) in the complex. We detect 510 sources with near or mid IR-excess, and we classify 87 Class I, and 165 Class II sources. The YSO are found in clusters surrounded by isolated YSO in a low-density distributed population. The ratio of clustered to total YSO is 0.8. We identify six new clusters in the complex. One of them, G192.63-00, is located around the ionizing star of the HII region S255. We hypothesize that the ionizing star of S255 was formed in this cluster. We also detect a southern component of the cluster in HII region S256. The cluster G192.54-0.15, located inside HII region S254 has a VLSR of 17 km/s with respect to the main cloud, and we conclude that it is located in the background of the complex. The structure of the molecular cloud is examined using 12CO and 13CO, as well as a near-IR extinction map. The main body of the molecular cloud has VLSR between 5 and 9 km/s. The arc-shaped structure of the molecular cloud, following the border of the HII regions, and the high column density in the border of the HII regions support the idea that the material has been swept up by the expansion of the HII regions.

The mysterious very high energy gamma-ray source, TeV J2032+4130, is coincident with the powerful Cygnus OB2 stellar association, though a physical association between the two remains uncertain. It is possible that the detected very high energy photons are produced via an overdensity of locally accelerated cosmic rays impinging on molecular clouds in the source region. In order to test this hypothesis, we used the Kitt Peak 12m, the Heinrich-Hertz Submillimeter Telescope (HH-SMT), and the Five College Radio Astronomy Observatory (FCRAO), to obtain observations in the J=1-->0 and J=2-->1 lines of both 12CO and 13CO. We report here on the detection of significant molecular material toward the TeV source region which could be acting as the target of locally accelerated CRs. We also find evidence of compact molecular clumps, showing large line widths in the CO spectra, possibly indicative of energetic processes in this region of Cygnus OB2.

The FCRAO Survey of the Taurus Molecular Cloud observed the 12CO and 13CO J=1-0 emission from 98 square degrees of this important, nearby star forming region. This set of data with 45" resolution comprises the highest spatial dynamic range image of an individual molecular cloud constructed to date, and provides valuable insights to the molecular gas distribution, kinematics, and the star formation process. In this contribution, we describe the observations, calibration, data processing, and characteristics of the noise and line emission of the survey. The angular distribution of 12CO and 13CO emission over 1 km/s velocity intervals and the full velocity extent of the cloud are presented. These reveal a complex, dynamic medium of cold, molecular gas.

We report the results of a 100 square degree survey of the Taurus Molecular Cloud region in the J = 1-0 transition of 12CO and 13CO. The image of the cloud in each velocity channel includes ~ 3 million Nyquist sampled pixels on a 20" grid. The high sensitivity and large linear dynamic range of the maps in both isotopologues reveal a very complex, highly structured cloud morphology. There are large scale correlated structures evident in 13CO emission having very fine dimensions, including filaments, cavities, and rings. The 12CO emission shows a quite different structure, with particularly complex interfaces between regions of greater and smaller column density defining the boundaries of the largest-scale cloud structures. The axes of the striations seen in the 12CO emission from relatively diffuse gas are aligned with the direction of the magnetic field. Using a column density-dependent model for the CO fractional abundance, we derive the mass of the region mapped to be 24,000 solar masses, a factor of three greater than would be obtained with canonical CO abundance restricted to the high column density regions. We determine that half the mass of the cloud is in regions having column density below 2.1x10^{21} per square cm. The distribution of young stars in the region covered is highly nonuniform, with the probability of finding a star in a pixel with a specified column density rising sharply for N(H2) = 6x10^{21} cm^{-2}. We determine a relatively low star formation efficiency (mass of young stars/mass of molecular gas), between 0.3 and 1.2 %, and an average star formation rate during the past 3 Myr of 8x10^{-5} stars yr^{-1}.

Velocity anisotropy induced by MHD turbulence is investigated using computational simulations and molecular line observations of the Taurus molecular cloud. A new analysis method is presented to evaluate the degree and angle of velocity anisotropy using spectroscopic imaging data of interstellar clouds. The efficacy of this method is demonstrated on model observations derived from three dimensional velocity and density fields from the set of numerical MHD simulations that span a range of magnetic field strengths. The analysis is applied to 12CO J=1-0 imaging of a sub-field within the Taurus molecular cloud. Velocity anisotropy is identified that is aligned within 10 degrees of the mean local magnetic field direction derived from optical polarization measurements. Estimated values of the field strength based on velocity anisotropy are consistent with results from other methods. When combined with new column density measurements for Taurus, our magnetic field strength estimate indicates that the envelope of the cloud is magnetically subcritical. These observations favor strong MHD turbulence within the low density, sub-critical, molecular gas substrate of the Taurus cloud.

We use 2MASS and MSX infrared observations, along with new molecular line (CO) observations, to examine the distribution of young stellar objects (YSOs) in the molecular cloud surrounding the halo HII region KR 140 in order to determine if the ongoing star-formation activity in this region is dominated by sequential star formation within the photodissociation region (PDR) surrounding the HII region. We find that KR 140 has an extensive population of YSOs that have spontaneously formed due to processes not related to the expansion of the HII region. Much of the YSO population in the molecular cloud is concentrated along a dense filamentary molecular structure, traced by C18O, that has not been erased by the formation of the exciting O star. Some of the previously observed submillimetre clumps surrounding the HII region are shown to be sites of recent intermediate and low-mass star formation while other massive starless clumps clearly associated with the PDR may be the next sites of sequential star formation.

We present the first results from a new 250, 350, and 500 micron Galactic Plane survey taken with the Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) in 2005. This survey's primary goal is to identify and characterize high-mass proto-stellar objects (HMPOs). The region studied here covers 4 sq. deg near the open cluster NGC 6823 in the constellation Vulpecula (l=59). We find 60 compact sources (<60'' diameter) detected simultaneously in all three bands. Their spectral energy distributions (SEDs) are constrained through BLAST, IRAS, Spitzer MIPS, and MSX photometry, with inferred dust temperatures spanning ~12-40K assuming a dust emissivity index beta=1.5. The luminosity-to-mass ratio, a distance-independent quantity, spans ~0.2-130 L_\odot M_\odot^{-1}. Distances are estimated from coincident 13CO (1->0) velocities combined with a variety of other velocity and morphological data in the literature. In total, 49 sources are associated with a molecular cloud complex encompassing NGC 6823 (distance ~2.3kpc), 10 objects with the Perseus Arm (~8.5kpc) and one object is probably in the outer Galaxy (~14kpc). Near NGC 6823, the inferred luminosities and masses of BLAST sources span ~40-10^4 L_\odot, and ~15-700 M_\odot, respectively. The mass spectrum is compatible with molecular gas masses in other high-mass star forming regions. Several luminous sources appear to be Ultra Compact HII regions powered by early B stars. However, many of the objects are cool, massive gravitationally-bound clumps with no obvious internal radiation from a protostar, and hence excellent HMPO candidates.

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

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

HI 21cm-line self-absorption (HISA) reveals the shape and distribution of cold atomic clouds in the Galactic disk. Many of these clouds lack corresponding CO emission, despite being colder than purely atomic gas in equilibrium models. HISA requires background line emission at the same velocity, hence mechanisms that can produce such backgrounds. Weak, small-scale, and widespread absorption is likely to arise from turbulent eddies, while strong, large-scale absorption appears organized in cloud complexes along spiral arm shocks. In the latter, the gas may be evolving from an atomic to a molecular state prior to star formation, which would account for the incomplete HISA-CO agreement.

The observational record of turbulence within the molecular gas phase of the interstellar medium is summarized. We briefly review the analysis methods used to recover the velocity structure function from spectroscopic imaging and the application of these tools on sets of cloud data. These studies identify a near-invariant velocity structure function that is independent of local the environment and star formation activity. Such universality accounts for the cloud-to-cloud scaling law between the global line-width and size of molecular clouds found by Larson (1981) and constrains the degree to which supersonic turbulence can regulate star formation. In addition, the evidence for large scale driving sources necessary to sustain supersonic flows is summarized.

We report the detection of molecular clouds around the X-ray bright interior feature in the Galactic supernova remnant (SNR) CTB 109 (G109.1-1.0). This feature, called the Lobe, has been previously suggested to be the result of an interaction of the SNR shock wave with a molecular cloud complex. We present new high resolution X-ray data from the Chandra X-ray Observatory and new high resolution CO data from the Five College Radio Observatory which show the interaction region with the cloud complex in greater detail. The CO data reveal three clouds around the Lobe in the velocity interval -57 < v < -52 km s^-1. The velocity profiles of 12CO at various parts of the east cloud are well fit with a Gaussian; however, at the position where the CO cloud and the Lobe overlap, the velocity profile has an additional component towards higher negative velocities. The molecular hydrogen density in this part of the cloud is relatively high (N_H2 = 1.9 x 10^20 cm^-2), whereas the foreground absorption in X-rays (N_H = 4.5 x 10^21 cm^-2), obtained from Chandra data, is lower than in other parts of the cloud and in the north and south cloud. These results indicate that this cloud has been hit by the SNR blast wave on the western side, forming the bright X-ray Lobe.

The role of turbulent fragmentation in regulating the efficiency of star formation in interstellar clouds is examined from new wide field imaging of 12CO and 13CO J=1-0 emission from the Rosette and G216-2.5 molecular clouds. The Rosette molecular cloud is a typical star forming giant molecular cloud and G215-2.5 is a massive molecular cloud with no OB stars and very little low mass star formation. The properties of the turbulent gas flow are derived from the set of eigenvectors and eigenimages generated by Principal Component Analysis of the spectroscopic data cubes. While the two clouds represent quite divergent states of star formation activity, the velocity structure functions for both clouds are similar. The sonic scale, lambda_S, defined as the spatial scale at which turbulent velocity fluctuations are equivalent to the local sound speed, and the turbulent Mach number evaluated at 1 pc, M_{1pc}, are derived for an ensemble of clouds including the Rosette and, G216-2.5 regions that span a large range in star formation activity. We find no evidence for the positive correlations between these quantities and the star formation efficiency, that are predicted by turbulent fragmentation models. A correlation does exist between the star formation efficiency and the sonic scale for a subset of clouds with L_{FIR}/M(H_2) > 1 that are generating young stellar clusters. Turbulent fragmentation must play a limited and non-exclusive role in determining the yield of stellar masses within interstellar clouds.

We have observed 12CO J = 2-1 and J = 1-0, and 13CO J = 1-0 emission in two regions of HI Self-Absorption (HISA) in Perseus: a small, isolated HISA feature called the globule and a more extended HISA cloud called the complex. Using both Large Velocity Gradient and Monte Carlo radiative transfer codes we found that, in the globule, N(12CO) < 6.0x10^15 cm-2 which, using PDR models, implies that N(H_2) < 9.9x10^20 cm-2. In the complex we found that the H_2 column densities ranged from 1.2 - 2.2 x 10^21 cm-2. By comparing the HISA and CO observations we are able to constrain the physical conditions and atomic gas fraction (f). In the globule, 8 K < T_spin < 22 K and 0.02 < f < 0.2 depending on whether the (unknown) gas density is 10^2, 10^3, or 10^4 cm-3. In the complex, 12 K < T_spin < 24 K, 0.02 < f < 0.05, and we were also able to constrain the gas density (100 < n < 1200 cm-3). These results imply that the gas in the HISA clouds is colder and denser than that usually associated with the atomic ISM and, indeed, is similar to that seen in molecular clouds. The small atomic gas fractions also imply that there is a significant molecular component in these HISA clouds, even when little or no 12CO is detected. The level of 12CO detected and the visual extinction due to dust is consistent with the idea that these HISA clouds are undergoing a transition from the atomic to molecular phase.

We describe algorithms that detect 21cm line HI self-absorption (HISA) in large data sets and extract it for analysis. Our search method identifies HISA as spatially and spectrally confined dark HI features that appear as negative residuals after removing larger-scale emission components with a modified CLEAN algorithm. Adjacent HISA volume-pixels (voxels) are grouped into features in (l,b,v) space, and the HI brightness of voxels outside the 3-D feature boundaries is smoothly interpolated to estimate the absorption amplitude and the unabsorbed HI emission brightness. The reliability and completeness of our HISA detection scheme have been tested extensively with model data. We detect most features over a wide range of sizes, linewidths, amplitudes, and background levels, with poor detection only where the absorption brightness temperature amplitude is weak, the absorption scale approaches that of the correlated noise, or the background level is too faint for HISA to be distinguished reliably from emission gaps. False detection rates are very low in all parts of the parameter space except at sizes and amplitudes approaching those of noise fluctuations. Absorption measurement biases introduced by the method are generally small and appear to arise from cases of incomplete HISA detection. This paper is the third in a series examining HISA at high angular resolution. A companion paper (Paper II) uses our HISA search and extraction method to investigate the cold atomic gas distribution in the Canadian Galactic Plane Survey.

We present a 21cm line HI self-absorption (HISA) survey of cold atomic gas within Galactic longitudes 75 to 146 degrees and latitudes -3 to +5 degrees. We identify HISA as spatially and spectrally confined dark HI features and extract it from the surrounding HI emission in the arcminute-resolution Canadian Galactic Plane Survey (CGPS). We compile a catalog of the most significant features in our survey and compare our detections against those in the literature. Within the parameters of our search, we find nearly all previously detected features and identify many new ones. The CGPS shows HISA in much greater detail than any prior survey and allows both new and previously-discovered features to be placed into the larger context of Galactic structure. In space and radial velocity, faint HISA is detected virtually everywhere that the HI emission background is sufficiently bright. This ambient HISA population may arise from small turbulent fluctuations of temperature and velocity in the neutral interstellar medium. By contrast, stronger HISA is organized into discrete complexes, many of which follow a longitude-velocity distribution that suggests they have been made visible by the velocity reversal of the Perseus arm's spiral density wave. The cold HI revealed in this way may have recently passed through the spiral shock and be on its way to forming molecules and, eventually, new stars. This paper is the second in a series examining HISA at high angular resolution. A companion paper (Paper III) describes our HISA search and extraction algorithms in detail.

The universality of interstellar turbulence is examined from observed structure functions of 27 giant molecular clouds and Monte Carlo modeling. We show that the structure functions, dv=v0 l^gamma, derived from wide field imaging of CO J=1-0 emission from individual clouds are described by a narrow range in the scaling exponent, gamma, and the scaling coefficient, v0. The similarity of turbulent structure functions emphasizes the universality of turbulence in the molecular interstellar medium and accounts for the cloud-to-cloud size-line width relationship initially identified by Larson (1981). The degree of turbulence universality is quantified by Monte Carlo simulations that reproduce the mean squared velocity residuals of the observed cloud-to-cloud relationship. Upper limits to the variation of the scaling amplitudes and exponents for molecular clouds are ~10-20%. The measured invariance of turbulence for molecular clouds with vastly different sizes, environments, and star formation activity suggests a common formation mechanism such as converging turbulent flows within the diffuse ISM and a limited contribution of energy from sources within the cloud with respect to large scale driving mechanisms.

(Modified) The scaling of velocity fluctuation, dv, as a function of spatial scale L in molecular clouds can be measured from size-linewidth relations, principal component analysis, or line centroid variation. Differing values of the power law index of the scaling relation dv = L^(g3D) in 3D are given by these different methods: the first two give g3D=0.5, while line centroid analysis gives g3D=0. This discrepancy has previously not been fully appreciated, as the variation of projected velocity line centroid fluctuations (dv_{lc} = L^(g2D)) is indeed described, in 2D, by g2D=0.5. However, if projection smoothing is accounted for, this implies that g3D=0. We suggest that a resolution of this discrepancy can be achieved by accounting for the effect of density inhomogeneity on the observed g2D obtained from velocity line centroid analysis. Numerical simulations of compressible turbulence are used to show that the effect of density inhomogeneity statistically reverses the effect of projection smoothing in the case of driven turbulence so that velocity line centroid analysis does indeed predict that g2D=g3D=0.5. Using our numerical results we can restore consistency between line centroid analysis, principal component analysis and size-linewidth relations, and we derive g3D=0.5, corresponding to shock-dominated (Burgers) turbulence. We find that this consistency requires that molecular clouds are continually driven on large scales or are only recently formed.

The multivariate tool of Principal Component Analysis (PCA) is applied to 23 fields in the FCRAO CO Survey of the Outer Galaxy. PCA enables the identification of line profile differences which are assumed to be generated from fluctuations within a turbulent velocity field. The variation of these velocity differences with spatial scale within a molecular region is described by a singular power law, delta v= c L^alpha which can be used as a powerful diagnostic to turbulent motions. For the ensemble of 23 fields, we find a mean value alpha = 0.62 +- 0.11. From a recent calibration of this method using fractal Brownian motion simulations (Brunt & Heyer 2001), the measured velocity difference-size relationship corresponds to an energy spectrum, E(k), which varies as k^-beta, where beta = 2.17 +- 0.31. We compare our results to both decaying and forced hydrodynamic simulations of turbulence. We conclude that energy must be continually injected into the regions to replenish that lost by dissipative processes such as shocks. The absence of large, widely distributed shocks within the targeted fields suggests that the energy is injected at spatial scales less than several pc.

We present new, high resolution 1420 and 408 MHz continuum images and HI and 12CO (J=1-0) spectral line maps of the diffuse supernova remnant CTB104A (G93.7-0.3). Analysis of the complex continuum emission reveals no significant spectral index variations across the remnant. Three prominences around CTB104A are found to be related to the SNR, while one extension to the east is identified as an HII region associated with a background molecular shell. Small scale polarization and rotation measure (RM) structures are turbulent in nature, but we find a well-ordered RM gradient across the remnant, extending from southeast to northwest. This gradient does not agree with the direction of the global Galactic magnetic field, but does agree with a large-scale RM anomaly inferred from rotation measure data by Cleg et al. (1992). We show that the observed morphology of CTB104A is consistent with expansion in a uniform magnetic field, and this is supported by the observed RM distribution. By modeling the RM gradient with a simple compression model we have determined the magnetic field strength within the remnant as Bo ~ 2.3 micro G. We have identified signatures of the interaction of CTB104A with the surrounding neutral material, and determined its distance, from the kinematics of the HI structure encompassing the radio emission, as 1.5 kpc. We also observed clear breaks in the HI shell that correspond well to the positions of two of the prominences, indicating regions where hot gas is escaping from the interior of the SNR.

We demonstrate the capability of Principal Component Analysis (PCA) as applied by Heyer & Schloerb (1997) to extract the statistics of turbulent interstellar velocity fields as measured by the energy spectrum, E(k)= k^-beta. Turbulent velocity and density fields with known statistics are generated from fBm simulations. These fields are translated to the observational domain, T(x,y,v), considering the excitation of molecular rotational energy levels and radiative transfer. Using PCA and the characterization of velocity and spatial scales from the eigenvectors and eigenimages respectively, a relationship is identified which describes the magnitude of line profile differences, deltav, and the scale, L, over which these differences occur, deltav=L^alpha. From a series of models with varying values of beta, we find, alpha = 0.33beta - 0.05 for 1 < beta < 3. This provides the basic calibration between the intrinsic velocity field statistics to observational measures and a diagnostic for turbulent flows in the interstellar medium. We also investigate the effects of noise, line opacity, and finite resolution on these results.