A&A 467, 93-106 (2007)
DOI: 10.1051/0004-6361:20066692
M. Boquien1,2 - P.-A. Duc1 - J. Braine3 - E. Brinks4 - U. Lisenfeld5 - V. Charmandaris6,7,8
1 - Laboratoire AIM,
CEA, CNRS et Université Paris Diderot
2 - CEA-Saclay, DSM/DAPNIA/Service d'Astrophysique, 91191 Gif-sur-Yvette Cedex, France
3 - Observatoire de Bordeaux, UMR 5804, CNRS/INSU, BP 89, 33270 Floirac, France
4 - Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield AL10 9AB, UK
5 - Dept. de Física Teórica y del Cosmos, Universidad de Granada, Granada, Spain
6 - Department of Physics, University of Crete, 71003 Heraklion, Greece
7 - IESL/Foundation for Research and Technology, Hellas, 71110 Heraklion, Greece
8 - Chercheur Associé, Observatoire de Paris, 75014 Paris, France
Received 3 November 2006 / Accepted 18 February 2007
Abstract
Context. Star formation (SF) takes place in unusual places such as way out in the intergalactic medium out of material expelled from parent galaxies.
Aims. Whether SF proceeds in this specific environment in a similar way than in galactic disks is the question we wish to answer. Particularly, we address the reliability of ultraviolet, H
and mid-infrared as tracers of SF in the intergalactic medium.
Methods. We have carried out a multiwavelength analysis of the interacting system NGC 5291, which is remarkable for its extended HI ring hosting numerous intergalactic HII regions. We combined new ultraviolet (GALEX) observations with archival H,
8
m (Spitzer Space Telescope) and HI (VLA B-array) images of the system.
Results. We have found that the morphology of the star forming regions, as traced by the ultraviolet, H,
and 8
m emission is similar. The 8.0
m infrared emission, normalised to emission from dust at 4.5
m, which is known to be dominated by PAH bands, is comparable to the integrated emission of dwarf galaxies of the same metallicity and to the emission of individual HII regions in spirals. The 8.0
m in the intergalactic environment is therefore an estimator of the star formation rate which is as reliable in that extreme environment as it is for spirals. There is a clear excess of ultraviolet emission compared to individual HII regions in spirals, i.e. the
and
SFR ratios are on average low although there are some large variations from one region to another, which cannot be explained by variations of the metallicity or the dust extinction along the HI structure. Comparing the observed SFR with a model of the evolution of
with time favours young, quasi-instantaneous though already fading starbursts. The total star formation rate measured in the intergalactic medium (which accounts for 80% of the total) surrounding NGC 5291 is up to 1.3
yr-1 - a value typical for spirals - assuming the standard SFR calibrations are valid. The SFR drops by a factor of 2 to 4 in case the star formation is indeed quasi-instantaneous.
Key words: galaxies: individual: NGC 5291 - stars: formation - HII regions - intergalactic medium
It has long been known that "vast and often very irregular swarms of stars and other matter exist in the spaces between the conventional spiral, elliptical, and irregular galaxies'' (Zwicky 1951). Those stars have been stripped off from their host galaxies during interactions, the intergalactic medium being too diffuse and hot to allow in situ star formation. However a fraction of the stellar population could nevertheless have been formed locally out of the gas which has been pulled out from peripheral regions of gas rich galaxies. This has been shown to occur in simulations of tidal encounters (Duc et al. 2004; Hibbard & Mihos 1995) as well as in the case of collisional rings (Horellou & Combes 2001; Appleton & Struck-Marcell 1996). This gas which has the metallicity of the parent galaxy (roughly solar metallicity in the case of normal star forming spirals, lower for dwarf galaxies, see Mendes de Oliveira et al. (2006) for instance) can then be reprocessed to form new stars and even new dwarf-like galaxies if the expelled gas reservoir was massive enough. How star formation (SF) proceeds in collisional debris between galaxies, in an environment which is isolated from the direct influence of the galactic disks, and how it compares with more conventional galactic star formation has so far not been studied in detail.
The number of reported "intergalactic'' star forming regions has grown rapidly in the last few years, especially in groups (Mendes de Oliveira et al. 2004) and clusters of galaxies (Cortese et al. 2006). The area around NGC 5291 - an interacting system located at the edge of the Abell 3574 cluster - is one of the most spectacular objects where numerous instances of intergalactic HII regions are found.
The system is composed of an early-type galaxy, NGC 5291, and a highly disturbed companion galaxy, "the Seashell''. The first detailed observation of the system was carried out by Longmore et al. (1979) who detected on each side of NGC 5291 blue clumps identified as giant extragalactic HII regions. The spectra of some of them were obtained thanks to an object prism survey of HII galaxies (Maza et al. 1991; Peña et al. 1991). Malphrus et al. (1997) used the Very Large Array to map the neutral gas of the system. They discovered a huge asymmetrical ring-like structure connected to NGC 5291. It turned out to be one of the two most gas rich objects in the Southern hemisphere (Koribalski et al. 2004). The HI structure consists of a series of sub-condensations which correspond to the previously identified intergalactic HII regions.
Duc & Mirabel (1998) published a detailed optical (imagery and spectra) and near-infrared study of the system. They found that the HII regions had moderately high metallicity (
)
and could thus exclude a primordial origin for the HI structure. They claimed that the ring structure resulted from the past collision of an extended HI disk with a companion. They also showed that the most recent starbursts along the ring are younger than
years without any evidence for an underlying old stellar population. Braine et al. (2001) obtained with the SEST antenna CO(1-0) and CO(2-1) spectra of two of the brightest condensations. Inferring the molecular gas content from the CO emission and comparing it with the SFR determined from H
,
they determined that their star formation efficiency (SFE) was closer to that measured on average in spiral galaxies than in dwarfs of the same luminosity. Incidently, the detection of the millimetre lines was a further indication that the gas around NGC 5291 had been previously enriched. Finally, Higdon et al. (2006) performed mid-infrared observations, both imaging and spectroscopy, using the IRAC and IRS instruments on the Spitzer Space Telescope. They unambiguously detected Polycyclic Aromatic Hydrocarbon (PAH) features towards the two brightest condensations of the system. Their spectral analysis also confirmed the young ages of the star forming episodes.
![]() |
Figure 1:
Composite image of the NGC 5291 system in pseudo-colours: 8.0 ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
In order to study in more detail the characteristics of the SF regions around NGC 5291, we collected ultraviolet images from the GALEX satellite, ground-based H
images, and high resolution HI maps obtained with the Very Large Array which we combined with optical and mid-infrared Spitzer images of the system. See Fig. 1 for a combination of those images.
In Sect. 2 we present the multiwavelength data set we have used. In particular, we describe how we defined and selected the individual intergalactic star forming regions and how we performed accurate aperture photometry on them. In Sect. 3, we describe their main properties: overall distribution and star formation rates inferred from different tracers. We then compare them with those measured in various samples of galaxies and individual star forming regions in spiral disks. In Sect. 4, we discuss the origin of the discrepancies between the various SF indicators and present the possible role of extinction, metallicity and starburst age. We finally speculate on the possible impact of intergalactic star formation on its environment.
Throughout this paper we use H0=72 km s-1 Mpc-1 which results in a distance for NGC 5291 of 62 Mpc. At this distance, 1
= 0.30 kpc.
We describe in this section the various observations used for this paper. We present for the first time an ultraviolet map obtained with the GALEX space telescope. The high resolution HI maps were taken from the paper by Bournaud et al. (2007). The ESO/3.6 m H
data and Spitzer/IRAC mid-infrared images were obtained from the literature or from the telescope archives.
On May 2005 GALEX was used to acquire a set of NUV (near-ultraviolet;
nm and
nm) broadband images of the target with a total exposure time of 2886 s. The field of view has a diameter of 1.24
,
much larger than the NGC 5291 system. The PSF (point spread function) has a width of
5.0
.
Due to technical difficulties at the time of observation, no far ultraviolet image (FUV) was obtained.
H
observations were carried out at the 3.6-m telescope of the European Southern Observatory (ESO) on April 2002, using the Fabry-Perot interferometer CIGALE. The seeing was
1
.
The zeroth moment maps were extracted from the datacubes. For further technical details about the observations and data reduction, see Bournaud et al. (2004).
Due to the small field of view of the camera, three individual images were acquired, covering respectively the area centred on the most luminous HII region to the south of NGC 5291, the central area close to NGC 5291 and the Seashell galaxies, and the area centred on the most luminous HII region to the north of NGC 5291. The southern and central images overlap partially but there is a gap between the central and northern images, as shown in Fig. 1.
The mid-infrared broadband images have been retrieved from the Spitzer public archives. The observations were carried out on February 2004 with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope at 3.6 m, 4.5
m, 5.8
m and 8.0
m. The resolution is
3
.
The data have been preprocessed using the Spitzer Science Center (SSC) pipeline version S11.0.2. Due to the limited field of view (about 5
), the final image is a combination of two different pointings (each using a Reuleaux dither pattern), separated by 180
(see Fig. 1). For further details about Spitzer data processing see the Spitzer Observing Manual
and the paper by Higdon et al. (2006).
We covered the entire velocity range of the system at 48.8 kHz (10.45 km s-1), after online Hanning smoothing. The AIPS
package was used to calibrate and Fourier transform the data. For this paper we used natural weighting which results in maps at
arcsec2 resolution (and a beam position angle of 13
). The
rms noise is 0.33 mJy (corresponding to 3.8 K).
As the emission in this system is extended over more than 10 arcmin, some flux due to a smooth, extended HI component might be missing from the interferometer data. It should be stressed that for spatial frequencies corresponding to scales smaller than 5
,
all HI is recovered in our maps and that any contribution from an extended component becomes negligible in the small regions studied in this paper.
To identify the intergalactic star-forming regions belonging to the NGC 5291 system and reject from the subsequent analysis unrelated objects we used a combination of optical (B, V and R bands from Duc & Mirabel 1998), H
Fabry-Perot, mid infrared (3.6, 4.5 and 8.0
m) and HI observations. We used the following criteria in this order:
As HII regions are irregularly shaped and the system is crowded, we have chosen to define polygonal apertures. We could thus simultaneously take into account the variation of the PSF from one image to another avoiding as much as possible pollution from stars and from nearby objects, which would not have been easily possible using circular apertures. The shape and angular sizes of the polygonal apertures are the same for all bands and thus cover exactly the same physical region.
Based on the above mentioned criteria, we have selected 29 regions across the field, which correspond to intergalactic star-forming regions. They are indicated in Fig. 3.
For reference, we have also measured the fluxes towards the parent galaxies (NGC 5291 and the Seashell) and towards an HI peak just North of NGC 5291 with only faint ultraviolet and mid-infrared counterparts (information on H
is lacking there).
Table 1:
Fluxes of selected regions in the ultraviolet, H,
and mid-infrared (3.6, 4.5, and 8.0
m). The "-'' symbol indicates that the region falls outside the field of view, "
'' means the aperture is heavily polluted by stars and the flux is therefore likely being overestimated and "
'' means the aperture does not contain all of the object, e.g., to avoid a foreground object, the actual flux then being underestimated.
![]() |
Figure 2:
The different tracers of star formation towards NGC 5291: VLA HI map ( top left), GALEX near-ultraviolet ( top right), Fabry-Perot H![]() ![]() ![]() ![]() |
Open with DEXTER |
First of all, we have substracted the background in each band. To do so, we have manually measured the background level on each image using the IRAF (Image Reduction and Analysis Facility) IMEXA procedure. Depending on the image, a single background was determined in case the background level was flat throughout the image (the ultraviolet image falls in this category as the standard pipeline already performs a background subtraction). Individually determined background levels had to be used for each of the regions in the H
and infrared images. Each level is evaluated as the mean of several (up to 30) measurements of the background at different locations around the selected regions. The uncertainty on the background is the standard deviation.
We have also tested background removal using SExtractor (Bertin & Arnouts 1996). However, the large scale diffuse luminosity made the determination of the level quite uncertain and a small variation of the background extraction parameters caused a large variation of the flux, especially for the faintest regions where the difference between two extractions could be up to a factor of 3.
The standard IRAC mid-infrared pipeline calibrates the images in such a way that the flux of a point source is correctly determined when using a 10 pixel radius circular aperture and calculating the background in an annulus with inner and outer radii of 10 and 20 pixels. This is due to the fact that the PSF is quite large, especially in the 8.0
m band where a significant part of the flux is scattered on arcminute size scales, which is much larger than 10 pixels. The calibrated units are set in such a way that the off aperture flux loss is exactly compensated when it is determined with a 10 pixel circular aperture. Therefore, when the aperture is larger (smaller) than 10 pixels of radius, the flux density is overestimated (underestimated). However in this paper we do not use a circular aperture for the determination of the fluxes. Instead, the areas of interest are enclosed using polygons. As a result another method needs to be applied in order to estimate the flux outside our apertures and perform the necessary corrections. To do so we use the standard correction values published in the IRAC manual for a point source correction, taking an equivalent radius of
,
S being the surface in pixel units. The equivalent background annulus was taken with a radius from 10 to 20 pixels. Corrections for aperture areas falling between the published values have been interpolated. The uncertainty on the aperture correction therefore depends on the aperture and the HII region shapes. The Poisson photon noise has also been taken into account. The absolute photometric calibration, according to the IRAC data handbook, is better than 5%.
As the Fabry-Perot observations were not flux calibrated, we made an indirect calibration of the H
data using the optical multi-object spectroscopic observations taken on July 1994 on the NTT ESO telescope at La Silla (Duc & Mirabel 1998). The calibration is accurate to only about 30% due to uncertainties on the slit positions. The overlap of the two southern fields of view has been used to check the consistency of the calibration between the individual Fabry-Perot pointings.
The absolute photometric uncertainty of the GALEX near ultraviolet image is 7%.
The fluxes derived in each band - when available - are presented in Table 1. The mid-infrared fluxes obtained have been compared whenever possible with the ones published by Higdon et al. (2006). The integrated fluxes of NGC 5291, the Seashell and PGC 048877 are consistent to within a few percent. Our photometry of individual star-forming regions agrees with the earlier values to within 30% for the case of bright isolated regions, but it deviates substantially for the fainter ones. Our analysis indicates that the differences mainly originate from the way apertures were selected and how the background substraction was performed: automatically by SExtractor in the Higdon et al. study, manually in this paper to make sure that regions with the same physical origin were enclosed.
The brightest regions (in H,
NUV and 8.0
m) are located near the apparent tip and points of maximum column density of the HI ring: they correspond to our regions 1, 3 and 21; the latter two, also known as NGC 5291S and NGC 5291N, have been discussed in the literature as the possible progenitors of dwarf galaxies, and are probably associated with gravitationally bound gas clouds (Bournaud et al. 2004). As shown in Fig. 2, the HI condensation just North of NGC 5291, which actually has one of the highest HI column density in the entire system (
cm-2, hereafter referred to as HI concentration) shows only faint diffuse ultraviolet and mid-infrared counterparts
. The UV/HI and IR/HI flux ratios measured within this condensation deviate by an order of magnitude from the values determined from the other regions. One possibility would be that the starburst is extremely young and whatever O and B stars formed there have not dispersed their surrounding dense cloud material yet. Having a CO spectrum of that region and a far-infrared/millimetre image would help in checking this hypothesis. Having said that, inspection of individual HI spectra across this cloud indicates that HI is found spread over a much wider range in velocities than at the other locations investigated. This could mean that what we are seeing in Fig. 2 is a blend of clouds along the line of sight, adding up to result in an apparent massive HI concentration whereas each contributing cloud is much less massive, one of which being just massive enough to have started to form stars.
Figure 3 presents combined NUV, H
and 8.0
m emission of the individual star-forming knots. The false-colour images, all displayed with the same intensity scale, indicate that although they have the same morphology in each band, they show large variations in their relative fluxes. Local variations of the properties of the star forming regions are thus important, despite the fact that they all belong to a single structure probably formed simultaneously during a collision.
Another region, labelled number 1 in Fig. 3, also shows a level of activity much higher than what would be expected from its HI column density. However the HI distribution of this object shows clear distortions: a gaseous tail without any 8.0 m counterpart is seen just to the West of the star-forming region (see Fig. 1). There is evidence that in addition to the gravitational effects of the collision, the objects in the system may have suffered from the interaction of their interstellar medium with the intergalactic medium belonging to the cluster through which they move. A systematic offset between UV/H
/MIR and HI peaks is observed with the HI being shifted to the opposite side of the cluster centre, which is located to the East of the galaxy. Why this process would be more efficient for object 1, which belongs to the same global structure as the other objects is however unclear.
![]() |
Figure 3:
Ultraviolet image of NGC 5291 with a mosaic of the selected regions. The individual regions and the main galaxies in the group are labelled. The thumbnails are zoomed images of each selected region using pseudo-colours (red representing 8.0 ![]() ![]() |
Open with DEXTER |
The ultraviolet, H
and infrared 8.0
m emission are all known to be tracers of the level of star forming activity. The reliability of measuring the star formation rate based on each of these indicators, or a combination of them, has been extensively debated in the literature. It seems to depend very much on the type of environments as many external factors may affect it: contamination by old stars, level of dust extinction, etc. These tracers have not yet been calibrated for the special environments of collisional debris and IGM. We will thus, as a first step, use the classical calibrations obtained for spiral disks to get estimates of the SFR. From a comparison between the various estimates, we will then both discuss the validity of such tracers and see whether some global properties, such as the starburst age may be constrained.
UV radiation is mainly emitted by massive stars younger than
years. The star formation rate has been calculated using the standard Kennicutt (1998) relation:
[W Hz-1]
yr-1. This estimator has been built under the assumptions of continuous star formation, with a Salpeter initial mass function (IMF; Salpeter 1955) and mass cut-offs from 0.1 to 100
.
It is sensitive to extinction which reprocesses a fraction of the UV emitted by young stars which is reemitted in the infrared.
Ionizing radiation (
nm) is mainly emitted by very massive, short-lived young stars. This radiation ionizes surrounding gas clouds, which reemit partly in H
,
thus tracing starbursts. We have used again Kennicutt (1998):
(H
)
[W]
yr-1, with the same IMF as used previously. This estimator is sensitive to absorption of the ionizing photons by the molecular clouds surrounding young stars and thus very dependent on geometry
.
At 8.0 m we observe a combination of PAH broad emission bands, continuum emission from hot VSG (very small grains) and stellar continuum. The photospheric stellar contribution at 8.0
m is lower than 10% in our special environment which is mostly devoid of an old stellar component (see below). The PAH are heated stochastically by far UV radiation from massive young stars and as such trace star formation. Wu et al. (2005) built a SFR estimator based on a correlation between the 8.0
m dust emission (stellar contribution being substracted using the 3.6
m band) and radio or H
luminosities. The star formation rate calibration derived by them using the 8.0
m/radio correlation is:
yr-1. The Chary-Elbaz templates (Chary & Elbaz 2001) and the Kennicutt (1998) estimator based on the total infrared luminosity produce consistent results.
The star formation rates have been calculated using the relations given above and the fluxes presented in Table 1. We have corrected the fluxes for galactic extinction both in UV and in H
using the Cardelli extinction curve (Cardelli et al. 1989), and taking
E(B-V)=0.064 from NED (NASA/IPAC Extragalactic Database). The systematic uncertainty of each estimator has not been taken into account. The UV SFR,
and
for each region is presented in Table 2.
We stress that those estimators are likely not very adapted for our case (we expect a very young starburst, see Higdon et al. 2006). However they give more intuitive values to study the variations of the ratio of the SFRs in two bands, from one region to another. We use SFR(UV) as the reference since it is the most direct estimator (unprocessed star light) and more importantly, this is the only estimator which samples the star formation emission along the whole HI structure.
Table 2:
Star formation rates deduced from the ultraviolet, H,
and mid infrared, in solar mass per year. The symbol ":'' means the estimation is uncertain.
The star formation rate derived from the NUV summed along the entire HI structure but excluding that of NGC 5291 and the Seashell galaxies amounts to
yr-1. The SFR of NGC 5291 and the Seashell is only 0.35
yr-1, or
20% of the SFR in the entire field of view. Thus, star formation in this system is currently mostly observed in the intergalactic medium.
The total star formation rate in the regions in which we have measured fluxes in all three bands is:
To compare the different SFR estimators, we plot
(Fig. 4) and
(Fig. 5) versus SFR(UV). On average, we get
and
.
![]() |
Figure 4:
Ratio of
![]() |
Open with DEXTER |
![]() |
Figure 5:
Ratio of
![]() |
Open with DEXTER |
Figures 4 and 5 indicate strikingly large variations of
and
from one region to an other, even among the most luminous ones. For instance, there is a factor
2 increase in
when comparing region 20 with 21, or a factor
6 decrease in
when going from region 1 to 3. Variations can also be important within a unique region, see for instance region 1 in Fig. 3. As discussed in Sect. 4, local variations of the star formation history, dust extinction and metallicity can in principle explain this scatter. This will be investigated further in the discussion.
The intergalactic star-forming regions we have studied are objects intermediate between isolated star-forming dwarf galaxies and individual HII regions in spiral disks. Indeed, the most luminous of them have simultaneously the structural properties and optical colours of typical primordial blue compact dwarf galaxies and the interstellar medium content, in particular the metallicity, of their parent disks. It therefore makes sense to compare the individual SF regions around NGC 5291 with both samples of individual galaxies, studied globally, and samples of individual HII regions within nearby galaxies (M 51, M 81 and Arp 82). Fortunately, several Spitzer based papers have recently provided the data required for these comparisons.
We first compare the properties of the intergalactic star-forming regions in NGC 5291 with those measured globally for star-forming galaxies. For this comparison, we used three samples: the one studied by Rosenberg et al. (2006), consisting of 19 galaxies with a B-band absolute magnitude MB>-18.0; the sample of galaxies, mostly spirals, in the Spitzer SINGS legacy survey (Dale et al. 2007); and the Engelbracht et al. (2005) sample consisting of a heterogeneous collection of star forming galaxies without strong active nuclei including metal deficient blue compact dwarf galaxies. In Fig. 6, we plot the 8.0 m fluxes, normalised to the 4.5
m fluxes (a tracer of the hot dust emission) as a function of the monochromatic luminosity at 8.0
m. The 4.5
m band was corrected for contamination by the photospheric stellar emission determined from the 3.6
m band (which we assumed not to be significantly polluted by hot dust):
,
where
is a scaling coefficient. We adopted
as used in Engelbracht et al. (2005). This is derived from Starburst99 models (Leitherer et al. 1999) which give typical values between 0.53 and 0.61; Engelbracht assigned
an uncertainty of 7%.
![]() |
Figure 6:
Ratio of the mid-infrared fluxes
![]() ![]() |
Open with DEXTER |
Table 3: Mean flux and SFR ratios for the different samples. Extremely discrepant points (if any) have been discarded before performing the calculation.
Most of the individual regions of NGC 5291 exhibit roughly the same range of monochromatic luminosities at 8.0 m as the Rosenberg et al. sample of dwarf galaxies except for the brightest ones (see Fig. 6). Based on their SFRs only, each of the 29 individual regions could by itself be considered a dwarf galaxy. Besides, they have on average a
ratio slightly, but not significantly, higher than the dwarf galaxy sample (see Table 3). The ratio for the SINGS is slightly larger but very similar. The Engelbracht et al. sample, in spite of its inhomogeneity, shows a similar mean value as for the NGC 5291 regions.
In Fig. 7 we compare
to
for the galaxies of the SINGS sample and the HII regions detected in the NGC 5291 system.
![]() |
Figure 7:
Star formation rate estimated from the IRAC 8.0 ![]() |
Open with DEXTER |
We notice that the least luminous galaxies from the SINGS sample present an important ultraviolet excess. We think it is due to the PAH defiency seen in the low metallicity galaxies, rather than an ultraviolet excess per se. On the other hand, the most luminous SINGS galaxies present a strong infrared excess, due to the higher metallicity and the large population of non-ionizing stars which can excite the PAH in addition to massive young stars.
The relative scatter (see Table 3) for all these samples are quite similar. This is surprising as the Rosenberg et al., the Englebracht et al. and the SINGS samples are composed of heterogeneous dwarf galaxies whereas the NGC 5291 HII regions lie in the same HI structure.
In Fig. 8 we perform a similar comparison as in Fig. 6 but with HII regions in the nearby spiral M 81 (for more details see Pérez-González et al. 2006) and the interacting system Arp 82 (Hancock et al. 2007).
![]() |
Figure 8:
Ratio of the mid-infrared fluxes
![]() ![]() |
Open with DEXTER |
The
ratio of the regions in M 81 is roughly similar to that of the most luminous HII regions of NGC 5291 even though their 8.0
m monochromatic luminosity is lower by one to two orders of magnitude. Beside, for the less luminous regions, the ratio is higher in M 81 which could indicate a deficiency in PAH emission in the HII regions of NGC 5291. The star forming regions in Arp 82, especially along the tidal arms of the interacting system, follow more closely what can be seen in NGC 5291 both in the 8.0
m flux and the
.
Therefore both systems have on average the same PAH properties compared to the dust continuum.
In Figs. 9 and 10, respectively, we present a comparison of the
to
and
to
SFRs for individual HII regions in the HI intergalactic structure of NGC 5291, for the HII regions in the disk of M 51 (Calzetti et al. 2005), M 81 (Pérez-González et al. 2006) and Arp 82 (Hancock et al. 2007).
![]() |
Figure 9:
Star formation rate estimated from the IRAC 8.0 ![]() ![]() |
Open with DEXTER |
We noticed in Sect. 3.2 that the SFRs of the NGC 5291 HII regions estimated from the UV fluxes were on average higher by a factor of 2 than those obtained from the H
and 8.0
m bands
. Figures 9 and 10 show that this systematic UV-excess is not observed for the individual HII regions of M 51 and M 81; the regions in M 51 present an IR excess, whereas those of M 81 seem to be too H
rich. No systematic UV excess can be detected in the Arp 82 interacting system.
![]() |
Figure 10:
Star formation rate estimated from the
![]() ![]() |
Open with DEXTER |
More quantitatively, when normalised to the 8.0 m band, the NUV flux in NGC 5291 is in excess by a factor of 8, 2 and 1.6 compared to M 51, M 81 and Arp 82 respectively. When normalised to H
,
the corresponding excess factors are 2, 6 and 2. Finally, one should note that NGC 5291 presents a scatter in the H
vs. UV relation which is smaller than for the comparison galaxies, while the MIR vs. UV is roughly comparable.
In summary, while globally the 8.0 m emission in the intergalactic star forming regions of NGC 5291 are not significantly different than those of star-forming dwarf galaxies, they seem to present with respect to individual HII regions in spirals, and with the H
and MIR SF indicators, a clear UV excess.
In the previous section, we presented the properties of the star forming regions along the huge HI structure around NGC 5291, a study based on three SF indicators: UV, H
and the mid-infrared. We found that the three of them were correlated, as for more classical galactic star-forming regions, but showed with the latter some differences, in particular a UV excess. We discuss here the origins of these differences and of the scatter in the correlations. The goal is to determine which SF indicator - or combination of them - may be used to estimate the total SFR in this special environment, bearing in mind that the previous calibrations were derived for galactic environments. Once a "total'' SFR is determined, we briefly discuss its impact on the IGM.
In most galactic star-forming regions, the dust emission around 8.0 m is largely dominated by PAH features. The latter are however known to be absent in metal-poor objects such as the Blue Compact Dwarf Galaxies (Madden et al. 2006; Wu et al. 2006; Galliano et al. 2005; Engelbracht et al. 2005). Thus for these systems, the IRAC based mid-IR-SFR calibration, as the one used in this study, could be fairly inaccurate. The intergalactic HII regions around NGC 5291, although they share many properties with those of BCDGs, lie however within an HI structure which had been previously pre-enriched. Duc & Mirabel (1998) measured an oxygen abundance of
to 8.6, with an average of
8.49+0.07-0.10. Although this is somewhat below the solar value, this is consistent with a PAH detection. Indeed, the IRS spectra obtained by Higdon et al. (2006) indicate that the IRAC 8.0
m emission of the two most luminous star forming regions is mainly due to PAHs.
![]() |
Figure 11:
Ratio of the mid-infrared fluxes
![]() ![]() |
Open with DEXTER |
In Fig. 11, we plot the PAH emission normalised to the dust continuum as a function of metallicity for the NGC 5291 star forming regions, the dwarf galaxies from the Rosenberg sample and the galaxies from the Engelbracht sample. The luminosity at 8.0 m appears to be consistent with that of dwarf galaxies with the same metallicity. Wu et al. (2006) directly measured a correlation between the equivalent widths of the PAH features and metallicity. We checked this for region 21. We found that the EWs measured in the IRS spectrum published in Higdon et al. (2006)
, using a method similar as in the Wu et al. (2006) paper, were consistent with the Wu et al. correlation.
Thus, on average, the relative strength of the 8.0 m emission in NGC 5291 appears to be normal, i.e. consistent with that expected for star-forming regions with the same metallicity. Therefore the calibration derived between the IRAC mid-infrared emission and the SFR so far only validated for spiral disks should also apply to the collisional debris.
The large scatter in the
flux ratio along the system (see Figs. 6 and 8) is however surprising. In any case, it cannot be simply explained by local variations of the metallicity as the latter was measured by Duc & Mirabel (1998) to be quite uniform all along the HI structure. For instance, very bright regions, such as 1 and 3, show discrepant
flux ratios although they lie very close on the sky and thus should have very similar metallicities. Unfortunately no optical spectrum of region 1 is available to check this. The IRS spectra of the two most luminous regions (regions 3 and 21 in this paper) also show variations in their PAH strengths and dust continuum shortwards of 5
m and above 9
m, which are not yet understood. Clearly, the metallicity is not the only parameter playing a role in the 8.0
m emission, as also shown in several studies of galactic star formation.
To quantify this HI based extinction, we used previous extinction estimates determined by Duc & Mirabel (1998) from the measurement of the Balmer decrement in several associated HII regions. Correlating the peak HI column densities and the optical extinction in the V band with 7 data points (after having excluded 2 which showed large discrepancies), we obtained a calibration of the HI-based extinction. We estimate the uncertainty to be about 0.4 dex with this method.
To convert the extinction from the V band to other bands, we used the average Large Magellanic Cloud
of Gordon et al. (2003). We find
and
.
We have corrected the NUV and H
fluxes for dust extinction based on the calibrated HI column densities. Results of this correction are shown in Fig. 12 for objects with uncorrected ultraviolet star formation rate greater than 0.02
yr-1, which have the most secure determination of the SFR.
![]() |
Figure 12:
![]() ![]() ![]() |
Open with DEXTER |
The relative scatter of
and
is not reduced by the extinction correction. Clearly, variations of the dust extinction alone cannot explain all the scatter in the SF properties found along the HI structure, unless very local variations dominate.
The so-called "UV-excess'' seems to be a determining characteristic of the intergalactic star forming regions around NGC 5291. This was also observed in a few other intergalactic SF regions (Werk et al. 2006). We discuss here whether this excess, as well as the large scatter of the
after the correction for dust extinction, could result from some starburst age effects.
Duc & Mirabel (1998) measured H
equivalent widths as large as 140 Å of several HII regions around NGC 5291, which is indicative of young starbursts (see for instance Terlevich et al. 2004). The detection of Wolf-Rayet features in the most luminous one shows that the current star formation episode there is less than
years old. More recently, Higdon et al. (2006) determined ages of
years from [Ne II]/[Ne III] emission line ratios for the objects for which mid-infrared IRS spectra were available. Are these indications for starbursts being young consistent with the UV excess observed in all intergalactic HII regions around NGC 5291? At first sight, one may think that a UV excess, or conversely, an H
deficit would rather suggest older, declining bursts since the UV time scale is a factor of ten larger than the H
one (see Sect. 3.2).
To pursue the analysis, we have modelled the HII regions using the spectral evolution code P´EGASE II (Fioc & Rocca-Volmerange 1997,1999). We have used a set of parameters differing by the star formation history (SFH), going from an instantaneous starburst to exponentially decreasing ones with a timescale between 1 and 40 million years. We used a Salpeter IMF from 0.1 to 100
and a constant metallicity Z=0.008 (65% solar metallicity, corresponding to
), in agreement with the metallicity determined by Duc & Mirabel (1998). Extinction from the gas has not been taken into account in the simulations in order to isolate age effects only. For comparison we therefore need to use the extinction corrected ratios. We have in particular studied how the H
equivalent width and the
ratio evolve with SFH (see Figs. 13 and 14). The first plot suggests that the ongoing star formation may not have lasted for more than
years for the regions presently showing the largest H
equivalent width, whatever the timescale assumed for the starburst. Other than that, very young, quasi-instantaneous starbursts ignited
years ago may well cohabit with older and more extended star-formation episodes (Neff et al. 2005), i.e. in regions showing smaller H
equivalent width. The second plot shows that the observed range for
may be obtained for various SFH, but that in any case, all bursts are well advanced. A constant or just slowly declining SFR is inconsistent with the data. However, our dataset may accommodate quasi-instantaneous bursts which occurred
years ago up to bursts with slightly longer time scales of
years which started about
years ago. Combining the two criteria, we conclude that the regions that are currently the most active ones (i.e. region 1, 3 and 21) are all very young but should have had their SF peak already
years ago while the more quiescent ones may be a few
years older, still sustaining star formation but at a lower rate than before.
![]() |
Figure 13:
Plot of the variation of the equivalent width of H![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 14:
Plot of the variation of the
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
One should note that quasi-instantaneous young starbursts for all HII regions around NGC 5291 are also fully compatible with the data, i.e. the H
equivalent width, the UV excess and the somehow large region-to-region variations of the
to UV ratio. The latter may be understood noting how this ratio changes quickly after the beginning of the starburst: a small age difference can produce a considerable difference in
(see Fig. 14). Probing the presence, or actually the absence, of an intermediate age stellar population of a few
years would be required to validate the scenario of a rapidly decreasing starburst. This requires a model of the full Spectral Energy Distribution of each region, from the UV to IR, which is beyond the scope of this study. This task is somehow easier in the case of the NGC 5291 intergalactic star-forming regions. Indeed, they should not contain a stellar component of a few
years or older coming from the parent galaxies.
As all the HII regions seem to have roughly the same age it is legitimate to wonder how star-formation may have started quasi simultaneously along a structure as extended as 180 kpc, and how likely we are to observe it so close to its onset. One should however remember that the HI ring which hosts all these regions was probably formed during a single violent event, likely the collision of a previously existing extended HI disk with a companion galaxy (Duc & Mirabel 1998). Given its ring-like shape and overall velocity field, the HI structure was most probably not formed by tidal forces but rather by an expanding density wave following a bull's eye collision. According to numerical simulations (Bournaud et al. 2007), the ring around NGC 5291 has an age between
and
years, very different from the star formation timescale. However, in the case of a collisional ring and unlike tidal tails, the star formation can be suppressed during the first few
years, i.e., until the HI ring expansion stalls after which star formation proceeds. If that scenario is indeed correct, this would bring the dynamical and the star-formation time scales together and explain the, at first sight, rather intriguing observation of quasi simultaneous star formation along the entire ring.
The integrated star formation rate along the HI structure - uncorrected for dust extinction - derived from the GALEX NUV emission, amounts to
yr-1,
yr-1 in H
and
yr-1 from the 8.0
m band emission (see Sect. 3.2).
Whether the uncorrected UV, H
and IR SFRs should be added to get the total, dust free SFR, or whether a simple extinction correction for each tracer is enough is a matter of active debate (Hirashita et al. 2003). Because our estimates of the extinction, based on the HI column density, are rather uncertain, we decided not to apply a dust correction for each individual region. Instead we decided to add the uncorrected SFR obtained from H
to the SFR based on the dust emission scaled by a coefficient which takes into account the amount of ionizing radiation (actually radiation emitted by the same stars as traced through H
)
absorbed by PAH. The latter was determined empirically by fitting SFR(H
) + k
SFR(IR) to the extinction corrected SFR(H
)
using only datapoints for which measures of AV obtained from the Balmer decrement were available. The advantage of using a combination of H
and IR rather than UV and IR (as commonly done) is that the former two indicators probe SF with relatively similar time scales of
years (rather than
years for UV). Doing so, we obtained k=0.27 and a total SFR all along the HI structure of 0.83
yr-1. This value is a factor of 1.7 smaller than the value determined from the uncorrected near UV emission, confirming again the earlier mentioned UV excess.
However, we should note at this stage that the SFRs were estimated using calibrations valid for a constant star formation over
years (Kennicutt 1998). As mentioned before, the star formation history in NGC 5291 may have been quite different. There are various independent indications pointing towards young, quasi-instantaneous, starbursts and the measure of the actual SFR may be different.
For young and rapidly decreasing starbursts, models (see Figs. 15 and 16) indicate that the Kennicutt UV and H
calibrations overestimate by factors of 2 to 4 the actual SFR. Thus if the quasi-instantaneous starburst hypothesis is correct, the real "intergalactic'' star formation rate around NGC 5291 should be 0.3-0.7
yr-1, still higher than for typical dwarf and irregular galaxies (Hunter & Elmegreen 2004). However, one should keep in mind that the UV-excess indicates that the star formation rate was higher a few 106 years ago, with values close to that observed in quiescent spiral galaxies.
![]() |
Figure 15:
Plot of the variation of
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 16:
Plot of the variation of the
![]() ![]() ![]() ![]() |
Open with DEXTER |
We thus find that around NGC 5291, star-formation proceeds at a rate typical of spirals but in an environment which is much less dense, and where the effect of stellar winds on the surrounding medium should be much more dramatic. As shown, among others, by Ferrara & Tolstoy (2000), the ability of an object to lose its heavy elements, and thus enrich the intergalactic medium, depends very much on the mass and shape of the halo it sits in. Collisional debris being dark matter deficient, they should have a reduced capability to retain their metals. How the particularly active star-forming regions in NGC 5291 manage to survive this feedback of mechanical energy still remains to be studied. Given their large gas reservoir, it is somewhat surprising that they are currently already fading, unless SF has ceased as a result of feedback from massive stars.
Clearly, a system like NGC 5291 is exceptional in the local Universe. Studying a sample of interacting systems (Hancock et al. 2007; Struck 2006; Smith et al. 2007) found that their external plumes and tails only account for at most 10% of the total star formation. In NGC 5291, from the UV band we can see that more than 80% of the stars are currently formed along the intergalactic ring. When using the 8.0 m band, removing the old stellar content from the NGC 5291 galaxy (this population is responsible for roughly 50% of the 8.0
m flux), the intergalactic SF still accounts for more than 50% of the SFR. In the distant Universe, however, such large gaseous structures may have been more common
and the star formation regions they hosted may have contributed significantly to the direct enrichment of the IGM.
We have performed a multi-wavelength study of the NGC 5291 HII regions using ground-based, GALEX, and Spitzer Space Telescope data of 29 intergalactic HII regions around NGC 5291. These regions lie along a previously identified huge HI ring-like structure of collisional origin. We have measured their fluxes in the near ultraviolet, H,
and mid-infrared in order to compare these different star formation tracers. The morphology observed in the three bands is remarkably similar, indicating that each can serve as a star formation tracer. We have shown that the intergalactic IR emission at 8.0
m, dominated by PAH, when normalised to the hot dust continuum, is comparable to the integrated emission of dwarf galaxies of the same metallicity and to the emission of individual HII regions in spirals. This indicates that the 8.0
m band emission is probably an estimator of the SFR that is as reliable for these regions as it is for spirals. Whereas the 8.0
m is "normal'', there is a clear excess of near ultraviolet emission compared to individual HII regions in spirals. The
and
are low although there are some strong variations from one region to an other. The variations of the metallicity sensitive
ratio cannot be explained by a metallicity effect since the oxygen abundance was found to be uniform along the entire HI structure. Moreover, correcting for the spatial variations of the dust extinction does not reduce the scatter in the
and
SFR ratios. Their variations are best explained by age effects. A model of the evolution
with time favours young instantaneous starbursts: the ultraviolet excess or deficit of H
plus the large scatter indicating rapid variations of the
ratio are difficult to explain with a constant SFR or moderately decreasing SFR. A young starburst is consistent with constraints from the H
equivalent widths and mid-infrared NeII fluxes. So far no indication for the presence of an old stellar component has been found. A detailed study of the full SED of the star-forming regions, including dust models, would be required to confirm it. Comparing the individual starburst ages with the dynamical timescale for the formation of the HI ring in which they lie, we concluded that the apparent observation of a quasi-simultaneous onset of star formation over the 180 kpc long gaseous structure can be understood if SF commences only after expansion of the ring has stalled, after several
yr.
Using the standard Kennicutt calibrations, we estimated from the UV, uncorrected for dust extinction, a total SFR around NGC 5291 of 1.3
yr-1. If the quasi-instantaneous starburst hypothesis is correct, then this current value of the SFR may have been overestimated by factors of 2 to 4. In that case, the SFR would have been much higher in the recent past. Finally it appears that about 80% of star formation in the NGC 5291 system takes place in the intergalactic HII regions, a value much higher than in nearby interacting systems.
Acknowledgements
This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.
We thank Pablo Pérez-González, Daniel Dale and Mark Hancock for providing additional data on the comparison galaxies. We also thank our referee, Claudia Mendes de Oliveira, as well as Véronique Buat, Frédéric Bournaud and Frédéric Galliano for useful comments and enlightening discussions.
U.L. acknowledges support by the DGI grants AYA 2005-07516-C02-01 and ESP 2004-06870-C02-02, and by the Junta de Andalucía (Spain).