Volume 521, October 2010
Herschel/HIFI: first science highlights
Article Number L24
Number of page(s) 7
Section Letters
Published online 01 October 2010

Online Material

HIFI has been designed and built by a consortium of institutes and university departments from across Europe, Canada and the United States under the leadership of SRON Netherlands Institute for Space Research, Groningen, The Netherlands and with major contributions from Germany, France and the US. Consortium members are: Canada: CSA, U.Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, MPS; Ireland, NUI Maynooth; Italy: ASI, IFSI-INAF, Osservatorio Astrofisico di Arcetri- INAF; Netherlands: SRON, TUD; Poland: CAMK, CBK; Spain: Observatorio Astronomico Nacional (IGN), Centro de Astrobiologia (CSIC-INTA). Sweden: Chalmers University of Technology - MC2, RSS & GARD; Onsala Space Observatory; Swedish National Space Board, Stockholm University - Stockholm Observatory; Switzerland: ETH Zurich, FHNW; USA: Caltech, JPL, NHSC. The work on star formation at ETH Zurich is partially funded by the Swiss National Science Foundation (grant nr. 200020-113556). This program is made possible thanks to the Swiss HIFI guaranteed time program. This work was supported by the German Deutsche Forschungsgemeinschaft, DFG project number Os 177/1-1. We thank the members of the Herschel key project ``Galactic Cold Cores: A Herschel survey of the source populations revealed by Planck'' lead by M. Juvela (KPOT_mjuvela$\_1$) for providing us with the results of the SPIRE 250 $\mu $m mapping and fruitful discussions. We would also like to acknowledge the use of the JCMT CO(3-2) archival data (PI M. Thompson, M08BU15). A portion of this research was performed at the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space administration. We would like to thank an anonymous referee for constructive comments.

Appendix A: Comparison of line profiles taken with different beams

The different lines discussed here were measured at various telescopes and at different frequencies, so that they all represent somewhat different spatial resolutions. To allow a comparison in terms of a physical interpretation, they have to be translated to a common resolution, so that they stem from the same area on the sky. As a reference, we use $80\hbox {$^{\prime \prime }$ }$, the resolution of the KOSMA observations in the 3-2 transition of the CO isotopes. In principle, all data taken at a finer resolution can be resampled to this beam if the mapped area is large enough.

Unfortunately, most HIFI observations were only single-point observations, not full maps, so that such a convolution is impossible. To derive scaling factors that describe the translation between the measured intensities and the intensity that would be obtained in an $80\hbox {$^{\prime \prime }$ }$ beam, we have to assume a source geometry of the emission. Instead of using any analytic geometry, we use the actually measured distribution of warm dust seen in the sub-mm. By assuming that the spatial distribution of all PDR tracers roughly follows the warm dust, we derive scaling factors for the line intensities at different beam widths, by convolving the sub-mm continuum map with the different beam sizes and picking ratios between the convolved intensities at the measured positions. As a continuum map at IRS1, we used the combination of the SPIRE 250 $\mu $m ( $18\hbox{$^{\prime\prime}$ }$ beam size[*] ) and SCUBA 450 $\mu $m image ( $8\hbox{$^{\prime\prime}$ }$beam size, Holland et al. 1999), because IRS1 is saturated in SPIRE 250 $\mu $m, whereas the observed area of SCUBA 450 $\mu $m is too small to make a convolution map with the beam size of $80\hbox {$^{\prime \prime }$ }$. We regrid the SCUBA 450 $\mu $m image to the same grid as SPIRE 250 $\mu $m, determine the scaling factor between the SCUBA 450 $\mu $m and the SPIRE 250 $\mu $m maps from the overlapping area, and replace the saturated pixels of SPIRE 250 $\mu $m with the scaled SCUBA data. As this combination implies some arbitrariness, we tested four different approaches. We derive the scaling factor in a least squares fit using either (1) all the valid overlapping pixels or (2) only overlapping pixels with SPIRE 250 $\mu $m >100 Jy/beam. When replacing SPIRE pixels, we replace either (1) only saturated pixels, or (2) the full square area ( $5 \times 5$ pixels) containing all saturated pixels. The combination of these provide 4 different beam scaling factors, which are consistent to within 3%. Taking the average of these 4 values, we derive the final factors as shown in Table A.1. A direct convolution to $80\hbox {$^{\prime \prime }$ }$ was possible for the ground-based maps that were observed with a smaller beam, such as the JCMT CO(3-2) map. All the resulting intensities are summarized in Table A.2.

Table A.1:   The scaling factors between the different beam size observations, which should be multiplied to the line intensity to estimate the one in a $80\hbox {$^{\prime \prime }$ }$ beam.

\end{figure} Figure A.1:

a) The spatial distribution of warm dust obtained by combining the SPIRE 250$\mu $m map with the SCUBA 450 $\mu $m map (see text). The HIFI beam is assumed to ``see'' convolved images at beam sizes of b) $42\hbox {$^{\prime \prime }$ }$, c) $55\hbox {$^{\prime \prime }$ }$, and d) $80\hbox {$^{\prime \prime }$ }$. The color scale is common to all the figures. The coordinates are relative to the position of IRS1. Thus, the flux ratios at (0, 0) in these images are the values in Table A.1.

Open with DEXTER

Table A.2:   List of the complementary data.

Appendix B: The clumpy PDR model

To model the far-infrared line emission from S140, we used a superposition of spherical clumps described by the KOSMA-$\tau$ PDR model (Röllig et al. 2006) that represent an ensemble of clumps with a fixed size-spectrum (Cubick et al. 2008). The KOSMA-$\tau$ PDR model simulates a spherical cloud with a radial density profile given by

$\displaystyle n(r) = n_{\rm s} \left\{
c^\delta & {\rm for}~ ...}& {\rm for}~ R/c< r \le R\\
0 & {\rm for}~ r > R .\\
\end{array}\right.$     (B.1)

The constant c determines the dynamic range that is covered by the power-law density decay. The spectrum of PDR clumps is characterized by the clump mass spectrum

\begin{displaymath}{{\rm d} N \over {\rm d} M} = a M^{-\alpha},
\end{displaymath} (B.2)

where the factor a is determined by the total mass of clumps within the beam, $M_{\rm ens}$, and the mass-size relation

\begin{displaymath}M = C R^{\gamma},
\end{displaymath} (B.3)

which implicitly defines the surface density of the individual clumps $n_{\rm s}$. The factor C is determined by the average ensemble density $n_{\rm ens}$. Guided by the clump-decomposition results from Heithausen et al. (1998), we fix the parameters of the spectrum to be $\alpha=1.8$, $\gamma=2.3$, $\delta=1.5$, and c=5. The boundaries of the clump size distribution $M_{\rm min}$ and $M_{\rm max}$ are used as free parameters with the constraint that the maximum clump mass cannot exceed half of the mass of the total ensemble.

Every individual clump is treated as a spherically symmetric configuration illuminated by an isotropic external UV radiation field, specified in terms of the average interstellar radiation field, $\chi_0=2.7\times
10^{-3}$ erg cm-2 s-1 in Draine units, and cosmic rays producing an average ionization rate, $\zeta_{\rm CR}= 5 \times 10^{-17}$ s-1. The internal velocity dispersion of molecules within the clumps is fixed to 1 km s-1. The model computes the stationary chemical and temperature structure by solving the coupled detailed balance of heating, line and continuum cooling, and the chemical network using the UMIST data base of reaction rates (Woodall et al. 2007) expanded by separate entries for the 13C chemistry (see Röllig et al. 2007, for details). The chemical network currently does not include 18O, so that C18O predictions can only be obtained by scaling the 13CO values ignoring fractionation between the two species. As this ignores the different self-shielding of 13CO and C18O, the model results for C18O are less reliable than for the other lines.

In the superposition of clumps, the line emission from the different clumps is simply added assuming that the velocity dispersion between the clumps is large enough, so that they do not shield each other in position-velocity space. This is valid for most species, except for the [O  I] emission, which is opticallyvery thick ( $\tau \ga 100$), so that the lines are much broader than the velocity distribution. Therefore, the model is unable to provide any reliable estimate of the [O  I] intensities. For the continuum extinction of the UV radiation, the situation is different. There, mutual shading of the clouds is relevant, leading to the concept of different clump ensembles that ``see'' different UV fields if the average ensemble extinction exceeds an $A_{\rm V}$ of about unity.

The parameters of the clump mass spectrum imply that most of the mass is actually contained in the largest clumps that also have the maximum column density or $A_{\rm V}$, respectively. The dependence of $A_{\rm V}$ and clump size on clump mass is given by

\begin{displaymath}R = 5.3 \times 10^{18} \left( \frac{M [M_{\odot}] }{n [{\rm cm}^{-3}]} \right)^{\frac{1}{3}}
\end{displaymath} (B.4)


\begin{displaymath}A_{\rm V} =1.6018 \times 10^{21} n[{\rm cm}^{-3}] R\rm [cm].
\end{displaymath} (B.5)

In terms of the total clump surface or the total solid angle of the different clumps, we find about equal contributions from each logarithmic mass bin in the ensemble (Eq. (16) in Cubick et al. 2008). Consequently, we find a non-trivial dependence of the intensity in the different HIFI lines on the clump mass. Figures 3 and 5 from Cubick et al. (2008) show that the [C II] emission is dominated by the largest clumps, while the high-J CO lines are dominated by the smallest clumps. A complex, non-monotonic behavior is observed for the lines from atomic oxygen and mid- to low-J CO isotopes.

As the clump spectrum is purely observationally based, it has no direct relation to a stability criterion. We indeed find that both the most massive clumps are unstable to gravitational collapse and that the smallest clumps will be dispersed on the timescale of a few million years. Therefore, the spectrum can only be considered as a snapshot of interstellar turbulence that reflects the density structure over a timescale of 106-107 years. The assumption of a steady-state chemistry and energy-balance is therefore only applicable if all rates are higher. This holds for the dense clumps in the S140 model fit with densities above 104 cm-3, but for lower densities, an explicitly time-dependent modeling would be required (Viti et al. 2006).

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.