Free Access
Issue
A&A
Volume 619, November 2018
Article Number A69
Number of page(s) 18
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201832876
Published online 08 November 2018

© ESO 2018

1. Introduction

The narrow-line Seyfert 1 galaxies (NLS1; Osterbrock & Pogge 1985) is a subclass of active galactic nuclei (AGN), characterised by the narrow permitted and forbidden emission lines (FWHM(Hβ) < 2000 km s−1, Goodrich 1989). In addition, [O III] λ5007 emission is relatively weak (F([O III])/F(Hβ) < 3, Osterbrock & Pogge 1985) and Fe II emission relatively strong in the majority of NLS1 sources (Osterbrock & Pogge 1985). Several NLS1 galaxies have been detected at gamma-rays (e.g. Abdo et al. 2009; Paliya et al. 2018), confirming the presence of powerful relativistic jets in them.

In contradiction to the usual view that only supermassive black holes (MBH >  108M) are able to launch powerful relativistic jets, (Laor 2000) the black holes in NLS1 galaxies have low or intermediate masses (MBH <  108M, Peterson et al. 2000) and accrete at extraordinarily high rates, from 0.1 Eddington ratio to even super-Eddington accretion (Boroson & Green 1992). However, it has been suggested that the unusually low black hole masses, when estimated using the emission line widths, are due to the orientation effects caused by a flattened broad-line region (BLR) seen face-on (Decarli et al. 2008; Rakshit et al. 2017). This would affect in particular the gamma-ray detected sources, known to have small inclinations. Alternative methods to estimate the black hole masses, for example, the MBH − Lbulge relation (Ryan et al. 2007; D’Ammando et al. 2018), spectropolarimetry (Baldi et al. 2016), and accretion disk modelling (Calderone et al. 2013), result in higher masses. On the other hand, mass estimates based on line dispersion, which is less affected by the orientation (Foschini et al. 2015; Berton et al. 2018), and reverberation mapping (Rafter et al. 2013; Wang et al. 2016) favour low black hole masses.

Only ∼10% of NLS1 sources are radio-loud (radio-loudness1R >  10; Komossa et al. 2006), ∼6% radio-quiet (R <  10), and the rest have not been detected in radio at all and are thus radio-silent. Particularly at lower frequencies radio-loudness does not necessarily indicate the presence of a jet. On the other hand, some radio-silent NLS1 sources have recently been detected at 37 GHz (Lähteenmäki et al. 2018), further highlighting the ambiguity of the radio-loudness parameter. Diverse radio properties and other disparities, for example, in the black hole masses and the large-scale environments, suggest that the NLS1 population is not homogeneous. It remains unclear which intrinsic properties account for the various characteristics observed in them, and how – if at all – these seemingly different NLS1 subpopulations are connected.

The connection between an AGN and its host galaxy has already been established (e.g. Storchi-Bergmann 2008; van de Ven & Fathi 2010; Pović et al. 2012; Fabian 2012; King & Pounds 2015). The host has a direct impact on the activity of the AGN by regulating the gas supply of the central black hole. In turn, various AGN feedback mechanisms, for example, radiation pressure, jets, and winds and outflows, affect the host galaxy and steer its evolution. Feedback can be either negative or positive, for example, quenching or enhancing star formation.

The nearby environment of a galaxy could affect its properties and subsequently the nuclear activity. A dense enough environment increases the probability of galaxy – galaxy interaction, and a galaxy might undergo multiple minor and major mergers, and close encounters, that distort and change its morphology. Interaction and mergers might replenish or strip the gas reservoirs of the galaxy and cause gas infall, which in turn triggers circumnuclear star formation and feeds the black hole. A clear connection between a merger and the level of AGN activity is yet to be established. Some studies argue that both minor (Taniguchi 1999; Barth et al. 2008) and major (Urrutia et al. 2008; Ellison et al. 2011) mergers can trigger nuclear activity, whereas others do not find mergers and nuclear activity connected at all (Corbin 2000; Cisternas et al. 2011; Kocevski et al. 2012). It has even been suggested that all radio-loud AGN are mergers (Chiaberge et al. 2015), predominantly major, or that the most luminous AGN are mostly triggered by major mergers (Treister et al. 2012; Villforth et al. 2014).

Host galaxy studies of NLS1 galaxies have so far mostly concentrated on the radio-quiet and radio-silent sources. They are preferably, but not exclusively, found in disk-like galaxies. Large-scale stellar bars are frequent (65%–80%) in NLS1 galaxies (Crenshaw et al. 2003; Ohta et al. 2007), as are nuclear dust spirals (83%) of which most are grand-design (80%; Deo et al. 2006). NLS1 galaxies also show enhanced star formation in their nuclear regions (Deo et al. 2006; Sani et al. 2010). Mathur (2000) suggested that NLS1 galaxies could be sources rejuvenated by a recent merger, but since only 8%–16% of NLS1 galaxies show signs of interaction or merging (Ohta et al. 2007; 20%–30% for Seyfert galaxies in general, Schmitt et al. 2001) this scenario seems insufficient to explain the nuclear activity in all NLS1 galaxies. The paucity of interaction and mergers (e.g. Ryan et al. 2007; Ohta et al. 2007) indicates that NLS1 galaxies mainly grow and evolve through secular processes. The predominance of pseudo-bulges in NLS1 host galaxies further supports this scenario (Orban de Xivry et al. 2011). A singular ensemble of properties such as bars, grand-design dust spirals, and stellar nuclear rings, might enable efficient fueling of the central black hole and in part explain extraordinarily high Eddington ratios and enhanced star formation in NLS1 galaxies.

So far the host galaxies of only four radio-loud NLS1 sources have been investigated, all of which have also been detected at gamma-rays. The morphology of 1H0323+342 is unclear. Zhou et al. (2007) claim that it is a one-armed spiral galaxy, whereas others find results indicating it is an elliptical galaxy and disturbed by merging (Antón et al. 2008; León Tavares et al. 2014). Olguín-Iglesias et al. (2017) found FBQS J1644+2619 to reside in a barred lenticular galaxy showing signs of a recent minor merger, but D’Ammando et al. (2017) argue that the host is an elliptical galaxy. PKS 2004-447 seems to reside in a barred disk-like galaxy (Kotilainen et al. 2016), and PKS 1502+036 in an elliptical galaxy (D’Ammando et al. 2018).

Comparing the radio-loud NLS1 host galaxy, multi-frequency, and environment data with the corresponding data for radio-quiet NLS1 galaxies and, eventually, other AGN classes, will help us to address several unsolved issues. We will have an additional perspective on the heterogeneity of the NLS1 population and can estimate how significant the host galaxy properties are for the nuclear activity. This helps us to unveil the mechanism responsible for the triggering and maintaining a relativistic jet. We can also investigate the intraclass evolution of NLS1 sources, especially whether the seemingly different NLS1 galaxies are connected through evolution, or have separate evolutionary lines. A better understanding of their diversity helps us to clarify their place in the AGN unification schemes and aids in the search for the parent population. Additionally, we get insight into the future evolution of NLS1 galaxies and their connection to other AGN classes, for example, flat-spectrum radio quasars, whose scaled-down young versions NLS1 galaxies are assumed to be. NLS1 sources offer us a view to the first stages of the evolution of the most powerful AGN.

In this paper we investigate the near-infrared morphologies of nine NLS1 galaxies. The paper is organised as follows: we explain the sample selection, observations and data reduction in Sect. 2, and the data analysis methods and the results for individual sources fitted reliably in Sect. 3. Compromised fits are presented in Appendix A. We discuss the results and their implications in Sect. 5, and conclude in Sect. 6. Throughout the paper we assume a cosmology with H0 = 73 km s−1 Mpc−1, Ωmatter = 0.27 and Ωvacuum = 0.73 Spergel et al. (2007).

2. Observations and data reduction

Our original sample consisted of 13 NLS1 galaxies included in the Metsähovi NLS1 observing programme (Lähteenmäki et al. 2017, 2018), with a redshift limit of z <  0.3 to ensure good enough quality of NIR data under favourable conditions. The sample is diverse: it includes radio-silent, radio-quiet, and radio-loud sources, both flat- and steep-spectrum sources, and sources detected and not detected at 37 GHz. Our purpose was to as widely as possible probe the host galaxy characteristics of a varied ensemble of NLS1 galaxies, and how they possibly affect the formation of jets in them. Due to challenging weather conditions we were able to observe nine sources. Basic properties of the sample and the observations are listed in Table 1. Exposure times of the sources vary because we were not able to complete the whole observing programme due to the weather. Details of individual sources and their modelling are provided in Sect. 3.

Table 1.

Basic properties of the sample and J-band observations.

J-band observations were performed with the Nordic Optical Telescope (NOT) during four nights, March 3–7, 2017, using the NOTCam with high resolution imaging. The NOTCam CCD is 1024  ×  1024 pixels, and with the high resolution imaging the field of view (FOV) is 80′′ × 80′′, giving a scale of 0.078′′ px−1. The weather and seeing during the observations were variable; we report them separately for each source. The total exposure times for each source are shown in Table 1. We used 9-point dithering with frame mode, with the dithering step between 10′′ and 15′′, depending on the size of the source, and skew of 1′′. Each J-band (λcentral = 1.247 μm) observation took 900s (9  ×  10  ×  10s). In addition to the targets we observed a bright (J-band magnitude ∼11) nearby star for each source to properly model the point spread function (PSF). This is needed to properly perform the deconvolution and to get rid of the AGN contamination to allow the modelling of the host galaxy. PSF stars were carefully observed so that they have very high quality, S/N >  1000, but not saturated. To perform the photometric calibration we observed the standard field AS16 (Hunt et al. 1998).

We reduced the data using the NOTCam data reduction package2 for IRAF3. The bad pixel mask provided was used in flat-field and actual images. Morning and evening flats were taken every day and combined to provide a masterflat. Dithered images were aligned using point-like sources in the FOV, and combined to a single image. Since we have multiple exposures of each source, these images were aligned and co-added to form the final image used in the data analysis. The data reduction and image combination were performed similarly for the PSF stars and standard field images.

Table 2.

Basic properties of the PSF stars.

3. Data analysis

3.1. Image analysis

We used the 2D fitting algorithm GALFIT version 3 (Peng et al. 2010) to perform a photometric decomposition of the images. The fitting procedure was similar for all sources. A crucial part of the fitting is to correctly model the PSF and thus be able to remove AGN contamination. For each source we used a separate PSF star with very high S/N, observed right after the source to have as similar conditions, for example, seeing and airmass, as possible. After manually centring the PSF image and subtracting the mean sky value – which we estimated with IRAF – GALFIT can extract the PSF directly from the image. This is the most accurate way to model the PSF, but it does not give an analytical form for the subtracted PSF. To ensure as good a spatial resolution as possible we observed the galaxies using the high resolution camera. However, the FOV of NOTCam is so small that generally there are no stars in the same FOV as our targets. Therefore we could not test the goodness of the PSF or estimate how using a separate PSF star affects the results. In most cases the conditions stayed similar during the PSF star and the NLS1 galaxy observations, but in some cases they slightly changed. These are mentioned for individual sources in Sect. 4.

Before fitting the galaxy the sky background was estimated using IRAF. For each source the sky value was measured in several separate, isolated regions of ∼150 × 150 px. The mean value was used in the fitting. The errors were estimated by fitting with ±1σ values. For magnitudes an additional zeropoint estimation error was added. We first fitted the source only with the PSF, and checked the residual image to see whether there was emission left and need for additional components. Then, if needed, we added more components, convolved with the PSF, one at a time, and after each addition checked the residual image for leftover flux and the subcomponent images to ascertain that the model components were physically reasonable. At each step we also varied the initial fitting parameters to make sure that the parameter values the fit converges to are correct and stable. We kept all the parameters free to have as unbiased a fit as possible, except in a few cases we had to freeze the central positions of some components. We used as large a fitting region as reasonably possible (but a smaller region is shown in the images for easier inspection) and tested whether the size of the fitting region affects the results. If there were any other sources inside the fitting region, we also fitted them to achieve better results. We varied the number of components and their parameters until the best, yet reasonable, fit was achieved. The goodness of fit was determined based on the χ ν 2 $ \chi_{\nu}^2$ of the fit as well as visual inspection of the model components. To model the various components of our sources we mainly used the Sérsic profile, given as: I ( r ) = I e exp [ κ n ( ( r r e ) 1 / n 1 ) ] , $$ \begin{aligned} I(r) = I_{{\mathrm{e}}}\mathtt{exp }\Bigg [-\kappa _{n}\Bigg (\bigg (\frac{r}{r_{{\mathrm{e}}}}\bigg )^{1/{n}} -1\Bigg )\Bigg ] \mathtt , \end{aligned} $$(1)

where I(r) is the surface brightness at radius r, κn is a parameter connected to the Sérsic index, n, so that Ie is the surface brightness at the half-light radius, re (Graham & Driver 2005). Sérsic profile can be used to model a variety of different light distributions in galaxies, for example, classical and pseudo-bulges, and elliptical and spiral morphologies, merely by changing the Sérsic index n. Setting n = 4 gives the traditional de Vaucouleurs profile, widely used to model classical bulges, n = 1 represents the exponential profile, useful for modelling disks, and n = 0.5 gives a Gaussian profile. In general, smaller values of n – a core flattening faster at r <  re and the intensity dropping faster at r >  re – are associated with galaxies with disk-like morphology and pseudo-bulges, and larger values of n with elliptical galaxies and classical bulges.

For the special case when n = 1, the functional form of the Sérsic profile is for historical reasons given as I ( r ) = I 0 exp ( r r s ) , $$ \begin{aligned} I(r) = I_{0}\mathtt{exp }\bigg (\frac{r}{r_{{\mathrm{s}}}}\bigg ) \mathtt , \end{aligned} $$(2)

where I(r) is the surface brightness at radius r, I0 is the central surface brightness, and rs is the scale length of the exponential disk. Explicitly when n = 1 the scale length is related to the effective radius by re = 1.678rs.

After achieving a sufficient fit we extracted the radial surface brightness profile from the observed image, the model image, and the separate component images using IRAF task ELLIPSE, which fits concentric elliptical isophotes to a 2D image. ELLIPSE is given initial values for the central x and y coordinates, ellipticity (ϵ), and the position angle (PA). The ellipticity and the position angle of the subsequent ellipses can vary. Output from this task includes the surface brightness (μ), semimajor axis length, PA, and ϵ of each isophote. We compared the radial surface brightness profiles of the galaxy and the best fit model using ELLIPSE with similar initial values, because the PSF does not have an analytical form to be plotted. The individual components in the surface brightness profiles were not necessarily concentric with the observed and model images. The errors were estimated taking into account the most significant error sources: the error in sky value estimation, and the error in the magnitude estimations of the stars used to calculate the magnitude zeropoint.

We used an additional diagnostic whenever the fit indicated the presence of a bar. When plotting the ellipticity and the position angle of the isophotes, given by ELLIPSE, as functions of radius, a bar can be identified as a region with increasing ellipticity and a rather constant position angle (ΔPA <  20, Wozniak et al. 1995; Menéndez-Delmestre et al. 2007).

3.2. Additional parameters

Due to the ambiguity of radio-loudness, we used additional characteristics to estimate the contribution of the jet in these sources. We assumed that sources detected at 37 GHz host a powerful, probably relativistic jet since no other known phenomena, for example, supernova remnants, are able to produce radio emission at the several hundred mJy level at such a high frequency. Seven out of our nine sources, of which three were formerly classified radio-silent, have been detected at 37 GHz. Additionally, we used the q22 parameter to estimate the jet contribution (Caccianiga et al. 2015) at low radio frequencies. It is defined as the logarithm of the flux density ratio of 22 μm mid-infrared (MIR) emission and 1.4 GHz radio emission (S22 μm); q 22 = log ( s 22 μ m / s 1.4 GHz ) $ \mathrm{q}22=\rm{log}({\mathrm{s}}_{22\,\mu\text{m}}/{\mathrm{s}}_{1.4\,\text{ GHz}}) $, and can be used to distinguish between sources in which the predominant source of the low frequency radio emission is the jet and sources in which the contribution of enhanced star formation is significant. It is based on the varied MIR-to-radio flux ratio in different sources; when q22 <  −0.8 the jet is presumed to be the main source of the radio emission, and when q22 >  1 the star formation may dominate the low frequency radio emission, especially in sources which have excess MIR emission, that is, the difference between the Wide-field Infrared Survey Explorer4 (WISE) magnitudes in W312 μm- and W422 μm-bands is > 2.5. Sources in between, with −0.8 <  q22 <  1, most likely have both, the star formation and the jet, contributing to the radio emission. The 1.4 GHz properties do not necessarily correlate with the properties at higher radio frequencies, and a source can host a powerful jet even if its low radio frequency emission is not dominated by it.

We used 1.4 GHz radio data from the Very Large Array (VLA) Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey5 and MIR data from WISE to calculate q22 and W3–W4 for our sample. Four out of six sources with sufficient data have q22 >  1, and the remaining two have −0.8 <  q22 <  1. For sources not detected at 1.4 GHz we calculated q22 assuming a flux density at the detection limit of the FIRST survey; 1 mJy. Since this is an upper limit, their q22 values are lower limits. Based on this these three radio-silent sources seem to be intermediate or star-formation dominated. Thus, according to the q22 parameter, the radio emission at 1.4 GHz is not dominated by the jet in any of these nine sources, and at least four out of nine can be assumed to exhibit enhanced star formation. Radio and infrared properties are shown in Table 3.

Table 3.

Radio and infrared properties of the sample, and the black hole mass estimates.

3.3. Black hole masses

The black hole masses were estimated using the second-order moment of Hβ (Foschini et al. 2015; Chen et al. 2018), which is less sensitive to the orientation than FWHM(Hβ). We used SDSS spectrum which was corrected for redshift using, when possible, the [S II]λ6731 line as a reference, to avoid wrong corrections due to blue outliers common among radio-loud NLS1 galaxies (Berton et al. 2016). If [S II] was weak or absent, we used the [O II]λ3727 line instead. The Galactic absorption correction was done using the values by Kalberla et al. (2005). The host galaxy was subtracted from the spectra following the procedure described in La Mura et al. (2007) and Chen et al. (2018). An example of the host subtraction is shown in Fig. 1. Then the power-law continuum of the AGN was subtracted. The Fe II lines in the Hβ region were subtracted using the online template6 developed by Kovačević et al. (2010) and described in Shapovalova et al. (2012). Hβ was fitted using three Gaussian components. One was used to reproduce the narrow component and two for the broad component. Figure 2 shows an example of the Hβ fitting. The wavelength of the narrow component was fixed to the restframe wavelength of Hβ. We forced its FWHM to be smaller than or equal to the FWHM of [O III], since the latter can be often turbulent and provide an overestimate of the narrow FWHM (e.g. Berton et al. 2015). After subtracting the narrow component, we estimated the second-order moment σ of Hβ as a velocity indicator. The σ is less biased for BLR geometry and inclination Peterson (2011). To derive the BLR radius, we exploited the relation between the Hβ flux and the BLR radius derived by Greene et al. (2010). For the virial factor f we adopted the value of 3.85, calculated by Collin et al. (2006). The black hole mass estimates are presented in Table 3.

thumbnail Fig. 1.

SDSS J091313.73+365817.2. An example of the host subtraction.

thumbnail Fig. 2.

SDSS J122749.14+321458.9. An example of the Hβ fitting.

4. Individual source analysis

4.1. SDSS J091313.73+365817.2

J091313.73+365817.2 is classified as a radio-quiet NLS1 galaxy and it has not been detected at 37 GHz. Based on parameter values q22 = 1.69 and W3 − W4 = 2.55 it is expected to exhibit enhanced star formation.

Seeing during the PSF star observation was 0.7′′, but varied between ∼0.7′′ − 1′′ during the four 900s exposures of the source itself. Visual inspection of the image clearly indicates that it is a spiral galaxy, probably barred. However, we did not try to fit the spiral arms since they are very faint. The optimal fit was achieved with two Sérsic functions, for the bulge and the bar, and an exponential function for the disk. The nearby galaxy was sufficiently fitted with a Sérsic function. The parameters of the best fit are listed in Table 4. The observed, model and residual images are shown in Fig. 3, and the radial surface brightness profile in Fig. 4. We tried to obtain the spectrum of the nearby galaxy using ALFOSC at NOT but were not successful; either it is too faint or has no strong emission or absorption lines. It remains unclear whether it is a satellite galaxy of J091313.73+365817.2, or a foreground or a background galaxy. Either way, with n = 0.65 it resembles a disk galaxy.

Table 4.

Best fit parameters of SDSS J091313.73+365817.2. χ ν 2 = 1 . 12 0.01 + 0.01 $ \chi_{\nu}^2 = 1.12^{+0.01}_{-0.01} $.

thumbnail Fig. 3.

SDSS J091313.73+365817.2. The field of view is 23.4′′/44.1 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

thumbnail Fig. 4.

SDSS J091313.73+365817.2. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

The fit of the NLS1 source ended up being better without the PSF, and in fits where it was included it was ∼3 mag fainter than other components, indicating that the J-band emission of this source is strongly dominated by the host galaxy also in the nuclear region. Sérsic 1 with n = 0.69 and re = 0.75 kpc depicts the bulge of the galaxy, which is a pseudo-bulge. The parameters of Sérsic 2 fit those of a bar with n = 0.51, re = 3.23 kpc and axial ratio of 0.40. To verify the presence of a bar we plotted ellipticity and position angle against radius. The plot is shown in Fig. 5, and a region matching the description of a bar is visible between approximately 1.5 − 4′′. The effective radius of the disk is 5.89 kpc. Since the spiral arms are not modelled they are distinct in the residual image in Fig. 3; there are two prominent, axisymmetric arms, connected to the ends of the bar. The other arm seems to be brighter at smaller and the other at larger radii. The spiral arms also cause the radial surface brightness profile of the model to deviate from the profile of the galaxy at radii larger than ∼5 kpc, as seen in Fig. 4. J091313.73+365817.2 appears to be a stereotype radio-quiet NLS1, hosted by a barred spiral with a pseudo-bulge.

thumbnail Fig. 5.

SDSS J091313.73+365817.2. The ellipticity and the position angle plotted against the major axis of the isophote. Symbols and colours are explained in the plot.

4.2. SDSS J111934.01+533518.7

J111934.01+533518.7 is classified as a radio-loud NLS1, but it has not been detected at 37 GHz. Its q22 = 0.5 and W3 − W4 = 2.5, indicating an intermediate case between jet-dominated and star-formation-dominated radio emission at lower frequencies.

The PSF star originally chosen for this sources had a contaminating source in the FOV, so we used the PSF star of J091313.73+365817.2 instead. The seeing during the observation was 0.7′′. Unfortunately the source itself has only two J-band exposures. The image shows that this is clearly an interacting source. To ascertain this we obtained spectra of the sources using ALFOSC at NOT. The observation was done on December 5, 2017, using grism #4 (resolution R  ∼  360) and 1.0′′ slit. Total exposure time was 3600s (3 × 1200s), and seeing 1′′. The arc lamp used was HeNe, the flat lamp Halogen, and the standard star HD19445. Also an additional bias subtraction was made. The spectra of J091313.73+365817.2 and the nearby galaxy are shown in Fig. 6. Based on emission lines visible in the spectra of both sources, Hα and [N II], they lie at the same redshift. The NLS1 nucleus resides in the brighter, northwest, source. We experimented with various combinations of components for this interacting pair, and the most reliable fit was achieved when the NLS1 was modelled with a PSF and a Sérsic, and the companion with two Sérsic functions. The best fit parameters are shown in Table 5. Other fits yielded slightly smaller χ ν 2 $ \chi_{\nu}^2$ and less residuals, but they were discarded based on the visual inspection of the subcomponents since they clearly did not correspond to actual physical components of the galaxy. Figure 7 shows the observed, model, and residual images, and Fig. 8 the radial surface brightness profiles. The r e of the Sérsic profile of the NLS1 galaxy is 0.96 kpc and its n = 1.19, compatible with a pseudo-bulge. Sérsic 3 is centred at the companion galaxy, and based on n it may be an elliptical. Sérsic 2 is a bar-like component, but it is slightly off-centre from the companion galaxy and fits the disturbed, elongated, structure between the two galaxies.

thumbnail Fig. 6.

Spectra of SDSS J111934.01+533518.7 and its companion galaxy.

Table 5.

Best fit parameters of SDSS J111934.01+533518.7. χ ν 2 = 1 . 18 0.00 + 0.00 $ \chi_{\nu}^2 = 1.18^{+0.00}_{-0.00} $.

thumbnail Fig. 7.

SDSS J111934.01+533518.7. The field of view is 23.4′′/44.1 kpc in all images. Left panel: observed image. NLS1 nucleus resides in the northwest component. Middle panel: model image, and right panel: residual image, smoothed over 3 px.

thumbnail Fig. 8.

SDSS J111934.01+533518.7. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

The modelling of this system is likely to have been affected by its noticeably disturbed morphology. The fit left considerable residuals, which is understandable taking into account the interaction which perturbs the morphology. The tidal structure is a prominent feature – also visible in the observed image – on the northeast side of the galaxy. It is the probable cause of the model surface brightness profile deviating from the observed surface brightness profile around 8 kpc in Fig. 8.

4.3. SDSS J151020.06+554722.0

J151020.06+554722.0 was not detected in FIRST survey and was therefore previously classified as radio-silent. However, it has been frequently detected at 37 GHz, with a maximum flux density of 0.83 Jy, indicating the presence of a relativistic jet (Lähteenmäki et al. 2018). Assuming a 1.4 GHz flux density at the detection limit of FIRST, 1 mJy, the star formation tracers give q22 = 0.89 and W3 − W4 = 2.42. These are lower limits and imply that the low radio frequency emission of this source is not jet-dominated.

The FWHM of the PSF star is 1.3′′, so the seeing was not optimal during this observation. We have in total four exposures, of which two unfortunately have even worse seeing (1.7′′). In addition to J151020.06+554722.0 we also fitted a nearby source, which was nicely modelled with a PSF. The best fit for J151020.06+554722.0 was achieved with a PSF, a Sérsic and en exponential disk functions. The parameters are listed in Table 6, the observed, model, and residual images are shown in Figs. 11 and 9 shows the radial surface brightness profiles. The Sérsic function, with n = 0.32 and an axial ratio of 0.30, accounts for the bar with r e = 4.18 kpc. The ϵ and PA versus radius plot is presented in Fig. 10, and it shows a region that can be interpreted as a bar around 1 − 3′′. The effective radius of the disk is 5.31 kpc. The radial surface brightness profiles of the galaxy and the model agree well and almost no residuals are left in Fig. 11. Thus, despite the non-optimal data quality, it is safe to assume that the host galaxy of J151020.06+554722.0 is a disk-like galaxy.

Table 6.

Best fit parameters of SDSS J151020.06+554722.0. χ ν 2 = 1 . 12 0.00 + 0.01 $ \chi_{\nu}^2=1.12^{+0.01}_{-0.00} $.

thumbnail Fig. 9.

SDSS J151020.06+554722.0. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

thumbnail Fig. 10.

SDSS J151020.06+554722.0. The ellipticity and the position angle plotted against the major axis of the isophote. Symbols and colours are explained in the plot.

thumbnail Fig. 11.

SDSS J151020.06+554722.0. The field of view is 23.4′′/58.8 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

4.4. SDSS J152205.41+393441.3

J152205.41+393441.3 is a radio-quiet NLS1 galaxy, detected at 37 GHz at an unusually high level with a maximum flux density of 1.43 Jy. Based on the values q22 = 0.6 and W3 − W4 = 1.95 it does not have particularly high star formation activity, but it should not be jet-dominated either at low radio frequencies.

We have 3600s of exposure time for this source, but the seeing was rather poor, ∼1.2′′. However, it is apparent from the original image in Fig. 13 that this is an interacting source. To ascertain this we obtained the spectra of both sources with the Asiago 1.82 m Copernico telescope. The observation was done on October 16, 2017. The exposure time was 3600s, and the seeing 1′′. We used grism #4 (resolution R ∼ 486). Arc lamps used were HgCd and Neon, and the standard star was HR7596. The spectra are shown in Fig. 12. The absorption lines of the companion galaxy are at the same redshift as the emission lines of the NLS1 galaxy. The NLS1 nucleus lies in the more extended, eastern, galaxy, resembling a barred spiral in visual inspection. This is a complicated source and almost equally good fits could be achieved with many combinations of functions. In addition to χ ν 2 $ \chi_{\nu}^2$ we inspected the subcomponents of the fits carefully and compared the radial surface brightness profiles of the models to the observed one. Based on these criteria the optimal, and also realistic, fit for the NLS1 galaxy was achieved with a PSF, Sérsic functions for the bulge and the bar, and an exponential function for the disk. The companion was best modelled with one Sérsic function. The model parameters are given in Table 7. Figure 13 shows the observed, model, and residual images, and the radial surface brightness profiles are shown in Fig. 14. The model parameters confirm that J152205.41+393441.3 is a late-type galaxy. The bulge component, Sérsic 1, has n = 0.84 suggesting that it is a pseudo-bulge with an effective radius of 1.62 kpc. Sérsic 2 component clearly resembles a bar with n = 0.96, re = 6.14 kpc and an axial ratio of 0.24. To study the bar more closely we plotted ϵ/PA versus radius in Fig. 15. A region fitting the description of a bar is seen between 0.75 and 2.5′′. The effective radius of the disk is 9.24 kpc. The bulge model is slightly extended towards the companion galaxy, explaining the rather small axial ratio and PA different compared to the bar and disk components. Based on the Sérsic 3 profile parameters, n = 2.83 and re = 0.51 kpc, the companion galaxy is likely to be a very small or a dwarf elliptical galaxy. All of the fit parameters are likely to be affected to some extent by the interaction. Considerable residuals are visible in the residual image in Fig. 13, owing to the disturbed morphology of the system. Large-scale residuals with S/N = 4 and possibly resembling tidal structures are seen on both sides of the galaxy.

thumbnail Fig. 12.

Spectra of SDSS J152205.41+393441.3 and its companion galaxy.

thumbnail Fig. 13.

SDSS J152205.41+393441.3. The field of view is 31.2′′/43.5 kpc in all images. Left panel: observed image. NLS1 nucleus resides in the eastern galaxy. Middle panel: model image, and right panel: residual image, smoothed over 3 px.

Table 7.

Best fit parameters of SDSS J152205.41+393441.3. χ ν 2 = 1 . 05 0.13 + 0.21 $ \chi_{\nu}^2=1.05^{+0.21}_{-0.13} $.

thumbnail Fig. 14.

SDSS J152205.41+393441.3. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

thumbnail Fig. 15.

SDSS J152205.41+393441.3. The ellipticity and the position angle plotted against the major axis of the isophote. Symbols and colours are explained in the plot.

4.5. SDSS J161259.83+421940.3

J161259.83+421940.3 is a radio-loud flat-spectrum NLS1 galaxy, also detected at 37 GHz at a maximum flux density level of 0.64 Jy. Based on the infrared characteristics, q22 = 1.9 and W3 − W4 = 2.6, this source exhibits enhanced star formation; Caccianiga et al. (2015) estimated the SFR to be 224 M yr−1, which is very high and indicates starburst activity. However, the flat radio spectrum and the 37 GHz detection suggests that also a relativistic jet is present. The kpc-scale radio morphology of J161259.83+421940.3 remains unclear due to a nearby contaminating source in JVLA observations (Berton et al. 2018).

Based on the PSF star the seeing was 0.8′′ during this observation. The four exposures we have are about the same quality, and the total exposure time was 3600s. Fitting a PSF and a Sérsic function left residuals (Fig. 16) showing an excess of flux around the nuclear region – possibly due to the circumnuclear star formation – and a prominent brightening on the northwest side of the galaxy. This structure could be a spiral arm or a morphological disturbance due to a recent merger. We experimented with different combinations of functions to achieve as good a fit as possible: the best combination was a PSF and two Sérsic functions. Parameters of the best fit are presented in Table 8, the observed, model, and residual images of J161259.83+421940.3 in Fig. 18, and the radial surface brightness profiles in Fig. 17. Sérsic 1 accounts for the bulge of the galaxy, which, based on the Sérsic index, is a pseudo-bulge. Based on the fit parameters Sérsic 2 resembles a bar-like structure, but it is very faint as seen in its surface brightness profile. However, when plotting the ellipticity and the position angle against radius, a region with bar-like properties can be seen between approximately 1.7 − 3.5′′. According to the best fit the bar component is quite extended with r e = 8.39 kpc, and there seems not to be a disk at all. However, the bar is very faint and the possible disk can be expected to be even fainter, in which case we probably would not have been able to detect it with these observations. An alternative option is that the morphology is perturbed, as suggested by Fig. 16, and instead of an undisturbed bar the feature seen in Fig. 19 traces a disturbed, elongated component caused by interaction. Both Sérsic functions are slightly off-centred from the PSF towards the bright feature, and together account for the excess flux. The contribution of the PSF to the total J-band flux is low, and the fit is almost as good without the PSF as it is with it. This indicates that the J-band emission is mostly dominated by the host galaxy, which seems plausible taking into account its enhanced star formation.

thumbnail Fig. 16.

SDSS J161259.83+421940.3. The residuals after fitting with a PSF and one Sérsic function. The field of view is 20.3′′/72.5 kpc.

Table 8.

Best fit parameters of SDSS J161259.83+421940.3. χ ν 2 = 1 . 11 0.00 + 0.01 $ \chi_{\nu}^2=1.11^{+0.01}_{-0.00} $.

thumbnail Fig. 17.

SDSS J161259.83+421940.3. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

thumbnail Fig. 18.

SDSS J161259.83+421940.3. The field of view is 20.3′′/72.5 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

thumbnail Fig. 19.

SDSS J161259.83+421940.3. The ellipticity and the position angle plotted against the major axis of the isophote. Symbols and colours are explained in the plot.

5. Discussion

We successfully modelled the host galaxy morphologies of five NLS1 sources; the results are summarised in Table 9. We were unable to model one source, J122749.14+321458.9, and the modelling of the remaining three galaxies, J123220.11+495721.8, J125635.89+500852.4 and J133345.47+414127.7, is probably compromised; these four sources suffered from poor seeing. Combined with some of the highest redshifts of the sample, the data obtained for these sources were not of good enough quality to properly model them. Their data are presented in Appendix A.

Table 9.

Summary of the results.

All five reliably fitted sources are disk-like galaxies with black hole masses well below 108M. Four of them have pseudo-bulges while the data of the remaining source were possibly not good enough for distinguishing all galaxy components from each other. The prevalence of pseudo-bulges supports the idea that NLS1 galaxies host black holes of low or intermediate masses (e.g. Järvelä et al. 2015; Cracco et al. 2016). Bars are also commonly seen, with four out of five sources showing a component that can be perceived as a bar.

The high fraction of sources that are interacting with another galaxy (J111934.01+533518.7 and J152205.41+393441.3) or showed signs of disturbed morphology during the fitting procedure (J161259.83+421940.3) is remarkable: three out of five of the properly modelled sources seem to be disturbed spiral galaxies. In addition one of the sources with worse quality data, J133345.47+414127.7, may show signs of interaction. However, since they in any case have the characteristics of late-type galaxies, for example, pseudo-bulges, it is improbable that they would previously have undergone major mergers. Two out of three have been detected at 37 GHz and we assume them to host powerful, most probably relativistic jets, and the third is a radio-loud NLS1 galaxy. This may indicate that morphological disturbances caused by galaxy – galaxy interaction indeed are connected to higher levels of nuclear activity and subsequently the launching of relativistic jets in NLS1 galaxies. If the scenario proposed by Mathur et al. (2001) that NLS1 galaxies are sources rejuvenated by mergers is correct, the previous mergers must have been minor for the spiral-like morphology of the host galaxy to have been maintained.

Compared to radio-quiet and radio-silent NLS1 galaxies the sources in our sample seem to be quite similar: late-type galaxies with pseudo-bulges and with a high fraction of bars. The divisive factor is the frequency of interaction. Among radio-quiet and radio-silent NLS1 sources only 8%–16% show signs of interaction (Ohta et al. 2007), whereas among NLS1 galaxies with probable jets, when including the previously studied sources and excluding the four in our sample with the poor quality data, four (or six if assuming that FBQS J1644+2619 and PKS 2004-447 have gone through recent minor mergers) out of eight show disturbed morphology. The sample size is still small, but assuming that 16% of NLS1 galaxies are interacting and we randomly choose eight sources, the probability that four of them are interacting is low (2.2%). Thus it seems safe to say that among NLS1 galaxies the presence of a jet is most likely connected to interaction, but since there is also proof of isolated NLS1 galaxies with jets, for example, J151020.06+554722.0 in this study, it is not the only factor that influences nuclear activity.

Our results highlight the disparity between NLS1 galaxies with and without jets (assuming that 37 GHz emission is an indicator of relativistic jets), and further support their heterogeneity. If interaction, mostly minor mergers, indeed proves to be one of the dividing factors, this suggests that the observed heterogeneity of the NLS1 population might be due to their development over time. This also explains and supports the findings in Järvelä et al. (2017) where NLS1 galaxies with jets were found to favour denser large-scale environments than NLS1 galaxies without jets. The probability of interaction, which might lead to the launching of a jet, is naturally higher in regions with more galaxies. One out of the five reliably modelled sources in our sample resides in a supercluster, three in an intermediate-density region, and one in a void (data from Järvelä et al. 2017). The only source in supercluster environment, J111934.01+533518.7, shows signs of interaction. One of the three in the intermediate-density regions, J161259.83+421940.3, shows disturbed morphology. J091313.73+365817.2 seems to be an undisturbed spiral galaxy, but then again it has not been detected at 37 GHz and therefore probably does not host a jet. J151020.06+554722.0 is an intriguing source: it does not seem to be disturbed, and according to the SDSS Main Galaxy sample it is an isolated galaxy – though it might have faint satellite galaxies not visible in the SDSS – but it is one of the sources that has the highest detection rate at 37 GHz. This raises the question of what has triggered the high-frequency radio emission that we assume to originate in the jet. However, the seeing during the J151020.06+554722.0 observation was not optimal so it is possible that the signs of interaction have passed unnoticed. The source in a void, J152205.41+393441.3, is clearly a component in an ongoing merger, suggesting that it has a smaller companion galaxy it is now merging with.

We find that q22 is not a good measure of the presence of a jet. In fact the two sources classified as exhibitingstarburst activity possess jets based on 37 GHz observations and also show disturbed morphology. This parameter may correctly estimate the contribution of the various emission components around 1.4 GHz, but it does not seem to correlate with higher radio frequencies. It might well be that the emission is dominated by star formation processes at 1.4 GHz, especially if the galaxy exhibits (post-)merger starburst activity, but it does not mean that the source would not host a jet at the same time. Vice versa, starbursts are often triggered in a merger, as seems to be the case for jets as well. Since q22 identifies starburst sources as the most star-formation-dominated, these sources might also have a higher probability of exhibiting a jet. q22 uses VLA FIRST and WISE data which in the worst case could have been acquired decades apart, and it does not take into account the high variability of the AGN and jet emission at low radio frequencies.

It is significant that none of the host galaxies of the NLS1 sources with assumed jets in our study is an elliptical galaxy (except the dubious case of J125635.89+500852.4). This result, combined with the other remarkable characteristics of NLS1 galaxies, for example, the unusually low black hole masses found also in this study, contradicts the usual AGN paradigm that only massive elliptical galaxies are able to launch and maintain relativistic jets.

6. Conclusions

We observed nine NLS1 galaxies in near-infrared to determine the morphologies of their host galaxies. We were able to reliably model five sources of which at least three have jets based on 37 GHz observations. Our main results are, firstly, that five out of nine sources in our sample are late-type galaxies, most with pseudo-bulges and bars, and three of them host jets according to 37 GHz observations. Interaction and mergers are seen in three sources. Secondly, the abundance of interacting systems may indicate that interaction is connected to nuclear activity and the launching of jets, and could explain the heterogeneity seen in the NLS1 population. Thirdly, we conclude that as proof against the conventional view of relativistic jets being exclusively launched from massive elliptical galaxies accumulates, the current evolution and unification models need to be revised.


1

Radio-loudness, R, is defined as the ratio between radio flux density, SR, and optical flux density, SO, originally at 5 GHz and in B-band, respectively (Kellermann et al. 1989).

3

IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.

Acknowledgments

We are grateful to J. Harmanen for his help with the NOT spectral data reduction and analysis, and to FINCA Science School for obtaining the spectra. We thank R. De Propris for his help with GALFIT. Based on observations made with the Nordic Optical Telescope, operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofísica de Canarias. This publication makes use of data obtained at the Metsähovi Radio Observatory, operated by Aalto University, Finland. This work is based on observations made with the Copernico and Schmidt Telescopes of the INAF-Asiago Observatory. The National Radio Astronomy Observatory 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 publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. Funding for the Sloan Digital Sky Survey (SDSS) has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Aeronautics and Space Administration, the National Science Foundation, the U.S. Department of Energy, the Japanese Monbukagakusho, and the Max Planck Society. The SDSS Web site is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium (ARC) for the Participating Institutions. The Participating Institutions are The University of Chicago, Fermilab, the Institute for Advanced Study, the Japan Participation Group, The Johns Hopkins University, the Korean Scientist Group, Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.

References

  1. Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, ApJ, 699, 976 [NASA ADS] [CrossRef] [Google Scholar]
  2. Antón, S., Browne, I. W. A., & Marchã, M. J. 2008, A&A, 490, 583 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Baldi, R. D., Capetti, A., Robinson, A., Laor, A., & Behar, E. 2016, MNRAS, 458, L69 [NASA ADS] [CrossRef] [Google Scholar]
  4. Barth, A. J., Bentz, M. C., Greene, J. E., & Ho, L. C. 2008, ApJ, 683, L119 [NASA ADS] [CrossRef] [Google Scholar]
  5. Berton, M., Foschini, L., Ciroi, S., et al. 2015, A&A, 578, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Berton, M., Foschini, L., Ciroi, S., et al. 2016, A&A, 591, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Berton, M., Congiu, E., Järvelä, E., et al. 2018, A&A, 614, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Boroson, T. A., & Green, R. F. 1992, ApJS, 80, 109 [NASA ADS] [CrossRef] [Google Scholar]
  9. Caccianiga, A., Antón, S., Ballo, L., et al. 2015, MNRAS, 451, 1795 [NASA ADS] [CrossRef] [Google Scholar]
  10. Calderone, G., Ghisellini, G., Colpi, M., & Dotti, M. 2013, MNRAS, 431, 210 [NASA ADS] [CrossRef] [Google Scholar]
  11. Chen, S., Berton, M., La Mura, G., et al. 2018, A&A, 615, A167 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Chiaberge, M., Gilli, R., Lotz, J. M., & Norman, C. 2015, ApJ, 806, 147 [NASA ADS] [CrossRef] [Google Scholar]
  13. Cisternas, M., Jahnke, K., Inskip, K. J., et al. 2011, ApJ, 726, 57 [NASA ADS] [CrossRef] [Google Scholar]
  14. Collin, S., Kawaguchi, T., Peterson, B. M., & Vestergaard, M. 2006, A&A, 456, 75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Corbin, M. R. 2000, ApJ, 536, L73 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cracco, V., Ciroi, S., Berton, M., et al. 2016, MNRAS, 462, 1256 [NASA ADS] [CrossRef] [Google Scholar]
  17. Crenshaw, D. M., Kraemer, S. B., & Gabel, J. R. 2003, AJ, 126, 1690 [NASA ADS] [CrossRef] [Google Scholar]
  18. D’Ammando, F., Acosta-Pulido, J. A., & Capetti, A. 2017, MNRAS, 469, L11 [NASA ADS] [CrossRef] [Google Scholar]
  19. D’Ammando, F., Acosta-Pulido, J. A., Capetti, A., et al. 2018, MNRAS, 478, L66 [NASA ADS] [Google Scholar]
  20. Decarli, R., Dotti, M., Fontana, M., & Haardt, F. 2008, MNRAS, 386, L15 [NASA ADS] [Google Scholar]
  21. Deo, R. P., Crenshaw, D. M., & Kraemer, S. B. 2006, AJ, 132, 321 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Ellison, S. L., Patton, D. R., Mendel, J. T., & Scudder, J. M. 2011, MNRAS, 418, 2043 [NASA ADS] [CrossRef] [Google Scholar]
  23. Fabian, A. C. 2012, ARA&A, 50, 455 [NASA ADS] [CrossRef] [Google Scholar]
  24. Foschini, L., Berton, M., Caccianiga, A., et al. 2015, A&A, 575, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Goodrich, R. W. 1989, ApJ, 342, 224 [NASA ADS] [CrossRef] [Google Scholar]
  26. Graham, A. W., & Driver, S. P. 2005, PASA, 22, 118 [NASA ADS] [CrossRef] [Google Scholar]
  27. Greene, J. E., Hood, C. E., Barth, A. J., et al. 2010, ApJ, 723, 409 [NASA ADS] [CrossRef] [Google Scholar]
  28. Hunt, L. K., Mannucci, F., Testi, L., et al. 1998, AJ, 115, 2594 [NASA ADS] [CrossRef] [Google Scholar]
  29. Järvelä, E., Lähteenmäki, A., & León-Tavares, J. 2015, A&A, 573, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Järvelä, E., Lähteenmäki, A., Lietzen, H., et al. 2017, A&A, 606, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Kellermann, K. I., Sramek, R., Schmidt, M., Shaffer, D. B., & Green, R. 1989, AJ, 98, 1195 [NASA ADS] [CrossRef] [Google Scholar]
  33. King, A., & Pounds, K. 2015, ARA&A, 53, 115 [Google Scholar]
  34. Kocevski, D. D., Faber, S. M., Mozena, M., et al. 2012, ApJ, 744, 148 [NASA ADS] [CrossRef] [Google Scholar]
  35. Komossa, S., Voges, W., Xu, D., et al. 2006, AJ, 132, 531 [NASA ADS] [CrossRef] [Google Scholar]
  36. Kotilainen, J. K., León-Tavares, J., Olguín-Iglesias, A., et al. 2016, ApJ, 832, 157 [NASA ADS] [CrossRef] [Google Scholar]
  37. Kovačević, J., Popović, L. Č., & Dimitrijević, M. S. 2010, ApJS, 189, 15 [NASA ADS] [CrossRef] [Google Scholar]
  38. Lähteenmäki, A., Järvelä, E., Hovatta, T., et al. 2017, A&A, 603, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Lähteenmäki, A., Järvelä, E., & Ramakrishnan, V. 2018, A&A, 614, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. La Mura, G., Popović, L. Č., Ciroi, S., Rafanelli, P., & Ilić, D. 2007, ApJ, 671, 104 [NASA ADS] [CrossRef] [Google Scholar]
  41. Laor, A. 2000, ApJ, 543, L111 [NASA ADS] [CrossRef] [Google Scholar]
  42. León Tavares, J., Kotilainen, J., Chavushyan, V., et al. 2014, ApJ, 795, 58 [NASA ADS] [CrossRef] [Google Scholar]
  43. Mathur, S. 2000, MNRAS, 314, L17 [NASA ADS] [CrossRef] [Google Scholar]
  44. Mathur, S., Kuraszkiewicz, J., & Czerny, B. 2001, New Ast., 6, 321 [Google Scholar]
  45. Menéndez-Delmestre, K., Sheth, K., Schinnerer, E., Jarrett, T. H., & Scoville, N. Z. 2007, ApJ, 657, 790 [NASA ADS] [CrossRef] [Google Scholar]
  46. Ohta, K., Aoki, K., Kawaguchi, T., & Kiuchi, G. 2007, ApJS, 169, 1 [NASA ADS] [CrossRef] [Google Scholar]
  47. Olguín-Iglesias, A., Kotilainen, J. K., León Tavares, J., Chavushyan, V., & Añorve, C. 2017, MNRAS, 467, 3712 [NASA ADS] [Google Scholar]
  48. Orban de Xivry, G., Davies, R., Schartmann, M., et al. 2011, MNRAS, 417, 2721 [NASA ADS] [CrossRef] [Google Scholar]
  49. Osterbrock, D. E., & Pogge, R. W. 1985, ApJ, 297, 166 [NASA ADS] [CrossRef] [Google Scholar]
  50. Paliya, V. S., Ajello, M., Rakshit, S., et al. 2018, ApJ, 853, L2 [NASA ADS] [CrossRef] [Google Scholar]
  51. Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2010, AJ, 139, 2097 [NASA ADS] [CrossRef] [Google Scholar]
  52. Peterson, B. M. 2011, ArXiv e-prints [arXiv:1109.4181] [Google Scholar]
  53. Peterson, B. M., McHardy, I. M., Wilkes, B. J., et al. 2000, ApJ, 542, 161 [NASA ADS] [CrossRef] [Google Scholar]
  54. Pović, M., Sánchez-Portal, M., Pérez García, A. M., et al. 2012, A&A, 541, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Rafter, S. E., Kaspi, S., Chelouche, D., et al. 2013, ApJ, 773, 24 [NASA ADS] [CrossRef] [Google Scholar]
  56. Rakshit, S., Stalin, C. S., Chand, H., & Zhang, X.-G. 2017, ApJS, 229, 39 [NASA ADS] [CrossRef] [Google Scholar]
  57. Ryan, C. J., De Robertis, M. M., Virani, S., Laor, A., & Dawson, P. C. 2007, ApJ, 654, 799 [NASA ADS] [CrossRef] [Google Scholar]
  58. Sani, E., Lutz, D., Risaliti, G., et al. 2010, MNRAS, 403, 1246 [NASA ADS] [CrossRef] [Google Scholar]
  59. Schmitt, H. R., Antonucci, R. R. J., Ulvestad, J. S., et al. 2001, ApJ, 555, 663 [NASA ADS] [CrossRef] [Google Scholar]
  60. Shapovalova, A. I., Popović, L. Č., Burenkov, A. N., et al. 2012, ApJS, 202, 10 [NASA ADS] [CrossRef] [Google Scholar]
  61. Spergel, D. N., Bean, R., Doré, O., et al. 2007, ApJS, 170, 377 [NASA ADS] [CrossRef] [Google Scholar]
  62. Storchi-Bergmann, T. 2008, Rev. Mex. Astron. Astrofis., 32, 139 [Google Scholar]
  63. Taniguchi, Y. 1999, ApJ, 524, 65 [NASA ADS] [CrossRef] [Google Scholar]
  64. Treister, E., Schawinski, K., Urry, C. M., & Simmons, B. D. 2012, ApJ, 758, L39 [NASA ADS] [CrossRef] [Google Scholar]
  65. Urrutia, T., Lacy, M., & Becker, R. H. 2008, ApJ, 674, 80 [NASA ADS] [CrossRef] [Google Scholar]
  66. van de Ven, G., & Fathi, K. 2010, ApJ, 723, 767 [NASA ADS] [CrossRef] [Google Scholar]
  67. Villforth, C., Hamann, F., Rosario, D. J., et al. 2014, MNRAS, 439, 3342 [NASA ADS] [CrossRef] [Google Scholar]
  68. Wang, F., Du, P., Hu, C., et al. 2016, ApJ, 824, 149 [NASA ADS] [CrossRef] [Google Scholar]
  69. Wozniak, H., Friedli, D., Martinet, L., Martin, P., & Bratschi, P. 1995, A&AS, 111, 115 [NASA ADS] [Google Scholar]
  70. Zhou, H., Wang, T., Yuan, W., et al. 2007, ApJ, 658, L13 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A Compromised fits

The following sources did not have good enough data to properly model them. These sources suffered a combination of some of the highest redshifts among our sample and worst weather conditions during the observations. Only the radial surface brightness profile of the galaxy are presented. The parameters of the fits were not stable but varied considerably depending on the sky value, and are therefore not reliable. The host galaxy morphologies of these sources remain unclear.

A.1. SDSS J122749.14+321458.9

J122749.14+321458.9 is classified as a radio-loud flat-spectrum NLS1 galaxy. It has been detected at 37 GHz but at a rather low flux density level of around 0.2 Jy. q22 of 1.19 and W3–W4 of 2.62 indicate strong star formation; 54 M yr−1 according to Caccianiga et al. (2015). In the JVLA radio map J122749.14+321458.9 has a compact, unresolved core and possibly some diffuse radio emission on one side of the source (Berton et al. 2018).

We have two 900s exposures of this source with seeing of 1.2′′. The fit with only a PSF left prominent axisymmetric residuals, probably due to compromised PSF modelling and poor quality data. We attempted to model the source with the PSF and a variety of other functions and their combinations, but in all cases the functions converged to parameters that are not physical; n around 15–20 and the axial ratio betwen 0.01–0.05. It seems that due to the poor seeing, combined with the short exposure time, J122749.14+321458.9 cannot be properly modelled. The fit is shown only with the PSF. The parameters are given in Table A.1, the observed, model, and residual images in Fig. A.1, and the radial surface brightness profile of the galaxy in Fig. A.2.

Table A.1.

Best fit parameters of SDSS J122749.14+321458.9. χ ν 2 = 1 . 22 0.01 + 0.02 $ \chi_{\nu}^2=1.22^{+0.02}_{-0.01} $.

thumbnail Fig. A.1.

SDSS J122749.14+321458.9. The field of view is 23.4′′/54.5 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

thumbnail Fig. A.2.

SDSS J122749.14+321458.9. The observed radial surface brightness profile. Symbols and colours are explained in the plot. The shaded area around the profile describes the associated errors.

A.2. SDSS J123220.11+495721.8

J123220.11+495721.8 is a source previously classified radio-silent. It has, however, been detected at 37 GHz multiple times with a maximum flux density of 0.56 Jy. Its star formation indicators are not particularly high, q22 = 0.94 (assuming S 1.4 GHz = 1 mJy) and W3 − W4 = 2.19, however, these are lower limits and in reality they can be higher.

The exposure time for this source was 3600s and the seeing ∼1′′. The source was fitted with a PSF and a Sérsic. The parameters of the fit are listed in Table A.2, Fig. A.3 shows the observed, model, and residual images, and Fig. A.4 the radial surface brightness profile of the galaxy. The Sérsic index is suspiciously high and depends on the sky value, making the fit physically unreliable. In fact, the fit seemed to be equally good with any Sérsic index between 1 and 20, adding to the inaccuracy of the result.

Table A.2.

Best fit parameters of SDSS J123220.11+495721.8. χ ν 2 = 1 . 17 0.00 + 0.00 $ \chi_{\nu}^2=1.17^{+0.00}_{-0.00} $.

thumbnail Fig. A.3.

SDSS J123220.11+495721.8. The field of view is 23.4′′/91.3 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

thumbnail Fig. A.4.

SDSS J123220.11+495721.8. The observed radial surface brightness profile. Symbols and colours are explained in the plot. The shaded area around the profile describes the associated errors.

A.3. SDSS J125635.89+500852.4

J125635.89+500852.4 is a source previously classified as radio-silent, however, it has been detected at 37 GHz with a flux density of S37 GHz, max = 0.61 Jy. Assuming a 1.4 GHz flux density of 1 mJy its q22 = 0.61 and W3 − W4 = 2.43, therefore the jet does not seem to be dominant at low radio frequencies.

Seeing during the PSF star observation was 0.9′′, but got worse (∼1.4′′) during the source observation. Since the source lies at a considerable redshift, the bad seeing is likely to affect the results of the modelling. We have four 900s exposures of this source. We modelled the source with a PSF and a Sérsic function, however, the fit parameters were not stable. Table A.3 lists the fit parameters, the observed, model, and residual images are shown in Fig. A.5, and the radial surface brightness profile of the galaxy in Fig. A.6. The Sérsic index of the model is 5.4, suggesting that J125635.89+500852.4 is an elliptical galaxy. However, the Sérsic index as well as other model parameters considerably change with the 1σ values of the sky, casting doubt on the credibility of the model. Interestingly, the nucleus of J125635.89+500852.4 is faint compared to the galaxy and indeed, the fit is equally good without the PSF.

Table A.3.

Best fit parameters of SDSS J125635.89+500852.4. χ ν 2 = 1 . 18 0.00 + 0.01 $ \chi_{\nu}^2=1.18^{+0.01}_{-0.00} $.

thumbnail Fig. A.5.

SDSS J125635.89+500852.4. The field of view is 23.4′′/87.0 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

thumbnail Fig. A.6.

SDSS J125635.89+500852.4. The observed radial surface brightness profile. Symbols and colours are explained in the plot. The shaded area around the profile describes the associated errors.

A.4. SDSS J133345.47+414127.7

J133345.47+414127.7 is a borderline source between radio-quiet and radio-loud. It was previously classified as a flat-spectrum source, but based on JVLA observations re-classified as a steep-spectrum source with compact radio morphology (Berton et al. 2018). However, detections at 37 GHz indicate that it has a strongly variable jet. Star formation tracers of J133345.47+414127.7 are high, q22 = 1.3 and W3 − W4 = 2.9. Caccianiga et al. (2015) estimate its SFR to be as high as 148 M yr−1.

Based on the PSF star the seeing was not optimal: ∼1.1′′. We have four 900s exposures of approximately the same quality. To increase the fitting area we also fitted two nearby sources, which were sufficiently modelled with PSFs. After several attempts the best fit of J133345.47+414127.7 was achieved with two PSFs and a Sérsic function. Trying to fit the galaxy with only one PSF, and a single or multiple Sérsic functions consistently yielded results with one of the Sérsic functions having an index of ∼20 that is not realistic, however, it indicates the presence of an additional compact component. Adding another PSF solves this problem. The separation of the centres of PSFs is 2.25 kpc. This is less than what the seeing allows to resolve so the fit should be treated with caution. For this reason the best fit parameters of the fit with a single and with two PSFs are shown in Tables A.4 and A.5, respectively. The observed, model, and residual images of the fit with one PSF and a Sérsic function are presented in Fig. A.8, and the radial surface brightness profile of the galaxy in Fig. A.7. Looking at the values and errors of the fits it is clear that they are not stable and most probably not reliable. However, in the residual image in Fig. A.8, even if the fit of the nuclear region is bad, it is evident that there is non-axisymmetric flux excess on the northern side of the galaxy that could be a tidal feature.

Table A.4.

Best fit parameters of SDSS J133345.47+414127.7, one PSF. χ ν 2 = 1 . 14 0.00 + 0.02 $ \chi_{\nu}^2=1.14^{+0.02}_{-0.00} $.

Table A.5.

Best fit parameters of SDSS J133345.47+414127.7, two PSFs. χ ν 2 = 1 . 13 0.00 + 0.01 $ \chi_{\nu}^2=1.13_{-0.00}^{+0.01} $.

thumbnail Fig. A.7.

SDSS J133345.47+414127.7. The observed radial surface brightness profile. Symbols and colours are explained in the plot. The shaded area around the profile describes the associated errors.

thumbnail Fig. A.8.

SDSS J133345.47+414127.7. The field of view is 23.4′′/81.5 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image,smoothed over 3 px.

All Tables

Table 1.

Basic properties of the sample and J-band observations.

Table 2.

Basic properties of the PSF stars.

Table 3.

Radio and infrared properties of the sample, and the black hole mass estimates.

Table 4.

Best fit parameters of SDSS J091313.73+365817.2. χ ν 2 = 1 . 12 0.01 + 0.01 $ \chi_{\nu}^2 = 1.12^{+0.01}_{-0.01} $.

Table 5.

Best fit parameters of SDSS J111934.01+533518.7. χ ν 2 = 1 . 18 0.00 + 0.00 $ \chi_{\nu}^2 = 1.18^{+0.00}_{-0.00} $.

Table 6.

Best fit parameters of SDSS J151020.06+554722.0. χ ν 2 = 1 . 12 0.00 + 0.01 $ \chi_{\nu}^2=1.12^{+0.01}_{-0.00} $.

Table 7.

Best fit parameters of SDSS J152205.41+393441.3. χ ν 2 = 1 . 05 0.13 + 0.21 $ \chi_{\nu}^2=1.05^{+0.21}_{-0.13} $.

Table 8.

Best fit parameters of SDSS J161259.83+421940.3. χ ν 2 = 1 . 11 0.00 + 0.01 $ \chi_{\nu}^2=1.11^{+0.01}_{-0.00} $.

Table 9.

Summary of the results.

Table A.1.

Best fit parameters of SDSS J122749.14+321458.9. χ ν 2 = 1 . 22 0.01 + 0.02 $ \chi_{\nu}^2=1.22^{+0.02}_{-0.01} $.

Table A.2.

Best fit parameters of SDSS J123220.11+495721.8. χ ν 2 = 1 . 17 0.00 + 0.00 $ \chi_{\nu}^2=1.17^{+0.00}_{-0.00} $.

Table A.3.

Best fit parameters of SDSS J125635.89+500852.4. χ ν 2 = 1 . 18 0.00 + 0.01 $ \chi_{\nu}^2=1.18^{+0.01}_{-0.00} $.

Table A.4.

Best fit parameters of SDSS J133345.47+414127.7, one PSF. χ ν 2 = 1 . 14 0.00 + 0.02 $ \chi_{\nu}^2=1.14^{+0.02}_{-0.00} $.

Table A.5.

Best fit parameters of SDSS J133345.47+414127.7, two PSFs. χ ν 2 = 1 . 13 0.00 + 0.01 $ \chi_{\nu}^2=1.13_{-0.00}^{+0.01} $.

All Figures

thumbnail Fig. 1.

SDSS J091313.73+365817.2. An example of the host subtraction.

In the text
thumbnail Fig. 2.

SDSS J122749.14+321458.9. An example of the Hβ fitting.

In the text
thumbnail Fig. 3.

SDSS J091313.73+365817.2. The field of view is 23.4′′/44.1 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

In the text
thumbnail Fig. 4.

SDSS J091313.73+365817.2. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

In the text
thumbnail Fig. 5.

SDSS J091313.73+365817.2. The ellipticity and the position angle plotted against the major axis of the isophote. Symbols and colours are explained in the plot.

In the text
thumbnail Fig. 6.

Spectra of SDSS J111934.01+533518.7 and its companion galaxy.

In the text
thumbnail Fig. 7.

SDSS J111934.01+533518.7. The field of view is 23.4′′/44.1 kpc in all images. Left panel: observed image. NLS1 nucleus resides in the northwest component. Middle panel: model image, and right panel: residual image, smoothed over 3 px.

In the text
thumbnail Fig. 8.

SDSS J111934.01+533518.7. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

In the text
thumbnail Fig. 9.

SDSS J151020.06+554722.0. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

In the text
thumbnail Fig. 10.

SDSS J151020.06+554722.0. The ellipticity and the position angle plotted against the major axis of the isophote. Symbols and colours are explained in the plot.

In the text
thumbnail Fig. 11.

SDSS J151020.06+554722.0. The field of view is 23.4′′/58.8 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

In the text
thumbnail Fig. 12.

Spectra of SDSS J152205.41+393441.3 and its companion galaxy.

In the text
thumbnail Fig. 13.

SDSS J152205.41+393441.3. The field of view is 31.2′′/43.5 kpc in all images. Left panel: observed image. NLS1 nucleus resides in the eastern galaxy. Middle panel: model image, and right panel: residual image, smoothed over 3 px.

In the text
thumbnail Fig. 14.

SDSS J152205.41+393441.3. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

In the text
thumbnail Fig. 15.

SDSS J152205.41+393441.3. The ellipticity and the position angle plotted against the major axis of the isophote. Symbols and colours are explained in the plot.

In the text
thumbnail Fig. 16.

SDSS J161259.83+421940.3. The residuals after fitting with a PSF and one Sérsic function. The field of view is 20.3′′/72.5 kpc.

In the text
thumbnail Fig. 17.

SDSS J161259.83+421940.3. The observed and model radial surface brightness profiles. Symbols and colours are explained in the plot. The shaded area around each profile describes the associated errors.

In the text
thumbnail Fig. 18.

SDSS J161259.83+421940.3. The field of view is 20.3′′/72.5 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

In the text
thumbnail Fig. 19.

SDSS J161259.83+421940.3. The ellipticity and the position angle plotted against the major axis of the isophote. Symbols and colours are explained in the plot.

In the text
thumbnail Fig. A.1.

SDSS J122749.14+321458.9. The field of view is 23.4′′/54.5 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

In the text
thumbnail Fig. A.2.

SDSS J122749.14+321458.9. The observed radial surface brightness profile. Symbols and colours are explained in the plot. The shaded area around the profile describes the associated errors.

In the text
thumbnail Fig. A.3.

SDSS J123220.11+495721.8. The field of view is 23.4′′/91.3 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

In the text
thumbnail Fig. A.4.

SDSS J123220.11+495721.8. The observed radial surface brightness profile. Symbols and colours are explained in the plot. The shaded area around the profile describes the associated errors.

In the text
thumbnail Fig. A.5.

SDSS J125635.89+500852.4. The field of view is 23.4′′/87.0 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image, smoothed over 3 px.

In the text
thumbnail Fig. A.6.

SDSS J125635.89+500852.4. The observed radial surface brightness profile. Symbols and colours are explained in the plot. The shaded area around the profile describes the associated errors.

In the text
thumbnail Fig. A.7.

SDSS J133345.47+414127.7. The observed radial surface brightness profile. Symbols and colours are explained in the plot. The shaded area around the profile describes the associated errors.

In the text
thumbnail Fig. A.8.

SDSS J133345.47+414127.7. The field of view is 23.4′′/81.5 kpc in all images. Left panel: observed image, middle panel: model image, and right panel: residual image,smoothed over 3 px.

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.