A&A 477, 491-501 (2008)
DOI: 10.1051/0004-6361:20077112

Star formation in the hosts of GHz peaked spectrum and compact steep spectrum radio galaxies

A. Labiano1,2 - C. P. O'Dea3 - P. D. Barthel2 - W. H. de Vries4,5 - S. A. Baum6


1 - Departamento de Astrofísica Molecular e Infrarroja, Instituto de Estructura de la Materia (CSIC), Madrid, Spain
2 - Kapteyn Astronomical Institute, Groningen, 9700 AV, The Netherlands
3 - Department of Physics, Rochester Institute of Technology, Rochester, NY 14623, USA
4 - University of California, Davis, CA 95616, USA
5 - Lawrence Livermore National Laboratory, IGPP, Livermore, CA 94550, USA
6 - Center for Imaging Science, Rochester Institute of Technology, Rochester, NY 14623, USA

Received 17 January 2007 / Accepted 29 October 2007

Abstract
Aims. We are searching for star formation regions in the hosts of potentially young radio galaxies (gigahertz peaked spectrum and compact steep spectrum sources).
Methods. We used near-UV imaging with the Hubble Space Telescope Advanced Camera for Surveys.
Results. We find near-UV light could be the product of recent star formation in five of the nine observed sources, though other explanations are not currently ruled out. An additional two sources show marginal detections. The UV luminosities of the GPS and CSS sources are similar to those of a sample of nearby large-scale radio galaxies. Stellar-population synthesis models are consistent with a burst of recent star formation occurring before the formation of the radio source. However, observations at other wavelengths and colors are needed to definitively establish the nature of the observed UV light. In the CSS source 1443+77, the near-UV light is aligned with and is co-spatial with the radio source. We suggest that the UV light in this source is produced by star formation triggered and/or enhanced by the radio source.

Key words: galaxies: active - galaxies: starburst - galaxies: evolution - galaxies: stellar content - galaxies: interactions - ultraviolet: galaxies

1 Introduction

The relationship between black hole mass and galaxy mass implies that the growth and evolution of black holes (therefore AGN) and their host galaxies must somehow be related (e.g., Gebhardt et al. 2000). Mergers and strong interactions can trigger AGN activity in a galaxy (e.g., Baum et al. 1992; Israel 1998; Heckman et al. 1986). These events can also produce instabilities in the ISM and trigger star formation (e.g., Ho 2005). Numerical simulations and models (e.g., Mellema et al. 2002; Rees 1989) suggest that the advancement of the jets through the host galaxy environment can also trigger star formation. Imaging studies in the ultraviolet (UV) light of large 3CR sources find evidence for episodes of star formation starting around the time when the radio source was triggered (i.e. $\la$ 107-108 yr, Koekemoer et al. 1999; Allen et al. 2002; O'Dea et al. 2001; Martel et al. 2002; Chiaberge et al. 2002; O'Dea et al. 2003), suggesting a possible link.


  \begin{figure}
\par\mbox{\includegraphics[width=4.5cm]{figs/7112fg1a.eps}\hspace...
...hspace*{1.2cm} \includegraphics[width=4.5cm]{figs/7112fg1o.eps} }\end{figure} Figure 1: ACS HRC F330W image ( left panels), GALFIT model ( center panels) and residuals ( right panels) of our sample. Arrows indicate the direction of the radio source. The dashed lines in 1443+77 sketch the arc of emission mentioned in Sect. 4.1.
Open with DEXTER


 \begin{figure}
\par\mbox{\includegraphics[width=4.5cm]{figs/7112fg1p.eps}\hspace...
...hspace*{1.2cm} \includegraphics[width=4.5cm]{figs/7112fg2z.eps} }\end{figure} Figure 1: continued.
Open with DEXTER

Gigahertz peaked spectrum (GPS) and compact steep spectrum (CSS) radio sources are apparently young, smaller (GPS $\la$ 1 kpc, CSS $\la$ 15 kpc, for a review see O'Dea 1998) versions of the large powerful radio sources, so they are expected to exhibit signs of more recent star formation. In addition, their subgalactic size makes them excellent probes of the interactions between the expanding lobes and the host. They have not completely broken through the host ISM, so these interactions are expected to be more important than in the larger sources.

Near UV observations are very sensitive to the presence of hot young stars and as such will trace recent star formation events. We have therefore obtained high-resolution HST/ACS near-UV images of these young compact sources to study the morphology and the extent of recent star formation.

This is the first time a sample of GPS and CSS sources has been imaged in the near-UV. It is also the first time that the relative sizes of radio sources in well-matched samples are used to study time evolution of merger-induced and jet-induced star formation.

2 Observations and data reduction

We have obtained high resolution near-UV images with the High Resolution Channel (HRC) of the Advanced Camera for surveys (ACS) on board the Hubble Space Telescope, through the F330W filter, with integration times of 1800 s. The objects observed are GPS and CSS galaxies 1117+146, 1233+418, 1345+125, 1443+77, 1607+268, 1814-637, 1934-638, 1946+708, 2352+495 (Table 1). Our sample is chosen to be representative of GPS and CSS sources with $z\la 0.5$, and nearby enough to eliminate strong effects due to evolution with cosmic time. The objects are drawn primarily from the well-defined samples of Fanti et al. (1990), Fanti et al. (2001), Stanghellini (1992) and Stanghellini et al. (1997). The comparison sample of large 3CR sources consists of FR I and FR II sources with redshifts less than 0.1 observed in the near-UV by Allen et al. (2002).

The standard ACS reduction pipeline was used to remove detector signatures such as bias, dark current, flat field and to perform flux calibration. Each target was observed in a three-point dither pattern. The frames were combined with Multidrizzle (Koekemoer et al. 2002) to correct for geometric distortions and cosmic rays. The resulting ACS images have a signal to noise of $\sim$100 for the brightest objects and $\sim$5 for the faintest.

The 2D fitting code GALFIT (Peng et al. 2002) was used to parameterize the UV emission. For each image we tested different combinations of point source and Sersic profiles, allowing the sky level, position and magnitudes of all components, as well the index and effective radii of the Sersic components, to vary. The final model was chosen according to the lowest  $\chi^2$ and best residuals (with the lowest number of components). The results are summarized in Table 2. Figure 1 shows the UV image, GALFIT model and residuals for the nine objects. These will be discussed in detail in Sect. 4.

A Sersic profile (Sersic 1968; Graham & Driver 2005; Sersic 1963, for a review), R1/n is described by:

\begin{displaymath}I(R) = I_{\rm e} ~ {\rm exp}\left[ -b_n \left( \left( \frac{R}{R_{\rm e}}\right)^{1/n} - 1 \right) \right]
\end{displaymath} (1)

where $I_{\rm e}$ is the intensity at the effective radius, $R_{\rm e}$ (distance that encloses half of the total emission), and bn is a constant coupled to the value of n. Special cases of the Sersic profile are those where n=4 (de Vaucouleurs profile), n=1 (exponential profile) and n=0.5 (Gaussian profile).

The TinyTim (Space Telescope Science Institute program to generate simulated HST point spread functions) models for the HRC PSF are not sufficient for our purposes. Instead, we used the calibration plan observations of Cycles 12 and 13 (programs 10054 and 10374). These programs contain observations of the spectrophotometric standard stars GD71, G191B2B, GD153 and HZ 14. To model our PSF, we compared the PSF created by each of these stars in each Cycle, as well as combinations of these. In general, the differences between each PSF are subtle and do not produce major differences in the models of the sources. However, the average PSF of the spectrophotometric stars of Cycle 13 produced better point source results (lower $\chi^2$ and residuals) so we chose to use it in our models. This PSF has a FWHM of $\sim$1.9 pixels, implying a resolution of $\sim$ $0.05\hbox{$^{\prime\prime}$ }$. However, the residuals for 1814-637 suggest that an even better PSF model is required (see Fig. 1).

GALFIT yields the coordinates of the modeled components in pixels. The conversion to RA and Dec was performed using the astrometry information stored in the header of each HST/ACS - pipeline reduced - image (model errors are listed in Table 2). The radio positions of the sources are from the literature. However, the registration between HST and radio frames can be off by as much as 1 $\hbox{$^{\prime\prime}$ }$.

Unless otherwise noticed, all presented fluxes, luminosities and magnitudes are for the F330W filter passband, corrected for galactic extinction and k-correction. The magnitudes are in the STMAG system. Galactic extinction was corrected using the Galactic de-reddening curve in Cardelli et al. (1989) and the measured Galactic extinction values of Schlegel et al. (1998). K-correction (typically $\sim$0.1, 0.2 mag) and conversion of the Allen et al. (2002) filters were done following the PEGASE (Fioc & Rocca-Volmerange 1997) templates and using the IRAF package SYNPHOT. SYNPHOT was also used for measuring F330W magnitudes of the PEGASE and GALAXEV (Bruzual & Charlot 2003) stellar population models. We use H0=71, $\Omega_{\rm M} = 0.27 , \Omega_\Lambda = 0.73$ (Spergel et al. 2003) throughout the paper.

Table 1: Redshift and wavelengths.

3 Nuclear and emission line contamination

Here we discuss several possible contributions to the near-UV light in these radio galaxies. Using ground-based, optical spectroscopy, Tadhunter et al. (2002) studied the nature of the UV excess in GPS, CSS and FR II sources at redshifts 0.15 < z < 0.7. They found that the UV continuum in these sources has contributions from (1) nebular continuum, (2) direct AGN light, (3) scattered AGN light, and a (4) starburst component. In addition to these continuum components we also consider contamination from emission line gas.

The observed objects are Narrow Line Radio Galaxies (NLRG) so we expect no contamination from direct light from the AGN. The main contribution from emission line gas to our observations would come from Mg II (see Table 1). It is usually only found in the nuclear Broad Line Region (BLR) of AGN hosts so we do not expect direct contamination from emission line gas.

To make a rough estimation of the extent of the possible contamination by direct light from the AGN, we use the STIS spectrum of 3C 277.1 (Labiano et al. 2005), a CSS QSO that could represent our worst case scenario: when we are looking into the nucleus. In 3C 277.1, all signs from nuclear contribution (broad lines and AGN continuum) and Mg II emission disappear at 0.8 kpc from the nucleus, roughly 1.5 times the FWHM of the PSF in the STIS observations. On average, 0.8 kpc correspond to $\sim$ $0.2\hbox{$^{\prime\prime}$ }$ in our ACS observations and we could expect the nuclear traces to disappear closer to the center in our ACS galaxies. The C II] line is usually fainter than Mg II (Peterson 1997) in AGN (it is also not present in 3C 277.1) so it is unlikely that it is affecting our observations.

There are UV polarization observations of only two objects, 1934-638 and 1345+125, so the presence of scattered nuclear light in the sample cannot be completely ruled out. Tadhunter et al. (2002) find AGN scattered light in $\sim$$37\%$ of their sample (including the GPS 1934-638), but in most cases it does not seem to dominate the UV emission. As we note below in Sect. 5.1, the UV light in the GPS sources tends not to be aligned with the radio source suggesting that it is not due to scattered nuclear light (assuming the scattered photons would escape along the radio axis). Observations of 1345+125 by Hurt et al. (1999) are consistent with the existence of polarized UV light (though with large uncertainty). However, Tadhunter et al. (2005) found evidence of recent star formation and strong jet cloud interactions in 1345+125 (see also Holt et al. 2003; Surace et al. 1998).

In general, the contribution from nebular continuum in radio galaxies varies between 3 and 40$\%$ (Tadhunter et al. 2002). However, 3 out of the 4 CSS they studied are dominated by young stellar populations. Furthermore, for up to 50$\%$ of their sources (including GPS, CSS and FR II), the UV excess is dominated by young stellar populations. Similar results have been found for FR I and FR II sources at a range of redshifts (Aretxaga et al. 2001; Wills et al. 2002,2004).

To sum up, the most probable sources of UV light in our images are star formation, scattered nuclear light and nebular continuum. However, it is likely that the contributions of the latter two are small.

   
4 UV morphology

Table 2: Galfit components.

We have modeled the UV emitting regions with GALFIT and found that 1117+146 and 1814-637 show one point source component. There is a marginal detection of a point source in 1233+418 and 1946+708. The hosts of 1345+125, 1443+77 and 1934-638 show a combination of at least one Sersic[*] component (with different indices) and one or several point sources (see Table 2 and Fig. 1). 1607+268 and 2352+495 show a combination of two and three point sources respectively. In the near-IR, de Vries et al. (2000) find Sersic indices $2 \la n \la 5$ and effective radii $2~{\rm kpc} \la R_{\rm e} \la 4$ kpc for GPS and CSS sources. Most of our data show point sources and/or Sersic profiles with indices $n \la 2$ and radii $R_{\rm e}\la 400$ parsec. The presence of these small clumps of near-UV emission is consistent with star forming regions in the host. Before addressing the nature of the UV emission, we describe the properties of this emission in relation to other properties, for the individual sample sources. We note that optical observations of the higher redshift Parkes half-Jansky sample suggest that the rest frame UV may contain a contribution from a young stellar population (de Vries et al. 2007).

For our observations, the 3$\sigma$ detection limit, for a point source ( ${\it FWHM} \sim3$ pixels) is 25.7 (no Galactic extinction applied), calculated from the noise in the images. The ACS Imaging Exposure Time Calculator gives the following 3$\sigma$ detection limits for our observations: Johnson U = 21/arcsec2 for an extended uniform source ($2\times 2$ pixel extraction region) and 25.2 for a point source[*]. For comparison purposes with this limit and published data, the magnitudes listed in this section and Table 2 have not been corrected for Galactic extinction.

   
4.1 Notes on individual sources

1117+146: Identified as a GPS by O'Dea et al. (1991). The counterpart of the radio source corresponds to a $m_{\rm R}=20.1$ galaxy (de Vries et al. 1995) at z = 0.362 (de Vries et al. 1998b). Radio observations (e.g., Fey & Charlot 1997) show a $\sim$100 mas double radio source. Our image shows an unresolved 24.06 mag source.

1233+418: A CSS galaxy with a photometric R-band redshift $z_{\rm R} = 0.25$ (Fanti et al. 2004). We have a marginal detection of a point source with magnitude 24.93. Our image (see Fig. 1, first panel) suggests faint extended emission towards Northeast but the errors in the GALFIT model for this component are too large to be sure of its existence. Checking the individual pre-drizzle images, it is not clear if these components are real or artifacts.

1345+125: 4C +12.50. A long known peaked spectrum radio source (e.g., Véron 1971). The counterpart of this GPS source is a $m_{\rm R}=15.5$ galaxy (e.g., Stanghellini et al. 1993) at z = 0.12174 (Holt et al. 2003), in a cluster of fainter galaxies (Stanghellini et al. 1993). IR images of this well known ULIRG show an extremely reddened source with two nuclei separated $\sim$ $1.8\hbox{$^{\prime\prime}$ }$ ($\sim$4 kpc) embedded in a common envelope and aligned roughly East-West (e.g., Surace & Sanders 2000; Scoville et al. 2000) suggesting an ongoing merger, which may have triggered the AGN (e.g., Xiang et al. 2002; Heckman et al. 1986). The western nucleus is the brightest and shows a Seyfert 2 spectrum (Gilmore & Shaw 1986). Veilleux et al. (1997) suggest that the source may have a hidden quasar, also supported by UV polarized continuum emission ( $p=16.4\% \pm2.6\%$, the polarization vector approximately perpendicular to the radio axis Hurt et al. 1999). VLBI imaging shows a complex, distorted $\sim$100 mas ($\sim$0.2 kpc) source (e.g., Lister et al. 2003), roughly oriented North-South ( ${\rm PA}\sim-20\deg$). Older observations of 1345+125 related the radio source to the East nucleus but improved astrometry showed that it is related to the western one (e.g. Fanti 2000; Axon et al. 2000; Stanghellini et al. 1997). Evans et al. (1999) study the molecular gas in 1345+125 and suggest that the molecular gas is fueling the AGN. Our UV image detects both nuclei separated by $\sim$ $1.8\hbox{$^{\prime\prime}$ }$. The East component was modeled with a Sersic profile with index $\sim$2.3 and an effective radius $\sim$ $ 0.3\hbox{$^{\prime\prime}$ }$ and magnitude 21.84. The West nucleus shows a more complex structure: an extended component with magnitude 21.20 and Sersic index 1.62 and a point source with magnitude 21.49. This more complex structure could be due to interaction with the radio source, the Eastern component or a combination of both. Optical emission line images (Batcheldor et al. 2006; Axon et al. 2000) show an arc of emission $\sim$ $1\hbox{$^{\prime\prime}$ }$ North of the Western nucleus and fainter emission at $2\hbox{$^{\prime\prime}$ }$. We do not detect UV emission associated with these emission line features. They also detect a faint tail of emission stretching from the West nucleus towards the west, present in our image and not modeled by GALFIT (see Fig. 1). Surace et al. (1998) detect compact blue knots in 1345+125 around the source which they attribute to star forming regions. We detect the southern knots $\sim$ $2\hbox{$^{\prime\prime}$ }$ south of the nuclei but not the northern one (see Fig. 1).

1443+77: 3C 303.1. The optical counterpart of this CSS corresponds to a galaxy of $m_v\sim20$ (e.g., Sanghera et al. 1995), at z = 0.267 (Kristian et al. 1978). The MERLIN map (Sanghera et al. 1995) shows a $\sim$ $1.8\hbox{$^{\prime\prime}$ }$ long ($\sim$7 kpc) double radio source aligned NW-SE ( ${\rm PA}\sim 47\deg$) aligned with the inner emission line gas (e.g. de Vries et al. 1999). XMM-Newton observations detect the ISM of the host galaxy as well as a second component which could be either Synchrotron Self Compton from the Southern radio lobe or hot gas shocked by the expansion of the radio source (O'Dea et al. 2006). The center of 1443+77 shows a complex structure in the optical, which could be due to an ongoing merger (e.g., Axon et al. 2000; de Vries et al. 1997) and it also shows up in our UV image. Axon et al. (2000) and de Vries et al. (1999) find an arc of emission south of the nucleus of the source and McCarthy et al. (1995) detects circumnuclear [O III] as far as 3 $\hbox{$^{\prime\prime}$ }$ from the center. Our image suggests the presence of an arc of emission $\sim$ $0.9\hbox{$^{\prime\prime}$ }$ ($\sim$4 kpc) from the nucleus. The integrated magnitude of this arc is $\sim$23. A region of the sky with the same area has a magnitude $\sim$25. The GALFIT model of the source consists of two point sources of magnitudes 25.21 (marginal detection) and 24.06, and an extended component with Sersic index 2.66 and effective radius 1.1 $\hbox{$^{\prime\prime}$ }$ ($\sim$4.5 kpc). The arc is not modeled by GALFIT so it is visible in the residuals of Fig. 1. 1443+77 is the best candidate in our sample to be undergoing jet induced star formation (see below).

1607+268: Also known as CTD93. VLBI maps (e.g. Dallacasa et al. 1998) show a two-component, $\sim$60 mas, GPS source. The counterpart is a galaxy at z = 0.473 (O'Dea et al. 1991) with $m_{\rm r} = 20.4$ (Stanghellini et al. 1993). We observe two point sources of magnitudes 24.38 and 23.92.

1814-637: VLBI imaging shows a two component CSS galaxy (e.g., Tzioumis et al. 2002). Although the radio source is small ($\sim$ $0.4\hbox{$^{\prime\prime}$ }$, $\sim$0.3 kpc) its radio spectrum does not peak at $\sim$1 GHz. The optical identification corresponds to a mv=18.0 galaxy at z = 0.063 (Wall & Peacock 1985). We observe a bright, 16.38 mag point source. Danziger & Goss (1979) describe 1814-637 as a narrow line radio galaxy with a bright star-like nucleus. They report H and H absorption lines at zero redshift and suggest that there is a foreground star superposed on the galaxy. The UV object is $\sim$2 arcsec to the Northeast of the radio source (Ma et al. 1998). However, the registration between HST and radio frames can be off by as much as 1 $\hbox{$^{\prime\prime}$ }$, so it is probably the counterpart.

1934-638: A long known GPS (Bolton et al. 1963). VLBI maps (e.g., Tzioumis et al. 1998) show two components separated by $\sim$40 mas (0.12 kpc at z = 0.183, Tadhunter et al. 1993). R-band observations (Jauncey et al. 1986) show a system of two galaxies, separated by $\sim$ $3\hbox{$^{\prime\prime}$ }$, consistent with our observations. The GALFIT model yields two extended components with magnitudes 21.08 and 21.75 and a fainter point source with magnitude 23.81. 1934-638 shows significant polarization in the UV (3.5$\%$, e.g., Morganti et al. 1997; Tadhunter et al. 1994), with the position angle of the electric vector perpendicular to the radio axis. Scattered AGN light probably makes a significant contribution to the UV excess in this source (Tadhunter et al. 2002). Optical polarization observations also suggest an anomalous environment (Morganti et al. 1997; Tadhunter et al. 1994).

1946+708: The counterpart is identified with a z = 0.10083, $m_{\rm R}=16.3$ galaxy (e.g., Snellen et al. 2003). VLBI observations of this GPS source show an elongated, $\sim$40 mas ($\sim$0.07 kpc) NE-SW structure (Taylor & Vermeulen 1997). Optical observations (Perlman et al. 2001) also suggest a NE-SW elongated source. Our near-UV image suggests a marginal detection of a point source with magnitude 24.53, present in two of the three individual pre-drizzle images. Perlman et al. (2001) suggest that 1946+708 may be part of a group of galaxies and the closest would be $\sim$ $1\hbox{$^\prime$ }$ far from 1946+708. These objects fall outside our field and they have unknown redshifts.

2352+495: Snellen et al. (2003) identify the counterpart of this GPS source (e.g., O'Dea 1998) with a z = 0.23790, $m_{\rm R} = 18.2$ galaxy. VLBI observations show a complex elongated $\sim$ $0.07\hbox{$^{\prime\prime}$ }$ ($\sim$0.25 kpc) source oriented NW-SE (Pollack et al. 2003). Our ACS observations show three point sources in a circle about 0.5 $\hbox{$^{\prime\prime}$ }$ wide, with magnitudes 24.79, 25.13 and 25.53 (slightly above the detection limit). The individual pre-drizzle show the two brightest sources while the third one is dubious.

5 UV and radio properties

Here we compare the UV and radio properties of our sample with those of the sample of 3CR FR I and FR II sources studied by Allen et al. (2002). The Allen et al. (2002) data consist of STIS near-UV-MAMA snapshots (exposure time of 1440 s) with filters F25SRF2 $\lambda_{\rm c} = 2270~\AA$) and F25CN182 ( $\lambda_{\rm c} = 1820~\AA$). The radio data have been collected from the papers listed in the individual notes above: O'Dea (1998), Morganti et al. (1993), Martel et al. (1999), Flesch & Hardcastle (2004) and the 3CRR on-line catalog (an update of the sample of Laing et al. 1983). The data are summarized in Table 3 and plotted in Figs. 2 to 10.

Table 3: UV and 5GHz radio properties.

We note that O'Dea & Baum (1997) found that the radio power of the bright GPS and CSS sample is independent of size and comparable to that of the most powerful FR II sources. However, the Allen et al. (2002) FR sources are low z 3CR, and therefore lower radio luminosity sources.

5.1 GPS sources

Most of the sources in our sample are very compact GPS sources, and therefore, are probably too small to strongly affect their environment on the scales resolved by these observations. However, some conclusions can be drawn from their integrated UV and radio properties.

Figure 2 suggests a trend between UV luminosity and radio power of GPS sources (consistent with the results of Raimann et al. 2005, for larger radio galaxies). This relationship could be produced by a larger reservoir of gas available in the more luminous sources for fueling the radio source and supplying the starburst.

Inspection of Fig. 3 shows that alignment between the UV and radio is found for sources $\ga$1 kpc, i.e., the compact GPS sources tend not to show systematic alignment. If the UV light was scattered nuclear light, it would be expected to align with the radio axis regardless of radio source size (assuming the scattered light would escape along the radio axis). The lack of alignment in the GPS sources suggests that the UV light is not scattered nuclear light. Instead, in the GPS sources, we may be detecting some clumps of star formation on scales larger than the radio source which are associated with the fueling of the radio activity.

5.2 CSS sources

The fact that the GPS sources show UV emission on larger scales than the radio source and which are not aligned with the radio source is consistent with the hypothesis that the GPS sources are too small to strongly affect their environment on the scales resolved by these observations. In large scale 3CR sources, the jets extend far beyond the host galaxy. Therefore, if there is jet induced star formation, CSS sources are the best candidates to reveal it. We also note that the two most luminous UV sources in our sample (1443+77 and 1814-637) are CSS.

Labiano et al. (2005) demonstrated the presence of gas ionized by the shocks from the expanding radio source in CSS sources. Furthermore, they found that 1443+77 shows the strongest contribution from shocks. It is possible these shocks are affecting the star formation in the host.

5.3 Is there jet-induced star formation?

Figure 4 suggests there is no correlation between UV luminosity and size of the radio source (i.e., the GPS, CSS and large 3CR sources have similar UV luminosities). Therefore, any UV emission which is related to source size does not seem to dominate the UV properties of the host. However, if the expansion of the radio lobes were enhancing the UV luminosity, it could have been unnoticed in our sample. The lack of a significant number of CSS sources in our sample (i.e. radio sizes comparable to the host) may cause us to be missing hosts with UV emission $\sim$ $10^9~L_{\odot}$.

If the jet is enhancing the star formation, we expect the radio source and UV emission to be aligned and have similar sizes. The high UV luminosity source 1443+77 shows alignment between the UV and radio source (Fig. 5). Observations of additional sources with sizes $\ga$1 kpc are needed to improve the statistics. Figure 6 shows that UV size and radio size are not correlated for most sources. However, two CSS and one GPS source (1443+77, 1814-637 and 2352+495) have similar UV and radio sizes. The radio source in 1233+418 is larger than the (marginal) UV emission. However, the UV magnitudes for 1233+418 are close to the detection limit so we may be missing extended and fainter UV emission in this object.

Table 3 shows the position angles of the radio and UV in the GPS and CSS sources. We see that the CSS 1443+77 and the GPS 2352+495[*] are aligned with the UV. In addition 1443+77 1814-637 and 2352+495 have a ratio of radio to UV size of order unity (Fig. 7). The UV emission in 2352+495 has two components aligned with the radio source (and a third which is not aligned), but the unknown offset between radio and HST reference frames prevents us from accurately overlaying them. Therefore, 1443+77 is currently the best candidate for jet induced star formation.

This is consistent with the hypothesis that the UV emission is dominated by extended regions of star formation rather than by point-like AGN.

We also note several interesting null results. The detected UV luminosity in the GPS, CSS, and large 3CR sources is independent of redshift (out to z=0.5) (Fig. 8). In the GPS and CSS sources, the alignment does not depend on either radio power (Fig. 9) or redshift (Fig. 10); though in both cases the statistics are dominated by the GPS sources which seem not to show alignment on the scales resolved by these observations.


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{figs/7112fg2.eps}
\end{figure} Figure 2: UV luminosity and radio power of the combined sample of GPS/CSS and FR I, FR II from Allen et al. (2002). The solid line shows the correlation between UV luminosity and radio power for GPS sources.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{figs/7112fg3.eps}
\end{figure} Figure 3: Difference between the UV and radio position angles (alignment) of the GPS and CSS galaxies versus radio size.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{figs/7112fg4.eps}
\end{figure} Figure 4: As Fig. 2, for radio size.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{figs/7112fg5.eps}
\end{figure} Figure 5: As Fig. 3, for UV luminosity.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{figs/7112fg6.eps}
\end{figure} Figure 6: UV size and radio size of GPS and CSS sources. The dashed line marks the locus of sources with equal radio and UV sizes.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{figs/7112fg7.eps}
\end{figure} Figure 7: Ratio of radio and UV sizes versus the alignment of the GPS and CSS sources. The dashed line marks the locus for sources with same UV and radio size and highest alignment i.e. where the radio source is most likely affecting the star formation in the host.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{figs/7112fg8.eps}
\end{figure} Figure 8: As Fig. 2, for redshift.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{figs/7112fg9.eps}
\end{figure} Figure 9: As Fig. 3, for radio power.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{figs/7112fg10.eps}
\end{figure} Figure 10: As Fig. 3, for redshift.
Open with DEXTER

Table 4: Mass from luminosity.

   
6 Stellar synthesis models

GPS and CSS sources are usually associated with massive ellipticals of ages $\sim$5 Gyr and solar metallicities (e.g., de Vries et al. 1998a; Snellen et al. 1998; de Vries et al. 2000; Snellen et al. 1999; de Vries et al. 2007). With this assumptions, we use stellar population synthesis models by Bruzual & Charlot (2003) using the Chabrier (Chabrier 2003) initial mass function and Padova evolutionary tracks, for a 1 Myr long continuous burst of star formation and a instantaneous single burst model for populations of 104, 103, 10 and 1 Myr to compare our UV luminosity measurements and to estimate the mass and age of stars producing the UV emission in our sources. We also compare the measurements with the elliptical galaxy templates from PEGASE.

Table 4 shows the mass of stars for different ages needed to reproduce the observed near-UV emission of each source. Our observations are generally consistent with models of a single instantaneous burst of 10 $^6{-}10^7~M_\odot$, $\la$10 Myr ago. However, some sources may be dominated by intermediate-age (0.1 to 1 Gyr) populations (with masses $\sim$ $10^9~M_\odot$, e.g., Tadhunter et al. 2005). These young and intermediate ages in our GPS/CSS sources are also consistent with stellar population ages measured in hosts of powerful radio galaxies (Raimann et al. 2005) and other CSS sources (e.g., Johnston et al. 2005). The models which require old and massive stellar populations to produce the UV light are inconsistent with the small measured sizes for the UV emitting regions.

Our observations are consistent with a scenario in which a single event is responsible for triggering the AGN and initiating a starburst. However, the dynamical ages of GPS and CSS sources are much smaller (between 103 and 106, e.g., Polatidis & Conway 2003; O'Dea 1998) than the estimated ages of the stellar populations produced in the starburst. This implies a minimum time delay between the start of the starburst and the start of the radio activity of 10 Myr. The delay could be as long as 1 Gyr if the 1 Gyr stellar population models are correct. The apparent delay between starburst and formation of the radio source has been found in objects with possible connections between the starburst and the formation of the AGN and radio source (Raimann et al. 2005; Tadhunter et al. 2005) and is predicted by theoretical calculations of the time needed for the gas to reach the center of the galaxy after a tidal interaction (e.g., Lin et al. 1988).

7 Summary

We have obtained HST/ACS/HRC near-UV high resolution images of compact GPS and CSS radio galaxies which are likely to be the young progenitors of large scale powerful radio galaxies. We detect near-UV emission in point sources and/or small clumps (tens to hundreds of pc scale) in seven of the sources, consistent with the presence of recent star formation. Two sources show very weak near-UV emission and may have not been detected. The UV luminosity in the GPS and CSS sources is similar to that measured in nearby large 3CR radio galaxies by Allen et al. (2002). In the GPS sources as a whole we do not see systematic alignment between the radio and UV. This may be because the GPS radio sources are smaller than the scales resolved by these HST observations. The lack of alignment in the GPS sources implies that the UV emission is not due to scattered nuclear light.

There is evidence for a correlation between radio and UV luminosity in the GPS sources. This relationship could be produced by having a larger reservoir of gas available in the more luminous sources for fueling the radio source and supplying the starburst.

In the CSS source 1443+77 the near-UV emission is aligned with and co-spatial with the the radio emission and we suggest that star formation has been triggered or at least enhanced by expansion of the radio source through the host.

We suggest that the starburst and AGN are triggered by the same event. Comparison with stellar population synthesis models suggests the data are consistent with star formation occurring in a burst $\la$10 Myr ago (though models as old as 1 Gyr are not ruled out). The radio source ages are much smaller (e.g., 103 to 106 yr) suggesting that there is a minimum delay of 10 Myr between the onset of the starburst and the onset of the radio activity (see also Raimann et al. 2005; Tadhunter et al. 2005). Observations at other wavelengths and measurement of the colors are needed to further asses the nature of the observed UV properties.

Acknowledgements
We thank the anonymous referee for valuable suggestions and comments on the paper. A.L. wishes to thank Dr. Peng (STScI) for his help with GALFIT and Dr. I. A. G. Snellen (Leiden) for fruitful scientific discussions. W.D.V.s work was performed under the auspices of the US Department of Energy, National Nuclear Security Administration by the University of California, Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48. These observations are associated with program 10117. Support for program 10117 was provided by NASA through a grant from the Space Telescope Science Institute. This research has made use of NASA's Astrophysics Data System Bibliographic Services and 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. This research is based on observations made with the NASA/ESA Hubble Space Telescope, obtained fro m the data archive at the Space Telescope Institute. STScI is operated by the association of Universities for Research in Astronomy, Inc. under the NASA contract NAS 5-26555.

References

 

Copyright ESO 2007