The initial conditions of isolated star formation – VII. Spitzer mapping of pre-stellar cores

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

A SCUBA survey of L1689 – the dog that didn't bark

We present submillimetre data for the L1689 cloud in the ρ Ophiuchi molecular cloud complex. We detect a number of starless and pre-stellar cores and protostellar envelopes. We also detect a number of filaments for the first time in the submillimetre continuum that are parallel both to each other, and to filaments observed in the neighbouring L1688 cloud. These filaments are also seen in the ¹³CO observations of L1689. The filaments contain all of the star-formation activity in the cloud. L1689 lies next to the well-studied L1688 cloud that contains the ρ Oph-A core. L1688 has a much more active star-formation history than L1689 despite their apparent similarity in ¹³CO data. Hence, we label L1689 as the dog that didn't bark. We endeavour to explain this apparent anomaly by comparing the total mass of each cloud that is currently in the form of dense material such as pre-stellar cores. We note firstly that L1688 is more massive than L1689, but we also find that when normalized to the total mass of each cloud, the L1689 cloud has a much lower percentage of mass in dense cores than L1688. We attribute this to the hypothesis of Loren that the star formation in the ρ Ophiuchi complex is being affected and probably dominated by the external influence of the nearby Upper Scorpius OB association and predominantly by σ Sco. L1689 is further from σ Sco and is therefore less active. The influence of σ Sco appears none the less to have created the filaments that we observe in L1689.     

Star-like activity from a very young 'isolated planet'

The discovery of isolated bodies of planetary mass has challenged the paradigm that planets form only as companions to stars. To determine whether 'isolated planets', brown dwarfs and stars can have a common origin, we have made deep submillimetre observations of part of the ρ Oph B star formation region. Spectroscopy of the 9-Jupiter-mass core Oph B-11 has revealed carbon monoxide line wings such as those of a protostar. Moreover, the estimated mass of outflowing gas lies on the force versus core-mass relation for protostars and protobrown dwarfs. This is evidence for a common process that can form any object between planetary and stellar masses in a molecular cloud. In a submillimetre continuum map, six compact cores in ρ Oph B were found to have masses presently below the deuterium-burning limit, extending the core mass function down to  0.01 M_⊙  with the approximate form  d N /d M ∝ M ^−3/2  . If these lowest-mass cores are not transient and can collapse under gravity, then isolated planets should be very common in ρ Oph in the future, as is the case in the Orion star formation region. In fact, the isolated planetary objects that may form from these cores would outnumber the massive planets that have been found as companions to stars.     

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.     

Hyperluminous infrared galaxies

39 galaxies are now known, from follow-up of faint IRAS sources and from submillimetre observations of high-redshift AGN, with far-infrared luminosities >10¹³ L_⊙. 13 of these, which have been found in 60- or 850-μm surveys, form an important unbiased subsample. 12 have been found by comparison of 60-μm surveys with quasar or radio galaxy catalogues, or from infrared surveys with colour selection biased towards AGN, while a further 14 have been found through submillimetre observations of known high-redshift AGN. In this paper I argue, on the basis of detailed modelling of the spectral energy distributions of hyperluminous galaxies with accurate radiative transfer models, and from evidence of high gas mass in several cases, that the bulk of the emission from these galaxies at rest frame wavelengths ≥50 μm is caused by star formation. Even after correction for the effects of lensing, hyperluminous galaxies with emission peaking at rest frame wavelengths ≥50 μm are therefore undergoing star formation at rates >10³ M_⊙ yr⁻¹ and are strong candidates for being primeval galaxies, in the process of a major episode of star formation.     

A shallow though extensive H2 2.122-μm imaging survey of Taurus–Auriga–Perseus – I. NGC 1333, L1455, L1448 and B1

We discuss wide-field near-infrared (near-IR) imaging of the NGC 1333, L1448, L1455 and B1 star-forming regions in Perseus. The observations have been extracted from a much larger narrow-band imaging survey of the Taurus–Auriga–Perseus complex. These H₂ 2.122-μm observations are complemented by broad-band K imaging, mid-IR imaging and photometry from the Spitzer Space Telescope , and published submillimetre CO   J = 3–2  maps of high-velocity molecular outflows. We detect and label 85 H₂ features and associate these with 26 molecular outflows. Three are parsec-scale flows, with a mean flow lobe length exceeding 11.5 arcmin. 37 (44 per cent) of the detected H₂ features are associated with a known Herbig–Haro object, while 72 (46 per cent) of catalogued HH objects are detected in H₂ emission. Embedded Spitzer sources are identified for all but two of the 26 molecular outflows. These candidate outflow sources all have high near-to-mid-IR spectral indices (mean value of  α∼ 1.4  ) as well as red IRAC 3.6–4.5 μm and IRAC/MIPS 4.5–24.0 μm colours: 80 per cent have [3.6]–[4.5] > 1.0 and [4.5]–[24] > 1.5. These criteria – high α and red [4.5]–[24] and [3.6]–[4.5] colours – are powerful discriminants when searching for molecular outflow sources. However, we find no correlation between α and flow length or opening angle, and the outflows appear randomly orientated in each region. The more massive clouds are associated with a greater number of outflows, which suggests that the star formation efficiency is roughly the same in each region.     

Discovery of molecular hydrogen line emission associated with methanol maser emission

We report the discovery of H₂ line emission associated with 6.67-GHz methanol maser emission in massive star-forming regions. In our UNSWIRF/AAT observations, H₂ 1–0 S(1) line emission was found associated with an ultracompact H  ii region IRAS 14567–5846 and isolated methanol maser sites in G318.95–0.20 , IRAS 15278–5620 and IRAS 16076–5134 . Owing to the lack of radio continuum in the latter three sources, we argue that their H₂ emission is shock excited, while it is UV-fluorescently excited in IRAS 14567–5846 . Within the positional uncertainties of 3 arcsec, the maser sites correspond to the location of infrared sources. We suggest that 6.67-GHz methanol maser emission is associated with hot molecular cores, and propose an evolutionary sequence of events for the process of massive star formation.     

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.     

Observation and Study of the CO, CO, HCO and N2H Molecular Lines Towards IRAS 02232+6138

Using the 13.7 m millimeter-wave telescope at the Qinghai Station of Purple Mountain Observatory, we have made observations of ¹³CO, C¹⁸O, HCO⁺ and N₂H⁺ molecular lines towards IRAS 02232+6138. As the excitation density of the probe molecule increases from ¹³CO to HCO⁺, the size of the cloud core associated with IRAS 02232+6138 decreases from 2.40 pc to 0.54 pc, and the virial mass of the cloud core decreases from 2.2 × 10³M to 5.1 × 10²M. A bipolar molecular outflow is found towards IRAS 02232+6138. Using the power function n(r) ∝ r⁻ to fit the spatial density structure of the cloud core, we obtain the power-law index  = 2.3 − 1.2; and we find that, as the probed density increases, the power function becomes more flat. The abundance ratio of ¹³CO to C¹⁸O is 12.4 ± 6.9, comparable with the values 11.8 ± 5.9 for dark clouds and the values 9.0–15.6 for massive cores. The abundance of N₂H⁺ molecules is 3.5 ± 2.5 × 10⁻¹⁰, consistent with the value 1.0 − 5.0 × 10⁻¹⁰ for dark cloud cores and the value 1.2 − 12.8 × 10⁻¹⁰ for massive cores. The abundance of HCO⁺ molecules is 0.9 ± 0.5 × 10⁻⁹, close to the value 1.6 − 2.4 × 10⁻⁹ for massive cores. An increase of HCO⁺ abundance in the outflow region was not found. Combining with the IRAS data, the luminosity-mass ratio of the cloud core is obtained in the range 37–163(L/M). Based on the IRAS luminosity, it is estimated that a main-sequence O7.5 star is probably embedded in the IRAS 02232+6138 cloud core.     

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.     

Photometric distances to nine dark globules

Distances to nine dark globules are determined by a method using optical ( VRI ) and near-infrared (near-IR) ( JHK ) photometry of stars projected towards the field containing the globules. In this method, we compute intrinsic colour indices of stars projected towards the direction of the globule by dereddening the observed colour indices using various trial values of extinction   A _V   and a standard extinction law. These computed intrinsic colour indices for each star are then compared with the intrinsic colour indices of normal main-sequence stars and a spectral type is assigned to the star for which the computed colour indices best match with the standard intrinsic colour indices. Distances ( d ) to the stars are determined using the   A _V   and absolute magnitude  ( M_V )  corresponding to the spectral types thus obtained. A distance versus extinction plot is made and the distance at which   A _V   undergoes a sharp rise is taken to be the distance to the globule. All the clouds studied in this work are in the distance range 160–400 pc. The estimated distances to dark globules LDN 544, LDN 549, LDN 567, LDN 543, LDN 1113, LDN 1031, LDN 1225, LDN 1252 and LDN 1257 are  180 ± 35, 200 ± 40, 180 ± 35, 160 ± 30, 350 ± 70, 200 ± 40, 400 ± 80, 250 ± 50  and 250 ± 50 pc, respectively. Using the distances determined, we have estimated the masses of the globules and the far-IR luminosity of the IRAS sources associated with them. The mass of the clouds studied are in the range  10–200 M_⊙  .     

ULTRACAM observations of the black hole X-ray transient XTE J1118+480 in quiescence

We present 450- and 850-μm maps of R Coronae Australis. We compare the maps with previous surveys of the region, and shed new light on the previously unknown nature of the protostellar sources at the centre of the cloud. We clarify the nature of two millimetre sources previously discovered in lower-resolution data. We identify one new Class 0 protostar, which we label SMM 1B, and we measure the envelope masses of a number of more evolved protostars. We identify two new pre-stellar cores, which we call SMM 1A and SMM 6.     

The SCUBA HAlf Degree Extragalactic Survey – VI. 350-μm mapping of submillimetre galaxies

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 μm has been carried out to map the regions around several 850-μm-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared (IR) luminosities and hence star formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350-μm flux densities for 24 850-μm sources, seven of which are detected at ≥2.5σ within a 10 arcsec search radius of the 850-μm positions. When results from the literature are included the total number of 350-μm flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-IR (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich  ( M _dust≃ 9 × 10⁸ M_⊙)  , luminous  ( L _FIR≃ 2 × 10¹² L_⊙)  star-forming galaxies with intrinsic dust temperatures of ≃35 K and star formation rates of  ≃400 M_⊙ yr⁻¹  . We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350-μm source counts, N ₃₅₀(>25 mJy) ∼ 200–500 deg⁻².     

The SCUBA Local Universe Galaxy Survey – II. 450-μm data: evidence for cold dust in bright IRAS galaxies

This is the second in a series of papers presenting results from the SCUBA Local Universe Galaxy Survey. In our first paper we provided 850-μm flux densities for 104 galaxies selected from the IRAS Bright Galaxy Sample and we found that the 60-, 100-μm ( IRAS ) and 850-μm (SCUBA) fluxes could be adequately fitted by emission from dust at a single temperature. In this paper we present 450-μm data for the galaxies. With the new data, the spectral energy distributions of the galaxies can no longer be fitted with an isothermal dust model – two temperature components are now required. Using our 450-μm data and fluxes from the literature, we find that the 450/850-μm flux ratio for the galaxies is remarkably constant, and this holds from objects in which the star formation rate is similar to our own Galaxy, to ultraluminous infrared galaxies (ULIRGs) such as Arp 220. The only possible explanation for this is if the dust emissivity index for all of the galaxies is ∼2 and the cold dust component has a similar temperature in all galaxies     . The 60-μm luminosities of the galaxies were found to depend on both the dust mass and the relative amount of energy in the warm component, with a tendency for the temperature effects to dominate at the highest L ₆₀. The dust masses estimated using the new temperatures are higher by a factor of ∼2 than those determined previously using a single temperature. This brings the gas-to-dust ratios of the IRAS galaxies into agreement with those of the Milky Way and other spiral galaxies which have been intensively studied in the submm.     

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.     

Infrared dark clouds as precursors to star clusters

Infrared dark clouds (IRDCs) are cold, dense molecular clouds identified as extinction features against the bright mid-infrared Galactic background. Our recent 1.2 mm continuum emission survey of IRDCs reveals many compact (<0.5 pc) and massive (10–2100 M_⊙) cores within them. These prestellar cores hold the key to understanding IRDCs and their role in star formation. Here, we present high angular resolution spectral-line and mm/sub-mm continuum images obtained with the IRAM Plateau de Bure Interferometer and the Sub-Millimeter Array towards three high-mass IRDC cores. The high angular resolution images reveal that two of the cores are resolved into multiple, compact protostellar condensations, while the remaining core contains a single, compact protostellar condensation with a very rich molecular spectrum, indicating that it is a hot molecular core. The derived gas masses for these condensations suggest that each core is forming at least one high-mass protostar, while two of the cores are also forming lower-mass protostars. The close proximity of multiple protostars of disparate mass indicates that these IRDCs are in the earliest evolutionary states in the formation of stellar clusters.     

Complex organic molecules in an early stage of protostellar evolution

We have detected the rotational lines of HCOOCH₃ toward a Class 0 low-mass protostar, NGC1333 IRAS4B, which is reported to be extremely young according to the dynamical age of the molecular outflow (a few 100 yr). This suggests that the complex organic molecules appear from the very early stage of protostellar evolution. On the other hand, the complex organic molecules are not detected in a more evolved protostar, L1527. We have also found a similar trend in a massive star forming region, NGC2264. The HCOOCH₃ emission is almost absent toward IRS1, whereas it is concentrated near MMS3, which is younger than IRS1. In addition, the HCOOCH₃ intensity peak is slightly shifted from the dust emission peak, as is seen in the Orion KL Compact Ridge, giving an important clue to solve its origin.     

Study of star formation in RCW 106 using far-infrared observations

High-resolution far-infrared observations of a large area of the star-forming complex RCW 106 obtained using the TIFR 1-m balloon-borne telescope are presented. Intensity maps have been obtained simultaneously in two bands centred around 150 and 210 μm. Intensity maps have also been obtained in the four IRAS bands using HIRES-processed IRAS data. From the 150- and 210-μm maps, reliable maps of dust temperature and optical depth have been generated. The star formation in this complex has occurred in five linear sub-clumps. Using the map at 210 μm, which has a spatial resolution superior to that of IRAS at 100 μm, 23 sources have been identified. The spectral energy distribution (SED) and luminosity of these sources have been determined using the associations with the IRAS maps. The luminosity distribution of these sources has been obtained. Assuming these embedded sources to be zero-age main-sequence stars and using the mass–luminosity relation for these, the power-law slope of the initial mass function is found to be −1.73±0.5 . This index for this very young complex is about the same as that for more evolved complexes and clusters. Radiation transfer calculations in spherically symmetric geometry have been undertaken to fit the SEDs of 13 sources with fluxes in both the TIFR and the IRAS bands. From this, the r ⁻² density distribution in the envelopes is ruled out. Finally, a correlation is seen between the luminosity of embedded sources and the computed dust masses of the envelopes.     

Imaging and integral-field spectroscopy of shocked H2 around G25.65+1.05

The results of B -band CCD imaging linear polarimetry obtained for stars from the Hipparcos catalogue are used to re-examine the distribution of the local interstellar medium towards the IRAS 100-μm emission void in the Lupus dark clouds. The analysis of the obtained parallax–polarization diagram assigns to the dark cloud Lupus 1 a distance between 130 and 150 pc and assures the existence of a low column density region coincident with the observed infrared void. Moreover, there are clear indications of the existence of absorbing material at distances closer than 60–100 pc, which may be associated with the interface boundary between the Local Bubble and its neighbourhood Loop I superbubble.     

Implications for unified schemes from submillimetre and far-infrared follow-up of radio-selected samples

We extend our previous analysis which used generalized luminosity functions (GLFs) to predict the number of quasars and galaxies in low-radio-frequency-selected samples as a function of redshift, radio luminosity, narrow-emission-line luminosity and type of unified scheme. Our extended analysis incorporates the observed submillimetre (850-μm) flux densities of radio sources, employs a new method which allows us to deal with non-detections, and focuses on the high-luminosity population. First, we conclude that the submillimetre luminosity L ₈₅₀ of low-frequency-selected radio sources is correlated with the bolometric luminosity L _bol of their quasar nuclei via an approximate scaling relation   L ₈₅₀∝ L ^0.7±0.2_bol  . Secondly, we conclude that there is quantitative evidence for a receding-torus-like physical process for the high-luminosity population within a two-population unified scheme for radio sources; this evidence comes from the fact that radio quasars are brighter in both narrow emission lines and submillimetre luminosity than radio galaxies matched in radio luminosity and redshift. Thirdly, we note that the combination of a receding-torus-like scheme and the assumption that the observed submillimetre emission is dominated by quasar-heated dust yields a scaling relation   L ₈₅₀∝ L ^1/2_bol  which is within the errors of that determined here for radio-selected quasars, and consistent with that inferred for radio-quiet quasars.