A publishing partnership

ArticlesThe following article is Free article

CHANG-ES. III. UGC 10288—AN EDGE-ON GALAXY WITH A BACKGROUND DOUBLE-LOBED RADIO SOURCE

, , , , , , , , , and

Published 2013 November 15 © 2013. The American Astronomical Society. All rights reserved.
, , Citation Judith Irwin et al 2013 AJ 146 164 DOI 10.1088/0004-6256/146/6/164

1538-3881/146/6/164

ABSTRACT

This third paper in the Continuum Halos in Nearby Galaxies—an EVLA Survey (CHANG-ES) series shows the first results from our regular data taken with the Karl G. Jansky Very Large Array. The edge-on galaxy, UGC 10288, has been observed in the B, C, and D configurations at L band (1.5 GHz) and in the C and D configurations at C band (6 GHz) in all polarization products. We show the first spatially resolved images of this galaxy in these bands, the first polarization images, and the first composed image at an intermediate frequency (4.1 GHz) which has been formed from a combination of all data sets. A surprising new result is the presence of a strong, polarized, double-lobed extragalactic radio source (CHANG-ES A) almost immediately behind the galaxy and perpendicular to its disk. The core of CHANG-ES A has an optical counterpart (SDSS J161423.28−001211.8) at a photometric redshift of zphot = 0.39; the southern radio lobe is behind the disk of UGC 10288 and the northern lobe is behind the halo region. This background "probe" has allowed us to do a preliminary Faraday rotation analysis of the foreground galaxy, putting limits on the regular magnetic field and electron density in the halo of UGC 10288 in regions in which there is no direct detection of a radio continuum halo. We have revised the flux densities of the two sources individually as well as the star formation rate (SFR) for UGC 10288. The SFR is low (0.4–0.5 M yr−1) and the galaxy has a high thermal fraction (44% at 6 GHz), as estimated using both the thermal and non-thermal SFR calibrations of Murphy et al. UGC 10288 would have fallen well below the CHANG-ES flux density cutoff, had it been considered without the brighter contribution of the background source. UGC 10288 shows discrete high-latitude radio continuum features, but it does not have a global radio continuum halo (exponential scale heights are typically ≈1 kpc averaged over regions with and without extensions). One prominent feature appears to form a large arc to the north of the galaxy on its east side, extending to 3.5 kpc above the plane. The total minimum magnetic field strength at a sample position in the arc is ∼10 μG. Thus, this galaxy still appears to be able to form substantial high latitude, localized features in spite of its relatively low SFR.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

This is the third paper in the series, "Continuum Halos in Nearby Galaxies—an EVLA Survey" (CHANG-ES), whose goals are to map a sample of 35 edge-on galaxies at L band (1.5 GHz) and C band (6.0 GHz) in all four Stokes parameters, in a variety of configurations of the Expanded Very Large Array (EVLA). The EVLA is now known as the Karl G. Jansky Very Large Array (JVLA) of the National Radio Astronomy Observatory10 (NRAO) and we will refer to it as the JVLA from now on. A complete description of the galaxy sample and goals of the survey can be found in Irwin et al. (2012a, hereafter Paper 1) and first results based on early test data of NGC 4631 can be found in Irwin et al. (2012b, hereafter Paper 2).

UGC 10288 (Table 1) is a member of the Lyon Group of Galaxies 404 (Garcia et al. 1993) which contains only two other, widely separated members: UGC 10290, 1fdg0 (0.61 Mpc; assuming the distance of UGC 10288) to the north, and NGC 6070, 1fdg4 (0.86 Mpc) to the northwest. UGC 10288, itself, is edge-on (inclination of 90°) and is listed in the Revised Flat Galaxy Catalog (Karachentsev et al. 1999).

Table 1. Properties of UGC 10288a

Property Value
R.A. (J2000) 16h14m24fs80
Decl. (J2000) −00°12m27fs1
Inclination (deg) 90
d25 (')b 4.9
V (km s−1)c 2046
D (Mpc)d 34.1
Morphological type Sc
S1.4 GHz (mJy)e 26.1
LFIR (L)f 2.55 × 109
SFRIRAS (M yr−1)g 1.3
SFRthis work (M yr−1)h 0.4–0.5
ρ (Mpc−3)i 0.23
W20 (km s−1)j 397.3
MT (M)k 2.2 × 1011
Sint (Jy km s−1)l 36.4
$M_{{\rm H\,\scriptsize{I}}}$ (M)m 9.9 × 109

Notes. aSee Paper 1 for original sources if not specified. bObserved blue diameter at the 25th mag arcsec−2 isophote. cHeliocentric systemic velocity. dDistance, assuming H0 = 73 km s−1 Mpc−1 and correcting for Virgo Cluster and Great Attractor perturbations. ePreviously measured flux density at a frequency of 1.4 GHz. fFar-IR luminosity. gStar formation rate based on IRAS photometry. hSee Section 3.1.2 for an explanation. iDensity of galaxies brighter than −16 mag in the vicinity of UGC 10288. jWidth of the H i profile at the 20% level from the H i Parkes All Sky Survey (HIPASS) Catalogue (Meyer et al. 2004). kTotal mass, using $M_T = {\rm W}_{20}^2\,{\rm d}_{25}/(8\,G)$ (cgs units). Note that this result is a factor of two higher than the quantity listed in Paper I which was taken from older H i data. lHIPASS integrated flux density of the H i line (Meyer et al. 2004). This value may be compared with Sint = 38.2Jy km s−1 and Sint = 34.7Jy km s−1 using the NRAO 43 m and Green Bank telescopes, respectively (Hogg et al. 2007). mH i mass, from $[{M}_{{\rm H\,\scriptsize{I}}}/{M_\odot }] = 2.35\,\times \,10^5 [D/{\rm Mpc}]^2\, [S_{{\rm int}}/({\rm Jy \, km\, s^{-1})}]$.

Download table as:  ASCIITypeset image

UGC 10288 has not previously been observed extensively, but it has been a target in a few searches for high-latitude emission. Rand (1996) detected a few Hα spurs extending to 1 kpc below the galaxy's disk at an emission measure level of 40 pc cm−6. Follow-up single-slit optical spectra at two positions (Collins & Rand 2001) then revealed emission lines to a vertical z height11 of z = 3 kpc and showed that the line ratios were best described by composite photoionization/shock models. An Hα observation by Rossa & Dettmar (2000) was too low in signal-to-noise (S/N) to detect more than the brightest disk emission. Alton et al. (2000) claim a dust absorption exponential scale height that is unresolved at 1farcs2 resolution (198 pc). Of these references, only Rand (1996) has published an original image of the galaxy.

Aside from the 1.4 GHz NRAO VLA Sky Survey (NVSS), whose rms noise is 450 μJy beam−1 in a 45'' beam (Condon et al. 1998), the only previous radio continuum observation of UGC 10288 was made by Hummel et al. (1991) at 5.0 GHz to an rms limit of 90 μJy beam−1 with a beam of 14farcs5. The latter image shows essentially a point source with a slight elongation to the north which has been interpreted as evidence for a radio halo (Rossa & Dettmar 2003). Because of these new JVLA observations, however, we now know that this elongation is due to a strong background double-lobed extragalactic radio source (EGRS). Our CHANG-ES observations achieve a best rms of 3 μJy beam−1 for a single array/frequency band combination (C array, C band), and provide a range of spatial resolution from 3farcs0 to 37'', significantly improving on these previous data.

2. OBSERVATIONS AND DATA REDUCTIONS

Observations and data reductions of UGC 10288 were similar to those described for the test data in Paper 2. However, there were also some differences12 and, since this is the first paper describing our "regular" (as opposed to "test") observations, we describe each step fully below. Note that these data reductions are consistent with current standard continuum tutorials,13 unless otherwise indicated.

2.1. Observations of UGC 10288

Observations were carried out using the JVLA during its commissioning period in L band (B, C, and D arrays) and at C band (C and D arrays). A log of the observational setup for each JVLA configuration and each frequency band is given in Table 2. A single array and frequency combination will be referred to as an "observation" in the following. As can be seen, the same calibrators were used for each observation, and the same correlator setup was used. We introduced a gap in the frequency band coverage at L band (from 1.503 to 1.647 GHz) because of persistent strong interference in that frequency range.

Table 2. UGC 10288 Observing and Calibration Informationa

Frequency 1.5 GHz (L band) 6.0 GHz (C band)
Array B C D C D
Date of observations 2011 Apr 5 2012 Mar 30 2011 Dec 30 2012 Feb 11 2011 Dec 10
        2012 Feb 16  
Frequency range (GHz)b 1.247 → 1.503 1.247 → 1.503 1.247 → 1.503 4.979 → 7.021 4.979 → 7.021
  1.647 → 1.903 1.647 → 1.903 1.647 → 1.903    
Total bandwidth (MHz) 512 512 512 2042 2042
No. of spectral windows 32 32 32 16 16
No. of channels per spectral window 64 64 64 64 64
Total no. channels 2048 2048 2048 1024 1024
Channel separation (MHz) 0.25 0.25 0.25 2.0 2.0
Spectral resolution (MHz)c 0.50 0.50 0.50 4.0 4.0
Integration time (s)d 10 10 10 10 10
Obs. time (minute)e 96.7 43.0 18.3 191 38.7
Primary calibratorf 3C 286 3C 286 3C 286 3C 286 3C 286
Secondary calibratorg J1557−0001 J1557−0001 J1557−0001 J1557−0001 J1557−0001
$S_{\nu _0}$(Jy)h 0.487 ± 0.004 0.461 ± 0.002 0.47 ± 0.01 0.429 ± 0.001 0.432 ± 0.002
        0.422 ± 0.004  
Pol.-leakage calibratori J1407+2827 J1407+2827 J1407+2827 J1407+2827 J1407+2827
$S_{\nu _0}$(Jy)h 0.953 ± 0.004 0.963 ± 0.003 0.943 ± 0.003 2.094 ± 0.002 2.103 ± 0.002

Notes. aThese values apply to the setup of the observations and do not take into account the flagging of bad data. bThe frequency range of 1.503–1.647 GHz was avoided due to interference. At C band, the upper and lower portions of the band overlap slightly at the center; consequently, the total bandwidth is slightly less than 1024 channels × 2 MHz channel−1 = 2048 MHz. cAfter Hanning smoothing. dSingle record measurement time. eOn-source observing time before flagging. fThis source was also used as the bandpass calibrator and for determining the absolute position angle for polarization. gThis source is a "primary" calibrator in the sense of its amplitude errors (<3% amplitude closure errors expected) in all arrays and in both bands. It is separated from UGC 10288 by 4fdg14 in the sky. hFlux density of the calibrator specified on the previous line, at a representative frequency, ν0, within the frequency band. For L band, we have taken ν0 = 1.495 GHz and for C band, ν0 = 5.94 GHz. Uncertainties are based on the extrapolation of amplitudes and phases from 3C 286 after phase, amplitude, and bandpass calibrations have been applied. Uncertainties at other frequencies within the respective bands are similar. For the C-array C-band data, the two rows correspond to the two dates of these observations, in the same order as the listed dates. iThis zero-polarization calibrator, also known as OQ208 and QSO B1404+2841, was used to determine the polarization leakage terms.

Download table as:  ASCIITypeset image

The observations were carried out in a standard fashion with scans of the source flanked by scans on a complex gain calibrator (J1557−0001) that was near it in the sky. We will refer to this source as the "secondary calibrator." In addition, a scan was made on a bright calibrator of known flux density and brightness distribution (3C 286), i.e., the "primary calibrator" (this source is also used for the bandpass and polarization angle calibration). Finally, a scan was also carried out on a calibrator that was known to have negligible polarization (J1407+2827 = OQ 208), the latter to determine the polarization leakage (see Section 2.4). The longest consecutive scan on the galaxy at any array, between scans of the secondary calibrator, was 25 minutes. Scans of 3C 286, J1557−0001, and J1407+2827 were typically 7.5–11 minutes, 4.5 minutes, and 8–11 minutes, respectively. For each observation, the galaxy was observed within a scheduling block (SB) that included other galaxies. This allowed us to observe any given galaxy near the beginning of the SB and then again near the end, so that the broadest possible uv14 coverage could be achieved. For UGC 10288, the gaps in time varied with array and frequency, ranging from about 1.5–4.8 hr.

2.2. Total Intensity Calibration

Data were reduced using the Common Astronomy Software Applications (CASA) package,15 versions 3.3, 3.4, and 4.0, and each observation was reduced in an identical fashion, as described below, unless otherwise indicated. Note that JVLA antennas detect right (R) and left (L) circularly polarized radiation. In the following description, R and L are treated separately in any antenna-based calibration and the parallel hands (RR and LL) are treated separately in any baseline-based calibration. Calibration of the cross-hands (RL and LR) is discussed in Section 2.4.

A copy of the data set was first made and then Hanning-smoothed so as to more easily detect bad data. An initial pass of flagging then took place to remove the worst of the radio frequency interference (RFI).16 An antenna position correction table was then produced, if needed, based on known positions as posted by NRAO.17 An antenna-based delay correction table was then formed using 1 minute of the scan on the primary calibrator, during which the antenna position correction table was applied on the fly. Delay errors produce phase slopes as a function of frequency; we found delay corrections of 6 ns or less, as determined with respect to a reference antenna in the same polarization (R or L). The same RFI flags determined in the first pass of the Hanning-smoothed data as well as the antenna position correction and delay correction tables were then applied to the original non-Hanning-smoothed data. This ordering of steps (which is not described in standard CASA tutorials) was needed because delay corrections must be applied to non-Hanning-smoothed data,18 whereas locating RFI is most effectively done on Hanning-smoothed data. The corrected data were then split off into a separate data set and Hanning-smoothed for further processing.

Another round of RFI flagging was then carried out until there were no obvious amplitude departures from typical values. In addition, each observing band contains contiguous spectral windows (Table 2), each with its own bandpass response. Since the bandpass response in each spectral window declines strongly at the upper and lower edges, we flagged five channels at both the upper and lower ends of each spectral window. JVLA observing parameters did not allow us to overlap the spectral windows without losing significant total bandwidth19 so that, after flagging, there are gaps of 10 channels between each spectral window. The flux density of the primary calibrator was then set for each channel in the band using a known model20 for Stokes I.

An initial calibration table of phase as a function of time for the primary calibrator was made for a small range of channels at the center of each spectral window (where the bandpass response is reasonably flat). A bandpass calibration table was then made using the same calibrator, applying the previous phase calibration table on the fly. Application of the initial phase calibration table ensures that no decorrelation occurs when vector averaging the data to determine bandpass solutions. Bandpasses were determined for each antenna and each spectral window. The initial phase calibration table is then no longer used.

Phase and amplitude calibration tables were then determined separately (first phase and then amplitude) for the primary calibrator, the secondary calibrator, and the polarization leakage calibrator. The bandpass calibration table was applied on the fly during these steps and the phase calibration table was also applied on the fly when the amplitude calibration table was formed. Phase solutions were determined for each integration time for the primary calibrator, and for each scan for the secondary and polarization leakage calibrators. Amplitude solutions were determined per scan for all calibrators. Flux densities for the secondary and polarization leakage calibrators were then found from the known primary calibrator values on a per spectral window basis. Results for a representative frequency near the center of each band are given in Table 2. The results at L band over the various independent observations agree to 1% for J1407+2827 and 3% for J1557−0001. At C band, the agreement is to within 0.2% for J1407+2827 and 1% for J1557−0001.

All calibration tables were then applied to the data (calibrators plus galaxy) to form a calibrated data set in total intensity, I. For the galaxy, the calibration solutions determined from the secondary calibrator were applied, interpolated linearly with time. At this point, low level RFI, which was not evident in the uncalibrated data, became much more obvious. Therefore, another round of flagging ensued and the above steps were all repeated, with new tables made and a new calibration applied. The flux density data for the calibrators presented in Table 2 refer to the final results.

2.3. Total Intensity Imaging and Self-calibration

Total intensity images were formed for each observation by first Fourier transforming the uv data into the sky plane (forming the "dirty map") and then deconvolving the "dirty beam" (i.e., the point-spread function) from the dirty map using the Clark CLEAN algorithm (Clark 1980). Further details regarding our imaging of wide-band data can be found in Paper 2 and the parameters used to make our maps are in Table 3.

Table 3. UGC 10288 Map Parameters (Individual Data Sets)

Frequency 1.5 GHz (L band) 6.0 GHz (C band)
Array B C D C D
uv weightinga Briggs Briggs Briggs Briggs Briggs
No. Self-calsb none 1 a&p 1 a&p 2 a&p + 1 p 1 a&p
  I images
Clean spatial scales ('')c 0, 7.5, 15 0, 20, 40, 80 0 0, 5, 10, 20 0, 10, 20, 40, 60
Synth. beamd          
('', '', °) 3.80, 3.58, 66.2 12.18, 10.71, −55.0 40.23, 34.34, −31.3 3.14, 2.87, −11.0 10.96, 9.46, −28.5
Map peak (mJy beam−1)e 9.57 12.6 19.2 7.62 8.94
rms (μJy beam−1)f 14 22 42 3.1 7.0
  Q and U imagesg
Synth. beamh          
('', '', °) 3.72, 3.54, 60.8 12.18, 10.71, −55.0 40.19, 34.33, −31.4 3.14, 2.87, −11.0 10.95, 9.47, −28.4
rms (μJy beam−1)i 13 22 35 2.9 5.4
  P imagesj
Map peak (μJy beam−1)e 177.9 178.8 272.1 149.5 146.7

Notes. These parameters represent the images presented in Figures 1-4. aSee Briggs (1995) for a description of Briggs weighting with various "robust" factors. Robust = 0 was used in each case, as employed in the CASA clean task. bNumber of self-calibration iterations, where "p" refers to phase and "a&p" refers to amplitude and phase together. cScales used for the multi-scale clean. The scale, 0, corresponds to a classic clean for which the emission is assumed to be described by point sources. dSynthesized beam FWHM of major and minor axis, and position angle. eMaximum value on the map. fRms map noise before primary beam correction. gStokes Q and U images. The clean spatial scales used were equivalent to the I images with the highest spatial scale dropped (except for D-array, L-band images for which only the zero spatial scale was used). hThe cross-hands (RL, LR) had slightly different sets of flags applied to them compared to the parallel hands (RR, LL), leading to small differences in the synthesized beams for Q and U compared to I. iThe rms noise is lower than for I because the lower and less extensive polarized emission results in fewer clean errors. jPolarized intensity images, where P = $\sqrt{Q^2\,+\,U^2}$, corrected for bias.

Download table as:  ASCIITypeset image

In brief, we used the wide-field imaging algorithm (Cornwell et al. 2008) and multi-scale, multi-frequency synthesis (ms-mfs) during imaging (Cornwell 2008; Rau & Cornwell 2011). Wide-field imaging takes into account the fact that all antennas may not be coplanar. The multi-scale clean assumes that the emission can be represented by a variety of spatial scales, rather than only point sources as has been largely done in the past. Multi-frequency synthesis makes use of the fact that, at different frequencies within the band, the uv spacings (measured in kλ) are also different and therefore the uv plane is more filled in than would be the case for monochromatic observations. We also fit the emission across the band in each pixel with a spectrum (Sault & Wieringa 1994) of the form Iν ∝ να, which, in practice, is expanded as a Taylor series with two terms (see Paper 2). At the end of the cleaning process, the image is constructed by reconvolving the clean components with a Gaussian of the same FWHM as the dirty beam.

Each data set was then self-calibrated (e.g., Pearson & Readhead 1984), with the exception of D-array L band for which self-calibration attempts did not improve the map. The number and type of self-calibration iterations is given in Table 3.

The resulting total intensity images are seen in Figures 1 and 2 for the L-band and C-band data, respectively. We additionally formed images (not shown) at spatial resolutions intermediate between those shown in these figures (e.g., by tapering the uv plane response). The images of Figures 1 and 2 are shown without a primary beam correction, but any flux density (or other numerical) measurements are made on primary beam-corrected images. The primary beam correction was carried out as described in Paper 2, except that the beam has now been constructed with frequency weighting that matches the frequency weighting of the data (for example, the weighting that necessarily results from flagging).21

Figure 1.

Figure 1. L-band images of UGC 10288 (both gray scale and contours). The beam size is shown as a black ellipse at the bottom left and corresponding map parameters are given in Table 3. (a) D-array map; contours are at 84 (2σ), 120, 200, 350, 600, 1200, 3000, 6000, 12,000, and 18,000 μJy beam−1. (b) C-array map; contours are at 44 (2σ), 60, 90, 150, 200, 400, 800, 1600, 3000, 6000, and 12,000 μJy beam−1. The background sources are labeled. (c) B-array map; contours are at 28 (2σ), 50, 100, 200, 400, 1200, 3600, and 8000 μJy beam−1.

Standard image High-resolution image
Figure 2.

Figure 2. C-band images of UGC 10288 (both gray scale and contours). The beam size is shown as a black ellipse at the bottom left and corresponding map parameters are given in Table 3. (a) D-array map; contours are at 14 (2σ), 20, 40, 80, 150, 300, 600, 1200, 3000, and 6000 μJy beam−1. Dashed lines denote the region used to estimate the fluxes of Table 5. (b) C-array map; contours are at 6.2 (2σ), 12, 20, 30, 60, 120, 200, 400, 1600, and 6000 μJy beam−1.

Standard image High-resolution image

2.4. Polarization Calibration and Imaging

Polarization calibration proceeded in the same fashion for all data sets. The cross-correlation data (RL and LR) were first inspected separately and additional flagging was carried out as needed. A known model for Stokes Q and U for the primary calibrator was introduced, based on formulae given in Paper 2. The relative antenna-based delay correction between R and L was then derived using the primary calibrator. The leakage terms between the R and L circularly polarized feeds were found from the polarization leakage calibrator (Table 2). Finally, the absolute polarization position angle on the sky was determined from the known position angle for the primary calibrator (χ = 33° for 3C 286). During each of these steps, all previously determined calibration tables were applied on the fly. Finally, all calibration tables, including the self-calibration table determined from the corresponding total intensity map, were applied to the RL and LR data.

Imaging of Stokes Q and U proceeded in the same fashion as for I but the largest spatial scale was dropped during cleaning (except for D-array L band) since the extent of the polarized emission is not as great as in total intensity. Maps of polarized intensity were then formed for each array and frequency, correcting for the bias introduced by the fact that P images do not obey Gaussian statistics (Simmons & Stewart 1985; Vaillancourt 2006). Maps of E vectors rotated by 90° were also formed after first blanking all emission less than 3σ, where σ is the rms noise of the respective Q and U maps, resulting in an uncertainty of ±10°. These vectors would be equivalent to the intrinsic direction of the regular magnetic field, B, if there were no Faraday rotation. The maps of polarized intensity (without primary beam corrections) superimposed with the rotated E vectors are shown in Figures 3 and 4.

Figure 3.

Figure 3. L-band polarization images of UGC 10288 (green contours) with E vectors rotated by 90° (red, formed where Q and U > 3σ). The gray scale represents the corresponding total intensity (I) image along with its 2σ contour (black) from Figure 1. Map parameters are given in Table 3, except that the D-array polarization image has been made from Q and U data that have been smoothed slightly to match the corresponding I data. (a) D-array map; contours are at 105 (3σ), 150, and 220 μJy beam−1. (b) C-array map; contours are at 65.1 (3σ), 110, and 175 μJy beam−1. (c) B-array map; contours are at 38.1 (3σ), 90, and 175 μJy beam−1.

Standard image High-resolution image
Figure 4.

Figure 4. C-band polarization images of UGC 10288 (green contours) with E vectors rotated by 90° (red, formed where Q and U > 3σ). The gray scale represents the corresponding total intensity (I) image along with its 2σ contour (black) from Figure 2. Map parameters are given in Table 3. (a) D-array map; contours are at 16.2 (3σ), 35, 70, and 120 μJy beam−1. (b) C-array map; contours are at 8.7 (3σ), 30, 75, and 140 μJy beam−1.

Standard image High-resolution image

2.5. Combined Array/Frequency Images and Spectral Index Map

The uv data at each frequency (B, C, and D arrays at L band and C and D arrays at C band) were combined in order to obtain higher S/N and more complete uv coverage for each of the bands. The results of these combinations are parameterized in Table 4 and each map was also corrected for the primary beam as described above. We will refer to them as needed.

Table 4. UGC 10288 Map Parameters (Combined Data Sets)

Parameter BCD L band CD C band BCD+CL (rob0) BCD+CL (uvtap)
Central frequency (GHz)a 1.575 5.998 4.134 4.134
uv weightingb Briggs + 10 kλ taper Briggs + 16 kλ taper Briggs Briggs + 16 kλ taper
Clean spatial scales ('') 0, 5, 10 0, 5, 10, 20 0, 7.5, 15, 30c 0, 15, 30, 60
Synth. beam        
('', '', °) 11.43, 10.11, −84.23 6.30, 6.01, −81.30 3.77, 3.40, −3.89 6.81, 6.46, −82.97
Map peak (mJy beam−1) 12.06 8.45 8.41 9.43
rms (μJy beam−1) 20 4.0 7.0 10

Notes. These parameters represent images shown in Figures 58. The meaning of each parameter is as specified in Table 3 unless otherwise indicated. aFrequency of the maps, intermediate between L and C bands. bRobust = 0 was used in each case. The term, "uvtap," means that a uv taper, of specified scale, was applied to the data. cThe largest spatial scale was dropped when this clean was very advanced.

Download table as:  ASCIITypeset image

The power and versatility of CASA, moreover, along with the multi-array and frequency observations that we have carried out, also allow us to form a single image from all arrays and frequencies. The self-calibrated uv data from all five observations are input into the clean algorithm and, with a point-by-point fit of the spectral index, an image can be made that corresponds to a frequency which is intermediate between L and C bands, i.e., at ν = 4.13 GHz. This process is time-consuming but produces a final image with improved sensitivity to spatial scales and, on average, lower rms. We show the result in Figure 5, where the two images differ only in the uv weighting function adopted.

Figure 5.

Figure 5. Total intensity combined array images of UGC 10288 (all arrays, both bands) corresponding to a frequency ν0 = 4.1 GHz, shown at two different uv weightings. (a) I contours of the higher resolution Briggs robust = 0 uv weighting at 14 (2σ), 20, 30, 56, 100, 200, 400, 800, 2000, and 4000 μJy beam−1 over the SDSS r-band color image. The beam size (∼3farcs5; Table 4) is shown at lower left. (b) I contours and color of the lower resolution Briggs robust = 0 + 16 kλ uv tapered image, with contours at 20 (2σ), 25, 30, 38, 55, 80, 160, 300, 600, 1200, 2400, 4800, and 8500 μJy beam−1. The beam size (∼6farcs5; Table 4) is shown at lower left.

Standard image High-resolution image

Since the spectral index, α (Iν ∝ να), is fitted during the imaging, a map of α has also been formed (see Paper 2 for details) and is shown in Figure 10 along with its error map. A 3σ cutoff has been applied to the data when forming the α map and the result has been corrected for the spectral index of the primary beam. The error map represents random errors at each point as described in Rau & Cornwell (2011; their Equation (39)). Note that similar spectral index maps were formed for each individual data set as well (Section 2.3), but we show only the combined array/frequency result because of its higher S/N. In principle, it is also possible to solve for curvature in the spectral index (see Paper 2 for examples); however, the uncertainties become large, given the S/N in the disk of UGC 10288, and we did not pursue this for UGC 10288.

The remaining issue is to consider whether there might be missing flux from a lack of short spacings. At L band, we detect all spatial scales less than 16' and, at C band, less than 4'. The optical galaxy is 4farcm9 in diameter (Table 1) but, as illustrated in Figures 1 and 2, the radio diameter does not extend beyond ≈3'. Moreover, we do not observe any "negative bowls" around the emission, as would be the case if there were missing spacings. Consequently, there should be no missing flux in these data.

3. RESULTS

3.1. The Total Intensity Emission

A remarkable result of these observations is the discovery of a strong, double-lobed EGRS which is most certainly in the background of the UGC 10288 field. Previous to these observations, the only radio continuum images obtained were of the galaxy and background source blended together (see Hummel et al. 1991; Condon et al. 1998). Indeed, the previously measured flux density from the combined sources, S1.4 GHz = 26.1 mJy (Condon et al. 1998), met our minimum flux density criterion (23 mJy) for membership in the CHANG-ES survey (Paper 1); however, it would have fallen far below that cutoff, had a measurement existed for the galaxy alone.

We now see both the galaxy and the background source in detail. There is an optical point source at the core of the EGRS (Figure 5), with the Sloan Digital Sky Survey (SDSS) DR9 identifer, J161423.28−001211.8, listed as a galaxy with a photometric redshift of zphot = 0.388 ± 0.026. We have named this background source CHANG-ES A and labeled it in Figure 1(b). Another strong source (CHANG-ES B, labeled in Figure 6) is blended with the disk of UGC 10288 in D-array L band (Figure 1(a)), making the UGC 10288 disk appear to be skewed toward the northwest in that figure. CHANG-ES B also has an optical background galaxy at its center, namely J161419.79−001155.6 at zphot = 0.39 ± 0.11, in agreement with the redshift of CHANG-ES A. Therefore, CHANG-ES A and CHANG-ES B possibly form a background pair.

Figure 6.

Figure 6. Combined C and D array C-band contours over a gray scale λ3.6 μm map from the Spitzer telescope. Contour levels are at 8 (2σ), 11, 15, 20, 30, 50, 100, 300, 600, 1200, and 3000 μJy beam−1, and the beam size (∼6''; Table 4) is shown at lower left. Several point sources showing radio emission that are distinct from the disk of UGC 10288 have been numbered (Table 6). The box encloses a region containing the radio continuum extension seen in Figure 5(b).

Standard image High-resolution image

3.1.1. Flux Densities and Global Spectral Index

Estimates of the flux densities of the various sources are given in Table 5. All flux densities were measured from the primary-beam-corrected C-array L-band (Figure 1(b)) and D-array C-band (Figure 2(a)) maps. For both bands, the total flux densities of UGC 10288 and CHANG-ES A were measured together. Then the flux density of CHANG-ES A was estimated from a measurement within the two dashed lines shown in Figure 2(a), a width of 34''. The flux density of UGC 10288 alone was then measured in a region of equivalent width to the east of CHANG-ES A along the galaxy's disk. Assuming that this measured flux density is approximately equivalent to the flux density of UGC 10288 within the box occupied by CHANG-ES A, it was then subtracted from the CHANG-ES A measurement as a correction for the contribution of the galaxy's disk within the dashed lines. At L band, this correction amounted to 4% of the flux density of CHANG-ES A, and 22% of the flux density of UGC 10288. At C band, the correction constituted 4% of the flux density of CHANG-ES A and 32% of the flux density of UGC 10288. In addition, at L band, the totals were then summed and compared to the combined, blended emission from all three sources, i.e., UGC 10288, CHANG-ES A, and CHANG-ES B as measured in the D-array L-band image (Figure 1(a)); the totals agree to within 0.5%.

Table 5. Derived Properties of UGC 10288 and CHANG-ES A

Property UGC 10288 CHANG-ES A
R.A. (J2000)a 16h14m24fs80 16h14m23fs28 ± 0fs04
Decl. (J2000)a −00°12m27fs1 −00°12m11.6 ± 0fs5
SL (mJy)b 4.4 ± 0.5 22 ± 2
SC (mJy)b 1.53 ± 0.05 12.0 ± 0.2
αtotalc −0.76 ± 0.11 −0.44 ± 0.08

Notes. aU10288 position from NED. For CHANG-ES A, the position has been measured from the highest resolution data. bFlux densities in L or C bands, measured from the primary-beam-corrected images as described in Section 3. Uncertainties represent variations from adjusting the box size containing the emission. cGlobal spectral index according to SL/SC = (νLC)α, where νL and νC are the L-band and C-band central frequencies, respectively.

Download table as:  ASCIITypeset image

The global spectral index of UGC 10288 was calculated from the flux densities as just described, and the result is given in Table 5. The two sources are both clearly dominated by non-thermal emission.

3.1.2. The Star Formation Rate of UGC 10288

It is remarkable that UGC 10288 reveals extra-planar emission (see the next two subsections) when the total flux density is only a few mJy. This leads us to consider the star formation rate (SFR) of the galaxy, since halo emission is often associated with such activity.

UGC 10288's SFRIRAS of 1.3 M yr−1 as given in Table 1, has been calculated from Infrared Astronomical Satellite (IRAS) data (Paper 1), but, like the radio continuum, the IR emission will certainly be contaminated by emission from CHANG-ES A.

A closer inspection of the IRAS flux densities, moreover, suggests foreground contamination as well. Given the observed IRAS 60–100 μm flux density ratio of UGC 10288 (f(60 μm)/f(100 μm) = 0.236), the galaxy appears to be extremely cold. In fact, this flux density ratio is ≈20% smaller (colder) than that for the most extreme (coldest) galaxy included in the Dale & Helou (2002) spectral energy distribution models of normal star-forming galaxies. If CHANG-ES A is indeed a background radio galaxy, we naively would not expect it to be colder than all galaxies in the local universe. If anything, we would expect a warmer temperature arising from hot dust being powered by the accreting black hole. Therefore, it seems that the most likely explanation for the cold dust temperature is significant contamination by difficult-to-remove foreground cirrus emission that is contributing to the IRAS flux densities.

Given these issues, an alternative approach is to use the calibrated, known radio continuum–far-IR (FIR) relation to estimate the SFR of UGC 10288 from its radio continuum emission as estimated in Section 3.1.1 with CHANG-ES A's flux density subtracted. For this, we adopt the calibrations of Murphy et al. (2011) which use a Kroupa initial mass function and a supernova cutoff of 8 M. The SFR estimated from the thermal radio continuum requires an electron temperature which we take to be 104 K. The SFR estimated from the non-thermal radio continuum requires a knowledge of the non-thermal spectral index, αNT, whereas we have, from the observed flux densities of Table 5, a global spectral index of α = −0.76, using the convention (Sν ∝ να).22

However, we can estimate αNT by combining the thermal and non-thermal calibrations of Murphy et al. (2011, their Equations (11) and (14)) and the constraint that the measured flux density is the sum of the non-thermal and thermal flux densities at any frequency, i.e.,

Equation (1)

where Sν is the total flux density at frequency ν, and these calibrations assume a linear dependence of SFR on both thermal and non-thermal luminosity. The use of this equation infers a thermal/non-thermal fraction which has a dependence on both frequency and αNT.

Using our measured flux densities of Table 5 at two frequencies, we can solve Equation (1) to find αNT = −1.04. This leads to a thermal to total and non-thermal to total flux density in L band of (ST/S)1.5 = 0.17 and (SNT/S)1.5 = 0.83, respectively. In C band, the results are (ST/S)6.0 = 0.44 and (SNT/S)6.0 = 0.56, respectively.

Finally, with the SFR calibration (Equation (15) of Murphy et al. 2011), we find SFR = 0.51 M yr−1.

We caution that, if the non-thermal luminosity, LNT ∝ SFRn, where n > 1 (e.g., n = 1.3; Niklas & Beck 1997), then there would be a weak dependence on SFR on the right-hand side of Equation (1); the adjustment, however, does not significantly perturb the above results.23

This newly determined SFR is considerably lower than the original SFR which was found from IRAS values and should be  more accurate since it corrects for the contribution of CHANG-ES A and also uses a spectral index determined from the data. If, instead, we use the global FIR–radio continuum relation (Equation (17) of Murphy et al. 2011) based on earlier calibrations (de Jong et al. 1985; Helou et al. 1985), then SFR = 0.39 M yr−1. In either case, the radio-derived SFRs are clearly much lower than the IRAS-based estimate of 1.3 M yr−1.

In summary, the SFR based on IRAS photometry (Table 1) neither corrects for the strong contribution from CHANG-ES A nor for contamination from foreground cirrus. We suggest that a value between 0.4 and 0.5 M yr−1 is a more accurate estimate of the true SFR in UGC 10288. Note that the range given here is dominated by the assumptions that have been made, rather than by a formal error determination. This galaxy is now one of the lower SFR galaxies in the CHANG-ES sample (Paper 1).

3.1.3. Foreground/Background Point-like Sources

Figure 6 shows the combined C+D array C-band map of UGC 10288 overlaid on the λ3.6 μm image from the Spitzer Space Telescope.24 These radio data have the lowest rms of all combined data sets (Table 4) as well as high resolution. Comparing with the Spitzer image, which was chosen because it should show both foreground stars as well as background quasars, should reveal how much radio emission (especially in regions of apparent halo features) is associated with point-like sources in the displayed field.

We searched for sources for which both the IR and radio emission are point-like and for which the radio emission is sufficiently separated from the disk of UGC 10288 to allow a discrete measurement of flux density. The process was then repeated for the combined B+C array L-band map (not shown), although the latter image has almost four times higher rms noise than at C band for approximately equal beam sizes. The resulting sources are labeled in Figure 6 and their properties listed, along with SDSS identifiers (where available) in Table 6.

Table 6. Point-like Sources with Radio Emission near UGC 10288

Source R.A. (J2000)a Decl. (J2000)a SDSS Identifierb SLc SCc αd
(h m s) (° ' '') (μJy) (μJy)
1 16 14 29.88 −00 12 01.2 J161429.87−001201.3 (S) 50 ± 10 9 ± 4 −1.3 ± 0.6
2 16 14 28.80 −00 12 13.6 ND <16 16 ± 4 >0
3 16 14 27.23 −00 13 21.0 ND 35 ± 8 17 ± 4 −0.5 ± 0.4
4 16 14 25.43 −00 12 12.1 J161425.40−001212.2 (S) <16 16 ± 4 >0
5 16 14 24.87 −00 13 18.5 J161424.84−001318.3 (S) 109 ± 5 53 ± 4 −0.5 ± 0.1
6 16 14 21.40 −00 12 36.6 ND 43 ± 8 17 ± 4 −0.7 ± 0.3
CHANG-ES B 16 14 19.80 −00 11 55.6 J161419.79−001155.6 (G) 3290 ± 50 1040 ± 20 −0.83 ± 0.02
7 16 14 19.04 −00 11 34.4 ND 28 ± 6 12 ± 4 −0.634 ± 0.4

Notes. aMeasured from the λ3.6 μm map at the maximum value (Figure 6). The uncertainty is approximately 0farcs8 (one-half of the spatial resolution). bSDSS identification number. "ND" means no detection in the SDSS. (S) = "Star"; (G) = "Galaxy." cC-band and L-band flux densities were measured from the primary-beam-corrected combined C+D array images and the combined B+C array images, respectively. Uncertainties represent a range resulting from varying the measurement box size. dSpectral index, determined from the previous two columns.

Download table as:  ASCIITypeset image

Let us consider the region that is encompassed by the box shown in Figure 6 (16h14m25fs9 ⩽ R.A. ⩽ 16h14m28fs9; −00°12'21'' ⩽ decl. ⩽−00°11'40''). This box contains significant "halo" emission in the form of a large extended arc-like feature visible in the all-data image of Figure 5(b). This arc extends to a height of 21'' (3.5 kpc) above the plane. The flux density contained within this box is S4.1 GHz = 450 ± 80 μJy.25 Only one identified point-like source (Source 2) is located within this region whose flux density, extrapolated to 4.1 GHz, is <16 μJy (Table 6). Consequently, although we cannot conclude that the arc-like feature represents a single coherent structure, its flux density does not result from background sources. Hereafter, we will refer to this features as the "arc" for simplicity.

There is also a feature extending to the north on the west side of the galaxy (at R.A. ≈ 16h14m21fs5, decl. ≈−00°12'12'') as shown in Figure 5(b), but there are no point sources meeting our criteria in this region at all.

We caution that if a background source displays double-lobed radio structure, then it could be displaced with respect to the IR source and might have been missed in this search, as would objects with weak IR flux and those that are blended with the disk. However, even if we increased our point-like source flux density estimate by a factor of two, its contribution to the flux density of halo features in UGC 10288 would still be negligible. By including other array configurations in the maps, in particular short spacings, the sensitivity to extended halo features is enhanced but there should be no change in sensitivity to point sources.

It is beyond the scope of this paper to examine in detail the nature of the point-like sources. SDSS identifiers, where available, indicate that three of the four identified sources are stars, which would be an interesting result given the detected radio emission. However, one of these sources (Source 5, SDSS J161424.84−001318.3), although labeled as a star, has SDSS colors (u − g = +0.18, g − r = −0.18, r − i = 0.45, and iz = −0.07) which are entirely consistent with those of a quasar over a range of redshift (Pâris et al. 2012, their Figure 28).

3.1.4. Comparison with Images at other Wavebands

It is apparent that halo radio emission is associated with UGC 10288. Here, we appeal to ancillary data, where available, to shed further light on the high-latitude features.

Few resolved images of UGC 10288 are available at other wavebands but, of those available, we show several overlays in Figures 7 and 8. Figure 7 shows our all-array and all-frequency radio continuum images of Figure 5 superimposed on the Hα image of Rand (1996). In Figure 8, we show contours from the Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010) at λ12 μm and λ22 μm over Hα emission and our L-band combined array image, respectively. In each overlay, the spatial resolutions of the two images are identical.

Figure 7.

Figure 7. Contours from Figure 5 over the Hα emission measure map of Rand (1996) smoothed to the same resolution as the radio data in each case. The beam size is shown at bottom left. (a) Contours are from Figure 5(a) and the resolution is ∼3farcs5 (Table 4). The arrows point to features discussed in Section 3.1.4. (b) Contours are from Figure 5(b) and the resolution is ∼6farcs5 (Table 4). A plus marks the center of UGC 10288 and two crosses mark the positions discussed in Section 3.3.1. The labeled strips are discussed in Section 3.1.4.

Standard image High-resolution image
Figure 8.

Figure 8. (a) Contours (in arbitrary units) of the WISE λ12 μm image over the Hα image of Rand (1996), both smoothed to 7farcs0 resolution (beam shown at lower left). (b) Contours (in arbitrary units) of the WISE λ22 μm image over the combined BCD array L-band data (original information given in Table 4) smoothed to the same 12''resolution as the WISE data (beam at lower left).

Standard image High-resolution image

Because radio emission from UGC 10288 is blended with that of CHANG-ES A, it is difficult to disentangle which source is responsible for specific features observed at locations close (in projection) to that background source. For example, a horizontal feature appears to emanate from the northern lobe of CHANG-ES A toward the east (see R.A. ∼ 16h14m25s, decl. ∼−00°12'06'' in Figure 7(b)). There is also a projection which extends to the south of CHANG-ES A's southern radio lobe (R.A. ∼ 16h14m23s, decl. ∼−00°12'43'' in Figure 7(b) and farther south as a narrower radio extension in Figure 8(b)). Neither of these have obvious counterparts at other wavebands, although the southern feature appears to lie between two Hα extensions (see below).

The IR emission does not distinguish between UGC 10288 and background sources either, since, for example, the bright background radio source, CHANG-ES B, to the west of CHANGES-A (marked in Figure 6) is also an IR emitter at both λ12 and 22 μm.

As we have no way to determine the origin of radio continuum features near CHANG-ES A, we will concentrate instead on other features that are farther from this region and are more likely to be associated with UGC 10288. The Hα emission is helpful in this regard, since it clearly has an origin in UGC 10288 itself.

We must also note that each of the images suffers from some low-level artifacts.

For the radio images, the map contains residual side lobes which, for the most part, we have included in measuring and quoting rms map values. However, residual side lobes from a source 6farcm3 to the east of UGC 10288 are likely responsible for the horizontal feature at the southeast side of Figure 7(b) (R.A. ∼ 16h14m30s, decl. ∼−00°12'50''), and residual side lobes from another point source 6farcm2 to the southwest (SW) may be responsible for the faint gray scale NE–SW feature visible on the west side of Figure 8(b) (R.A. ∼ 16h14m17s, decl. ∼−00°12'15'').

Similarly, in the Hα image, it has been noted that the galaxy's bulge has been oversubtracted (Rand 1996). And both WISE images show broad (arcmin scale) NE–SW striping which may be producing or contributing to the broad southward extension at R.A. ∼ 16h14m22.5s, decl. ∼−00°13'15'' (Figure 8(a)).

In any following discussion, therefore, we will keep these cautions in mind and discuss results that appear to be robust to such artifacts. Figure 7(a) shows the highest resolution (∼3farcs5) radio/Hα overlay with arrows pointing to specific features. Discrete in-disk Hα emission reveals specific sites of unobscured star formation.

Feature 1 is a discrete star-forming region in the disk above which there is a vertical plume of Hα emission (to the north). The Hα plume is more readily seen in the smoothed map of Figure 7(b) and has a λ12 μm counterpart (Figure 8(a)). No radio emission can be confirmed in this plume, but a slight extension visible in the L-band gray scale image of Figure 8(b) is close to its location, slightly displaced to the east.26 That is, the radio extension's location is more closely aligned with the gap27 immediately to the east of star-forming region #1.

Feature 2 refers to several star-forming complexes that are located below the strong northward radio plume seen in Figure 7. The two arrows point to gaps in the Hα in-disk emission, the easternmost one also corresponding to a gap in the in-disk radio emission. Two "prongs" in the radio emission close to the disk (Figure 7(b)) are immediately above these gaps.

The southern Hα extension of feature 3 was originally pointed out by Rand (1996) and is in a region of low brightness in comparison to the star-forming regions to the east and west of it. This feature shows a small northward elongation in the radio at high resolution (Figure 7(a)).

Feature 4 points to another gap in the Hα emission between two star-forming regions. There is a small radio extension to the north that is visible in the high-resolution image (Figure 7(a)) but a much larger corresponding radio continuum feature can be seen at lower resolution (Figure 7(b)).

Finally, feature 5 refers to two southward prongs which were originally identified by Rand (1996) who also noted a "hemispherical" shape formed by these two features together. We have not detected any high-latitude radio emission associated with this feature.

In summary, our high resolution data suggest that radio continuum features are more likely to be seen above in-disk gaps in the Hα emission.

It should be pointed out, of course, that apparent correlations above an edge-on disk could apply to any location along the line of sight (although such indeterminacy is less severe at large galactocentric radii). This means that such apparent correlations do not necessarily represent physical correlations. Nevertheless, the observed trend is worth pointing out, since it is rare that an edge-on galaxy has adequate data to compare Hα and radio continuum results at the same resolution with good sensitivity. A similar anti-correlation has been observed before between in-disk Hα emission and high-latitude radio continuum features in NGC 5775 (Lee et al. 2001). We discuss this further in Section 4.

On broader scales, it is more difficult to associate high-latitude radio emission with in-disk activity. Intriguing features are a southward Hα arc seen in Figure 7(b) at R.A. ∼ 16h14m24fs5, decl. ∼−00°12'48'', and a possible companion feature at R.A. ∼ 16h14m22fs7, decl. ∼−00°12'42'' between which the southern radio continuum extension, mentioned above, lies. In Figure 8(a), the λ12 μm emission shows several other northward extensions as well as a southward arc-like feature on the east side of the galaxy, and in Figure 8(b), a southward radio extension seen at L band on the far east end of the disk.

We note two other observations regarding the distributions.

The first is the north–south radio asymmetry observed on the east side of the major axis, best seen in Figure 7(b). That is, there is a large arc-like extension toward the north but nothing equivalent to the south.28 This may reflect N/S asymmetries in star-forming regions and/or N/S density asymmetries in the disk. Such asymmetries can also be produced by ram pressure as a galaxy moves through an intergalactic medium; this is possible for UGC 10288, given that it is in a sparse group (Section 1), but we have no independent observations that would argue for or against this speculation.

The second observation is that UGC 10288 does not appear to have a global radio halo, but rather it has discrete high-latitude radio extensions. Table 7 compares beam-corrected exponential scale heights between the 4.1 GHz all-data map (Figure 5(b)) and both the Hα and WISE λ12 μm data at the same (7'') resolution according to the method of Dumke et al. (1995). The fitting has been carried out for data that have been averaged in strips, each 26farcs4 = 1.76 s wide, that are labeled in Figure 7(b).29 The values are from single-exponential fits that apply to the disk (for example, for Strip #2 North, the horizontal feature at decl. ∼ −00°12'06'' was avoided). Note that the strips that do not have significant extensions (mainly #2 South but also #1 South and #2 North) have low to modest scale heights in radio continuum (0.46, 0.85, and 0.86 kpc, respectively), whereas the strips containing discrete extensions (#1 North and #4 North/South) are wider in z on average (1.3, 1.4, and 1.1 kpc, respectively), as would be expected (see Section 4 for more discussion).

Table 7. Disk Exponential Scale Heights

Strip # BCD+CL WISE
(4.1 GHz) λ12 μm
1 (North) 7.9 (1.3) 2.3 (0.38) 6.4 (1.1)
  (South) 5.1 (0.85) 5.1 (0.84) 7.7 (1.3)
2 (North) 5.2 (0.86) 3.2 (0.54) 6.8 (1.1)
  (South) 2.8 (0.46) 4.0 (0.66) 6.8 (1.1)
4 (North) 8.6 (1.4) 2.9 (0.47) 7.5 (1.2)
  (South) 6.7 (1.1)a 4.0 (0.66) 7.1 (1.2)

Notes. Comparison of beam-corrected exponential scale heights between the radio emission from combined array, combined frequency data (4.1 GHz), the Hα emission measure map, and the WISE λ12 μm map, all at 7'' resolution. Corrections for non-zero baselines have been made. Values are in arcsec with kpc in parentheses. aThis value reduces to 4.6 (0.76) if the fit is cut at the first "shoulder" of emission.

Download table as:  ASCIITypeset image

In order to give an overview of the relative spatial distributions of the various sources of emission, we finally display a multi-frequency image (Figure 9) in which several data sets have been overlaid with different colors (see caption). Our radio data show CHANG-ES A and CHANG-ES B in cyan; our masking technique (Rector et al. 2007) allows us to show their optical counterparts (yellowish) at their cores. The disk of UGC 10288 looks quite wide vertically in this picture because of the WISE λ12 μm image. For more information as to the techniques used to combine the data sets, see Rector et al. (2007).

Figure 9.

Figure 9. Multi-frequency color image of the UGC 10288 field showing our C-array C-band image (Figure 2(b)) in cyan, the combined all-array, all-frequency image in darker cyan (Figure 5(b)), the WISE λ12 μm image (Figure 8(a)) in orange, the Spitzer λ3.6 μm image (Figure 6) in yellow, the Hα image of Rand (1996; Figure 7) in rose, the SDSS r-band image (Figure 5(a)) in blue, and the SDSS g-band image in purple. Smoothing, color choice, and masking follow the techniques described in Rector et al. (2007). Spatial resolutions vary and have been chosen for visual effect. Both CHANG-ES A and CHANG-ES B (cyan and darker cyan at lower resolution) are quite obvious, and the large vertical width of the disk which appears salmon colored at the edges, is due to the large width seen in the WISE λ12 μm image.

Standard image High-resolution image

3.2. Spectral Index

The spectral index map along with a point-by-point error map (see Section 2.5) are shown in Figure 10. The disk of UGC 10288 is of course blended with CHANG-ES A which shows the classic signature of a double-lobed radio source: a flatter spectral index at the center and steeper in the lobes.

Figure 10.

Figure 10. Spectral index map, α (a), and uncertainty in the spectral index (b) (random errors only), made from all data (all arrays, both bands) corresponding to the uv weighting image of Figure 5(b) (Sν ∝ να) with a 3σ intensity cutoff applied. The three brightest star-forming regions from Figure 7(a) have been marked with crosses and the beam is shown at lower left. Several points from the α map have uncertainties that appear white because they have values that are >1. The beam size (∼6farcs5; Table 4) is shown at lower left.

Standard image High-resolution image

The spectral index in the disk of UGC 10288 shows complex structure but, as also shown by the global spectral index (Section 3.1.1), non-thermal spectra dominate (Table 5). At the locations of the three brightest star-forming regions in the disk (denoted with crosses), as identified from the Hα map, the spectral indices are α = −0.7 ± 0.3, α = −0.4 ± 0.2, and α = −0.1 ± 0.3 from east to west, respectively. If averaged over a beam rather than measured at the Hα peak, these values become α = −0.8 ± 0.3, α = −0.4 ± 0.2, and α = −0.3 ± 0.2 from east to west, respectively. Thus, only the westernmost star-forming region is consistent with a purely thermal spectrum at its peak. However, if we insist that the SFR calibrations discussed in Section 3.1.2 hold at every point, then a flatter observed spectral index actually implies a flatter non-thermal spectral index, presumably reflecting a younger cosmic ray electron population. Such details will be explored in future CHANG-ES papers.

Note that these three star-forming regions do not align exactly with the flattest spectral indices. The region with the flattest spectral index is directly north of the central star-forming region and there is another region of flatter spectral to the southwest of the westernmost star-forming region. This could be a result of outflows in these regions; such spectral index flattenings have been observed before in regions of outflows (e.g., NGC 5775; Lee et al. 2001). However, it is important to keep in mind the uncertainties on these quantities. Generally, the errors also increase toward the edges of the map.

3.3. Polarization

3.3.1. Polarization Images and Sample Magnetic Field

Figures 3 and 4 show the polarization maps at L and C bands, respectively. We do not detect polarization in the disk of UGC 10288 in the individual observations. Rather, it is the background radio source, CHANG-ES A, that shows polarization. Although a more complete analysis of the polarization-related properties awaits a future paper, a few results are presented here.

For both B-array, L band (resolution ∼ 3farcs6) and C-array, C band (resolution ∼ 3farcs0), the highest polarization occurs right at the core of CHANG-ES A (peak values and resolutions provided in Table 3). At this location, the percentage polarization at L and C bands is 1.9% and 2.0%, respectively.

North and south of the core, although the polarized intensity is lower, the percentage polarization increases. At the peak polarization position of the C-band data in the northern lobe, for example, the values are 17% and 30% at L and C bands, respectively (R.A. = 16h14m23fs7, decl. = −00°11'48''). At the southern peak position, the corresponding values are 6.3% and 31% at L and C bands, respectively (R.A. = 16h14m23fs47, decl. = −00°12'21'').30

The L-band percentage polarization is weaker than at C band, as expected, since more Faraday depolarization occurs at the lower frequency, either intrinsic to the source or contributed via UGC 10288 or both. Moreover, the L-band value is lower still for the southern lobe, compared to the northern, likely because of stronger Faraday depolarization in the (possibly turbulent) disk of UGC 10288. The fractional polarization in the core is much lower than in either lobe; this is not unusual for double-lobed radio sources (e.g., Rudnick et al. 1986; Saikia & Kulkarni 1998) suggesting that much of the core/lobe difference is intrinsic to CHANG-ES A.

UGC 10288 does not lend itself well to global magnetic field calculations, given the interfering presence of CHANG-ES A. However, we have taken two sample positions, one in the disk, and one in the large northern arc on the east side of UGC 10288 visible in the 4.1 GHz map, to estimate the total minimum energy magnetic field strength, B. We use the revised minimum energy formulae given in Beck & Krause (2005) for isotropic fields and adopt a proton-to-electron number density ratio of 100. The positions are marked with crosses in Figure 7(b), the disk position being immediately below the arc.

At the position in the eastern disk, a projected distance of x = 39'' from the center, α = −0.44 ± 0.18 and I4.1 GHz = 70.4 μJy beam−1. Then αNT = −0.54 and INT = 55.4 μJy beam−1 by Equation (1) and the arguments presented in Section 3.1.2. For a line-of-sight distance of 6.0 kpc (estimated from a 1/e radial scale length) then the minimum energy magnetic field strength is B = 10 ± 3 μG. Here, the error represents measurement error only, which is dominated by the uncertainty in α; it does not include errors that may be associated with the assumptions.

At the position in the northern extension, 15'' above the plane (2.5 kpc in projection), α = −1.00 ± 0.50 (Figure 10), I4.1 GHz = 36.0 μJy beam−1, we assume that the emission is dominantly non-thermal based on the absence of Hα emission, and take l = 1.5 kpc, which is approximately the width of the arc at this location (9''). Then, provided that the conditions for minimum energy are met for this region of steep spectral index, B = 12 ± 5μG, where, again, the error represents measurement error.

As Figure 10 illustrates, there is considerable variation in spectral index (and its error) and the magnetic field strength will also vary with position. It is also not clear whether the thermal/non-thermal emission ratios follow the expectations of the calibrations of Murphy et al. (2011), point by point. Nevertheless, the approach in which the thermal/non-thermal ratio is tied to an SFR calibration should be an improvement upon previous estimates for other galaxies which usually assume a constant and singular thermal/non-thermal ratio for every point in the galaxy. These preliminary calculations, then, suggest a total magnetic field strength of order 10 μG at the two sample positions selected.

3.3.2. Rotation Measures and the Regular Field

A powerful tool, which is potentially available for UGC 10288, is the technique of using the background radio galaxy, CHANG-ES A, to probe the foreground magnetic field. To this end, we computed the classical Faraday rotation measure (RM) map (Figure 11) between the B-array L-band and the C-array C-band data, at 5'' resolution. The maximum 1σ error is ±10 rad m−2 at the outer edge of each patch, while the 1σ error of the mean RM of Figure 11 is ±2 rad m−2. The RM ambiguity is ±93 rad m−2.

Figure 11.

Figure 11. Rotation measure map, made from the L-band B-array and the C-band C-array data, at 5''resolution, computed at all pixels where the polarized intensities in both bands exceed twice the rms noise. Several contours (red) are marked, with all values in units of rad m−2. The two regions with no contours marked have average values of 33 rad m−2 (northern region) and 62 rad m−2 (southern). Letters denote the regions discussed in Section 3.3.2. Black contours represent the 5'' C-band C-array total intensity image at 6.2, 30, 150, 300, 1200, and 6000 μJy beam−1. The beam is marked at lower left.

Standard image High-resolution image

The observed RM values are a superposition from CHANG-ES A, the galaxy UGC 10288, and the foreground of the Milky Way. The lobes of radio galaxies can have high RM values with a large dispersion, so that, when observed with low resolution, the average RM is small. The central sources of radio galaxies (jets) generally do not reveal significant intrinsic RM. The RM of CHANG-ES A's northern lobe (at about 7 kpc projected distance above the plane of UGC 10288, marked "A" in Figure 11) of +30 ± 4 rad m−2 (average of the two visible regions) agrees well with the foreground RM of the Milky Way (Taylor et al. 2009; Oppermann et al. 2012) at the Galactic coordinates of our radio maps (l = +12°, b = +34°). This indicates that no correction for the RM ambiguity is needed.

The southernmost blob in Figure 11 (marked "C") is located behind the disk of UGC 10288, about 1.5 kpc above the plane, with an observed average value of RM ≃ +62 ± 4 rad m−2. We interpret the (foreground subtracted) value of RMi ≃ +32 rad m−2 as indication of a regular magnetic field in the galaxy disk, pointing toward us. A line-of-sight component of about +2 μG, a thermal electron density of 0.02 cm−3, and a path length of 1 kpc can provide the observed RMi. If the path length is more than 1 kpc in the halo, the field strength and/or the electron density lowers by the same factor. However, as noted earlier (Section 3.1.4), when high-latitude radio emission is observed, the features have generally been discrete with typical sizes of the order of a kpc. In addition, disk exponential scale heights are also typically ≈1 kpc. Consequently, at the location of "C" at 1.5 kpc above mid-plane, we might expect the effective path length of the magnetic field to be approximately that of a typical high-latitude discrete feature.

RMi ≃ −7 rad m−2 (observed average value of RM ≃ +23 ± 2 rad m−2 for both parts) toward the central region of the radio galaxy, at about 3.5 kpc above the plane of UGC 10288 ("B"), indicates that the field direction is reversed with respect to the plane. A thermal density in the halo of 0.01 cm−3 requires a regular halo field with a line-of-sight component of about −1 μG, pointing away from us. Again, if the path length is more than 1 kpc in the halo, the field strength and/or the electron density lowers by the same factor. Between disk and halo, no regular field from the galaxy is detected, so that the RM is solely due to the Galactic foreground.

The conclusion that a field reversal is present depends, of course, on the accuracy of the RM measurements, especially region "C" which is present in the total power maps only as diffuse emission between the core of CHANG-ES A and its southern radio lobe. If we use a higher cutoff of 3σ to form the RM map, then the value of "C" becomes 64 rad m−1, rather than the 62 rad m−2 seen in Figure 11, i.e., within the quoted 1σ uncertainty of ±10 rad m−2 on the RM at the map edges (see above). A 3σ cutoff also reduces the 1σ uncertainty in RM to ±7 rad m−2; however, such accuracy is not required for this conclusion.

A reversal between the azimuthal components of the disk field and halo field may be quite common among galaxies; it has been found from modeling the magnetic field in the Milky Way from extragalactic RMs (e.g., Sun et al. 2008; Jansson & Farrar 2012) and also in the northern part of the galaxy M51 (Fletcher et al. 2011). Dynamo models can hardly explain such reversals if the same dynamo process operates in disk and halo, unless wind outflows allow such "mixed parities" (Moss et al. 2010).

These results are necessarily preliminary. The fact that many spectral channels are available in each of our observations (Table 2) lends itself well to the application of the RM Synthesis technique (Brentjens & de Bruyn 2005) to UGC 10288. This requires that (R.A.–decl.–frequency) cubes be formed for all data sets, and is an approach that will be attempted in the future.

4. DISCUSSION

Our observations have revealed some new surprises in UGC 10288 and the multi-array, multi-frequency approach has been essential in understanding this galaxy and disentangling its emission from the background radio galaxy. Not only have we been able to correct for the emission of CHANG-ES A, but CHANG-ES B would also have been assumed to be part of UGC 10288 at low resolution (see Figure 1(a)), had multiple arrays (i.e., higher resolutions) not been used.

Moreover, the combination of all arrays and both frequencies has allowed us to detect high-latitude emission in UGC 10288. However, rather than having a broad scale halo, the high-latitude features in UGC 10288 are discrete and, in the case of the arc-like feature seen extending toward the north on the eastern side of UGC 10288's disk (Figure 5(b)), also large. This feature extends to 21'' (∼3.5 kpc) above the disk.

It is interesting to explore some properties of this arc. For example, taking the estimate of the total magnetic field strength (B ∼ 10 μG; Section 3.3.1) at the position marked in Figure 7(b) (15'', or 2.5 kpc above the plane) and using standard formulae (e.g., Heesen et al. 2009), the particle lifetime to synchrotron losses is tsynch ∼ 1.6 × 107 yr. We do not know the outflow velocity, but it must be substantial in order to detect synchrotron radiation at this z height. For example, if there were no additional losses, this lifetime implies (for constant velocity) a vertical flow of vz ∼ 134 km s−1, and higher still if adiabatic or other losses are significant. The largest uncertainty, apart from those introduced by assumptions, belongs to α which affects B by ± a few μG. Since tsynch ∝ B−3/2, then the uncertainty on vz due to the uncertainty on α alone is ∼30%. What we see, therefore, is that outflow speeds can be substantial in localized regions, even in galaxies with modest SFRs.

As pointed out in Section 3.1.4, the disk exponential scale heights in regions outside of the discrete high-latitude features is less than 1 kpc (Section 3.1.4, Table 7) and, even when one considers the scale height averaged over regions with and without radio extensions, the average is only ≈1 kpc. This lack of a global radio halo, we suggest, is consistent with a rather modest SFR in UGC 10288 which we have modified downward to be 0.4–0.5 M yr−1. That is, discrete extensions may be formed in relation to underlying disk activity, but not globally. A comparison could be made to M 31 which has a similar (slightly lower) SFR of 0.3 M yr−1 (Arshakian et al. 2011) and for which the synchrotron exponential scale height is only ∼300 pc (Fletcher et al. 2004).

It is clear that galaxies with similar global SFRs can be quite different in other ways. For example, although the SFRs of M 31 and M 33 are similar, their star formation efficiencies, SFE = SFR/M(H2), where M(H2) is the molecular gas mass, differ by at least a factor of three (Tabatabaei & Berkhuijsen 2010; Gardan et al. 2007, their Figures 17 and 11, respectively).31 Similarly, galaxies with comparably low SFRs might show differences in the presence or absence of discrete high-latitude radio continuum features. Possible reasons might be the degree of clustering of star-forming regions, or the relative strength and direction of the magnetic fields.

When the various scale heights are averaged, irrespective of location, we find beam-corrected scale heights of he(4.1 GHz) ∼ 1.0 kpc, he(ne) ∼ 1.2 kpc, where we have doubled the observed average Hα scale height to take into account the fact that the observed emission measure $\propto n_e^2$, and he(12 μm) ∼ 1.2 kpc. Consequently the warm dust and thermal electron density scale heights are comparable and, arguably, so is the radio scale height, although the latter is far more irregular as Figure 7(b) (from which the values of Table 7 were obtained) readily shows.

The high thermal fraction in UGC 1028 (44% in C band; Section 3.1.2) also has been seen before, again in galaxies with low or modest SFRs. For example, Dumke et al. (2000) find a C-band thermal fraction as high as 54% for the galaxy, NGC 5907, and Chyzy et al. (2007) find thermal fractions of 50%–70% in three late-type galaxies at 4.85 GHz compared to what is normally measured (≈20%). These authors suggest that, in galaxies with higher SFRs, the non-thermal fraction increases because of the non-linear dependence of non-thermal luminosity on SFR, as noted in Section 3.1.2. Our results depend on an assumption of a linear dependence, but as already noted, we do not anticipate a significant change at the low SFR of UGC 10288 unless the non-thermal calibration is also markedly different. Consequently, our results of a high thermal fraction are consistent with previous estimates for galaxies with similar SFRs.

A significant new result is, of course, the detection of a double-lobed radio galaxy, CHANG-ES A, behind the disk of UGC 10288, which has allowed us to make some preliminary estimates of the regular magnetic field in UGC 10288 prior to a full RM synthesis analysis. Aside from the Milky Way, there have been only a few limited examples in which polarized background sources have been used to probe foreground galaxies (Han et al. 1998; Mao et al. 2012) and these were for galaxies of large angular size (M 31 and the Large Magellanic Cloud). In the case of UGC 10288, which is only 4farcm9 in diameter, the 90° alignment of the double-lobed radio galaxy is so fortuitous that both the disk and halo regions of the foreground galaxy can be probed.

For an adopted electron density of ne = 0.02 cm−3, a projected distance 1.5 kpc above the plane (position "C" in Figure 11), and a path length of 1 kpc, we find a regular field of +2 μG parallel to the line of sight. At 3.5 kpc above the plane (position "B" in Figure 11) with ne = 0.01 cm−3, B|| ∼ −1 μG, suggesting a reversal of the azimuthal field above the plane.

These values are still preliminary. Nevertheless, CHANG-ES A has the potential to be a powerful probe of both B|| and ne in regions in which neither are measured directly via radio emission. At the lower position ("C"), for example, the Hα emission has already dropped to the level of the noise (Figure 7(a)). This map has a 3σ emission measure noise limit of 13 pc cm−6 (Rand 1996) which implies that any electron density less than ne = 0.1 cm−3 over the same path length would not have been detected (adjusting as $1/\sqrt{l}$ for a different path length, l).

At the same time, our sampling of the total magnetic field strength in the discrete northern arc, 2.5 kpc above the plane, reveals a rather strong total magnetic field B ∼ 10 μG (Section 3.3.1). Although not at the same position as RMs were measured, the implication is that field strengths in discrete high-latitude features could be similar to those in the disk, possibly because of field compression in the arc. A similar result was found for NGC 5775 (Soida et al. 2011) for which X-shaped extensions have fields almost as strong as in the disk, although the latter galaxy has, by contrast, a high SFR and global radio halo. NGC 4631 also has a relatively strong halo field (Hummel et al. 1991). It is unlikely that discrete strong high-latitude features could be a result of a mean-field dynamo alone. Localized star-forming activity in the underlying disk is likely contributing to this feature.

5. CONCLUSIONS

We have obtained JVLA wide-band radio continuum data of UGC 10288 in both total power as well as polarized intensity as part of the CHANG-ES program (Paper 1). Three different array configurations (B, C, and D) at two frequencies, 1.5 GHz (L band) and 6 GHz (C band), have been used and careful attention has been paid to the techniques needed to reduce wide-band data. Moreover, we have combined all total intensity data to form a single map along with its spectral index, applicable to the intermediate frequency of 4.1 GHz. The combined data image has revealed new features in UGC 10288 that were too faint to be seen in the individual data sets.

Our main conclusions are as follows.

  • 1.  
    We have discovered a background double-lobed EGRS, which we have named CHANG-ES A (zphot = 0.39) behind and perpendicular to the disk of UGC 10288. In previous observations, this source had been blended with UGC 10288. Approximately 50'' to the west of this source is another background source, CHANG-ES B, at the same redshift.
  • 2.  
    We have disentangled the flux densities of UGC 10288 and CHANG-ES A and find, for the foreground galaxy, that it is only 17% and 11% of the total blended flux density at L and C bands, respectively (Table 5). UGC 10288 would not have been in the CHANG-ES survey had it not been blended with CHANG-ES A, since its flux density is a factor of five below the minimum cutoff for the survey.
  • 3.  
    We have revised the SFR of UGC 10288 downward from the IRAS value of 1.3 M yr−1 to our new value of 0.4–0.5 M yr−1. In the process, we have employed the SFR calibrations of Murphy et al. (2011) for both thermal and non-thermal radio emission, which also allows us to solve for the global non-thermal spectral index (αNT = −1.0).
  • 4.  
    The global thermal fractions tend to be high, as has been observed before in galaxies which have lower SFRs, i.e., (ST/S)1.5 = 0.17 and (ST/S)6.0 = 0.44, at 1.5 GHz and 6 GHz, respectively.
  • 5.  
    UGC 10288 shows discrete high-latitude radio continuum features. In particular, a large radio continuum arc-like feature is observed to the north of the major axis on the east side of the galaxy (Figure 5(b)). This arc extends to 3.5 kpc above the plane.
  • 6.  
    At high resolution, there is a trend for high-latitude radio emission to occur above gaps in the underlying Hα emission; however, we cannot confirm whether this relationship is a physical one because of uncertainty in the line-of-sight location of the features.
  • 7.  
    UGC 10288's radio continuum emission does not form a disk-wide global halo. The beam-corrected exponential scale heights are typically ∼1 kpc when averaged over regions both with and without discrete high-latitude radio continuum features (smaller in regions which do not show such features). High-latitude radio emission therefore appears to occur in localized discrete features in this galaxy.
  • 8.  
    The fortuitous placement of CHANG-ES A has allowed us to do a preliminary RM analysis of UGC 10288 at several positions. Between 1.5 and 3.5 kpc above mid-plane, the azimuthal components of magnetic field appear to reverse directions.
  • 9.  
    The minimum energy total magnetic field strength in the large northeastern arc at z = 2.5 kpc is ∼10 μG at a location 2.5 kpc above the plane, indicating that the magnetic field is substantial in this high-latitude feature.

We thank Jeremy MacHattie for retrieving previous radio flux density information at the time of proposal submission. J.A.I. would like to thank the Natural Sciences and Engineering Research Council of Canada for grant. This research has used the Karl G. Jansky Very Large Array operated by the National Radio Astronomy Observatory (NRAO). NRAO is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This work is based in part on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the US Department of Energy Office of Science. The SDSS-III Web site is http://www.sdss3.org/. R.B. acknowledges support by DFG FOR1254.

Facility: VLA - Very Large Array

Footnotes

  • 10 

    The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.

  • 11 

    The coordinate, z, is taken to be the perpendicular distance from the galaxy's plane.

  • 12 

    For example, our correlator setup is different to avoid some strong, persistent interference, our observations were staggered in time, and we have made delay corrections which were not standardly available previously.

  • 13 
  • 14 

    In this paper, we use "uv" to describe the projected plane within which the antennas are distributed on the ground, and "UV" to specify the ultraviolet part of the spectrum.

  • 15 
  • 16 

    All flagging has been done manually.

  • 17 
  • 18 

    If there were antennas with large delays, then delay errors, which are a function of frequency, lead to loss of amplitude which in turn introduces non-closing errors. As a precaution, we ensured that delay corrections were carried out on data with the full spectral resolution.

  • 19 

    The exception is at the center of the entire observing band; at C band, the upper and lower halves of the band overlap slightly.

  • 20 

    The Perley–Butler 2010 model, available in the CASA routine, setjy, was used.

  • 21 

    Available in CASA 4.0.0 and higher.

  • 22 

    This is opposite to that of Murphy et al. (2011).

  • 23 

    If, for example, the constant of proportionality for the Niklas & Beck non-thermal relation is the same as for the linear one used here, then one would replace both "ones" on the right-hand side of Equation (1) with SFR0.3. This would result in αNT = −1.14, or a change of 9% using SFR = 0.5 M yr−1. The SFR, in turn, would change by 3%.

  • 24 

    The "post-Basic Calibrated Data" (post-BCD) were used.

  • 25 

    Measured from the primary-beam-corrected image. The error is measured as $\Delta \,S = {\rm rms} * \sqrt{N_B}$, where NB is the number of independent beams in the region.

  • 26 

    This feature is just below point source #1 (Figure 6 and Table 6).

  • 27 

    By "gap," we do not mean the absence of Hα emission but rather a sharp drop in the emission compared to adjacent locations.

  • 28 

    This is not a result of side lobes or negative bowls in the emission.

  • 29 

    Strip #3 has been avoided since it contains CHANG-ES A.

  • 30 

    This southern peak position is not as far from the core as is the northern peak position noted above. Since the southern emission is blended with the disk of UGC 10288, it is not clear whether the southern peak position is in the jet or the lobe.

  • 31 

    No molecular gas measurements have yet been made for UGC 10288.

Please wait… references are loading.
10.1088/0004-6256/146/6/164