Issue |
A&A
Volume 507, Number 3, December I 2009
|
|
---|---|---|
Page(s) | 1467 - 1473 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/200911894 | |
Published online | 01 October 2009 |
A&A 507, 1467-1473 (2009)
7 mm continuum observations of ultra compact HII regions
P. Leto1,2 - G. Umana1 - C. Trigilio1 - C. S. Buemi1 - S. Dolei3 - P. Manzitto3 - L. Cerrigone4 - C. Siringo3
1 - INAF - Osservatorio Astrofisico di Catania, via S. Sofia 78, 95123 Catania, Italy
2 -
INAF - Istituto di Radioastronomia, Sezione di Noto, C.P. 161, Noto (SR), Italy
3 -
Università di Catania, Dipartimento di Fisica e Astronomia, via S. Sofia 78, 95123 Catania, Italy
4 -
Max-Planck-Institut für Radioastronomie, Bonn, Germany
Received 20 February 2009 / Accepted 8 September 2009
Abstract
Aims. Ultra compact HII (UCHII) regions are indicators of
high-mass star formation sites and are distributed mainly in the
Galactic plane. They exhibit a broad band spectrum with significant
emission between near-IR and radio wavelengths. We intend to
investigate the possible contribution of the forthcoming ESA Planck
mission to the science of UCHII regions by evaluating the
possibility of detecting UCHIIs that are bright in the radio regime.
Methods. We performed new 7 mm observations of a sample of
UCHII regions. The observations were designed to acquire
high-frequency radio spectra. For each source in our sample, the
free-free radio spectrum has been modeled. Along with far-IR
measurements, our spectra allow us to estimate the flux densities of
the sources in the millimeter and sub-millimeter bands. We extrapolated
and summed the ionized-gas (free-free radio emission) and dust (thermal
emission) contributions in the afore mentioned wavelength ranges. The
possibility of Planck detecting the selected sources can be assessed by
comparing the estimated flux densities to the expected sensitivity in
each Planck channel. To obtain a realistic estimation of the noise
produced by the Galactic emission, the Planck sky model software
package was used.
Results. For each target source, from our new 7 mm data and
other radio measurements from the literature, important physical
parameters such as electron density and their spatial distribution,
source geometry and emission measure were derived. We conclude that, in
the case of the present sample, located close to the Galactic center,
Planck will have a very low detection rate. In contrast, assuming that
our sample is representative of the whole UCHII-region population, we
derive a very high probability of detecting this kind of source with
Planck if located instead close to the anticenter. From the analysis of
the ionized-gas properties, we suggest that the selected sample could
also be contaminated by other kinds of Galactic objects.
Key words: H II regions - radio continuum: ISM - stars: formation - ISM: dust, extinction - stars: circumstellar matter
1 Introduction
The Planck mission of the European Space Agency (ESA) will operate over
a wide range of frequencies, from 30 to 900 GHz
(from 1 cm to about 300 m).
The frequency coverage of the Planck channels was selected to
correspond to the so-called cosmic microwave background (CMB) window,
where temperature fluctuations due to Galactic emission are minimal.
Although Planck was designed for cosmological studies, its
multi-frequency nearly-full-sky maps will have a profound impact on
fundamental physics as well as Galactic and extragalactic astrophysics
and several Galactic projects are indeed included in its core program
(see The Planck Collaboration 2005, for details about the scientific case for the Planck mission).
Table 1: Observations.
In the science of compact Galactic objects, Planck will be sensitive to the millimeter emission from dusty envelopes of young and old stars and to the thermal radio emission from the ionized gas surrounding many classes of stellar objects.
HII regions of size about 0.1 pc, emission measure higher than 107 pc cm-6, electron density higher than 104 cm-3, and ionized mass of around
are
classified as
ultra compact HII regions (UCHIIs) and are sites of high-mass star
formation. In these regions, newly born stars are still embedded in
their natal molecular clouds. The ambient dust thermally reradiates the
entire luminosity of the central star in the infrared (IR), making
UCHIIs among the brightest IR sources in the Galaxy. The radiation
fields of the massive hot stars inside UCHIIs (OB spectral type)
ionize the circumstellar material. The signature of ionization is the
typical free-free flux-density distribution at radio wavelengths (
).
Therefore, we expect the Galactic population of UCHIIs to be a strong foreground contamination source to the Planck experiment.
As part of the preparation of a Planck pre-launch catalog, we started an observing program to measure the 7 mm flux density of samples of Galactic sources with the 32 m INAF-IRA Noto Radiotelescope. In a previous paper (Umana et al. 2008a), we published the flux density measurements of planetary nebulae (PNe).
In this paper, we present new 7 mm single-dish observations of UCHIIs. Main goal was to obtain reliable estimates of the flux densities expected in Planck channels by building and modeling the spectral energy distribution (SED) of our targets. The observing frequency (43 GHz) is within the cosmic microwave background (CMB) window and represents one of the channels of the forthcoming Planck mission. Therefore, in this frequency band, we obtain a direct measurement of the expected flux density with minimal frequency extrapolation.
As added value, the present work provides a sizeable dataset of 7 mm measurements of UCHII regions. These data provide strong constraints on the observed SEDs in the spectral region where free-free emission dominates.
The SED modeling in the radio-millimeter spectral range is a crucial step for the study of the physics of ionized envelopes around UCHII regions. The 7 mm measurements presented in this paper provide the possibility of constraining parameters such as the electron density and the radial profile, the source geometry, and the emission measure.
2 Sample selection
Our sample was selected from the catalog of Giveon et al. (2005a), where a cross-correlation between
the Midcourse Space Experiment (MSX6C) catalog (Egan et al. 2003) and the Multi-Array Galactic Plane Imaging Survey (MAGPIS) (White et al. 2005) was performed. MAGPIS was carried out at 6 cm
and covered Galactic regions with longitude
and latitude
.
A comparison of the two catalogs (radio and infrared) allowed the
authors to find a sample of 687 radio sources with highly-reliable
infrared counterparts. Only 228 of 687 targets were detected
in the new MAGPIS at 20 cm (Helfand et al. 2006).
All sources with both 6 and 20 cm data have red MSX colors and flat or inverted spectral indices, as is typical of thermal radio sources. The average spectral index between 20 and 6 cm is equal to 0.28. Giveon et al. (2005a) classified these sources as compact HII regions on the basis of their radio and infrared properties. For most sources, this was their first classification. Starting from this subsample of 228, we first inspected all the corresponding MAGPIS 6 cm images.
To minimize the source-structure effect on the flux density measurements, we selected only sources
compact enough to be considered point-like at the resolution of the Noto Radiotelescope at 7 mm (HPBW about
)
and rejected those located in high-confusion regions. In our sample, we included single sources of 6 cm angular size
and sources that, although mapped as multiple at 6 cm, showed separations
and
were associated with only one infrared counterpart. Finally, from the
resulting sample we selected only sources with integrated flux
densities higher than 100 mJy at 6 cm and declination higher
than
.
In the conservative hypothesis of an optically thin nebula at 6 cm, a cut-off at both 100 mJy and
guarantees 7 mm detectability and observability with the Noto Radiotelescope.
The expected 7 mm flux density is a lower limit in the case of an optically thick nebula at 6 cm.
All the above selection criteria reduced our sample to a total of 51 UCHIIs.
The names and positions of the selected targets are reported in Cols. 1-3 and 6-8 of Table 1. When observing multiple sources associated with the same infrared counterpart, we pointed the telescope at the MSX coordinates.
3 Observations
The observations reported in this paper were carried out at
different epochs between January 2007 and July 2008, using
the 32 m Radiotelescope in Noto (Italy). The telescope has an
active surface system that accounts for the gravitational deformation
of the primary mirror, providing high quality observations at high
frequencies (see Orfei et al. 2004,2002,
for details of the active surface system mounted on the Noto
radiotelescope). The 7 mm receiver used in our observations is a
cooled superheterodyne receiver of typical zenith system
temperature (
)
about 100 K. The observations were performed with a 400 MHz instantaneous band centered on the sky frequency
MHz, the beam size being
.
The gain ranges from 0.04 to 0.08 K/Jy, depending on the air
mass as a function of the elevation. Our targets were observed using
the On The Fly (OTF) scan technique, which consists of driving the beam
of the telescope across the sources in the RA direction. The
typical scan duration was approximately 20 s, which is short
enough to remain close to the white-noise regime of the radiometer. The
scan length was
,
corresponding to about three telescope beams. To achieve a good
signal-to-noise ratio, each source was observed several times, for a
total integration time of 30 min. Multiple OTF scans were
then added together.
Daily gain curves were obtained and the flux scale fixed using NGC 7027 as a primary calibrator and 3C 286 and DR21 as secondary calibrators. The adopted flux density for NGC 7027 is 5.07 Jy. This value was obtained by correcting that reported in Ott et al. (1994) for the observed evolution of the spectral flux density in NGC 7027 (Zijlstra et al. 2008).
All measurements were corrected for atmospheric extinction. The
correction factor was calculated by performing a series of tipping
scans at elevations between
and
,
covering a range of
from 1.01
to 3. Tipping scans were performed several times (once per hour)
during each observing run. The opacity correction was applied to each
flux measurement after interpolating the results of contiguous tipping
scans. The sky-opacity measurement averaged over all of the observing
sessions was about 0.14. The opacity fluctuation ranges from
about 7% (measured during observing sessions with very stable sky
conditions) to about 30% (observing sessions characterized by
worse weather).
Table 2: Sources parameters.
The telescope pointing was checked before each observing session by
measuring its accuracy for a large sample of pointing calibrators:
SiO maser sources and bright sources at 7 mm. The estimated
pointing accuracy is about
.
Taking into account the telescope pointing accuracy, the error
in the flux density of the primary
gain calibrator, the errors in the opacity measurements, and the
gain-curve error, we estimate the overall calibration accuracy to be
around 10-15.
We detected all targets with a 3
confidence threshold, with a detection rate of 100
.
The results of these observations are reported in Table 1, where the measured 7 mm flux densities (
)
are listed with their errors (Cols. 4, 5 and 9, 10).
4 The SEDs
4.1 The dust contribution
![]() |
Figure 1:
SED template of the dust contribution. The data collected from the
far-IR to the mm range have been normalized to 1000 Jy at 100 |
Open with DEXTER |
To determine the dust contribution for our sample of UCHIIs, we firstly
searched both the literature and online catalogs for data in the range
from 8 m to about 1 mm, where the dust contribution is relevant. As a consequence of the selection criteria, MSX data (Egan et al. 2003) are available for all of the 51 sources of the sample. We found IRAS data (Helou & Walker 1988) for 45 sources, Scuba data for 12 sources at 450
m and for 16 at 850
m (Thompson et al. 2006; Wu et al. 2005; Williams et al. 2004), and Calthech Submillimeter Observatory (CSO) data for 4 sources at 350
m (Motte et al. 2003; Hunter et al. 2000). For 8 sources, we found measurements at 1.2 mm
carried out with the IRAM 30 m Radiotelescope (Beuther et al. 2002; Altenhoff et al. 1994) and SIMBA (Faúndez et al. 2004). Ten sources were also observed at 1.3 mm at the IRTF on Mauna Kea (Chini et al. 1986b,a), at the SEST (Launhardt & Henning 1997), and with the MAMBO bolometer array installed on the IRAM 30 m Radiotelescope (Motte et al. 2003).
To summarize, we have enough measurements to compile SEDs for
34 sources in our sample. From the SEDs, we derive an SED template
by normalizing the 100
m high-quality IRAS data to 1000 Jy and then averaging together all the individual SEDs.
The SED template obtained can be modeled by using two dust components
at different temperatures. We assume that each of the two dust
components emits as a blackbody modified by the frequency-dependent
dust opacity. This gives
,
where
is the Planck function for a dust temperature of
and where
is the emissivity index. The index
strongly depends on the mineralogical composition
of the grains and on their physical shapes. Following the results in Chini & Krügel (1986)
about modeling dust emission from UCHII regions, we find that for
the cold dust component a good fit is achieved for a dust temperature
of 30 K and an emissivity index
.
The hot component has been fitted for a dust temperature of 120 K and by assuming that
(Faúndez et al. 2004).
In Fig. 1, we show the averaged data and overplot their best-fit solution. The deviation of the point at 8 m
could be related to the contribution of either a very hot dust
component or the stellar component, whereas the deviation of the two
measurements at 1.2 and 1.3 mm can be explained as evidence
of an ionized gas contribution.
For 17 sources, we do not have enough measurements at
m
to compile their individual SEDs. To estimate their SED
contribution from cold dust, we assume that all sources in our sample
have the same dust properties. Our SED template is therefore a good
assessment of their emission curves (Fig. 1) and we can obtain an individual SED for each of them by just scaling the template to their MSX 21
m measurement. In the SED template, the flux density at 21
m is found to be a factor of 0.04 of its value at 100
m.
4.2 The ionized-gas contribution
We estimate the contribution from ionized gas to the SEDs of the selected sources by combining our 7 mm measurements with the flux densities at 6 cm in Giveon et al. (2005a) and those at 20 cm, the observations at the lower frequency being available for each source in our sample from the new MAGPIS catalog reported by Giveon et al. (2005b).
In the present work, we need to assess the contribution from free-free
emission in an angular element of size comparable to the Noto
Radiotelescope beam (HPBW about
).
For those sources with multiple peaks in their 6 cm and 20 cm
high-resolution observations, we use as an estimation of their
6 cm and 20 cm flux densities the sum of the integrated
fluxes of each peak.
Literature data at radio wavelengths shorter than 6 cm were also
collected for 18 targets in our sample (Martin-Hernandez et al. 2003; Wood et al. 1988b; Fey et al. 1995; Longmore et al. 2007; Walsh et al. 1998; Kurtz et al. 1994; Wink et al. 1982; Wood et al. 1988a; Gehrz et al. 1995).
The spectral index
between 6 cm (Giveon et al. 2005a) and 7 mm (this paper) was calculated for each source. The values obtained are reported in Table 2. The free-free emission from the ionized gas was derived by modeling all the available radio measurements for every source.
As described by Hoare et al. (1991),
a typical UCHII region can be modeled using spherical symmetry. The
source can be modeled as a central cavity surrounded by a shell of
ionized gas, the newly born high-mass star being located inside the
cavity center. The source geometry is characterized by the linear
diameter and the ratio of the internal to external radii (
/
).
The shell of ionized material is located in the inner region of the
dust cloud
responsible for the infrared emission. In this simple
spherical-geometry model, the spatial density structure of the ionized
gas is described by the power law
.
Adopting the model software developed by Umana et al. (2008b),
the free-free radio spectrum was simulated for each source by assuming
a fixed value of the source linear diameter. In the following, we
consider two cases of the source linear diameter (equal to 2
): 0.1 pc and 0.05 pc. The free-free model spectrum was calculated for different values of the ratio
/
,
the index of the power law that describes the ionized-gas spatial
distribution, and the corresponding average electron density.
The distance to the source was treated as a free parameter in obtaining the best-fit SED to both the observations and simulations. The source angular size derived from the assumed linear diameter and distance needs to be comparable to the angular size measured at 6 cm (Giveon et al. 2005a).
![]() |
Figure 2: Comparison between the derived model parameters for the two analyzed source diameters (0.1 pc and 0.05 pc). Left panel: index of the power law that describes the density gradient inside the source. Right panel: ratio of internal to external source radii. |
Open with DEXTER |
The new 7 mm measurements allow us to characterize well the
free-free radio spectrum of each source in our sample. The
distance-independent parameters are, in fact, univocally defined. The
values of
and
/
for the two sizes do not show any significant difference, as shown in Fig. 2. Using our model we found that the index
covers wide range of model solutions, many sources being reproduced
well by assuming homogeneous electron-density radial distribution,
while others have distributions that are highly inhomogeneous. In all
cases, the geometry of most of the sources could be described by the
ratio
/
below 0.2. The derived electron density
is on the order of 104 cm-3.
Following Terzian & Dickey (1973), it is possible to derive the mean emission measure (EM)
if the observed radio optically thin flux density (S7, in mJy) and the corresponding
source angular size (,
in arcsec) are known,
To derive the mean emission measure from our 7 mm measurements, we assume that the source angular size at 7 mm equals the angular size at 6 cm measured by Giveon et al. (2005a). The derived EM is reported in Table 2. For sources with simulated spectra that are optically thick up to 7 mm, the derived EM is considered to be a lower limit.
Note that the spectral index between 6 cm and 7 mm on its own, is unable to reproduce successfully the emission at 7 mm. In many cases, sources with positive spectral index are characterized by optically thin emission at 7 mm.
Our modeling also allows us to estimate the emission measure. In Fig. 3, the emission measure obtained from the model fit is plotted as a function of the EM calculated from Eq. (1).
![]() |
Figure 3: Comparison of the Emission Measure deduced from the simulation with those calculated from the 7 mm measurements. Arrows refer to those sources optically thick ut to 7 mm. |
Open with DEXTER |
The model used to simulate the radio emission from UCHII regions assumes the highly simplified hypothesis of spherical geometry, whereas radio images at high angular resolution show that these sources are characterized by various morphologies (Kurtz et al. 1994; Wood & Churchwell 1989a; De Pree et al. 2005), i.e., shell-like, bipolar, cometary, irregular, and spherical. In spite of the possible difference between the real and the modeled source geometry, the clear correlation between the values of EM deduced from the model and those calculated from the observations indicates that our model parameters allow us to reproduce well the overall properties of the targets.
Our calculations exhibit a large spread in EM (Fig. 3). HII regions characterized by
pc cm-6 are defined as ultra compact. Sources with EM values
at the edges of the range do not show any significant difference in
their dust emission compared to the other sources of the sample.
No relation was observed between infrared colors and EM. This physical parameter is also randomly distributed as a function of the Galactic longitude.
The infrared colors of the selected sample follow a typical trend for UCHII regions, while in the radio-millimeter regime not all sources have typical UCHII characteristics, as illustrated by their EM values. Our study in the radio range demonstrates that a sample of UCHIIs selected only on the basis of their infrared colors may be contaminated by other classes of Galactic sources. Planetary nebulae or high-mass evolved stars could have infrared colors similar to those of UCHII, but different ionized gas properties.
5 Estimated detectability with Planck
Planck is equipped with a Low Frequency Instrument (LFI), operating at 30, 44, and 70 GHz and a High Frequency Instrument (HFI) operating at 100, 143, 217, 353, 545, and 857 GHz. To evaluate the fraction of UCHIIs in our sample that is detectable with Planck, we need to compare the modeled flux densities of the sources with the expected sensitivity in each observational band.
We derive the modeled flux densities at each Planck channel from the source SEDs by calculating the dust and gas emission as described in Sects. 4.1 and 4.2 and then summing both contributions. In Fig. 4, four of the derived SEDs are shown as an example; the shaded area represents the spectral range covered by the Planck satellite. The dust contribution was obtained by extrapolating the modified-blackbody fit of the cold dust component (dashed lines in Fig. 4), while the free-free contribution was obtained by modeling the data points as described in Sect. 4.2 and extending the derived synthetic spectra to the submillimeter range (dotted line in Fig. 4).
![]() |
Figure 4: Examples of the source SEDs, different symbols refer to different components. Diamonds: hot dust component; squares: cold dust component; triangles: ionized gas. The dust component is deduced from the template SED, the free-free emission of the ionized gas derived from its model. The spectral range covered by Planck is represented by the shaded area. |
Open with DEXTER |
Following Umana et al. (2006), the expected total sensitivity is assumed to be the sum in quadrature of all the sources of confusion noise, namely (i) the nominal instrumental Planck sensitivity per resolution element; (ii) the extragalactic-foreground confusion noise (Toffolatti et al. 1998); (iii) the CMB anisotropy (Bennett et al. 2003a); and (iv) the Galactic-foreground confusion noise.
As pointed out by Wood & Churchwell (1989b), the Galactic distribution of candidate massive-star formation sites is mainly confined to the Galactic plane, and our sample consists indeed of UCHIIs located in the Galactic plane, where we also expect the Galactic emission to be the major source of confusion noise in all Planck channels because of the emission of the dust and gas located in the Galactic disk.
The Galactic noise reported by Umana et al. (2006)
was derived for regions of the sky outside of
the Galactic plane, therefore we cannot use this value for
UCHII regions. For a realistic estimation of the noise related to
the Galactic emission, we simulate the sky emission in the range
between 30 and 900 GHz with the Planck sky model (PSM)
software package (version 1.6.3). PSM is a set of programs
and ancillary data developed by the Planck ESA Working Group 2
that can simulate the sky emission at Planck frequencies and angular
resolution, and reproduces the thermal free-free, synchrotron, and the
dust components, i.e., thermal dust and spinning dust. The free-free component is modeled following Dickinson et al. (2003)
and the spatial structure estimated from the dust-extinction corrected
emission.
The synchrotron component is calculated by assuming a perfect power-law
with different values of the spectral index in different sky
directions, as derived from the 408 MHz (Haslam et al. 1982) and WMAP 23 GHz (Bennett et al. 2003b) maps. The thermal emission from interstellar dust
is assessed following Finkbeiner et al. (1999). The prediction of the rotating grain emission is obtained by multiplying the column densities of the cold neutral medium (Schlegel et al. 1998) with the emissivities calculated using the spinning dust model developed by Draine & Lazarian (1998).
![]() |
Figure 5:
Expected Galactic emission fluctuation for
|
Open with DEXTER |
For the two Planck channels of lower frequency (30 and 44 GHz), the Galactic emission is dominated by the free-free component, whereas for the third channel (70 GHz) the free-free and the thermal dust components are comparable; for the other channels the emission of the thermal dust dominates. The other emission mechanisms provide negligible large-scale contributions.
To obtain an idea of the trend of the Galactic emission
fluctuations as a function of the Galactic longitude, the standard
deviation of the Galactic emission for each PSM channel map was
calculated in boxes of size 3
(
is the Planck beam at the considered frequency), centered on Galactic latitude
and spaced every
in longitude.
We selected the size of the region in which the standard deviation of
the Galactic emission was to be calculated as a compromise between the
necessity both to have a large number of pixels in the region and to
remain well inside the plane of the Galaxy. We show the rms of the
Galactic emission for each channel map at different values of Galactic
longitude for the two Planck instruments in Fig. 5.
The Galactic emission fluctuation strongly depends on the longitude.
The first and fourth Galactic quadrants (center) are characterized by
higher noise with respect to the second and third quadrants
(anticenter), as an obvious consequence of the larger amount of dust in
the direction of the center of the Galaxy.
The correctness of the PSM predictions was tested by comparing them to the observed all-sky WMAP temperature map for Stokes I at 41 GHz (Q band). The WMAP data were retrieved from the NASA Legacy Archive for Microwave Background Data Analysis (LAMBDA), while a sky temperature map at the same frequency and angular resolution was constructed with PSM. We found very good agreement between them, allowing us to use our results on the Galactic emission fluctuation with confidence.
In Fig. 6, we show the percentage of objects from our sample that could be detected above the 3threshold
for each channel of the instruments onboard Planck. Unfortunately, all
the UCHIIs in our sample are located in the first and fourth quadrants
and because of the high Galactic emission fluctuation a very small
fraction of them should be detected by Planck.
![]() |
Figure 6: Expected Planck capability of detecting the sample of UCHII analyzed here. Dashed line: true source Galactic position. Dot dashed line: results for the same population but located in the anticenter region. |
Open with DEXTER |
However, UCHIIs have been detected in any Galactic direction (Wood & Churchwell 1989b). Therefore, assuming that our sample is representative of the entire class of UCHIIs, we can also calculate the detectability of those located in the region of the Galactic anticenter. Figure 6 shows how UCHIIs in the direction of the anticenter should be easily seen by Planck.
Unfortunately, the high-mass star-formation rate is expected to strongly decrease far away from the center of the Galaxy. Bronfman et al. (2000) reported that the OB star-formation sites in the anticenter direction reduced to about 10
those detected in directions close to the center of the Galaxy. The
expected increased Planck capability to detect UCHIIs in the anticenter
direction has to
be compared with the decreasing star formation rate. The true expected
number of detectable UCHIIs must be weighted by the radial distribution
of star formation in the Galaxy.
The UCHII sample analyzed in this paper represents only a subsample of the 228 UCHII regions reported by Giveon et al. (2005a), which are distributed in a region of the Galactic disk
that cover about
of Galactic longitude. Taking into account the estimation of Bronfman et al. (2000),
the expected UCHIIs that Planck could detect for a similar range of
Galactic longitudes in the anticenter direction could be tens. Note
that the Planck channels that have the higher capability of detecting
UCHIIs are the channels at the lower wavelengths (Fig. 6), where the dust emission contribution is dominant. However, the sample selected by Giveon et al. (2005a)
has a bias towards sources above the threshold detectability at the
radio wavelengths. Therefore, it is reasonable to expect that the
estimation reported here is only a lower limit to the true number of
UCHII regions detectable by the Planck satellite,
as there may be many sources with very low radio emission levels, for
actual sensitivities, but with detectable IR dust emission. We
emphasize that this conclusion is supported by many detections of
UCHIIs in the region of the Galactic anticenter obtained with the
balloon experiment Archeops at mm and sub-mm wavelengths (Desert et al. 2008).
In many cases UCHIIs have been observed on the borders of larger and fainter HII regions. In higher wavelength channels, the Planck beam may be large enough to detect emission from these associated normal HII regions. Considering their contribution, the Planck capability to detect UCHIIs here reported could be an underestimate.
We conclude that we expect the Planck mission to provide an important contribution to studies of UCHIIs, as it will provide the scientific community with measurements in a still largely-unexplored spectral region.
6 Summary
We have presented 7 mm (43 GHz) observations of a sample of radio-bright ultra-compact HII regions carried out with the INAF-IRA Noto Radiotelescope. The results of this survey provide a sizeable database of new millimeter measurements that allow us to characterize the high-frequency free-free contribution from the ionized-gas components of the selected targets. The free-free spectrum of each source were derived by modeling all the available radio and millimeter observations and the new 7 mm data.
We used our 7 mm measurements along with infrared and sub-millimeter data to compile SEDs from the radio to the far-IR and to estimate the expected flux densities in the spectral region between the cm and the sub-mm ranges.
The study of properties of Galactic UCHIIs is one of the secondary science objectives of the Planck core program (The Planck Collaboration 2005) therefore, it appears to be important to test the capability of Planck to detect UCHII regions. We, thus, estimated the fluctuation of the Galactic emission as a function of the Galactic longitude by using the Planck Sky Model software package. In this paper, we note that we report for the first time a quantitative estimation of the emission fluctuation of the Galactic disk as seen by Planck. As expected, the Galactic emission fluctuation is higher near the center of the Galaxy. When comparing the modeled radio/sub-mm flux densities to the expected total sensitivity of the Planck satellite, we evaluate that few sources are likely to be detected by Planck in the direction of the Galactic center, where the sources of our sample are located. If we assume that the sample analyzed in the present paper is representative of the whole Galactic UCHII population, we infer a higher rate of detectability by Planck in the direction of the Galactic anticenter, where the Galactic emission fluctuation is significantly lower.
Therefore, even though designed for cosmological studies, Planck should also contribute to the UCHII science. The rich source database provided by Planck will allow follow-up observations aimed at the investigating morphological details with future instruments working in the same spectral region, such as ALMA.
The source parameter emission measures were derived both directly from the 7 mm measurements and, independently, from a model of the ionized-gas continuum emission. The derived EM values for each sources are clearly correlated, supporting the goodness of the model fit. The analysis of this parameter leads us to conclude that our sample could be contaminated by other types of Galactic sources, although all of the selected targets can be classified as UCHII regions because of their infrared colors.
Based on observations with the Noto Radiotelescope operated by INAF-Istituto di Radioastronomia. This work has been partially supported by ASI through contract Planck LFI Activity of Phase E2. The authors acknowledge the use of the Planck Sky Model, developed by the Component Separation Working Group (WG2) of the Planck Collaboration. We thank the referee for his/her constructive criticism which enabled us to improve this paper.
References
- Altenhoff, W. J., Thum, C., & Wendker, H. J. 1994, A&A, 281, 161 [NASA ADS]
- Bennett, C. L., Halpern, M., Hinshaw, G., et al. 2003a, ApJS, 148, 1 [NASA ADS] [CrossRef]
- Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2003b, ApJS, 148, 97 [NASA ADS] [CrossRef]
- Beuther, H., Schilke, P., Menten, K. M., et al. 2002, ApJ, 566, 945 [NASA ADS] [CrossRef]
- Bronfman, L., Casassus, S., May, J., & Nyman, L.-Å. 2000, A&A, 358, 521 [NASA ADS]
- Chini, R., & Krügel, E. 1986, A&A, 167, 315 [NASA ADS]
- Chini, R., Kreysa, E., Mezger, P. G., & Gemund, H.-P. 1986a, A&A, 154, L8 [NASA ADS]
- Chini, R., Kreysa, E., Mezger, P. G., & Gemund, H.-P. 1986b, A&A, 157, L1 [NASA ADS]
- Delabrouille, J., et al. 2009, in preparation
- Desert, F.-X., Macias-Perez, J. F., Mayet, F., et al. 2008, A&A, 481, 411 [NASA ADS] [CrossRef] [EDP Sciences]
- De Pree, C. G., Wilner, D. J., Deblasio, J., Mercer, A. J., & Davis, L. E. 2005, ApJ, 624, L101 [NASA ADS] [CrossRef]
- Dickinson, C., Davies, R. D., & Davis, R. J. 2003, MNRAS, 341, 369 [NASA ADS] [CrossRef]
- Draine, B. T., & Lazarian, A. 1998, ApJ, 508, 157 [NASA ADS] [CrossRef]
- Egan, M. P., Price, S. D., Kraemer, K. E., et al. 2003, VizieR On-line Data Catalog, originally published in: Air Force Research Lab. Technical Rep. AFRL-VS-TR-2003-1589
- Faúndez, S., Bronfman, L., Garay, G., et al. 2004, A&A, 426, 97 [NASA ADS] [CrossRef] [EDP Sciences]
- Fey, A. L., Gaume, R. A., Claussen, M. J., & Vrba, F. J. 1995, ApJ, 453, 308 [NASA ADS] [CrossRef]
- Finkbeiner, D. P., Davis, M., & Schlegel, D. J. 1999, ApJ, 524, 867 [NASA ADS] [CrossRef]
- Gehrz, R. D., Hayward, T. L., Houck, J. R., et al. 1995, ApJ, 439, 417 [NASA ADS] [CrossRef]
- Giveon, U., Becker, R. H., Helfand, D. J., & White, R. L. 2005a, AJ, 129, 348 [NASA ADS] [CrossRef]
- Giveon, U., Becker, R. H., Helfand, D. J., & White, R. L. 2005b, AJ, 130, 156 [NASA ADS] [CrossRef]
- Haslam, C. G. T., Salter, C. J., Stoffel, H., & Wilson, W. E. 1982, A&AS, 47, 1 [NASA ADS]
- Helfand, D. J., Becker, R. H., White, R. L., Fallon, A., & Tuttle, S. 2006, AJ, 131, 2525 [NASA ADS] [CrossRef]
- Helou, G., & Walker, D. W. 1988, IRAS Catalogs and Atlases Explanatory Supplement, NASA RP-1190, 1
- Hoare, M. G., Roche, P. F., & Glencross, W. M. 1991, MNRAS, 251, 584 [NASA ADS]
- Hunter, T. R., Churchwell, E., Watson, C., et al. 2000, AJ, 119, 2711 [NASA ADS] [CrossRef]
- Kurtz, S., Churchwell, E., & Wood, D. O. S. 1994, ApJS, 91, 659 [NASA ADS] [CrossRef]
- Launhardt, R., & Henning, T. 1997, A&A, 326, 329 [NASA ADS]
- Leach, S. M., Cardoso, J.-F., Baccigalupi, C., et al. 2008, A&A, 491, 597 [NASA ADS] [CrossRef] [EDP Sciences]
- Longmore, S. N., Burton, M. G., Barnes, P. J., et al. 2007, MNRAS, 379, 535 [NASA ADS] [CrossRef]
- Martin-Hernandez, N. L., Van Der Hulst, J. M., & Tielens, A. G. G. M. 2003, A&A, 407, 957 [NASA ADS] [CrossRef] [EDP Sciences]
- Motte, F., Schilke, P., & Lis, D. C. 2003, ApJ, 582, 277 [NASA ADS] [CrossRef]
- Orfei, A., Morsiani, M., Zacchiroli, G., et al. 2002, Proceedings of the 6th European VLBI Network Symposium, 13
- Orfei, A., Morsiani, M., Zacchiroli, G., et al. 2004, SPIE, 5495, 116 [NASA ADS]
- Ott, M., Witzel, A., Quirrenbach, A., et al. 1994, A&A, 284, 3310
- Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [NASA ADS] [CrossRef]
- Terzian, Y., & Dickey, J. 1973, AJ, 78, 875 [NASA ADS] [CrossRef]
- The Planck Collaboration 2005, ESA-SCI [arXiv:astro-ph/0604069]
- Thompson, M. A., Hatchell, J., Walsh, A. J., MacDonald, G. H., & Millar, T. J. 2006, A&A, 453, 1003 [NASA ADS] [CrossRef] [EDP Sciences]
- Toffolatti, L., Argueso Gomez, F., de Zotti, G., et al. 1998, MNRAS, 297, 117 [NASA ADS] [CrossRef]
- Umana, G., Burigana, C., & Trigilio, C. 2006, MSAIS, 9, 279 [NASA ADS]
- Umana, G., Leto, P., Trigilio, C., et al. 2008a, A&A, 482, 529 [NASA ADS] [CrossRef] [EDP Sciences]
- Umana, G., Trigilio, C., Cerrigone, L., Buemi, C. S., & Leto, P. 2008b, MNRAS, 386, 1404 [NASA ADS] [CrossRef]
- Walsh, A. J., Burton, M. G., Hyland, A. R., & Robinson, G. 1998, MNRAS, 301, 640 [NASA ADS] [CrossRef]
- White, R. L., Becker, R. H., & Helfand, D. J. 2005, AJ, 130, 586 [NASA ADS] [CrossRef]
- Williams, S. J., Fuller, G. A., & Sridharan, T. K. 2004, A&A, 417, 115 [NASA ADS] [CrossRef] [EDP Sciences]
- Wink, J. E., Altenhoff, W. J., & Mezger, P. G. 1982, A&A, 108, 227 [NASA ADS]
- Wood, D. O. S., & Churchwell, E. 1989a, ApJS, 69, 831 [NASA ADS] [CrossRef]
- Wood, D. O. S., & Churchwell, E. 1989b, ApJ, 340, 265 [NASA ADS] [CrossRef]
- Wood, D. O. S., Churchwell, E., & Salter, C. J. 1988a, ApJ, 325, 694 [NASA ADS] [CrossRef]
- Wood, D. O. S., Handa, T., Fukui, Y., et al. 1988b, ApJ, 326, 884 [NASA ADS] [CrossRef]
- Wu, Y., Zhu, M., Wei, Y., et al. 2005, ApJ, 628, L57 [NASA ADS] [CrossRef]
- Zijlstra, A. A., van Hoof, P. A. M., & Pearly, R. A. 2008, ApJ, 681, 1296 [NASA ADS] [CrossRef]
Footnotes
- ... dust
- See Delabrouille et al. (2009) and Leach et al. (2008) for applications and capabilities of the PSM package.
All Tables
Table 1: Observations.
Table 2: Sources parameters.
All Figures
![]() |
Figure 1:
SED template of the dust contribution. The data collected from the
far-IR to the mm range have been normalized to 1000 Jy at 100 |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Comparison between the derived model parameters for the two analyzed source diameters (0.1 pc and 0.05 pc). Left panel: index of the power law that describes the density gradient inside the source. Right panel: ratio of internal to external source radii. |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Comparison of the Emission Measure deduced from the simulation with those calculated from the 7 mm measurements. Arrows refer to those sources optically thick ut to 7 mm. |
Open with DEXTER | |
In the text |
![]() |
Figure 4: Examples of the source SEDs, different symbols refer to different components. Diamonds: hot dust component; squares: cold dust component; triangles: ionized gas. The dust component is deduced from the template SED, the free-free emission of the ionized gas derived from its model. The spectral range covered by Planck is represented by the shaded area. |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Expected Galactic emission fluctuation for
|
Open with DEXTER | |
In the text |
![]() |
Figure 6: Expected Planck capability of detecting the sample of UCHII analyzed here. Dashed line: true source Galactic position. Dot dashed line: results for the same population but located in the anticenter region. |
Open with DEXTER | |
In the text |
Copyright ESO 2009
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.