Offsets of molecular line emission peaks from continuum peaks are very common but frequently difficult to explain with a single spherical cloud chemical model. We propose that the spatial projection effects of an irregular three dimensional (3D) cloud structure can be a solution. This work shows that the idea can be successfully applied to the Planck cold clump G224.4-0.6 by approximating it with four individual spherically symmetric cloud cores whose chemical patterns overlap with each other to produce observable line maps. With the empirical physical structures inferred from the observation data of this clump and a gas-grain chemical model, the four cores can satisfactorily reproduce its 850 $\mu$m continuum map and the diverse peak offsets of CCS, HC$_3$N and N$_2$H$^+$ simultaneously at chemical ages of about $8\times 10^5\sim 3\times 10^6$ yrs. The 3D projection effects on chemistry has the potential to explain such asymmetrical distributions of chemicals in many other molecular clouds.

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation.

In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt.

Many species of complex organic molecules (COMs) have been observed in several astrophysical environments but it is not clear how they are produced, particularly in cold, quiescent regions. One process that has been proposed as a means to enhance the chemical complexity of the gas phase in such regions is the explosion of the ice mantles of dust grains. In this process, a build up of chemical energy in the ice is released, sublimating the ices and producing a short lived phase of high density, high temperature gas. The gas-grain chemical code UCLCHEM has been modified to treat these explosions in order to model the observed abundances of COMs towards the TMC-1 region. It is found that, based on our current understanding of the explosion mechanism and chemical pathways, the inclusion of explosions in chemical models is not warranted at this time. Explosions are not shown to improve the model's match to the observed abundances of simple species in TMC-1. Further, neither the inclusion of surface diffusion chemistry, nor explosions, results in the production of COMs with observationally inferred abundances.

Iram 30-m Observations towards eight protostellar outflow sources were taken in the 96-\SI{176}{\giga\hertz} range. Transitions of CH$_3$OH and CH$_3$CHO were detected in seven of them. The integrated emission of the transitions of each species that fell into the observed frequency range were measured and fit using RADEX and LTE models. Column densities and gas properties inferred from this fitting are presented. The ratio of the A and E-type isomers of CH$_3$OH indicate that the methanol observed in these outflows was formed on the grain surface. Both species demonstrate a reduction of terminal velocity in their line profiles in faster outflows, indicating destruction in the post-shock gas phase. This destruction, and a near constant ratio of the CH$_3$OH and CH$_3$CHO column densities imply it is most likely that CH$_3$CHO also forms on the grain surface.

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.

Recent observations made by the Rosetta/ROSINA instrument have detected molecular oxygen in the coma of comet 67P/Churyumov-Gerasimenko with abundances at the 1-10% level relative to H2O. Previous studies have indicated that the likely origin of the O2 may be surface chemistry of primordial (dark cloud) origin, requiring somewhat warmer, denser and extreme H-atom poor conditions than are usually assumed. In this study we propose a primordial gas-phase origin for the O2 which is subsequently frozen and effectively hidden until the ice mantles are sublimated in the comet's coma. Our study presents results from a three-phase astrochemical model that simulates the chemical evolution of ices in the primordial dark cloud phase, its gravitational collapse, and evolution in the early protosolar nebula. We find that the O2 abundance can be produced and is fairly robust to the choice of the free parameters. Good matches for the O2:H2O ratio and, to a lesser extent, the N2:CO and CO:H2O ratios are obtained, but the models significantly over-produce N2. We speculate that the low value of N2:O2 that is observed is a consequence of the specific thermal history of the comet.

At the time of formation, protoplanetary discs likely contain a comparable mass to their host protostars. As a result, gravitational instabilities (GIs) are expected to play a pivotal role in the early phases of disc evolution. However, as these young objects are heavily embedded, confirmation of GIs has remained elusive. Therefore, we use the radiative transfer code LIME to produce line images of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME and explore methods to improve upon the default gridding. We synthesise noiseless observations to determine the sensitivities required to detect the spiral flux, and find that the line flux distribution does not necessarily correlate to the abundance density distribution; hence performing radiative transfer calculations is imperative. Moreover, the spiral features are seen in absorption, due to the GI-heated midplane and high extinction, which could be indicative of GI activity. If a small beamsize and appropriate molecular line are used then spatially resolving spirals in a protosolar-like disc should be possible with ALMA for an on-source time of 30 hr. Spectrally resolving non-axisymmetric structure takes only a tenth as long for a reasonable noise level, but attributing this structure to GI-induced activity would be tentative. Finally, we find that identifying finger-like features in PV diagrams of nearly edge-on discs, which are a direct indicator of spirals, is feasible with an on-source time of 19 hr, and hence likely offers the most promising means of confirming GI-driven spiral structure in young, embedded protoplanetary discs.

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.

Glycolaldehyde (HOCH$_2$CHO) and ethylene glycol ((CH$_2$OH)$_2$) are two complex organic molecules detected in the hot cores and hot corinos of several star-forming regions. The ethylene glycol/glycolaldehyde abundance ratio seems to show an increase with the source luminosity. In the literature, several surface-chemistry formation mechanisms have been proposed for these two species. With the UCLCHEM chemical code, we explored the different scenarios and compared the predictions for a range of sources of different luminosities with the observations. None of the scenarios reproduce perfectly the trend. A better agreement is, however, found for a formation through recombination of two HCO radicals followed by successive hydrogenations. The reaction between HCO and CH$_2$OH could also contribute to the formation of glycolaldehyde in addition to the hydrogenation pathway. The predictions are improved when a trend of decreasing H$_2$ density within the core region with T $\geq$ 100 K as a function of luminosity, is included in the model. Destruction reactions of complex organic molecules in the gas phase would also need to be investigated, since they can affect the abundance ratios once the species have desorbed in the warm inner regions of the star-forming regions.

In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H$_{\rm 2}$O, H$_{\rm 2}$S, HNO, N$_{\rm 2}$, NH$_{\rm 3}$, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH$_{\rm 4}$, CN, CS, H$_{\rm 2}$CO), and species which are abundant in the circumfragmentary material (e.g. HCO$^{\rm +}$). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material.

Gravitational instabilities (GIs) are most likely a fundamental process during the early stages of protoplanetary disc formation. Recently, there have been detections of spiral features in young, embedded objects that appear consistent with GI-driven structure. It is crucial to perform hydrodynamic and radiative transfer simulations of gravitationally unstable discs in order to assess the validity of GIs in such objects, and constrain optimal targets for future observations. We utilise the radiative transfer code LIME to produce continuum emission maps of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME as is and explore methods to improve upon the default gridding. We use CASA to produce synthetic observations of 270 continuum emission maps generated across different frequencies, inclinations and dust opacities. We find that the spiral structure of our protosolar-like disc model is distinguishable across the majority of our parameter space after 1 hour of observation, and is especially prominent at 230$\,$GHz due to the favourable combination of angular resolution and sensitivity. Disc mass derived from the observations is sensitive to the assumed dust opacities and temperatures, and therefore can be underestimated by a factor of at least 30 at 850$\,$GHz and 2.5 at 90$\,$GHz. As a result, this effect could retrospectively validate GIs in discs previously thought not massive enough to be gravitationally unstable, which could have a significant impact on the understanding of the formation and evolution of protoplanetary discs.

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.

Methanol (CH$_3$OH) is found to be abundant and widespread towards the Central Molecular Zone, the inner few hundred parsecs of our Galaxy. Its origin is, however, not fully understood. It was proposed that the high cosmic ray ionisation rate in this region could lead to a more efficient non-thermal desorption of this species formed on grain surfaces, but it would also mean that this species is destroyed in a relatively short timescale. In a first step, we run chemical models with a high cosmic ray ionisation rate and find that this scenario can only reproduce the lowest abundances of methanol derived in this region ($\sim$10$^{-9}$-10$^{-8}$). In a second step, we investigate another scenario based on episodic explosions of grain mantles. We find a good agreement between the predicted abundances of methanol and the observations. We find that the dominant route for the formation of methanol is through hydrogenation of CO on the grains followed by the desorption due to the grain mantle explosion. The cyclic aspect of this model can explain the widespread presence of methanol without requiring any additional mechanism. We also model silicon monoxide (SiO), another species detected in several molecular clouds of the Galactic Centre. An agreement is found with observations for a high depletion of Si (Si/H $\sim$ 10$^{-8}$) with respect to the solar abundance.

This paper presents $\textit{Herschel}$/SPIRE spectroscopic observations of several fundamental rotational ortho- and para-water transitions seen in three Jupiter-family comets and one Oort-cloud comet. Radiative transfer models that include excitation by collisions with neutrals and electrons, and by solar infrared radiation were used to produce synthetic emission line profiles originating in the cometary coma. Ortho-to-para ratios (OPRs) were determined and used to derived water production rates for all comets. Comparisons are made with the water production rates derived using an OPR of 3. The OPR of three of the comets in this study are much lower than the statistical equilibrium value of 3, however they agree with observations of comets 1P/Halley and C/2001 A2 (LINEAR), and the protoplanetary disc TW Hydrae. These results provide evidence suggesting that OPR variation is caused by post-sublimation gas-phase nuclear-spin conversion processes. The water production rates of all comets agree with previous work and, in general, decrease with increasing nucleocentric offset. This could be due to a temperature profile, additional water source, or OPR variation in the comae, or model inaccuracies.

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

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

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

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

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

To date, most simulations of the chemistry in protoplanetary discs have used 1+1D or 2D axisymmetric $\alpha$-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical and chemical evolution of a 0.17 M$_{\odot}$ self-gravitating disc over a period of 2000 yr. The 0.8 M$_{\odot}$ central protostar is likely to evolve into a solar-like star, and hence this Class 0 or early Class I young stellar object may be analogous to our early Solar System. Shocks driven by gravitational instabilities enhance the desorption rates, which dominate the changes in gas-phase fractional abundances for most species. We find that at the end of the simulation, a number of species distinctly trace the spiral structure of our relatively low-mass disc, particularly CN. We compare our simulation to that of a more massive disc, and conclude that mass differences between gravitationally unstable discs may not have a strong impact on the chemical composition. We find that over the duration of our simulation, successive shock heating has a permanent effect on the abundances of HNO, CN and NH$_3$, which may have significant implications for both simulations and observations. We also find that HCO$^+$ may be a useful tracer of disc mass. We conclude that gravitational instabilities induced in lower mass discs can significantly, and permanently, affect the chemical evolution, and that observations with high-resolution instruments such as ALMA offer a promising means of characterising gravitational instabilities in protosolar discs.

We separate the extragalactic radio source population above ~50 uJy into active galactic nuclei (AGN) and star-forming sources. The primary method of our approach is to fit the infrared spectral energy distributions (SEDs), constructed using Spitzer/IRAC and MIPS and Herschel/SPIRE photometry, of 380 radio sources in the Extended Chandra Deep Field-South. From the fitted SEDs, we determine the relative AGN and star-forming contributions to their infrared emission. With the inclusion of other AGN diagnostics such as X-ray luminosity, Spitzer/IRAC colours, radio spectral index and the ratio of star-forming total infrared flux to k-corrected 1.4 GHz flux density, qIR, we determine whether the radio emission in these sources is powered by star formation or by an AGN. The majority of these radio sources (60 per cent) show the signature of an AGN at some wavelength. Of the sources with AGN signatures, 58 per cent are hybrid systems for which the radio emission is being powered by star formation. This implies that radio sources which have likely been selected on their star formation have a high AGN fraction. Below a 1.4 GHz flux density of 1 mJy, along with finding a strong contribution to the source counts from pure star-forming sources, we find that hybrid sources constitute 20-65 per cent of the sources. This result suggests that hybrid sources have a significant contribution, along with sources that do not host a detectable AGN, to the observed flattening of the source counts at ~1mJy for the extragalactic radio source population.

We present the first Swift Ultraviolet/Optical Telescope Serendipitous Source Catalogue (UVOTSSC). The catalogue was compiled from 23,059 Swift datasets taken within the first five years of observations with the Swift UVOT. A purpose-built processing pipeline, based around the standard Swift processing tools, was employed. The catalogue contains positions, photometry in three UV and three optical bands, morphological information and data quality flags. In total, the catalogue contains 6,200,016 unique sources of which more than 2 million have multiple observations in the catalogue.

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

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

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

Low-mass stars are generally understood to form by the gravitational collapse of the dense molecular clouds known as starless cores. Continuum observations have not been able to distinguish among the several different hypotheses that describe the collapse because the predicted density distributions are the almost the same, as they are for all spherical self-gravitating clouds. However, the predicted contraction velocities are different enough that the models can be discriminated by comparing the velocities at large and small radii. This can be done by observing at least two different molecular line transitions that are excited at different densities. For example, the spectral lines of the H2O (110 - 101) and C18O (1-0) have critical densities for collisional de-excitation that differ by 5 orders of magnitude. We compare observations of these lines from the contracting starless core L1544 against the spectra predicted for several different hypothetical models of contraction including the Larson-Penston flow, the inside-out collapse of the singular isothermal sphere, the quasi-equilibrium contraction of an unstable Bonnor-Ebert sphere, and the non-equilibrium collapse of an over-dense Bonnor-Ebert sphere. Only the model of the unstable quasi-equilibrium Bonnor-Ebert sphere is able to produce the observed shapes of both spectral lines. This model allows us to interpret other observations of molecular lines in L1544 to find that the inward velocities seen in observations of CS(2-1) and N2H+ are located within the starless core itself, in particular in the region where the density profile follows an inverse square law. If these conclusions were to hold in the analysis of other starless cores, this would imply that the formation of hydrostatic clouds within the turbulent interstellar medium is not only possible but not exceptional and may be an evolutionary phase in low-mass star formation.

A rich variety of molecular species has now been observed towards hot cores in star forming regions and in the interstellar medium. An increasing body of evidence from millimetre interferometers suggests that many of these form at the interfaces between protostellar outflows and their natal molecular clouds. However, current models have remained unable to explain the origin of the observational bias towards wide-angled "brandy snifter" shaped outflows over narrower "champagne flute" shapes in carbon monoxide imaging. Furthermore, these wide-angled systems exhibit unusually high abundances of the molecular ion HCO$^+$. We present results from a chemo-dynamic model of such regions where a rich chemistry arises naturally as a result of turbulent mixing between cold, dense molecular gas and the hot, ionized outflow material. The injecta drives a rich and rapid ion-neutral chemistry in qualitative and quantitative agreement with the observations. The observational bias towards wide-angled outflows is explained naturally by the geometry-dependent ion injection rate causing rapid dissociation of CO in the younger systems.

The Herschel Space Observatory's recent detections of water vapor in the cold, dense cloud L1544 allow a direct comparison between observations and chemical models for oxygen species in conditions just before star formation. We explain a chemical model for gas phase water, simplified for the limited number of reactions or processes that are active in extreme cold ($<$ 15 K). In this model, water is removed from the gas phase by freezing onto grains and by photodissociation. Water is formed as ice on the surface of dust grains from O and OH and released into the gas phase by photodesorption. The reactions are fast enough with respect to the slow dynamical evolution of L1544 that the gas phase water is in equilibrium for the local conditions thoughout the cloud. We explain the paradoxical radiative transfer of the H$_2$O ($1_{10}-1_{01}$) line. Despite discouragingly high optical depth caused by the large Einstein A coefficient, the subcritical excitation in the cold, rarefied H$_2$ causes the line brightness to scale linearly with column density. Thus the water line can provide information on the chemical and dynamical processes in the darkest region in the center of a cold, dense cloud. The inverse P-Cygni profile of the observed water line generally indicates a contracting cloud. This profile is reproduced with a dynamical model of slow contraction from unstable quasi-static hydrodynamic equilibrium (an unstable Bonnor-Ebert sphere).

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

Low mass star-forming regions are more complex than the simple spherically symmetric approximation that is often assumed. We apply a more realistic infall/outflow physical model to molecular/continuum observations of three late Class 0 protostellar sources with the aims of (a) proving the applicability of a single physical model for all three sources, and (b) deriving physical parameters for the molecular gas component in each of the sources. We have observed several molecular species in multiple rotational transitions. The observed line profiles were modelled in the context of a dynamical model which incorporates infall and bipolar outflows, using a three dimensional radiative transfer code. This results in constraints on the physical parameters and chemical abundances in each source. Self-consistent fits to each source are obtained. We constrain the characteristics of the molecular gas in the envelopes as well as in the molecular outflows. We find that the molecular gas abundances in the infalling envelope are reduced, presumably due to freeze-out, whilst the abundances in the molecular outflows are enhanced, presumably due to dynamical activity. Despite the fact that the line profiles show significant source-to-source variation, which primarily derives from variations in the outflow viewing angle, the physical parameters of the gas are found to be similar in each core.

Complex organic molecules, such as propylene (CH3CHCH2), are detected in molecular clouds (such as TMC1) with high fractional abundances (~2.10^{-9}, relative to hydrogen) that cannot be explained by gas-phase chemical reactions under normal dark cloud conditions. To obtain such high abundances requires an efficient grain-driven chemistry to be operating, coupled with effective desorption of the complex organics back into the gas-phase. We propose that the mechanism that we have previously described; rapid high density gas-phase chemistry in the gas released following sudden, total, ice mantle sublimation - can explain the high abundances, without recourse to ill-defined surface chemical pathways. Our model provides a natural explanation for why it is that some sources harbour high abundances of propylene, whilst others do not; based on the age and level of dynamical activity within the source (which affects the ice composition) and the chemical composition of the ambient gas.

The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. ALMA will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here we present simulated ALMA observations based on radiative transfer modelling of a relatively massive (0.39 M_solar) self-gravitating disc embedded in a 10 M_solar dense core, with structure similar to the pre-stellar core L1544. We focus on simple species and conclude that C17O 3-2, HCO+ 3-2, OCS 26-25 and H2CO 404-303 lines can be used to probe the disc structure and kinematics at all scales.

We have carried out an investigation of the chemical evolution of gas in different carbon-rich circumstellar environments. Previous studies have tended to invoke terrestrial flame chemistries, based on acetylene (C2H2) combustion to model the formation of carbon dust, via Polycyclic Aromatic Hydrocarbons (PAHs). In this work we pay careful attention to the accurate calculation of the molecular photoreaction rate coefficients to ascertain whether there is a universal formation mechanism for carbon dust in strongly irradiated astrophysical environments. A large number of possible chemical channels may exist for the formation of PAHs, so we have concentrated on the viability of the formation of the smallest building block species, C2H2, in a variety of carbon-rich stellar outflows. C2H2 is very sensitive to dissociation by UV radiation. This sensitivity is tested, using models of the time-dependent chemistry. We find that C2H2 formation is sensitive to some of the physical parameters and that in some known sources of dust-formation it can never attain appreciable abundances. Therefore multiple (and currently ill-defined) dust-formation channels must exist.

We present a model for the formation of large organic molecules in dark clouds. The molecules are produced in the high density gas-phase that exists immediately after ice mantles are explosively sublimated. The explosions are initiated by the catastrophic recombination of trapped atomic hydrogen. We propose that, in molecular clouds, the processes of freeze-out onto ice mantles, accumulation of radicals, explosion and then rapid (three-body) gas-phase chemistry occurs in a cyclic fashion. This can lead to a cumulative molecular enrichment of the interstellar medium. A model of the time-dependent chemistries, based on this hypothesis, shows that significant abundances of large molecular species can be formed, although the complexity of the species is limited by the short expansion timescale in the gas, immediately following mantle explosion. We find that this mechanism may be an important source of smaller organic species, such as methanol and formaldehyde, as well as precursors to bio-molecule formation. Most significantly, we predict the gas-phase presence of these larger molecular species in quiescent molecular clouds and not just dynamically active regions, such as hot cores. As such the mechanism that we propose complements alternative methods of large molecule formation, such as those that invoke solid-state chemistry within activated ice mantles.

We present the mid-infrared spectra of seven of the most powerful radio-galaxies known to exist at 1.5 < z < 2.6. The radio emission of these sources is dominated by the AGN with 500 MHz luminosities in the range 10^27.8 - 10^29.1 W/Hz. The AGN signature is clearly evident in the mid-infrared spectra, however, we also detect polycyclic aromatic hydrocarbons emission, indicative of prodigious star formation at a rate of up to ~1000 Msun/yr. Interestingly, we observe no significant correlation between AGN power and star formation in the host galaxy. We also find most of these radio galaxies to have weak 9.7 um silicate absorption features (tau_{9.7 um} < 0.8) which implies that their mid-IR obscuration is predominantly due to the dusty torus that surrounds the central engine, rather than the host galaxy. The tori are likely to have an inhomogeneous distribution with the obscuring structure consisting of individual clouds. We estimate that these radio galaxies have already formed the bulk of their stellar mass and appear to lie at a stage in their evolution where the obscured AGN dominates the energy output of the system but star formation is also prevalent.

We investigate the importance of the shielding of chemical photorates by molecular hydrogen photodissociation lines and the carbon photoionization continuum deep within models of photon dominated regions. In particular, the photodissociation of N2 and CN are significantly shielded by the H2 photodissociation line spectrum. We model this by switching off the photodissociation channels for these species behind the HI to H2 transition. We also model the shielding effect of the carbon photoionization continuum as an attenuation of the incident radiation field shortwards of 1102\AA. Using recent line and continuum cross section data, we present calculations of the direct and cosmic ray induced photorates for a range of species, as well as optically thick shielding factors for the carbon continuum. Applying these to a time dependent PDR model we see enrichments in the abundances of N2, N2H+, NH3 and CN by factors of roughly 3-100 in the extinction band Av=2.0 to Av=4.0 for a range of environments. While the precise quantitative results of this study are limited by the simplicity of our model, they highlight the importance of these mutual shielding effects, neither of which has been discussed in recent models.

We present a detailed study of the infrared spectral energy distribution of the high-redshift radio galaxy MRC 1138-26 at z = 2.156, also known as the Spiderweb Galaxy. By combining photometry from Spitzer, Herschel and LABOCA we fit the rest-frame 5-300 um emission using a two component, starburst and active galactic nucleus (AGN), model. The total infrared (8 - 1000 um) luminosity of this galaxy is (1.97+/-0.28)x10^13 Lsun with (1.17+/-0.27) and (0.79+/-0.09)x10^13 Lsun due to the AGN and starburst components respectively. The high derived AGN accretion rate of \sim20% Eddington, and the measured star formation rate (SFR) of 1390pm150 Msun/yr, suggest that this massive system is in a special phase of rapid central black hole and host galaxy growth, likely caused by a gas rich merger in a dense environment. The accretion rate is sufficient to power both the jets and the previously observed large outflow. The high SFR and strong outflow suggest this galaxy could potentially exhaust its fuel for stellar growth in a few tens of Myr, although the likely merger of the radio galaxy with nearby satellites suggest bursts of star formation may recur again on time scales of several hundreds of Myr. The age of the radio lobes implies the jet started after the current burst of star formation, and therefore we are possibly witnessing the transition from a merger-induced starburst phase to a radio-loud AGN phase. We also note tentative evidence for [CII]158um emission. This paper marks the first results from the Herschel Galaxy Evolution Project (Project HeRGE), a systematic study of the evolutionary state of 71 high redshift, 1 < z < 5.2, radio galaxies.

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

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

Until now, axisymmetric, alpha-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surrounded by a massive (0.39 Msun), self-gravitating disc, similar to those that may be found around Class 0 and early Class I protostars, in a study of disc chemistry. We find that for each of a number of species, e.g. H2O, adsorption and desorption dominate the changes in the gas-phase fractional abundance; because the desorption rates are very sensitive to temperature, maps of the emissions from such species should reveal the locations of shocks of varying strengths. The gas-phase fractional abundances of some other species, e.g. CS, are also affected by gas-phase reactions, particularly in warm shocked regions. We conclude that the dynamics of massive discs have a strong impact on how they appear when imaged in the emission lines of various molecular species.

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

We suggest that three-body chemistry may occur in warm high density gas evaporating in transient co\textendash desorption events on interstellar ices. Using a highly idealised computational model we explore the chemical conversion from simple species of the ice to more complex species containing several heavy atoms, as a function of density and of adopted three body rate coefficients. We predict that there is a wide range of densities and rate coefficients in which a significant chemical conversion may occur. We discuss the implications of this idea for the astrochemistry of hot cores.

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We have coupled a chemical model with two dynamical models of collapsing low mass star-forming cores to predict abundances across the core of the commonly used infall tracers, CS and HCO$^+$, at various stages of the collapse. The models investigated are a new ambipolar diffusion model and the `inside-out' collapse model. We have then used these results as an input to a radiative transfer model to predict the line profiles of several transitions of these molecules. For the inside-out collapse model, we predict significant molecular depletion due to freeze-out in the core centre, which prevents the formation of the blue asymmetry (believed to be the `signature' of infall) in the line profiles. Molecular depletion also occurs in the ambipolar diffusion model during the late stages of collapse, but the line profiles still exhibit a strong blue asymmetry due to extended infall. For the inside-out collapse model to exhibit the blue asymmetry it is necessary to impose a negative kinetic temperature gradient on the core and suppress freeze-out. Since freeze-out is observed in several class 0 protostars which are thought to be collapsing, this presents a major inconsistency in the inside-out collapse model of star formation.

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

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

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.

Photon-dominated regions (PDRs) are powerful molecular line emitters in external galaxies. They are expected in galaxies with high rates of massive star formation due to either starburst (SB) events or starburst coupled with active galactic nuclei (AGN) events. We have explored the PDR chemistry for a range of physical conditions representing a variety of galaxy types. Our main result is a demonstration of the sensitivity of the chemistry to changes in the physical conditions. We adopt crude estimates of relevant physical parameters for several galaxy types and use our models to predict suitable molecular tracers of those conditions. The set of recommended molecular tracers differs from that which we recommended for use in galaxies with embedded massive stars. Thus, molecular observations can in principle be used to distinguish between excitation by starburst and by SB+AGN in distant galaxies. Our recommendations are intended to be useful in preparing Herschel and ALMA proposals to identify sources of excitation in galaxies.

We investigate the physical and chemical conditions necessary for low-mass star formation in extragalactic environments by calculating various characteristic timescales associated with star formation for a range of initial conditions. The balance of these timescales indicates whether low-mass star formation is enhanced or inhibited under certain physical conditions. In this study, we consider timescales for free-fall, cooling, freeze-out, desorption, chemistry and ambipolar diffusion and their variations with changes in the gas density, metallicity, cosmic ray ionisation rate and FUV radiation field strength. We find that extragalactic systems with high FUV radiation field strengths and high cosmic ray fluxes considered at a range of metallicities, are likely to have enhanced low-mass star formation unless the magnetic pressure is sufficient to halt collapse. Our results indicate that this is only likely to be the case for high-redshift galaxies approaching solar metallicities. Unless this is true for all high-redshift sources, this study finds little evidence for a high-mass biased IMF at high redshifts.

We present the first detection of CS in the Antennae galaxies towards the NGC 4038 nucleus, as well as the first detections of two high-J (5-4 and 7-6) CS lines in the center of M 82. The CS(7-6) line in M 82 shows a profile that is surprisingly different to those of other low-J CS transitions we observed. This implies the presence of a separate, denser and warmer molecular gas component. The derived physical properties and the likely location of the CS(7-6) emission suggests an association with the supershell in the centre of M 82.

The evolution of star-forming core analogues undergoing inside-out collapse is studied with a multi-point chemodynamical model which self-consistently computes the abundance distribution of chemical species in the core. For several collapse periods the output chemistry of infall tracer species such as HCO+, CS, and N2H+, is then coupled to an accelerated Lambda-iteration radiative transfer code, which predicts the emerging molecular line profiles using two different input gas/dust temperature distributions. We investigate the sensitivity of the predicted spectral line profiles and line asymmetry ratios to the core temperature distribution, the time-dependent model chemistry, as well as to ad hoc abundance distributions. The line asymmetry is found to be strongly dependent on the adopted chemical abundance distribution. In general, models with a warm central region show higher values of blue asymmetry in optically thick HCO+ and CS lines than models with a starless core temperature profile. We find that in the formal context of Shu-type inside-out infall, and in the absence of rotation or outflows, the relative blue asymmetry of certain HCO+ and CS transitions is a function of time and, subject to the foregoing caveats, can act as a collapse chronometer. The sensitivity of simulated HCO+ line profiles to linear radial variations, subsonic or supersonic, of the internal turbulence field is investigated in the separate case of static cores.

Hot core molecules should be detectable in external active galaxies out to high redshift. We present here a detailed study of the chemistry of star-forming regions under physical conditions that differ significantly from those likely to be appropriate in the Milky Way Galaxy. We examine, in particular, the trends in molecular abundances as a function of time with respect to changes in the relevant physical parameters. These parameters include metallicity, dust:gas mass ratio, the H$_{2}$ formation rate, relative initial elemental abundances, the cosmic ray ionization rate, and the temperature of hot cores. These trends indicate how different tracers provide information on the physical conditions and on evolutionary age. We identify hot core tracers for several observed galaxies that are considered to represent spirals, active galaxies, low-metallicity galaxies, and high-redshift galaxies. Even in low-metallicity examples, many potential molecular tracers should be present at levels high enough to allow unresolved detection of active galaxies at high redshift containing large numbers of hot cores.

CO isotopes are able to probe the different components in protostellar clouds. These components, core, envelope and outflow have distinct physical conditions and sometimes more than one component contributes to the observed line profile. In this study we determine how CO isotope abundances are altered by the physical conditions in the different components. We use a 3D molecular line transport code to simulate the emission of four CO isotopomers, 12CO J=2-1, 13CO J=2-1, C18O J=2-1 and C17O J=2-1 from the Class 0/1 object L483, which contains a cold quiescent core, an infalling envelope and a clear outflow. Our models replicate JCMT (James Clerk Maxwell Telescope) line observations with the inclusion of freeze-out, a density profile and infall. Our model profiles of 12CO and 13CO have a large linewidth due to a high velocity jet. These profiles replicate the process of more abundant material being susceptible to a jet. C18O and C17O do not display such a large linewidth as they trace denser quiescent material deep in the cloud.

The desorption of molecular species from ice mantles back into the gas phase in molecular clouds results from a variety of very poorly understood processes. We have investigated three mechanisms; desorption resulting from H_2 formation on grains, direct cosmic ray heating and cosmic ray induced photodesorption. Whilst qualitative differences exist between these processes (essentially deriving from the assumptions concerning the species-selectivity of the desorption and the assumed threshold adsorption energies, E_t) all three processes are found to be potentially very significant in dark cloud conditions. It is therefore important that all three mechanisms should be considered in studies of molecular clouds in which freeze-out and desorption are believed to be important. Employing a chemical model of a typical static molecular core and using likely estimates for the quantum yields of the three processes we find that desorption by H_2 formation probably dominates over the other two mechanisms. However, the physics of the desorption processes and the nature of the dust grains and ice mantles are very poorly constrained. We therefore conclude that the best approach is to set empirical constraints on the desorption, based on observed molecular depletions - rather than try to establish the desorption efficiencies from purely theoretical considerations. Applying this method to one such object (L1689B) yields upper limits to the desorption efficiencies that are consistent with our understanding of these mechanisms.

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.

Stars form in the densest, coldest, most quiescent regions of molecular clouds. Molecules provide the only probes which can reveal the dynamics, physics, chemistry and evolution of these regions, but our understanding of the molecular inventory of sources and how this is related to their physical state and evolution is rudimentary and incomplete. The Spectral Legacy Survey (SLS) is one of seven surveys recently approved by the JCMT Board. Starting in 2007, the SLS will produce a spectral imaging survey of the content and distribution of all the molecules detected in the 345 GHz atmospheric window (between 332 GHz and 373 GHz) towards a sample of 5 sources. Our intended targets are: a low mass core (NGC1333 IRAS4), 3 high mass cores spanning a range of star forming environments and evolutionary states (W49, AFGL2591, and IRAS20126), and a PDR (the Orion Bar). The SLS will use the unique spectral imaging capabilities of HARP-B/ACSIS to study the molecular inventory and the physical structure of these objects, which span different evolutionary stages and physical environments, to probe their evolution during the star formation process. As its name suggests, the SLS will provide a lasting data legacy from the JCMT that is intended to benefit the entire astronomical community. As such, the entire data set (including calibrated spectral datacubes, maps of molecular emission, line identifications, and calculations of the gas temperature and column density) will be publicly available.

We investigate the chemical and observational implications of repetitive transient dense core formation in molecular clouds. We allow a transient density fluctuation to form and disperse over a period of 1 Myr, tracing its chemical evolution. We then allow the same gas immediately to undergo further such formation and dispersion cycles. The chemistry of the dense gas in subsequent cycles is similar to that of the first, and a limit cycle is reached quickly (2 - 3 cycles). Enhancement of hydrocarbon abundances during a specific period of evolution is the strongest indicator of previous dynamical history. The molecular content of the diffuse background gas in the molecular cloud is expected to be strongly enhanced by the core formation and dispersion process. Such enhancement may remain for as long as 0.5 Myr. The frequency of repetitive core formation should strongly determine the level of background molecular enhancement. We also convolve the emission from a synthesised dark cloud, comprised of ensembles of transient dense cores. We find that the dynamical history of the gas, and therefore the chemical state of the diffuse inter-core medium, may be determined if a sufficient sample of cores is present in an ensemble. Molecular ratios of key hydrocarbons with SO and SO2 are crucial to this distinction. Only surveys with great enough angular resolution to resolve individual cores, or very small groupings, are expected to show evidence of repetitive dynamical processing. The existence of non-equilibrium chemistry in the diffuse background may have implications for the initial conditions used in chemical models. Observed variations in the chemistries of diffuse and translucent regions may be explained by lines of sight which intersect a number of molecular cloud cores in various stages of evolution.

New molecular line observations of the Bok globule Barnard 68 in HCO+ irrefutably confirm the complex pattern of red and blue asymmetric line profiles seen across the face of the cloud in previous observations of CS. The new observations thus strengthen the previous interpretation that Barnard 68 is undergoing peculiar oscillations. Furthermore, the physical chemistry of B68 indicates that the object is much older than the sound crossing time and is therefore long-lived. A model is presented for the globule in which a modest external pressure perturbation is shown to lead to oscillations about a stable equilibrium configuration. Such oscillations may be present in other stable starless cores as manifested by a similar signature of inward and outward motions.

The ionization fraction is an important factor in determining the chemical and physical evolution of star forming regions. In the dense, dark starless cores of such objects, the ionization rate is dominated by cosmic rays; it is therefore possible to use simple analytic estimators, based on the relative abundances of different molecular tracers, to determine the cosmic ray ionization rate. This paper uses a simple model to investigate the accuracy of two well-known estimators in dynamically evolving molecular clouds. It is found that, although the analytical formulae based on the abundances of H3+,H2,CO,O,H2O and HCO+ give a reasonably accurate measure of the cosmic ray ionization rate in static, quiescent clouds, significant discrepancies occur in rapidly evolving (collapsing) clouds. As recent evidence suggests that molecular clouds may consist of complex, dynamically evolving sub-structure, we conclude that simple abundance ratios do not provide reliable estimates of the cosmic ray ionization rate in dynamically active regions.

We construct models of molecular clouds that are considered as ensembles of transient cores. Each core is assumed to develop in the background gas of the cloud, grow to high density and decay into the background. The chemistry in each core responds to the dynamical state of the gas and to the gas-dust interaction. Ices are deposited on the dust grains in the core's dense phase, and this material is returned to the gas as the core expands to low density. The cores of the ensemble number typically one thousand and are placed randomly in position within the cloud, and are assigned a random evolutionary phase. The models are used to generate molecular line contour maps of a typical dark cloud. These maps are found to represent extremely well the characteristic features of observed maps of the dark cloud L673, which has been observed at both low and high resolutions. The computed maps are found to exhibit the general morphology of the observed maps, and to generate similar sizes of emitting regions, molecular column densities, and the separations between peaks of emissions of various molecular species. The models give insight into the nature of molecular clouds and the dynamical processes occurring within them, and significantly constrain dynamical and chemical processes in the interstellar medium.

Lynds 1498 is a pre-protostellar core (PPC) and was one of the initial objects toward which molecular depletion and differentiation was detected. Despite the considerable scrutiny of L1498, there has not been a extensive study of the density and temperature structure as derived from radiative transfer modeling of dust continuum observations. We present deep SCUBA observations of L1498 at 850 and 450 micron, high resolution BEARS maps of the N2H+ 1-0 transition, CSO observations of the N2H+ 3-2 transition, and GBT observations of the C3S 4-3 transition. We also present a comparison of derived properties between L1498 and nearby PPCs that have been observed at far-infrared and submillimeter wavelengths. We present a more realistic treatment of PPC heating which varies the strength of the ISRF, Sisrf, and includes attenuation of the ISRF due to dust grains at the outer radius of the core, Av. The best-fitted model consists of a Bonner-Ebert sphere with a central density of 1 - 3 x 10^4 cm-3, R_o ~ 0.29 pc, 0.5 <= Sisrf <= 1, Av ~ 1 mag, and a nearly isothermal temperature profile of ~ 10.5 K for OH8 opacities. C3S emission shows a central depletion hole while N2H+ emission is centrally peaked. The observed depletions of C3S and H2CO, the modest N2H+ abundance, and a central density that is an order of magnitude lower than other modeled PPCs suggests that L1498 may be a forming PPC. Our derived temperature and density profile will improve modeling of molecular line observations that will explicate the core's kinematical and chemical state. (abridged)

Nova Cassiopeiae 1993 (V705 Cas) was an archetypical dust-forming nova. It displayed a deep minimum in the visual light curve, and spectroscopic evidence for carbon, hydrocarbon and silicate dust. We report the results of fitting the infrared spectral energy distribution with the DUSTY code, which we use to determine the properties and geometry of the emitting dust. The emission is well described as originating in a thin shell whose dust has a carbon:silicate ratio of ~2:1 by number (1.26:1 by mass) and a relatively flat size distribution. The 9.7micron and 18micron silicate features are consistent with freshly-condensed dust and, while the lower limit to the grain size distribution is not well constrained, the largest grains have dimensions \~0.06micron; unless the grains in V705 Cas were anomalously small, the sizes of grains produced in nova eruptions may previously have been overestimated in novae with optically thick dust shells. Laboratory work by Grishko & Duley may provide clues to the apparently unique nature of nova UIR features.

The very high rates of second generation star formation detected and inferred in high redshift objects should be accompanied by intense millimetre-wave emission from hot core molecules. We calculate the molecular abundances likely to arise in hot cores associated with massive star formation at high redshift, using several independent models of metallicity in the early Universe. If the number of hot cores exceeds that in the Milky Way Galaxy by a factor of at least one thousand, then a wide range of molecules in high redshift hot cores should have detectable emission. It should be possible to distinguish between independent models for the production of metals and hence hot core molecules should be useful probes of star formation at high redshift.

We present observations of 25 transitions of 17 isotopologues of 9 molecules toward B335. With a goal of constraining chemical models of collapsing clouds, we compare our observations, along with data from the literature, to models of chemical abundances. The observed lines are simulated with a Monte Carlo code, which uses various physical models of density and velocity as a function of radius. The dust temperature as a function of radius is calculated self-consistently by a radiative transfer code. The gas temperature is then calculated at each radius, including gas-dust collisions, cosmic rays, photoelectric heating, and molecular cooling. The results provide the input to the Monte Carlo code. We consider both {\it ad hoc} step function models for chemical abundances and abundances taken from a self-consistent modeling of the evolution of a star-forming core. The step function models can match the observed lines reasonably well, but they require very unlikely combinations of radial variations in chemical abundances. Among the self-consistent chemical models, the observed lines are matched best by models with somewhat enhanced cosmic-ray ionization rates and sulfur abundances. We discuss briefly the steps needed to close the loop on the modeling of dust and gas, including off-center spectra of molecular lines.

Recent observations suggest that the behaviour of tracer species such as N_2H+ and CS is significantly different in regions of high and low mass star formation. In the latter, N_2H+ is a good tracer of mass, while CS is not. Observations show the reverse to be true in high-mass star formation regions. We use a computational chemical model to show that the abundances of these and other species may be significantly altered by a period of accelerated collapse in high mass star forming regions. We suggest these results provide a potential explanation of the observations, and make predictions for the behaviour of other species.

The search for the onset of star formation in pre-stellar cores has focussed on the identification of an infall signature in the molecular line profiles of tracer species. The classic infall signature is a double peaked line profile with an asymmetry in the strength of the peaks such that the blue peak is stronger. L1689B is a pre-stellar core and infall candidate but new JCMT HCO+ line profile data, presented here, confirms that both blue and red asymmetric line profiles are present in this source. Moreover, a dividing line can be drawn between the locations where each type of profile is found. It is argued that it is unlikely that the line profiles can be interpreted with simple models of infall or outflow and that rotation of the inner regions is the most likely explanation. A rotational model is developed in detail with a new 3D molecular line transport code and it is found that the best type of model is one in which the rotational velocity profile is in between solid body and Keplerian. It is firstly shown that red and blue asymmetric line profiles can be generated with a rotation model entirely in the absence of any infall motion. The model is then quantitively compared with the JCMT data and an iteration over a range of parameters is performed to minmize the difference between the data and model. The results indicate that rotation can dominate the line profile shape even before the onset of infall.

In cold molecular clouds submillimetre emission lines are excited by the ambient radiation field. The pumping is dominated by the cosmic microwave background (CMB). It is usual in molecular line radiative transfer modelling to simply assume that this is the only incident radiation field. In this paper, a molecular line transport code and a dust radiative transfer code are used to explore the effects of the inclusion of a full interstellar radiation field (ISRF) on a simple test molecular cloud. It is found that in many galactic situations, the shape and strength of the line profiles that result are robust to variations in the ISRF and thus that in most cases, it is safe to adopt the CMB radiation field for the molecular line transport calculations. However, we show that in two examples, the inclusion of a plausible radiation field can have a significant effect on the the line profiles. Firstly, in the vicinity of an embedded massive star, there will be an enhanced far infared component to the radiation field. Secondly, for molecular clouds at large redshift, the CMB temperature increases and this of course also alters the radiation field. In both of these cases, the line profiles are weakened significantly compared to a cloud exposed to a standard radiation field. Therefore this effect should be accounted for when investigating prestellar cores in massive star forming regions and when searching for molecular clouds at high redshift.

A plausible model is proposed for the enhancement of the abundance of molecular species in bipolar outflow sources. In this model, levels of HCO+ enhancement are considered based on previous chemical calculations, that are assumed to result from shock-induced desorption and photoprocessing of dust grain ice mantles in the boundary layer between the outflow jet and the surrounding envelope. A radiative transfer simulation that incorporates chemical variations within the flow shows that the proposed abundance enhancements in the boundary layer are capable of reproducing the observed characteristics of the outflow seen in HCO+ emission in the star forming core L1527. The radiative transfer simulation also shows that the emission lines from the enhanced molecular species that trace the boundary layer of the outflow exhibit complex line profiles indicating that detailed spatial maps of the line profiles are essential in any attempt to identify the kinematics of potential infall/outflow sources. This study is one of the first applications of a full three dimensional radiative transfer code which incorporates chemical variations within the source.

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.

We present results from the observations and modeling of seventeen Class I cores with the Submillimetre Common Users Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). By modeling the transfer of radiation through the envelope for nine cores, we find, for a power law distribution n(r)=n_f(r/r_f)^-p, the average and standard deviation p=1.6 +/- 0.4 and a median of p=1.8. However, the inclusion of a disk or other point-like component can cause the derived p to be shallower by as much as 0.5. In addition, we test the Shu collapse model for our sources and discuss the application of simpler analyses that derive a density power law distribution directly from the slope of the intensity radial profile. The total mass of the envelope in our sample has a range from 0.04 to 5.0 M_sun, but these masses disagree with the virial masses derived from molecular line observations. Finally, we discuss the nature of these sources in light of various evolutionary indicators and find that T_bol and L_obs/L_smm are often inconsistent in distinguishing Class 0 from Class I cores.

Seven Class 0 sources mapped with SCUBA at 850 and 450 micron are modeled using a one dimensional radiative transfer code. The modeling takes into account heating from an internal protostar, heating from the ISRF, realistic beam effects, and chopping to model the normalized intensity profile and spectral energy distribution. Power law density models, n(r) ~ r^{-p}, fit all of the sources; best fit values are mostly p = 1.8 +/- 0.1, but two sources with aspherical emission contours have lower values (p ~ 1.1). Including all sources, <p> = 1.63 +/- 0.33. Based on studies of the sensitivity of the best-fit p to variations in other input parameters, uncertainties in p for an envelope model are \Delta p = +/- 0.2. If an unresolved source (e.g., a disk) contributes 70% of the flux at the peak, p is lowered in this extreme case and \Delta p = ^{+0.2}_{-0.6}. The models allow a determination of the internal luminosity (<L_{int}> = 4.0 \lsun) of the central protostar as well as a characteristic dust temperature for mass determination (<T_{iso}> = 13.8 +/- 2.4 K). We find that heating from the ISRF strongly affects the shape of the dust temperature profile and the normalized intensity profile, but does not contribute strongly to the overall bolometric luminosity of Class 0 sources. There is little evidence for variation in the dust opacity as a function of distance from the central source. The data are well-fitted by dust opacities for coagulated dust grains with ice mantles (Ossenkopf & Henning 1994). The density profile from an inside-out collapse model (Shu 1977) does not fit the data well, unless the infall radius is set so small as to make the density nearly a power-law.

We have modeled the emission from dust in pre-protostellar cores, including a self-consistent calculation of the temperature distribution for each input density distribution. Model density distributions include Bonnor-Ebert spheres and power laws. The Bonnor-Ebert spheres fit the data well for all three cores we have modeled. The dust temperatures decline to very low values ($\Td \sim 7$ K) in the centers of these cores, strongly affecting the dust emission. Compared to earlier models that assume constant dust temperatures, our models indicate higher central densities and smaller regions of relatively constant density. Indeed, for L1544, a power-law density distribution, similar to that of a singular, isothermal sphere, cannot be ruled out. For the three sources modeled herein, there seems to be a sequence of increasing central condensation, from L1512 to L1689B to L1544. The two denser cores, L1689B and L1544, have spectroscopic evidence for contraction, suggesting an evolutionary sequence for pre-protostellar cores.

We have obtained 850 and 450 micron continuum maps of 21 low mass cores with SED's ranging from pre-protostellar to Class I (18K < Tbol < 370K), using SCUBA at the JCMT. In this paper we present the maps, radial intensity profiles, and photometry. Pre-protostellar cores do not have power-law intensity profiles, whereas the intensity profiles of Class 0 and Class I sources can be fitted with power laws over a large range of radii. A substantial number of sources have companion sources within a few arcminutes (2 out of 5 pre-protostellar cores, 9 out of 16 Class 0/I sources). The mean separation between sources is 10800 AU. The median separation is 18000 AU including sources without companions as a lower limit. The mean value of the spectral index between 450 and 850 microns is 2.8 +- 0.4, with PPCs having slightly lower spectral indices (2.5 +- 0.4). The mean mass of the sample, based on the dust emission in a 120 arcsecond aperture, is 1.1 +- 0.9 solar masses. For the sources fitted by power-law intensity distributions (I_nu(b)/I_nu(0) = (b/b_0) ^ m), the mean value of m is 1.52 +- 0.45 for Class 0 and I sources at 850 microns and 1.44 +- 0.25 at 450 microns. Based on a simple analysis, assuming the emission is in the Rayleigh-Jeans limit and that T_d(r) ~ r ^ {-0.4}, these values of m translate into power-law density distributions (n ~ r ^ {-p}) with p ~ 2.1. However, we show that this result may be changed by more careful consideration of effects such as beam size and shape, finite outer radii, more realistic T_d(r), and failure of the Rayleigh-Jeans approximation.