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.     

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.     

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.     

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.     

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.     

Discovery of a dense bipolar outflow from a new class 0 protostar in NGC 2068/LBS 17

We report the discovery of high-velocity dense gas from a bipolar outflow source near NGC 2068 in the L1630 giant molecular cloud. CO and HCO⁺ J =3→2 line wings have a bipolar distribution in the vicinity of LBS 17-H with the flow orientated roughly east–west and perpendicular to the elongation of the submillimetre dust continuum emission. The flow is compact (total extent ∼0.2 pc) and contains of the order of 0.1 M_⊙ of swept-up gas. The high-velocity HCO⁺ emission is distributed over a somewhat smaller area <0.1 pc in extent.
A map of C¹⁸O J =2→1 emission traces the LBS 17 core and follows the ambient HCO⁺ emission reasonably well, with the exception of the direction towards LBS 17-H where there is a significant anticorrelation between the C¹⁸O and HCO⁺. A comparison of beam-matched C¹⁸O and dust-derived H₂ column densities suggests that CO is depleted by up to a factor of ∼50 at this position if the temperature is as low as 9 K, although the difference is substantially reduced if the temperature is as high as 20 K. Chemical models of collapsing clouds can account for this discrepancy in terms of different rates of depletion on to dust grains for CO and HCO⁺.
LBS 17-H has a previously known water maser coincident with it but there are no known near-infrared, IRAS or radio continuum sources associated with this object, leading to the conclusion that it is probably very young. A greybody fit to the continuum data gives a luminosity of only 1.7 L_⊙ and a submillimetre-to-bolometric luminosity ratio of 0.1, comfortably satisfying the criteria for classification as a class 0 protostar candidate.     

The Gunn–Peterson effect and the Lyman alpha forest

We show that spatial correlations in a stochastic large-scale velocity field in an otherwise smooth intergalactic medium (homogeneous comoving density) superposed on the general Hubble flow may cause a 'line-like' structure in QSO spectra similar to the population of unsaturated Lyα forest lines which usually are attributed to individual clouds with 10¹¹ ≲ N _{H i}  5 × 10¹³ cm⁻². Therefore there is no clear observational distinction between a diffuse intergalactic medium and discrete intergalactic clouds. It follows that the H  i density in the diffuse intergalactic medium might be substantially underestimated if it is determined from the observed intensity distribution near the apparent continuum in high-resolution spectra of QSOs. Our tentative estimate implies a diffuse neutral hydrogen opacity τ_GP ∼ 0.3 at z  ∼ 3 and a current baryon density Ω_IGM ≃ 0.08, assuming a Hubble constant H ₀ = 70 km s⁻¹ Mpc⁻¹.     

Exciting molecular hydrogen in the central galaxies of cooling flows

The origin of rovibrational H₂ emission in the central galaxies of cooling flow clusters is poorly understood. Here we address this issue using data from our near-infrared spectroscopic survey of 32 of the most line-luminous such systems, presented in the companion paper by Edge et al.
We consider excitation by X-rays from the surrounding intracluster medium (ICM), ultra-violet (UV) radiation from young stars, and shocks. The   v = 1–0  K -band lines with upper levels within  10⁴ K  of the ground state appear to be mostly thermalized (implying gas densities  ≳10⁵ cm⁻³  ), with the excitation temperature typically exceeding 2000 K, as found earlier by Jaffe, Bremer & van der Werf. Together with the lack of strong   v = 2–0  lines in the H -band, this rules out UV radiative fluorescence.
Using the cloudy photoionization code, we deduce that the H₂ lines can originate in a population of dense clouds, exposed to the same hot  ( T ∼ 50 000 K)  stellar continuum as the lower density gas which produces the bulk of the forbidden optical line emission in the Hα-luminous systems. This dense gas may be in the form of self-gravitating clouds deposited directly by the cooling flow, or may instead be produced in the high-pressure zones behind strong shocks. Furthermore, the shocked gas is likely to be gravitationally unstable, so collisions between the larger clouds may lead to the formation of globular clusters.     

On the detection of the linear C5 molecule in the interstellar medium

An upper limit of the column density of the C₅ linear molecule in translucent interstellar clouds is estimated from high-resolution ( R =80 000) and very high signal-to-noise ratio (∼1000) echelle spectra. It is 10¹² cm⁻² per E ( B − V )=1 (two orders of magnitude lower than that of C₂).     

Multifrequency study of the ring nebula SG 13

We investigate the morphology and kinematics of the interstellar medium in the environs of the open cluster Mrk 50, which includes the Wolf–Rayet star WR 157 and a number of early B-type stars. The analysis was performed using radio continuum images at 408 and 1420 MHz, and H  i 21-cm line data taken from the Canadian Galactic Plane Survey, molecular observations of the ¹²CO  ( J = 1 → 0)  line at 115 GHz from the Five College Radio Astronomy Observatory and available mid- and far-infrared (FIR) observations obtained with the Midcourse Space Experiment and IRAS satellites, respectively.
This study allowed the identification of the radio continuum and molecular counterpart of the ring nebula SG 13, while no neutral atomic structure was found to be associated. The nebula is also detected in the images in the mid- and FIR, showing the existence of dust well mixed with the ionized gas. We estimate the main physical parameters of the material linked to the nebula.
The interstellar gas distribution in the environs of Mrk 50 is compatible with a stellar wind bubble created by the mass loss from WR 157.
The distribution of young stellar object candidates in the region shows that the stellar formation activity may be present in the molecular shell that encircles the ring nebula.     

Interstellar CO towards X Persei (HD 24534) – II. Two-component model

The observations made by the Goddard High Resolution Spectrograph (GHRS) aboard the Hubble Space Telescope ( HST ) of molecular CO in absorbing gas towards X Persei are reported. The two-component statistical equilibrium model incorporating radiative excitation of CO by line emission at the same velocity that originates in nearby molecular clouds has been used to reproduce high-resolution GHRS spectra. Earlier analysis indicates that the cloud has a complex structure and at least a two-component model should be used to obtain accurate results. The spectra obtained from the International Ultraviolet Explorer ( IUE ) were used to complement GHRS data and constrain the space of possible solutions. The new oscillator strengths recommended by Eidelsberg et al. for A–X bands have been used. The results show that one of the components may be attributed to the Perseus OB2 molecular cloud, and the other component to an extension of the Taurus dark cloud. The total CO column density N (CO)=(1.0±0.2)×10¹⁶ cm⁻² has been determined. According to the results about 85 per cent of the observed CO belongs to an extension of the Taurus dark cloud. The CO radiation that originates in nearby molecular clouds may be the dominant excitation mechanism of the observed CO. The early results of ¹³CO line analysis indicate a ¹³CO/¹²CO ratio of about 40.     

Accounting for the foreground contribution to the dust emission towards Kepler's supernova remnant

Whether or not supernovae contribute significantly to the overall dust budget is a controversial subject. Submillimetre (sub-mm) observations, sensitive to cold dust, have shown an excess at 450 and 850 μm in young remnants Cassiopeia A (Cas A) and Kepler. Some of the sub-mm emission from Cas A has been shown to be contaminated by unrelated material along the line of sight. In this paper, we explore the emission from material towards Kepler using sub-mm continuum imaging and spectroscopic observations of atomic and molecular gas, via H  i , ¹²CO( J = 2–1) and ¹³CO( J = 2–1). We detect weak CO emission (peak   T *_A  = 0.2–1 K, 1–2 km s⁻¹ full width at half-maximum) from diffuse, optically thin gas at the locations of some of the sub-mm clumps. The contribution to the sub-mm emission from foreground molecular and atomic clouds is negligible. The revised dust mass for Kepler's remnant is  0.1–1.2 M_⊙  , about half of the quoted values in the original study by Morgan et al., but still sufficient to explain the origin of dust at high redshifts.     

A SCUBA survey of compact dark Lynds clouds

We present the first results of a submillimetre continuum survey of Lynds dark clouds. Submillimetre surveys of star-forming regions are an important tool with which to obtain representative samples of the very first phases of star formation. Maps of 24 small clouds were obtained with SCUBA, the bolometer array receiver at the James Clerk Maxwell Telescope, and 19 clouds were detected. The total dark cloud area surveyed was ∼130 arcmin², and a total gas mass of 90 M_⊙ was detected. The dust emission is in general in good agreement with the extinction of optical starlight. The observed clouds contain a newly discovered protostar in L944, and a previously known protostar IRAS 23228+4320 in L1246. Another eight starless cores, either gravitationally unbound or pre-stellar in nature, were also detected. All starless cores and protostars were detected in only seven clouds, and the remaining 17 clouds seem quiescent and do not show any signs of recent star formation activity. The 850-μm images of all detected clouds are presented, as well as 450-μm images of L328, L944, L1014 and L1262. The outflows of the protostars in L944 and L1246 were also discovered and were mapped in ¹²CO J =2→1. The detection of the young protostar in L944, which is not present in the IRAS Point Source Catalog, shows the capacity of submillimetre surveys to detect unknown protostars.     

A disturbed molecular disc in the Circinus galaxy

We have used the Swedish ESO Submillimeter Telescope to observe the molecular gas in the Circinus galaxy using the CO(1 → 0) transition as a tracer. The central region and major axis have been mapped and several other points were also observed. The gas in the galaxy is concentrated towards the nucleus, the peak being coincident with the radio/optical core. The inclination of the molecular galactic disc is more comparable to that of the radio continuum than to that of the large-scale H  i emission. Evidence for an anomalous spur structure pointing radially away from the galactic centre is presented, and may indicate a causal link between it and similar features seen in optical lines and radio continuum. Our data suggest the presence of a central molecular ring or disc with radius 300 ± 50 pc and a rotation velocity of about 200 km s⁻¹ (assuming i  = 73°). The dynamical mass of the nucleus is estimated to be no greater than 3.9 × 10⁹ M_⊙. Assuming that the distribution of gas varies smoothly in the outer regions, we calculate the mass of molecular gas in the galaxy to be at least M _mol = 1.1 × 10⁹ M_⊙, and the star-forming efficiency to be 11 ± 2 L_⊙ M_⊙⁻¹. These results imply that Circinus is undergoing a massive central starburst which may be, at least partially, responsible for its extended minor axis emission seen in several wavebands.     

Molecular hydrogen emission in Cygnus A

We present J , H and K -band spectroscopy of Cygnus A, spanning 1.0–2.4 μm in the rest-frame and hence several rovibrational H₂, H recombination and [Fe  ii ] emission lines. The lines are spatially extended by up to 6 kpc from the nucleus, but their distinct kinematics indicate that the three groups (H, H₂ and [Fe  ii ]) are not wholly produced in the same gas. The broadest line, [Fe  ii ] λ 1.644, exhibits a non-Gaussian profile with a broad base (FWHM≃1040 km s⁻¹), perhaps because of the interaction with the radio source. Extinctions to the line-emitting regions substantially exceed earlier measurements based on optical H recombination lines.
Hard X-rays from the quasar nucleus are likely to dominate the excitation of the H₂ emission. The results of Maloney, Hollenbach & Tielens are thus used to infer the total mass of gas in H₂ v=1–0 S(1)-emitting clouds as a function of radius, for gas densities of 10³ and 10⁵ cm⁻³, and stopping column densities N _H=10²²–10²⁴ cm⁻². Assuming azimuthal symmetry, at least 2.3×10⁸ M_⊙ of such material is present within 5 kpc of the nucleus, if the line-emitting clouds see an unobscured quasar spectrum. Alternatively, if the bulk of the X-ray absorption to the nucleus inferred by Ueno et al. actually arises in a circumnuclear torus, the implied gas mass rises to ∼10¹⁰ M_⊙. The latter plausibly accounts for 10⁹ yr of mass deposition from the cluster cooling flow, for which within this radius.     

An absorption event in the X-ray light curve of NGC 3227

We have monitored the Seyfert galaxy NGC 3227 with the Rossi X-ray Timing Explorer ( RXTE ) since 1999 January. During late 2000 and early 2001 we observed an unusual hardening of the 2–10 keV X-ray spectrum which lasted several months. The spectral hardening was not accompanied by any correlated variation in flux above 8 keV. We therefore interpret the spectral change as transient absorption by a gas cloud of column density 2.6 × 10²³ cm⁻² crossing the line of sight to the X-ray source. A spectrum obtained by XMM–Newton during an early phase of the hard-spectrum event confirms the obscuration model and shows that the absorbing cloud is only weakly ionized. The XMM–Newton spectrum also shows that ∼10 per cent of the X-ray flux is not obscured, but this unabsorbed component is not significantly variable and may be scattered radiation from a large-scale scattering medium. Applying the spectral constraints on the cloud ionization parameter and assuming that the cloud follows a Keplerian orbit, we constrain the location of the cloud to be   R ∼ 10–100  light-days from the central X-ray source, and its density to be   n _H∼ 10⁸ cm⁻³  , implying that we have witnessed the eclipse of the X-ray source by a broad line region cloud.     

Detection of an H92α recombination line in the starburst galaxy NGC 660

We have used the Very Large Array (VLA) to search for the H92α radio recombination line (RRL) in four starburst galaxies. In NGC 660, the line was detected over a 17Å‐8 arcsec² region near its starburst nucleus. The line and continuum emission indicate that the RRL-emitting gas is most likely in the form of a cluster of H ii regions with a small filling factor. Using a simple model we find that the total ionized mass in the nuclear region is in the range 2–8Å‐10⁴ M⊙ and the rate of production of UV photons N _Lyc∼1–3Å‐10⁵³ s⁻¹. The ratio of H92α and Brγ line intensities in NGC 660 indicates that extinction is significant even at λ=2 μm. The velocity field of the ionized gas is consistent with a rotating disc with an average velocity gradient of ∼15 km s⁻¹ arcsec⁻¹. The dynamical mass within the central 500 pc is ∼4Å‐10⁸ M⊙ and may be about ∼6Å‐10⁷ M⊙ within the central 120 pc. No line was detected in the other galaxies (NGC 520, NGC 1614 and NGC 6946) to a 3σ limit of 300 μJy. In the starburst galaxies in which RRLs have been detected, we find that there is a rough correlation between the integrated H92α line flux density and both the total far-infrared flux density and the radio continuum emission from the central region.     

IRAM observations of JVAS/CLASS gravitational lenses

We have searched for molecular absorption lines at millimetre wavelengths in 11 gravitational lens systems discovered in the JVAS/CLASS surveys of flat spectrum radio sources. Spectra of only one source 1030+074 were obtained in the 3-, 2- and 1.3-mm bands at the frequencies corresponding to common molecular transitions of CO and HCO⁺ as continuum emission was not found in any of the other sources. We calculated upper limits to the column density in molecular absorption for 1030+074, using an excitation temperature of 15 K, to be N _CO<6.3×10¹³ cm⁻² and N _HCO+<1.3×10¹¹ cm⁻² , equivalent to hydrogen column density of the order N _H<10¹⁸ cm⁻² , assuming standard molecular abundances. We also present the best upper limits of the continuum at the lower frequency for the other 10 gravitational lenses.     

Radiation pressure acceleration by X-rays in active galactic nuclei

We present calculations of the dynamics of highly ionized gas (HIG) clouds that are confined by external pressure, and are photoionized by AGN continuum. We focus on the gas that is seen, in absorption, in the X-ray spectrum of many AGN and show that such gas can reach hydrostatic equilibrium under various conditions. The principal conclusion is that the clouds can be accelerated to high velocities by the central X-ray source. The dynamical problem can be reduced to the calculation of a single parameter, the average force multiplier, 〈 M 〉. The typical value of 〈 M 〉 is ∼10 suggesting that radiation pressure acceleration by X-rays is efficient for L / L _Edd≳0.1 . The terminal velocity scales with the escape velocity at the base of the flow and can exceed it by a large factor. The typical velocity for a HIG flow that originates at R =10¹⁷ cm in a source with L _x =10⁴⁴ erg s⁻¹ is ∼1000 km s⁻¹, i.e. similar to the velocities observed in several X-ray and UV absorption systems.
Highly ionized AGN clouds are driven mainly by bound–free absorption, and bound–bound processes are less important unless the lines are significantly broadened or the column density is very small. Pressure laws that result in constant or outward decreasing ionization parameters are most effective in accelerating the flow.     

The high-velocity molecular stream in W51B: a ring structure with compact H ii regions

We present ¹³ CO J  = 1 − 0 line observations of the H  ii region complex W51B located in the high-velocity (HV) stream. These observations reveal a filamentary and clumpy structure in the molecular gas. The mean local standard of rest (LSR) velocity ∼ + 65 km s⁻¹ of the molecular gas in this region is greater than the maximum velocities allowed by kinematic Galactic rotation curves. The size and mass of the molecular clouds are ∼ 48 × 17 pc² and ∼ 2.4 × 10⁵ M⊙ respectively. In a position–velocity diagram, molecular gas in the southern part comprises a redshifted ring structure with v _LSR=+ 60 to +73 km s⁻¹. The velocity gradient of this ring is ∼ 0.5 km s⁻¹ pc⁻¹, and the mass is ∼ 6.2 × 10⁴ M⊙. If we assume that the ring is expanding with a uniform velocity, the expansion velocity, radius and kinetic energy are ∼ 7 km s⁻¹, ∼ 13 pc and ∼ 3.0 × 10 ⁴⁹ erg respectively. The kinetic energy and mass spectrum of the ring could be explained by an expanding cylindrical cloud with a centrally condensed mass distribution. The locations of two compact H  ii regions, G49.0−0.3 and G48.9−0.3, coincide with the two molecular clumps in this ring. We discuss star formation, and the mechanism that produced the ring structure.