Issue |
A&A
Volume 521, October 2010
Herschel/HIFI: first science highlights
|
|
---|---|---|
Article Number | L25 | |
Number of page(s) | 5 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201015129 | |
Published online | 01 October 2010 |
Herschel/HIFI: first science highlights
LETTER TO THE EDITOR
Gas morphology and energetics at the surface of PDRs: New insights with Herschel
observations of NGC 7023
C. Joblin1,2 -
P. Pilleri1,2 -
J. Montillaud1,2 -
A. Fuente3 -
M. Gerin4 -
O. Berné5 -
V. Ossenkopf6,7 -
J. Le Bourlot8 -
D. Teyssier9 -
J. R. Goicoechea10 -
F. Le Petit8 -
M. Röllig6 -
M. Akyilmaz6 -
A. O. Benz11 -
F. Boulanger12 -
S. Bruderer11 -
C. Dedes11 -
K. France13 -
R. Güsten14 -
A. Harris15 -
T. Klein14 -
C. Kramer16 -
S. D. Lord17 -
P. G. Martin13 -
J. Martin-Pintado10 -
B. Mookerjea18 -
Y. Okada6 -
T. G. Phillips19 -
J.R. Rizzo10 -
R. Simon6 -
J. Stutzki6 -
F. van der Tak7,20 -
H. W. Yorke21 -
E. Steinmetz22 -
C. Jarchow22 -
P. Hartogh22 -
C. E. Honingh6 -
O. Siebertz6 -
E. Caux1,2 -
B. Colin17
1 - Université de Toulouse, UPS, CESR, 9 avenue du colonel Roche, 31028 Toulouse Cedex 4, France
2 -
CNRS, UMR 5187, 31028 Toulouse, France
3 -
Observatorio Astronómico Nacional (OAN), Apdo. 112, 28803 Alcalá de Henares (Madrid), Spain
4 -
LERMA, Observatoire de Paris, 61 Av. de l'Observatoire, 75014 Paris, France
5 -
Leiden Observatory, Universiteit Leiden, PO Box 9513, 2300 RA Leiden, The Netherlands
6 -
I. Physikalisches Institut der Universität
zu Köln, Zülpicher Straße 77, 50937 Köln, Germany
7 -
SRON Netherlands Institute for Space Research, PO Box 800, 9700 AV
Groningen, The Netherlands
8 -
Observatoire de Paris, LUTH and Université Denis Diderot, Place J. Janssen, 92190 Meudon, France
9 -
European Space Astronomy Centre, Urb. Villafranca del Castillo, PO Box 50727, 28080 Madrid, Spain
10 -
Centro de Astrobiología, CSIC-INTA, 28850 Madrid, Spain
11 -
Institute for Astronomy, ETH Zürich, 8093 Zürich, Switzerland
12 -
Institut d'Astrophysique Spatiale, Université Paris-Sud, Bât. 121, 91405 Orsay Cedex, France
13 -
Department of Astronomy and Astrophysics, University of Toronto, 60 St. George Street, Toronto, ON M5S 3H8, Canada
14 -
Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121, Bonn, Germany
15 -
Astronomy Department, University of Maryland, College Park, MD 20742, USA
16 -
Instituto de Radio Astronomía Milimétrica (IRAM), Avenida Divina Pastora 7, Local 20, 18012 Granada, Spain
17 -
IPAC/Caltech, MS 100-22, Pasadena, CA 91125, USA
18 -
Tata Institute of Fundamental Research (TIFR), Homi Bhabha Road, Mumbai 400005, India
19 -
California Institute of Technology, 320-47, Pasadena, CA 91125-4700, USA
20 -
Kapteyn Astronomical Institute, University of Groningen, PO box 800, 9700 AV Groningen, The Netherlands
21 -
Jet Propulsion Laboratory, Caltech, Pasadena, CA 91109, USA
22 -
MPI für Sonnensystemforschung, 37191 Katlenburg-Lindau, Germany
Received 1 June 2010 / Accepted 19 July 2010
Abstract
Context. We investigate the physics and chemistry of the gas
and dust in dense photon-dominated regions (PDRs), along with their
dependence on the illuminating UV field.
Aims. Using Herschel/HIFI observations, we study the gas
energetics in NGC 7023 in relation to the morphology of this
nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V
star and is one of the key targets of Herschel.
Methods. Our approach consists in determining the energetics of
the region by combining the information carried by the mid-IR spectrum
(extinction by classical grains, emission from very small dust
particles) with that of the main gas coolant lines. In this letter, we
discuss more specifically the intensity and line profile of the
158 m (1901 GHz) [C II] line measured by HIFI and provide information on the emitting gas.
Results. We show that both the [C II]
emission and the mid-IR emission from polycyclic aromatic hydrocarbons
(PAHs) arise from the regions located in the transition zone between
atomic and molecular gas. Using the Meudon PDR code and a simple
transfer model, we find good agreement between the calculated and
observed [C II] intensities.
Conclusions. HIFI observations of NGC 7023 provide the
opportunity to constrain the energetics at the surface of PDRs. Future
work will include analysis of the main coolant line [O I] and use of a new PDR model that includes PAH-related species.
Key words: ISM: structure - ISM: kinematics and dynamics - ISM: molecules - submillimeter: ISM
1 Introduction
One main goal of the guarantee time key programme ``Warm and dense
interstellar medium'' (WADI) of the HIFI heterodyne spectrometer (de Graauw et al. 2010)
onboard Herschel (Pilbratt et al. 2010)
is to investigate the physics and chemistry
of the gas and dust in dense photon-dominated regions (PDRs), as well as
their dependence on the illuminating UV field. As part of this programme,
we observed a prototype PDR, NGC 7023.
The region is illuminated by the B2Ve HD 200775
[RA(2000) = 21h01m36.9s ; Dec(2000
],
and has been shaped by the star formation process leading to the formation of a cavity. NGC 7023 has been widely studied
at many wavelengths. It has been shown that this region hosts structures at different gas densities:
cm-3 in the cavity,
104 cm-3 in the PDRs
that are located north-west (NW), south (S) and east, and
105-106 cm-3 in dense filaments and clumps that are observed in the mm (Gerin et al. 1998; Fuente et al. 1996
and references therein) and near-IR (Martini et al. 1997; Lemaire et al. 1996).
NGC 7023 has been mapped by the instruments PACS and SPIRE of Herschel to study the emission of large cold grains (Abergel et al. 2010). We present here some observations of the gas at the surface of this nebula, taking advantage of the very high spectral resolution of HIFI. By combining these observations with previous mid-IR observations, we study the geometry and energetics of the NW and S PDRs.
![]() |
Figure 1:
The NW and S PDRs of NGC 7023 observed by Spitzer-IRAC at 8 |
Open with DEXTER |
Table 1: Summary of the HIFI data.
2 Observations and results
2.1 HIFI observations
The HIFI observations presented here consist in (offset positions are relative to the star; see Fig. 1):
- -
- single pointing using frequency switch mode towards the NW PDR in bands 1a and 3b.
In band 1a, the frequency ranges covered by the wide band spectrometer (WBS) were
[535-539] GHz (LSB) and [547-551] GHz (USB).
In band 3b, covered ranges were [879-883] GHz (LSB) and [891-895] GHz (USB).
The observed position (
,
, called H2 peak) corresponds to the peak intensity of the H2ro-vibrational emission associated to the near-IR filaments (Lemaire et al. 1996);
- -
- on-the-fly (OTF) mapping of the [C II]
1901 GHz emission line in band 7b. The line was covered by both
the WBS and the high-resolution spectrometer (HRS) in USB.
Two cuts were performed: a cut from the star to the NW PDR (star-NW
cut) and a south-north cut (hereafter S-NW cut) covering the
NW PDR, the cavity and the S PDR (
,
, see Fig. 1). The pixel size after regriding is 6.5
.




![]() |
Figure 2:
Left: 13CO, C18O, and HCO+high-J transitions observed with HIFI toward the H2 peak. Right: examples of
[C II] emission profiles at the H2 peak and at the three positions shown by crosses on Fig. 1:
north (
|
Open with DEXTER |
2.2 Gas kinematics
Figure 2 (left panel) shows the 13CO 5-4, 13CO 8-7, C18O 5-4, and HCO+ 6-5 lines observed by HIFI towards the H2 peak. All the lines have a central velocity of about 2.2 km s-1, comparable to previous ground-based
observations in several molecular lines (Fuente et al. 1993).
Figure 2 (right panel) shows the [C II] line profiles at the H2 peak and at different positions along the S-NW
cut.
The line is much broader than molecular lines and its profile shows a complex multi-component structure.
The observations towards the PDRs show that the emission peaks at intermediate velocities (
km s-1),
which have already been observed towards the NW PDR
in several molecular tracers (see, for example, Fuente et al. 1996). There is also a contribution
from higher velocity components (
km s-1), which dominate the emission in the cavity.
2.3 Emission from very small dust particles and C+
![]() |
Figure 3:
Comparison between the [C II] 158 |
Open with DEXTER |
The 158 m (1901 GHz) [C II] and 63
m [O I] lines are the major coolants of the gas at the surface of PDRs
(Hollenbach & Tielens 1999).
In these regions, photoelectric effect dominates the heating, while H2
formation provides a minor contribution. Since the smallest dust
particles, polycyclic aromatic hydrocarbons (PAHs) and very small
grains (VSGs) contribute to a large fraction to this process (Bakes & Tielens 1994; Habart et al. 2001), and these particles emit in the mid-IR most of the energy they absorb in the UV, we expect that the mid-IR and [C II] emissions arise in the same regions.
We used mid-IR spectro-imagery data of the NW PDR of NGC 7023 that were obtained in the 5.5-14.5 m range with the Infrared Spectrograph onboard Spitzer
(Werner et al. 2004). For the S PDR, we used ISOCAM highly-processed data products (Boulanger et al. 1996) from the ISO data archive. To analyse the mid-IR spectra, we followed the method explained in Rapacioli et al. (2005) and Berné et al. (2007),
in which the mid-IR emission is decomposed into three aromatic IR band (AIB) spectra whose carriers are neutral PAHs (PAH0), cationic PAHs (PAH+), and evaporating VSGs.
The fitting procedure was recently improved by including the convolution of the composite spectrum by extinction (Pilleri et al. 2010).
Figure 3 displays the AIB flux,
,
obtained by summing the fluxes of the PAH0, PAH+, and VSG components that were derived from the fit. It shows that the AIB intensity correlates well with the [C II] line intensity, strongly suggesting that both emissions arise from the same regions.
The fit of the mid-IR emission provides two independent tracers of the total gas column density N(H) along the line of sight as explained below.
- (i)
- Owing to the excitation mechanism, the AIB intensity can be written as
where NCAIB is the column density of carbon in the AIB carriers and G0 the UV flux in Habing units (Habing 1968). Assuming that
stays constant at the PDR surface,
can therefore be used as a tracer of N(H) if G0 is known (Pilleri et al. 2010).
- (ii)
- If the column density N(H) is high enough, the
effect of extinction by silicates can be seen on the AIB spectrum. In
the mid-IR fit, the extinction is derived from a simple correction
term, assuming that the emitting and absorbing materials are fully
mixed:
where
is the optical depth in the line of sight. The extinction cross-section per nucleon
is taken from Weingartner & Draine (2001) for RV = 5.5 and N(H) is a free parameter of the fit.
Method (ii) is precise for column densities higher than
cm-2.
Method (i) can probe lower column densities but has two limitations.
The AIB emission needs to be corrected for the variation in the UV
field G0to retrieve the value of N(H). This was done assuming that G0 scales as the inverse squared distance to
the illuminating star HD200775 and a value of G0=2600 at 42
from this star (Pilleri et al. 2010). We used the projected distance as an estimate of the true distance.
This introduces an error that can be especially strong at positions close to the star in the plane of the sky.
Figure 3 shows that the AIB emission stays almost constant at d < 16
,
therefore we used this value as the minimum effective distance of the
NW PDR to the star.
Method (i) also needs to be calibrated since the local emissivity of
the AIB carriers is not known precisely. Our approach was therefore to
derive a calibration factor using the values obtained by
method (ii) around position 42
on the star-NW cut. The same calibration factor was used for all positions along the two cuts. Figure 4 shows that the column densities
that were derived on the two cuts correlate quite well with the [C II] line intensity.
![]() |
Figure 4:
Comparison between the [C II] 158 |
Open with DEXTER |
3 Modelling C+ emission
Table 2:
Summary of the PDR modelling of the [C II] emission for 5 points along the HIFI S-NW cut (
).
The critical density for the [C II] 158 m line is
cm-3 for collisions with H and therefore
the line emissivity depends mainly on temperature for n>ncrit. We selected a few positions on
the HIFI S-NW [C II] cut, three points on the NW PDR and two on the South PDR (cf. Table 2).
The values of G0 were determined as explained in Sect. 2.3,
and we assumed a constant average density with two different values:
cm-3 that is characteristic
of the molecular cloud (Gerin et al. 1998) and
cm-3 that was derived by Rapacioli et al. (2006) in their study of
PAH-related species.
![]() |
Figure 5:
S-NW cut in [C II] emission measured with HIFI at a beam size of 11
|
Open with DEXTER |
We used the 1D Meudon PDR code (Le Petit et al. 2006) to compute the gas temperature T at the cloud surface for all the
selected positions (cf. Table 2). The values of T are used to calculate the C+ level populations.
Line intensities are then derived by integrating along the line-of-sight (perpendicular to model results) and by
assuming uniform excitation conditions.
The thickness of the observed regions leads to an optical depth
,
which implies that transfer effects must be taken into account.
If we assume constant excitation conditions and gas properties along the line of sight, then
and the line intensity can be computed.
The agreement between calculated and observed flux values is very good when using
cm-3.
In the NW PDR, the ratio is 1.0 for NW3 (16) and NW2 (12), and 1.4 for NW1 (-3).
For the S PDR, a value of 2.3 is derived for the two positions, suggesting that systematic effects are causing
the deviation between observed and calculated values of the [C II] flux.
There are several parameters that are not precise in our model but looking at Table 2, it seems the local [C II] emissivity
is mainly affected by the local density and not by the value of G0. For N(H), we assumed the same regions emit in PAHs and [C II] , in agreement with the profiles shown in Fig. 3.
There is also an error for N(H) due to our method (cf. Sect. 2.3), but this error is expected to be the same for both PDRs.
Dividing N(H) by a factor of two leads to lower values of the ratio of the calculated over the observed [CII] flux: 0.7-0.8
for the NW PDR and 1.6-1.7 for the S PDR.
One step further in the model would consist in studying the effect of the grain charge on the photoelectric efficiency (Bakes & Tielens 1994). The relative abundances of PAH+, PAH0, and evaporating VSGs vary significantly over the nebula (Fig. 5). Regions in the cavity appear mainly populated by PAH+ (cf. NW1 (-3) in Table 2).
Since the ionization potential of PAH+ is much higher than that of PAH0 (10 eV compared to
6 eV; Malloci et al. 2007), PAH+ should contribute less to the photoelectric heating than PAH0,
leading to a decrease in the heating rate, hence in the gas cooling. In
its current version, the PDR code uses classical grains with an MRN
distribution (Mathis et al. 1977) and absorption and scattering cross-sections from Laor & Draine (1993).
We have used grains of sizes from 15 Å to 3000 Å with a
dust-to-gas mass ratio of 1%. As a result, the ionization parameter
that quantifies the grain charge (cf. Table 2) does not reflect well the variations of the PAH charge observed in Fig. 5. An upgraded version, in which the PDR code is coupled to the code DUSTEM (Compiègne et al. 2010),
is under development (Gonzalez et al., in prep.) and will allow
for including PAHs. NGC 7023 is clearly a template region that could be
used for these studies.
4 Conclusion
By using HIFI and complementary mid-IR data, we have shown that the [C II] cooling line and the AIB emission arise from the same regions, in the transition zone between atomic and molecular gas. The prototype PDR NGC 7023 was found to be a good object for comparison with PDR models. Further progress on the energetics of this region awaits for the coming [O I] data from the PACS instrument and a PDR model that treats the photophysics of PAHs consistently.
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 Astronímico Nacional (IGN), Centro de Astrobiología (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. This work was supported by the German Deutsche Forschungsgemeinschaft, DFG project number Os 177/1-1. 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.
References
- Abergel, A., Arab, H., Compiègne, M., et al. 2010, A&A, 518, L96 [Google Scholar]
- Bakes, E. L. O., & Tielens, A. G. G. M. 1994, ApJ, 427, 822 [NASA ADS] [CrossRef] [Google Scholar]
- Berné, O., Joblin, C., Deville, Y., et al. 2007, A&A, 469, 575 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Boulanger, F., Reach, W. T., Abergel, A., et al. 1996, A&A, 315, L325 [NASA ADS] [Google Scholar]
- Compiègne, M., Verstraete, L., Jones, A., et al. 2010, ApJ, submitted [Google Scholar]
- de Graauw, Th., Helmich, F. P., Phillips, T. G., et al. 2010, A&A, 518, L6 [Google Scholar]
- Fuente, A., Martin-Pintado, J., Cernicharo, J., & Bachiller, R. 1993, A&A, 276, 473 [NASA ADS] [Google Scholar]
- Fuente, A., Martin-Pintado, J., Neri, R., Rogers, C., & Moriarty-Schieven, G. 1996, A&A, 310, 286 [NASA ADS] [Google Scholar]
- Gerin, M., Phillips, T. G., Keene, J., Betz, A. L., & Boreiko, R. T. 1998, ApJ, 500, 329 [NASA ADS] [CrossRef] [Google Scholar]
- Habart, E., Verstraete, L., Boulanger, F., et al. 2001, A&A, 373, 702 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Habing, H. J. 1968, Bull. Astron. Inst. Netherlands, 19, 421 [Google Scholar]
- Hollenbach, D. J., & Tielens, A. G. G. M. 1999, Rev. Mod. Phys., 71, 173 [Google Scholar]
- Laor, A., & Draine, B. T. 1993, ApJ, 402, 441 [NASA ADS] [CrossRef] [Google Scholar]
- Le Petit, F., Nehmé, C., Le Bourlot, J., & Roueff, E. 2006, ApJS, 164, 506 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lemaire, J. L., Field, D., Gerin, M., et al. 1996, A&A, 308, 895 [NASA ADS] [Google Scholar]
- Malloci, G., Joblin, C., & Mulas, G. 2007, A&A, 462, 627 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martini, P., Sellgren, K., & Hora, J. L. 1997, ApJ, 484, 296 [NASA ADS] [CrossRef] [Google Scholar]
- Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 [NASA ADS] [CrossRef] [Google Scholar]
- Ott, S. 2010, in Astronomical Data Analysis Software and Systems XIX, ed. Y Mizumoto, K. I. Morita, & M. Ohishi, ASP Conf. Ser., in press [Google Scholar]
- Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
- Pilleri, P., Montillaud, J., Berné, O., & Joblin, C. 2010, A&A, submitted [Google Scholar]
- Rapacioli, M., Joblin, C., & Boissel, P. 2005, A&A, 429, 193 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rapacioli, M., Calvo, F., Joblin, C., et al. 2006, A&A, 460, 519 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296 [NASA ADS] [CrossRef] [Google Scholar]
- Werner, M. W., Uchida, K. I., Sellgren, K., et al. 2004, ApJS, 154, 309 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
- ... 7023
- Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.
All Tables
Table 1: Summary of the HIFI data.
Table 2:
Summary of the PDR modelling of the [C II] emission for 5 points along the HIFI S-NW cut (
).
All Figures
![]() |
Figure 1:
The NW and S PDRs of NGC 7023 observed by Spitzer-IRAC at 8 |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Left: 13CO, C18O, and HCO+high-J transitions observed with HIFI toward the H2 peak. Right: examples of
[C II] emission profiles at the H2 peak and at the three positions shown by crosses on Fig. 1:
north (
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Comparison between the [C II] 158 |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Comparison between the [C II] 158 |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
S-NW cut in [C II] emission measured with HIFI at a beam size of 11
|
Open with DEXTER | |
In the text |
Copyright ESO 2010
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.