The Close AGN Reference Survey (CARS). Discovery of a global [CII] 158 $\mu$m line excess in AGN HE1353-1917

The [CII]158$\mu$m line is one of the strongest far-infrared (FIR) lines and an important coolant in the interstellar medium of galaxies that is accessible out to high redshifts. The excitation of [CII] is complex and can best be studied in detail at low redshifts. Here we report the discovery of the highest global [CII] excess with respect to the FIR luminosity in the nearby AGN host galaxy HE1353-1917. This galaxy is exceptional among a sample of five targets because the AGN ionization cone and radio jet directly intercept the cold galactic disk. As a consequence, a massive multiphase gas outflow on kiloparsec scales is embedded in an extended narrow-line region. Because HE1353-1917 is distinguished by these special properties from our four bright AGN, we propose that a global [CII] excess in AGN host galaxies could be a direct signature of a multiphase AGN-driven outflow with a high mass-loading factor.


Introduction
The [C ii] 157.74 µm emission line arises from the fine-structure transition 2 P 3/2 → 2 P 1/2 of the ground state of singly ionized carbon C + (ionization potential of 11.2 eV).Working as a coolant in multiple phases of the interstellar medium (ISM), the [C ii] line is one of the brightest emission lines in the far-infrared (FIR); it contributes 0.1-0.3% to the FIR luminosity.
The [C ii] line has been calibrated as a probe for the cold gas content and associated star formation rates (SFR) in galaxies (Stacey et al. 1991;Boselli et al. 2002;Herrera-Camus et al. 2015).However, using the [C ii] line as an SFR tracer is complex because of the multiple mixed excitation mechanisms of the line.In local star-forming galaxies, 66-82 % of [C ii] arises from the neutral gas of photodissociation regions (PDRs), and the rest comes from the ionized phase (Croxall et al. 2017).

Lyman Spitzer Jr. Fellow
The mechanism of dust infrared emission, on the other hand, is rather simple: dust preferentially absorbs UV radiation from the stellar population and therefore is sensitive to the bright young stars.The infrared luminosity has been well calibrated as an SFR tracer at SFR >1 M yr −1 (Hirashita et al. 2003;Murphy et al. 2011) so that the [C ii] and FIR luminosity are expected to be correlated.
At the highest SFRs, luminous and ultra-luminous infrared galaxies (U/LIRGs) exhibit the so-called [C ii] line deficit (Helou et al. 2001;Malhotra et al. 2001;Luhman et al. 2003), where [C ii] becomes unreliable as an SFR indicator (Díaz-Santos et al. 2013).The origin of the line deficit is still debated and is directly connected to physical processes that are crucial for understanding the [C ii] excitation mechanisms.
Active galactic nuclei (AGN) are able to affect the [C ii]/FIR ratio in several ways: they can increase the infrared luminosity through dust heating (Herrera-Camus et al. 2018a); act as an ad-  3+, etc. with their hard radiation field (Langer & Pineda 2015).Herschel surveys of nearby galaxies such as KINGFISH (Smith et al. 2017) and SHINING (Herrera-Camus et al. 2018b) found no link between the [C ii] line emission and AGN luminosity, but these AGN may be not luminous enough to outshine the star formation.Given its brightness, the [C ii] line is the most important ISM diagnostic at high redshifts that can be observed with unprecedented spatial resolution and depth on submillimeter interferometers.The sample of the observed high-redshift objects includes starburst and AGN-dominated systems, with the [C ii]/FIR ratios spanning a wide range from 0.02 % to 5 % (e.g., Gullberg et al. 2015;Brisbin et al. 2015;Decarli et al. 2018).In order to provide an interpretation, we need to investigate the [C ii] emission in local galaxies and determine the effect of luminous AGN.
In this letter, we present [C ii] line observations with the Stratospheric Observatory For Infrared Astronomy (SOFIA, Temi et al. 2014) for five nearby (0.024 < z < 0.040) luminous Seyfert 1 AGN host galaxies from the Close AGN Reference Survey (CARS; Husemann et al. (2017); www.cars-survey.org)to investigate the impact of AGN on the global [C ii] luminosity.

SOFIA/FIFI-LS observations
We observed five CARS objects with the Far Infrared Field-Imaging Line Spectrometer (FIFI-LS; Klein et al. (2014)) on board SOFIA.The objects were picked to cover a broad range of SFRs (1-11 M yr −1 ) and avoid the strong atmospheric absorption in the redshifted [C ii] line wavelength region.The SFR esti-mates were initially based on the predictions calculated from the AGN-subtracted extinction-corrected Hα line from observations obtained with the Mulit-Unit Spectroscopic Explorer (MUSE, Bacon et al. 2010, section 2.2).The properties of the observed galaxies are listed in Table 1.
The observations were performed during SOFIA Cycle 4 (plan ID 04_0056, PI: Husemann) and Cycle 5 (plan ID 05_0077, PI: Husemann).FIFI-LS is a double-beam spectrometer that covers 1 × 1 in the red channel (105-200 µm) and 0.5 ×0.5 in the blue channel (50-125 µm), split up into 5×5 spatial pixels.We tuned the setups to cover the [C ii] line in the red channel with spectral resolution of R∼ 1200 (250 km s −1 ) and either [O iii] 88 µm or [O i] 63 µm in the blue channel, depending on atmospheric transmission.
The pipeline-processed data are provided by the FIFI-LS team.We used LEVEL_3 science-ready data, which consist of a number of ∼30 second exposures, to apply an additional selection (see Appendix A for details) and background subtraction, and constructed data cubes with 6 sampling using the Drizzle algorithm (Fruchter & Hook 1997).To derive total [C ii] line fluxes from the FIFI-LS cubes, we summed the spectra within an aperture with 36 diameter and fit the line shape with Gaussian profiles.The HE 1108−2813and HE 2211−3903  used in the [O iii] upper limit estimation, taking into account the spectral resolution of R∼ 670 (450 km s −1 ) for the observed wavelengths.

VLT/MUSE observations
All five CARS targets were observed with MUSE at the Very Large Telescope (VLT) under ESO programs 094.B-0345(A) and 095.B-0015(A).The MUSE data cover 1 × 1 FoV at a 0 . 2 spatial sampling and a wavelength coverage of 4650-9300 Å with R ∼ 2500.Integration times range from 400 s-900 s split up into two or three exposures, which are rotated by 90 • against each other for cosmic-ray rejection and better image cosmetics.The data were reduced with the standard MUSE pipeline (version 1.6.0,Weilbacher et al. 2012Weilbacher et al. , 2014)).
As a first step in the post-processing, we subtracted the bright point-like AGN emission from the reconstructed datacube using QDeblend 3D , as described in Husemann et al. (2013Husemann et al. ( , 2014Husemann et al. ( , 2019)).Afterward we modeled the stellar continuum and ISM emission lines in the AGN-subtracted datacube with PyParadise (see Husemann et al. 2016;Weaver et al. 2018, for details).
In Fig. 2 we show the classical Baldwin-Philips-Terlevich (BPT, Baldwin et al. 1981) diagrams and corresponding spatial maps for the targeted galaxies after binning 8×8 pixels.Essentially, the BPT diagnostic diagrams highlight different dominating ionization mechanisms.The face-on galaxies display line ratios that are mainly consistent with star formation or intermediate line ratios between star formation and AGN excitation.The edge-on galaxy HE 1353−1917 is instead dominated by AGN photoionization.The fraction of AGN excitation to the extinction-corrected Hα for HE 1353−1917 is >80%, whereas a fraction of <40% is found for all face-on galaxies.

SED fitting and FIR luminosities
We constructed spectral energy distributions (SEDs) for all targets using publicly available broad-band photometry from the following catalogs: the Herschel/SPIRE Point Source Catalogue 1.17 ± 0.41 1.94 ± 0.25 (*) Integrated over 42.5-122.5µm (Schulz et al. 2017), the AKARI/FIS Bright Source Catalogue (Yamamura et al. 2009), the 2MASS Extended Source Catalog (Jarrett et al. 2000), and the GALEX Source Catalog (Bianchi et al. 2017).For the following datasets we applied aperture measurements on the survey images: Herschel/PACS images from the Herschel Science Archive2 , WISE Image Atlas (Cutri et al. 2011), optical grizy photometric images3 from the Pan-STARRS DR1 (Chambers et al. 2016), and Swift/UVOT images from the Barbara A. Mikulski Archive for Space Telescopes4 .
We obtained the FIR (42.5-122.5 µm) luminosities through SED fitting with AGNfitter (Calistro Rivera et al. 2016).AG-Nfitter uses an MCMC approach and various template libraries to model the SED as a super-position of emission from an AGN accretion disk, a torus of AGN-heated dust, stellar light from the galaxies, and cold dust in star-forming regions.The AGNfitter output SED models are shown in Appendix B and the FIR (42.5-122.5 µm) luminosities are listed in the Table 1.

Results and discussion
The [C ii] and FIR luminosities of our five targets are shown in Fig. 3  Based on the FIR and AGN-subtracted extinction-corrected Hα luminosity, we expected HE 1353−1917 to be the faintest target in [C ii], but it turned out to be the brightest source of our sample.While all galaxies have similar metallicities, stellar masses, and AGN bolometric luminosities, the obvious dissimilarity is the edge-on orientation of HE 1353−1917.When we consider the unobscured nature of this AGN, this means that the AGN ionization cone directly pierces the gas-rich disk of the galaxy.This leads to a large biconical extended narrow-line region (ENLR) that is oriented almost along the disk axis on kiloparsec scales, as clearly shown in Fig. 2. In addition, a massive multiphase outflow on kiloparsec scales with a mass outflow rate of Ṁout ∼10-100 M yr −1 is detected in this galaxy, as discussed in detail by Husemann et al. (2019), while less prominent outflows are detected in the other four galaxies (Singha et al. in prep.).The global SFR of ∼2M yr −1 implies an integrated mass-loading factor of ∼10 or more that has a similar scale as the observed [C ii] excess.
Several powering sources may contribute to the [C ii] line excess in HE 1353−1917, but the challenge is to distinguish the dominant [C ii] line excitation mechanism.As highlighted in Fig. 3 2010), we estimate an XDR contribution of only 10% given an X-ray luminosity of L 2−10 keV = 1.69 × 10 43 erg s −1 (Husemann et al. 2019).
[C ii] emission can originate in any gas phase that is illuminated by UV photons.Using other line and continuum diagnostics, we can obtain some idea of which phases contribute to the [C ii] emission, for example, 66-82 % for the neutral phase and the rest for the ionized phase in local star-forming galaxies (Croxall et al. 2017).However, these fractions do not necessarily remain the same in the ENLR of an AGN.How much of the [C ii] emission originates from an ENLR has not been systematically explored so far.Dissipation of the kinetic energy of shocks and outflows can also have important consequences on the [C ii] emission, as theoretically explored by Lesaffre et al. (2013).[C ii]/FIR ratios of 3-7% are detected in between merging galaxies (Appleton et al. 2013;Peterson et al. 2018) and locally within galaxies (Appleton et al. 2018).On the global scale, the radio galaxy 3C 326 emits around 3% of FIR in [C ii] line for which jet-driven turbulence is likely responsible (Guillard et al. 2015).Nesvadba et al. ( 2010 loading factor of >10 similar to that found in HE 1353−1917.HE 1353−1917 has a global [C ii]/FIR ratio of 4.3 ± 0.4% and a radio-jet powering the multiphase gas outflow (Husemann et al. 2019).The [C ii] line width of 760 ± 60 km s −1 is broader than the CO(1-0) line for cold molecular gas outflow ∼210 km s −1 , but narrower than the ionized gas outflow ∼1020 km s −1 (see Fig. 1).
[C ii] therefore likely traces the interface between the warm and cold gas phase within the outflow.We can relate the region where the [C ii] excess occurs to the center of the galaxy (Fig. 4) where the brightest region of the ENLR and the kpc-scale multi-phase outflow are indeed located.The [O iii] upper limit implies a line ratio limit of [O iii]/[C ii] < 0.96, which suggests that the pure AGN ionization in the ENLR is not the primary cause for the [C ii] excess.Hence, the observations are firm evidence that the [C ii] line excess in HE 1353−1917 is related to the multiphase outflow initiated by the jet.
High-redshift galaxies span an entire range between [C ii] deficient and normal regimes with [C ii]/FIR ratios from 0.02 % to 5 %.The large scatter in the ratios may be caused by the orderof-magnitude difference in spatial resolutions, for instance, 18 for CSO, Brisbin et al. (2015) with average ratio of 1.8 % and 1 for ALMA, Decarli et al. (2018)

Conclusion
We presented the discovery of a global [C ii] line excess in one out of five AGN host galaxies.Based on ancillary information from an extensive multiwavelength analysis of this galaxy as part of the CARS survey, we can directly connect the [C ii] line excess in HE 1353−1917 to the impact of the AGN that drives a massive multiphase outflow on kiloparsec scales embedded in an ENLR (Husemann et al. 2019).The detection of such a global [C ii] excess in AGN host galaxies is of crucial importance for the interpretation of [C ii] line observations and the detection of massive gas outflows in luminous high-redshift AGN host galaxies.Given the evidence of HE 1353−1917 and the similar outflow seen in 3C 326, we propose that a significant [C ii] line excess in luminous AGN, if detected, can be used as an inference for a multiphase AGN outflow with a high mass-loading factor even at high redshifts.

Appendix A: FIFI-LS time window selection
Although SOFIA reduces the water vapor absorption to 99% of the ground level, the atmospheric variations still play a perceptible role and affect the quality of the FIFI-LS data.To probe this effect, we used the science-ready LEVEL_3 pipeline output files.Each file corresponds to only ∼ 30 seconds of exposure time and cannot be examined independently due to the high levels of noise.Calculating [C ii] line flux signal-to-noise ratio (S/N) of cumulatively summed files allows us to trace the trends of the S/N and therefore judge the atmospheric conditions.Ideally, the trend should rise as the square root function.If the S/N remains the same or even decreases when exposure time is added, we assume that an atmospheric variation degrades the quality and exclude those time windows from the analysis, as shown with the shaded areas in Fig. A.1.
To calculate the S/N as a function of the exposure time, we fit a fixed Gaussian shape taken from the fitting result before the additional selection to the spectra, summed within the 36 diameter aperture.In the case of HE 1353−1917, when the line shape is not Gaussian, the procedure still works similarly because only relative trends are important in this analysis.Masking the spectral regions with strong atmospheric absorption lines is, conversely, very important for the flux estimations.The first few values of the S/N as a function of the coadded files are unreliable, but the overall trends remained when we performed the reverse-order coadding test.
The rapid altitude shifts also affect the S/N trends as seen in the case of HE 1353−1917and HE 2211−3903.In the last 15 files for HE 2211−3903 the altitude is higher and of much better quality, but we decided not to exclude the majority of the files from the analysis and kept all of them.
Even though the selection due to atmospheric conditions is performed in the reduction pipeline, our additional selection technique helps to increase S/N of the detection by excluding the poor quality artifact-rich files.

Fig. 1 .
Fig. 1.SOFIA/FIFI-LS [C ii] spectra (black lines) for each object integrated within 36 diameter apertures.The CO(1-0) and best-fit Gaussian or double-Gaussian profiles are shown as red dashed and solid green lines, respectively.The CO(1-0) and [O iii]5007Å line profiles of the HE 1353−1917 outflow are shown as orange dashed and dotted lines.The comparison line profiles are degraded to the spectral resolution of SOFIA.The atmospheric transmission curves are shown in blue, and the shaded regions of low transmission are excluded from the analysis.

Fig. 2 .
Fig. 2. Upper panel: MUSE BPT diagnostic diagrams with empirical dividing lines.The solid line shows data from Kauffmann et al. (2003) and the dashed line from Kewley et al. (2001).The colors represent AGN ionization (red), ionization by H ii regions (blue), and mixed regime (green).Lower panel: MUSE emission line 1.6" × 1.6" binned maps color-coded according to the BPT diagram regions; white-light contours are overplotted in black.
Fig. 3. [C ii] line luminosity as a function of FIR luminosity.Our five CARS targets are shown as red triangles compared to the literature compilation of Herrera-Camus et al. (2018a): normal star-forming and star-burst galaxies (green circles), AGN host galaxies (purple triangles), LINER galaxies (orange upside-down triangles), and high-redshift galaxies including the samples fromBrisbin et al. (2015) andDecarli et al. (2018) (brown diamonds).The blue line is the mean of the low-redshift galaxies distribution, and the blue-shaded areas correspond to 1, 2, and 3 σ.Pink and purple shades represent LIRG L FIR > 10 11 and ULIRG L FIR > 10 12 regimes.The luminosities for 3C 326 are taken fromGuillard et al. (2015).
, the [C ii] line luminosity of HE 1353−1917 cannot be powered by the star formation alone.The SFRs calculated from the extinction-corrected Hα (SFR Hα = 1.23±0.03M yr −1 ) and 42.5-122.5µm luminosity (SFR FIR = 2.3 ± 0.1 M yr −1 ) can account for only about 25% of the observed [C ii] line luminosity (SFR [C ii] = 2.286 • 10 −43 × L 1.034 [C ii] , Herrera-Camus et al. 2015).X-ray dominated regions (XDRs) produced by the hard X-ray photons from an AGN may contribute and even dominate the PDR [C ii] emission.From the scaling relation L [C ii], XDR = 2 × 10 −3 L 2−10 keV by Stacey et al. ( If the ionized gas phase produces more than 20-40 % of [C ii] as in PDR paradigm, then it can explain the observed [C ii] excess.In HE 1353−1917 only ∼20 % of Hα originates from star-forming regions and ∼80 % comes from the AGN-ionized regions.If we naively assume that 20 % of the [C ii] flux originates from star formation, the [C ii] SF of HE 1353−1917 falls within 3σ of the [C ii]-FIR relation. Fig. 4. SOFIA/FIFI-LS [C ii] line flux map of HE 1353 − 1917.The contours represent Herschel PACS 70 (solid green lines) and 160 µm (blue dashed line) photometric observations.
with average ratio of 0.08 %, as low-resolution observations can include a larger fraction of intergalactic [C ii] emission.Several high-redshift QSO have shown a [C ii] enhancement (Maiolino et al. 2009; Wagg et al. 2010), which was interpreted as due to a low metallicity given the high [C ii]/FIR and [C ii]/CO(1 − 0) ratios in local dwarf galaxies.We measured [C ii]/FIR ∼ 4% and [C ii]/CO(1 − 0) ∼ 10 4 for HE 1353−1917, which can neither be explained by intergalactic [C ii] emission nor by a low metallicity.

Fig
Fig. A.1.[C ii] line S/N as a function of the cumulatively coadded exposures (black circles).The trends (gray lines) are purely subjective and are shown to guide the eye.The shaded areas cover the time windows and associated files that were excluded from the further analysis.The corresponding altitudes of SOFIA are shown as blue lines.

Table 1 .
Measured galaxy properties