arXiv:2010.12595v1 [astro-ph.GA] 23 Oct 2020 ★ © rdt earlal intr ftepeec fa G ngal in AGN an of presence the of ( signature ies reliable a be a to they ered potential, ionization high very their of Because eV). eoatc n pc Administration. Space t and with Aeronautics NNH14CK55B contract under Hawaii of University the by .Prieto A. MNRAS ∼ (SNRs; 2003 al. et Reunanen potent (ionization species ionized highly transitio in emission collisions structure by ionization fine forbidden high from forbidden originate (FHILs), or (CLs), lines Coronal INTRODUCTION 1 † n tr eyatvl.Mroe,i NsCssrnl vary strongly CLs SNRs time. in Moreover, be actively. to need very sources if stars the case ing emission that in CL and out, detectable ruled produce is origin to AGN necessary be might them 1999 ( 2009 ooa ieFrs G :pyia rpriso h emissi the of properties regions physical I: AGN Forest Line Coronal cetdXX eevdYY noiia omZZZ form original in YYY; Received XXX. Accepted 1 2 4 3 otshe l 2009 al. et Pottasch .C Cerqueira-Campos, C. F. iiã eAtoíia nttt ainld eqia Es Pesquisas de Nacional Instituto Astrofísica, de Divisão aoaói ainld srfsc u o sao Unido Estados dos Rua - Astrofísica de Nacional Laboratório nttt eAtoíiad aais /VaLce,SN Sa S/N, Láctea, Via C/ Canarias, de Astrofísica de Instituto 10 eatmnod srnma nttt eFsc,Universi Física, de Instituto Astronomia, de Departamento iiigAtooe tteIfae eecp aiiy wh Facility, Telescope Infrared the at Astronomer Visiting 09TeAuthors The 2019 -al [email protected] E-mail: 31 eso ta.1984 al. et Penston .TpclC uioiisi hs atrojcsaelow, are objects latter these in luminosities CL Typical ). remnants supernova in produced also are they However, ). oos tal. et Komossa − 33 000 lv ta.1999 al. et Oliva r s erg , 1 – 4 − 20 1 n .G Dahmer-Hahn G. L. and (2019) ( oiae l 2018 al. et Dopita ,adWl–ae tr ( stars Wolf–Rayet and ), ( ; 2008 ayple l 2008 al. et Satyapal ; reo&Vea 2000 Viegas & Prieto ; ), mt ta.2009 al. et Smith oos tal. et Komossa 6 on lcrndniisi h ag 3 range the narr in the spec densities of a electron conditions found physical there indeed derive that We are AGNs. conclude they non-CLiF we and if ratios unveiling emission-line and coronal objects paring these AGN most of pr in physical (NLR) found the understanding is at c what aimed analysis to are wavelength contrast AGN) which (CLiF lines, Nuclei ionization Galactic Active Forest Coronal-Line ABSTRACT tla eoiydseso hwta oto h bet ha objects the of most cen that the show matter-boun of dispersion by mass velocity stellar dominated The observed. NLR spectrum a line suggest high-ionization We AGN. the by ization u eut ml htCi G sntasprt aeoyo A words: of Key class. category AGN the separate within a distribution represe not the they of is analyzed, AGN properties CLiF emission-line infrared that imply results Our .Tu,svrltosnsof thousands several Thus, ). . 3 × 1 10 ; ★ olig&Alexander & Goulding ( 4 2009 carr&Stasinska & Schaerer ,paeaynebulae planetary ), .Rodríguez-Ardila, A. ,sgetn htteinzto ehns sascae p associated is mechanism ionization the that suggesting K, ; reoe l 2002 al. et Prieto and ) aais cie–glxe:Syet–ifae:galaxies infrared: – Seyfert galaxies: – active galaxies: aii,AeiadsAtoats15,SoJs o Campos, dos José São 1758, Astronautas dos Avenida paciais, c soperated is ich rsoa eL aua eeie Spain Tenerife, Laguna, La de Cristobal n aeFdrld i rned u,A.BnoGnavs 9500. Gonçalves, Bento Av. Sul, do Grande Rio do Federal dade a,IP ial, age al. et Wang 2 5,Bir a aõs E 70-6,Iauá G Brazil MG, Itajubá, 37504-364, CEP Nações. das Bairro 154, s eNational he sexcited ns econsid- re ≥ form- lines with rpit2 coe 00Cmie sn NA L MNRAS using Compiled 2020 October 27 Preprint 100 ax- an ; ( 0 V r h otpoietCsi h pia region optical the in CLs prominent most the are ( eV) 100 = ly [Ne ally, 2006 fojcsi eun n o rvnb estvt detectio sensitivity [Si faintest by the driven analysis, not their i and In 3 CL remaining genuine sues. the one is in objects, CLs objects of the of lack of the that 67% showed They in identified. that found (NIR) near-infrared up of scales ( to parsecs extended hundred be few can region a emission their and ar p lines energetic uous these of highly Some activity. of nuclear existence the with the associated to attributed is spectra ( ( ∼ stesrnetC ntena-nrrd(I) oee,CL However, (NIR). near-infrared the in CL to up strongest the is lythem. play oe bu he reso antd,fo 10 from magnitude, of orders three about cover years. few orge-riae l 2011 al. et Rodríguez-Ardila 2010 al. et Mazzalay 1998 Taniguchi & Murayama 2011 10 lhuhCsaefeun nANseta o l G dis- AGN all not spectra, AGN in frequent are CLs Although ; (AGN nuclei galactic active in lines coronal of presence The 38 hwdta nteeojcs hyfd ntmsae fa of timescales on fade they objects, these in that showed ) . 1 6 j azlye l 2010 al. et Mazzalay , × r s erg 55e ([S eV =505 2 v 10 † ] ie tal. et Riffel − 2 _ .Riffel, R. 1 45Å(P=10e)ad[Fe and eV) 100 = (IP Å 3435 1 - esrdi G 01a redshift a at 4051 NGC in measured , . 7 × 10 xii .Teosre uioiiso uhlines such of luminosities observed The ). ( 4 ] 2006 reoe l 2005 al. et Prieto ta xeso otelwhg ends low/high the to extension an nt cm _ 69Å aebe eotdi AGN in reported been have Å) 7609 evle nteitra 10 interval the in values ve 3 cse nteNro ieRegion Line Narrow the in ocesses ; ). ,uigasml f4 G nthe in AGN 47 of sample a using ), − r odffrne ewe CLiF between differences no are hl [Si while ) .Marinello, M. rlbakhl,drvdfo the from derived hole, black tral 3 ülrSnhze l 2011 al. et Müller-Sánchez n eprtrso 3 of temperatures and .Hr,w ar u multi- a out carry we Here, s. wln ein(L)gsand gas (NLR) region line ow aatrzdb toghigh- strong by haracterized a ls fAN ycom- By AGN. of class ial 22-1,S,Brazil SP, 12227-010, ot lge S Brazil. RS, Alegre, Porto e lust xli the explain to clouds ded N.I l optical/near- all In GNs. vi iaiywt photoion- with rimarily vi ooa uioiywas luminosity coronal ] ; 1.963 ] orge-riae al. et Rodríguez-Ardila 38 vii pt 10 to up A ` T on-line ] I 6 eV) 166 = (IP m E tl l v3.0 file style X 2 _ I 07Å(IP Å 6087 = conspic- e 41 . 7 7 0.00234. rocesses − .Usu- ). × r s erg 8 10 M is- n 3% 3 ⊙ − to 1 s - ) . 2 F. C. Cerqueira-Campos et al.
More distant sources, such as Ark 564 (I =0.0247) or PG 1448+273 2.1 Optical Spectroscopy (I =0.065) display a [Si vi] luminosity nearly 2 dex higher. All ob- Optical spectra for SDSS J124134.25+442639.2, MRK 1388, jects were observed with similar signal to noise ratio (S/N). Similar III Zw 77 and SDSS J164126.91+432121.5 were obtained from the results are also found by Lamperti et al. (2017). In the optical re- Sloan Digital Sky Survey Data Release 7 (SDSS DR7) database gion, CLs of iron (i.e. [Fe vii] and [Fe x]) are usually observed, but (Abazajian et al. 2009). These spectra have already been published no quantitative studies of its frequency have been made yet. Despite by RET15. Below, we detailed the observations and data reduction its importance, these lines tend to be weak. [Fe vii] 6087 Å, the procedure for ESO 138-G001 and NGC 424. brightest CL in the 3500 - 7500 Å region, is usually 1-10% the strength of [O iii] 5007 Å (Murayama & Taniguchi 1998). Recently, a new class of AGN, named AGNs with coronal line forest (CLiF AGNs), were introduced in the literature (Rose et al. 2.1.1 Goodman/SOAR Observations 2015a; Rose et al. 2015b; Glidden et al. 2016). Among several inter- Optical spectra of ESO 138-G001 and NGC 424 were obtained esting properties, CLiF AGN are characterized by: (8) the emission with the 4.1 m Southern Observatory for Astrophysical Research vii _ V > . 88 x _ V > line flux ratio [Fe ] 6087/H 0 25; ( ) [Fe ] 6374/H (SOAR) Telescope at Cerro Pachon, Chile. The observations were 0.2; (888) [Ne v] _3426/HV > 1; (8E) The CLs are not blueshifted carried out using the Goodman Spectrograph (Clemens et al. 2004), with respect to the low-ionization lines (Δv . 100 km s−1); (E) The equipped with a 400 l/mm grating and a 0.8 arcsec slit width, giving velocity widths of the CLs are narrow and similar to that of low- a resolution R∼1500. In addition to the science frames, standard ionization lines (FWHM < 300 km s−1) using single Gaussian fits. stars (Baldwin & Stone 1984) were observed for flux calibration. Optical luminosities of these CLs are found to be between 1040 to 41 −1 HgAr arc lamps were taken after the science frames for wavelength 10 erg s (see Sect. 4.4). These values are well within the interval calibration. Daytime calibrations include bias and flat field images. of CL luminosities reported by Gelbord et al. (2009) in a sample 63 The data were reduced using standard IRAF tasks. It includes AGN. subtraction of the bias level and division of the science and flux stan- Another peculiar characteristic identified by Rose et al. dard star frames by a normalized master flat-field image. Thereafter, (2015a, hereafter RET15) in the CLiF AGNs is the high value of the spectra were wavelength calibrated by applying the dispersion U V the H /H ratio, which in all objects studied was between 3.8 and solution obtained from the arc lamp frames. Finally, the spectra of 6.6 (taking into account the limits of the error bars). That is, higher standard stars were extracted and combined to derived the sensi- than the Balmer decrement assuming case B for low density gas. tivity function, later applied to the science 1D spectra. The final They interpreted this result as due to high density gas in the narrow products are wavelength and flux calibrated optical spectra. line region instead of being attributed to dust extinction because In all cases above, the final spectra were corrected for Galac- the values found for HW/HV are consistent with the intrinsic case B tic extinction using the extinction maps of Schlafly & Finkbeiner value of 0.47 (Osterbrock & Ferland 2006). (2011) (see column 6 of Table 1) and the Cardelli et al. (1989) law. The results obtained by RET15 suggest that most of the HU The final reduced spectra are shown in Fig. 1. emission is produced in the same region of the coronal line forest, enhancing the flux ratio HU/HV. The gas density they obtained by modelling the NLR is in the range 106 − 107.5cm−3, much higher than that expected for typical NLR conditions (< 104 cm−3 2.2 NIR Spectroscopy Osterbrock & Ferland 2006). NIR of the sample listed in Table 1 were obtained at three different RET15 classified only 7 CLiF AGNs from an initial sample of observatories: Cerro Tololo Inter-American Observatory (CTIO), ∼5000 AGN spectra from the Sloan Digital Sky Survery (SDSS). Gemini North, and the NASA Infrared Telescope Facility (IRTF). Here, we carry out the first near-infrared study of CLiF AGN In all cases, a cross-dispersed spectrograph, providing simultaneous in the literature along with a re-analysis of the optical spectra of coverage of the 0.8 – 2.4 `m region was employed. Below we these objects. We detail in Sect. 2 the sample and observations. provide basic description of the observations and data reduction Our additional data widen the number and species of coronal lines procedure for each set of data. suitable to derive the physical conditions of the coronal and narrow line region in these interesting sources (Sects 3 & 4). Stellar and gas kinematics properties are derived in Sect. 5. We also derive 2.2.1 GNIRS/Gemini spectroscopy additional properties of CLiF AGN such as the mass of the SMBH in order to construct a more complete picture of these sources. Main Near-infrared spectra of SDSS J124134.25+442639.2 and conclusions are drawn in Sect. 6. SDSS J164126.91+432121.5 were obtained in queue mode with the 8.1 m Gemini North telescope atop Mauna Kea (PROGRAM ID=GN2017A-Q40). The Gemini Near-IR spectrograph (GNIRS, Elias et al. 2006) in the cross-dispersed mode was employed. It al- lows simultaneous z+J, H and K band observations, covering the 2 SAMPLE, OBSERVATIONS AND DATA REDUCTION spectral range 0.8 − 2.5`m in a single exposure. GNIRS science de- The sample chosen for this analysis consists of the seven CLiF AGN tector consist of an ALADDIN 1k × 1k InSb array. The instrument already identified by RET15, observed using optical (3400 ≤ _ ≤ setup includes a 32 l/mm grating and a 0.8×7 arcsec slit, giving a 7500 Å) along with NIR (7800 ≤ _ ≤ 25000 Å) spectroscopy. Table 1 spectral resolution of R∼1300 (or 320 km s−1 FWHM). Individual list basic properties of the objects as well as the log of observations. exposures were taken, nodding the source in a ABBA pattern along Details about telescope setup and data reduction according to the the slit. Right after the observation of the science frames, an A0V wavelength interval are given below. Note that we did not collect star was observed at a similar airmass, with the purpose of flux optical spectra for 2MASX J113111.05+162739 nor NIR spectra calibration and telluric correction. for SDSS J124134.25+442639.2. The NIR data were reduced using the XDGNIRS pipeline
MNRAS 000, 1–20 (2019) A study of CLiF AGNs 3
Table 1. Sample of CLiF AGN as defined by Rose et al. (2015a) and employed in this work. Column 9 lists the classification into AGN of Type I or II according to Rose et al. (2015b) while the last column list the classification based on results obtained from this work (see Sect. 4.1).
◦ Galaxy Redshift Telescope/Instrument Date E(B-V)G Airmass Exp. Time PA Classification Classification (z) dd.mm.yyyy mag (s) (E of N) Rose et al. (2015b) this work ESO138-G001 0.0091 SOAR/Goodman 11.03.2017 0.176 1.13 3×900 279 Type-2 Type-1 Blanco/ARCoIRIS 07.04.2017 1.16 16×180 0.0 SDSS J164+43 0.2210 Gemini/GNIRS 15.02.2017 0.01 1.24 16×180 90 Type-2 Type-1 III Zw 77 0.0342 IRTF/SpeX 12.04.2015 0.09 1.08 18×180 338 Type-1 Type-1 MRK 1388 0.0213 IRTF/SpeX 04.05.2015 0.03 1.00 8×200 300 Type-2 Type-1 SDSS J124+44 0.0420 ...... 0.02 ...... Type-2 Type-2 2MASX J113+16 0.1740 Gemini/GNIRS 20.02.2017 0.03 1.25 12×160 90 Type-2 Type-1 NGC 424 0.0118 SOAR/Goodman 07.08.2014 0.01 1.02 4×900 60 Type-2 Type-1 Blanco/ARCoIRIS 19.09.2016 1.02 12×180 270
(v2.0) 1, which delivers a full reduced, wavelength and flux cali- angle. Differential refraction is, indeed, the main source of uncer- brated, 1D spectrum with all orders combined (Mason et al. 2015). tainty in flux calibration and it was minimised following the above Briefly, the pipeline cleans the 2D images from radiative events and procedure. Finally, the spectra were corrected for Galactic extinc- prepares a master flat constructed from quartz IR lamps to remove tion using the Cardelli et al. (1989) law and the extinction maps of pixel to pixel variation. Thereafter, The s-distortion solution is ob- Schlafly & Finkbeiner (2011). The spectra are shown in Fig. 3. tained from daytime pinholes flats and applied to the science and telluric images to rectify them. Argon lamp images are then used to find the wavelength dispersion solution, followed by the extraction 2.2.3 Blanco/ARCoIRIS data of 1D spectra from the combined individual exposures. The telluric NIR spectra of ESO 138-G001 and NGC 424 were obtained using features from the science spectrum are removed using the spectrum the ARCoIRIS spectrograph attached to the 4.1 m Blanco Tele- of a A0V star. Finally, the flux calibration is achieved assuming a scope (Schlawin et al. 2014). The science detector employed is a black body shape for the standard star (Pecaut & Mamajek 2013) 2048 × 2048 Hawaii-2RG HgCdTe array with a sampling of 0.41 scaled to its K-band magnitude (Skrutskie et al. 2006). The differ- arcsec/pixel. The slit assembly is 1.1 arcsec wide and 28 arcsec ent orders are combined in to a single 1D spectrum and corrected long. The delivered spectral resolution R is ∼3500 across the six for Galactic extinction using the Cardelli et al. (1989) law and the science orders. Similar to the IRTF/SpeX data, telluric standards extinction maps of Schlafly & Finkbeiner (2011). The spectra are were observed close to the science targets to warrant good telluric shown in Fig. 2. band cancellations and flux calibration. Data reduction was carried out using Spextool V4.1 with some modifications specifically designed for the data format and charac- 2.2.2 IRTF/SpeX data teristics of ARCoIRIS. The procedure follows the same steps and NIR spectra of III Zw 77 and Mrk 1388 were obtained at the NASA recipes as for SpeX. The final data product consist of a 1D spectrum, 3 m Infrared Telescope Facility (IRTF) using the SpeX spectrograph wavelength and flux calibrated with all individual orders merged (Rayner et al. 2003) in the short cross-dispersed mode (SXD, 0.7- forming a continuous 0.9 – 2.4 `m spectrum. We then corrected 2.4 `m). The detector consists of a Teledyne 2048 × 2048 Hawaii- these data for Galactic extinction using the Cardelli et al. (1989) 2RG array with a spatial scale of 0.10 arcsec/ pixel. A 0.8 arcsec x 15 law and the extinction maps of Schlafly & Finkbeiner (2011). arcsec slit, providing a spectral resolution, on average, of320km s−1 In order to verify the quality of the flux calibrations in the was employed. NIR spectra, we compared our data with the 2MASS photometric Observations were done nodding in two positions along the points. The only data that required a calibration adjustment were slit. Right before or after the science target, a telluric star, close those taken with ARCoIRIS (ESO 138-G001 and NGC 424). In in airmass to the former, was observed to remove telluric features order to re-scale the flux level to that of 2MASS, a third-degree and to perform the flux calibration. CuHgAr frames were also ob- polynomial function was applied using the Python LMFIT routine served at the same position as the galaxies for wavelength calibra- (Newville et al. 2016). The spectra are shown in Fig. 4. tion. The spectral reduction, extraction and wavelength calibration procedures were performed using Spextool V4.1, an IDL-based software developed and provided by the SpeX team for the IRTF 3 CONTINUUM EMISSION ANALYSIS community (Cushing et al. 2004). Telluric features removal and flux calibration were done using XTELLCOR (Vacca et al. 2003), an- In this section we describe the modelling of the continuum emission other software available by the SpeX team. The different orders in the optical and NIR in order to remove the effects of both the were merged into a single 1D spectrum from 0.7 to 2.4 `m using underlying stellar population of the host galaxy and the continuum the XMERGEORDERS routine. After this procedure, the IDL rou- emission attributed to dust heated by the AGN. The first component tine Xlightloss, also written by the SpeX team, was employed. It is observed in both spectral regions while the second is restricted to corrects an input spectrum for slit losses relative to the standard star the NIR. used for flux calibration. This program is useful if either the object or the standard were not observed with the slit at the parallactic 3.1 Optical continuum A rapid inspection to the optical spectra of the CLiF sample (see 1 Based on the Gemini IRAF packages Fig. 1 of RET15 and Fig. 1 of this manuscript) reveals the presence
MNRAS 000, 1–20 (2019) 4 F. C. Cerqueira-Campos et al.
Figure 1. Optical spectra obtained by Goodman/SOAR. In the bottom panel for both spectra a zoom around the continuum level is displayed to show the faintest lines identified. The ions belonging to coronal line transitions are represented in red.
MNRAS 000, 1–20 (2019) A study of CLiF AGNs 5
Figure 2. NIR spectra obtained by GNIRS/Gemini. The shaded area represents the region of low atmospheric transmission. Lines located in these regions were not employed in the analysis. In the bottom panel, for both spectra, a zoom is displayed to identify the faintest lines. Coronal lines are identified with red labels.
MNRAS 000, 1–20 (2019) 6 F. C. Cerqueira-Campos et al.
Figure 3. Same as Fig. 2 for IRTF/Spex spectra
MNRAS 000, 1–20 (2019) A study of CLiF AGNs 7
Figure 4. Same as Fig. 2 for Blanco/ARCoIRIS spectra
MNRAS 000, 1–20 (2019) 8 F. C. Cerqueira-Campos et al.
Galaxy Galaxy 10 Fit 6 Galaxy - Fit Galaxy - Fit Fit 5 8 Stellar Power law 4 6 3
4 2
ESO 138 ESOG1 138 2 1 SDSS J164+43 0 0 7 Galaxy 10 Galaxy 6 Fit Fit Galaxy - Fit Galaxy - Fit 5 8
4 6 flux 3 4 2 III ZW 77 2 MRK 1388 1
0 0
Normalized 7 4 Galaxy Galaxy Fit 6 Fit Galaxy - Fit Galaxy - Fit 3 5
4 2 3
1 2 NGC 424 NGC
SDSS J124+44 1 0 0 4000 4500 5000 5500 6000 6500 7000 4000 4500 5000 5500 6000 6500 7000 Wavelength (Å)
Figure 5. Modelling of the stellar contribution for CLiF AGN in the optical region. In all panels, the observed spectrum is in black, the stellar continuum is in red and the nebular emission obtained after subtraction of the stellar continuum is in blue.
of strong absorption lines, evidencing the contribution of stellar 3.2 The NIR continuum population from the host galaxy to the observed integrated spectra. It is not of interest in this paper to study the stellar population, Accepting the unified model to describe the different types of AGNs, but to model and subtract this contribution to obtain suitable flux the presence of a dusty torus leaves a spectroscopic signature that measurements for the gas in the central region. This procedure starts to show up in the NIR. Dust absorbs a fraction of the op- particularly impact the measurements of the Balmer decrement and tical/ultraviolet continuum from the central source and radiates it the flux of lines that are partially absorbed by the stellar continuum. back at longer wavelengths, from 1 `m up to the far-infrared. Observational evidence corroborates models describing the NIR continuum in AGNs from 1 to 10 `m as being predominantly or entirely dominated by hot dust (Edelson & Malkan 1986; Barvainis 1987). In order to remove the stellar contribution, the Penalized Rodriguez-Ardila & Mazzalay (2006) confirmed this scenario Pixel-Fitting (pPXF) software developed by Cappellari & Emsellem in the Seyfert 1 AGN Mrk 1239 by adjusting the excess of NIR con- (2004) and updated by Cappellari (2017) was employed. This rou- tinuum over the power-law extrapolated from the optical region, tine uses maximum penalized likelihood estimation with the princi- as due to dust emission. They employed a simple blackbody func- ple of searching for the best spectral fit by combining stellar spectra tion at a temperature of 1200 K. This is close to the sublimation from a library of simple stellar populations. Here, the templates that temperature of dust grains (∼1500 K). were used to generate the pPXF fits are those of Bruzual & Charlot Thus, to fit the continuum of our sample of CLiF AGN in (2003). Fig. 5 shows the fits carried out to the galaxy sample and the the NIR, in addition to the stellar population templates, a function resulting nebular emission after subtraction of the stellar compo- that represents the hot dust emission should be employed. For the nent. In addition, pPXF also estimates the stellar velocity dispersion, former component, the IRTF library of stellar spectra Rayner et al. which will be employed in the determination of the mass of the black (2009) was used to fit the data. It is composed of 296 spectra in hole (see Sect. 5.2). the wavelength range 0.8 to 5.0 `m, with a resolving power of
MNRAS 000, 1–20 (2019) A study of CLiF AGNs 9