A&A 418, 429-443 (2004)
DOI: 10.1051/0004-6361:20034486

Further clues to the nature of composite LINER/H II galaxies

M. E. Filho 1,2,3 - F. Fraternali 1,4 - S. Markoff5,[*] - N. M. Nagar 1 - P. D. Barthel 1 - L. C. Ho 6 - F. Yuan 7


1 - Kapteyn Astronomical Institute, PO Box 800, 9700 AV Groningen, The Netherlands
2 - Istituto di Radioastronomia, CNR, via P. Gobetti, 101, 40129 Bologna, Italy
3 - Centro de Astrofísica da Universidade do Porto, Rua das Estrelas, 4150 - 762 Porto, Portugal
4 - ASTRON, PO Box 2, 7990 AA Dwingeloo, The Netherlands
5 - Massachusetts Institute of Technology, Center for Space Research, 77 Massachusetts Av., R. NE80 - 6035, Cambridge, Massachusetts 02139, USA
6 - The Observatories of the Carnegie Institution of Washington, 813 Santa Barbara Street, Pasadena, California 91101, USA
7 - Purdue University, Department of Physics, 525 Northwestern Av., West Lafayette, Indiana, 47907 - 2036, USA

Received 9 October 2003 / Accepted 14 January 2004

Abstract
We have analyzed new, archival and published high resolution radio and X-ray observations of a sample of composite LINER/H  II galaxies known to exhibit AGN-like properties. Five of the 16 AGN candidates have milliarcsecond-scale detections and are found to display a compact, flat spectrum, high brightness temperature radio core, four of which also exhibit extended radio emission. Five of the eight AGN candidates with available high resolution X-ray observations were found to possess a hard X-ray nuclear source, two of which have no milliarcsecond radio detection. The combined high resolution radio and X-ray data yield a 50% detection rate of low luminosity AGN among the AGN candidates, which translates into a 12% detection rate for the entire composite LINER/H  II sample. In the sources where the AGN has been unambiguously detected, the ionizing power of the AGN is not sufficient to generate the observed emission lines, unless the hard X-rays are heavily obscured. We attempt to apply a canonical advection-dominated accretion flow (ADAF) and jet model to the sample sources in order to explain the observed radio and X-ray emission. While ADAFs may be responsible for the observed emission in submillijansky radio cores like NGC 7331, they do not appear consistent with the radio emission observed in the milliarcsecond-scale radio detected cores; the latter sources are more likely to have an energetically important contribution from a radio-emitting jet.

Key words: galaxies: active - galaxies: nuclei - radiation mechanisms: general

1 Introduction

The radio emission observed in nearby galaxies can be attributed to stellar and interstellar processes - non-thermal synchrotron radiation from electrons accelerated by supernova remnants (SNRs) and thermal emission from H  II regions - or to the presence of accreting black holes. In the former case, the radio emission coincides with the infrared and optical emission-line sources, which are themselves tracers of young stars and supernovae. The radio emission originates in diffuse, low surface brightness regions that do not extend beyond the stellar light distribution. On the other hand, radio sources produced by accreting black holes are generally characterized by compact, nuclear radio cores, sometimes associated with jets and lobes. In many instances circumnuclear starburst regions and black holes coexist spatially, making it difficult to disentangle the two components in the optical regime and also in the radio. Furthermore, at high resolution, even star-forming galaxies may show compact radio sources, associated with radio supernovae (RSN) and SNRs (e.g. Weiler et al. 1986). In an active galactic nucleus (AGN), the presence of a flat radio spectrum source ( $\alpha\la0.5$; $F_{\nu}\propto\nu^{-\alpha}$) indicates synchrotron self-absorption in components whose brightness temperatures are comparable to the kinetic temperatures (mc2/k; $T_{B}\ga10^5$ K) of the relativistic electrons producing the synchrotron radiation. Steep spectrum emission in star-forming galaxies ( $\alpha\ga0.5$) arises from the non-thermal emission of SNRs, whereas flat spectra are consistent with free-free absorption of optically thin synchrotron radiation and/or thermal emission ( $T_{B}\la10^4$ K) originating from the electrons in H  II regions. So, while the spectra of thermal emitters may mimic those of partially opaque synchrotron sources, a combination of brightness temperature and radio spectral index of the strong nuclear radio components can indicate whether the emission is thermal ( $T_{B}\la10^4$ K; Condon 1992) or non-thermal ( TB>105 K) and whether optically thin thermal ( $\alpha\sim0.1$), self-absorbed synchrotron ( $\alpha<0.5$), or optically thin synchrotron ( $\alpha\sim0.5$-0.8) radiation dominates.

We have therefore applied this radio technique to a sample of mildly active galaxies - composite LINER/H  II galaxies - known to exhibit both starburst and AGN characteristics in their optical spectra. The motivation for studying this class of galaxies is to study the demographics of black hole accretion in cases suspected of weak AGN activity. This is an important issue because, if composite sources do indeed harbour a genuine active nucleus, then their considerable number (13% of emission-line nuclei with BT< 12.5 mag; Ho et al. 1997) will imply a significant contribution to the AGN population.

2 Sample selection

In Filho et al. (2000, 2002) we have presented a study of the radio characteristics of a complete sample of composite LINER/H  II galaxies, i.e., galaxies that are hypothesized to be LINERs spatially contaminated by circumnuclear star-forming regions (Ho 1996). The first part of our project consisted in gathering arcsecond-scale radio information (published or new) of the 65 composite sources contained in the magnitude-limited Palomar survey (Ho et al. 1995, 1997a). Their radio properties (Filho et al. 2000, 2002) indicate composite sources come in two types: AGN-like sources with compact cores and starburst-like sources dominated by diffuse emission, co-spatial with the projected galactic disk. The complete composite source sample contained 14 ($\sim$25%) AGN candidates as revealed by compact, flat/inverted radio spectrum cores on arcsecond scales, with peak flux densities above 1 mJy. NGC 660 and NGC 7331 have been included in the AGN candidate sample although they have sub-mJy radio peaks as explained in Filho et al. (2002). These sixteen sources comprise the present sample of candidate AGN-driven composite LINER/H  II galaxies.

3 New and published radio observations

Radio imaging at sub-arcsec and mas resolution of several of the sample sources have already been published. Our observations, using the Very Large Array (VLA) and the Very Large Baseline Array (VLBA), complement these observations and yield a complete picture of the occurrence of compact nuclei in the sample.

In this section we report on high resolution VLA observations of 12 (four new and eight published) and multiwavelength VLBA observations of 15 (13 new and two published) out of the 16 sample sources. Regarding the former, of the possible eight galaxies with no published sub-arcsec-scale radio information, four were chosen to fit the observational window assigned to us by the VLA. As for the latter, NGC 7331 was not observed with the VLBA because its radio core flux density at 1 $.\!\!^{\prime\prime}$5 resolution is well below 1 mJy at both 1.4 and 5 GHz (Cowan et al. 1994) and NGC 5866 and NGC 4552 because they had been previously observed and detected with the VLBA by Falcke et al. (2000) and Nagar et al. (2002), respectively. Following a description of our observations, we will combine the new and published data in a coherent analysis of the sample.

3.1 Observations and data reduction

VLA observations of four sources in the sample were performed on 1999 September 5 with the A-array, X-band system (8.4 GHz). Two IFs of 50 MHz each were combined to give a total bandwidth of 100 MHz. In its A-array configuration, the VLA yields typical resolutions of 0 $.\!\!^{\prime\prime}$25 at 8.4 GHz (which corresponds to 12 pc at D=10 Mpc). Observations of the sample sources were interspersed with observations of nearby phase calibrators, for a total integration time of about 10 min per galaxy. The primary flux calibrator was 3C 286 (1328+307) with adopted 8.4 GHz flux densities of 5.1915 Jy and 5.1702 Jy for IF1 and IF2, respectively. The uncertainty associated with the flux calibration procedure is mainly due to the uncertainty in the absolute flux density of 3C 286, which is conservatively set to 5%. Antenna gains were found to behave well throughout the observations.

VLBA observations of 13 out of the 16 sample galaxies were obtained in five sessions - 2000 June 22, 2001 September 1, 2001 September 8, 2001 September 17 and 2001 October 6 - with the standard ten VLBA antennas. The standard observing frequency was 5 GHz, although some sources have multi-wavelength observations. The four 8 MHz IFs yielded a total bandwidth of 32 MHz. Typical resolutions for the VLBA at 5 GHz are 2 mas (which corresponds to 0.1 pc at D=10 Mpc). Phase referencing was performed, which includes several minute on-source scans alternated with scans of nearby phase calibrators, yielding generally 180 (in some cases 45) min total integration time on each target source.

Reduction of both the VLA and VLBA data was performed using standard NRAO AIPS (version 15Oct99) image processing routines. After initial deleting of bad data points and calibration, the data were "cleaned''. The AIPS "cleaning'' task IMAGR was employed to Fourier transform the data and remove sidelobes from the maps in an interactive and iterative mode. Some sources were observed in multiple VLBA runs separated by several days. For these sources, maps were made for each run separately, since combining the data did not result in a better map. Full resolution images, having a 0 $.\!\!^{\prime\prime}$25 synthesized beam, were obtained for the VLA data and 2 mas synthesized beam at 5 GHz for the VLBA data. The image noise level reached the theoretical level of $\sim$0.07 mJy/beam (uniform weighting; Perley et al. 1989) to within a factor of two for the VLA data and $\sim$0.05 mJy/beam for the VLBA data. Phase self-calibration was employed on the strongest sources, leading to some reduction in the image noise level. The 8.4 GHz VLA calibrated data provided by J. Ulvestad for NGC 660 have been re-analysed, yielding slightly different radio parameters than those presented in Filho et al. (2002).

The image noise levels were measured with the AIPS task IMSTAT in a source-free region. Using the AIPS task IMFIT, the brightness peaks of the radio sources were fitted with single bi-dimensional Gaussians.

3.2 Results of the radio observations

Table 1: Map parameters (published or new) of the sample sources. Column 1: source name. Column 2: observing frequency. Column 3: integration time. Column 4: rms noise level of the image. Column 5: restoring beam or resolution. Column 6: position angle of the beam. Column 7: peak radio flux density. Column 8: integrated radio flux density. Columns 9 and 10: radio position. Column 11: deconvolved source size. Column 12: position angle of the source. Value in parenthesis is an uncertain value. Column 13: note. Letter is the observing epoch and number is the reference (published or new) from which the radio data were taken.


  \begin{figure}
\par\mbox{\hspace*{3.4cm}\subfigure[]{\includegraphics[width=6.3c...
...}
\subfigure[]{\includegraphics[width=6.4cm,clip]{0486fig1c.eps} }}
\end{figure} Figure 1: Radio emission contours with the following multiples of the rms noise (Table 1): (-3, -2.1, 2.1, 3, 4.2, 6, 8.4, 12, 16.8, 24). The size of the restoring beam in milliarcseconds, observing epoch, integration time in minutes, and frequency in GHz are given in parentheses after each object name (see Table 1). a) NGC 524 ( $4.5 \times 2.0$; 2001 Sep. 01; 180; 5); b) NGC 5354 ( $7.6 \times 3.8$; 2001 Sep. 17; 40; 2.3); c) NGC 5354 ( $2.3 \times 1.6$; 200 June 22; 105; 5).
Open with DEXTER


  \begin{figure}
\par\subfigure[]{\includegraphics[width=6.3cm,clip]{0486fig2a.eps...
...par\subfigure[]{\includegraphics[width=6.3cm,clip]{0486fig2c.eps} }
\end{figure} Figure 2: Same as Fig. 1. a) NGC 5846 (9.2 $\times $ 3.7; 2001 Oct. 06; 48; 2.3); b) NGC 5846 (3.1 $\times $ 1.4; 2000 June 22; 112; 5); c) NGC 5846 ( $1.7 \times 0.7$; 2001 Oct. 06; 48; 15).
Open with DEXTER

Four of four VLA and three of 13 newly observed VLBA sample sources were detected. Table 1 lists image parameters for the new data, where we have also included relevant published radio data from Nagar et al. (2000, 2002), Cowan et al. (1994; NGC 7331) and Falcke et al. (2000; NGC 5866). We refer to Filho et al. (2000, 2002) for an extensive discussion on the radio properties of the sample sources. Undetected sources have been given 5$\sigma$ upper limits for the VLA and VLBA data. We estimate the radio peak positions to be accurate to within $\la$0 $.\!\!^{\prime\prime}$2 for the VLA and $\la$1 mas for the VLBA data, due to the good phase solutions and accurate positions of the phase calibrators. Sources are considered unresolved if their deconvolved source sizes are less than half the beamwidth in any component.

Figures 1 and 2 show the highest signal-to-noise radio maps of the newly detected VLBA sources. The contour levels are the following multiples of the rms noise (Table 1) in the map: -3, -2.1, 2.1, 3, 4.2, 6, 8.4, 12, 16.8, 24.

3.2.1 Radio core properties

In Table 2 we list the radio properties of the detected mas-scale sample sources (from new or published data). Unless otherwise mentioned, we have used the longest integration time VLBA 5 GHz radio flux densities to calculate the spectral indices and brightness temperatures. Radio spectral indices should be taken with caution, due to non-simultaneity and mismatch in resolution. For the multiple radio components in NGC 5846, we have given 5$\sigma$ limits to any 5 GHz emission in order to obtain upper limits for the spectral indices. We calculate the brightness temperature as:

\begin{eqnarray*}T_{B}=\displaystyle 7.8\times10^6\left(\frac{F_{\nu}}{\rm
mJy}...
...{2.5~\rm mas}\right)^{-2}\left(\frac{\nu}{5~\rm GHz}\right)^{-2}
\end{eqnarray*}


where $F_{\nu}$ is the peak flux density at 5 GHz in the VLBA maps (Table 1) and $\theta$ is the FWHM of the Gaussian beam, taken to be typically 2.5 mas at 5 GHz. For unresolved sources, the quoted brightness temperatures should be considered lower limits.

Table 2: Nuclear radio properties (published or new). Column 1: source name. Column 2: high resolution radio flux density. Column 3: observing frequencies. Column 4: radio spectral index. Column 5: radio variability. Column 6: peak brightness temperature. Column 7: note.

In a distance-limited sample of low luminosity Seyferts, LINERs and composite LINER/H  II galaxies imaged with the VLA (150 mas resolution), 50% of the sources showed flat spectrum radio cores (Nagar et al. 2002, 2000). Subsequent VLBA 5 GHz observations of a flux-limited subsample (16 sources with flux densities >3 mJy) yielded a 100% detection rate of low luminosity AGN (LLAGN; Nagar et al. 2002; see also Falcke et al. 2000). In comparison, the present sample of AGN candidates yielded a VLA radio core detection rate of 83% (ten of 12 sources with VLA data, Table 1). Follow-up VLBA observations (including NGC 4552; Nagar et al. 2002 and NGC 5866; Falcke et al. 2000) have revealed that 33% (five of 15) of the sources have mas-scale radio cores. These results translate to an overall radio detection rate of 8% (five of 65) or 12% if we include NGC 660 and NGC 7331 (see Sect. 5), among the entire composite LINER/H  II sample. The mas-scale non-detections were sources with VLA flux densities <3 mJy; thus, our results are consistent with those of Falcke et al. (2000) and Nagar et al. (2002). Moreover, as has been pointed out in these papers and discussed also in Filho et al. (2002), scaling arguments suggest that the radio non-detections could well be lower power versions of the radio detections and thus harder to detect. The true AGN fraction in composite LINER/H  II galaxies is likely to be higher.

Of the four VLA observed sample sources, only NGC 660 was slightly resolved at VLA A-array resolution. The present new and published VLBA data show unresolved cores in NGC 524 (Fig. 1a), NGC 4552 (Nagar et al. 2002) and NGC 5354 (Figs. 1b,c). Some sources also show evidence of radio variability, as seems the case for NGC 4552 (see discussion in Filho et al. 2000). Furthermore, all VLBA detected sample sources (including NGC 4552; Nagar et al. 2002 and NGC 5866; Falcke et al. 2000) show flat spectrum radio cores, although the inverted spectrum in NGC 524 is most probably an effect of resolution mismatch. The origin of the mas-scale radio emission is unlikely to be free-free emission. The detected radio cores show brightness temperatures several orders of magnitude above the thermal limit (104 K; Condon 1992). Moreover, as argued in Falcke et al. (2001), a thermal origin for the radio emission would imply much higher X-ray fluxes than those observed (see Sect. 4). Therefore, the combination of brightness temperatures in excess of 107 K and spectral indices $\alpha\la0.5$ shows that the radio emission can be identified with (multi-component) non-thermal, self-absorbed synchrotron emission related to an accreting black hole.

On the other hand, the nature of the radio components associated with NGC 5846 is more complex. Our 2.3, 5, and 15 GHz VLBA maps of NGC 5846 show a multiple radio source roughly aligned in the North-South direction (Figs. 2a,b,c) - components NGC 5846A, NGC 5846B, NGC 5846C and NGC 5846D. Apparently we have resolved the single core observed with the VLA A-array (0 $.\!\!^{\prime\prime}$25 resolution) and VLA C-array (2 $.\!\!^{\prime\prime}$5 resolution; Filho et al. 2002). At 5 GHz, 11 mas separate the Southernmost (component A) from the Northern ones (components B and C). Component A appears in both the 2.3 GHz and 5 GHz maps (Figs. 2a,b), remaining unresolved at the latter frequency. However, it is not detected at 15 GHz (Fig. 2c), showing the source to be steep spectrum and/or variable. The steepness of the spectrum and high brightness temperature (Table 2) make it a candidate RSN (e.g. Tarchi et al. 2000). However, the compactness of the source, associated with the high brightness temperature and given its linear radio structure, suggests that component A is the variable radio core of NGC 5846.

3.2.2 Extended radio nuclear components

There are three (NGC 5354, NGC 5846 and also NGC 4552; Nagar et al. 2002), perhaps four (see discussion on NGC 5866 in Falcke et al. 2000) sources in the sample which show extended structure on mas-scales.

Although unresolved on arcsec-scales (Filho et al. 2000; Nagar et al. 2000), and in our VLBA maps, a recent high signal-to-noise VLBA image of NGC 4552 (Nagar et al. 2002) revealed a core plus two-sided emission in the East and West direction.

Our 0 $.\!\!^{\prime\prime}$2 resolution VLA map of NGC 5354 shows an unresolved core, whereas our longer integration VLBA map of this source shows about 4 mJy in symmetric double-sided emission to the East and West (Fig. 1c), suggesting the presence of jets cradling an AGN core. This extended emission is not seen, however, on our other VLBA shorter exposure maps, most probably due to low signal-to-noise.

Our first epoch 5 GHz VLBA map shows a triple source which we have designated as components NGC 5846A, B, and C. Furthermore, at 15 GHz we detect another component (NGC 5846D; Fig. 2c), about $\la$6 mas from the 5 GHz radio position of NGC 5846A and which cannot be identified with any of the components in the 5 GHz map. NGC 5846A has been discussed in the previous subsection, but the nature of the latter three components is more difficult to constrain. Components B, C and D appear only at one frequency each, and most probably have an inverted (NGC 5846D) or steep spectrum (NGC 5846B and NGC 5846C), suggestive of optically thin synchrotron emission. The apparent inverted nature of the spectrum of NGC 5846D could also be due to radio variability. All three components also appear to have a somewhat amorphous, resolved morphology and high brightness temperatures typical for non-thermal sources (Table 2). Furthermore, they are more or less aligned with component A in the North-South direction. One scenario is that these components are compact SNRs near the center of the galaxy (e.g. Tarchi et al. 2000) or they may be emission associated with a radio jet.

4 Archival and published X-ray observations

For most of our sample sources, published X-ray data are available - they mainly refer to low resolution observations not suited to our purposes. Therefore, in this section we present the results of the analysis of archival and/or published high resolution X-ray Chandra data for the sample sources.

4.1 Observations and data reduction

Public Chandra ACIS-S archival data are available for seven of the sample sources, three of which have already been published (NGC 1055 and NGC 3627; Ho et al. 2001 and NGC 5846; Trinchieri & Goudfrooij 2002). NGC 5866 Chandra data have been published (Terashima & Wilson 2003) but are not as yet publicly available.

For reasons of self-consistency and homogeneity, all publicly available Chandra data were reduced and analyzed using the CIAO (version 2.2) and XSPEC (version 11.2) packages. We removed high background time by excluding events exceeding $\pm$3$\sigma$ of the mean image count rate. The energy range was restricted to the 0.4-10 keV band. We searched for X-ray nuclei coincident with the radio core positions (Filho et al. 2000, 2002; Nagar et al. 2000, 2002; Sect. 3) to within the positional accuracy of Chandra ($\sim$1''). When found, spectra were extracted from a circular region with a diameter of 2.5'' ($\sim$5 pixels on the ACIS image) in order to spatially match both the VLA C-array radio maps (Filho et al. 2000, 2002) and the regions from where the H$\alpha $ luminosities were obtained (Ho et al. 1997). The background was extracted from an annulus around the source with an inner diameter of 5''.

Possible difficulties in pinpointing a hard X-ray nucleus can arise from the presence of ultraluminous X-ray sources (ULXs). Such confusion is minimized by the high resolution of Chandra. However, this possibility can not be totally excluded for galaxies like NGC 5846 (see Sect. 4.2), where several X-ray sources are detected in the nuclear region. With the possible exception of NGC 4552, pile-up effects should not be a problem, given the relatively low count rates. Even in the case of NGC 4552, a pile-up correction with a subsequent factor 2 difference in flux would not alter the following results.

For three of the X-ray sources (NGC 4552, NGC 5846, NGC 7331), the longer exposure times ( $\tau\simeq30$ ks), allowed a spectral analysis. We considered a power law model plus, for NGC 4552, a thermal plasma component (Raymond & Smith 1977). The column densities have been fixed to the Galactic values (Dickey & Lockman 1990). The fitted photon indices (1.5-2.3) are consistent with values obtained for LLAGN (Terashima & Wilson 2003). For sources observed with shorter exposure times, the count rates were converted to X-ray fluxes assuming a simple power-law spectrum with a photon index of 1.7 and Galactic absorption. Upper limits for non-detections were calculated after running a source detection procedure (wavdetect) on the full resolution Chandra image, with a threshold of 10-6 (4.7$\sigma$ level). Assuming a power-law ( $\Gamma=1.7$) and Galactic absorption, the upper limits were calculated to be equal to the flux of the weakest detected source on the image. Our results are consistent with those obtained by Ho et al. (2001) for NGC 1055 and NGC 3627 and Trinchieri & Goudfrooij (2002) for NGC 5846 (see Sect. 4.2).

4.2 Results of the X-ray observations

The X-ray parameters as given by Chandra archival or published data are specified in Table 3.

Five of eight (63%) sample sources with available Chandra data were detected as having hard X-ray nuclear cores. Given the hypothetical weakness of the hard X-ray core, the three non-detections are most likely related to the low integration time compared to the observations of the other sources; these have $\tau$ = 1-2 ks (Table 3).

Below we give a brief description of the X-ray properties of the AGN candidates.

Table 3: X-ray properties (published or archival). Column 1: source name. Column 2: exposure time. Column 3: photon index (and Raymond-Smith plasma temperature) of the power law fit and assumed. Column 4: unabsorbed X-ray flux in the 0.5-2 keV band. Column 5: unabsorbed X-ray flux in the 2-10 keV band. Column 6: logarithm of the soft X-ray luminosity (0.5-2 keV). Column 7: logarithm of the hard X-ray luminosity (2-10 keV). Column 8: note.

NGC 660 The soft X-ray region is about 0 $.\mkern-4mu^\prime$5 in extent in the Northeast-Southwest direction, about the size and orientation of the radio emitting region (see Filho et al. 2002). The hard X-ray nuclear region, spatially coincident with the radio core (Filho et al. 2002), comprises only about 10 counts, such that no spectral fit was possible.

NGC 1055 No obvious X-ray source - hard or soft - was found on the Chandra image (see also Ho et al. 2001). A source detection procedure (wavdetect) was run in order to obtain an upper limit to the X-ray emission.

NGC 3245 There is a hard nuclear X-ray source coincident with the optical nucleus. The soft X-ray emission is resolved and extended ($\sim$6'') in the North-South direction.

NGC 3627 The soft X-ray emission extends over about 2' in a Northwest-Southeast direction, similar to the extent and orientation of the triple-radio source (Filho et al. 2000). No obvious nuclear X-ray source - hard or soft - was found on the Chandra image (see also Ho et al. 2001). Although the source detection procedure (wavdetect) found a source positionally coincident with the VLA core (Filho et al. 2000; Nagar et al. 2000, 2002), it is the weakest source on the image and we prefer to consider it an upper limit (see also Ho et al. 2001).

NGC 4552 The soft X-ray emission extends almost to the limits of the optical galaxy ($\sim$2'), exhibiting a somewhat hourglass structure. We find a strong hard X-ray nuclear source, coincident with the radio core position (Nagar et al. 2002; Filho et al. 2000). The spectrum of the nuclear source (1200 net counts) yields a satisfactory fit with a two component model - power law ( $\Gamma=1.51$) plus Raymond-Smith thermal plasma (kT=0.95 keV).

NGC 5846 The X-ray morphology of this source is complex (see also Trinchieri & Goudfrooij 2002), extending over more than 3', almost to the limits of the optical galaxy. It exhibits two "spiral arms'', directed to the Sortheast and Southwest. There is a weak hard X-ray nucleus, coincident with the VLA radio core in NGC 5846, whose spectrum was fit with a power law ( $\Gamma=2.29$). The source we have identified with the nucleus coincides with the third most luminous X-ray source in NGC 5846, listed as source number 12 in Trinchieri & Goudfrooij (2002). Their quoted flux for this source is six times higher then our estimation because they convert counts to fluxes assuming a constant conversion factor. If we apply the same formulation, our results are consistent with theirs by a factor of 1.4.

NGC 5866 Recent published Chandra observations of this source failed to detect any hard X-ray emission from the nucleus (Terashima & Wilson 2003).

NGC 7331 The Chandra soft X-ray image of this galaxy shows several blotchy regions of emission, along about 3', in a North-South direction, similar to the morphology of the optical galaxy and also the radio (Cowan et al. 1994). Previous observations with ROSAT (Stockdale et al. 1998) have detected a luminous X-ray nucleus (few 1040 erg s-1). The Chandra observation has resolved the ROSAT emission into several X-ray sources and revealed a hard X-ray nucleus positionally coincident with the radio core detected in Cowan et al. (1994). The spectrum of the nuclear source is fit with a one component power law ( $\Gamma=1.46$) model and the luminosity is two magnitudes lower ($\sim$1038 erg s-1) than the ROSAT value.

5 Radio, optical and X-ray relations

In Table 4 we summarize the relevant properties for the sample sources from published, archival or new observations.

Table 4: Summary of the AGN candidate properties (published, archival or new). Column 1: source name. Column 2: adopted distance for $H_{\rm0}=75$ km s-1Mpc-1 (Tully 1988 or Ho et al. 1997a). Columns 3-5: extinction-corrected [OI]$\lambda $6300, H$\alpha $, and H$\beta $ line luminosity (Ho et al. 1997a, 2003). Column 6: published or new high resolution VLBA 5 GHz (or VLA; see footnote) radio core luminosity. When multiple 5 GHz, VLBA observations are available, the longest integration time observations is quoted. Column 7: archival and/or published Chandra hard X-ray (2-10 keV) nuclear luminosity. Column 8: nuclear B-band magnitude. Column 9: logarithm of the radio to nuclear B band ratio. Column 10: logarithm of the radio to hard X-ray ratio. Column 11: logarithm of the X-ray to H$\alpha $ ratio. Column 12: central velocity dispersion (Lyon/Meudon Extragalactic Database or HYPERCAT - NGC 660, 3627, 4419, 4527, and 5354). Column 13: estimated black hole mass obtained from the velocity dispersion. Column 14: logarithm of the Eddington luminosity. Column 15: logarithm of the X-ray to Eddingtion luminosity ratio.

With the exception of NGC 3245 (Barth et al. 2001), there are no reliable dynamical black hole mass estimates for the sample sources. The next best method for determining black hole mass utilizes the central velocity dispersion. The correlation between black hole mass and central velocity dispersion is tight (Ferrarese & Merritt 2000; Gebhardt et al. 2000; Tremaine et al. 2002):

\begin{eqnarray*}\displaystyle \log \left(\frac{M_{\rm BH}}{M_\odot}\right) = 8....
...; \log \left[ \frac{\sigma_{\rm e}}{200~{\rm km s^{-1}}} \right]
\end{eqnarray*}


where $\sigma _{\rm e}$ is the luminosity-weighted stellar velocity dispersion measured at the effective radius. However, in this case we will simply use the value $\sigma _{\rm c}$, the central velocity dispersion, since Gebhardt et al. (2000) have shown that the difference between  $\sigma _{\rm e}$ and  $\sigma _{\rm c}$ is only $\sim$10%.

As can be seen in Table 4, the dynamical black hole mass estimate for NGC 3245 ( $2\times10^8~M_{\odot}$; Barth et al. 2001) is in excellent agreement with the velocity dispersion estimate.

Nuclear B-band magnitudes have been calculated using the correlation found for low luminosity Seyfert 1s between $M^{\rm nuc}_{B}$ and the H$\beta $ (narrow plus broad) line luminosity (Ho & Peng 2001):

\begin{eqnarray*}\log~ (L_{\rm H\beta}) = -0.34 ~ {M^{\rm nuc}_{B}} + 35.1.
\end{eqnarray*}


With the exception of the Type 1.9 (broad H$\alpha $emission present) source NGC 1161, all the sample sources are Type 2 (no broad-line emission). Therefore, we have used only the narrow H$\beta $ line luminosities taken from Ho et al. (1997a) to estimate $M^{\rm nuc}_{B}$. In the case of composite galaxies, the derived $M^{\rm nuc}_{B}$ should be taken as an upper limit, due to contamination of the nucleus by stellar light and obscuration of the H$\beta $ line.

The radio-loudness parameter $R_{\rm o}$, was calculated as the ratio of the 5 GHz, mas-scale radio (except for NGC 7331) to nuclear B-band luminosity.

In comparison, the radio-loudness parameter $R_{\rm x}$ was calculated as the ratio of the 5 GHz, mas-scale radio (except for NGC 7331) to hard X-ray luminosity, as formulated in Terashima & Wilson (2003).

The Eddington luminosity - the maximum luminosity output of an accreting object in an isotropic homogeneous system - was calculated as (in units of erg s-1):

\begin{eqnarray*}L_{\rm Edd}=1.3 \times 10^{38}\left( \frac{M_{\rm BH}}{M_{\odot}} \right)\cdot
\end{eqnarray*}


5.1 Radio-loudness

Objects with log $R_{\rm o}>1$ are commonly said to be radio-loud (e.g. Kellermann et al.1994). Ho & Peng (2001) have recently challenged the long-thought notion that Seyfert galaxies are mainly radio-quiet galaxies. By using nuclear B-band magnitudes and radio powers, they found that as many as 60% of Seyferts are radio-loud.

An alternative method for parameterizing radio-loudness involves using the radio to hard X-ray ratio (Terashima & Wilson 2003). This method has the advantage of suffering less from extinction than the optical and allows the identification of the nucleus by means of high spatial resolution observations (e.g. Chandra). The equivalent radio-loud/radio-quiet boundary is established at log $R_{\rm x}=-4.5$(Terashima & Wilson 2003).

In Fig. 3 we have plotted the logarithm of the $R_{\rm x}$ versus $R_{\rm o}$ parameters, where the solid lines denote the radio-loud/radio-quiet boundaries.


  \begin{figure}
\par\includegraphics[width=8.6cm,clip]{0486fig3.eps}\end{figure} Figure 3: The radio to X-ray ratio ($R_{\rm x}$) versus radio to optical ($R_{\rm o}$) relation in logarithmic units. Upper limits in X-ray and radio luminosity are illustrated by arrows. The solid lines define the radio-loud/radio-quiet boundary for the optical and X-ray diagnostic.
Open with DEXTER

NGC 660, NGC 3245, NGC 3627 and NGC 7331 are found to be optically radio-quiet. The remaining sources are radio-loud, according to the log $R_{\rm o}>1$ criterion (e.g. Kellermann et al.1994). In particular, all of the mas-scale radio cores, as shown by the 5 GHz VLBA observations, are radio-loud. Our results support those of Ho & Peng (2001), whereby radio-loud objects are also found in disk-dominated hosts (see Filho et al. 2000, 2002).

As can be seen from the plot as well as Table 4, with the possible exception of NGC 660, NGC 3245 and NGC 7331, all sources are radio-loud according to the log $R_{\rm x}>-4.5$criterion. In particular, all the detected mas-scale cores are radio-loud.

The $R_{\rm o}$ and $R_{\rm x}$ results are consistent, in that all mas-scale radio cores are radio-loud. However, differences in the remaining sources most likely arise from stellar contribution. It is to be noted, for example, the high H$\alpha $ luminosities associated with NGC 660, NGC 3245, and NGC 3627.

5.2 Radio vs. H$\alpha $ and [O I] luminosity

It is well known that classical active galaxies like Seyferts, radio galaxies and LINERs display correlations between their radio and emission-line properties (Ho & Ulvestad 2001; Ulvestad & Ho 2001a; Zirbel & Baum 1995; Nagar et al. 2000, 2002). In particular, Nagar et al. (2000, 2002) also show that there is a correlation between the radio core luminosity and [O  I]$\lambda $6300 line luminosity in LINERs and low luminosity Seyferts.

The radio luminosity of the present sample does not seem to correlate with H$\alpha $ or with [O  I]$\lambda $6300 luminosity, but because of the small sample and large intrinsic scatter in the correlation we are not able to draw any conclusions. The absence of a radio and emission-line correlation for composite sources had already been noted in Filho et al. (2002) using lower resolution, VLA radio observations.

5.3 X-Ray vs. H$\alpha $ luminosity

By examining the hard X-ray to H$\alpha $ luminosity ratio, we wish to investigate whether the observed X-ray emission is related to the continuum source which powers the emission lines.

Such a correlation between the X-ray and H$\alpha $ luminosities is observed in classical AGN such as Seyferts and quasars (Ward et al. 1988) and more recently in LINER 1s (Terashima et al. 2000; Ho et al. 2001), indicating that the AGN is capable of powering the optical emission lines. Typical ratios for LLAGN are in the range:

\begin{eqnarray*}\displaystyle \log \left(\frac{L_{\rm X}}{L_{\rm H\alpha}}\right) \sim 1{ -} 2.
\end{eqnarray*}


In Fig. 4 we plot the relation for our sample sources, where the X-ray luminosity is taken from the hard X-ray band. The shaded area defines ratios that are typical for AGN-powered sources.


  \begin{figure}
\par\includegraphics[width=8.6cm,clip]{0486fig4.eps}\end{figure} Figure 4: The relation between hard X-ray core and H$\alpha $ emission line luminosity. Upper limits in X-ray luminosity are illustrated by arrows. The shaded area corresponds to typical values for LLAGN - ${\log \left(\frac{L_{\rm X}}{L_{\rm H\alpha}}\right) \sim 1 {-} 2}$.
Open with DEXTER

All of the sources fall short of the empirical X-ray to H$\alpha $ ratio observed in classical AGN, indicating that the observed X-ray luminosities are too low to generate the H$\alpha $ lines. There are three possible scenarios - either the emission lines are powered by a source other than the AGN, the ionizing spectrum is different than that inferred, or the X-ray emission is heavily obscured at energies above 2 keV. The first possibility has already been investigated for example in the work of Maoz et al. (1998), whereby circumnuclear stellar populations may provide enough ionizing photons to explain the observed optical emission-line flux. However, as investigated by Ho et al. (2003; see also Ho et al. 2001), there is no evidence for the presence of young stellar clusters. Alternatively, it may be that the underlying ionizing spectrum is more complicated then assumed - most of the ionizing photons may occur in the extreme ultraviolet (EUV) and soft X-ray region, not in the hard X-rays. Powering of the emission lines from a central AGN could still occur if the hard X-ray flux were being heavily absorbed (Terashima et al. 2000; Terashima & Wilson 2003). If intrinsic X-ray luminosities are 1-2 orders of magnitude higher, then the X-ray to H$\alpha $ ratio would be typical for LLAGN (Terashima & Wilson 2003). Although scattering is still a viable explanation, observations at energies above 10 keV are needed in order to observe any Compton-thick component.

5.4 Radio vs. X-ray luminosity

In Fig. 5 we plot the mas-scale radio core luminosity versus the hard X-ray core luminosity for the sample sources.


  \begin{figure}
\par\includegraphics[width=8.6cm,clip]{0486fig5.eps}\end{figure} Figure 5: The relation between the mas-scale radio core and hard X-ray core luminosity. Upper limits in radio and X-ray luminosity are illustrated by arrows. The solid line defines the radio-loud/radio-quiet boundary for the X-ray diagnostic.
Open with DEXTER

It can be seen from the plot that the detected mas-scale radio cores have $L_{\rm R}>10^{36.5}$ erg s-1, with two sources (NGC 4552 and NGC 5846) being firm Chandra detections (Table 3). Three of the sources with no mas-scale radio detection, however, show hard X-ray Chandra detections - NGC 660, NGC 3245 and NGC 7331. The Chandra observations of NGC 1055 and NGC 3627 (see also Ho et al. 2001) failed to detect any hard X-ray nuclear source, most likely due to the low integration time for the X-ray observations ($\tau<2$ ks; Table 3). Combining both the radio and X-ray data yields a $\sim$50% LLAGN detection rate among the sample sources or equivalently a 12% detection rate for the entire composite LINER/H  II galaxy sample.

NGC 7331 has a very low luminosity VLA radio core (Cowan et al. 1994), 10 times that of Sgr A*, and also the lowest hard X-ray core luminosity of all the sample sources, 104 times that of Sgr A*. Similarities in the arcsec-scale radio properties - radio morphology and low radio flux density - of NGC 7331 and the cores that were not detected on mas-scales, indicate that NGC 7331 may actually be a prototype. These sources may all harbour submillijansky radio and low luminosity hard X-ray AGN cores, difficult to detect unless deep observations are obtained. Nevertheless, the results demonstrate the combined power of high resolution X-ray and radio imaging in LLAGN detection. We will return to this issue in the next section.


  \begin{figure}
\par\includegraphics[width=8.6cm,clip]{0486fig6.eps}\end{figure} Figure 6: The relation between the mas-scale radio core/hard X-ray core luminosity ratio ( $L_{\rm 5~GHz}$/ $L_{\rm
2{-}10~keV}$) and the hard X-ray core luminosity. Upper limits in radio and X-ray luminosity are illustrated by arrows. Numbers in parentheses refer to the logarithm of the black hole masses and the respective error, considering an error in log $\sigma _{\rm c}=0.06$, 0.02 and 0.01 plus a 10% error in the conversion of $\sigma _{\rm e}$ to $\sigma _{\rm c}$. The solid lines refer to the ADAF model prediction - ${L_{\rm R} \propto M_{\rm BH}~ L_{\rm X}^{0.1}}$- for several black hole masses - ${\log ~ M_{\rm BH}=7.0, 8.0,
9.0}$ - in logarithmic units of solar mass.
Open with DEXTER

6 ADAFs or jets?

We wish to use our observations to test models of radio and X-ray emission in low luminosity sources, namely advection-dominated accretion flows (ADAFs) and jet models. These models have been extensively applied to several other low luminosity sources, of which Sgr A* (Yuan et al. 2002a, 2003) and NGC 4258 (Yuan et al. 2002b) are just two examples.

ADAFs (see reviews by Narayan et al. 1998 and Quataert 2001) occur in the low accretion rate regime ( $m_{\rm
acc}\rm (crit)<10^{-1.6}$, where $ m_{\rm acc}$ is given in dimensionless units of Eddington accretion rate). They are characterized by a low radiative efficiency, due to the advection of a significant amount of the energy into the black hole. ADAF theory has developed quickly in recent years with the inclusion of outflows and the possible existence of nonthermal electrons (see Yuan et al. 2003 for a review). However, the analysis and therefore conclusions of this paper are based on the canonical ADAF model.

ADAFs may be able to account for correlations observed between the radio and X-ray emission. Particularly, this type of low radiative efficiency flow is supported by the highly sub-Eddington luminosities observed in our sources ( $\frac{ L_{\rm X}}{L_{\rm Edd}}<10^{-5}$; Table 4). Furthermore, ADAFs are naturally bright in radio due to radio synchrotron cooling and dim in the UV/optical regime due to the absence of a "blue bump'', the signature of a standard, optically thick but geometrically thin accretion disk (Shakura & Sunyaev 1973). Both of these ADAF characteristics conspire to produce "radio-loud'' cores.

In ADAFs the hard X-ray emission arises from bremsstrahlung and Compton scattering of synchrotron photons subject to self-absorption. As a result, the ADAF model predicts the correlation (Yi & Boughn 1998, 1999):

\begin{eqnarray*}L_{\rm R} \propto {M_{\rm BH}}~{L_{\rm X}^{0.1}}.
\end{eqnarray*}


In Fig. 6 we plot the mas-scale radio core to hard X-ray ratio versus the hard X-ray core luminosity. In parentheses are the logarithms of the black hole mass plus error estimations, which include the error in velocity dispersion plus a 10% error associated with the conversion from $\sigma _{\rm e}$ to $\sigma _{\rm c}$. The lines refer to the ADAF model prediction for different black hole masses, appropriately scaled to 5 GHz.

If the sources with upper limits to their mas-scale radio emission contain accreting black holes as is the case for NGC 7331 (due to a hard X-ray detection), then it is possible to explain their radio emission using the ADAF model, as originating in underfed black holes ( $m_{\rm acc}<10^{-3}$). However, for the mas-scale detected radio cores, the ADAF model overpredicts the black hole masses by a factor 10 or so (Table 4) even when the errors associated with the black hole mass determination are considered (Fig. 6). This suggests that the ADAF model is not entirely consistent with the observed radio to X-ray ratio unless the hard X-rays are heavily obscured (see Sect. 5.3).

In Fig. 7 we plot the mas-scale radio core luminosity versus the black hole mass. The solid lines refer to the ADAF prediction of:

\begin{eqnarray*}L_{\rm R} \propto {m_{\rm acc}}^{1.2}~ {M_{\rm BH}^{1.6}}
\end{eqnarray*}


where the relation has been appropriately scaled to 5 GHz. Typical errors in black hole mass are given.


  \begin{figure}
\par\includegraphics[width=8.4cm,clip]{0486fig7.eps}\end{figure} Figure 7: The relation between the mas-scale radio core luminosity and the black hole mass. Upper limits in radio luminosity are illustrated by arrows. We give typical error bars for the black hole mass, considering an error in log $\sigma _{\rm c}=0.06$, 0.02 and 0.01 plus a 10% error in the conversion of $\sigma _{\rm e}$ to $\sigma _{\rm c}$. The solid lines refer to the ADAF model prediction - ${L_{\rm R} \propto
m_{\rm acc}^{1.2}~M_{\rm BH}^{1.6}}$ - for different mass accretion rates - ${m_{\rm acc}=10^{-2}, 10^{-3}, 10^{-4}}$ - given in dimensionless units of Eddington accretion rate.
Open with DEXTER

The results show that composite sources with mas-scale radio cores have a radio luminosity to black hole mass ratio similar to low luminosity Seyferts and LINERs (log ($L_{\rm R}$/ $M_{\rm BH}$$\sim$ 1.14; Nagar et al. 2002). On the other hand, if the sources with no mas-scale radio detection harbour an AGN, then on average they should have lower radio to black hole ratios than the LLAGN in the Nagar et al. (2002) sample.

Figure 7 also shows the ADAF model prediction for several accretion rates. Sources with upper limits to their mas-scale radio emission are consistent with low mass accretion rates ( $m_{\rm acc}<10^{-2}$). However, the detected mas-scale cores fall very close to the ADAF permitted accretion limit of $m_{\rm acc}({\rm crit})=10^{-1.6}$. These results suggest that the mas-scale radio cores show more radio emission than can perhaps be accounted for from an ADAF.

Plotted in Fig. 8 is the relation between the hard X-ray luminosity and the black hole mass. The solid lines refer to the predicted correlation in an ADAF model for various mass accretion rates:

\begin{eqnarray*}L_{\rm X}\propto M_{\rm BH} m_{\rm acc}^{\rm x}
\end{eqnarray*}


where x=2 for thermal emission. Again, typical errors in black hole mass determination are given in the plot.


  \begin{figure}
\par\includegraphics[width=8.4cm,clip]{0486fig8.eps}\end{figure} Figure 8: The relation between the hard X-ray core luminosity and black hole mass. Upper limits in X-ray luminosity are illustrated by arrows. We give typical error bars for the black hole mass, considering an error in log $\sigma _{\rm c}=0.06$, 0.02 and 0.01 plus a 10% error in the conversion of $\sigma _{\rm e}$ to $\sigma _{\rm c}$. The solid lines refer to the ADAF model prediction - ${L_{\rm X} \propto M_{\rm BH} m_{\rm acc}^{2}}$ - for various accretion rates - ${m_{\rm acc}=10^{-2}, 10^{-3}, 10^{-4}}$ - given in dimensionless units of Eddington accretion rate.
Open with DEXTER

The X-ray data are consistent with the ADAF model for low accretion rates ( $m_{\rm acc}<10^{-3}$), both for the sources with mas-scale upper limits and for those with mas-scale cores. This could suggest the same underlying mechanism for X-ray production in all sample sources - an ADAF. In fact, in two of the mas-scale radio cores where an X-ray spectral fit was possible (NGC 4552 and NGC 7331; see Table 3), the X-ray photon index is consistent with the predicted ADAF value of $\sim$1.4.

In summary, if we take NGC 7331 as a prototypical submillijansky AGN core, then scaling arguments suggest that an ADAF could simultaneously explain both the radio and X-ray emission in the sources with no mas-scale radio detection. In fact, these sources may in fact be "bare'' ADAFs, where the radio and X-ray contribution is coming mainly from the accretion flow. The absence of a jet or of an energetically significant contribution of a jet to the radio emission is not, however, clear. It may be that, like in M 81, the jet is present but very compact and difficult to detect (Bietenholz et al. 2000).

On the other hand, thermal ADAF models predict inverted radio spectra ( $\alpha \sim {-} 0.4$; Narayan et al. 1998; see also Yuan et al. 2003) while the mas-scale sources possess flat spectrum radio cores (Table 2). Modified ADAF models that include convection (Ball et al. 2001) or outflows (Quataert & Narayan 1999) do not provide satisfactory solutions, since the predicted spectral slopes are roughly the same as in the canonical ADAF models (Quataert & Narayan 1999)[*]. Therefore, while the X-ray emission in mas-scale radio cores may be ADAF-powered, the radio emission may require another emission mechanism.

The most likely scenario is that the radio emission in the mas-scale detected radio cores is a result of an ADAF-type accretion (or variants thereof) plus a compact jet or outflow (Quataert et al. 1999; Yuan 2000; Falcke et al. 2001; Nagar et al. 2001; Ulvestad & Ho 2001b; Nagar et al. 2002). The consideration of a jet is not ad hoc since we detect extended emission in 4 out of the 5 mas-scale radio cores. Furthermore, the addition of a self-absorbed jet component to an ADAF core may "flatten'' out the spectral index to the observed values. In fact, flat spectrum radio emission in "classical'' AGN are commonly attributed to the presence of jets.

In a recent paper, Falcke et al. (2004; see also Merloni et al.2003) investigated an accretion rate-dependent scheme for black hole-powered sources, whereby sub-Eddington systems are predicted to be dominated by emission from a relativistic jet. Such an idea was suggested already in Falcke & Biermann (1995), and a jet model was shown to analytically predict the radio/X-ray correlation seen in the low/hard state of the black hole candidate X-ray binary (XRB) GX 339-4 (Markoff et al. 2003). This correlation has now been seen in many other low/hard state XRBs, and is likely universal (Gallo et al. 2003). Falcke et al. (2004; see also Merloni et al.2003) show that this correlation can be scaled with black hole mass to extend to certain classes of AGN as well. This leads to a radio/X-ray/black hole mass fundamental plane which unifies diverse sources such as FR I radio galaxies, LINERs and XRBs as jet-dominated (see Fig. 1 in Falcke et al. 2004).

This scaling relation is given by (Eq. (7); Falcke et al. 2004):

\begin{eqnarray*}\displaystyle L^{\rm cor}_{\rm X}\propto L_{\rm X'} \left(\frac...
...)^{-\alpha_{\rm X}} M_{\rm BH}^{-m\alpha_{\rm R}+\alpha_{\rm X}}
\end{eqnarray*}


where $\alpha_{\rm R}$ and $\alpha_{\rm X}$ are the radio and X-ray spectral indices and m is a function of both $\alpha_{\rm R}$ and $\alpha_{\rm X}$. We here apply the black hole mass scaling relation in order to investigate the presence of an energetically dominant jet in our sample sources. We have extrapolated our 2-10 keV luminosities (Table 3) to 3-9 keV assuming an average $\alpha_{\rm X}=0.6$ ( ${\equiv} \Gamma-$1) and $\alpha_{\rm R}=-0.15$ (Falcke et al. 2004). The assumption of a spectral index of $\Gamma=1.6$ ( $\alpha_{\rm X}=0.6$) is consistent within 10% of values obtained using our actual tabulated spectral indices (Table 3) and therefore does not alter the general result.

In Fig. 9 we have plotted the scaled correlation normalized to the $6~M_{\odot}$ black hole mass of GX 339-4. Typical errors in the "black hole mass-corrected'' X-ray luminosity are conservatively estimated to be of the order of 0.08 in logarithmic units of luminosity.


  \begin{figure}
\par\includegraphics[width=8.6cm,clip]{0486fig9.eps}\end{figure} Figure 9: The relation between radio and "black hole mass-corrected'' X-ray luminosity, scaled to the $6~M_{\odot}$ black hole in XRB GX 339-4. The solid line is the predicted radio/X-ray correlation from the jet model and the ellipse details the broadband observations of GX 339-4 (Corbel et al. 2003; Markoff et al. 2003).
Open with DEXTER

In addition to the scenario discussed above, a jet model alone seems to interpret the present data as well. It is noteworthy that the mas-scale detected sources fall very close to the jet model predicted relation. The deviation of the sources NGC 4552, NGC 5846 and NGC 5866 is in the direction of excess radio emission, which is entirely consistent with the jet model. The scaling relation is valid for the compact inner jet core while these sources are known to possess jet emission on mas-scales; hence, the radio flux density likely includes some extended jet emission.

7 Conclusions

We have embarked on a high resolution X-ray and radio imaging study of a sample of composite LINER/H  II galaxies known to exhibit AGN-like properties. Five of the 16 sample sources possess milliarcsecond radio cores, four of which also show extended radio emission. Two of the sources with no milliarcsecond radio detection were found to possess hard X-ray nuclear cores. The 50% detection rate of radio and/or hard X-ray nuclei support the view that these AGN candidate composite H  II galaxies host weak AGN central engines. When considering the entire composite LINER/H  II sample, these results yield an overall 12% detection rate. Modeling of the radio and X-ray emission suggests that the mas-scale radio detected sources exhibit an energetically important contribution from a jet, associated with a sub-Eddingtion accreting black hole. The sources that were not detected on mas-scales, however, appear consistent with the radio and X-ray emission coming mainly from the low radiative efficiency accretion flow (ADAF). Nonetheless, even in the cases where the presence of an AGN is irrefutable, this AGN cannot be responsible solely for powering the observed emission lines unless the X-ray flux is heavily obscured.

Acknowledgements

M.E.F. acknowledges support from the Fundação para a Ciência e Tecnologia, Ministério da Ciência e Ensino Superior, Portugal through the grant PRAXIS XXI/BD/15830/98 and SFRH/BDP/11627/2002. We would like to thank Heino Falcke for the extensive discussions on ADAF and jet models and Jim Ulvestad for providing us with calibrated VLA data for NGC 660. We wish also to thank the anonymous referee for his helpful comments.

This research has made extensive use of NED (NASA/IPAC Extragalactic Database), which is operated by the Jet Propulsion Laboratory, California Institute of Tecnology, under contract with the National Aeronautics and Space Administration, LEDA (Lyon/Meudon Extragalactic Database), HYPERCAT and the Chandra data archive. The VLA and the VLBA are facilities of the National Radio Astronomy Observatory (NRAO) which is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.

References



Copyright ESO 2004