A submillimetre survey of the kinematics of the Perseus molecular cloud – I. Data

We present submillimetre observations of the   J = 3 → 2  rotational transition of ¹²CO, ¹³CO and C¹⁸O across over 600 arcmin² of the Perseus molecular cloud, undertaken with the Heterodyne Array Receiver Programme (HARP), a new array spectrograph on the James Clerk Maxwell Telescope. The data encompass four regions of the cloud, containing the largest clusters of dust continuum condensations: NGC 1333, IC348, L1448 and L1455. A new procedure to remove striping artefacts from the raw HARP data is introduced. We compare the maps to those of the dust continuum emission mapped with the Submillimetre Common-User Bolometer Array (SCUBA; Hatchell et al.) and the positions of starless and protostellar cores (Hatchell et al.). No straightforward correlation is found between the masses of each region derived from the HARP CO and SCUBA data, underlining the care that must be exercised when comparing masses of the same object derived from different tracers. From the ¹³CO/C¹⁸O line ratio the relative abundance of the two species  ([¹³CO]/[C¹⁸O]∼ 7)  and their opacities (typically τ is 0.02–0.22 and 0.15–1.52 for the C¹⁸O and ¹³CO gas, respectively) are calculated. C¹⁸O is optically thin nearly everywhere, increasing in opacity towards star-forming cores but not beyond  τ₁₈∼ 0.9  . Assuming the ¹²CO gas is optically thick, we compute its excitation temperature, T _ex (around 8–30 K), which has little correlation with estimates of the dust temperature.     

Iron needles in supernova remnants?

It has been suggested by Dwek that iron needles could explain the submillimetre emission from the Cas A supernova remnant (SNR) with only a very small total mass. We investigate whether a similar model holds for the Kepler SNR, and find that its emission could indeed be explained by a dust mass of less than  10⁻² M_⊙  , dependent on the axial ratio l / a of the needles – which we constrain to be less than 700. But the implied needle model for Kepler is inconsistent with that suggested for Cas A since either the needles would have to have a resistivity one or two orders of magnitude greater than those in Cas A or the electron density in Kepler's shocked plasma must be 40 times greater than suggested by X-ray observations. An additional problem with the needle model is that the implied thickness of the needles seems to be implausibly small, if the emission properties are calculated under the usual approximations.     

The JCMT Legacy Survey of the Gould Belt: a first look at Orion B with HARP

The Gould Belt Legacy Survey will survey nearby star-forming regions (within 500 pc), using Heterodyne Array Receiver Programme (HARP), Submillimetre Common-User Bolometer Array 2 and Polarimeter 2 on the James Clerk Maxwell Telescope. This paper describes the initial data obtained using HARP to observe ¹²CO, ¹³CO and C¹⁸O   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 the 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 ¹²CO emission, and compare the relative energetics of the high- and low-velocity material in the two clouds. We present a clumpfind analysis of the ¹³CO condensations. The slope of the condensation mass functions, at the high-mass ends, is similar to the slope of the initial mass function.     

CO abundances in a protostellar cloud: freeze-out and desorption in the envelope and outflow of L483

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, ¹²CO   J = 2 → 1, ¹³CO J = 2 → 1  , C¹⁸O   J = 2 → 1  and C¹⁷O   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 James Clerk Maxwell Telescope (JCMT) line observations with the inclusion of freeze-out, a density profile and infall. Our model profiles of ¹²CO and ¹³CO have a large linewidth due to a high-velocity jet. These profiles replicate the process of more abundant material being susceptible to a jet. C¹⁸O and C¹⁷O do not display such a large linewidth as they trace denser quiescent material deep in the cloud.     

Gas in early-type galaxies: cross-fuelling in late-type–early-type pairs?

We present  ¹²CO ( J = 1–0)  and  ¹²CO ( J = 2–1)  observations of eight early-type galaxies, forming part of a sample of interacting galaxies, each consisting of one late- and one early-type system. All of the early-type galaxies observed are undetected in CO to low levels, allowing us to place tight constraints on their molecular gas content. Additionally, we present H  i absorption data for one system. The implications for possible gas transfer from the late- to the early-type galaxy during the interaction are discussed.     

Molecular gas in the Perseus cooling flow galaxy, NGC 1275

The central arcminute of the Perseus cooling flow galaxy, NGC 1275, has been mapped with the JCMT in ¹²CO(2–1) at 21-arcsec resolution, with detections out to at least 36 arcsec (12 kpc). Within the limits of the resolution and coverage, the distribution of gas appears to be roughly east–west, consistent with previous observations of CO, X-ray, Hα and dust emission. The total detected molecular hydrogen mass is ∼ 1.6 × 10¹⁰ M_⊙, using a Galactic conversion factor. The inner central rotating disc is apparent in the data, but the overall distribution is not one of rotation. Rather, the line profiles are bluewards-asymmetric, consistent with previous observations in H  i and [O  iii ]. We suggest that the blueshift may be due to an acquired mean velocity of ∼ 150 km s⁻¹ imparted by the radio jet in the advancing direction. Within the uncertainties of the analysis, the available radio energy appears to be sufficient, and the interpretation is consistent with that of Bo¨hringer et al. for displaced X-ray emission. We have also made the first observations of ¹³CO(2–1) and ¹²CO(3–2) emission from the central 21-arcsec region of NGC 1275 and combined these data with IRAM data supplied by Reuter et al. to form line ratios over equivalent, well-sampled regions. An LVG radiative transfer analysis indicates that the line ratios are not well reproduced by single values of kinetic temperature, molecular hydrogen density and abundance per unit velocity gradient. At least two temperatures are suggested by a simple two-component LVG model, possibly reflecting a temperature gradient in this region.     

The 'Carina Flare' supershell: probing the atomic and molecular ISM in a Galactic chimney

The 'Carina Flare' supershell, GSH 287+04−17, is a molecular supershell originally discovered in  ¹²CO( J = 1–0)  with the NANTEN 4 m telescope. We present the first study of the shell's atomic ISM, using H  i 21-cm line data from the Parkes 64-m telescope Southern Galactic Plane Survey. The data reveal a gently expanding,  ∼230 × 360  pc H  i supershell that shows strong evidence of Galactic Plane blowout, with a break in its main body at   z ∼ 280  pc and a capped high-latitude extension reaching   z ∼ 450  pc. The molecular clouds form comoving parts of the atomic shell, and the morphology of the two phases reflects the supershell's influence on the structure of the ISM. We also report the first discovery of an ionized component of the supershell, in the form of delicate, streamer-like filaments aligned with the proposed direction of blowout. The distance estimate to the shell is re-examined, and we find strong evidence to support the original suggestion that it is located in the Carina Arm at a distance of  2.6 ± 0.4 kpc  . Associated H  i and H₂ masses are estimated as   M _H I≈ 7 ± 3 × 10⁵ M_⊙  and     , and the kinetic energy of the expanding shell as   E _K ∼ 1 × 10⁵¹  erg. We examine the results of analytical and numerical models to estimate a required formation energy of several 10⁵¹ to  ∼10⁵²  erg, and an age of  ∼10⁷ yr  . This age is compatible with molecular cloud formation time-scales, and we briefly consider the viability of a supershell-triggered origin for the molecular component.     

Ram-pressure stripping of the interstellar medium in NGC 4485

This paper describes submm, ¹²CO (J = 2–1) observations of the interacting pair of galaxies NGC 4490 and 4485, and together with high resolution H  i and multifrequency radio continuum data we investigate the evolution of the ISM in this system. We find the following. (i) The smaller member of the pair, NGC 4485, has had the atomic, molecular and dust components of its ISM stripped via ram pressure during its recent passage through the extended H  i distribution of NGC 4490. A bow-shock is identified in the H  i ahead of the stripped gas. (ii) Within the disc of NGC 4490 we find a very low H₂-to-H  i ratio as well as a strong correlation between thermal emission and mass of H₂ suggesting that the star formation rate is limited in this case by the conversion of H  i to H₂. (iii) ¹²CO emission from an H  i and radio-continuum bridge between the two galaxies is detected.     

Comparison of 13CO line and far-infrared continuum emission as a diagnostic of dust and molecular gas physical conditions – I. Motivation and modelling

Determining temperatures in molecular clouds from ratios of CO rotational lines or from ratios of continuum emission in different wavelength bands suffers from reduced temperature sensitivity in the high-temperature limit. In theory, the ratio of far-infrared (FIR), submillimetre or millimetre continuum to that of a ¹³CO (or C¹⁸O) rotational line can place reliable upper limits on the temperature of the dust and molecular gas. Consequently, FIR continuum data from the COBE /Diffuse Infrared Background Experiment (DIRBE) instrument and Nagoya 4-m  ¹³CO  J = 1 → 0  spectral line data were used to plot  240 μm/¹³CO  J = 1 → 0  intensity ratios against 140/240 μm dust colour temperatures, allowing us to constrain the multiparsec-scale physical conditions in the Orion A and B molecular clouds.
The best-fitting models to the Orion clouds consist of two components: a component near the surface of the clouds that is heated primarily by a very large scale (i.e. ∼1 kpc) interstellar radiation field and a component deeper within the clouds. The former has a fixed temperature and the latter has a range of temperatures that vary from one sightline to another. The models require a dust–gas temperature difference of 0 ± 2 K and suggest that 40–50 per cent of the Orion clouds are in the form of dust and gas with temperatures between 3 and 10 K. The implications are discussed in detail in later papers and include stronger dust–gas thermal coupling and higher Galactic-scale molecular gas temperatures than are usually accepted, and an improved explanation for the N (H₂)/ I (CO) conversion factor. It is emphasized that these results are preliminary and require confirmation by independent observations and methods.     

Hyperfine splitting of [Al vi] 3.66 μm and the Al isotopic ratio in NGC 6302

The core of planetary nebula NGC 6302 is filled with high-excitation photoionized gas at low expansion velocities. It represents a unique astrophysical situation in which to search for hyperfine structure (HFS) in coronal emission lines from highly ionized species. HFS is otherwise blended by thermal or velocity broadening. Spectra containing  [Al  vi ] 3.66 μm ³P₂←³P₁  , obtained with Phoenix on Gemini South at resolving powers of up to 75 000, resolve the line into five hyperfine components separated by 20–60 km s⁻¹ as a result of the coupling of the   I = 5/2  nuclear spin of ²⁷Al with the total electronic angular momentum J . The isotope ²⁶Al has a different nuclear spin of   I = 5  , and a different HFS, which allows us to place a 3σ upper limit on the ²⁶Al/²⁷Al abundance ratio of 1/33. We measure the HFS magnetic dipole coupling constants for [Al  vi ], and provide the first estimates of the electric quadrupole HFS coupling constants obtained through astronomical observations of an atomic transition.     

Wind-swept clouds and possible triggered star formation associated with the supernova remnant G357.7+0.3

We present evidence for interaction between the supernova remnant (SNR) G357.7+0.3 and nearby molecular clouds, leading to the formation of wind-swept structures and bright emission rims. These features are not observed at visual wavelengths, but are clearly visible in mid-infrared mapping undertaken using the Spitzer Space Telescope . Analysis of one of these clouds, the bright cometary structure G357.46+0.60, suggests that it contains strong polycyclic aromatic hydrocarbon emission features in the 5.8 and 8.0 μm photometric bands, and that these are highly variable over relatively small spatial scales. The source is also associated with strong variations in electron density; a far-infrared continuum peak associated with dust temperatures of ∼30 K; and has previously been observed in the 1720 MHz maser transition of OH, known to be associated with SNR shock excitation of interstellar clouds. This source also appears to contain a young stellar object (YSO) within the bright rim structure, with a steeply rising spectrum between 1.25 and 24 μm. If the formation of this star has been triggered recently by the SNR, then YSO modelling suggests a stellar mass  ∼5–10 M_⊙  , and luminosity   L _YSO∼10²–2 × 10³ L_⊙  .
Finally, it is noted that a further, conical emission region appears to be associated with the Mira V1139 Sco, and it is suggested that this may represent the case of a Mira outflow interacting with a SNR. If this is the case, however, then the distance to the SNR must be ∼half of that determined from CS   J = 2–1  and 3–2 line radial velocities.     

Extended dust emission in NGC 7465

We present SCUBA 850-μm, JCMT  CO( J =2→1)  , B -band imaging and VLA H  i observations of the NGC 7465/4/3 group of galaxies. The 850-μm emission associated with NGC 7465 extends to at least ∼2 R ₂₅ and is well correlated with the H  i . We investigate a range of possible mechanisms by which dust beyond R ₂₅ may be heated to give the observed extended submillimetre emission. By modelling the dust heating by stars in two extreme geometries, we fail to find any reasonable star formation scenario that is consistent with both the 850-μm and optical data. Furthermore, we do not detect any  CO( J =2→1)  emission coincident with the extended dust and atomic gas as would be expected if significant star formation were occurring. We show that shock-heating of dust via cloud–cloud collisions in the stripped interstellar medium of NGC 7465 could be sufficient to explain the extended 850-μm emission and lack of optical emission in the stripped gas, and suggest that cloud–cloud collisions may be an important dust heating mechanism in gas-rich systems.     

On the possibility of observing H2 emission from primordial molecular cloud kernels

We study the prospects for observing H₂ emission during the assembly of primordial molecular cloud kernels. The primordial molecular cloud cores, which resemble those at the present epoch, can emerge around  1+ z ∼20  according to recent numerical simulations. The kernels form inside the cores, and the first stars will appear inside the kernels. A kernel typically contracts to form one of the first generation stars with an accretion rate that is as large as ∼0.01 M_⊙ yr⁻¹. This occurs owing to the primordial abundances, which result in a kernel temperature of order 1000 K, and the collapsing kernel emits H₂ line radiation at a rate ∼10³⁵ erg s⁻¹. Predominantly   J =5-3   ( v =0)  rotational emission of H₂ is expected. At redshift  1+ z ∼20  , the expected flux is ∼0.01 μJy for a single kernel. While an individual object is not observable by any facilities available in the near future, the expected assembly of primordial star clusters on subgalactic scales can result in fluxes at the sub-mJy level. This is marginally observable with ASTRO-F and ALMA. We also examine the rotational   J =2-0   ( v =0)  and vibrational   δv =1  emission lines. The former may possibly be detectable with ALMA.     

Shocked molecular gas towards the supernova remnant G359.1–0.5 and the Snake

We have found a bar of shocked molecular hydrogen (H₂) towards the OH(1720 MHz) maser located at the projected intersection of supernova remnant (SNR)  G359.1–0.5  and the non-thermal radio filament known as the Snake. The H₂ bar is well aligned with the SNR shell and almost perpendicular to the Snake. The OH(1720 MHz) maser is located inside the sharp western edge of the H₂ emission, which is consistent with the scenario in which the SNR drives a shock into a molecular cloud at that location. The spectral line profiles of ¹²CO, HCO⁺ and CS towards the maser show broad-line absorption, which is absent in the ¹³CO spectra and most probably originates from the pre-shock gas. A density gradient is present across the region and is consistent with the passage of the SNR shock, while the H₂ filament is located at the boundary between the pre-shock and post-shock regions.     

Unresolved emission and ionized gas in the bulge of M31

We study the origin of unresolved X-ray emission from the bulge of M31 based on archival Chandra and XMM–Newton observations. We demonstrate that three different components are present. (i) Broad-band emission from a large number of faint sources – mainly accreting white dwarfs and active binaries, associated with the old stellar population, similar to the Galactic ridge X-ray emission of the Milky Way. The X-ray to K -band luminosity ratios are compatible with those for the Milky Way and for M32; in the 2–10 keV band, the ratio is  (3.6 ± 0.2) × 10²⁷ erg s⁻¹ L⁻¹_⊙  . (ii) Soft emission from ionized gas with a temperature of about ∼300 eV and a mass of  ∼2 × 10⁶ M_⊙  . The gas distribution is significantly extended along the minor axis of the galaxy, suggesting that it may be outflowing in the direction perpendicular to the galactic disc. The mass and energy supply from evolved stars and Type Ia supernovae is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the 'vertical' direction. (iii) Hard extended emission from spiral arms, most likely associated with young stellar objects and young stars located in the star-forming regions. The   L _X/SFR  (star formation rate) ratio equals  ∼9 × 10³⁸ (erg s⁻¹)(M_⊙ yr⁻¹)⁻¹  , which is about ∼1/3 of the high-mass X-ray binary contribution, determined earlier from Chandra observations of other nearby galaxies.     

Cosmic ray generation by quasar remnants: constraints and implications

We report the first detection of CO in the bulge of M31. The ¹²CO (1–0) and (2–1) lines are both detected in the dust complex D395A/393/384, at 1.3 arcmin (∼0.35 kpc) from the centre. From these data and from visual extinction data, we derive a CO luminosity to reddening ratio (and a CO luminosity to H₂ column density ratio) quite similar to that observed in the local Galactic clouds. The (2–1) to (1–0) line intensity ratio points to a CO rotational temperature and a gas kinetic temperature of >10 K. The molecular mass of the complex, inside a 25-arcsec (100 pc) region, is 1.5×10⁴ M_⊙.     

Cassiopeia A: dust factory revealed via submillimetre polarimetry

If Type II supernovae – the evolutionary end points of short-lived, massive stars – produce a significant quantity of dust  (>0.1 M_⊙)  then they can explain the rest-frame far-infrared emission seen in galaxies and quasars in the first Gyr of the Universe. Submillimetre (submm) observations of the Galactic supernova remnant, Cas A, provided the first observational evidence for the formation of significant quantities of dust in Type II supernovae. In this paper, we present new data which show that the submm emission from Cas A is polarized at a level significantly higher than that of its synchrotron emission. The orientation is consistent with that of the magnetic field in Cas A, implying that the polarized submm emission is associated with the remnant. No known mechanism would vary the synchrotron polarization in this way and so we attribute the excess polarized submm flux to cold dust within the remnant, providing fresh evidence that cosmic dust can form rapidly. This is supported by the presence of both polarized and unpolarized dust emission in the north of the remnant where there is no contamination from foreground molecular clouds. The inferred dust polarization fraction is unprecedented  ( f _pol∼ 30 per cent)  which, coupled with the brief time-scale available for grain alignment (<300 yr), suggests that supernova dust differs from that seen in other Galactic sources (where   f _pol= 2−7  per cent) or that a highly efficient grain alignment process must operate in the environment of a supernova remnant.     

The far-infrared–submillimetre spectral energy distribution of high-redshift quasars

We combine photometric observations of high-redshift     quasars, obtained at submillimetre to millimetre wavelengths, to obtain a mean far-infrared (rest-frame) spectral energy distribution (SED) of the thermal emission from dust, parametrized by a single temperature ( T ) and power-law emissivity index ( β ). The best-fitting values are     and     . Our method exploits the redshift spread of this set of quasars, which allows us to sample the SED at a larger number of rest wavelengths than is possible for a single object: the wavelength range extends down to ∼60 μm, and therefore samples the turnover in the greybody curve for these temperatures. This parametrization is of use to any studies that extrapolate from a flux at a single wavelength, for example to infer dust masses and far-infrared luminosities.
We interpret the cool, submillimetre component as arising from dust heated by star formation in the host galaxy of the quasar, although we do not exclude the presence of dust heated directly by the active galactic nucleus (AGN). Applying the mean SED to the data, we derive consistent star formation rates ∼1000 M_⊙ yr⁻¹ and dust masses ∼10⁹ M_⊙, and investigate a simple scheme of AGN and host galaxy co-evolution to account for these quantities. The time-scale for formation of the host galaxy is     , and the luminous quasar phase occurs towards the end of this period, just before the reservoir of cold gas is depleted. Given the youth of the Universe at     (1.6 Gyr), the coexistence of a massive black hole and a luminous starburst at high redshifts is a powerful constraint on models of quasar host galaxy formation.     

Sakurai's Object: characterizing the near-infrared CO ejecta between 2003 and 2007

We present observations of Sakurai's Object obtained at 1–5 μm between 2003 and 2007. By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption features around  4.7 μm  , we determine the excitation conditions in the line-forming region. We find  ¹²C/¹³C = 3.5^+2.0_−1.5  , consistent with CO originating in ejecta processed by the very late thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence of  2.2 × 10⁻⁶≤ M _CO≤ 2.7 × 10⁻⁶ M_⊙  of CO ejecta outside the dust, forming a high-velocity wind of  500 ± 80 km s⁻¹  . We find evidence for significant weakening of the CO band and cooling of the dust around the central star between 2003 and 2005. The gas and dust temperatures are implausibly high for stellar radiation to be the sole contributor.     

Submillimetre continuum images of the NGC 2024 star-forming ridge

We present 450- and 800-μm images, made with the James Clerk Maxwell Telescope, of the NGC 2024 molecular ridge. The seven previously known compact cores, FIR1–7, have been detected, and FIR5 has been resolved into a compact object and an associated extended source to the east. The estimated masses of the dense cores vary between 1.6 and 5.1 M⊙ per 14-arcsec beam, assuming a dust temperature of 30 K and a dust opacity of κ_800 μm = 0.002 m² kg⁻¹. A spectral index map made from the 450- and 800-μm images shows spatial variations, with the spectral index, α ( F _ν ∝ ν^α), being systematically lower towards the dense cores. We interpret this as evidence for a lower value of the frequency dependence of the dust opacity, β, towards the denser cores relative to the surrounding molecular material. This may indicate that grain growth is occurring in the cores, prior to planetesimal formation. By comparing the high-resolution 450-μm image with interferometer maps of the integrated CS(2–1) emission, the previously reported discrepancy between dust continuum emission and molecular line emission is found to be very localized. Depletion and temperature variations are discussed as possible explanations.