A&A 418, 429-443 (2004)
DOI: 10.1051/0004-6361:20034486
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
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 (
;
)
indicates synchrotron self-absorption in
components whose brightness temperatures are comparable to the kinetic
temperatures (mc2/k;
K) of the relativistic
electrons producing the synchrotron radiation. Steep spectrum
emission in star-forming galaxies (
)
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 (
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 (
K; Condon 1992) or non-thermal
(
TB>105 K) and whether optically thin thermal
(
), self-absorbed synchrotron (
), or
optically thin synchrotron (
-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.
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
(
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.
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
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.
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
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
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
0.07 mJy/beam (uniform weighting; Perley et al. 1989) to within a factor of two for the VLA data and
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.
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.
![]() |
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
(
|
| Open with DEXTER | |
![]() |
Figure 2:
Same as Fig. 1. a) NGC 5846
(9.2 |
| Open with DEXTER | |
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.
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
limits to any 5 GHz emission in order to obtain upper
limits for the spectral indices. We calculate the brightness
temperature as:
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
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
25 resolution) and VLA C-array
(2
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.
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
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
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.
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.
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
3
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 (
1''). When
found, spectra were extracted from a circular region with a diameter
of 2.5'' (
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
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
(
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
level). Assuming a
power-law (
)
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).
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
= 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
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 (
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 (
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 (
)
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 (
).
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 (
)
model and the luminosity
is two magnitudes lower (
1038 erg s-1) than the ROSAT value.
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
km s-1Mpc-1 (Tully 1988 or Ho et al. 1997a).
Columns 3-5: extinction-corrected [OI]
6300, H
,
and
H
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
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):
where
is the luminosity-weighted stellar velocity
dispersion measured at the effective radius. However, in this case we
will simply use the value
,
the central velocity dispersion,
since Gebhardt et al. (2000) have shown that the difference between
and
is only
10%.
As can be seen in Table 4, the dynamical black hole mass estimate
for NGC 3245 (
;
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
and
the H
(narrow plus broad) line luminosity (Ho & Peng 2001):
With the exception of the Type 1.9 (broad H
emission present) source NGC 1161, all the sample sources are Type 2 (no
broad-line emission). Therefore, we have used only the narrow
H
line luminosities taken from Ho et al.
(1997a) to estimate
.
In the case of composite
galaxies, the derived
should be taken as an
upper limit, due to contamination of the nucleus by stellar light and
obscuration of the H
line.
The radio-loudness parameter
,
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
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):
Objects with log
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
(Terashima & Wilson 2003).
In Fig. 3 we have plotted the logarithm of the
versus
parameters, where the solid lines denote the
radio-loud/radio-quiet boundaries.
![]() |
Figure 3:
The radio to X-ray ratio ( |
| 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
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
criterion. In particular, all the detected mas-scale cores are
radio-loud.
The
and
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
luminosities associated
with NGC 660, NGC 3245, and NGC 3627.
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]
6300 line luminosity in
LINERs and low luminosity Seyferts.
The radio luminosity of the present sample does not seem to correlate
with H
or with [O I]
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.
By examining the hard X-ray to H
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
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:
![]() |
Figure 4:
The relation between hard X-ray core and H |
| Open with DEXTER | |
All of the sources fall short of the empirical X-ray to H
ratio observed in classical AGN, indicating that the observed
X-ray luminosities are too low to generate the H
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
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.
In Fig. 5 we plot the mas-scale radio core luminosity versus the hard X-ray core luminosity for the sample sources.
![]() |
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
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 (
ks; Table 3). Combining
both the radio and X-ray data yields a
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.
![]() |
Figure 6:
The relation between the mas-scale radio
core/hard X-ray core luminosity ratio (
|
| Open with DEXTER | |
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 (
,
where
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
(
;
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):
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
(
). 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:
where the relation has been appropriately scaled to 5 GHz. Typical errors in black hole mass are given.
![]() |
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
|
| 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 (
/
)
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 (
). However, the
detected mas-scale cores fall very close to the ADAF
permitted accretion limit of
.
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:
where x=2 for thermal emission. Again, typical errors in black hole mass determination are given in the plot.
![]() |
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
|
| Open with DEXTER | |
The X-ray data are consistent with the ADAF model for low accretion
rates (
), 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
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
(
;
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):
In Fig. 9 we have plotted the scaled correlation normalized to the
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.
![]() |
Figure 9:
The relation between radio and "black hole mass-corrected''
X-ray luminosity, scaled to the
|
| 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.
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.