EDP Sciences
Free Access
Issue
A&A
Volume 603, July 2017
Article Number A32
Number of page(s) 5
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201730616
Published online 04 July 2017

© ESO, 2017

1. Introduction

Narrow-line Seyfert 1 galaxies (NLS1s) are a puzzling class of active galactic nuclei (AGN), which were first classified by Osterbrock & Pogge (1985) according to their full width at half maximum (FWHM) of Hβ< 2000 km s-1. However, despite the narrowness of Hβ, their ratio [O iii]λ5007/Hβ< 3 and the presence of strong Fe ii multiplets in the optical and UV spectrum indicate that these objects are type 1 AGN.

Radio-quiet1 NLS1s (RQNLS1s) constitute 93% of the total population up to redshift 0.8 (Komossa et al. 2006) and 96.5% at z< 0.35 (Cracco et al. 2016). Radio-loud NLS1s (RLNLS1s) are relatively uncommon. They can be divided into two different classes according to their radio spectrum in the cm range. Flat-spectrum RLNLS1s (F-NLS1s) probably have a relativistic jet pointed toward Earth and can produce γ-rays (Abdo et al. 2009a,b), while steep-spectrum RLNLS1s (S-NLS1s) often show an extended radio morphology and are likely misaligned F-NLS1s.

One of the most interesting possibilities concerning the nature of NLS1s is that they are young and evolving objects (Mathur 2000). In particular, this appears to be true for RLNLS1s: F-NLS1s might be young flat-spectrum radio quasars (FSRQs) with a small black hole mass and S-NLS1s young radio galaxies (Foschini et al. 2015; Berton et al. 2016a). However, a preference for low inclination might also play a role (e.g., Shen & Ho 2014; Peterson 2011). Thus NLS1s are a somewhat heterogeneous group.

S-NLS1s have often been associated with compact steep-spectrum objects (CSS; Oshlack et al. 2001; Komossa et al. 2006; Gallo et al. 2006; Yuan et al. 2008; Caccianiga et al. 2014, 2017; Gu et al. 2015; Schulz et al. 2015; Berton et al. 2016a), which are usually believed to be young and growing radio galaxies (Fanti et al. 1995). Only a handful of S-NLS1s were investigated in radio (Whalen et al. 2006; Antón et al. 2008; Doi et al. 2012, 2015; Richards & Lister 2015; Gu et al. 2015; Caccianiga et al. 2017). RLNLS1s indeed have a lower observed jet power than FSRQs (Foschini et al. 2015) because of their low black hole mass (Heinz & Sunyaev 2003; Foschini 2014). Therefore, while F-NLS1s are relatively easy to find because their luminosity is enhanced by relativistic beaming, S-NLS1s are not as easily detectable.

To study the radio properties of NLS1s, we carried out a survey with the Karl G. Jansky Very Large Array (JVLA) at 5 GHz in A configuration. Our sample consists of 60 sources drawn from the papers by Foschini et al. (2015) and Berton et al. (2015), and it contains radio-quiet (but not radio-silent) NLS1s, F-NLS1s, and S-NLS1s. In this paper we report the detection of extended emission in one S-NLS1s, Mrk 783. This source is one of the few NLS1 showing such an extended emission at z< 0.1. In Sect. 2 we describe the source according to results published in the literature, in Sect. 3 we describe the data reduction we performed, in Sect. 4 we present our results, in Sect. 5 we discuss them and, finally, in Sect. 6 we provide a brief summary. Throughout this work, we adopt a standard ΛCDM cosmology, with a Hubble constant H0 = 70 km s-1 Mpc-1, and ΩΛ = 0.73 (Komatsu et al. 2011). Spectral indexes are specified with flux density Sννα at frequency ν.

2. Mrk 783

Mrk 783 (RA = 13h02m58.8s Dec = + 16d24m27s) is a NLS1 galaxy first classified by Osterbrock & Pogge (1985) at z = 0.0672 (Hewitt & Burbidge 1991) with a bolometric luminosity of the AGN Lbol = 3.3 × 1044 erg s-1 (Berton et al. 2015). Its host galaxy was classified as a lenticular galaxy (Petrosian et al. 2007), but the SDSS image clearly shows the presence of a tidal tail, or a spiral arm, extended in the east direction.

The mass of the central black hole inferred from the Hβ broad component line width is about 4.3 × 107M (Berton et al. 2015). Hβ shows a prominent red wing in the broad component, indicating a receding outflow with a velocity of ~500 km s-1. This broad component is clearly visible in all the permitted lines of the optical spectrum. Conversely, narrow lines, and particularly [O iii]λ5007, do not show any outflowing component and are well reproduced by a single Gaussian profile (Berton et al. 2016b).

Mrk 783 is a strong X-ray emitter that has been detected by ROSAT (Schwope et al. 2000), INTEGRAL (Krivonos et al. 2007), and Swift/XRT (Panessa et al. 2011). Panessa et al. (2011) reported a luminosity of 9.33 × 1043 erg s-1 between 20 and 100 keV and a photon index of 1.7 ± 0.2 between 0.3 and 100 keV. This is consistent with nonsaturated comptonization, which occurs in the accretion disk corona and not in relativistic jets.

In the last 30 yrs, the galaxy was observed several times in several radio bands, for example, the WSRT at 1.4 GHz (Meurs & Wilson 1981), VLA at 5 GHz (Ulvestad & Wilson 1984; Ulvestad et al. 1995), and Green Bank telescope at 1.4 GHz (Bicay et al. 1995). However, no extended emission was found. Recently, Doi et al. (2013) observed the galaxy nucleus with the Very Long Baseline Array (VLBA) looking for extended emission near the core of the AGN. The image only shows a compact core, but the flux density recovered by the authors at 1.7 GHz is only 4% of the NRAO VLA Sky Survey (NVSS) flux density at 1.4 GHz (Sν = 33.2 mJy; Condon et al. 1998). This discrepancy means that the vast majority of the flux emitted by the galaxy is distributed in structures with relatively low brightness temperature, which could not be seen by the instrument. Another hint of the extended emission can be found in the FIRST image of the galaxy (Becker et al. 1995). The source is elongated along position angle (PA) 131° and shows a peak and a total flux density of 18.5 mJy and 28.72 mJy, respectively. At low frequencies, the TIFR Giant Metrewave Radio Telescope Sky Survey (TGSS; Intema et al. 2017) at 147 MHz reports a flux density of 89.2 ± 10.9 mJy.

Mrk 783 was classified as moderately radio-loud (Berton et al. 2015) or radio-quietm (Doi et al. 2013). The R parameter is indeed close to 10. Therefore, a different estimate of the optical magnitude or optical variability in the source could have provided two different classifications. This is not uncommon, as has been clearly shown by Ho & Peng (2001) and Kharb et al. (2014). However the radio emission does not appear to be dominant over the optical magnitude as in classical radio galaxies.

3. Data reduction

The galaxy was observed on 2015 September 6 with the JVLA at 5 GHz in A configuration with a bandwidth of 2 GHz, for a total exposure time of 10 min. We reduced and analyzed the data using the Common Astronomy Software Applications (CASA) version 4.5, the standard Expanded VLA (EVLA) data reduction pipeline, and the Astronomical Image Processing System (AIPS). The main calibrator was 3C 286. We split off the measurement set of the object from the main dataset, averaging over the 64 channels of each spectral window. After that, the object was cleaned using all the spectral windows and a natural weighting to create a first image. To improve the quality of the final map, we proceeded with iterative cycles of phase only self-calibration of the visibilities. The results of the CASA self-calibration were not satisfactory because this self-calibration caused a general increase of the noise level of the maps. Therefore we tried to redo the self-calibration process using AIPS, which, indeed, significantly improved the quality of the final images. Once the presence of extended emission was confirmed, we returned to CASA and we proceeded with a second cleaning of the data to obtain the final images. In addition to the high resolution image we produced another image, using a taper of 50 kλ, to recover the highest possible fraction of the extended emission flux density. Figure 1 shows the maps of the object before and after the application of the taper. The lower panels show the radio contours superimposed on a Sloan Digital Sky Survey (SDSS) optical image.

thumbnail Fig. 1

Top: VLA images with natural weighting of Mrk 783 at 5 GHz. The scale is 1.3 kpc arcsec-1. Left: the high resolution image is shown; the beam size is 0.45′′ × 0.40′′. The contour levels are separated by a factor starting from 3σ with σ = 0.011 mJy. The dashed contours represent a − 2σ value. Right: image of the source with a Gaussian taper of 50 kλ is shown; the beam size is 2.00′′ × 2.33′′. The contour levels are separated by a factor starting from 3σ with σ = 0.039 mJy. The dashed contours represent a − 2σ value. Bottom: the same contours superimposed on a SDSS image of the galaxy in the g band.

Open with DEXTER

4. Results

The images in Fig. 1 clearly show a compact core and extended emission at PA = 131°, which is in agreement with the FIRST data. In the high resolution image the extended emission is observed only in the southeast region of the galaxy up to a projected distance of ~8 kpc. In the tapered image, instead, we observe it on both sides of the nucleus up to a projected distance of 14 kpc in the southeast and 12.5 kpc in the northwest direction.

We fitted the central component in the high resolution image with a 2D Gaussian, using the fitting algorithm in the CASA viewer. We measured a peak flux density of 3.32 ± 0.04 mJy beam-1, an integrated flux density Sν,c = 4.03 ± 0.08 mJy, and the size of deconvolved core axes 216 ± 21 (~280) and 173 ± 24 (~225) mas (pc) for the major and minor axis, respectively. The beam size is 0.45′′ × 0.40′′. With these values we calculated the core luminosity Lc = 2.3 × 1039 erg s-1 and we used the following equation to estimate the brightness temperature (Doi et al. 2013): (1)where z is the redshift, φmaj is the core major axis, and φmin is the minor axis. This TB is only a lower limit to the peak value because our measurement includes some extended emission. The real core emission as detected with VLBI is smaller and Doi et al. (2013) quote a peak value of TB> 7.7 × 107 K at 1.7 GHz. From the tapered image we measured a total flux density Sν,t = 18.85 ± 0.03 mJy (Ltot = 1.03 × 1040 erg s-1) and we obtained the flux density of the extended emission (Sν,e = 14.82 ± 0.1 mJy, Lext = 8.5 × 1039 erg s-1), after subtracting the JVLA core flux density. To evaluate the relative importance of the core emission with respect to the total emission we measured the core dominance parameter (CD) as follows: (2)The core flux is likely an upper limit, therefore the CD parameter is also an upper limit, meaning that the galaxy is not core dominated.

To better characterize the emission we measured the spectral indexes of the two components. We measured these components by dividing the spectral windows into two bins: the first centered at 4.7 GHz and the second at 5.7 GHz. Images were made at the two frequencies, adjusting the weighting of the visibilities using the ROBUST parameter in the AIPS task IMAGR so as to obtain similar beams; they were finally convolved with identical circular beams of size 0.45 arcsec. The spectral index images were made using task COMB after blanking total intensity values below 1.5 sigma. We obtained the following spectral indexes for the core component and the southern lobe, respectively: αc ~ 0.67 ± 0.13, αext ~ 2.02 ± 0.74. The αext was estimated for the southeast lobe, avoiding the noisiest regions toward the north and east of the core.

Table 1

Fluxes of the galaxy in several radio and infrared bands.

5. Discussion

Our data revealed emission in a compact core as well as an extended component, which can be observed on both sides of the nucleus in the tapered image. The luminosity at 5 GHz is Lν = 2.1 × 1030 erg s-1 Hz-1, which is below the nominal ~7 × 1031 erg s-1 Hz-1, the Fanaroff-Riley (FR) cutoff luminosity. The CD (0.27) indicates that the emission is not core dominated and the steep spectral index of the core is different from those of highly beamed radio quasars and BL Lacs (α ~ 0), but it is more similar to that of CSS sources (α ~ 0.5). The extended component does not show any sign of an highly collimated jet, but it has a double-lobe shape that resembles the radio emission of some Seyfert galaxies, such as NGC 6764 (Hota & Saikia 2006) and Mrk 6 (Kharb et al. 2006).

Extended radio surveys of normal Seyfert galaxies and LINERs (e.g., Baum et al. 1993; Gallimore et al. 2006; Singh et al. 2015) have discovered that kiloparsec-scale radio emission is not uncommon in these objects. The emission usually seems to have an AGN origin (Gallimore et al. 2006; Singh et al. 2015), but in some cases star formation (SF) can contribute significantly (e.g., Baum et al. 1993). This could be true, indeed, in the case of NLS1 galaxies where the presence of SF (both in RQ and RL objects) has often been reported in the literature (e.g., Sani et al. 2011; Caccianiga et al. 2015). The very red mid-IR colors (W3−W4 = 2.6) measured by the Wide-field Infrared Survey Explorer (WISE) of Mrk 783 and its strong emission at 60 μm measured by the InfraRed Astronomical Satellite (IRAS) (310 mJy) seem to support this hypothesis. In order to evaluate the impact of this possible SF on the observed radio emission, we computed the parameter q22, defined as(3)where S22 μm and S1.4 GHz are the WISE 22 μm flux density and the NVSS flux density, respectively. The resulting value (q22 ~ 0.3) is significantly lower than that usually observed in SF galaxies (q22 > 1; Caccianiga et al. 2015). This means that, even if all the observed mid-IR emissions were produced by the SF, the expected radio flux at 1.4 GHz would be much lower than the observed flux (~20 per cent of the observed flux). Considering that part of the IR emission is likely due to the AGN, we conclude that, even if present, the SF alone cannot explain the majority of the observed radio emission.

We compared our spectral indexes with results in the literature. Using TGSS, NVSS, and our flux densities (Table 1) we recovered a spectral index α ~ 0.44. This is significantly flatter with respect to our measurements. A possible explanation for this discrepancy might be the presence of a break or a cutoff at 5 GHz. Furthermore, biases of the surveys toward extended emissions, for example, due to short integration time, might also cause an underestimate of the spectral index.

A possible explanation for the very steep spectral indexes observed in our VLA images might be that we are observing relic emission. This could be supported by the absence of collimated structure in our images and in the high resolution images from Doi et al. (2013). Relic emission has already been observed in some Seyfert galaxies (e.g., NGC 4235, Kharb et al. 2016) and they also show very steep spectra (e.g., Jamrozy et al. 2004; Kharb et al. 2016).

Czerny et al. (2009) found that in young radio sources with high accretion rates, radiation pressure instabilities of the accretion disk can result in intermittent activity of the radio jet. In their model, the activity phases last for 103−104 yr and they are separated by periods of 104−106 yr in which the radio jet is switched off. During the cycle, the previous activity period should manifest in the form of a relic extended emission that continues its expansion until the emitting cocoon cools down and recollapses. The maximum extension and cooling time of the cocoon depend on many factors, such as jet power and duration of the activity phase. In the case of NLS1, both quantities are considered relatively small, resulting in a very low detection rate of such emission in this class of AGN (Czerny et al. 2009; Foschini et al. 2015). Such extended emission might appear only if the central black hole mass is in the higher part of the mass distribution (which spans between 106 to 108M, e.g., Cracco et al. 2016), more precisely if MBH> 107M the BH should produce the necessary jet power to make it escape from the central regions of the AGN (Doi et al. 2012). The mass of Mrk 783 black hole is 4.3 × 107M (Berton et al. 2015), therefore it belongs to the objects that in principle could produce the extended emission. Some of the properties of the emission, such as its size and the absence of a collimated jet, might suggest that we are observing the galaxy in one of its quiescent phases not long after the switching off of the radio jet.

It is worth noting that the large scale structure of the radio emission in the tapered image has an S-like shape, which is typical of precessing radio jets (Ekers et al. 1978; Parma et al. 1985). Precessing jets could arise in binary black hole systems (e.g., Roos et al. 1993; Romero et al. 2000; Rubinur et al. 2017) or result from accretion disk instabilities (Pringle 1996; Livio & Pringle 1997). In the latter case, it might be consistent with scenario of episodic activity described previously. The resemblance to the case of Mrk 6 is strong (Kharb et al. 2006).

Another possible explanation of why we observe such steep spectral indexes might be that there is a strong interaction between a jet and interstellar medium of the galaxy. This could cause shocks and magnetic field amplification, leading to greater radiative losses and steep spectra.

This source might represent one of the few examples of the elusive parent population of F-NLS1s, i.e., S-NLS1s, Mrk 783, which lies at the edge of the RQ/RL division. The same is true for another S-NLS1 with extended emission, Mrk 1239 (Doi et al. 2015). This might indicate that more sources of this kind could be found among RQ objects. Very few RQNLS1s show detectable jets or diffuse radio emission (Berton et al. 2016b). This is partly because they are by definition faint in the radio regime, but some objects that fit the RQ definition show nonthermal core emission when observed with high sensitivity (Giroletti & Panessa 2009).

Another hypothesis regarding the parent population of F-NLS1s is based on a different assumption about NLS1s nature. Some authors believe that the narrowness of permitted lines is not due to the low black hole mass, but instead to a flat BLR observed pole-on (e.g., Decarli et al. 2008). Mrk 783, however, shows extended emission on both sides of the core, hence it likely has a non-negligible inclination. In this case then the low black hole mass estimate should not be significantly affected by any BLR flattening, although this might happen in other objects (Shen & Ho 2014).

To better understand the nature of this object, simultaneous radio observations at different frequencies are needed to study the spectral energy distribution (SED) of the emission. In particular, it is fundamental to investigate the reason of such a difference between our spectral indexes and what we found from data in the literature. Also images with an intermediate spatial resolution between our data and the data from Doi et al. (2013) (e.g., from e-MERLIN) could be useful to investigate the presence of an intermediate scale radio jet.

6. Summary

In this paper we present the first result of a survey of NLS1s carried out with the JVLA at 5 GHz in A configuration. In particular we report the detection of extended emission in the S-NLS1 Mrk 783. We found a compact core and extended emission that, in the tapered image, is observed on both the southeast and northwest sides of the galaxy nucleus, up to a maximum projected distance of 14 kpc. We excluded star formation as the dominating cause of the extended emission owing to the low value of the IR-to-radio flux ratio (q22). The latter, together with the morphology of the emission, indeed suggests an AGN origin. At the same time, in the high resolution image (Fig. 1) we could not find any sign of a large scale collimated jet and Doi et al. (2013) did not find any small scale jets. These facts, and the very steep spectral indexes that we measured, led us to hypothesize that the extended emission might be a relic and that the source might be in a quiescent period of its activity cycle, which is in agreement with the theoretical model by Czerny et al. (2009) and the young age scenario for NLS1s (Mathur 2000). Further observations are needed both in radio and in other bands to confirm our results and to investigate the nature of this object in more detail.


1

The radio loudness is defined by the parameter R, the ratio between the 5 GHz flux and the optical B-band flux (Kellermann et al. 1989). A source is considered to be radio-loud if R> 10 and radio-quiet if R< 10.

Acknowledgments

The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 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 research has made use of the NASA/IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. 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. We thank the staff of the GMRT, that made these observations possible. GMRT is run by the National Centre for Radio Astrophysics of the Tata Institute of Fundamental Research. Funding for the Sloan Digital Sky Survey has been provided by the Alfred P. Sloan Foundation and the US Department of Energy Office of Science. The SDSS web site is http://www.sdss.org. SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University.

References

All Tables

Table 1

Fluxes of the galaxy in several radio and infrared bands.

All Figures

thumbnail Fig. 1

Top: VLA images with natural weighting of Mrk 783 at 5 GHz. The scale is 1.3 kpc arcsec-1. Left: the high resolution image is shown; the beam size is 0.45′′ × 0.40′′. The contour levels are separated by a factor starting from 3σ with σ = 0.011 mJy. The dashed contours represent a − 2σ value. Right: image of the source with a Gaussian taper of 50 kλ is shown; the beam size is 2.00′′ × 2.33′′. The contour levels are separated by a factor starting from 3σ with σ = 0.039 mJy. The dashed contours represent a − 2σ value. Bottom: the same contours superimposed on a SDSS image of the galaxy in the g band.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.