Large-Scale Environment Properties of Narrow-Line Seyfert 1 Galaxies at z < 0.4

The large-scale environment is believed to affect the evolution and intrinsic properties of galaxies. It offers a new perspective on narrow-line Seyfert 1 galaxies (NLS1) which have not been extensively studied in this context before. We study a large and diverse sample of 960 NLS1 galaxies using a luminosity-density field constructed using Sloan Digital Sky Survey. We investigate how the large-scale environment is connected to the properties of NLS1 galaxies, especially their radio loudness. Furthermore, we compare the large-scale environment properties of NLS1 galaxies with other active galactic nuclei (AGN) classes, for example, other jetted AGN and broad-line Seyfert 1 (BLS1) galaxies, to shed light on their possible relations. In general NLS1 galaxies reside in less dense large-scale environments than any of our comparison samples, thus supporting their young age. The average luminosity-density and distribution to different luminosity-density regions of NLS1 sources is significantly different compared to BLS1 galaxies. This contradicts the simple orientation-based unification of NLS1 and BLS1 galaxies, and weakens the hypothesis that BLS1 galaxies are the parent population of NLS1 galaxies. The large-scale environment density also has an impact on the intrinsic properties of NLS1 galaxies; the radio loudness increases with the increasing luminosity-density. However, our results suggest that the NLS1 population is indeed heterogeneous, and that a considerable fraction of them are misclassified. We support a suggested description that the traditional classification based on the radio loudness should be replaced with the division to jetted and non-jetted sources.

The large-scale environment is believed to affect the evolution and intrinsic properties of galaxies. It offers a new perspective on narrow-line Seyfert 1 galaxies (NLS1) which have not been extensively studied in this context before. We study a large and diverse sample of 960 NLS1 galaxies using a luminosity-density field constructed using Sloan Digital Sky Survey. We investigate how the large-scale environment is connected to the properties of NLS1 galaxies, especially their radio loudness. Furthermore, we compare the large-scale environment properties of NLS1 galaxies with other active galactic nuclei (AGN) classes, for example, other jetted AGN and broad-line Seyfert 1 (BLS1) galaxies, to shed light on their possible relations. In general NLS1 galaxies reside in less dense large-scale environments than any of our comparison samples, thus supporting their young age. The average luminosity-density and distribution to different luminosity-density regions of NLS1 sources is significantly different compared to BLS1 galaxies. This contradicts the simple orientation-based unification of NLS1 and BLS1 galaxies, and weakens the hypothesis that BLS1 galaxies are the parent population of NLS1 galaxies. The largescale environment density also has an impact on the intrinsic properties of NLS1 galaxies; the radio loudness increases with the increasing luminosity-density. However, our results suggest that the NLS1 population is indeed heterogeneous, and that a considerable fraction of them are misclassified. We support a suggested description that the traditional classification based on the radio loudness should be replaced with the division to jetted and non-jetted sources.

INTRODUCTION
An intriguing class of active galactic nuclei (AGN) are narrow-line Seyfert 1 (NLS1) galaxies. They are characterized by narrow permitted emission lines (Goodrich, 1989), relatively weak [O III] (Osterbrock and Pogge, 1985) and strong Fe II (Osterbrock and Pogge, 1985). They host low or intermediate mass black holes (Peterson et al., 2000) accreting close to the Eddington limit (Boroson and Green, 1992). NLS1 sources are preferably, but not exclusively, hosted by disklike galaxies. These remarkable properties indicate that they are evolutionarily young, or possibly rejuvenated, AGN (Mathur et al., 2001).
Only ∼7% of these sources are radio-loud (RL) (Komossa et al., 2006). They usually show only a compact radio core, but recent studies suggest that pc-and kpc-scale structures are not that rare either (e.g., Gliozzi et al., 2010;Doi et al., 2013;Richards and Lister, 2015;Lister et al., 2016). Moreover, Fermi detected gamma-ray emission in some NLS1 galaxies (e.g., Abdo et al., 2009), confirming the presence of powerful relativistic jets. NLS1 population is seemingly heterogeneous including radioand gamma-loud sources as well as totally radio-silent (RS) sources. This divergence complicates the studies of the whole NLS1 population. An alternative take on this issue is to study the environments of NLS1 galaxies. The environment affects the evolution of galaxies, and consequently their nuclear activity. The environment can be divided into several scales: (1) the host galaxy; (2) the local environment, including the closest galaxies and the group or cluster the galaxy belongs to; and (3) the largescale environment, tracing the cosmic web of voids, filaments and superclusters.
Here we present the main results of a study of the large-scale environments of NLS1 galaxies. The large-scale environment position of the galaxy affects its properties; the galaxies residing in the densest environments are preferably ellipticals, and the fraction of late-type galaxies increases as the density decreases (Lietzen et al., 2012;Einasto et al., 2014;Chen et al., 2017;Kuutma et al., 2017;Pandey and Sarkar, 2017). Using a diverse and large, statistically significant sample we examine the large-scale environments of NLS1 galaxies. We compare the environment properties of NLS1 galaxies and other types of AGN, as well as study the trends within the NLS1 population. In addition, we investigate whether the large-scale environment is connected to the intrinsic properties of NLS1 galaxies. We assume a cosmology with H 0 = 73 km s −1 Mpc −1 , matter = 0.27 and vacuum = 0.73 (Spergel et al., 2007).

SAMPLE AND DATA
Our sample includes 960 sources from Zhou et al. (2006), Komossa et al. (2006), Whalen et al. (2006), Yuan et al. (2008), andFoschini (2011) at a redshift range of z = 0.0726 − 0.3996, and with a mean redshift of z = 0.2340. The radio-loudness (R, the ratio between the radio flux density, F R , and the optical flux density, F O ) was computed for sources with radio (1.4 GHz) and optical (468.6 nm) data, obtained from ASI Science Data Center (ASDC 1 ). The archival data are non-simultaneous.
K-correction was not applied since the variability of NLS1 sources in radio and optical bands affects R significantly more than the K-correction. Even though R is not an ideal way to estimate the radio characteristics of AGN, it is commonly used and enables us to compare our results with earlier studies. Accordingly, 74 of our sources are RL (10< R <100), 13 are very radio-loud (VRL, R >100) and 73 RQ (R <10). 161 sources are radio-detected (RD), leaving 799 seemingly RS sources. Some RS NLS1 galaxies have been detected at Metsähovi Radio Observatory at 37 GHz (Lähteenmäki et al. in preparation.), indicating that a fraction of them is misclassified. The arbitrary boundary between RL and RQ sources, their variability, and misclassifications might affect the data analysis. Additionally we use a subsample of 15 jetted NLS1 galaxies in which the radio 1 www.asdc.asi.it emission is known to be dominated by a relativistic jet (see section 5.3).

FWHM(Hβ)
The λL 5100 and FWHM(Hβ) values were from Zhou et al. (2006). We were able to estimate M BH for 937 sources. This method does not take into account possible inclination effects caused by the geometry of the BLR (Decarli et al., 2011), and for jetted sources it overestimates the M BH (Wu et al., 2004). Regardless, this method can be used as an order of magnitude estimate in statistical studies.

SDSS LRG LDF
SDSS LRG LDF is a three-dimensional low-resolution luminosity-density field (LDF) limited to 225-1,000 h −1 Mpc and constructed using a sample of luminous red galaxies (LRG) in SDSS DR7 (Abazajian et al., 2009). The mean luminositydensity around each source is calculated in a volume of 3 h −1 Mpc, using a 16 h −1 Mpc smoothing. Detailed description is available in Lietzen et al. (2011) and Liivamägi et al. (2012).
We study the large-scale environments of the subsamples and the total NLS1 sample and compare them to AGN samples studied in Lietzen et al. (2011) using the same LDF. Table 1 lists average luminosity-densities and the distributions into different luminosity-density regions. For the average luminosity-densities the error is the standard error of the mean and for the percentages it is the Poissonian error. Seyfert 1 (Sy1) galaxies from Lietzen et al. (2011) are also listed. They used a cut of [O III]/Hβ > 3 when selecting their Sy1 sample thus it is a pure BLS1 sample. The regions are defined as in Lietzen et al. (2011); in a void the luminosity-density is <1.0, in an intermediate luminosity-density region it is between 1.0 and 3.0, and in a supercluster >3.0. The intermediate luminosity-density regions might correspond to either physical filaments or other areas of intermediate density, e.g., the outer parts of superclusters. It should be kept in mind that whereas superclusters are the densest structures at these scales this study does not probe smaller-scale galaxy clusters which are the densest regions in the Universe.
The large-scale environment properties of NLS1 galaxies are distinct when compared to samples in Lietzen et al. (2011); they have smaller average luminosity-density than any comparison sample, and a different spatial distribution. Of the samples in Lietzen et al. (2011) NLS1 galaxies most resemble BLS1 and Sy2 galaxies; this was expected since Seyfert galaxies in general are young, late-type galaxies, whereas other AGN classes in Lietzen et al. (2011) are evolved, older galaxies, for example, BL Lac objects and Fanaroff-Riley I and II galaxies. Nevertheless, the difference is significant and supported by the divergence in their  spatial distributions; compared to BLS1 and Sy2 galaxies a larger fraction of NLS1 galaxies reside in voids and a smaller fraction in intermediate-density regions. In superclusters the percentage is almost the same. Welch's t-test for the average luminositydensities of the NLS1 and BLS1 samples gives a t-value of 3.96 which corresponds to the p-value of 0.0001, so it is highly unlikely that the averages are different by chance. Among the NLS1 subsamples the average luminosity-density seems to increase with increasing R. Statistically the differences in the averages are not so significant and suffer from the small sample sizes; for RS and RQ, RS and RL, and RQ and RL samples the p-values according to Welch's t-test are 0.98, 0.11, and 0.23, respectively. Interestingly, the average luminosity-density of RL subsample is similar to the average luminosity-density of BLS1 sources in Lietzen et al. (2011). Increasing R also seems to induce changes in the spatial distribution, but confirming this would require a larger sample.

PRINCIPAL COMPONENT ANALYSIS
Principal component analysis (PCA, e.g., Abdi and Williams, 2010) is a useful tool in AGN studies. It can be utilized to find underlying connections and the most dominant variables in a data set, and helps to identify the physical properties connected with the eigenvectors (EV). Previous PCA studies of AGN (Boroson and Green, 1992;Boroson, 2002;Grupe, 2004;Xu et al., 2012;Järvelä et al., 2015) have led to the discovery of the 4DE1 parameter space, defined by FWHM(Hβ), F([O III]) λ5007, R4570, and C IV λ1549, explaining a host of differences observed in AGN, and possibly tracing the general AGN evolution (Marziani et al., 2006;Sulentic et al., 2007).
We performed weighted PCA using the pca 2 function in MATLAB Statistics and Machine Learning Toolbox. Including a large-scale environment parameter to the PCA allows us to investigate its connection to the AGN correlation space. In PCA we used  Zhou et al. (2006). The results are presented in Table 2. The coefficients have been grouped together based on their sign, i.e., whether they correlate or anticorrelate with the EV and each other. The sign of the coefficient is insignificant as such and only their respective direction matters. , and R4570. EV1 correlates with the Eddington ratio 3 to some extent (Spearman rank correlation coefficient ρ = −0.506, and probability value p = 0.000) and weakly with the M BH (ρ = 0.291, p = 0.000).

EV3
EV3 explains 20% of the variance. Density has the strongest impact on EV3 and all the other parameters are negligible. EV3 does not correlate strongly with the Eddington ratio (ρ = -0.157, p = 0.000) or M BH (ρ = 0.100, p = 0.002). These results indicate that the large-scale environment density does not directly affect any singular properties, but rather the overall evolution of the galaxies.

NLS1s Compared to Other AGN
Earlier studies have found NLS1 galaxies to be morphologically young, late-type sources. In the framework of the densitymorphology relation (Hubble and Humason, 1931;Dressler, 1980)-more evolved galaxies are found in denser regions-our results support this and confirm the young nature of the NLS1 population.
There is a clear difference in the average luminosity-densities of the NLS1 and BLS1 sources, contradicting the simple orientation-based unification model in which they can be unified when the geometrical effects are taken into account (Decarli et al., 2008;Rakshit et al., 2017). However, the diversity of BLS1 galaxies was not taken into account in Lietzen et al. (2011) and might affect the results; for example, BLS1 galaxies with pseudobulges might be alike to NLS1 galaxies and possibly part of their parent population. Alternatively, NLS1 and BLS1 galaxies might be unified through evolution; if the narrower broad lines of NLS1 galaxies compared to BLS1 galaxies are due to an undermassive black hole, growing NLS1 galaxies would eventually evolve to have higher mass black holes and broader broad lines, thus, BLS1 galaxies.
Recently, a few higher-inclination NLS1 sources have been found (Doi et al., 2015;Congiu et al., 2017), proving that the narrow lines in these sources are due to the undermassive black hole and indicating that no additional broad-line parent population is necessarily needed. However, the parent population studies are in the early stage and include a variety of scenarios, for example, compact steep-spectrum sources and disk-hosted radio galaxies have been proposed as possible parent populations of jetted NLS1 sources (Berton et al., 2015(Berton et al., , 2016. PCA supports the idea that NLS1 sources are a part of the AGN continuum; NLS1 EVs 1 and 2 are a manifestation of the 4DE1 AGN correlation space. Shen and Ho (2014) propose that the variations in FWHM(Hβ), R4570, and [O III], and their correlations can be explained with the Eddington ratio and the orientation; increasing Eddington ratio leads to increasing R4570 and decreasing [O III], and the scatter in FWHM(Hβ) is mainly a geometrical effect. NLS1 galaxies follow the same trend as broad-line AGN and complete the high R4570, low [O III] extreme of the continuum. We can not confirm or disprove the role of the orientation, but if the FWHM(Hβ) was orientationdependent BLS1 sources should form the parent population of NLS1 galaxies, which, based on our large-scale environment results, seems improbable.

NLS1s as a Class
The incidence of RL sources increases toward the denser regions, indicating that the large-scale environment holds the capacity to alter the radio characteristics of NLS1 sources. The average M BH of the subsamples at fixed redshift are similar, as can be seen in Table 1, and R does not correlate with M BH , suggesting that M BH and the radio properties do not have a strong connection. The average M BH within the different luminosity-density regions are similar. For example, RL sources have similar average M BH in voids, intermediate-density regions, and superclusters, as well as do RQ and RS sources. This might indicate that, at least up to this point, the accretion history of these sources has been similar in all density regions.
There is some indication that the spatial distributions of subsamples change with the changing R, which would be a natural consequence of the changing average luminosity-density, but a larger, cafefully classified, sample will be needed to verify this result. Heterogeneous samples are one of the caveats in extensive AGN population studies, and a possible reason why large statistical studies might give complex, or even inconsistent, results. Our study is no exception; some faint radio sources might be incorrectly in the RS subsample, R might be too high for sources with enhanced star formation, and the variability of NLS1 sources causes R to be time-dependent. Indeed, according to the two-sample Kolmogorov-Smirnov test our subsamples are drawn from the same luminosity-density distribution, further indicating that they are, in some degree, mixed.

Special Case of Jetted NLS1s
Padovani (2016) suggests that the main difference between RQ and RL sources is the presence of a powerful relativistic jet, and that definitions RQ and RL should be replaced with non-jetted and jetted. While this is true for most sources the situation is more complicated, for example, we know that radio-loudness can be induced by enhanced star formation (Caccianiga et al., 2015) and some RS sources do host relativistic jets (Lähteenmäki et al. in preparation.). Instead of using radio-loudness as a proxy for the nuclear activity the real physical properties, specifically the presence of a jet, should be used, but this is often impossible due to the lack of data.
We can begin to address this issue by concentrating on the sources known to have jets. We implement this by selecting a sample of NLS1 galaxies whose radio emission is dominated by the AGN; gamma-detected NLS1 galaxies, NLS1 sources with extended jets, sources detected at Metsähovi (Lähteenmäki et al. (2017) and Lähteenmäki et al. in preparation), and sources in which the radio emission is jet dominated based on the q22 parameter [log(F 22 µm /F 1.4 GHz ), Caccianiga et al., 2015].
After excluding two sources from the Metsähovi sample, initially selected for the observing programme due to their extraordinarily dense environments, we have 15 almost certainly jetted sources. For these sources the average luminosity-density is 1.75 ± 0.31 (N = 15). The average luminosity-densities for the mixed and star formation (SF) dominated subsamples (based on q22) are 1.68 ± 0.15 (N = 82) and 1.44 ± 0.14 (N = 57), respectively.
Despite the rather small subsamples it seems clear that the jetted NLS1 galaxies tend to reside in denser environments. Interestingly, the average luminosity-density of these is almost similar to the traditionally defined RL NLS1 sources (1.71 ± 0.14) and the average luminosity-density of the SF dominated sources is similar to the RS sources.
These results are consistent with the LRG LDF results; the increasing R increases the probability that the radio emission originates from the jet, and in SF dominated sources the contribution of the probably faint AGN to the radio emission is negligible. This does not exclude the possibility that jetted sources could have enhanced SF.
It is noteworthy that the average large-scale environment density of jetted NLS1 galaxies is similar to that of BLS1 galaxies (Lietzen et al., 2011). Jetted AGN samples in Lietzen et al. (2011) include BL Lac objects (2.50 ± 0.20), flat-spectrum radio galaxies (2.60 ± 0.07), FR I radio galaxies (3.01 ± 0.07), and Frontiers in Astronomy and Space Sciences | www.frontiersin.org 4 November 2017 | Volume 4 | Article 54 FR II radio galaxies (3.20 ± 0.04), all with significantly higher average luminosity-densities than the jetted NLS1 galaxies (1.75 ± 0.31). This pronounced difference confirms that jets are not exclusive to evolved galaxies in high density regions, but it is possible to trigger a jet in diverse environments. The largescale environment affects the incidence of jets, but the triggering mechanism remains unclear.

SUMMARY
We examined the large-scale environment of NLS1 galaxies and its effect on their intrinsic properties. The environment provides an interesting new perspective on the AGN phenomena and helps us to understand how changes at the largest cosmic scales affect the smaller-scale properties of AGN and their evolution. Our main conclusions are the following.
1. Our results further establish the young, unevolved nature of NLS1 galaxies, and show that their radio properties are affected by the environment; RL and jetted NLS1 galaxies are more frequent in denser regions compared to RS, RQ and non-jetted sources. 2. As a class NLS1 galaxies are diverse and heterogeneous, and a considerable amount of them are misclassified. Our study supports the proposition to start using definitions non-jetted and jetted instead of RQ and RL (Padovani, 2016), and we further suggest that this division should be based on observed physical properties instead of radio-loudness parameter. 3. Our study suggests that NLS1 and BLS1 galaxies can not be unified only based on orientation. It is improbable that BLS1 galaxies would constitute the parent population of NLS1 galaxies. However, an evolutionarily link might exist. 4. NLS1 sources are part of the AGN 4DE1 continuum, but clearly distinctive in their nature when compared to the other gamma-ray emitting AGN.
The large-scale environment has an impact on the properties and evolution of AGN, amongst them NLS1 galaxies. However, it likely affects at very long time-scales and thus its connection to specific characteristics of variable and intermittent nuclear activity is weak. To effectively exploit this novel perspective the NLS1 environment studies should be extended to concern smaller-scale environment; the group and cluster scales, and host galaxy properties, both of which are likely to induce changes in AGN properties at shorter time-scales. The local environment studies (e.g., Dultzin-Hacyan et al., 1999;Krongold et al., 2001;Koulouridis et al., 2006Koulouridis et al., , 2013Koulouridis, 2014;Villarroel and Korn, 2014) indicate that Sy1 and Sy2 galaxies differ in their smaller-scale environments and host galaxies, but comprehensive studies including NLS1 galaxies are still needed. Extensive studies at various scales will clarify how the environment affects the AGN evolution, and help us establish the big picture of AGN. In addition, multiwavelength observations will be crucial to correctly classify the various NLS1 sources. This is essential to properly perform population-wide studies and address the heterogeneity issue of the NLS1 class. .

AUTHOR CONTRIBUTIONS
All of the authors mentioned below fulfil the criteria for an authorship. EJ: data analysis and interpretation, main writer, final acceptance. AL: original research idea, writing and revision of the draft, final acceptance. HL: construction of the luminositydensity field, revision of the draft, final acceptance.

FUNDING
HL is funded by PUT1627 grant from Estonian Research Council.