---
author:
- 'E. Congiu [^1]'
- 'M. Berton[^2]'
- 'M. Giroletti'
- 'R. Antonucci'
- 'A. Caccianiga'
- 'P. Kharb'
- 'M. L. Lister'
- |
    \
    L. Foschini
- 'S. Ciroi'
- 'V. Cracco'
- 'M. Frezzato'
- 'E. Järvelä'
- 'G. La Mura'
- 'J. L. Richards'
- 'P. Rafanelli'
bibliography:
- './biblio.bib'
title: |
    Kiloparsec-scale emission in the\
    narrow-line Seyfert 1 galaxy Mrk 783
---

Introduction
============

Narrow-line Seyfert 1 galaxies (NLS1s) are a puzzling class of active galactic nuclei (AGN), which were first classified by @Osterbrock85 according to their full width at half maximum (FWHM) of ${\rm H\beta }< 2000\,{\rm km\,s^{-1}}$. However, despite the narrowness of ${\rm H\beta }$, their ratio \[\]$\lambda5007/{\rm H\beta }< 3$ and the presence of strong multiplets in the optical and UV spectrum indicate that these objects are type 1 AGN.

Radio-quiet[^3] NLS1s (RQNLS1s) constitute $93\%$ of the total population up to redshift $0.8$ [@Komossa06] and $96.5\%$ at $\rm z<0.35$ [@Cracco16]. 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 $\gamma$-rays [@Abdo09a; @Abdo09b], 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 [@Mathur00]. 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 [@Foschini15; @Berton16c]. However, a preference for low inclination might also play a role [e.g., @Shen14; @Peterson11]. Thus NLS1s are a somewhat heterogeneous group.

S-NLS1s have often been associated with compact steep-spectrum objects [CSS; @Oshlack01; @Komossa06; @Gallo06a; @Yuan08; @Caccianiga14; @Gu15; @Schulz15; @Berton16c; @Caccianiga17], which are usually believed to be young and growing radio galaxies [@Fanti95]. Only a handful of S-NLS1s were investigated in radio [@Whalen06; @Anton08; @Doi12; @Richards15; @Doi15; @Gu15; @Caccianiga17]. RLNLS1s indeed have a lower observed jet power than FSRQs [@Foschini15] because of their low black hole mass [@Heinz03; @Foschini14]. 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 @Foschini15 and @Berton15a, 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 $\rm z<0.1$. In Sect.\[sec:mrk783\] we describe the source according to results published in the literature, in Sect.\[sec:datared\] we describe the data reduction we performed, in Sect.\[sec:results\] we present our results, in Sect.\[sec:discussion\] we discuss them and, finally, in Sect.\[sec:summary\] we provide a brief summary. Throughout this work, we adopt a standard $\rm \Lambda CDM$ cosmology, with a Hubble constant $H_0 = 70\,{\rm km\,s^{-1}}Mpc^{-1}$, and $\Omega_\Lambda = 0.73$ [@Komatsu11]. Spectral indexes are specified with flux density $S_{\nu} \propto \nu^{-\alpha}$ at frequency $\nu$.

Mrk 783 {#sec:mrk783}
=======

Mrk783 (R.A. = $13$h $02$m $58.8$s Dec=$+16$d $24$m $27$s) is a NLS1 galaxy first classified by @Osterbrock85 at $z = 0.0672$ [@Hewitt91] with a bolometric luminosity of the AGN $L_{bol} = 3.3\times10^{44}\,{\rm erg\,s^{-1}}$ [@Berton15a]. Its host galaxy was classified as a lenticular galaxy [@Petrosian07], 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 ${\rm H\beta }$ broad component line width is about $4.3\times10^7$ M$_{\odot}$ [@Berton15a]. ${\rm H\beta }$ shows a prominent red wing in the broad component, indicating a receding outflow with a velocity of $\sim500\,{\rm km\,s^{-1}}$. This broad component is clearly visible in all the permitted lines of the optical spectrum. Conversely, narrow lines, and particularly \[\]$\lambda5007$, do not show any outflowing component and are well reproduced by a single Gaussian profile [@Berton16b].

Mrk783 is a strong X-ray emitter that has been detected by ROSAT [@Schwope00], INTEGRAL [@Krivonos07], and Swift/XRT [@Panessa11]. @Panessa11 reported a luminosity of $9.33\times10^{43}\,{\rm erg\,s^{-1}}$ between $20$ and $100$ keV and a photon index of $1.7\pm0.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 years, the galaxy was observed several times in several radio bands, for example, the WSRT at $1.4$ GHz [@Meurs81], VLA at $5$ GHz [@Ulvestad84; @Ulvestad95], and Green Bank telescope at $1.4$ GHz [@Bicay95]. However, no extended emission was found. Recently, @Doi13 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_{\nu}=33.2\,$mJy; @Condon98]. 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 [@Becker95]. The source is elongated along position angle (PA) $\ang{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; @Intema17] at $147$ MHz reports a flux density of $89.2\pm10.9$ mJy.

Mrk783 was classified as moderately radio-loud [@Berton15a] or radio-quietm [@Doi13]. 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 @Ho01 and @Kharb14. However the radio emission does not appear to be dominant over the optical magnitude as in classical radio galaxies.

Data reduction {#sec:datared}
==============

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$ minutes. 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$\lambda$, to recover the highest possible fraction of the extended emission flux density. Fig.\[fig:map\] 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.

![image](Mrk783_MG.eps){width="45.00000%"} ![image](Mrk783_tape.eps){width="45.00000%"}\
-1.5cm ![image](SDSS_radio.pdf){width="32.00000%"} ![image](SDSS_radio_tape.pdf){width="32.00000%"}\

Results {#sec:results}
=======

The images in Fig.\[fig:map\] clearly show a compact core and extended emission at PA$=\ang{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 $\sim 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\pm0.04$ mJy beam$^{-1}$, an integrated flux density $S_{\nu,c}=4.03\pm0.08$ mJy, and the size of deconvolved core axes $216\pm21$ ($\sim 280$) and $173\pm24$ ($\sim225$) mas (pc) for the major and minor axis, respectively. The beam size is $0.45"\times0.40"$. With these values we calculated the core luminosity $L_{c}=2.3\times 10^{39}\,{\rm erg\,s^{-1}}$ and we used the following equation to estimate the brightness temperature [@Doi13]: $$\label{eq:Tb}
T_B=1.8\times 10^9 (1+z)\frac{S_{\nu,c}}{\nu^2\phi_{maj}\phi_{min}}\sim7600 {\rm\, K},$$ where $z$ is the redshift, $\phi_{maj}$ is the core major axis, and $\phi_{min}$ is the minor axis. This $T_B$ 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 @Doi13 quote a peak value of $T_B>7.7\times 10^7$ K at $1.7$ GHz. From the tapered image we measured a total flux density $S_{\nu,t} = 18.85 \pm 0.03$ mJy ($L_{tot}=1.03\times10^{40}\,{\rm erg\,s^{-1}}$) and we obtained the flux density of the extended emission ($S_{\nu,e} = 14.82\pm0.1$ mJy, $L_{ext}=8.5\times10^{39}\,{\rm 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: $$\label{eq:cd}
{\rm CD}=S_{\nu,c}/(S_{\nu,t}-S_{\nu,c})\sim0.27.$$ 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: $\alpha_{c} \sim 0.67\pm0.13$, $\alpha_{ext} \sim 2.02\pm0.74$. The $\alpha_{ext}$ was estimated for the southeast lobe, avoiding the noisiest regions toward the north and east of the core.

  Band                 Flux density (mJy)    Notes     Reference
  ------------------- -------------------- ---------- -----------
  $150$ MHz                  $89.2$          total         1
  $1.4$ GHz                  $33.2$          total         2
  $1.4$ GHz                  $28.7$          total         3
  $5$ GHz                    $4.0$            core         4
  $5$ GHz                    $14.8$         extended       4
  $5$ GHz                    $18.9$          total         4
  $60{\,\rm \mu m}$         $310.0$          total         5
  $22{\,\rm \mu m}$          $65.9$          total         6

  : Fluxes of the galaxy in several radio and infrared bands.[]{data-label="tab:flux"}

Discussion {#sec: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_{\nu}=2.1\times10^{30}\,{\rm {\rm erg\,s^{-1}}Hz^{-1}}$, which is below the nominal $\sim 7\times10^{31}\,{\rm {\rm erg\,s^{-1}}Hz^{-1}}$, the Fanaroff-Riley (FR) cutoff luminosity. The CD ($\leq 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 ($\alpha \sim 0$), but it is more similar to that of CSS sources ($\alpha \sim 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 NGC6764 [@Hota06] and Mrk6 [@Kharb06].

Extended radio surveys of normal Seyfert galaxies and LINERs [e.g., @Baum93; @Gallimore06; @Singh15] have discovered that kiloparsec-scale radio emission is not uncommon in these objects. The emission usually seems to have an AGN origin [@Gallimore06; @Singh15], but in some cases star formation (SF) can contribute significantly [e.g., @Baum93]. 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., @Sani11; @Caccianiga15]. 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\,{\rm\mu 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$$\label{eq:q22}
q22 = \log (S_{22\,{\rm \mu m}}/S_{1.4\,{\rm GHz}}),$$ where $S_{22\,{\rm \mu m}}$ and $S_{1.4\,{\rm GHz}}$ are the WISE $22\,{\rm \mu m}$ flux density and the NVSS flux density, respectively. The resulting value ($q22\sim0.3$) is significantly lower than that usually observed in SF galaxies [$q22>1$; @Caccianiga15]. 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 ($\sim20$ 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 (Tab.\[tab:flux\]) we recovered a spectral index $\alpha\sim 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 @Doi13. Relic emission has already been observed in some Seyfert galaxies [e.g., NGC 4235, @Kharb16] and they also show very steep spectra [e.g., @Jamrozy04; @Kharb16].

@Czerny09 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 $10^3-10^4$ yr and they are separated by periods of $10^4-10^6$ 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 [@Czerny09; @Foschini15]. Such extended emission might appear only if the central black hole mass is in the higher part of the mass distribution [which spans between $10^6$ to $10^8$ M$_{\odot}$, e.g., @Cracco16], more precisely if M$_{BH}>10^7$ M$_{\odot}$ the BH should produce the necessary jet power to make it escape from the central regions of the AGN [@Doi12]. The mass of Mrk 783 black hole is $4.3\times10^7$ M$_{\odot}$ [@Berton15a], 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 [@Ekers78; @Parma85]. Precessing jets could arise in binary black hole systems [e.g., @Roos93; @Romero00; @Rubinur17] or result from accretion disk instabilities [@Pringle96; @Livio97]. 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 [@Kharb06].

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 [@Doi15]. 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 [@Berton16b]. 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 [@Giroletti09].

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., @Decarli08]. 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 [@Shen14].

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 @Doi13 (e.g., from e-MERLIN) could be useful to investigate the presence of an intermediate scale radio jet.

Summary {#sec: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 Mrk783. 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.\[fig:map\]) we could not find any sign of a large scale collimated jet and @Doi13 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 @Czerny09 and the young age scenario for NLS1s [@Mathur00]. 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.

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 U.S. 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.

[^1]: enrico.congiu@phd.unipd.it

[^2]: marco.berton@unipd.it

[^3]: The radio loudness is defined by the parameter R, the ratio between the $5$ GHz flux and the optical B-band flux [@Kellermann89]. A source is considered to be radio-loud if $\rm R>10$ and radio-quiet if $\rm R <10$.
