A&A 374, 42-65 (2001)
DOI: 10.1051/0004-6361:20010639
M. Dahlem1 - J. S. Lazendic2,3 - R. F. Haynes3 - M. Ehle4,5 - U. Lisenfeld6
1 - Sterrewacht Leiden, Postbus 9513, 2300 RA Leiden,
The Netherlands
2 -
Astrophysics Department, School of Physics A28, University of
Sydney, NSW 2006, Australia
3 -
Australia Telescope National Facility, PO Box 76, Epping, NSW
2121, Australia
4 -
XMM-Newton Science Operations Centre, Apartado 50727, 28080
Madrid, Spain
5 -
Astrophysics Division, Space Science Department of ESA, ESTEC,
2200 AG Noordwijk, The Netherlands
6 -
IRAM, Avenida Divina Pastora 7, NC, 18012 Granada, Spain
Received 20 December 2000 / Accepted 1 May 2001
Abstract
We present radio continuum observations conducted with the VLA
and ATCA of a sample of 15 edge-on spiral galaxies. 11 of these
galaxies, with inclination angles of
and
neither active galactic nuclei nor nearby interaction partners,
are suitable for studies of halo properties in relation to the
level of star formation in their disks.
In 6 of these 11 galaxies radio halos were detected at the angular
resolution of the current data. In the remaining cases the presence
of halo emission could not be proven unambiguously, partly due to
relatively low angular resolution.
A clear trend was found that galaxies with radio halos are those
with the highest far-infrared 60
m to 100
m flux ratios.
This shows the suitability of high
f60/f100 ratios of
0.4 as a reliable tracer of galaxies with high star
formation rates and related disk-halo interactions, leading to
the presence of extraplanar emission, e.g. from cosmic ray
electrons. The measured exponential scale heights of those 6
radio halos that were clearly detected range from about 1.4 to
3.1 kpc.
All 4 physically small galaxies in our sample do show extraplanar
synchrotron radio emission, indicating that their more shallow
gravitational potential compared to normal-sized spirals might
facilitate the escape of cosmic-ray electrons from the sites of
star formation in their disks.
Although the galaxies with the highest energy input rates into
the ISM of their disks are those that have the most prominent
radio halos, there is no direct relation between the halo scale
heights and the energy input rates. Instead, the scale heights
of the radio halos are dominated by the energy losses of the
cosmic ray electrons on their way out of the galaxy disks.
Key words: ISM: general - galaxies: evolution - galaxies: halos - galaxies: starburst - radio continuum: galaxies
Despite accumulating evidence for the existence of gaseous halos around a number of nearby spiral galaxies (see, e.g., Dettmar 1992; Dahlem 1997 for overviews), it is currently a matter of controversy how many and exactly which galaxies have such halos. While the existence of gaseous halos (including radio halos, i.e. cosmic ray electrons) seems to be in general a rather rare phenomenon (Hummel et al. 1991a found that only about 5.5% of the 181 galaxies in their initial optically selected sample show signs of extraplanar radio emission), their existence appears to be quite common among nearby galaxies with high star formation rates (SFRs; e.g. Dahlem 1997). We find that this apparent discrepancy is caused by two factors.
First, different selection criteria were used to define the sample of candidate galaxies. Hummel et al.'s objects were selected by their optical (blue) magnitude. However, especially when studying edge-on systems (in searches for halos), the blue magnitude is not always a good indicator of the level of star formation (SF) because of dust extinction in the disks. The fraction of galaxies considered by Hummel et al. (1991a) as good candidates for follow-up searches of radio halos, based on extended emission found in radio synthesis images, is 10 out of 18, i.e. 56%.
In previous studies (Dahlem et al. 1995, 1998; hereafter DLG95
and D98, respectively) we have started using a far-infrared
criterion, originally defined by Heckman et al. (1990) to
characterise starburst galaxies, for selecting galaxies with
warm dust as targets for searches of gaseous halos: their 60
m to 100
m far-infrared (FIR) IRAS flux ratio.
Objects with
have warm dust, heated
by massive stars, with average dust temperatures,
K. It turns out that all edge-on galaxies
selected by us in this way (with 60
m far-infrared [FIR]
fluxes of
Jy) show clear signs of halo emission.
Emission is detected not only in the radio, but also in other
wavebands (D98). Galaxies with known luminous AGNs (that might
contribute significantly to the heating of the dust) were excluded
from our sample.
Second, contrary to searches of halo gas in very large samples, such as conducted by Hummel et al. (1991a), our present approach of carefully selecting target galaxies is combined with the employment of high-sensitivity observations to detect low surface brightness emission.
Irwin et al. (1999) conducted a search for radio halos in a
quite heterogeneous sample of target galaxies. Despite this
heterogeneity, these authors achieved a high detection rate of
halo emission of >
.
A discussion of the interpretation
of their data in the context of the new results presented here
will be provided in Sect. 4.4.
Here we present the results from our radio observations of a
number of FIR-warm edge-on galaxies with 3 Jy
< f60 < 30
Jy FIR flux densities conducted with the Australia Telescope
Compact Array (ATCA
) and the Very Large Array
(VLA
).
These galaxies are fainter than the galaxies studied by us
previously, because there are too few nearby systems (with
Jy) to derive general properties of their
halos. Most objects were selected to be FIR-warm, with
,
except NGC1055, NGC2820 and
NGC7090, which have slightly lower FIR flux ratios. The
primary goal is to investigate whether indeed all FIR-warm
galaxies (without AGNs) do have gaseous halos and an attempt
at detecting these at distances beyond that of the Virgo
cluster.
| Galaxy | Arraya | Date | Int. time |
| [h:mm] | |||
| NGC1511 | 750A | 1998, May 03 | 9:00 |
| 750B | 1997, Aug. 10 | 9:30 | |
| 750C | 1997, Oct. 20 | 9:30 | |
| 1.5D | 1998, Oct. 19 | 4:50 | |
| 1.5D | 1998, Oct. 20 | 6:00 | |
| 1.5C | 1999, Apr. 10 | 8:50 | |
| NGC7090 | 750A | 1998, May 02 | 9:30 |
| 750B | 1997, Aug. 06 | 9:30 | |
| 750C | 1997, Oct. 18 | 9:30 | |
| 1.5D | 1998, Oct. 17 | 5:20 | |
| 1.5D | 1998, Oct. 20 | 3:00 | |
| 1.5B | 1999, Apr. 02 | 8:40 | |
| NGC7462 | 750A | 1998, May 04 | 9:00 |
| 750B | 1997, Aug. 08 | 9:30 | |
| 750C | 1997, Oct. 19 | 9:30 | |
| 1.5D | 1998, Oct. 18 | 2:30 | |
| 1.5D | 1998, Oct. 20 | 1:00 |
| Note to Table 1: |
| a) Full aperture synthesis for each ATCA array is obtained by combining data from 4 configurations, named A, B, C and D. Three configurations already provide good uv coverage. |
The galaxies observed by us with the ATCA are NGC1511, NGC7090 and NGC7462. Each was observed with three different configurations of the 750 m array and two different configurations of the 1.5 km array, see Table 1. The individual observing runs were 13 hours long (including time for calibration), providing almost full 12 hour aperture syntheses in each configuration. The total on-source integration times range from 32 h (NGC7462) to 47 h (NGC1511). The ATCA is capable of observing at two frequencies simultaneously. We observed the radio continuum at 2.45 GHz and 1.43 GHz (13 cm and 21 cm wavelength, respectively). All data are polarisation calibrated. Due to technical problems part of the 1.5 km array observations of NGC7462 could not be used.
1934-638 was used as the primary flux calibrator, 0407-658 (for NGC1511) and 2106-413 (for both NGC7090 and NGC7462) as polarisation and phase calibrators. The adopted flux of 1934-638 is 11.14 (14.94) Jy at 2.45 (1.43) GHz. The data reduction was performed in a standard fashion, using the software package MIRIAD.
With angular extents of their radio emission of
4' (see
below), the observed galaxies are so small that they fit easily
into the primary beam of the ATCA's 22-m antennae (Full Width at
Half Maximum FWHM = 20'/34' at 2.45/1.43 GHz, respectively)
and no primary beam correction is necessary. With a shortest
spacing of 45.9 m and good uv coverage, flux losses due to missing
short spacings are negligible.
| Galaxy | Arraya | Date | Int. time |
| [h:mm] | |||
| NGC1055 | D | 1999, Mar. 25/29 | 3:30 |
| NGC1406 | DnC | 1999, Feb. 16 | 2:45 |
| NGC1421 | D | 1999, Apr. 02 | 2:45 |
| NGC2748 | D | 1999, Apr. 04 | 3:30 |
| NGC2820 | D | 1999, Apr. 04 | 3:05 |
| NGC3175 | DnC | 1999, Feb. 18 | 2:00 |
| NGC3437 | D | 1999, Mar. 29 | 3:00 |
| NGC3717 | DnC | 1999, Feb. 18 | 2:00 |
| NGC4527 | D | 1999, May 28 | 1:30 |
| NGC4700 | D | 1999, May 28 | 1:35 |
| NGC5073 | D | 1999, May 28 | 1:35 |
| NGC7541 | D | 1999, Apr. 02 | 2:05 |
| Note to Table 2: | |
| Open with DEXTER | |
| a) Southern galaxies (
i.e., with an extended northern arm of the array. | |
Those galaxies observed with the VLA are listed in Table 2. Two IFs of 50 MHz bandwith were used, centred at 1.465 GHz and 1.385 GHz, respectively. Primary flux calibrators are either 0137+331 (3C 48; 15.62 Jy at 1.465 GHz and 16.32 Jy at 1.385 GHz) or 1331+305 (3C 286; 14.55/14.94 Jy, respectively) or both, if available, adopting the flux scale by Baars et al. (1977). The combined final maps have a centre frequency of 1.425 GHz. The data reduction was performed in a standard way, utilizing the NRAO software package AIPS.
Due to the shortness of most of the observing runs, no polarisation calibration could be obtained (owing to insufficient parallactic angle coverage).
For galaxies with an extent of 5' or larger a primary beam correction was performed. This affects NGC2820 (and its partners), NGC3717 and NGC4527. For NGC3717, one of the most extended objects in our sample, a test was performed on the influence of the primary beam correction on the total flux measurement, which turned out to be negligible (0.1 mJy with respect to a total flux of 235 mJy).
For investigations of the z height of the observed radio emission it is important to have the highest possible resolution in the direction perpendicular to the disk plane of the sample galaxies, while normally interpretation of the results is easiest with a circular beam. In those cases where the clean beam was elongated roughly along the minor axis of the observed galaxy it was restored circular by slightly re-weighting the visibilities. In a few cases, where the beam was elongated in a direction close to the major axis of the target galaxy, it was left elliptical, but restored with its major axis oriented exactly along the galaxy disk. This leaves the beam's minor axis, and thus the highest resolution, perpendicular to the galaxy disk, while keeping the data analysis in the z direction straight-forward.
| Galaxy | D | i a |
|
Beamb | S1.43 | Pos. Angle | |||
| (2000) | (2000) | (Mpc) | ( |
FWHM ('') | (mJy beam-1) | (mJy) | PA ( |
||
| NGC1055 | 02 41 43.7 | +00 25 54 | 16.0 | 67 | 0.34 | 54 | 0.08 |
|
|
| NGC1406 | 03 39 22.5 | -31 19 19 | 14.9 | 90 | 0.44 | 43 | 0.08 |
|
|
| NGC1421 | 03 42 29.5 | -13 29 23 | 31.1 | 90 | 0.40 |
|
0.10 |
|
|
| NGC1511 | 03 59 35.7 | -67 38 07 | 17.5 | 72 | 0.62 | 17.7c/33.0d | 0.070c/0.05d |
|
|
| NGC2748 | 09 13 44.6 | +76 28 32 | 28.7 | 71 | 0.39 | 52 | 0.10 |
|
|
| NGC2820 | 09 21 44.1 | +64 15 26 | 29.9 | 90 | 0.33e |
|
0.06 |
|
|
| NGC3175 | 10 14 42.3 | -28 52 20 | 15.9 | 85 | 0.46 |
|
0.07 |
|
|
| NGC3437 | 10 52 34.9 | +22 56 04 | 25.5 | 75 | 0.57 |
|
0.10 |
|
|
| NGC3717 | 11 31 31.8 | -30 18 32 | 27.1 | 90 | 0.44 | 41 | 0.08 |
|
|
| NGC4527 | 12 34 08.5 | +02 39 11 | 31.9 | 68 | 0.47 |
|
0.15 |
|
|
| NGC4700 | 12 49 07.3 | -11 24 46 | 25.5 | 90 | 0.57 |
|
0.10 |
|
|
| NGC5073 | 13 19 20.8 | -14 50 35 | 43.9 | 87 | 0.69 |
|
0.15 |
|
|
| NGC7090 | 21 36 28.6 | -54 33 26 | 11.7 | 90 | 0.33 | 17.7c/34.5d | 0.075c/0.05d |
|
|
| NGC7462 | 23 02 46.5 | -40 50 07 | 15.1 | 90 | 0.49 | 23.5c/46.5d | 0.085c/0.05d |
|
|
| NGC7541 | 23 14 43.0 | +04 32 05 | 42.4 | 75 | 0.48 | 47 | 0.15 |
|
|
The 12 galaxies observed with the VLA are displayed in Figs. 1-3, 5-12 and 15. For these objects only 1.425 GHz data are currently available.
ATCA radio continuum images of the three galaxies NGC1511, NGC7090 and NGC7462 at 2.45 GHz and 1.43 GHz are presented in Figs. 4, 13 and 14. Each figure displays the 2.45 GHz image on the left side and the 1.43 GHz map on the right.
In the following we briefly describe the radio images of each galaxy
in our small sample of 15 systems. Note that the first positive
contour always displays the
confidence level and further
contours increase by factors of
from there.
The angular resolution and sensitivities (1
)
of the
maps are listed in Table 3. We also include in
Table 3 the measured total 1.43 GHz flux densities,
S1.43, and the position angles of the radio continuum emission
distribution, as determined from fitting two-dimensional ellipsoidal
Gaussians to the radio continuum emission distributions, using the
AIPS task JMFIT.
NGC1055: The 1.425 GHz radio continuum emission of NGC1055 is displayed
in Fig. 1. It is centrally peaked, with resolved
disk emission out to radii of about
.
As in the subsequent
figures, the full width at half maximum (FWHM) of the beam is
displayed in the lower left corner. There are several nearby point
sources, which are most likely unrelated.
NGC1406: Like NGC1055, NGC1406 (Fig. 2) exhibits centrally peaked emission from its disk, however with little substructure at the angular resolution of the current data.
NGC1421: NGC1421 is resolved by our observations, which were combined with data by Irwin et al. (1999). The emission is centrally peaked, but the central maximum is not very bright (Fig. 3). Several nearby unrelated point sources had to be removed for investigations of its structure perpendicular to the disk plane.
NGC1511: NGC1511 (Fig. 4) has an unusual radio continuum distribution in its disk. There are two prominent maxima, about 30'' apart. The surface brightness of the emission is high at both 2.45 GHz and 1.43 GHz. Perpendicular to the disk plane the emission distribution is very thick. The total height of the emission distribution at 2.45 GHz is about 5 beamwidths (about 4 at 1.43 GHz). The axial ratio in the optical is 2.9 (NASA Extragalactic Database; NED), while the radio axial ratio is approximately 1.4 (at 2.45 GHz) and even 1.2 at 1.43 GHz.
None of the radio continuum emission from NGC1511 is associated with SN1935C (or "Nova'' Hyi 1935; van den Bergh & Hazen 1988), which occurred further to the north-east. Also the two sources about 2' south of NGC1511 are probably unrelated.
NGC2748: The radio continuum emission distribution of NGC2748, displayed in Fig. 5, is centrally peaked but shows almost no significant substructure at the current resolution. The source to the west of NGC2748 is unrelated.
NGC2820: Figure 6 shows a complex emission distribution, comprised of contributions from NGC2820 (east of the centre) and NGC2814 (to the west) as the major peaks and NGC2820A as a tertiary maximum in between these two. Fitting Gaussians to the emission peaks of these sources, a fourth component remains that is associated with the companion of NGC2814. Another member of this galaxy group, NGC2805 (outside the displayed field of view) was also detected.
This proximity of the different emission sources makes it difficult to disentangle the contributions from each of the systems. Using Gaussian ellipsoids to fit the individual emission distributions, we succeeded in separating the radio emitters. However, the same problem of flux contamination also applies to the FIR data, based on which the target was selected (see Table 4), where the sources cannot be spatially separated due to the lower angular resolution of IRAS compared with our data.
NGC3175: The radial extent of NGC3175's radio continuum emission distribution (Fig. 7) is small compared to its optical extent (Condon et al. 1996). The emission is centrally peaked and marginally resolved along the disk plane, with no directly visible substructure. The source to the north-west is unrelated.
NGC3437: NGC3437 (Fig. 8) shows little substructure in its radio continuum emission distribution at the present resolution.
NGC3717: The radio continuum emission distribution of NGC3717 is complex (Fig. 9). Besides a strong central peak extended disk emission is visible as an emission ridge along the major axis, with a secondary emission peak about 3' south-west of the centre and another, weaker source about 3' north-east of the nucleus. The total radial extent of the emission distribution is about 7'. South-east of the centre, further emission is visible away from the disk plane, which is not easily removable by subtracting a point source profile and might therefore be genuinely extended. The weak point source north-west of the centre is probably unrelated.
![]() |
Figure 2:
VLA 1.43 GHz radio continuum contour map of NGC1406, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.23, -0.16, 0.16 (= |
| Open with DEXTER | |
![]() |
Figure 3:
VLA C+D array 1.465 GHz radio continuum contour map of
NGC1421, superimposed on a Digital Sky Survey (DSS) image.
The contour levels displayed are -0.4, -0.28, -0.2, 0.2
(= |
| Open with DEXTER | |
NGC4527: The radio continuum emission of NGC4527, as displayed in
Fig. 10, is bright, centrally peaked and extended
along its major axis, with an extent of about
.
The
ratio of major-to-minor axis extent is less than 2, which is low
for a highly inclined galaxy. The fact that the minor axis extent
of the emission distribution is about 4 beamwidths indicates that
it is resolved in this direction, despite the fact that there is
little visible substructure. The central source, which is
classified as an H II/LINER nucleus (NED), was subtracted for
further processing (see Sect. 3.2).
NGC4700: Despite its small angular extent, the radio continuum emission distribution of NGC4700 (Fig. 11) is resolved in both the major and minor axis direction. The axial ratio of the emission distribution is only 1.25, which is extremely low. This is a clear indication of extraplanar emission, especially when the elliptical shape of the beam, with its major axis aligned with that of the galaxy, is taken into account.
Two strong point sources east of NGC4700 (one just outside the
displayed field of view) are unrelated. A third, weaker point
source about
south-west of the galaxy centre might
possibly be related with NGC4700, being located near its disk
plane, however beyond the radial limit of the optical emission
distribution. This point source was not included in the total
flux measurement and when determining the position angle, PA (Table 3).
NGC5073: The emission of NGC5073 is unresolved (Fig. 12), with no immediately visible sign of extraplanar emission at the angular resolution of the current data.
NGC7090: NGC7090 (Fig. 13) exhibits centrally peaked
radio continuum emission from its disk. The radial extent of the
disk emission is about 4' at 2.45 GHz and
at 1.43 GHz. The disk emission is clearly resolved, with secondary maxima
most clearly visible in the 2.45 GHz map.
There are also clear indications of extraplanar emission from cosmic
ray electrons (CRs) in the halo. The halo emission is detected and
spatially resolved at both observing frequencies, 2.45 GHz and 1.43
GHz. The 1.43 GHz halo emission is apparently more extended and
clearly brighter than at 2.45 GHz.
![]() |
Figure 4:
ATCA radio continuum maps of NGC1511. Left panel: 2.45 GHz map;
right panel: 1.43 GHz map; both are overlaid on a DSS optical
image. The contours display the -4, -2.8, -2, 2, 2.8, 4, ...
|
| Open with DEXTER | |
![]() |
Figure 5:
VLA 1.43 GHz radio continuum contour map of NGC2748, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.4, -0.28, -0.2, 0.2 (= |
| Open with DEXTER | |
![]() |
Figure 6:
VLA 1.43 GHz radio continuum contour map of NGC2820, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.17, -0.12, 0.17, 0.24, ..., 30.7 mJy beam-1, with
a contour spacing of a factor of |
| Open with DEXTER | |
| Galaxy | Obs. | Radiala | Resolution | Disk | Halo | ||||
| Freq. | Range |
|
Rel. |
|
Rel. | ||||
| (GHz) | (''/kpc) | (''/kpc) | (''/kpc) | Ampl. | Side1b | Side2b | Average | Ampl. | |
| NGC1055 | 1.425 | 54.0/4.19 | 19.2/1.49 | 1.00 | -- | -- | -- | 0.00 | |
| NGC1406 | 1.425 | 43.0/5.28 | 4.8/0.35 | 0.85 | 27.3/1.97 | 18.6/1.34 | 23.0/1.66 | 0.15 | |
| NGC1421 | 1.465 | 37.5/5.66 | 13.3/2.01 | 1.00 | -- | -- | -- | 0.00 | |
| NGC1511 | 1.43 | 33.0/2.80 | 10.5/0.89 | 0.96 | 14.9/1.27 | 26.6/2.26 | 20.7/1.76 | 0.04 | |
| 2.45 | 17.7/1.50 | 7.9/0.67 | 0.90 | 17.9/1.52 | 18.1/1.54 | 18.0/1.53 | 0.10 | ||
| NGC2748 | 1.425 | 52.0/7.24 | 12.5/1.74 | 1.00 | -- | -- | -- | 0.00 | |
| NGC2820 | 1.425 | 41.0/5.95 | 10.8/1.57 | 1.00 | -- | -- | -- | 0.00 | |
| NGC3175 | 1.425 | 33.0/2.55 | 10.3/0.80 | 0.94 | 20.3/1.57 | 0.0/ 0.0 | 10.2/0.79 | 0.06 | |
| 1.49c | 15.0/1.16 | 5.3/0.41 | 0.84 | 27.6/2.13 | 14.4/1.11 | 21.0/1.62 | 0.16 | ||
| NGC3437 | 1.425 | 41.0/5.07 | 9.8/1.21 | 1.00 | -- | -- | -- | 0.00 | |
| NGC3717 | 1.425 | 41.0/5.39 | 14.2/1.87 | 1.00 | -- | -- | -- | 0.00 | |
| NGC4527 | 1.425 | 50.0/7.74 | 20.8/3.22 | 1.00 | -- | -- | -- | 0.00 | |
| NGC4700 | 1.425 | 46.0/5.69 | 0.1/0.01 | 0.33 | 26.9/3.33 | 22.5/2.78 | 24.7/3.06 | 0.67 | |
| NGC5073 | 1.425 | 44.0/9.37 | 14.7/3.13 | 1.00 | -- | -- | -- | 0.00 | |
| NGC7090 | 1.43 | 34.5/1.96 | 7.0/0.40 | 0.58 | 24.5/1.39 | 38.0/2.16 | 31.2/1.77 | 0.42 | |
| 2.45 | 17.7/1.00 | 2.5/0.14 | 0.59 | 21.0/1.19 | 29.0/1.65 | 25.0/1.42 | 0.41 | ||
| NGC7462 | 1.43 | 46.5/3.41 | 10.0/0.73 | 0.65 | 25.3/1.85 | 24.4/1.79 | 24.8/1.82 | 0.35 | |
| 2.45 | 23.5/1.72 | 2.5/0.18 | 0.66 | 26.2/1.92 | 20.6/1.51 | 23.4/1.71 | 0.34 | ||
| NGC7541 | 1.425 | 47.0/9.67 | 9.8/2.02 | 1.00 | -- | -- | -- | 0.00 | |
![]() |
Figure 7:
VLA 1.43 GHz radio continuum contour map of NGC3175, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.2, -0.14, 0.14 (= |
| Open with DEXTER | |
![]() |
Figure 8:
VLA 1.43 GHz radio continuum contour map of NGC3437, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.28, -0.2, 0.2 (= |
| Open with DEXTER | |
NGC7462: Figure 14 displays centrally peaked emission in
the disk of NGC7462 at both 2.45 GHz and 1.43 GHz. The 2.45 GHz
map exhibits a secondary emission maximum about
east
of the centre and a weaker corresponding maximum to the west of
the nucleus.
The 2.45 GHz map also shows some emission south of the centre, outside the disk of NGC7462. At 1.43 GHz one can discern extended emission from the halo. Despite the different angular resolutions of the two maps it appears that the halo is more prominent at 1.43 GHz than at higher frequencies. The 2.45 GHz map may also lack sensitivity and thus not show the full extent of the emission at this frequency yet. Note that the total flux densities of this galaxy at 2.45 GHz and 1.43 GHz are quite low (Table 3).
NGC7541: NGC7541 (Fig. 15) was observed in bad weather conditions, which caused phase problems that could not be removed entirely from the visibility data using the self-calibration technique. The 1.43 GHz emission distribution is centrally peaked and marginally resolved along the major axis, without visible substructure at the current resolution. The source to the south-west of NGC7541 is NGC7537. The other nearby point sources are probably also unrelated.
In order to investigate the z structure of the sample galaxies (i.e., the emission distribution perpendicular to the disk plane), the following technique was applied (as done by us earlier; see e.g. Hummel et al. 1991b; Dahlem et al. 1994; Dahlem et al. 1997):
![]() |
Figure 9:
VLA 1.43 GHz radio continuum contour map of NGC3717, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.16, 0.16 (= |
| Open with DEXTER | |
![]() |
Figure 10:
VLA 1.43 GHz radio continuum contour map of NGC4527, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.42, -0.3, 0.3 (= |
| Open with DEXTER | |
![]() |
Figure 11:
VLA 1.43 GHz radio continuum contour map of NGC4700, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.28, -0.2, 0.2 (= |
| Open with DEXTER | |
NGC7090: The relatively high resolution of our data of NGC7090 allows
for a clear distinction between disk and halo emission in this
galaxy.
Therefore, we use it here as a showcase to demonstrate the
process of modeling the z distribution of the emission.
The observed emission distribution is reproduced on a linear
scale as the solid line in the left panel of Fig. 16. An exponential disk model was created
(narrow peaked distribution, solid line) and convolved with the
beam (bold line), leading to the smoother profile (dotted line).
Its amplitude is initially fixed at 1 to see whether it can fit
the observed emission distribution.
However, deviations from the disk model set in at about 25-30%
of the peak surface brightness. In such cases, when residuals
remain, a second, wider exponential is added and beam-convolved.
The two components are then summed and the resulting function
(dashed line) fitted to the data.
For easier inspection of the exponential wings, the data and
model, including its individual components, are then reproduced
on a logarithmic scale (right panel of Fig. 16).
One can see that the halo emission is asymmetric, with
brighter emission towards positive offsets (the north-east)
than south-west of the disk. To save space, only logarithmic
plots are shown for all other galaxies.
Deviations from a purely exponential decay are visible at 2.45
GHz (Fig. 17). These might be caused by the
spurs with relatively high surface brightness detected by
Harnett & Reynolds (1985). Such details are visible, because
the 2.45 GHz map of NGC7090 is not only one of two maps with
the highest angular resolution in the sample, but-because of
the proximity of NGC7090-also the spatial resolution is good
(1 kpc).
The peak surface brightness of the averaged z profile at
2.45 GHz is only 4 mJy beam-1 and this map is accordingly the
one with the smallest dynamic range (Fig. 17).
Thus, the data used for calculating
go
down only to the 5%-level (log = -1.3) and the influence
of noise is visible in the wings of the emission distribution.
![]() |
Figure 12:
VLA 1.43 GHz radio continuum contour map of NGC5073, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.6, -0.42, -0.3, 0.3 (= |
| Open with DEXTER | |
![]() |
Figure 13:
ATCA radio continuum maps of NGC7090. Left panel: 2.45 GHz
map; right panel: 1.43 GHz map; both are overlaid on a DSS
optical image. The contours display the -4, -2.8, -2, 2,
2.8, 4, ... |
| Open with DEXTER | |
![]() |
Figure 14:
ATCA radio continuum maps of NGC7462. Left panel: 2.45 GHz
map; right panel: 1.43 GHz map; both are overlaid on a DSS
optical image. The contours display the -4, -2.8, -2, 2,
2.8, 4, ... |
| Open with DEXTER | |
![]() |
Figure 15:
VLA 1.43 GHz radio continuum contour map of NGC7541, superimposed
on a Digital Sky Survey (DSS) image. The contour levels displayed
are -0.6, -0.42, -0.3, 0.3 (= |
| Open with DEXTER | |
The values of the "disk scale height'',
,
in
Table 4 represent the apparent exponential scale
height of the thin disk emission component. At the resolution
of our data, and given the uncertainties in the values of the
inclination angles, i, it is not clear what causes the apparent
thickness of the thin disk component in each galaxy. The total
width could be due to either an inclination of
or intrinsic disk thickness, or a superposition of both.
The
values are thus not in all cases a direct
measure of the intrinsic disk thickness.
A comparison of the best-fitting model of a thin disk with the observed cut perpendicular to the disk plane tells immediately whether there is any residual emission that cannot be explained by the single-component model or not. Any emission at high projected z offsets that is not well represented by the best-fitting single exponential (dotted line) must arise from outside the galaxy disk. In the presence of significant amounts of residual emission, especially in the form of broad wings in the emission distribution perpendicular to the disk plane, it is assumed that a radio halo exists. For all practical purposes, it is assumed that the galaxy disks are not warped in those parts from where the observed radio continuum emission arises. In a number of cases there is only marginal or no evidence for the existence of extraplanar synchrotron radio continuum emission. Some of our maps do not have sufficient resolution to draw firm conclusions yet. This affects in particular the most distant objects in the sample (NGC5073 and NGC7541), but also a few others.
In those cases where there is excess emission beyond the disk,
the scale height of the second, wider exponential component
is interpreted as the exponential scale height of the halo
emission. It has been noted by different authors before that
the surface brightness of radio halos declines exponentially
away from the galaxy disks (see, for example, the high-resolution
z profile through NGC891 by Dahlem et al. 1994).
Therefore, the halo thickness can be quantified quite naturally
by the exponential scale height of the halo emission,
,
as listed in Cols. 7-9 of Table 4.
All other relevant parameters, including the relative amplitudes
of the exponential components fitted and the radial ranges
(Col. 3) over which the fits to the disk and halo emission
components were performed, are tabulated there as well.
Note that the relative amplitude of the two exponential components
is no indication of the actual total flux percentages, because the
z profiles were performed only over a limited radial range,
excluding the outer parts of the galaxy disks.
Data for the fits from which the exponential z scale heights were determined (as listed in Table 4) were used only above the noise level of the images (Figs. 18-33). The threshold in all cases lies in the range from 1-5% of the peak surface brightness, depending on data quality and source brightness. In the following we will present a short description of the observed z structure of all galaxies in the sample (except NGC7090, which was already introduced above).
![]() |
Figure 16: Cut through the 1.43 GHz emission of NGC7090, perpendicular to the disk plane, on a linear scale (left) and logarithmic (right). In the linear graph the observed data (solid line), beam profile (bold line), exponential disk model and exponential halo model (solid lines, peaked distributions), the two beam-convolved exponential components (dotted lines) and the total model (sum of the two exponential components; dashed line) are presented. The logarithmic plot shows the same components, except the initial unconvolved exponentials. |
| Open with DEXTER | |
NGC1055: There is no evidence in Fig. 18 for the presence of emission beyond the thin disk of NGC1055. The beam-convolved exponential disk model is a good representation of the data. Therefore, the dotted line of the single-exponential component model and the dashed one of the total model overlap exactly, resulting in a dash-dotted line that in turn conincides in large parts with the solid line representing the data.
NGC1406:
Although relatively weak, the emission beyond the best-fitting
disk model in the logarithmic plot in Fig. 19 is
statistically significant and can be fitted only by adding a
second exponential component. Higher angular resolution is
required to confirm this preliminary result and to allow
for a more reliable calculation of
.
NGC1421: The z profile of NGC1421 (Fig. 20) does not give any indication of the presence of wings that cannot be fitted by a single-exponential model.
NGC1511: Figures 21 and 22 display the z profiles at 2.45 GHz and 1.43 GHz, respectively. Both provide evidence for the existence of extraplanar radio continuum emission. The slopes of the halo emission are about linear in the logarithmic plots, indicating the expected exponential decay. Although the relative amplitude of the halo emission at 1.43 GHz is only 4% of the peak, it is clearly detected (Table 4). The significantly lower amplitude of the halo component at 1.43 GHz compared to 2.45 GHz is most likely caused by beam smearing of disk emission.
NGC2748: Figure 23 illustrates how precisely the emission profile of NGC2748 can be approximated by a single beam-smeared exponential distribution.
NGC2820: No proof of excess emission beyond the best-fitting disk model is visible in our data of NGC2820 at the present resolution (Fig. 24).
NGC3175: Although not very extended and at a low level, deviations from the best-fitting disk model were found on the south-west side of NGC3175, which are modeled by an additional exponential component displayed in Fig. 25. In Sect. 4.2.1 radio data with a higher angular resolution are presented that clearly confirm the presence of a radio halo.
![]() |
Figure 17: Cut through the 2.45 GHz emission of NGC7090, perpendicular to the disk plane (same line types as the logarithmic plot in Fig. 16). |
| Open with DEXTER | |
![]() |
Figure 18: Cut perpendicular to the disk plane of NGC1055, perpendicular to the disk plane (same line types as in Fig. 17). |
| Open with DEXTER | |
NGC3437: There is no evidence for the existence of extraplanar emission from NGC3437 from the z profile in Fig. 26. Some low surface brightness excess emission is visible on the north side (positive offsets). However, the result is not clear enough to be considered a secure detection because of residual phase uncertainties in the data. Observations with higher resolution and better signal-to-noise ratio are warranted.
NGC3717: NGC3717 exhibits no deviations from the best-fitting disk model in Fig. 27, but note the large scale height of the emission distribution compared to those galaxies where disk and halo emission could be separated.
NGC4527: In the case of NGC4527 insufficient angular resolution is not the primary problem in the detection of halo emission (see Fig. 28). However, it has the widest apparent disk emission distribution of all objects in the sample, probably due to its relatively low inclination (see Col. 5 in Table 3). No significant emission beyond the disk model is detected. The slope of the measured data is partly even steeper than that of the model (especially on the north side), indicating that there is either a residual data problem or that the emission distribution falls off more steeply with increasing z-distance than an exponential.
NGC4700: The most convincing new detection of a radio halo amongst the
galaxies observed with the VLA D array has been made in NGC4700,
a starburst galaxy that had scarcely been studied prior to our
observations. The logarithmic z profile in Fig. 29
exhibits wide linear wings on both sides of the galaxy disk,
indicating the presence of an exponentially decaying CR halo
that can be traced reliably out to the 1.25% level, up to
(10 kpc) away from the disk plane, with an unusually
high scale height of 3.1 kpc. Deviations from the best fit to
the disk emission start at a very high surface brightness. In
fact, NGC4700 is the only galaxy in the present sample in
which the halo exponential has a higher amplitude than the
disk component. The scale height of the disk component is not
well-constrained because it is unresolved.
NGC5073: This is the most distant object in our sample. No halo emission could be detected from this galaxy and the z profile of the disk is also almost unresolved by the beam (Fig. 30). The fact that the observed profile crosses the beam profile might indicate residual phase problems in the data or a slope that is steeper than exponential.
![]() |
Figure 19: Cut perpendicular to the disk plane of NGC1406 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 20: Cut perpendicular to the disk plane of NGC1421 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 21: Cut through the 2.45 GHz emission of NGC1511, perpendicular to the disk plane (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 22: Cut through the 1.43 GHz emission of NGC1511, perpendicular to the disk plane (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 23: Cut perpendicular to the disk plane of NGC2748 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 24: Cut perpendicular to the disk plane of NGC2820 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 25: Cut perpendicular to the disk plane of NGC3175 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 26: Cut perpendicular to the disk plane of NGC3437 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 27: Cut perpendicular to the disk plane of NGC3717 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 28: Cut perpendicular to the disk plane of NGC4527 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 29: Cut perpendicular to the disk plane of NGC4700 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 30: Cut perpendicular to the disk plane of NGC5073 (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 31: Cut through the 2.45 GHz emission of NGC7462, perpendicular to the disk plane (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 32: Cut through the 1.43 GHz emission of NGC7462, perpendicular to the disk plane (same line types as in Fig. 17). |
| Open with DEXTER | |
![]() |
Figure 33: Cut perpendicular to the disk plane of NGC7541 (same line types as in Fig. 17). |
| Open with DEXTER | |
NGC7462: Clear signatures of extraplanar emission are detected in the z profiles through NGC7462 at both observing frequencies, 2.45 GHz (Fig. 31) and 1.43 GHz (Fig. 32). Also at both frequencies an asymmetry of the halo is apparent, the brighter emission being observed on the south side (towards negative offsets).
NGC7541: This galaxy, at almost the same distance as NGC5073, is also unresolved along its minor axis by our observations, as displayed in Fig. 33. No excess emission is observed beyond the exponential disk model.
| Galaxy | 2.45 GHz |
|
| flux (mJy) | (
|
|
| NGC1511 |
|
|
| NGC7090 |
|
|
| NGC7462 |
|
From the VLA observations, and thus for 12 of the 15 galaxies,
no spectral index information is available. However, the ATCA
observations at two frequencies simultaneously do provide as
a by-product spectral indices,
.
We have
derived mean spectral indices,
,
over
the entire galaxies, as tabulated in Table 5.
Mean values around
,
as observed in the three
galaxies NGC1511, NGC7090 and NGC7462, are quite normal
for late-type spiral galaxies (e.g., Lisenfeld & Völk 2000
and references therein).
Images of the spectral index distribution within the three galaxies
are displayed in Figs. 35-37.
All three displays are presented in the same fashion, i.e., with
the same grey scales and contours (see figure captions), making
them directly comparable. The angular resolution of the spectral
index maps is that of the 1.43 GHz data. Spectral indices are
displayed only in those areas where the signal-to-noise ratio
in both input images (1.43 and 2.45 GHz) is ![]()
.
The boundaries surrounding the spectral index distributions
indicate that the maps were blanked outside this area.
The spectral index distribution in NGC1511 (as displayed in Fig. 35) is quite homogeneous, with a few edge effects (especially at the southern and northern end of the distribution) that should not be taken at face value.
In NGC7090 (Fig. 36) the flattest spectrum (i.e., most visibly influenced by thermal emission) is observed in the central part of the disk, with a slight steepening towards the outer disk and the halo regime.
NGC7462 (Fig. 37) exhibits a spectral index steepening slightly away from the disk plane, with little substructure at the resolution of the current data.
A general trend is visible that the spectral indices steepen slightly away from the disk planes, indicating that energy losses of the CRs play an important role.
The maps in Figs. 1-15 show only total emission because, with the current sensitivity and angular resolution, no polarised radiation was detected (ATCA data) or calibration not performed (VLA observations). Thus, no information is available on magnetic field configurations in the observed galaxies. This implies that any possible effect that the magnetic field configuration might have on the observed CR distribution in the halos cannot be investigated here.
![]() |
Figure 34: Cut perpendicular to the disk plane of NGC3175 (same line types as in Fig. 17). Based on VLA C array data by Condon et al. (1996). |
| Open with DEXTER | |
![]() |
Figure 35: 2.45 GHz vs. 1.43 GHz spectral index distribution within NGC1511. The grey-scale ranges from -2 to 1; the contours start at -2, increasing by 0.2 each. The central oval contour represents a spectral index of -0.8. |
| Open with DEXTER | |
![]() |
Figure 36: 2.45 GHz vs. 1.43 GHz spectral index distribution within NGC7090. The grey-scale ranges from -2 to 1; the contours start at -2, increasing by 0.2 each. The central closed contour represents a spectral index of -0.6. |
| Open with DEXTER | |
![]() |
Figure 37: 2.45 GHz vs. 1.43 GHz spectral index distribution within NGC7462. The grey-scale ranges from -2 to 1; the contours start at -2, increasing by 0.2 each. The central contours closest to the disk mid-plane represent a spectral index of -0.6. |
| Open with DEXTER | |
| Galaxy | Radio | H |
| NGC1055 | no | -- |
| NGC1406 | yes | -- |
| NGC1421 | no | no |
| NGC1511 | yes | yes |
| NGC2748 | no | -- |
| NGC3175 | yes | no |
| NGC3437 | no | -- |
| NGC3717 | no | -- |
| NGC4700 | yes | yesa |
| NGC7090 | yes | yesa |
| NGC7462 | yes | yesa |
The absence of clear evidence for the presence of halos in the remaining galaxies can, as indicated above, have different reasons:
In addition to the synchrotron radio continuum from relativistic
CR electrons there are other tracers of gaseous halos around
spiral galaxies, of which the most important are H
recombination radiation of warm ionised gas and soft X-ray
emission from a hot thermal plasma (cf. Dahlem 1997).
In the following we briefly summarise those previously published radio continuum images that are relevant for searches of extraplanar emission in the galaxies presented here.
NGC1421: Irwin et al. (1999) found some evidence for extraplanar emission in NGC1421 from their VLA C array 1.4 GHz radio continuum image (see their Fig. 4). The authors provided us with their data, which were combined with our D array visibilities to obtain a map that contains the full flux and has improved angular resolution (as displayed in Fig. 3).
NGC2748: A 15'' resolution 1.465 GHz map of NGC2748 by Hummel et al. (1985) might indicate the presence of radio continuum emission away from its disk plane (their Fig. 1c). However, these authors quote a lower total flux than measured by us. Due to missing short spacings (and thus missing extended flux) these data cannot be used for the present study.
NGC2820: The VLA C array image by Hummel & van der Hulst (1989) has a
higher angular resolution than ours (15'' [their Fig. 1]
compared to
;
Fig. 6, above),
but lower sensitivity. Based on a z profile as used by us here
(their Fig. 2), Hummel & van der Hulst (1989) found that
NGC2820 has a thick radio disk or halo at 1.465 GHz. Since
these authors did not report an exponential scale height for
the halo emission, we rely here on our own measurement (note
that this galaxy has been removed from the sample anyway).
The fact that Hummel & van der Hulst find evidence for the
presence of a halo, while there is none in our data,
demonstrates the need for sufficient angular resolution.
NGC3175: A 1.49 GHz radio continuum map of NGC3175 by Condon et al. (1996)
displays three maxima in the galaxy disk, one at the centre and
one each about 45'' on either side of it, thus coinciding with
the most prominent H II regions in the H
image by Ryder
& Dopita (1994). The contour map displayed in Condon et al. (1996) might
also indicate halo emission above its central region. This VLA C
array map has a higher angular resolution (15'' compared to
)
than ours, but lower sensitivity. It exhibits
a low axial ratio of the radio continuum emission, suggesting that
there might be a relatively small outflow from the nuclear region.
Since it recovers the total flux (the same within the error margins
as our measurement; see Table 3), Condon et al.'s image
can be used for a study of the z structure of the emission
distribution.
Because of the very different integration times, it is not possible to combine Condon et al.'s data with ours. However, the image being available via NED, we have retrieved it and used it to produce a second z profile. This higher-resolution z profile is displayed separately in Fig. 34. The plot reassuringly reproduces the major features from Fig. 25, however in more detail. The distribution of excess emission beyond the thin disk is slightly asymmetric, being brighter on the south-east side (negative offset values), as already visible in our DnC array data. But now it is much clearer that this excess emission, in particular also that on the north-west side, is indeed significant. The results from the fitting process are collected in Table 4.
NGC3437: Condon et al. (1990) presented a high-resolution 1.49 GHz image of NGC3437. However, since it does not recover the total flux of the galaxy, it is not suitable for investigations of the extended emission distribution, including the z structure perpendicular to the disk plane.
NGC5073: We note here that a 1.49 GHz radio image of NGC5073 by Condon
et al. (1990) exhibits a compact central source that is unresolved
at a resolution of about 5''. This is relevant in the context
of the H
data of this galaxy mentioned below. Since it
does not recover the full extended flux of the galaxy, the radio
map cannot be used for a study of the large-scale emission
distribution along the minor axis.
NGC7090: Two "spurs'' of radio emission above the disk of NGC7090 were detected by Harnett & Reynolds (1985) at 843 MHz (see their Plate 3). Our maps are more sensitive and have better angular resolution thus showing more extended halo emission. The brightest spur north-east of the central region, now embedded in extended emission, is visible in our 1.43 GHz map (cf. Fig. 13) as kinks in the isophotal contours, which can be traced from close to the disk plane out to the detection limit of the halo.
NGC7462: There is a tentative detection of extraplanar radio continuum emission from NGC7462 in a 1.49 GHz VLA map by Condon et al. (1987). However, the angular resolution (of 1') and sensitivity of these snapshot observations were insufficient to make a firm statement. Our detection is much more conclusive, in particular because the result comes from observations at two independent frequencies.
NGC7541: Condon et al. (1990) published a snapshot 1.49 GHz map of NGC7541, with a resolution of 15''. It exhibits a dominant unresolved central source, which might hint at the presence of an active nucleus. Since this radio image does not recover all extended source flux, it is not used as part of the present investigation.
5 of the 15 galaxies studied here are also in the sample by Lehnert
& Heckman (1995): NGC1511, NGC2820, NGC4527, NGC5073 and
NGC7541. Three of these, namely NGC1511, NGC4527 and NGC5073,
were imaged in H
line emission. In addition, H
images of a few more individual galaxies can be found in the
literature.
NGC1421: An unpublished H
image by Dettmar (priv. comm.) exhibits
numerous H II regions in the disk of NGC1421, indicating a
two-arm spiral structure and an inclination of
.
No extraplanar H
line emission was detected. However,
this might be due to a lack of sensitivity.
NGC1511: The H
image of NGC1511 by Lehnert & Heckman (1995) shows a
disturbed structure, with H II regions that are located outside the
thin disk plane. There are also indications of diffuse H
,
which apparently does not arise from the central disk plane, but
there is no direct evidence for halo emission. The disturbance of
its disk might arise from interactions with one of its partners,
NGC1511A or B, or just be its intrinsic structure. With an absolute R magnitude of -18.6 (Lehnert & Heckman 1995) NGC1511 is a low
optical luminosity, low-mass galaxy, in which disturbances are known
to be quite normal (compare, e.g., with NGC1569, NGC4449 or
NGC2188).
NGC3175: An H
image by Ryder & Dopita (1994) shows bright H II regions at the position of the nucleus and what appears to be a
starburst ring with a radius of 3.5 kpc around it. The H II regions at the tangential points of this ring are apparently
the brightest. These might also be "thick'', i.e. extended in
the direction perpendicular to the disk plane. However, no
extended diffuse H
emission is visible, especially
away from the disk plane. Note that NGC3175 does not exhibit
any signs of ongoing star formation outside the starburst ring.
Its outer disk has a very red colour (Ryder & Dopita 1994)
and is H I deficient (Dahlem et al. 2001).
NGC4527: The H
image of NGC4527 by Lehnert & Heckman (1995)
shows no obvious signs of extraplanar emission. Instead, it
indicates that the disk is not inclined exactly edge-on. From
its optical axial ratio of 3 Tully (1988) derived an inclination
angle of about
,
which has been taken into account
in our representation of the radio emission from the disk in
Fig. 28 and which is the reason for removing
the object from the present sample.
NGC4700: NGC4700 was observed by us (Rossa 2001) and the resulting image
will be presented separately. Clear evidence for extraplanar
H
emission is visible in our frame; therefore the
corresponding entry was made in Table 6.
NGC5073: In the H
image of NGC5073 by Lehnert & Heckman (1995)
one can see little diffuse emission above the nuclear area and a
strong point source in the centre. Together with the unresolved
nuclear radio continuum source, this might hint at the presence
of a previously unclassified active nucleus, which could explain
why a high
f60/f100 flux ratio is observed and the radio
emission is unresolved by our data, while there is no indication
of a halo.
NGC7090 and NGC7462: Both galaxies were observed and extraplanar H
emission
clearly detected (Rossa 2001). These results will be presented
at a later time, together with those on NGC4700.
Resume:
The small number of H
detections is probably in most cases
due to sensitivity problems, because the images available in the
literature were not taken for dedicated searches of low surface
brightness extraplanar emission. Long integration times are
required for attempts to detect such emission.
Data existing in archives. ROSAT archival data exist for NGC1406, NGC2820, NGC4527; however, they do not contribute anything new to the current investigation. NGC1406 is detected, but unresolved; a by-eye inspection of the data reveals no detected emission from either NGC2820 or NGC4527.
Einstein observations of NGC1421 show weak emission from its disk, but also from a "tail'' to the south of the optically visible disk and also west of the nuclear area, thus perpendicular to the disk plane, where the extreme colour stretch used by Fabbiano et al. (1992) for the underlying optical image might suggest the presence of a stellar tidal arm.
Future observations. Of the present sample, Chandra observations of NGC1055 have been approved. We will obtain XMM-Newton Guaranteed Time observations of NGC1511 and NGC5073, which will be presented elsewhere in due course.
Considering all observations listed above of tracers of gaseous
halos, 6 out of the 11 galaxies presented here show signs of
extraplanar emission in at least one waveband. Different
components of the halo ISM in the target galaxies were detected
in four cases, as previously done in more nearby systems (see
e.g. Dahlem 1997). A list of the detection/non-detection of
extraplanar radio and H
emission is provided in
Table 6.
Thus, it appears that the existence of gaseous halos around late-type spiral galaxies is indeed not as uncommon as it appeared in the past. Rather, the search criteria are now optimised. Also the high detection rate by Irwin et al. (1999) of >90% fits into this scenario.
FIR brightness and in particular "warm'' FIR colours, i.e. high f60/f100 flux ratios, have turned out to be a very effective search criterion (Heckman et al. 1990; Lehnert & Heckman 1995). This is also born out by the fact that all FIR-warm nearby edge-on starburst galaxies studied by us previously in the soft X-ray regime (D98) exhibit diffuse extraplanar X-ray emission. A somewhat lower detection rate in the present sample can be expected, because we include here objects at larger distances and thus with lower flux densities and surface brightnesses than before. With values of in some cases 5-9 kpc (Table 4) also the spatial resolution of our data is not always sufficient yet.
As outlined by Dahlem (1997), several processes can contribute to the creation of gaseous halos around late-type spirals. In order to make sure that only star-formation related energy input is taken into account, galaxies with AGNs and closely interacting systems have been removed from the original sample above. Many starburst galaxies are interacting, or reversely, there is an overabundance of starbursts in interacting galaxies (cf. Hummel et al. 1990; Lutz 1992). In the case of distant encounters one can hope that the interaction has only caused minor disturbances of the gas in the galaxies, leading to the accretion of gas near their inner Lindblad resonance (Combes 1987) or turnover point of galactic rotation (Lesch et al. 1990) and subsequently to the onset of the observed starbursts, while the starbursts themselves dominate the energy balance at the time of the observation. This leaves galaxies that have high f60/f100 flux ratios because of SF-related heating of their dust, with approximately unperturbed gravitational potentials. SF-related heating sources are photo-ionisation, shocks and/or turbulent mixing layers, but in principle all energy comes from processes relating to the evolution of high-mass stars (winds and type II supernovae; Leitherer & Heckman 1995; Leitherer et al. 1999). The above points out that careful selection of candidates and interpretation of the results is crucial; too many free parameters would leave doubts about the source of the energy driving the gaseous halos. Based on the above selection, we are confident that halo emission detected in the remaining 11 galaxies of our sample arises from energy input into the ISM from high-mass stars and type II SNe.
Having ensured that SF dominates the energy input into the disk ISM, the next step in an investigation of a potential connection between halo properties and SF activity in the underlying disk is to quantify the level of energy input and the halo properties.
In earlier studies of halo properties (e.g. Hummel et al. 1991b; Dahlem et al. 1994; Dahlem et al. 1997) we have used z profiles averaged over a broad radial range to quantify the properties of radio halos. This has the advantage of increasing the signal-to-noise ratio compared to a single-pixel profile. At the same time, data are only averaged over regions where emission in the halo has clearly been detected, adapting the measurements to the individual target's geometry.
Ellipsoidal fitting to radio images, as performed by Irwin et al. (1999) makes sense under the assumption that radio halos are
engulfing their host galaxies. However, this is not observed
in images with sufficient angular resolution. Images that might
suggest such behaviour (like, for example, our 1.43 GHz map of
NGC1511 in Fig. 4) often suffer from insufficient
angular resolution. Instead, with good resolution, a behaviour
as described by us earlier (DLG95)
is observed in most galaxies: radio halos exist only above the
brightest emission regions of the underlying galaxy disks, i.e.
the active SF regions. The radial extent of the halos is similar
to the radius out to which SF is observed,
,
and thus
smaller than that of the entire disks. This is evident in our new
1.43 GHz and 2.45 GHz images of NGC7090 (Fig. 13).
Supporting evidence for this result comes from recent H
observations by Rossa & Dettmar (2000). Ellipses do not fit the
geometry of the halo plus disk emission properly. We find that
the most reliable technique to detect and quantify extraplanar
emission is by producing averaged z profiles.
Based on radio data of the nearby edge-on galaxies NGC891 and
NGC4631, and - with less detailed information - on a few other
galaxies, we started investigating whether such a dependence
exists (DLG95). As a measure of the energy input into the disk ISM
per unit time,
,
we used the nonthermal radio continuum
emission produced by type II SNe. In order to obtain a good
measure of the energy density in a given volume within the
galaxy disk, the area over which energy input by active SF
takes place,
,
was quantified by us. Determining
from radial profiles of the radio continuum emission distribution
in the disk a cut-off radius of SF processes,
,
we calculated, adopting circular symmetry,
,
as
the circular area within the disk of the galaxies over which
SF-related energy production occurs, following the simple
relation
.
can be determined equally well from the radial distribution
of radio continuum emission in the disk.
Data of in total 9 galaxies indicated a trend that those galaxies
with the highest energy injection rates,
,
per unit time and unit surface area, have the most prominent
radio halos (DLG95). Below a certain energy input level no outflows
are enabled and galaxies thus do not exhibit extraplanar radio
emission.
Similar results were reached by Rand (1996) and Meurer et al.
(1995, 1997), based on data from other wavebands. A possible
link between high energy input rates into the disk ISM and the
existence of soft X-ray halo emission is discussed by Ehle
et al. (1998) and Dahlem et al. (1998).
The remaining 11 galaxies in the present sample can also be
studied for such a dependence. The relevant properties are
listed in Table 7.
We use as a measure of
the beam-deconvolved
radial extent of the 1.43 GHz radio continuum emission from
our current observations. We limit our studies here to the
global energy input as a measure of which we will use the
total FIR luminosity
,
as done by Rossa &
Dettmar (2000).
Based on the total FIR flux, FIR, following the relation
FIR =
), we calculate
.
is the FIR luminosity normalised per unit surface area of the
star-forming part of the disk, representing the energy input
rate (
), as defined similarly by us based
on radio surface brightnesses (DLG95). This quantity is used
as a measure of the global rate of energy injection into the
ISM via SF-related processes (stellar winds and type II
supernovae), see Col. 7 in Table 7.
,
finally, is the exponential scale height of
the radio halos from Table 4.
Rossa & Dettmar (2000) have made further progress in another
direction, plotting the
f60/f100 flux ratio vs. the
total FIR luminosity (as a measure of
)
normalised by
the disk surface area,
(in particular,
their Fig. 13), where they find a clear distinction between
quiescent and starburst galaxies.
All galaxies classified as starbursts by their location in the
f60/f100 vs.
plot (their Fig. 13)
have detections of H
emission from extraplanar diffuse
ionised gas (eDIG). Compared to DLG95, the FIR luminosity,
,
substitutes the radio luminosity to represent
the level of energy injection,
.
As a first approximation, Rossa & Dettmar use in their
determination of the surface area of the energy input the
25th magnitude isophotal diameter, D25. For the
galaxies with wide-spread SF over their disks as used in
their sample, this approximation is acceptable.
However, in other cases-especially the classic starburst
galaxies with circumnuclear activity-the use of D25 can
lead to serious errors: systems like, e.g., NGC1808 actively
form stars over a much smaller area than their disks out to
D25 (e.g., Dahlem et al. 1990). In the present sample
NGC3175 is such a case. While its D25 is 5', the
diameter of the star-forming part of its disk can be determined
(from H
and radio continuum images; see Condon et al.
1996; Ryder & Dopita 1994) very reliably to be only 90''
(7 kpc). Thus, using D25 would lead to an over-estimate of
the surface area over which its energy input is distributed, by
a factor of 11.1 and thus underestimate the energy density in
the disk by the same factor. In these cases one must determine
the area over which SF takes place more accurately, by using
and
instead.
Figure 38 corresponds with Figs. 11 and 12 by Rossa
& Dettmar (2000). It displays the
f60/f100 FIR flux
ratio (thus a measure of the mean dust temperature) vs. the
total FIR luminosity, normalised with the area in which active
SF is taking place. In this logarithmic plot all galaxies occupy
the upper right-hand quadrant where, according to Rossa & Dettmar
(2000) the starburst galaxies fall. Based on our selection criteria
the bottom left-hand corner of the plot is empty; there are no
quiescent galaxies in the sample.
Amongst the galaxies plotted in Fig. 38 a trend is
visible that the galaxies with the highest
values have the highest
f60/f100 flux ratios and thus
mean dust temperatures.
![]() |
Figure 38:
Plot of the
f60/f100 FIR flux density ratios
vs. the FIR luminosity normalised by the area in which star
formation and thus energy input into the ISM takes place,
|
| Open with DEXTER | |
![]() |
Figure 39:
Plot of the
f60/f100 FIR flux density ratios
vs. the 1.43 GHz radio luminosity normalised by the area in
which star formation and thus energy input into the ISM takes
place,
|
| Open with DEXTER | |
| Galaxy |
|
|
|
FIR |
|
|
|
| (''/kpc) | (103 kpc2) | (10-14 W m-2) | (1043 erg s-1) | (1040 erg s-1 kpc-2) | (kpc) | ||
| NGC1055 | 0.34 | 236/18.3 | 1.06 | 135.85 | 4.16 | 3.94 | -- |
| NGC1406 | 0.44 | 110/ 7.9 | 0.20 | 72.74 | 1.93 | 9.75 | 1.66 |
| NGC1421 | 0.40 | 145/21.9 | 1.50 | 54.43 | 6.30 | 4.20 | -- |
| NGC1511 | 0.62 | 73/ 6.2 | 0.12 | 124.81 | 4.57 | 38.08 | 1.76 |
| NGC2748 | 0.39 | 103/14.3 | 0.65 | 45.54 | 4.49 | 6.96 | -- |
| NGC3175 | 0.46 | 45/ 3.5 | 0.04 | 78.02 | 2.36 | 62.11 | 1.62 |
| NGC3437 | 0.57 | 96/11.9 | 0.48 | 62.77 | 4.88 | 10.92 | -- |
| NGC3717 | 0.44 | 177/23.3 | 1.71 | 60.24 | 5.29 | 3.09 | -- |
| NGC4700 | 0.57 | 54/ 6.6 | 0.14 | 16.67 | 1.30 | 9.35 | 3.06 |
| NGC7090 | 0.33 | 131/ 7.4 | 0.17 | 41.76 | 0.68 | 3.93 | 1.77 |
| NGC7462 | 0.49 | 103/ 7.5 | 0.18 | 18.07 | 4.93 | 27.54 | 1.82 |
| NGC1569 | 0.96 | 72/0.77 | 0.00186 | 207.20 | 0.12 | 64.50 | -- |
The new result of our investigation is that in 6 out of 11 objects
radio halo emission has been detected (open squares). There is a
clear trend that halos were detected in those galaxies with the
highest
and
f60/f100 values, with the
notable exception of NGC3437, which was discussed above.
We have selected highly inclined galaxies only for technical
reasons. The fact that most, if not all FIR-warm edge-on
galaxies have radio halos should apply to galaxies with any
inclination angle. We just need to develop tools to find
them in more face-on systems.
We have produced the same kind of plot as Fig. 38 based on the total 1.43 GHz radio continuum flux density, as
listed in Table 3, which is presented in
Fig. 39. The total 1.43 GHz radio luminosity
was calculated using the equation
,
neglecting the term
,
because it
is not relevant for the low-redshift galaxies studied here
(Condon et al. 1990).
Because of the strong radio-FIR correlation the two plots look
almost identical, although the positions of individual galaxies
can shift with respect to their location in Fig. 38.
This similarity of the two plots-and thus relations between
luminosity and dust temperature-makes our present study and
the results by Rossa & Dettmar (2000) comparable to our
earlier investigations (DLG95). As in Fig. 38,
the data points in Fig. 39 also indicate a
trend of the galaxies with the highest
values having the highest mean dust temperatures.
A second important result is that all physically small galaxies
(with
kpc; in the present sample NGC1511,
NGC4700, NGC7090 and NGC7462) have prominent radio halos.
Their total mass being low, it appears to be easier for CR
electrons to escape from the thin disk than in galaxies with
higher total masses,
.
Even NGC7090, with a moderate
mean dust temperature and also a correspondingly moderate energy
injection rate (
;
note that NGC1055 and
NGC7090 have almost equal values, see Table 7),
can produce a radio halo. Radio halos are visible, because these
galaxies are no dwarfs (in which case they might be too light to
retain their magnetic fields so that CRs might escape almost
loss-free and thereby remain invisible), but intermediate-mass
systems.
We note here that the simplifying assumption has been made that all energy contained in CR electrons emitting synchrotron radiation originally arises from the galaxy disks. Lacking images showing details, such an assumption must be made, without any distinction as to where the radio continuum emission comes from. The same assumption is made for the FIR data too (not only by us, but also by Rossa & Dettmar 2000 and other investigators).
Based on a larger sample, Irwin et al. (1999) did not find
a connection between the properties of gaseous halos and the
global level of underlying SF activity, using radio continuum
images. However, these authors did not distinguish
between different object classes contained in their sample due to the small number of galaxies in total.
From their sample, NGC3735 and NGC4388 should be excluded,
because they host Sy-2 nuclei (NED). The existence of a
radio halo in NGC3432, despite its low energy injection rate,
can be explained because it is a magellanic irregular with low
total mass and thus shallow gravitational potential, similar to
the small galaxies in our sample. In addition, NGC3432 is
closely interacting and should therefore be removed from the
sample for studies of the dependence of halo properties on the
level of underlying SF activity in the galaxy disk. NGC5433
and IC562 from the sample by Irwin et al. (1999) should be
excluded, as we did above with NGC5073 and NGC7541, because
their emission distribution is not resolved by the data. Once
this target selection is applied, we do see a trend in the data
by Irwin et al. (1999) suggesting that the galaxies with the
highest
values have halos, while those with lower
energy input levels do not always have halos. Therefore, there
is no contradiction in the data by these authors to our hypothesis
from the paper by DLG95. Instead, their data support our claim,
as do our new results presented here, that late-type spiral
galaxies with warm dust do have gaseous halos and that those
with the most active SF tend to have the most prominent halos.
However, trying to go one step further, we find that there is
clearly no direct quantitative relationship between the
halo scale height,
,
and the normalised energy
input,
,
as listed in Table 7.
In our earlier sample (DLG95) we presented only L* galaxies,
i.e. systems near the "knee'' of the galaxy luminosity function,
with total masses somewhere on the order of
.
There, we found a trend that galaxies with the highest SF rates
in their disks have the most prominent halos and that the halos
are brightest above the most actively star-forming regions.
The present sample, as indicated above, contains a number of physically small galaxies as well as L* galaxies. This implies that, compared to our earlier study (DLG95), an additional free parameter has entered the equation, namely the different resistence that particles meet when expelled from their birth sites in the disks.
We have therefore tried to normalise the energy input rate
not only by the area over which it is distributed, but also
by the gravitational potential of the galaxy as measured by
the total mass. Also for this quantity we have found no
convincing correlation with
.
This implies
that yet more, or other, factors affect the propagation of
CRs out of galaxy disks.
In particular, the energy losses that CR electrons undergo
on their way out of the disk influence the properties of
the observed radio halos:
Although predicting different kinds of behaviour, both
static and dynamic models of CR radio halos (Lerche &
Schlickeiser 1981a-c; Werner 1988) predict a dependence
of the extent of radio halos on the mean lifetime of
the CR electrons. In case of pure diffusive CR propagation
the dependence is directly proportional to the diffusion
index of the CRs. The influence of convection is expected
to increase with increasing energy input. However, at the
same time inverse Compton losses become more important
with increasing energy input too. This is the reason why
no direct dependence is found between the measured scale
heights and the level of energy input in the underlying
disks.
Significant energy losses of the CR electrons leaving the disk are in agreement with the observed slight steepening of the radio spectral indices (Sect. 3.3).
Based on images with sufficient spatial resolution a direct
morphological comparison can be made, as done by us earlier
(DLG95), of the radial extent of the radio halo with respect
to the radial extent of the actively star-forming part of
the disk,
.
We had found in our previous sample
that the radio halo emission breaks off radially near the
point where the surface brightness in the disk drops rapidly,
i.e.
.
This implies that halo emission is observed only above
the most actively star-forming parts of the galaxy disks. In
the present sample, the ATCA data have the highest angular
resolution. The effect is visible most clearly in NGC7090,
and also-though to a lesser degree-in NGC7462 (here
especially at 2.45 GHz, where the angular resolution is highest).
Higher resolution is required to study the same effect in the
galaxies observed with the VLA D array, especially in cases
like NGC4700 (cf. Fig. 11). In the case of
NGC3175 the data by Condon et al. (1996) suggest the same
behaviour. The extraplanar radio emission is restricted to the
radial regime of the central starburst. More details on this
will be presented elsewhere.
Despite the relatively good spatial resolution of our data, NGC1511, with its very smooth emission distribution, shows little of this effect. It comes closest to the ellipsoidal emission distribution that might be expected in case of purely diffusive particle propagation.
The above results prove that, after exclusion of host galaxies of AGNs, their f60/f100 FIR colour (and thus mean dust temperature) is a very successful criterion to select candidate galaxies for searches of radio halos (gaseous halos in general). Sensitive radio continuum imagery is a powerful tool to detect the halo emission.
For the 11 suitable objects which are presented here, we have
a detection rate of 55%. Thus, again a substantial fraction
of FIR-warm objects, although fainter than those studied by us
earlier (DLG95; D98), exhibit halo emission. In some of the
remaining cases, especially the galaxies with high
f60/f100
ratios but no detected radio halos, the angular resolution of
our data might not be sufficient to separate thin disk and halo
emission.
All physically small galaxies in our sample, even some with
,
have radio halos.
The measured z scale heights of the radio halos range from
about 1.4 to 3.1 kpc.
Our results suggest that the galaxies with the highest energy input rates into their disk ISM are the ones that have the most prominent radio halos. However, there is no direct relationship between the halo scale heights and the level of energy input in the underlying disks, because the halo properties far away from the disk planes are dominated by the energy losses of the CR electrons.
To confirm the results obtained here and to improve the angular resolution of the datasets where it limits our investigation, we plan to extend our VLA observations at 1.43 GHz by adding longer baselines from the C array. This will enable us to study the dependence of the intrinsic radio halo properties on the local level of SF in the underlying disks, as already done for NGC891 and NGC4631 (DLG95), in a larger sample.
Acknowledgements
We thank both the ATNF and NRAO for the generous allocation of observing time and the referee, Dr. J. Condon, for fruitful discussions, which helped improve the paper considerably. We thank J. Rossa very much for communicating to us his results of Hobservations prior to publication. M. D. thanks Drs. R. Sault, R. Wark and H. May at the ATNF for their software and computing support and Dr. F. Israel for his kind hospitality at Sterrewacht Leiden, where most of this work was done. Thanks are also due to the VLA analysts for their support during the observations and data reduction. Many thanks to Dr. M. Gotzens-Petr, J. Krist and K. Kraiberg-Knudsen for helping us with the production of the z profiles and to Dr. G. Meurer for the calculation of the distance values. We thank Dr. J. Irwin for making available to us the VLA C array observations of NGC1421 and Dr. J. Condon for the C array data of NGC3175. Part of J. L.'s work was performed as part of an ATNF Vacation Scholarship. This research has made use of the NASA Extragalactic Database (NED), whose contributions to this paper are gratefully acknowledged. The Digitized Sky Survey was produced at the Space Telescope Science Institute under U.S. Government grant NAG W-2166. The National Geographic Society - Palomar Observatory Sky Atlas (POSS-I) was made by the California Institute of Technology with grants from the National Geographic Society.