Testing the circumstellar disc hypothesis: a search for H2 outflow signatures from massive young stellar objects with linearly distributed methanol masers

The results of a survey searching for outflows using near-infrared imaging are presented. Targets were chosen from a compiled list of massive young stellar objects associated with methanol masers in linear distributions. Presently, it is a widely held belief that these methanol masers are found in (and delineate) circumstellar accretion discs around massive stars. If this scenario is correct, one way to test the disc hypothesis is to search for outflows perpendicular to the methanol maser distributions. The main objective of the survey was to obtain wide-field near-infrared images of the sites of linearly distributed methanol masers using a narrow-band 2.12-μm filter. This filter is centred on the  H₂ v = 1–0 S(1)  line; a shock diagnostic that has been shown to successfully trace CO outflows from young stellar objects. 28 sources in total were imaged of which 18 sources display H₂ emission. Of these, only two sources showed emission found to be dominantly perpendicular to the methanol maser distribution. Surprisingly, the H₂ emission in these fields is not distributed randomly, but instead the majority of sources are found to have H₂ emission dominantly parallel to their distribution of methanol masers. These results seriously question the hypothesis that methanol masers exist in circumstellar discs. The possibility that linearly distributed methanol masers are instead directly associated with outflows is discussed.     

ISO observations of molecular hydrogen in the DR 21 bipolar outflow

We demonstrate that a wide range of molecular hydrogen excitation can be observed in protostellar outflows at wavelengths in excess of 5 μm. Cold H₂ in DR 21 is detected through the pure rotational transitions in the ground vibrational level (0–0). Hot H₂ is detected in pure rotational transitions within higher vibrational levels (1–1, 1–2, etc.). Although this emission is relatively weak, we have detected two 1–1 lines in the DR 21 outflow with the ISO SWS instrument. We thus investigate molecular excitation over energy levels corresponding to the temperature range 1015–15 722 K, without the uncertainty introduced by differential extinction when employing near-infrared data.
This gas is thermally excited. We uncover a rather low H₂ excitation in the DR 21 West Peak. The line emission cannot be produced from single C-shocks or J-shocks; a range of shock strengths is required. This suggests that bow shocks and/or bow-generated supersonic turbulence is responsible. We are able to distinguish this shock-excited gas from the fluoresced gas detected in the K band, providing support for the dual-excitation model of Fernandes, Brand & Burton.     

A textbook case of bow shock entrainment in a YSO outflow

Near-infrared images in H₂ line emission and submillimetre maps in CO J  = 3–2 emission illustrate the remarkable association between a molecular bow shock and the redshifted molecular outflow lobe in W75N. The flow lobe fits perfectly into the wake of the bow, as one would expect if the lobe represented swept-up gas. Indeed, these observations strongly support the 'bow shock' entrainment scenario for molecular outflows driven by young stars.   The characteristics of the bow shock and CO outflow lobe are compared with those of numerical simulations of jet-driven flows. These models successfully reproduce the bulge and limb-brightening in the CO outflow, although the model H₂ bow exhibits more structure extending back along the flow axis. We also find that the size of the flow, the high mass fraction in the flow at low outflow velocities (low γ values) and the high CO/H₂ luminosity ratio indicate that the system is evolved. We also predict a correlation, in evolved systems, between outflow age and the CO/H₂ luminosity ratio.     

IRAC photometric analysis and the mid-IR photometric properties of Lyman-break galaxies

We present photometric analysis of deep mid-infrared (mid-IR) observations obtained by Spitzer /IRAC covering the fields Q1422+2309, Q2233+1341, DSF2237a,b, HDFN, SSA22a,b and B20902+34, giving the number counts and the depths for each field. In a sample of 751 Lyman-break galaxies (LBGs) lying in those fields, 443, 448, 137 and 152 are identified at 3.6-, 4.5-, 5.8-, 8.0-μm IRAC bands, respectively, expanding their spectral energy distribution to rest-near-IR and revealing that LBGs display a variety of colours. Their rest-near-IR properties are rather inhomogeneous, ranging from those that are bright in IRAC bands and exhibit  [ R ]−[3.6] > 1.5  colours to those that are faint or not detected at all in IRAC bands with  [ R ]−[3.6] < 1.5  colours and these two groups of LBGs are investigated. We compare the mid-IR colours of the LBGs with the colours of star-forming galaxies and we find that LBGs have colours consistent with star-forming galaxies at   z ∼ 3  . The properties of the LBGs detected in the 8-μm IRAC band (rest-frame K band) are examined separately, showing that they exhibit redder  [ R ]−[3.6]  colours than the rest of the population and that although in general, a multiwavelength study is needed to reach more secure results, IRAC 8-μm band can be used as a diagnostic tool, to separate high z , luminous AGN-dominated objects from normal star-forming galaxies at   z ∼ 3  .     

Interaction between the north-eastern boundary of Sgr A East and giant molecular clouds

We have detected the   v = 1 → 0 S(1) (λ= 2.1218 μm)  and   v = 2 → 1 S(1) (λ= 2.2477 μm)  lines of H₂ in the Galactic Centre, in a  90 × 27 arcsec²  region between the north-eastern boundary of the non-thermal source Sgr A East, and the giant molecular cloud (GMC)  M−0.02 − 0.07  . The detected  H₂ v = 1 → 0  S(1) emission has an intensity of  1.6–21 × 10⁻¹⁸ W m⁻² arcsec⁻²  and is present over most of the region. Along with the high intensity, the large linewidths  (FWHM = 40–70 km s⁻¹)  and the  H₂ v = 2 → 1 S(1)  to   v = 1 → 0 S(1)  line ratios (0.3–0.5) can be best explained by a combination of C-type shocks and fluorescence. The detection of shocked H₂ is clear evidence that Sgr A East is driving material into the surrounding adjacent cool molecular gas. The H₂ emission lines have two velocity components at ∼+50 and  ∼0 km s⁻¹  , which are also present in the NH₃(3, 3) emission mapped by McGary, Coil & Ho. This two-velocity structure can be explained if Sgr A East is driving C-type shocks into both the  GMC M−0.02 − 0.07  and the northern ridge of McGary et al.     

Unlocking the Keyhole: H2 and PAH emission from molecular clumps in the Keyhole Nebula

To better understand the environment surrounding CO emission clumps in the Keyhole Nebula, we have made images of the region in H₂ 1–0 S(1) (2.122-μm) emission and polycyclic aromatic hydrocarbon (PAH) emission at 3.29 μm. Our results show that the H₂ and PAH emission regions are morphologically similar, existing as several clumps, all of which correspond to CO emission clumps and dark optical features. The emission confirms the existence of photodissociation regions (PDRs) on the surface of the clumps. By comparing the velocity range of the CO emission with the optical appearance of the H₂ and PAH emission, we present a model of the Keyhole Nebula whereby the most negative velocity clumps are in front of the ionization region, the clumps at intermediate velocities are in it and those which have the least negative velocities are at the far side. It may be that these clumps, which appear to have been swept up from molecular gas by the stellar winds from η  Car, are now being overrun by the ionization region and forming PDRs on their surfaces. These clumps comprise the last remnants of the ambient molecular cloud around η Car.     

Dust and hydrogen molecules in the extremely metal-poor dwarf galaxy SBS 0335–052

During the early stages of galaxy evolution, the metallicity is generally low and nearby metal-poor star-forming galaxies may provide templates for primordial star formation. In particular, the dust content of such objects is of great importance, because early molecular formation can take place on grains. To gain insight into primeval galaxies at high redshift, we examine the dust content of the nearby extremely low-metallicity galaxy SBS  0335–052  which hosts a very young starburst (≲10⁷ yr). In young galaxies, the dust formation rate in Type II supernovae governs the amount of dust, and by incorporating recent results on dust production in Type II supernovae we model the evolution of dust content. If the star-forming region is compact (≲100 pc), as suggested by observations of SBS  0335–052  , our models consistently explain the quantity of dust, far-infrared luminosity, and dust temperature in this low-metallicity object. We also discuss the H₂ abundance. The compactness of the region is important to H₂ formation, because the optical depth of dust for UV photons becomes large and H₂ dissociation is suppressed. We finally focus on implications for damped Ly α systems.     

The excitation and kinematical properties of H2 and [Fe ii] in the HH 46/47 bipolar outflow

Long-slit spectra of the molecular outflow Herbig–Haro (HH) 46/47 have been taken in the J and K near-infrared bands. The observed H₂ line emission confirms the existence of a bright and extended redshifted counter-jet outflow south-west of HH 46. In contrast with the optical appearance of this object, we show that this outflow seems to be composed of two different emission regions characterized by distinct heliocentric velocities. This implies an acceleration of the counter-jet.
The observed [Fe  ii ] emission suggests an average extinction of 7–9 visual magnitudes for the region associated with the counter-jet.
Through position–velocity diagrams, we show the existence of different morphologies for the H₂ and [Fe  ii ] emission regions in the northern part of the HH 46/47 outflow. We have detected for the first time high-velocity (−250 km s⁻¹) [Fe  ii ] emission in the region bridging HH 46 to HH 47A. The two strong peaks detected can be identified with the optical positions B8 and HH 47B.
The H₂ excitation diagrams for the counter-jet shock suggest an excitation temperature for the gas of T _ex≈2600 K . The lack of emission from the higher energy H₂ lines, such as the 4–3 S(3) transition, suggests a thermal excitation scenario for the origin of the observed emission. Comparison of the H₂ line ratios with various shock models yielded useful constraints about the geometry and type of these shocks. Planar shocks can be ruled out whereas curved or bow shocks (both J- and C-type) can be parametrized to fit our data.     

Spitzer IRAC infrared colours of submillimetre-bright galaxies

High-redshift submillimetre-bright galaxies identified by blank field surveys at millimetre and submillimetre wavelengths appear in the region of the Infra Red Array Camera (IRAC) colour–colour diagrams previously identified as the domain of luminous active galactic nuclei (AGNs). Our analysis using a set of empirical and theoretical dusty starburst spectral energy distribution (SED) models shows that power-law continuum sources associated with hot dust heated by young (≲100 Myr old), extreme starbursts at z > 2 also occupy the same general area as AGNs in the IRAC colour–colour plots. A detailed comparison of the IRAC colours and SEDs demonstrates that the two populations are distinct from each other, with submillimetre-bright galaxies having a systematically flatter IRAC spectrum (≳1 mag bluer in the observed [4.5]–[8.0] colour). Only about 20 per cent of the objects overlap in the colour–colour plots, and this low fraction suggests that submillimetre galaxies powered by a dust-obscured AGN are not common. The red infrared colours of the submillimetre galaxies are distinct from those of the ubiquitous foreground IRAC sources, and we propose a set of infrared colour selection criteria for identifying SMG counterparts that can be used even in the absence of radio or Multiband Imaging Photometer for Spitzer (MIPS) 24 μm data.     

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.     

Molecular hydrogen in damped Lyα systems: clues to interstellar physics at high redshift

In order to interpret H₂ quasar absorption-line observations of damped Lyα systems (DLAs) and subDLAs, we model their H₂ abundance as a function of dust-to-gas ratio, including H₂ self-shielding and dust extinction against dissociating photons. Then, we constrain the physical state of the gas by using H₂ data. Using H₂ excitation data for DLAs with H₂ detections, we derive a gas density  1.5 ≲ log n (cm⁻³) ≲ 2.5  , temperature  1.5 ≲ log T (K) ≲ 3  , and an internal ultraviolet (UV) radiation field (in units of the Galactic value)  0.5 ≲ log χ≲ 1.5  . We then find that the observed relation between the molecular fraction and the dust-to-gas ratio of the sample is naturally explained by the above conditions. However, it is still possible that H₂ deficient DLAs and subDLAs with H₂ fractions less than  ∼10⁻⁶  are in a more diffuse and warmer state. The efficient photodissociation by the internal UV radiation field explains the extremely small H₂ fraction  (≲10⁻⁶)  observed for  κ≲ 1/30  (κ is the dust-to-gas ratio in units of the Galactic value); H₂ self-shielding causes a rapid increase in, and large variations of, H₂ abundance for  κ≳ 1/30  . We finally propose an independent method to estimate the star formation rates of DLAs from H₂ abundances; such rates are then critically compared with those derived from other proposed methods. The implications for the contribution of DLAs to the cosmic star formation history are briefly discussed.     

The emission and kinematic structures of the bipolar planetary nebula M 1-8

We have undertaken echelle spectroscopy and narrow-band line imaging of the bipolar planetary nebula M 1-8. This has permitted us to map the outflow in [N  ii ]λλ 6548+6583 Å, Hα, and in the v = 1–0 S(1) transition of H₂ at λ 2.122 μm. It has also permitted us to acquire high-resolution spectra for [N  ii ]λ 6583 Å, Hα and He  ii λ 6560 Å. Our observations support the results of a previous 2MASS analysis by two of the authors (J. P. Phillips and G. Ramos-Larios), and confirm that there is strong H₂ emission outside of the ionized zone, as well as along the major axis of the outflow. Finally, we have investigated the spatial structure of the outflow in low and high excitation lines, and noted evidence for strong ionization stratification within the envelope of the source. We also note that major axis spectra show asymmetries attributable to outflow along the lobes, oriented at an angle i ∼ 35°–40° to the line of sight. Asymmetries along the minor axis, by contrast, appear to be associated with the central collimating disc, and may be interpretable in terms of asymmetries in disc structure, or rotation at an angular velocity of Ω∼ 1.4 10⁻¹² rad s⁻¹. If the disc arises due to common-envelope evolution, then it seems that angular momentum constraints must be relatively tight, and can only be satisfied given fairly extreme physical assumptions (such as low disc mass, high primary star mass, a low distance to the source and so forth).     

A GLIMPSE of new and possible planetary nebulae at mid-infrared wavelengths

We have undertaken a mid-infrared (MIR) search for new planetary nebulae (PNe) using the Spitzer Space Telescope GLIMPSE Galactic plane survey. This has involved searching extant GLIMPSE data products for morphologically appropriate structures, and investigating sources having IRAS colours similar to those of Galactic PNe. We have found 12 sources which have a high probability of being high-extinction PNe, and which possess MIR and IRAS colours, and shell morphologies similar to those of previously identified Galactic nebulae. Calibrated mapping of these structures and profiles in all four of the IRAC bands (3.6, 4.5, 5.8 and  8.0 μm  ) suggests that many (if not all) of the nebulae possess at least two primary structures: an interior high surface brightness shell, corresponding to what is probably the primary ionized zone, and a much weaker halo extending to very much greater distances from the nucleus. These latter regimes are particularly evident at longer MIR wavelengths (5.8 and  8.0 μm  ), and it is probable that they trace the nebular photodissociative regimes, where emission derives from small-grain continua and/or polycyclic aromatic hydrocarbon molecular bands. This latter behaviour has also been noted in previous analyses of Galactic PNe.     

Further constraints on the evolution of Ks-selected galaxies in the GOODS/CDFS field

We have selected and analysed the properties of a sample of  2905 K_s < 21.5  galaxies in  ∼131 arcmin²  of the Great Observatories Origins Deep Survey (GOODS) Chandra Deep Field South (CDFS), to obtain further constraints on the evolution of K_s -selected galaxies with respect to the results already obtained in previous studies. We made use of the public deep multiwavelength imaging from the optical B through the infrared (IR) 4.5-μm bands, in conjunction with available spectroscopic and COMBO17 data in the CDFS, to construct an optimized redshift catalogue for our galaxy sample. We computed the K_s -band luminosity function and determined that its characteristic magnitude has a substantial brightening and a decreasing total density from   z = 0  to  〈 z 〉= 2.5  . We also analysed the colours and number density evolution of galaxies with different stellar masses. Within our sample, and in contrast to what is observed for less massive systems, the vast majority (∼85–90 per cent) of the most massive  ( M > 2.5 × 10¹¹ M_⊙)  local galaxies appear to be in place before redshift   z ∼ 1  . Around 65–70 per cent of the total assemble between redshifts   z = 1  and 3 and most of them display extremely red colours, suggesting that plausible star formation in these very massive systems should mainly proceed in obscured, short-time-scale bursts. The remaining fraction (up to ∼20 per cent) could be in place at even higher redshifts   z = 3–4  , pushing the first epoch of formation of massive galaxies beyond the limits of current near-IR surveys.     

Global velocity field and bubbles in the blue compact dwarf galaxy Mrk 86

We have studied the velocity field of the blue compact dwarf galaxy Mrk 86 (NGC 2537) using data provided by 14 long-slit optical spectra obtained in 10 different orientations and positions. This kinematical information is complemented with narrow-band ([O  iii ]5007 Å and H α ) and broad-band ( B , V , Gunn r and K ) imaging. The analysis of the galaxy global velocity field suggests that the ionized gas could be distributed in a rotating inclined disc, with projected central angular velocity of Ω=34 km s⁻¹ kpc⁻¹. The comparison between the stellar, H  i and modelled dark matter density profile indicates that the total mass within its optical radius is dominated by the stellar component. Peculiarities observed in its velocity field can be explained by irregularities in the ionized gas distribution or local motions induced by star formation.
Kinematical evidences for two expanding bubbles, Mrk 86–B and Mrk 86–C, are given. They show expanding velocities of 34 and 17 km s⁻¹, H α luminosities of 3×10³⁸ and 1.7×10³⁹ erg s⁻¹, and physical radii of 374 and 120 pc, respectively. The change in the [S  ii ]/H α , [N  ii ]/H α , [O  ii ]/[O  iii ] and [O  iii ]/H β line ratios with the distance to the bubble precursor suggests a diminution in the ionization parameter and, in the case of Mrk 86–B, an enhancement of the shock-excited gas emission. The optical–near-infrared colours of the bubble precursors are characteristic of low‐metallicity star‐forming regions (∼0.2 Z_⊙) with burst strengths of about 1 per cent in mass.     

The mid-infrared colours of Galactic bulge, disc and Magellanic planetary nebulae

We present mid-infrared (MIR) photometry for 367 Galactic disc, bulge and Large Magellanic Cloud (LMC) planetary nebulae (PNe), determined using data acquired with the Spitzer Space Telescope , and through the Legacy Programs GLIMPSE II (Galactic Legacy Infrared Mid-plane Survey Extraordinaire II) and SAGE (Surveying the Agents of the Galaxy's Evolution). This has permitted us to make a comparison between the luminosity functions of bulge and LMC PNe, and between the MIR colours of all three categories of source. It is determined that whilst the  3.6 μm  luminosity functions of the LMC and bulge sources are likely to be closely similar, the [3.6]–[5.8] and [5.8]–[8-0] indices of LMC nebulae are different from those of their disc and bulge counterparts. This may arise because of enhanced  6.2 μm  polycyclic aromatic hydrocarbon emission within the LMC sources, and/or as a result of further, and more radical differences between the spectra of LMC and Galactic PNe. We also determine that the more evolved disc sources listed in the Macquarie/AAO/Strasbourg (MASH) catalogues of Parker et al. and Miszalski et al. have similar colours to those of the less evolved (and higher surface brightness) sources in the catalogue of Acker et al., a result which appears at variance with previous studies of these sources.     

A burst of outflows from the Serpens cloud core: wide-field submillimetre continuum, CO J=2–1 and optical observations

Wide-field mapping of Serpens in submillimetre continuum emission and CO J =2–1 line emission is here complemented by optical imaging in [S  ii ] λλ 6716, 6731 line emission. Analysis of the 450- and 850-μm continuum data shows at least 10 separate sources, along with fainter diffuse background emission and filaments extending to the south and east of the core. These filaments describe 'cavity-like' structures that may have been shaped by the numerous outflows in the region. The dust opacity index, β , derived for the identifiable compact sources is of the order of 1.0±0.2, with dust temperatures in excess of 20 K. This value of β is somewhat lower than for typical class I YSOs; we suggest that the Serpens sources may be 'warm', late class 0 or early class I objects.
With the combined CO and optical data we also examine, on large scales, the outflows driven by the embedded sources in Serpens. In addition to a number of new Herbig–Haro flows (here denoted HH 455–460), a number of high-velocity CO lobes are observed; these extend radially outwards from the cluster of submillimetre sources in the core. A close association between the optical and molecular flows is also identified. The data suggest that many of the submillimetre sources power outflows. Collectively, the outflows traced in CO support the widely recognized correlation between source bolometric luminosity and outflow power, and imply a dynamical age for the whole protostellar cluster of ∼3×10⁴ yr. Notably, this is roughly equal to the proposed duration of the 'class 0' stage in protostellar evolution.     

Excitation and kinematics in H2 bow shocks: near-infrared observations of HH 99 and VLA 1623A (HH 313)

We present a comprehensive near-infrared study of two molecular bow shocks in two protostellar outflows, HH 99 in R Coronae Australis and VLA 1623A (HH 313) in Rho Ophiuchi. New, high-resolution, narrow-band images reveal the well-defined bow shock morphologies of both sources. These are compared with two-dimensional MHD modelling of molecular bows from which we infer flow inclination angles, shock speeds and the magnetic field in the pre-shock gas in each system. With combined echelle spectroscopy and low-resolution K -band spectra we further examine the kinematics and excitation of each source. Bow shock models are used to interpret excitation (CDR) diagrams and estimate the extinction and, in the case of VLA 1623, the ortho–para ratio associated with the observed H₂ population. For the first time, morphology, excitation and kinematics are fitted with a single bow shock model.
Specifically, we find that HH 99 is best fitted by a C-type bow shock model (although a J-type cap is probably responsible for the [Fe  ii ] emission). The bow is flowing away from the observer (at an angle to the line of sight of ∼45°) at a speed of roughly 100 km s⁻¹. VLA 1623A is interpreted in terms of a C-type bow moving towards the observer (at an angle to the line of sight of ∼75°) at a speed of ∼80 km s⁻¹. The magnetic field associated with HH 99 is thought to be orientated parallel to the flow axis; in VLA 1623A the field is probably oblique to the flow axis, since this source is clearly asymmetric in our H₂ images.     

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.