Open Access
Issue
A&A
Volume 659, March 2022
Article Number A125
Number of page(s) 38
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202142011
Published online 17 March 2022

© I. Smirnova-Pinchukova et al. 2022

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Open Access funding provided by Max Planck Society.

1. Introduction

The evolution of galaxies is driven by many different physical processes, which we can indirectly study by probing various galactic properties across cosmic time. One key parameter is the star-formation rate (SFR), which changes with redshift and diskloses essential information about the evolution of galaxies (e.g., Madau & Dickinson 2014). Furthermore, at lower redshifts the galaxy population shows a clear bi-modality in the color-magnitude diagram (e.g., Strateva et al. 2001; Baldry et al. 2004; Wyder et al. 2007; Brammer et al. 2009) of which the specific SFR versus stellar mass (sSFR–M) diagram (e.g., Brinchmann et al. 2004; Noeske et al. 2007; Karim et al. 2011; Renzini & Peng 2015) is a more physical representation. Such diagrams separate star-forming from passive galaxies where the so-called star-forming main sequence (SFMS, Noeske et al. 2007) indicates the average location of star-forming galaxies as a function of stellar mass. Moreover, the SFMS tends to be present across a wide range of redshifts as a tight correlation (e.g., Elbaz et al. 2007; Speagle et al. 2014). How this galaxy bi-modality is established is one of the open key questions in our current galaxy evolution model. One option is that accreting super-massive black holes (SMBHs) release enough energy to suppress the star formation in their host galaxies as they pass quickly through from the blue cloud to the red sequence (e.g., Faber et al. 2007; Schawinski et al. 2014), but also opposite tracks via rejuvenation of star forming have been proposed (e.g., Hasinger 2008; Mancini et al. 2019; Chauke et al. 2019). Furthermore, environmental processes (e.g., Balogh et al. 2004; Peng et al. 2010; Schaefer et al. 2019) and gas consumption (e.g., Colombo et al. 2020) may facilitate this transformative process.

The ability both to rapidly quench and maintain low levels of star formation in red-sequence massive galaxies has often been attributed to the impact of active galactic nuclei (AGN). Radiatively-driven winds, mechanical energy, and/or thermal pressure release from the AGN are physical processes routinely incorporated into numerical simulations of galaxy evolution to solve the overcooling problem (e.g., Vogelsberger et al. 2014a; Schaye et al. 2015; Pillepich et al. 2018; Bassini et al. 2019; Wittor & Gaspari 2020), so that the characteristics of the observed galaxy population can be reproduced (e.g., Vogelsberger et al. 2014b; Nelson et al. 2018). Often two types of the AGN feedback mode are being distinguished, which are called the quasar mode or radiative mode (e.g., Di Matteo et al. 2005; Hopkins et al. 2008; Hopkins & Elvis 2010) and the radio mode also known as kinetic or maintenance mode (e.g., Bower et al. 2006; Croton et al. 2006; Fabian 2012; Gaspari et al. 2020). The quasar mode is thought to initially quench star formation through powerful radiatively driven winds that remove the gas from the galaxy (e.g., Nesvadba et al. 2008; Feruglio et al. 2010; Maiolino et al. 2012) and thereby suppress star formation. The radio mode is a heating mechanism affecting the gaseous halo that prevents continuous condensation of cold and warm gas, thereby quenching star formation over long timescales and maintaining a low level of star formation in red sequence galaxies (e.g., Brüggen & Kaiser 2002; McNamara et al. 2005; Gaspari et al. 2019; McDonald et al. 2021). Strong observational evidence for radio mode feedback has been collected from galaxy clusters hosting powerful radio galaxies (e.g., Fabian et al. 2006; McNamara & Nulsen 2007). The role of quasar mode feedback is much more controversial from an observational perspective. Various studies have reported negative (Ho 2005; Nandra et al. 2007; Schawinski et al. 2009; Farrah et al. 2012; Page et al. 2012; Mullaney et al. 2015; Shimizu et al. 2015; Wylezalek & Zakamska 2016; Kakkad et al. 2017; Catalán-Torrecilla et al. 2017; Bing et al. 2019; Bluck et al. 2020; Ramos Padilla et al. 2020; Brownson et al. 2020; Smith et al. 2020), positive (Kim et al. 2006; Cresci et al. 2015a,b; Bernhard et al. 2016; Santoro et al. 2016; Maiolino et al. 2017; Koss et al. 2021), or no effect from feedback (Elbaz et al. 2011; Bongiorno et al. 2012; Harrison et al. 2012; Husemann et al. 2014; Balmaverde et al. 2016; Leung et al. 2017; Woo et al. 2017; Shangguan et al. 2018; Scholtz et al. 2020) on the star formation rates of observed AGN host galaxies.

One common difficulty in observational studies is the measurement of the basic AGN host galaxy parameters such as SFR and stellar mass, because the AGN light can dominate the total emission of the host galaxy across nearly the entire electro-magnetic spectrum from radio to X-rays. Hence, it can be challenging to measure the SFR or stellar mass from commonly used proxies and various strategies have been used in the literature to mitigate the issue. The stellar mass can be easily determined from optical spectra or NIR-optical photometry in obscured AGN (e.g., Kauffmann et al. 2003; Silverman et al. 2009; Aird et al. 2012; Mignoli et al. 2013) as the light from the central AGN is blocked by dust along our line of sight. However, the drawback is that AGN parameters such as black hole (BH) mass can usually not be directly inferred for obscured AGN, except for few very nearby systems using the narrow-line region (NLR) kinematics (e.g., Barth et al. 2001; Walsh et al. 2013).

The SFR of a galaxy is usually determined either from the UV radiation of young stars (e.g., Salim et al. 2007), the ionizing photons from H II regions using Hα or [O II] (e.g., Kennicutt & Kent 1983; Kewley et al. 2004), the reheated warm and cold dust at MIR and FIR wavelengths (e.g., Kennicutt 1998a), or the radio emission produced in supernova explosions (e.g., Condon 1992). All those tracers can be significantly contaminated by AGN emission even in obscured AGN. The UV or radio emission can be spatially separated through high-angular observations (Fernández-Ontiveros et al. 2018; Rosario et al. 2021). Alternatively, the overall spectral energy distribution (SED) is decomposed into various emitting components of a galaxy including AGN, star-formation heated dust and stars. Examples of such SED fitting codes are (X-)Cigale (Noll et al. 2009; Yang et al. 2020), Magphys (da Cunha et al. 2008), AGNfitter (Calistro Rivera et al. 2016) or FortesFit (Rosario 2019), which can be used to infer stellar mass, IR-based SFR and AGN luminosity simultaneously even for unobscured type 1 AGN. One potential drawback of IR-based SFR tracers for AGN feedback studies is that it traces the SFR on relatively long timescales of ∼100 Myr and therefore the SFR can be significantly over or underestimated for a strongly declining or rising SFR within that timescale (Hayward et al. 2014).

Considering that the quasar-mode feedback is expected to be a fast process, the Hα line may be a more viable tracer for the feedback, considering that it probes the current SFR on < 5 Myr timescales. Unfortunately, the Hα line can be strongly boosted by AGN ionization. In order to use Hα as a SFR tracer, optical emission lines diagnostic diagrams such as the Baldwin-Phillips-Terlevich diagrams (Baldwin et al. 1981) have been extensively used to distinguish between AGN and H II region excitation. A comprehensive classification scheme with such diagnostic diagrams has been developed Kewley et al. (2006). In particular with integral-field unit (IFU) spectroscopy, the usage of BPT diagrams has been extended to map the ionization conditions across galaxies to exclude regions from the SFR estimation dominated by AGN ionization using such BPT demarcation lines (e.g., Husemann et al. 2014; Nascimento et al. 2019). However, the superposition of different ionization conditions along the line-of-sight is a common scenario. The observed AGN mixing sequence (MS) arises from such a superposition, and an analytic approach to separate the relative contribution was presented in Davies et al. (2014a,b, 2016, 2017).

In this paper, we infer and compare the SFR estimated from Hα and the IR for a sample of 41 galaxies as part of the Close AGN Reference Survey (CARS1). CARS aims to shed light on the relationship between the central AGN and its host galaxy by following a spatially resolved multiwavelength approach. As presented in Husemann et al. (2022), the 41 CARS targets are a subsample of 99 luminous unobscured AGN at z < 0.06 previously identified with the flux-limited Hamburg/ESO AGN survey (Wisotzki et al. 2000). The AGN targeted by CARS therefore correspond to the most luminous unobscured AGN in the nearby Universe allowing a detailed investigation of AGN – host galaxies interactions at relatively high spatial resolution.

In this work we use both panchromatic SEDs (to estimate stellar masses and FIR-based SFRs), and analyze the CARS IFU data to infer Hα based SFRs (adopting a new algorithm to separate the contribution to the Hα-line from other excitation mechanisms). In addition, we explore the use of different prescriptions to predict the expected SFR from the non-AGN galaxy population as control samples. This allows us to investigate links with AGN parameters depending on the actual SFR tracers and control sample properties being used as well as to distinguish between potentially declining and increasing star formation histories amongst the AGN host galaxies. Throughout the paper, we assume a flat cosmological model with H0 = 70 km s−1 Mpc−1, ΩM = 0.3, and ΩΛ = 0.7.

2. Observational data set

2.1. Archival broad-band photometry

We collect publicly available broad-band photometry from archival sources combined with dedicated observations for CARS targets to construct their panchromatic spectral energy distributions (SEDs). Where catalogs are not available or do not properly handle extended sources, we extract consistent aperture photometry from the public survey images. All photometric measurements are listed in Table 1 and their origin is described below (see Fig. 1 for photometric band compilation).

thumbnail Fig. 1.

Photometric broad bands from a number of surveys and instruments used for the SED modeling. The y-axis is normalized roughly by the year of observational campaign to visually compare bands with similar coverage. When more than one option was available for certain bands, the preference was given to the band sets, which are temporally closer to each other (especially in near infrared and optical ranges) to minimize the impact of AGN variability.

Table 1.

Multiband aperture photometry measurements.

2.1.1. Catalog data

We directly took data from the Herschel/SPIRE point source catalog (250 μm, 350 μm, and 500 μm bands, Schulz et al. 2017), the AKARI/FIS bright source catalog (N60, N160, WIDE-L, and WIDE-S bands, Yamamura et al. 2009), and the 2MASS extended source catalog (J, H, and Ks bands, Jarrett et al. 2000) without further processing. The GALEX source catalog (FUV and NUV bands, Bianchi et al. 2017) were used only to compare the fluxes to our own aperture photometry measurements.

2.1.2. Aperture photometry

For the following data sets we performed customized aperture measurements on the survey images: Herschel/PACS 70, 100, and 160 μm images from the Herschel science archive2, WISE atlas images (Cutri et al. 2011), 2MASS J, H and Ks images from the interactive image service3, VISTA/Paranal J, H, and Ks band images from the Cambridge astronomical survey unit4, optical grizy photometric images5 from the Pan-STARRS DR1 (Chambers et al. 2016) and ugriz images6 from the SDSS-III DR12 (Alam et al. 2015), and GALEX and Swift/UVOT ultraviolet images from the Barbara A. Mikulski Archive for Space Telescopes7.

The elliptical apertures are defined for each instrument to cover the entire galaxy in the sequence of the bands. The aperture fluxes are computed using the AstroPy Python package PHOTUTILS8 (Bradley et al. 2020). Since the images of CARS targets are extended at these wavelengths the apertures are not affected by beam smearing and PSF matching as done for higher redshift targets is not a concern for our sample.

2.2. Galactic extinction and AGN broad-line correction

The observed photometry for external galaxies is attenuated by the dust in the Milky Way depending on the specific line-of-sight. We apply an extinction correction to the NIR, optical, and UV photometry of the SED using the Fitzpatrick (1999) Milky Way attenuation curve. The normalization of the attenuation curve is anchored to the measured V band extinction (AV) along our line-of-sight toward the target is taken from Schlafly & Finkbeiner (2011) and retrieved from the NASA/IPAC Extragalactic Database (NED)9. All corresponding photometric data in Table 1 are provided after extinction correction.

The broad emission lines of unobscured AGN can significantly contribute to broad-band photometry. Because the SED model library used for the fitting assumes smooth power-law for the optical-UV accretion disk component we need to subtract the broad emission line contribution from the photometry. We used the available IFU spectroscopy of the CARS sample to measure Hβ and Hα fluxes, and subtract them from the optical bands (mainly g and r). The values presented in Table 1 are presented after correction for the BLR contribution.

2.3. Dedicated broad-band imaging

Given the lack of available data in some bands for some of the CARS targets, we obtained dedicated observations to fill the missing gaps in our photometry. Optical observations were mainly missing in the u band which Pan-STARRS does not cover and HE 0108−4743 and HE 2211−3903 are not covered by SDSS or Pan-STARRS, given their low declination. Therefore, we obtained deep wide-field optical imaging with the Dark Energy Camera (DECam, DePoy et al. 2008) mounted to the 4 m Blanco telescope of the Cerro Tololo Inter-American Observatory (CTIO) and with the Wide-Field Imager (WFI, Baade et al. 1999) mounted to the 2.2 m telescope of the La Silla observatory. Missing NIR imaging were obtained with the PAnoramic Near-Infrared Camera (PANIC, Cárdenas Vázquez et al. 2018) mounted to the 2.2 m telescope at the Calar Alto Observatory. Finally, deeper FIR observations were obtained for most of the CARS targets undetected with AKARI using the High-resolution Airborne Wideband Camera Plus (HAWC+, Harper et al. 2018) aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA) and the Submillimetre Common-User Bolometer Array 2 (SCUBA-2, Holland et al. 2013) camera mounted to the James Clark Maxwell Telescope (JCMT). The data characteristics and data reduction of all those new observations are described in more detail below.

2.3.1. DECam optical observations

DECam imaging in the u band was obtained for 14 CARS targets during the nights of 1−3 April 2017 as part of the program 2017A-0914 (P.I.: Grant Tremblay). DECam covers almost 3 sq. deg using a mosaic of 62 2 k × 4 k CCDs with a pixel size of 0.27″. The total exposure times for each target ranged from 1200 s to 4800 s. The data were automatically reduced by the DECam community pipeline (Valdes et al. 2014) and retrieved from the NOAO archive. Absolute flux calibration for individual images was achieved by comparing detected stars in the field with the photometrically-calibrated SkyMapper DR1 (Wolf et al. 2018) catalog magnitudes.

2.3.2. WFI optical observations

WFI imaging in the UBV bands was obtained for HE 0108−4743 and HE 2211−3903 to provide optical photometry given the missing Pan-STARRS and SDSS coverage. The observations were taken during a larger observing run from 18 to 28 October under program 0100.A-9003 (P.I.: Bernd Husemann). WFI covers a 34′×33′ FoV using a mosaic of 8 2 k × 4 k CCDs. At the time of observations 2 of the CCDs were broken so that a two-pointing offset scheme plus dithering was used to cover the same field. A total exposure time of 3600 s was obtained for the U band.

The images were fully reduced and combined with the latest version of the THELI pipeline10 (version 3) (Schirmer 2013). The absolute magnitude zero-point for each combined image was obtained by comparing stellar counts with catalog values from the SkyMapper DR1 and APASS DR1 (Henden et al. 2009) for the UBV bands.

2.3.3. PANIC NIR observations

Additional NIR imaging was also obtained with the Panoramic Near-Infrared Camera (PANIC, Cárdenas Vázquez et al. 2018) mounted to the 2.2 m telescope at the Calar Alto observatory. Although PANIC consists of 4 Hawaii-2RG detectors only 2 were operational at the time of observations of which only one had nominal performance. 20 CARS targets were observed as part of programmes F17-2.2-014 and H17-2.2-008 (P.I.: Bernd Husemann). For the first program only the best performing chip was used covering a FoV of 15′×15′. For the second program a mosaic strategy was used to cover the nominal 30′×30′ with the two operational detectors. The sampling of PANIC is per pixel. For the analysis in this paper we only use a small part of the images covering the AGN host galaxies so that the full FoV is not relevant.

Total exposure times ranged from 1−2 h for the J, H, and Ks bands, but not all NIR bands could be obtained for all the targets. A spiral dither pattern was used along all observations for background subtraction, excluding satellite tracks and image cosmetics. The PANIC exposures were reduced and combined with the dedicated PANIC pipeline11. The astrometric registration and absolute photometry was further optimized with the PhotometryPipline (PP, Mommert 2017) software package using 2MASS as photometric reference.

2.3.4. SCUBA-2 FIR observations

SCUBA-2 is a 10 000 pixel bolometer camera simultaneously operating at 450 μm and 850 μm and covering an area of 41 arcmin2 on the sky (Holland et al. 2013). Effective beam sizes are 10″ and 15″ at 450 μm and 850 μm, respectively.

Given the angular size of our AGN hosts on the sky, we observe using the ‘constant velocity (CV) Daisy’ mapping mode which provides uniform sensitivity in the central 3′ of the observation. The ‘CV Daisy’ is a circular scanning pattern designed so that the target is always within the field-of-view of the array throughout the integration while moving at a constant 155″ s−1. The observation provides usable coverage out to ∼6.0′ in radius, but beyond a radius of 1.5′ the map sensitivity decreases rapidly.

The observations were performed in service mode between February and December 2019 during grade 3 weather conditions (0.08 < τ225 < 0.12) as part of project M19AP019 (P.I.: Timothy Davis). On-source exposure times ranged from 1 h to 7.5 h. The resulting maps were reduced using the Dynamic Iterative Map Maker (DIMM) within the sub-mm user reduction facility (SMURF, Chapin et al. 2013). For a full overview of the procedure see Chapin et al. (2013).

The data were reduced using the blank field reduction in order to detect point sources within the map. The maps were then calibrated using the standard flux calibration factors of FCF450 = 491 Jy beam−1 pW−1 and FCF850 = 537 Jy beam−1 pW−1 (Dempsey et al. 2013). To improve the point source detection, a matched filter was applied to the maps. This matched filtering caused a ∼10% loss in flux (determined by inserting artificial sources into the map and comparing the flux before and after), so an additional 10% was applied to the FCFs to account for this.

2.3.5. SOFIA HAWC+ FIR observations

The HAWC+ FIR camera (Harper et al. 2018) aboard SOFIA was used to obtain FIR photometry for CARS targets that could not be detected with all-sky AKARI FIR survey and that were not targeted with Herschel. Observations were taken for 8 out of 11 proposed targets as part of a SOFIA survey program (Plan ID: 07_193, P.I.: Bernd Husemann) with exposure times ranging from 15 min to 60 min. HE 0227−0913 was observed in the E band (214 μm) while all other targets were observed in the D band (154 μm). The FWHM of the HAWC+ beam is and for the E and D band, respectively. All data were reduced by the HAWC+ SOFIA team using a dedicated instrument pipeline. The pipeline was run specifically with the faint object flag to account for the faint signal in our targets.

2.4. Integral-field spectroscopy observations

IFU observations for the CARS sample were mainly obtained with the MUSE instrument at the Very Large Telescope (VLT; Bacon et al. 2010) for 37 targets under programs 094.B-0345(A), 095.B-0015(A), 099.B-0242(A), and 099.B-0249(B), with the VIMOS instrument at the VLT (Le Févre et al. 2003) for 2 targets under program 083.B-0801(A), and with the PMAS instrument at the 3.5 m Calar Alto telescope (Roth et al. 2005) for 2 targets under program H18-3.5-010. MUSE covers a 1′×1′ FoV at a sampling of , VIMOS covers a 27″ × 27″ FoV at a sampling of and PMAS covers a 16″ × 16″ FoV with a sampling of 1″. All details of the IFU observations and their data reduction are presented in Husemann et al. (2022) to which the interested reader is referred.

For the present paper, we are working with emission line maps that were created after the AGN point-like emission was removed from the reconstructed IFU datacube as described in Husemann et al. (2022). Emission-line fluxes are extracted by fitting Gaussian profiles after modeling and subtracting the stellar continuum with PYPARADISE (see Walcher et al. 2015; Weaver et al. 2018; Husemann et al. 2022, for details) from the AGN-subtracted data. All the details of this process and its application to the CARS IFU data are presented in Husemann et al. (2022). For the analysis presented in this paper we use the 2D emission line maps of [O III] λ5007, Hβ, [N II] λ6583, Hα, [S II] λλ6716,6731 and their associated errors. All those data are also accessible from the CARS data release 112.

3. Analysis

3.1. modeling the panchromatic spectral energy distribution

In order to infer primary physical parameters such as the stellar mass M or the SFR, we need to model the panchromatic SED constructed for all the sources. We decided to perform the SED fitting with the publicly available AGNfitter package (Calistro Rivera et al. 2016). AGNfitter models the observed photometry as a super-position of various template libraries for the stellar continuum, the AGN accretion disk emission, AGN-heated dust emission from the torus and the cold dust emission excited by star formation.

The stellar library is taken from Bruzual & Charlot (2003) and parametrized by a stellar age of the galaxy and an exponential decay time of the starburst τ which can be further modified by reddening and normalized in absolute flux. The accretion disk emission is represented by an empirical composite spectrum generated from 259 type-1 AGN from the Sloan Digital Sky Survey (SDSS, York et al. 2000), which can be reddened using the Small Magellanic Cloud reddening law of Prevot et al. (1984). This leaves the flux normalization and the reddening as the only two free parameters for this component. The hot dust component of the AGN is reconstructed from the empirical SED library collected by Silva et al. (2004) based on photometric observations of AGN that were modeled with the radiative transfer code GRASIL Silva et al. (1998) to obtain the full NIR SEDs. The library is divided in a range of absorbing neutral hydrogen column densities NH through the torus along our line-of-sight as the main free parameter in addition to the normalization. Semi-empirical starburst template spectra from Chary & Elbaz (2001) and Dale & Helou (2002) are used to model the cold dust component from star formation which are denoted by a luminosity parameter and can be normalized.

AGNfitter uses the emcee package (Foreman-Mackey et al. 2013) for an affine-invariant Markov-chain Monte Carlo (MCMC) ensemble sampling to explore the 10-dimensional parameter space associated with the underlying template library of spectra. The advantage of AGNfitter is that it natively accepts upper limits in the photometry, allows to set priors on parameters and additional constraints can be incorporated easily to break degeneracies. Particularly, one big issue in the SED modeling of luminous type 1 AGN is that the accretion disk model is competing with the stellar population model in the UV–optical rest-frame wavelength regime. This leads to a degeneracy between the accretion disk and stellar continuum model when fitting the observed integrated broad-band photometry as shown in Fig. 2 (left panel).

thumbnail Fig. 2.

SED for HE 0114−0015 and the best-fit model determined by AGNfitter (Calistro Rivera et al. 2016). 20 MCMC realizations of the SED model are shown for the broadband photometric data (black points with error bars) where the red points are the corresponding predictions from the model. The red lines represent the total model, which consists of four components: the cold and warm dust in star-forming regions (green lines), the torus of AGN-heated dust (purple lines), the stellar continuum (yellow lines), and the AGN accretion disk (blue lines). On the right panel is an additional constraint to the accretion disk model (e.g. g and i AGN photometry is included).

This degeneracy could be solved by linking the UV radiation field with the corresponding re-emission in the IR as an additional constraint implemented for example in the CIGALE or MAGPHYS SED fitting codes. However, this method requires an assumption of the UV escape fractions and dust content of the galaxies, which can be very different for starburst galaxies and AGN-dominated galaxies. Due to the low redshift of the CARS sample, it is possible to spatially separate the AGN and host galaxy light using 2D image synthesis modeling even in ground-based images. This allows us to obtain direct measurements for the AGN brightness in the optical bands. The AGN magnitudes in the g and i bands are obtained by modeling the ground-based images as a super-position of two Sersic profiles and a point source for the AGN with the galfit package (Peng et al. 2002, 2010) as described in Husemann et al. (2022). Hence, we adjusted the likelihood function within the AGNfitter code to consider these additional photometric data as constraints solely for the AGN accretion disk component model. The effectiveness of this approach is shown in Fig. 2 in comparison to the modeling without those additional AGN constraints.

With the additional photometric constraints we reduced the uncertainties of the various physical quantities that are directly computed by AGNfitter from the posterior parameter distribution. For our work we are most interested in the stellar mass M that are based on the stellar population model and its associated luminosity normalization. In addition we report the integrated luminosity of the torus model (Ltor), and the luminosities of the cold dust component integrated between entire IR 8−1000 μm (L8 − 1000 μm) wavelength range and a restricted range of 42.5−122.5 μm (L42.5 − 122.5 μm) for literature comparisons. AGNfitter estimates a SFR from the integrated IR luminosity using the calibration established by Murphy et al. (2011):

(1)

In cases where the cold dust SED is not well constrained due to upper limits in the FIR photometry, we determine 5σ upper limits from the posterior distribution function of the integrated luminosities and SFR. We list all the inferred parameters from SED fitting as described above in Table 2 and all the corresponding SED models with AGNfitter are shown in Fig. A.1.

Table 2.

AGN host galaxy properties inferred from SED modeling.

3.2. Measuring Hα-based star formation rates

The IFU observations for the CARS sample as presented in Husemann et al. (2022) are well suited to measure the current star formation rate of the host galaxies in comparison to the FIR-based SFR. Various calibrations have been determined to convert the extinction-corrected Hα luminosity into a SFR from which we use the calibration of (Calzetti et al. 2007)

(2)

in this work to estimate the integrated SFR based on the integrated Hα. However, not all of the Hα emission is associated with star-forming H II regions in AGN host galaxies. A detailed spatially resolved emission line diagnostic analysis is therefore essential for AGN host galaxies to separate the contribution from AGN photoionization and star-forming H II regions to the excitation of Hα.

In order to provide reliable Hα-based SFR estimates, we use a combination of line maps from pixel-by-pixel and binned IFU cubes after subtracting the AGN point-source emission as described in Husemann et al. (2022) in more detail. While the pixel-by-pixel maps provide high-spatial resolution information for bright emission line regions, the binned maps provide higher signal-to-noise ratio (S/N) information on the more extended ionized gas. We only consider the diagnostic lines Hα, Hβ, [N II] and [O III] as detected above a 3σ level. For the maps we also include the weaker spaxels where only Hα S/N is required to be more than 3.

Here, we consider several ways to correct the Hα flux by the AGN contribution as described in more detail in the following subsections. All those methods leverage the potential of the IFU data to perform flexible spatially resolved emission-line diagnostics that open possibilities for de-contamination compared to traditional narrow-band imaging or long-slit spectroscopy.

3.2.1. Demarcation lines

The simplest approach to correct for the AGN affected Hα emission is to ignore those spaxels displaying AGN ionization and integrate only the Hα flux from regions clearly associated with H II regions. Specific BPT demarcation lines are often used as hard boundaries to distinguish between different ISM ionization processes in galaxies even though physically (due to line-of-sight blending and other physical processes) those demarcations cannot be hard boundaries.

The theoretical “maximum starburst” demarcation line proposed by Kewley et al. (2001) divides the region where ionization could still be explained by pure star-formation from the region where AGN ionization is required to contribute to the emission. Similarly, the empirical ‘pure star-forming’ BPT demarcation line proposed by Kauffmann et al. (2003) defines a part in the BPT assumed to be powered purely by H II regions whereas other mechanisms can play a role in the remaining part. The gap in between the stricter “pure star-forming” and the “maximum starburst” demarcation lines is often referred to as the mixed or composite BPT region.

The AGN contamination to Hα across the BPT and, in particular, through the composite region can be tackled in different ways. Nascimento et al. (2019) considers the pure star-forming region and the composite region to be entirely associated with star formation assuming that the AGN contribution in the composite region compensates for the ‘hidden’ star formation. Wylezalek et al. (2018) assigns a fixed fraction of 80% of AGN ionization contribution to Hα in the BPT region above the “maximum starburst” demarcation line and only a 20% contribution to the composite region. We use the former approach (illustrated in Fig. 3, left panel) to be able to consistently compare our results (Table 3) with previous works and the more sophisticated methods as described below.

thumbnail Fig. 3.

BPT diagram of HE 0853+0102 as an example for a mixing sequence between H II regions and AGN photoionization across the host galaxy. We highlight three ways of determining the SF fraction from the BPT diagram. Left panel: a simple demarcation line approach where everything above the Kewley et al. (2001) line is classified as pure AGN and below as pure star formation. Middle panel: simple linear combination with a χ2 fit using a fixed line vector as a reference point for the characteristic AGN and H II line ratios. The colored halos around the points represent the difference between the calculated SF fraction adopting two different different reference pairs (triangle and square symbols) to highlight the systematic uncertainties introduced by a certain choice. Right panel: MCMC fitting results of the star formation fraction adopting a large cloud of potential reference points for the AGN (red points) and H II region (magenta points) which are explored during the MCMC sampling.

Table 3.

Hα-based SFRs and metallicity estimated from the IFU emission lines.

3.2.2. BPT mixing sequence

In contrast to the fixed contribution assumed before, there is likely a smooth transition of ionization contributions across galaxies. AGN host galaxies therefore often show a mixing sequence which appears as an elongated structure in the BPT diagram. This mixing sequence (MS) usually exhibits a prominent dependency on the distance from the AGN. As proposed by Davies et al. (2014a,b, 2016, 2017), one can pick basis vectors of emission line ratios to characterize pure AGN and star-forming ionization and treat the composite data points of the MS as a linear combination of the basis vectors:

(3)

Here, the parameters fAGN and fSF represent the non-negative linear coefficient for the AGN and star-forming ionization fractions, correspondingly. An important additional basic constraint of fAGN + fSF = 1 ensures that flux is preserved. In principle more excitation mechanisms could be added to this equation, but this will lead to degeneracies in the basis vector as they are not fully orthogonal.

The emission line columns are normalized to either Hα or Hβ, which minimizes the impact of internal extinction within the host on the results. Here, we include the [S II]/Hα line ratio from the other classical BPT diagram and more line ratios could in principle be considered if S/N are high enough to be informative.

The idea behind the method is to pick emission-line basis vectors at the extreme ends of the mixing sequence, which should best reflect the assumption of two completely independent excitation mechanisms. The process of picking such emission-line ratio basis vectors allows considerable freedom and uncertainty. For example, the selected pair of basis vectors may not fully capture the different ISM conditions across the galaxies such as metallicity or ionization parameter. In the Fig. 3, middle panel, the example of two different bases pairs – triangle and square symbols – leads to a 10% difference in the resulting SF fraction.

For more meaningful results, the uncertainties caused by the choice of the basis vector need to be estimated and included into the final SFR uncertainty. With Monte Carlo Markov chain (MCMC) algorithms it is possible to use the Bayesian approach and also take into account the uncertainty of the underlying model. So, our solution to the basis-picking problem is to define a large number of basis ratio sets and let the MCMC sampling to choose basis vectors for a certain spaxel and, more importantly, take the uncertainty of this choice into account. We use the MCMC python package EMCEE to develop a python package called rainbow13 designed to handle the mixing sequence fitting. For each given spaxel it tries various combinations of basis vectors and determines the posterior probability distribution of the parameters through maximizing the likelihood function:

(4)

where y represents the vector of emission line ratios of the fitted spaxel, σy is the corresponding error vector, xAGN and xSF represent basis vectors, and the model is ymodel(xAGN, xSF, fSF) = (1 − fSF)xAGN + fSFxSF such as the vector Eq. (3). The likelihood function here is a Gaussian where the variance is underestimated by the fractional amount parameter ℱ, that represents an uncertainty of the model. The star-forming ionization fraction fSF and its uncertainty then is inferred from the probability distribution. An example of the MCMC fitting result from rainbow is shown in the right panel of Fig. 3.

Our criteria of picking SF emission-line ratio basis vectors from the mixing sequence are: first, SF bases should be below the “pure star-forming line” (Kauffmann et al. 2003); second, the signal-to-noise ratio of all the involved lines should be 3σ or more (> 3σ threshold is used to reduce the number of the basis vectors which speeds-up the fitting process); third, a certain radius may be specified to bound the basis vectors closer to the center and represent a narrower range of metallicities and ionization parameters. The criteria of picking AGN (or shock/evolved star, as described in Sect. 3.2.3) photoionization basis vector from the mixing sequence are the following: first, [O III]/Hβ is above a certain threshold, this threshold cannot be same for all galaxies, as it depends on the strength of the AGN, and, therefore, is set to a value between 0.5−0.8 (0.25−0.5 for LINERs) depending on a galaxy, so that only the top spaxels are selected; second, a certain distance from the AGN may be specified to separate the AGN basis in the central spatial region and extended AGN cloud. The adopted selection is shown on a galaxy-by-galaxy case in Fig. B.1 for the entire sample.

3.2.3. BPT morphology

The idea of treating the mixing sequence as a signature of physical mixing of H II regions and AGN ionized regions suggests that different populations of spaxels on a BPT diagram can be differentiated based on their properties. We name this approach BPT morphology, as the idea is to distinguish different shapes on the diagram (see Fig. 4):

  • (1)

    The star-forming cloud is an extended structure located under the “maximum starburst line” of Kewley et al. (2001) on a BPT diagram. This emission arises from across the star forming body of a galaxy and, therefore, contains emission from regions with different metallicities, densities, and ionization parameters. The variance of these ISM parameters results in the extensive shape that covers large area in the BPT diagram. Although it spreads out to the composite area between the “pure star-forming” and the “maximum starburst” demarcation lines on the BPT diagram, we assume that there is no need to clean this group of spaxels from the AGN contamination. See HE 0150−0344 (Fig. 4 middle row, left) as an example of a galaxy without AGN but with the spaxels spreading to the composite BPT area.

  • (2)

    The AGN-ionized cloud is a structure which spans the upper region of a BPT diagram and spreads to the left side toward lower [N II]/Hα line ratios. This corresponds to the extended narrow-line region or AGN-illuminated gas. Spatially, this region can reach from the AGN across the entire galaxy and even outside of the galaxy’s main body. This group of spaxels does not contribute to the galaxy’s SFR, although we can infer the upper limits considering hidden star formation as highlighted below.

  • (3)

    The mixing sequence is already described in the section above; it is an elongated structure on a BPT diagram which tends to be located in the central region of a galaxy. However, not all mixing sequences end in the AGN ionized area of the BPT diagram, where some are extending from the H II regions area to the LINER area likely due to different excitation mechanisms for the LINER emission. Nevertheless, the LINER mixing sequence can be treated in the same way as the AGN mixing sequence, following the same assumption of the two main excitation sources – star formation and LINER excitation. The model (Eq. (3)) will then have fLINER instead of fAGN with the same flux preservation constrain fLINER + fSF = 1.

thumbnail Fig. 4.

Schematic sketch of the BPT morphology populations. Upper row: on the right there is a representation of a galaxy with a central region in green, the rest of the galactic body in blue, and an AGN ionization cone in yellow. These regions populate specific areas on the BPT diagram: left diagram: star-forming cloud with two different sizes; middle diagram: AGN-ionized cloud with two sizes; right diagram: mixing sequence elongated toward the AGN area and a mixing sequence elongated toward the LINER area. Middle row: BPT diagram schematic examples for the SF cloud, AGN cloud, mixing sequence and LINER mixing sequence dominated galaxies. Lower row: BPT diagram and spatial map examples for SF cloud and AGN cloud combined with a mixing sequence in one galaxy.

An AGN host galaxy can have one of the described populations or even a combination of those (see the galaxy-by-galaxy cases in Appendix B). We separate star-forming clouds and AGN-ionized clouds with a simple yet efficient method: manually introducing fixed spatial radii for each individual galaxy (listed on the Table B.1) and separating the regions with the circles as on the bottom examples of Figs. 4 and B.1. After separating different populations we assume that the SF cloud and AGN-ionized cloud have 100% and 0% contribution to the star-formation rate, correspondingly. The mixing sequence is modeled, out to a certain radius, with RAINBOW as described in the previous section to compute SF fractions along the sequence.

However, we attempt to correct all populations for the AGN contribution (similarly as with the mixing sequence). Given that some small contribution to the Hα emission from star formation may be hidden in galaxies which display only an AGN cloud morphology, we estimate an upper limit for the SF contribution adopting AGN bases or SF bases from the other galaxies with a prominent mixing sequence and apply our rainbow analysis.

3.2.4. Gas-phase metallicity

Another important characteristic of the ionized gas is the metallicity. The gas-phase metallicity tracks the immediate enrichment history of the ISM due to the evolution of stars and their metal yields across the galaxy. Undisturbed disk galaxies typically have negative metallicity gradients (Sánchez et al. 2014) and their central metallicity correlates with the stellar mass of a galaxy (Tremonti et al. 2004; Kewley & Ellison 2008). As stellar evolution is a long-term process, outflows and inflows of gas on short timescales can significantly affect the observed metallicity distribution and therefore can be used as a key diagnostic to understand the origin and motions of gas on galactic scales. For example a flattening and dilution of gas-phase metallicity have been observed during galaxy mergers and interactions (e.g., Ellison et al. 2008; Kewley et al. 2010; Thorp et al. 2019), and in barred galaxies (Martin & Roy 1994; Sánchez et al. 2014) where low-metallicity gas from the outskirts is efficiently transported toward the center. On the other hand, gas outflows can enrich the circum-galactic medium with metals from the galaxy center which has been observationally confirmed from absorption line studies (e.g., Bordoloi et al. 2011; Tumlinson et al. 2011; Bouché et al. 2012; Nielsen et al. 2015; Schroetter et al. 2019).

Measuring gas-phase metallicities across AGN host galaxies is more complicated because all strong-line metallicity diagnostics are calibrated for H II regions. Their application to regions photoionized by the AGN is invalid in the majority of cases, but specific calibrations have been developed (e.g., Storchi-Bergmann et al. 1998). The metallicity measurements can in principle be restricted to H II regions ionized by star formation (as identified from BPT diagnostics) in AGN host galaxies (e.g., Husemann et al. 2014), but this approach greatly limits suitable targets and radial coverage. The [N II]/Hα (N2 index) line ratio is one of the prominent strong-line metallicity calibrators (Pettini & Pagel 2004; Marino et al. 2013) for H II regions. Photoionization models show that the N2 index should also trace the metallicity in AGN photoionized regions (e.g., Groves et al. 2006) with a different scale. That the situation is more complex was shown by Stern & Laor (2013) as they found a secondary dependence of the N2 index with AGN luminosity. At the same time they diskovered that the [N II]/[S II] (N2S2 index) is well recovering the mass–metallicity relation of galaxies for AGN hosts independent of AGN luminosity. Indeed, the N2S2 index is generally a good metallicity calibrator (Dopita et al. 2016) also for SF ionized H II regions, but not as widely used given that the diagnostic lines are significantly fainter.

In Husemann et al. (2019), we established a N2S2 index calibration based on the SDSS galaxy sample, and used it to map the gas-phase metallicity across the CARS AGN galaxy HE 1353−1917. The metallicity pattern revealed that the ENLR of the edge-on galaxy follows a very similar radial metallicity than the star forming disk. This let us conclude that the diffuse extra-planar gas was expelled by SN-driven winds rather than an outflow from the central AGN. Here, we apply the N2S2 index calibration established in Husemann et al. (2019) for all CARS targets to infer the metallicity gradient and absolute scale of the central metallicity. In order to directly use the metallicity determined for SDSS based on Tremonti et al. (2004) we determined an N2S2 = log ([N II]λ 6583/[S II]λλ 6716,6731) calibration for this SDSS metallicity scale in the similar way as described in Husemann et al. (2019):

(5)

Based on this calibration we estimated the oxygen abundance across the galaxies and reconstruct the radial metallicity gradient, as shown for the face-on disk galaxy HE2211−3903 in Fig. 5. For this galaxy a negative radial metallicity gradient is clearly recovered with a slope of αO/H = −0.055 dex with a zero-point central metallicity of 12 + log(O/H) = 8.732. As expected from the work of Stern & Laor (2013), the N2S2 index indeed recovers a matching absolute metallicites even in the case of strong AGN ionization without significant offset in metallicity as shown in Fig. 5. Hence, the 2D map of the estimated oxygen abundance does not show any features related to the varying ionization conditions throughout its disk. While the [N II] and [S II] emission lines may not exactly originate from the same location within H II region, due to the ionization structure of the nebulae on 100 pc scales (e.g., Sanders et al. 2020; Mannucci et al. 2021), our calibration based on SDSS spectra (3″ apertures) consistently include also part of the surrounding diffuse gas similar to our MUSE observation with typical resolutions of 07-10 covering 300–1000 pc depending on redshift. Because we are unable to isolate individual H II regions with CARS, any calibration based on isolated H II would be invalid.

thumbnail Fig. 5.

Example of the 2D and radial metallicity distribution for HE 2211−3903. The oxygen abundance was measured using the N2S2 index which is almost insensitive to the excitation. On the left panel the full 2D distribution is shown and on the right panel a radial projection in units of the effective radius taken from the single Sersic model. A linear fit is performed and shown as the red line.

This methodology therefore allows us an initial characterization of the metallicity for the entire sample, as listed in Table 3 and maps shown in Fig. C.1. Notably, we exclude a close star-forming companion in case of HE 0203−0031 and HE 1017−0305 from the radial fitting and fixed the center to the AGN position in all cases. As the absolute metallicity scale is strongly dependent on the underlying calibration we also transform our oxygen abundance estimate from the original O3N2 calibration scale as used in Husemann et al. (2019) to the one of Tremonti et al. (2004) so that we can more easily apply calibrations from the non-AGN SDSS sample as discussed below. Uncertainties are determined by bootstrapping; re-fitting the linear relation after randomly sub-selecting 80% of the data points.

4. Results and discussions

4.1. Finding SFR dependencies on AGN parameters

In order to investigate the potential impact of AGN on the total SFR of their host galaxies we collect the derived CARS host galaxies parameters such as AGN-corrected integrated Hα and IR luminosity and associated SFRs together with the stellar mass from our analysis presented above. From the SFRs we compute the specific SFR (sSFR) by dividing with the stellar mass, that is sSFR = SFR/M and compare it against the stellar mass as shown in the Fig. 6 for the Hα and IR-based SFRs, respectively.

thumbnail Fig. 6.

Comparison of the sSFR against stellar mass for Hα- and IR-based SFR tracers. The CARS objects are shown as red points with error bars, while xCOLD GASS (xCG) are in small gray dots, KINGFISH are in blue and Herschel Reference Survey (HRS) in green. Specific SFR is calculated from Hα luminosity (left panel) derived after BPT morphology analysis and RAINBOW fitting using the SFR formula SFRHα = 5.3 × 10−42 × LHα (Calzetti et al. 2007). IR specific SFR (right panel) derived from the AGNFITTER modeling using the Murphy et al. (2011) SFR formula. The stellar mass is taken from the AGNFITTER modeling parameters. The background on the left plot is the density map of SDSS galaxies (Brinchmann et al. 2004). Linear models to the star-forming main sequence are shown as black solid lines from Renzini & Peng (2015) and Shimizu et al. (2015), respectively. The best-fit linear relation from the xCG sample is shown as the blue dashed line in both cases. Lower panel: residuals around the adopted star-forming main sequence where the colored bands (red, blue, and purple) emphasize three bins used further in Fig. 8.

We compare the CARS host galaxies with the non-AGN galaxy population at low redshifts using the data from SDSS (Brinchmann et al. 2004; Abazajian et al. 2009), xCOLD GASS (Saintonge et al. 2017), KINGFISH (Skibba et al. 2011), and the Herschel Reference Survey (HRS, Smith et al. 2012; Boselli et al. 2015). Those reference samples highlight the position of the so-called star-forming main-sequence (SFMS) as determined by various authors and is indicated on Fig. 6, with blue dashed lines. These were derived from SDSS galaxies in the Hα diagram by Renzini & Peng (2015) and Shimizu et al. (2015) for the IR diagram.

Those SFMS relations are derived from different data sets, and as our goal is to achieve consistency between the two diagrams, we consistently determine the SFMS from the xCOLD GASS (xCG) sample. The xCG is a subsample from the SDSS sample, therefore Hα, IR, stellar mass, and also CO(1−0) and metallicity data are available from the literature sources listed above. As there is also a BPT classification of the objects, we exclude AGN and LINER galaxies to define a clean star-forming subsample of galaxies (hereafter training sample) resulting in 197 objects for the Hα diagram and 86 objects for the IR diagram. The linear SFMS fits to the training sample data and the corresponding R2 scores (defined below in Eq. (7)) are:

(6)

Our SFMS linear relations are consistent with the literature determinations shown in Fig. 6 within the reported scatter.

In order to assess the differences between the CARS hosts and the training sample we compute the residuals along the SFMS relation (Fig. 6 bottom panel). The CARS objects exhibit a scatter around the SFMS of 0.69 dex based on Hα and 0.60 dex based on IR SFR, respectively, ignoring the upper limits. Moreover, the residuals reveal a clear negative correlation with a Pearson correlation coefficient of r = −0.69 and r = −0.61 for the Hα and IR SFR, respectively. This negative trend might be naïvely interpreted as a result of star-formation quenching where higher mass galaxies appear to be more passive potentially due to the effects of the AGN (e.g., Cicone et al. 2014; Saintonge et al. 2017; Lacerda et al. 2020) considering our AGN selection.

However, the stellar mass is not the only parameter which controls the SFR of a galaxy. The cold gas content was shown to be another fundamental parameter linked to the SFR and explaining part of the scatter in SFR perpendicular to the SFMS (Tacconi et al. 2018; Colombo et al. 2020; Piotrowska et al. 2020; Ellison et al. 2020). Furthermore, the gas-phase metallicity is being discussed to be linked to the SFR and stellar mass, the so-called fundamental mass–metallicity relation (Lequeux et al. 1979; Mannucci et al. 2010; Yates et al. 2012; Curti et al. 2020), but see also Sánchez et al. (2017, 2019). In particular, the metallicity could be related to the metallicity-dependent conversion of LCO to the total cold gas mass (e.g., Genzel et al. 2012; Bolatto et al. 2013; Carleton et al. 2017; Utomo et al. 2017) or the general ability of the gas to cool to form star-forming gas clouds.

In order to quantitatively compare the previous model with the ones introduced further we use the coefficient of determination or the R2 score (Draper 1998). It provides an indication of how well the fit is for the given dataset with the best possible score of 1.0. Unsuitable models can result in negative R2 scores and an R2 score of 0.0 would result from a constant model that always predicts the expected value. For a dataset y1yn with fitted values y1  modelyn  model and a mean value the R2 score is calculated as follows:

(7)

As our model above does not take into account the cold gas content nor the gas-phase metallicity of the host galaxy, we expand our one-dimensional linear model into a multidimensional linear model for the SFR (hereafter SFMS+gas model) with stellar mass (M), CO(1−0) luminosity (LCO) and metallicity (12 + log(O/H)) as the independent parameters. Here, we use the CO(1−0) luminosity from Saintonge et al. (2017) as the main cold gas mass proxy and the oxygen abundance 12 + log(O/H) determined from the ionized gas emission lines in the SDSS spectra as determined by Tremonti et al. (2004) as the main gas-phase metallicity parameter. This leads to the following formulas for the linear multidimensional model

(8)

(9)

with R2 scores of R2 = 0.64 and R2 = 0.89, respectively.

As gas content is by far the most dominant driver compared to metallicity and metallicity may not be always measurable, we also determine relations that only rely on stellar mass and gas content:

(10)

(11)

The quality of these models can be represented by R2 scores of R2 = 0.64 and R2 = 0.84, respectively. As expected the quality of the model does not become significantly worse and can be securely used for all samples where metallicity is not available. Here, we use the full model including metallicity because this parameter is available to us and so we can use it in our analysis. In any case the choice of the models does not alter our final results given the similarity of both models and low impact of metallicity.

After building the new model to predict SFR using the xCG training sample, we apply it to the CARS objects with the additional information of LCO(1 − 0) taken from Bertram et al. (2007) and 12 + log(O/H) inferred from the CARS IFU observations as described above. With the new SFMS+gas model we obtain new predictions for the expected SFR of individual CARS objects and present the corresponding residuals in Fig. 7. In comparison with the previous residuals (Fig. 6 bottom panel), the residuals are significantly smaller. The scatter is reduced by a factor of two as the trend with the gas content is taken into account. The initially observed correlation with stellar mass is significantly flattened, with a correlation coefficient consistent with no remaining mass dependence. One potential issue may arise for AGN hosts as the abundance of cold gas has been suspected to be linked with the AGN either by expelling or heating the cold gas (e.g., Carniani et al. 2017; Kakkad et al. 2017; Perna et al. 2018; Bischetti et al. 2021; Circosta et al. 2021). Considering that the residual sSFR scatter is significantly reduced when incorporating cold gas mass in the SFR prediction indicates that this potential effect is negligible for our sample.

thumbnail Fig. 7.

Difference between measured and expected Hα (upper panels) and FIR luminosity (lower panels) as a function of stellar mass, CO(1−0) luminosity and 12 + log(O/H) ionized gas metallicity. Here we used our more complex model (Eqs. (1) and (2)) to predict the SFR-sensitive luminosities. The black dots are the data from our xCG training sample and the red symbols are the CARS data with triangles indicating upper limits. The colored bands (red, blue and purple) define three bins in residual luminosity as discussed in the text and used for Fig. 8.

In order to understand if there are any remaining link between the residuals in sSFR and various AGN parameters we divide the residual space into three bins (negative, central, and positive residual bins) and compare the mean values of the AGN parameters for each bin. The bins are created such, that the number of objects in the bins is approximately equal (with a small excess of objects in the central bin). Objects with upper limits are treated separately as we cannot necessarily associate them to a specific bin. The central residual bin for IR luminosity is shifted slightly toward higher values because of the larger number of upper limits that are too high to be informative. Those objects with upper limits are missing to populate the central and negative residuals bins, but are mostly included in the case of Hα. As we are only interested in the relative changes along the residuals we focussed to have sufficient objects in each bin for statistical reasons. In Fig. 8, we compare bins in residual SFR for the SFMS and the new SFMS+gas models for both Hα and IR SFR tracers allowing us to explore trends with BH mass, AGN luminosity, Eddington ratio and b/a host galaxy axis ratio as reported by Husemann et al. (2022) and Singha et al. (2022) for both proxies.

thumbnail Fig. 8.

Comparison of SFR residuals from the SFMS and SFMS+gas model against various AGN parameters. Upper left panel: histograms of the residuals for SFMS and SFMS+gas for both Hα and IR models, with the red, blue, and purple colors corresponding to the previous plots (Figs. 6 and 7, background colors). The other panels show the comparison of the AGN parameters taken from Husemann et al. (2022) and Singha et al. (2022): logarithm of BH mass (upper center), logarithm of AGN bolometric luminosity (upper right), number of the objects, containing extended outflows (lower left), logarithm of Eddington ratio (lower center), and b/a morphological parameter (lower right). Each of the four models is represented with three bins (from left to right: negative red, central blue, positive purple) and the forth bin with the upper limits (gray). The error bars of the bins also represent the distribution of the parameters within the bins. The black dots and lines highlight the trend of the mean values that are written over the bins.

We find no trends of the residuals to systematically change with the log MBH and log LAGN given the uncertainties. While a putative correlation is seen with the fraction of objects containing extended outflows in the SFMS+gas model using the Hα proxy, it is not reflected in the IR tracer and therefore not robust. While the logarithm of the Eddington ratio (log λEdd) has a slight increasing trend (e.g. higher AGN Eddington ratio in SFR excess systems) for the SFMS model, such a trend vanishes for the SFMS+gas model. It is important to note that most of the previous studies that investigated the impact of AGN on the SFR only used stellar mass as a proxy for the expected SFR (e.g., Page et al. 2012; Husemann et al. 2014; Shimizu et al. 2015; Balmaverde et al. 2016; Bernhard et al. 2016; Catalán-Torrecilla et al. 2017; Circosta et al. 2018; Scholtz et al. 2020). The inclusion of the cold gas mass to predict the expected SFR in addition to the stellar mass can significantly alter the derived conclusions about potential positive, negative or no AGN feedback on their host galaxies. It is therefore crucial for future studies to take more parameters into account than just stellar mass to properly characterize the parent population of non-AGN galaxies.

In Fig. 8 the geometrical parameter b/a (tracing the galaxy inclination) shows a clear trend for both Hα and IR SFRs for the SFR model. When comparing the b/a distribution for the negative and positive residual bins with Kolmogorov–Smirnov and Anderson–Darling tests we find a highly significant difference between them for the IR SFR with a p-value of 0.003, but less significant for Hα. This suggests that dust extinction or hidden star formation is not the driver for this trend. The trend significantly flattens for the SFMS+gas model considering that the overall scatter in the residuals is greatly reduced, so that the difference become insignificant for our sample size. A potential explanation for the tentative trend might be the orientation of the AGN. The CARS objects are unobscured AGN hosts, which means that the ionization cone of the AGN is pointed to the observer. Hence, the cross-section of the AGN ionization cone which can interact with the ISM is increased when the host galaxies are closer to an edge-on geometry. Such misalignments of the galaxies rotation axis and the central AGN engine can have an important influence on the impact of the AGN on the host galaxy, as has already been proposed for various individual galaxies (e.g., García-Burillo et al. 2014; Gallimore et al. 2016; Cresci et al. 2015b; Alonso-Herrero et al. 2018; Husemann et al. 2019; Smirnova-Pinchukova et al. 2019). With the CARS sample we start to see a potential systematic trend across the population where the potential of the AGN to suppress SFR might also dependent on the relative orientation of the central engine with respect to the galaxy. However, our CARS sample is strongly limited by the low number of more edge-on systems and a bigger sample in needed to gain more insight into this process. Our results highlight that the AGN luminosity may not be the only factor determining the ability of the AGN to impact star formation. The efficiency with which the released energy can couple to the host galaxy adds complexity to the picture, where geometry is certainly just one of many additional parameters to be considered.

4.2. Comparison of SFR tracers for AGN host galaxies

Another complexity in investigating the impact of AGN on star formation is caused by the fact that different SFR indicators trace different timescales of star formation. Hα emission is related to the most recent star formation (∼5 Myr) as excited necessarily by hot O stars, whereas IR traces the dust heated by a wider range of stars, and, therefore traces SF over longer timescales (∼100 Myr; e.g. Kennicutt 1998b; Hayward et al. 2014; Flores Velázquez et al. 2021). Both SFR tracers are expected to be implicitly linked to the CO luminosity as a molecular gas mass tracer, because the cold gas is the necessary seed reservoir of which new stars can be formed.

By combining the SFR calibration for the IR and Hα luminosity (Eqs. (1) and (2)) we can make a prediction for their expected relation. Furthermore, we can link the SFR to the CO luminosity given the calibration of Lin et al. (2019). This leads to the following relation between the Hα luminosity and the IR and CO luminosity, respectively.

(12)

(13)

Comparing the three luminosities with each other in Fig. 9 we find that the CARS sample and the training sample from the xCG are significantly offset with respect to the predicted relation between Hα and IR luminosity. Despite this, the relation between LCO and LHα as predicted from the combined Calzetti et al. (2007) and Lin et al. (2019) calibrations agree well with our data. This points to the notion that either the theoretical assumption made by Murphy et al. (2011) is inaccurate or that some of the ongoing star formation is completely obscured in Hα which cannot be recovered with an extinction correction. Nevertheless, the xCG non-AGN sample can be used as a reference because the impact of extinction on the FIR/Hα ratio should be comparable to our AGN host galaxies. We therefore performed a linear fit directly on the xCG data and obtained the following relation which are shown as the black lines in Fig. 9:

(14)

(15)

thumbnail Fig. 9.

8−1000 μm IR luminosity LIR (left panel) and the CO (1−0) luminosity LCO (right panel) plotted against Hα luminosity. The training sample from the xCG (gray dots) is compared to the CARS objects (red circles and triangles for the 5σ upper limits). Galaxy mergers in the xCG sample as classified by GalaxyZoo2 are additionally marked (black crosses) on the left panel. The blue dashed line corresponds to Eqs. (12) and (13) are derived by the combination of literature calibrations as described in the main text. The black lines represent a linear relation fitted to the XCOLD GASS reference sample, see Eqs. (14) and (15). The red dotted line is fitted to the CARS objects with the same slope as the black line, therefore highlight the shift in the normalization.

Even though the CARS sample scatter is close to the comparison sample scatter, the CARS galaxies are on average located ∼0.2 dex above the fitted LIR to LHα trend. The question here is whether this offset has a physical meaning or is just caused by low-number statistics and unknown biases of the CARS sample. Here, we claim that this offset is physical, based on two sets of evidence: (1) the CARS data are in a good agreement with the non-AGN reference sample on the LCO to LHα plot; (2) the training sample also has a few galaxies, located at similarly high LIR/LHα ratios which drive the offset in the CARS sample. The xCG objects in this specific area of the LIR/LHα diagram are almost exclusively populated with galaxy mergers according to the GalaxyZoo2 classification (Willett et al. 2013; Hart et al. 2016). Those mergers seem to have a lower instantaneous SFRHα compared to the slightly longer timescale IR SFR, which can be expected from the rapid star formation history evolution and bursts of star formation on 100 Myr timescale that are observed and predicted by detailed galaxy simulation (e.g., Mihos et al. 1992; Barnes 2004; Springel et al. 2005; Di Matteo et al. 2008; Hopkins et al. 2013). Considering that the offset for the CARS sample is mainly caused by individual galaxies with similarly high LIR/LHα ratio, as in the non-AGN merger sample, some AGN host galaxies potentially had higher SFR in the past (over a ∼100 Myr timescale).

Whether the potential difference in the star formation history is caused by a recent excess of star formation, which might actually be linked to a delayed BH growth (e.g., Wild et al. 2010), or by ongoing suppression of star formation as expected from AGN feedback remains unclear. To shed light on these different possibilities we discuss the difference between current and recent star formation as probed by Hα and FIR emission in the following for individual objects where more information is available and can be interpreted.

4.3. Individual comparison of current and recent SFRs

In the previous subsection we compared different SFR tracers for the CARS sample with respect to the non-AGN comparison sample. Now we look at the SFR difference derived from the different SFR proxies for individual CARS galaxies and try to discuss the potential origins. We visualize the SFR results together with host galaxy morphologies and an indication for the AGN variability in Fig. 10. The AGN variability is set here by the brightness difference in the NIR photometry between old 2MASS and more recent VISTA or PANIC observations separated by a few years. The Hα-, IR- and CO-based SFR are again defined according to Eqs. (2) and (1), respectively. We plot objects with upper limits in either of the SFR traces separately as any difference between the SFR cannot be inferred on an individual basis.

thumbnail Fig. 10.

Comparison of different SFR measurements for individual objects. Left column: Hα-, IR-, and CO(1−0)-based SFRs for the CARS objects sorted by the SFRIR − SFRHα difference. Central column: near infrared variability; the objects with no AGN or known changing look AGN are also marked here. Right column: morphology of the host galaxies. Lower panel: contains the objects with either Hα and/or IR SFR uncertainties (upper limits or missing data).

The sign of the relative difference of Hα- and IR-based SFR is highlighted in green and blue color when the IR-based SFR is higher or lower compared to Hα, respectively. Assuming that the SFR timescale is the prime physical explanation for the difference (e.g., Hayward et al. 2014; Davies et al. 2015), where Hα probes more recent star formation than the IR, we can in principle differentiate between declining or enhanced star formation. In Fig. 10 we order the objects from the strongest decline to the highest enhancement in SFR. The objects HE 0253−1641, HE 0119−0118, HE 1126−0407, HE 0045−2145, HE 2233+0124 and HE 1029−1831 belong to the declining SFR group and HE 0150−0340, HE 0212−0059, HE 0853+0102, HE 0108−4743, HE 0227−0913, HE 1338−1423, and HE 1310−1051 belong to the galaxies with potentially enhanced star formation. Below we discuss different scenarios for those individual cases that may provide physical explanations, but certainly require further tests and observations to verify.

All of these galaxies are disk-dominated and none are bulge-dominated or strongly interacting systems. Strong interactions therefore seem to play a minor role, but we still identify cases where this might be important but not obvious. In addition, bars might suppress the star formation dynamically due to bar quenching (e.g., Khoperskov et al. 2018; Fraser-McKelvie et al. 2020), but we identified nonstar forming and star forming bars in the decreasing SF group, based on the analysis of Neumann et al. (2019). While the majority of the galaxies with declining SFR are indeed barred galaxies, most of them appear star-forming contrary to expectation if they would suppress star formation. AGN or starburst-driven winds may play a role in individual cases of HE 1126−0407 and HE 0045−2145, which have both declining SFR. HE 1126−0407 (aka PG 1126−041) is known to have a powerful ultra-fast outflow driven by the AGN (Wang et al. 1999; Giustini et al. 2011) which may couple more efficiently on kpc scales (Marasco et al. 2020) as the galaxy is significantly inclined with respect to our line-of-sight and central AGN engine orientation. Furthermore, HE 0045−2145 was misclassified in the Hamburg/ESO survey because of the broader lines caused by a starburst-driven outflow (Nevin et al., in prep.). In both cases the outflow might be related to the anticipated decline in SFR. The targets HE 0253−1641 and HE 0119−0118 show high gas dispersion in the AGN-ionized region on kpc scales (see Fig. B.1 in Husemann et al. 2022) which also points to a past or ongoing galactic outflow event in these systems. Such outflow signatures are much weaker in HE 2233+0124, but the galaxy is also more edge on with a misaligned ionization cone similar to HE 1126−0407 so that the impact of the AGN may be amplified. The cause of the potential decline in the SFR of HE 1029−1831 is less clear, but the stellar population modeling in high-angular NIR IFU spectroscopy with SINFONI revealed a recent circumnuclear starburst about 100 Myr ago that is rapidly declining (Busch et al. 2015). Although the cause of this decline is hard to directly link to AGN feedback or simple gas consumption, it clearly supports the timescale interpretation of the difference in the Hα and FIR-based SFR.

For the objects with increased SFR, we identified HE 1310−1051 to exhibit a strong interaction with a minor companion (Husemann et al. 2014) and HE 0150−0344 to be a strongly interacting non-AGN system. In those two cases, the interactions are most likely responsible for a recent burst of star formation. We cannot pin-point any obvious reasons why the star formation may be enhanced in HE 0212−0059, HE 0853+0102, HE 0108−4743, HE 0227−0913 and HE 1338−1423, but we note that the bolometric luminosity is close to the Eddington luminosity for the last three sources as reported in Husemann et al. (2022) due to their rather low BH mass with respect to the AGN luminosity. This is consistent with the scenario discussed in Husemann et al. (2022) that low BH mass AGN are likely to be observed in earlier phase of the AGN cycle. In this case, the circumnuclear starburst would be observed much closer than 100 Myr in time with respect to its peak activity. Indeed, narrow-line Seyfert 1 galaxies as high-Eddington ratio AGN show systematically higher SFR based on PAH emission detected with Spitzer (Sani et al. 2010) compared to the broad-line Seyfert 1 counterparts with lower Eddington ratios. This highlights the intriguing connection between the AGN and the circumnuclear starburst, where AGN with increasing SFR are potentially young AGN in a fueling mode powered by the starburst.

Lastly, we obtained some information on AGN variability as a side-product of our analysis when we compared old 2MASS NIR photometry with the more recent VISTA or PANIC observations. Interestingly, we see that the targets with declining SFR show systematically less variability on a few tens of years timescale. It is still speculative whether this points to a much more stable energy output of the AGN on longer timescales and thereby enhancing the impact as a cumulative effect of energy release over time. However, it shows that the process of AGN feedback is complex and the time evolution of the AGN phase with respect to the galaxy needs to be considered as well to get a comprehensive picture. Indeed, in current self-regulated feedback models (Gaspari et al. 2020 for a review) the AGN is expected to flicker on-off with rapid variability that increases toward low-mass systems often due to the chaotic cold accretion feeding the SMBH (e.g., Gaspari et al. 2015; Tremblay et al. 2018; Rose et al. 2019).

5. Summary and conclusions

In this paper, we presented a complete census of the integrated SFR properties across the entire CARS sample using the entire multiwavelength data set available for this local AGN host galaxy sample. In particular, we inferred robust stellar masses and IR-based SFR from panchromatic SED modeling with AGNFITTER and Hα-based SFR from a careful analysis of spatially resolved optical emission-line diagnostics with our new rainbow algorithm. Using the large xCG sample of local non-AGN galaxies as a control and training set we built different models to predict the SFR for our AGN host galaxies to investigate the role of the AGN in terms of star formation feedback. Our main conclusions from this analysis can be summarized as follows:

  • We find that stellar mass alone is an insufficient proxy for the expected SFR in AGN host galaxies. The cold gas content and possibly the metallicity are crucial to consider in order to avoid artificial trends with AGN parameters that could mimic expected AGN feedback trends.

  • No systematic suppression of SFR could be detected with respect to the non-AGN galaxy reference sample and there is also no trend with AGN luminosity.

  • A potential link between lower than expected SFRs and the axis ratio b/a of the AGN host galaxies was identified. As our sample contains only unobscured AGN, the central engine (and thus the ionization cones) of low b/a systems must be mis-aligned with the galaxy rotation axis for inclined disks, leading to a higher cross-section of the AGN radiation field which can interact with the galaxy disk. This tentative trend is much weaker and currently insignificant for the SFR+gas model given the low-number of strongly inclined systems, but it may imply that such mis-alignments could amplify the coupling of the released AGN energy with the cold gas disk of the galaxy and thereby impact star formation more efficiently. This should be confirmed with bigger AGN host galaxies samples.

  • Interpreting the IR and Hα SFR tracers as proxies for the recent star formation history on ∼100 Myr and ∼5 Myr timescales respectively, we identified systems with decreasing or increasing SFR. The declining SFR cases might often be associated with galactic outflows while the increasing SFR cases can be associated with interaction or potentially with a young AGN phase.

That we cannot find any strong evidence for a global positive or negative AGN feedback on the SFR across the entire CARS sample is in agreement with various recent studies reporting no immediate impact of AGN on the star formation (e.g., Scholtz et al. 2020). We emphasize that predicting the expected SFR of galaxies is difficult, and using stellar mass alone may not be sufficient. The cold gas mass is fundamental, and its use in such relations allows one to capture the significant scatter in the star-forming main sequence (as already demonstrated in several works for the non-AGN population; e.g., Colombo et al. 2020; Popesso et al. 2020). Despite the lack of a global and obvious impact of AGN on star formation, we diskover subtle effects that should be investigated in the future. Most importantly, the relative orientation of the AGN central engine and associated ionization cones may be relevant for the cross-section of release AGN energy and the cold gas of the galaxy. This is most prominent in disk galaxies and previously studied in several individual cases (e.g., Cecil et al. 2001; Morganti et al. 2015; Mahony et al. 2016; Mukherjee et al. 2018; Husemann et al. 2019), but CARS reveals a potential systematic trend that should be explored with larger samples in the future.

The nondetection of a relation of AGN luminosity, BH mass and Eddington ratio with the global SFR may be related to the different time scales of the AGN phase and star formation in galaxies. In case the AGN phase is short, there would not be enough time passed to see the impact on the global star formation when selecting AGN samples rather than post-starburst system as discussed in the review of Alexander & Hickox (2012) and in Hickox et al. (2014). Indeed, the CARS sample suggest a potential correlation with the duration of a luminous AGN phase as a function of BH mass (Husemann et al. 2022), which is further corroborated by current models of AGN feedback self-regulated via chaotic cold accretion (e.g., Gaspari et al. 2020). The predicted durations are on the order of 1 Myr for a single AGN phase which would be clearly too short to suppress the star formation in the entire host galaxy and can explain our observations. Still, the circumnuclear SFR could be affected on these timescales. We observed interesting patterns of increasing and decreasing SFR by comparing the IR and Hα SFR tracers among the sample that we can partially attribute to outflows, circumnuclear star formation and galaxy interaction. This highlights the potential of this approach, and also the complexity in the galaxy properties to be considered. We plan to expand the diagnostics of the star formation history determination by inferring radio-based SFR for the CARS sample which probe intermediate SFR timescale of a few tens of Myr, filling the gap in the Hα and IR-based SFRs. Furthermore, we expand our SFR investigation by zooming into the circumnuclear region of the galaxies in the CARS sample in the future. This requires the construction of an appropriate control sample with similar resolution which was beyond the scope of this paper.

Overall, we identify cold gas content, relative AGN engine orientation with respect to the host galaxy, as well as the time domain variability as potential key parameters that need to be explored in the future to understand the impact of AGN on their galaxy on a population wide basis. This leads to obvious challenges in the sample selection, sample size and parameter space to be measured to gain more insights into the putative AGN feedback process.


6

Available at https://dr12.sdss.org/

13

Publicly available open-source at https://gitlab.com/SPIrina/rainbow

Acknowledgments

We thank the anonymous referee for helpful comments that improved the quality of the manuscript. We thank Mischa Schirmer for his support with the usage of the updated Theli pipeline v3. I.S.P. greatly appreciates financial support from the DLR through grant 50OR2006. B.H. is financially supported through DFG grant GE625/17-1. I.S.P. and B.H. also acknowledge travel support from the DAAD via grant 57509925. T.A.D. acknowledges support from the UK Science and Technology Facilities Council through grant ST/S00033X/1. M.G. acknowledges partial support by NASA Chandra GO8-19104X/GO9-20114X and HST GO-15890.020/023-A, and the BlackHoleWeather program. S.B., C.O., and M.S. acknowledge support from the Natural Sciences and Engineering Research Council (NSERC) of Canada. M.P.T. acknowledges financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709) and through grants PGC2018-098915-B-C21 and PID2020-117404GB-C21 (MCI/AEI/FEDER, UE). The work of T.C. was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA. Based on observations collected at the European Organization for Astronomical Research in the Southern Hemisphere under ESO programme 094.B-0345(A) and 095.B-0015(A). Based on observations collected at the Centro Astronómico Hispano-Alemán (CAHA) at Calar Alto, operated jointly by Junta de Andalucía and Consejo Superior de Investigaciones Científicas (IAA-CSIC). This project used data obtained with the Dark Energy Camera (DECam), which was constructed by the Dark Energy Survey (DES) Collaboration. Funding for the DES Projects has been provided by the US Department of Energy, the US National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute for Cosmological Physics at the University of Chicago, Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Enérgeticas, Medioambientales y Tecnológicas–Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciències de l’Espai (IEEC/CSIC), the Institut de Física d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, NSF’s NOIRLab, the University of Nottingham, the Ohio State University, the OzDES Membership Consortium, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University. Based on observations at Cerro Tololo Inter-American Observatory, NSF’s NOIRLab (NOIRLab Prop. ID 2017A-0914; PI: Grant Tremblay), which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. The Pan-STARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation. The James Clerk Maxwell Telescope is operated by the East Asian Observatory on behalf of The National Astronomical Observatory of Japan; Academia Sinica Institute of Astronomy and Astrophysics; the Korea Astronomy and Space Science Institute; Center for Astronomical Mega-Science (as well as the National Key R&D Program of China with No. 2017YFA0402700). Additional funding support is provided by the Science and Technology Facilities Council of the United Kingdom and participating universities and organizations in the United Kingdom and Canada. 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. Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA. 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. This research has made use of the NASA/IPAC Infrared Science Archive, which is funded by the National Aeronautics and Space Administration and operated by the California Institute of Technology. This research has made use of the APASS database, located at the AAVSO website. Funding for APASS has been provided by the Robert Martin Ayers Sciences Fund. Supported by the international Gemini Observatory, a program of NSF’s NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation, on behalf of the Gemini partnership of Argentina, Brazil, Canada, Chile, the Republic of Korea, and the United States of America. The Science, Technology and Facilities Council is acknowledged by JN for support through the Consolidated Grant Cosmology and Astrophysics at Portsmouth, ST/S000550/1.

References

  1. Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543 [Google Scholar]
  2. Aird, J., Coil, A. L., Moustakas, J., et al. 2012, ApJ, 746, 90 [CrossRef] [Google Scholar]
  3. Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12 [Google Scholar]
  4. Alexander, D. M., & Hickox, R. C. 2012, New Astron. Rev., 56, 93 [Google Scholar]
  5. Alonso-Herrero, A., Pereira-Santaella, M., García-Burillo, S., et al. 2018, ApJ, 859, 144 [Google Scholar]
  6. Baade, D., Meisenheimer, K., Iwert, O., et al. 1999, The Messenger, 95, 15 [NASA ADS] [Google Scholar]
  7. Bacon, R., Accardo, M., Adjali, L., et al. 2010, Proc. SPIE, 7735, 773508 [Google Scholar]
  8. Baldry, I. K., Glazebrook, K., Brinkmann, J., et al. 2004, ApJ, 600, 681 [Google Scholar]
  9. Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 [Google Scholar]
  10. Balmaverde, B., Marconi, A., Brusa, M., et al. 2016, A&A, 585, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Balogh, M. L., Baldry, I. K., Nichol, R., et al. 2004, ApJ, 615, L101 [Google Scholar]
  12. Barnes, J. E. 2004, MNRAS, 350, 798 [NASA ADS] [CrossRef] [Google Scholar]
  13. Barth, A. J., Sarzi, M., Rix, H.-W., et al. 2001, ApJ, 555, 685 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bassini, L., Rasia, E., Borgani, S., et al. 2019, A&A, 630, A144 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Bernhard, E., Mullaney, J. R., Daddi, E., Ciesla, L., & Schreiber, C. 2016, MNRAS, 460, 902 [NASA ADS] [CrossRef] [Google Scholar]
  16. Bertram, T., Eckart, A., Fischer, S., et al. 2007, A&A, 470, 571 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Bianchi, L., Shiao, B., & Thilker, D. 2017, ApJS, 230, 24 [Google Scholar]
  18. Bing, L., Shi, Y., Chen, Y., et al. 2019, MNRAS, 482, 194 [NASA ADS] [CrossRef] [Google Scholar]
  19. Bischetti, M., Feruglio, C., Piconcelli, E., et al. 2021, A&A, 645, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Bluck, A. F. L., Maiolino, R., Sánchez, S. F., et al. 2020, MNRAS, 492, 96 [Google Scholar]
  21. Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207 [CrossRef] [Google Scholar]
  22. Bongiorno, A., Merloni, A., Brusa, M., et al. 2012, MNRAS, 427, 3103 [Google Scholar]
  23. Bordoloi, R., Lilly, S. J., Knobel, C., et al. 2011, ApJ, 743, 10 [CrossRef] [Google Scholar]
  24. Boselli, A., Fossati, M., Gavazzi, G., et al. 2015, A&A, 579, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Bouché, N., Hohensee, W., Vargas, R., et al. 2012, MNRAS, 426, 801 [CrossRef] [Google Scholar]
  26. Bower, R. G., Benson, A. J., Malbon, R., et al. 2006, MNRAS, 370, 645 [Google Scholar]
  27. Bradley, L., Sipőcz, B., Robitaille, T., et al. 2020, https://doi.org/10.5281/zenodo.4044744 [Google Scholar]
  28. Brammer, G. B., Whitaker, K. E., van Dokkum, P. G., et al. 2009, ApJ, 706, L173 [Google Scholar]
  29. Brinchmann, J., Charlot, S., Heckman, T. M., et al. 2004, ArXiv e-prints [arXiv:astro-ph/0406220] [Google Scholar]
  30. Brownson, S., Belfiore, F., Maiolino, R., Lin, L., & Carniani, S. 2020, MNRAS, 498, L66 [NASA ADS] [CrossRef] [Google Scholar]
  31. Brüggen, M., & Kaiser, C. R. 2002, Nature, 418, 301 [CrossRef] [Google Scholar]
  32. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  33. Busch, G., Smajić, S., Scharwächter, J., et al. 2015, A&A, 575, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Calistro Rivera, G., Lusso, E., Hennawi, J. F., & Hogg, D. W. 2016, ApJ, 833, 98 [Google Scholar]
  35. Calzetti, D., Kennicutt, R. C., Engelbracht, C. W., et al. 2007, ApJ, 666, 870 [Google Scholar]
  36. Cárdenas Vázquez, M. C., Dorner, B., Huber, A., et al. 2018, PASP, 130, 025003 [CrossRef] [Google Scholar]
  37. Carleton, T., Cooper, M. C., Bolatto, A. D., et al. 2017, MNRAS, 467, 4886 [Google Scholar]
  38. Carniani, S., Marconi, A., Maiolino, R., et al. 2017, A&A, 605, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Catalán-Torrecilla, C., Gil de Paz, A., Castillo-Morales, A., et al. 2017, ApJ, 848, 87 [Google Scholar]
  40. Cecil, G., Bland-Hawthorn, J., Veilleux, S., & Filippenko, A. V. 2001, ApJ, 555, 338 [NASA ADS] [CrossRef] [Google Scholar]
  41. Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, ArXiv e-prints [arXiv:1612.05560] [Google Scholar]
  42. Chapin, E. L., Berry, D. S., Gibb, A. G., et al. 2013, MNRAS, 430, 2545 [Google Scholar]
  43. Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 [Google Scholar]
  44. Chauke, P., van der Wel, A., Pacifici, C., et al. 2019, ApJ, 877, 48 [NASA ADS] [CrossRef] [Google Scholar]
  45. Cicone, C., Maiolino, R., Sturm, E., et al. 2014, A&A, 562, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Circosta, C., Mainieri, V., Padovani, P., et al. 2018, A&A, 620, A82 [EDP Sciences] [Google Scholar]
  47. Circosta, C., Mainieri, V., Lamperti, I., et al. 2021, A&A, 646, A96 [EDP Sciences] [Google Scholar]
  48. Colombo, D., Sanchez, S. F., Bolatto, A. D., et al. 2020, A&A, 644, A97 [EDP Sciences] [Google Scholar]
  49. Condon, J. J. 1992, ARA&A, 30, 575 [Google Scholar]
  50. Cresci, G., Mainieri, V., Brusa, M., et al. 2015a, ApJ, 799, 82 [Google Scholar]
  51. Cresci, G., Marconi, A., Zibetti, S., et al. 2015b, A&A, 582, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Croton, D. J., Springel, V., White, S. D. M., et al. 2006, MNRAS, 365, 11 [Google Scholar]
  53. Cutri, R. M., Wright, E. L., Conrow, T., et al. 2011, Explanatory Supplement to the WISE Preliminary Data Release Products, Tech. Rep. [Google Scholar]
  54. Curti, M., Mannucci, F., Cresci, G., & Maiolino, R. 2020, MNRAS, 491, 944 [Google Scholar]
  55. da Cunha, E., Charlot, S., & Elbaz, D. 2008, MNRAS, 388, 1595 [Google Scholar]
  56. Dale, D. A., & Helou, G. 2002, ApJ, 576, 159 [Google Scholar]
  57. Davies, R. L., Rich, J. A., Kewley, L. J., & Dopita, M. A. 2014a, MNRAS, 439, 3835 [Google Scholar]
  58. Davies, R. L., Kewley, L. J., Ho, I. T., & Dopita, M. A. 2014b, MNRAS, 444, 3961 [NASA ADS] [CrossRef] [Google Scholar]
  59. Davies, L. J. M., Robotham, A. S. G., Driver, S. P., et al. 2015, MNRAS, 452, 616 [NASA ADS] [CrossRef] [Google Scholar]
  60. Davies, R. L., Groves, B., Kewley, L. J., et al. 2016, MNRAS, 462, 1616 [NASA ADS] [CrossRef] [Google Scholar]
  61. Davies, R. L., Groves, B., Kewley, L. J., et al. 2017, MNRAS, 470, 4974 [NASA ADS] [CrossRef] [Google Scholar]
  62. Dempsey, J. T., Friberg, P., Jenness, T., et al. 2013, MNRAS, 430, 2534 [Google Scholar]
  63. DePoy, D. L., Abbott, T., Annis, J., et al. 2008, SPIE Conf. Ser., 7014, 70140E [NASA ADS] [Google Scholar]
  64. Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604 [NASA ADS] [CrossRef] [Google Scholar]
  65. Di Matteo, P., Bournaud, F., Martig, M., et al. 2008, A&A, 492, 31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Dopita, M. A., Kewley, L. J., Sutherland, R. S., & Nicholls, D. C. 2016, Ap&SS, 361, 61 [NASA ADS] [CrossRef] [Google Scholar]
  67. Draper, N. R. 1998, Applied Regression Analysis (New York: John Wiley& Sons) [CrossRef] [Google Scholar]
  68. Elbaz, D., Daddi, E., Le Borgne, D., et al. 2007, A&A, 468, 33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Elbaz, D., Dickinson, M., Hwang, H. S., et al. 2011, A&A, 533, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Ellison, S. L., Patton, D. R., Simard, L., & McConnachie, A. W. 2008, AJ, 135, 1877 [NASA ADS] [CrossRef] [Google Scholar]
  71. Ellison, S. L., Thorp, M. D., Lin, L., et al. 2020, MNRAS, 493, L39 [CrossRef] [Google Scholar]
  72. Faber, S. M., Willmer, C. N. A., Wolf, C., et al. 2007, ApJ, 665, 265 [Google Scholar]
  73. Fabian, A. C. 2012, ARA&A, 50, 455 [Google Scholar]
  74. Fabian, A. C., Sanders, J. S., Taylor, G. B., et al. 2006, MNRAS, 366, 417 [NASA ADS] [CrossRef] [Google Scholar]
  75. Farrah, D., Urrutia, T., Lacy, M., et al. 2012, ApJ, 745, 178 [NASA ADS] [CrossRef] [Google Scholar]
  76. Fernández-Ontiveros, J. A., Tristram, K. R. W., Hönig, S., Gandhi, P., & Weigelt, G. 2018, A&A, 611, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Feruglio, C., Maiolino, R., Piconcelli, E., et al. 2010, A&A, 518, L155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Fitzpatrick, E. L. 1999, PASP, 111, 63 [Google Scholar]
  79. Flores Velázquez, J. A., Gurvich, A. B., Faucher-Giguère, C.-A., et al. 2021, MNRAS, 501, 4812 [CrossRef] [Google Scholar]
  80. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
  81. Fraser-McKelvie, A., Merrifield, M., Aragón-Salamanca, A., et al. 2020, MNRAS, 499, 1116 [CrossRef] [Google Scholar]
  82. Gallimore, J. F., Elitzur, M., Maiolino, R., et al. 2016, ApJ, 829, L7 [Google Scholar]
  83. García-Burillo, S., Combes, F., Usero, A., et al. 2014, A&A, 567, A125 [Google Scholar]
  84. Gaspari, M., Brighenti, F., & Temi, P. 2015, A&A, 579, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Gaspari, M., Eckert, D., Ettori, S., et al. 2019, ApJ, 884, 169 [NASA ADS] [CrossRef] [Google Scholar]
  86. Gaspari, M., Tombesi, F., & Cappi, M. 2020, Nat. Astron., 4, 10 [Google Scholar]
  87. Genzel, R., Tacconi, L. J., Combes, F., et al. 2012, ApJ, 746, 69 [Google Scholar]
  88. Giustini, M., Cappi, M., Chartas, G., et al. 2011, A&A, 536, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Groves, B. A., Heckman, T. M., & Kauffmann, G. 2006, MNRAS, 371, 1559 [NASA ADS] [CrossRef] [Google Scholar]
  90. Harper, D. A., Runyan, M. C., Dowell, C. D., et al. 2018, J. Astron. Instrum., 7, 1840008 [Google Scholar]
  91. Harrison, C. M., Alexander, D. M., Mullaney, J. R., et al. 2012, ApJ, 760, L15 [NASA ADS] [CrossRef] [Google Scholar]
  92. Hart, R. E., Bamford, S. P., Willett, K. W., et al. 2016, MNRAS, 461, 3663 [NASA ADS] [CrossRef] [Google Scholar]
  93. Hasinger, G. 2008, A&A, 490, 905 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  94. Hayward, C. C., Lanz, L., Ashby, M. L. N., et al. 2014, MNRAS, 445, 1598 [NASA ADS] [CrossRef] [Google Scholar]
  95. Henden, A. A., Welch, D. L., Terrell, D., & Levine, S. E. 2009, Am. Astron. Soc. Meet. Abstr., 214, 407.02 [NASA ADS] [Google Scholar]
  96. Hickox, R. C., Mullaney, J. R., Alexander, D. M., et al. 2014, ApJ, 782, 9 [NASA ADS] [CrossRef] [Google Scholar]
  97. Ho, L. C. 2005, ApJ, 629, 680 [CrossRef] [Google Scholar]
  98. Holland, W. S., Bintley, D., Chapin, E. L., et al. 2013, MNRAS, 430, 2513 [Google Scholar]
  99. Hopkins, P. F., & Elvis, M. 2010, MNRAS, 401, 7 [NASA ADS] [CrossRef] [Google Scholar]
  100. Hopkins, P. F., Hernquist, L., Cox, T. J., & Kereš, D. 2008, ApJS, 175, 356 [Google Scholar]
  101. Hopkins, P. F., Cox, T. J., Hernquist, L., et al. 2013, MNRAS, 430, 1901 [Google Scholar]
  102. Husemann, B., Jahnke, K., Sánchez, S. F., et al. 2014, MNRAS, 443, 755 [NASA ADS] [CrossRef] [Google Scholar]
  103. Husemann, B., Scharwächter, J., Davis, T. A., et al. 2019, A&A, 627, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  104. Husemann, B., Singha, M., Scharwächter, J., et al. 2022, A&A, 659, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Jarrett, T. H., Chester, T., Cutri, R., et al. 2000, AJ, 119, 2498 [Google Scholar]
  106. Kakkad, D., Mainieri, V., Brusa, M., et al. 2017, MNRAS, 468, 4205 [Google Scholar]
  107. Karim, A., Schinnerer, E., Martínez-Sansigre, A., et al. 2011, ApJ, 730, 61 [Google Scholar]
  108. Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346, 1055 [Google Scholar]
  109. Kennicutt, R. C., Jr. 1998a, ARA&A, 36, 189 [Google Scholar]
  110. Kennicutt, R. C., Jr. 1998b, ApJ, 498, 541 [Google Scholar]
  111. Kennicutt, R. C., Jr., & Kent, S. M. 1983, AJ, 88, 1094 [NASA ADS] [CrossRef] [Google Scholar]
  112. Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183 [Google Scholar]
  113. Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 121 [Google Scholar]
  114. Kewley, L. J., Geller, M. J., & Jansen, R. A. 2004, AJ, 127, 2002 [Google Scholar]
  115. Kewley, L. J., Groves, B., Kauffmann, G., & Heckman, T. 2006, MNRAS, 372, 961 [Google Scholar]
  116. Kewley, L. J., Rupke, D., Jabran Zahid, H., Geller, M. J., & Barton, E. J. 2010, ApJ, 721, L48 [NASA ADS] [CrossRef] [Google Scholar]
  117. Khoperskov, S., Haywood, M., Di Matteo, P., Lehnert, M. D., & Combes, F. 2018, A&A, 609, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  118. Kim, M., Ho, L. C., & Im, M. 2006, ApJ, 642, 702 [NASA ADS] [CrossRef] [Google Scholar]
  119. Koss, M. J., Strittmatter, B., Lamperti, I., et al. 2021, ApJS, 252, 29 [CrossRef] [Google Scholar]
  120. Lacerda, E. A. D., Sánchez, S. F., Cid Fernandes, R., et al. 2020, MNRAS, 492, 3073 [Google Scholar]
  121. Le Févre, O., Saisse, M., Mancini, D., et al. 2003, Proc. SPIE, 4841, 1671 [Google Scholar]
  122. Lequeux, J., Peimbert, M., Rayo, J. F., Serrano, A., & Torres-Peimbert, S. 1979, A&A, 500, 145 [NASA ADS] [Google Scholar]
  123. Leung, G. C. K., Coil, A. L., Azadi, M., et al. 2017, ApJ, 849, 48 [NASA ADS] [CrossRef] [Google Scholar]
  124. Lin, L., Pan, H.-A., Ellison, S. L., et al. 2019, ApJ, 884, L33 [NASA ADS] [CrossRef] [Google Scholar]
  125. Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 [Google Scholar]
  126. Mahony, E. K., Oonk, J. B. R., Morganti, R., et al. 2016, MNRAS, 455, 2453 [NASA ADS] [CrossRef] [Google Scholar]
  127. Maiolino, R., Gallerani, S., Neri, R., et al. 2012, MNRAS, 425, L66 [Google Scholar]
  128. Maiolino, R., Russell, H. R., Fabian, A. C., et al. 2017, Nature, 544, 202 [Google Scholar]
  129. Mancini, C., Daddi, E., Juneau, S., et al. 2019, MNRAS, 489, 1265 [Google Scholar]
  130. Mannucci, F., Cresci, G., Maiolino, R., Marconi, A., & Gnerucci, A. 2010, MNRAS, 408, 2115 [NASA ADS] [CrossRef] [Google Scholar]
  131. Mannucci, F., Belfiore, F., Curti, M., et al. 2021, MNRAS, 508, 1582 [NASA ADS] [CrossRef] [Google Scholar]
  132. Marasco, A., Cresci, G., Nardini, E., et al. 2020, A&A, 644, A15 [EDP Sciences] [Google Scholar]
  133. Marino, R. A., Rosales-Ortega, F. F., Sánchez, S. F., et al. 2013, A&A, 559, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  134. Martin, P., & Roy, J.-R. 1994, ApJ, 424, 599 [NASA ADS] [CrossRef] [Google Scholar]
  135. McDonald, M., McNamara, B. R., Calzadilla, M. S., et al. 2021, ApJ, 908, 85 [NASA ADS] [CrossRef] [Google Scholar]
  136. McNamara, B. R., & Nulsen, P. E. J. 2007, ARA&A, 45, 117 [NASA ADS] [CrossRef] [Google Scholar]
  137. McNamara, B. R., Nulsen, P. E. J., Wise, M. W., et al. 2005, Nature, 433, 45 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  138. Mignoli, M., Vignali, C., Gilli, R., et al. 2013, A&A, 556, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  139. Mihos, J. C., Richstone, D. O., & Bothun, G. D. 1992, ApJ, 400, 153 [NASA ADS] [CrossRef] [Google Scholar]
  140. Mommert, M. 2017, Astron. Comput., 18, 47 [NASA ADS] [CrossRef] [Google Scholar]
  141. Morganti, R., Oosterloo, T., Oonk, J. B. R., Frieswijk, W., & Tadhunter, C. 2015, A&A, 580, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  142. Mukherjee, D., Wagner, A. Y., Bicknell, G. V., et al. 2018, MNRAS, 476, 80 [NASA ADS] [CrossRef] [Google Scholar]
  143. Mullaney, J. R., Alexander, D. M., Aird, J., et al. 2015, MNRAS, 453, L83 [Google Scholar]
  144. Murphy, E. J., Condon, J. J., Schinnerer, E., et al. 2011, ApJ, 737, 67 [Google Scholar]
  145. Nandra, K., Georgakakis, A., Willmer, C. N. A., et al. 2007, ApJ, 660, L11 [NASA ADS] [CrossRef] [Google Scholar]
  146. Nascimento, J. C. D., Storchi-Bergmann, T., Mallmann, N. D., et al. 2019, MNRAS, 486, 5075 [NASA ADS] [CrossRef] [Google Scholar]
  147. Nelson, D., Pillepich, A., Springel, V., et al. 2018, MNRAS, 475, 624 [Google Scholar]
  148. Nesvadba, N. P. H., Lehnert, M. D., De Breuck, C., Gilbert, A. M., & van Breugel, W. 2008, A&A, 491, 407 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  149. Neumann, J., Gadotti, D. A., Wisotzki, L., et al. 2019, A&A, 627, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  150. Nielsen, N. M., Churchill, C. W., Kacprzak, G. G., Murphy, M. T., & Evans, J. L. 2015, ApJ, 812, 83 [NASA ADS] [CrossRef] [Google Scholar]
  151. Noeske, K. G., Weiner, B. J., Faber, S. M., et al. 2007, ApJ, 660, L43 [CrossRef] [Google Scholar]
  152. Noll, S., Burgarella, D., Giovannoli, E., et al. 2009, A&A, 507, 1793 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  153. Page, M. J., Symeonidis, M., Vieira, J. D., et al. 2012, Nature, 485, 213 [NASA ADS] [CrossRef] [Google Scholar]
  154. Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, AJ, 124, 266 [Google Scholar]
  155. Peng, Y.-J., Lilly, S. J., Kovač, K., et al. 2010, ApJ, 721, 193 [Google Scholar]
  156. Perna, M., Sargent, M. T., Brusa, M., et al. 2018, A&A, 619, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  157. Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 [NASA ADS] [CrossRef] [Google Scholar]
  158. Pillepich, A., Springel, V., Nelson, D., et al. 2018, MNRAS, 473, 4077 [Google Scholar]
  159. Piotrowska, J. M., Bluck, A. F. L., Maiolino, R., Concas, A., & Peng, Y. 2020, MNRAS, 492, L6 [Google Scholar]
  160. Popesso, P., Concas, A., Morselli, L., et al. 2020, MNRAS, 496, 2531 [NASA ADS] [CrossRef] [Google Scholar]
  161. Prevot, M. L., Lequeux, J., Maurice, E., Prevot, L., & Rocca-Volmerange, B. 1984, A&A, 132, 389 [Google Scholar]
  162. Ramos Padilla, A. F., Ashby, M. L. N., Smith, H. A., et al. 2020, MNRAS, 499, 4325 [NASA ADS] [CrossRef] [Google Scholar]
  163. Renzini, A., & Peng, Y.-J. 2015, ApJ, 801, L29 [Google Scholar]
  164. Rosario, D. J. 2019, Astrophysics Source Code Library [record ascl:1904.011] [Google Scholar]
  165. Rosario, D. J., Alexander, D. M., Moldon, J., et al. 2021, MNRAS, 505, 5283 [NASA ADS] [CrossRef] [Google Scholar]
  166. Rose, T., Edge, A. C., Combes, F., et al. 2019, MNRAS, 489, 349 [NASA ADS] [CrossRef] [Google Scholar]
  167. Roth, M. M., Kelz, A., Fechner, T., et al. 2005, PASP, 117, 620 [Google Scholar]
  168. Saintonge, A., Catinella, B., Tacconi, L. J., et al. 2017, ApJS, 233, 22 [Google Scholar]
  169. Salim, S., Rich, R. M., Charlot, S., et al. 2007, ApJS, 173, 267 [NASA ADS] [CrossRef] [Google Scholar]
  170. Sánchez, S. F., Rosales-Ortega, F. F., Iglesias-Páramo, J., et al. 2014, A&A, 563, A49 [CrossRef] [EDP Sciences] [Google Scholar]
  171. Sánchez, S. F., Barrera-Ballesteros, J. K., Sánchez-Menguiano, L., et al. 2017, MNRAS, 469, 2121 [CrossRef] [Google Scholar]
  172. Sánchez, S. F., Barrera-Ballesteros, J. K., López-Cobá, C., et al. 2019, MNRAS, 484, 3042 [CrossRef] [Google Scholar]
  173. Sanders, R. L., Jones, T., Shapley, A. E., et al. 2020, ApJ, 888, L11 [Google Scholar]
  174. Sani, E., Lutz, D., Risaliti, G., et al. 2010, MNRAS, 403, 1246 [Google Scholar]
  175. Santoro, F., Oonk, J. B. R., Morganti, R., Oosterloo, T. A., & Tadhunter, C. 2016, A&A, 590, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  176. Schaefer, A. L., Croom, S. M., Scott, N., et al. 2019, MNRAS, 483, 2851 [NASA ADS] [CrossRef] [Google Scholar]
  177. Schawinski, K., Virani, S., Simmons, B., et al. 2009, ApJ, 692, L19 [NASA ADS] [CrossRef] [Google Scholar]
  178. Schawinski, K., Urry, C. M., Simmons, B. D., et al. 2014, MNRAS, 440, 889 [Google Scholar]
  179. Schaye, J., Crain, R. A., Bower, R. G., et al. 2015, MNRAS, 446, 521 [Google Scholar]
  180. Schirmer, M. 2013, ApJS, 209, 21 [NASA ADS] [CrossRef] [Google Scholar]
  181. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [Google Scholar]
  182. Scholtz, J., Harrison, C. M., Rosario, D. J., et al. 2020, MNRAS, 492, 3194 [NASA ADS] [CrossRef] [Google Scholar]
  183. Schroetter, I., Bouché, N. F., Zabl, J., et al. 2019, MNRAS, 490, 4368 [CrossRef] [Google Scholar]
  184. Schulz, B., Marton, G., Valtchanov, I., et al. 2017, ArXiv e-prints [arXiv:1706.00448] [Google Scholar]
  185. Shangguan, J., Ho, L. C., & Xie, Y. 2018, ApJ, 854, 158 [NASA ADS] [CrossRef] [Google Scholar]
  186. Shimizu, T. T., Mushotzky, R. F., Meléndez, M., Koss, M., & Rosario, D. J. 2015, MNRAS, 452, 1841 [NASA ADS] [CrossRef] [Google Scholar]
  187. Silva, L., Granato, G. L., Bressan, A., & Danese, L. 1998, ApJ, 509, 103 [Google Scholar]
  188. Silva, L., Maiolino, R., & Granato, G. L. 2004, MNRAS, 355, 973 [Google Scholar]
  189. Silverman, J. D., Lamareille, F., Maier, C., et al. 2009, ApJ, 696, 396 [NASA ADS] [CrossRef] [Google Scholar]
  190. Singha, M., Husemann, B., Urrutia, T., et al. 2022, A&A, 659, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  191. Skibba, R. A., Engelbracht, C. W., Dale, D., et al. 2011, ApJ, 738, 89 [NASA ADS] [CrossRef] [Google Scholar]
  192. Smirnova-Pinchukova, I., Husemann, B., Busch, G., et al. 2019, A&A, 626, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  193. Smith, M. W. L., Gomez, H. L., Eales, S. A., et al. 2012, ApJ, 748, 123 [NASA ADS] [CrossRef] [Google Scholar]
  194. Smith, K. L., Koss, M., Mushotzky, R., et al. 2020, ApJ, 904, 83 [NASA ADS] [CrossRef] [Google Scholar]
  195. Speagle, J. S., Steinhardt, C. L., Capak, P. L., & Silverman, J. D. 2014, ApJS, 214, 15 [Google Scholar]
  196. Springel, V., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 776 [Google Scholar]
  197. Stern, J., & Laor, A. 2013, MNRAS, 431, 836 [NASA ADS] [CrossRef] [Google Scholar]
  198. Storchi-Bergmann, T., Schmitt, H. R., Calzetti, D., & Kinney, A. L. 1998, AJ, 115, 909 [NASA ADS] [CrossRef] [Google Scholar]
  199. Strateva, I., Ivezić, Ž., Knapp, G. R., et al. 2001, AJ, 122, 1861 [CrossRef] [Google Scholar]
  200. Tacconi, L. J., Genzel, R., Saintonge, A., et al. 2018, ApJ, 853, 179 [Google Scholar]
  201. Thorp, M. D., Ellison, S. L., Simard, L., Sánchez, S. F., & Antonio, B. 2019, MNRAS, 482, L55 [Google Scholar]
  202. Tremblay, G. R., Combes, F., Oonk, J. B. R., et al. 2018, ApJ, 865, 13 [NASA ADS] [CrossRef] [Google Scholar]
  203. Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898 [Google Scholar]
  204. Tumlinson, J., Thom, C., Werk, J. K., et al. 2011, Science, 334, 948 [CrossRef] [Google Scholar]
  205. Utomo, D., Bolatto, A. D., Wong, T., et al. 2017, ApJ, 849, 26 [NASA ADS] [CrossRef] [Google Scholar]
  206. Valdes, F., Gruendl, R., & DES Project 2014, ASP Conf. Ser., 485, 379 [NASA ADS] [Google Scholar]
  207. Vogelsberger, M., Genel, S., Springel, V., et al. 2014a, MNRAS, 444, 1518 [Google Scholar]
  208. Vogelsberger, M., Genel, S., Springel, V., et al. 2014b, Nature, 509, 177 [Google Scholar]
  209. Walcher, C. J., Coelho, P. R. T., Gallazzi, A., et al. 2015, A&A, 582, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  210. Walsh, J. L., Barth, A. J., Ho, L. C., & Sarzi, M. 2013, ApJ, 770, 86 [NASA ADS] [CrossRef] [Google Scholar]
  211. Wang, T. G., Brinkmann, W., Wamsteker, W., Yuan, W., & Wang, J. X. 1999, MNRAS, 307, 821 [NASA ADS] [CrossRef] [Google Scholar]
  212. Weaver, J., Husemann, B., Kuntschner, H., et al. 2018, A&A, 614, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  213. Wild, V., Heckman, T., & Charlot, S. 2010, MNRAS, 405, 933 [NASA ADS] [Google Scholar]
  214. Willett, K. W., Lintott, C. J., Bamford, S. P., et al. 2013, MNRAS, 435, 2835 [Google Scholar]
  215. Wisotzki, L., Christlieb, N., Bade, N., et al. 2000, A&A, 358, 77 [NASA ADS] [Google Scholar]
  216. Wittor, D., & Gaspari, M. 2020, MNRAS, 498, 4983 [NASA ADS] [CrossRef] [Google Scholar]
  217. Wolf, C., Onken, C. A., Luvaul, L. C., et al. 2018, PASA, 35, e010 [Google Scholar]
  218. Woo, J.-H., Son, D., & Bae, H.-J. 2017, ApJ, 839, 120 [Google Scholar]
  219. Wyder, T. K., Martin, D. C., Schiminovich, D., et al. 2007, ApJS, 173, 293 [Google Scholar]
  220. Wylezalek, D., & Zakamska, N. L. 2016, MNRAS, 461, 3724 [NASA ADS] [CrossRef] [Google Scholar]
  221. Wylezalek, D., Zakamska, N. L., Greene, J. E., et al. 2018, MNRAS, 474, 1499 [NASA ADS] [CrossRef] [Google Scholar]
  222. Yamamura, I., Makiuti, S., Ikeda, N., et al. 2009, ASP Conf. Ser., 418, 3 [NASA ADS] [Google Scholar]
  223. Yang, G., Boquien, M., Buat, V., et al. 2020, MNRAS, 491, 740 [Google Scholar]
  224. Yates, R. M., Kauffmann, G., & Guo, Q. 2012, MNRAS, 422, 215 [NASA ADS] [CrossRef] [Google Scholar]
  225. York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 120, 1579 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: SED models

thumbnail Fig. A.1.

Overview of the AGNfitter SED modeling results for the entire CARS sample. Black data points represent the measured panchromatic photometry collected for each object. Upper limits are highlighted as black arrows. The red dots are the predicted photometry for XX MCMC relations from the superposition of individual SED components. The red, yellow, blue, purple and green lines correspond to the total SED, the stellar component, the AGN component, the hot dust component and the cold dust component, respectively.

thumbnail Fig. A.1.

continued.

thumbnail Fig. A.1.

continued.

Appendix B: CARS sample BPT morphology

thumbnail Fig. B.1.

BPT morphology and RAINBOW results for the sample. The BPT and a spatial map for each galaxy show SF fraction. The red dotted line shows the threshold for the AGN basis spaxels. The contours on the maps represent MUSE whitelight contours. The gray pixels on the maps represent those pixels, where Hα S/N is more than 3σ, but the other lines are weaker so the datapoint cannot be shown on a BPT.

thumbnail Fig. B.1.

continued.

thumbnail Fig. B.1.

continued.

thumbnail Fig. B.1.

continued.

thumbnail Fig. B.1.

continued.

thumbnail Fig. B.1.

continued.

thumbnail Fig. B.1.

continued.

thumbnail Fig. B.1.

continued.

thumbnail Fig. B.1.

continued.

Table B.1.

Separation radii for galaxies with star-forming and/or AGN-ionised cloud.

Appendix C: Metallicity maps

thumbnail Fig. C.1.

Same as Fig. 5 for the entire sample.

thumbnail Fig. C.1.

continued.

thumbnail Fig. C.1.

continued.

thumbnail Fig. C.1.

continued.

All Tables

Table 1.

Multiband aperture photometry measurements.

Table 2.

AGN host galaxy properties inferred from SED modeling.

Table 3.

Hα-based SFRs and metallicity estimated from the IFU emission lines.

Table B.1.

Separation radii for galaxies with star-forming and/or AGN-ionised cloud.

All Figures

thumbnail Fig. 1.

Photometric broad bands from a number of surveys and instruments used for the SED modeling. The y-axis is normalized roughly by the year of observational campaign to visually compare bands with similar coverage. When more than one option was available for certain bands, the preference was given to the band sets, which are temporally closer to each other (especially in near infrared and optical ranges) to minimize the impact of AGN variability.

In the text
thumbnail Fig. 2.

SED for HE 0114−0015 and the best-fit model determined by AGNfitter (Calistro Rivera et al. 2016). 20 MCMC realizations of the SED model are shown for the broadband photometric data (black points with error bars) where the red points are the corresponding predictions from the model. The red lines represent the total model, which consists of four components: the cold and warm dust in star-forming regions (green lines), the torus of AGN-heated dust (purple lines), the stellar continuum (yellow lines), and the AGN accretion disk (blue lines). On the right panel is an additional constraint to the accretion disk model (e.g. g and i AGN photometry is included).

In the text
thumbnail Fig. 3.

BPT diagram of HE 0853+0102 as an example for a mixing sequence between H II regions and AGN photoionization across the host galaxy. We highlight three ways of determining the SF fraction from the BPT diagram. Left panel: a simple demarcation line approach where everything above the Kewley et al. (2001) line is classified as pure AGN and below as pure star formation. Middle panel: simple linear combination with a χ2 fit using a fixed line vector as a reference point for the characteristic AGN and H II line ratios. The colored halos around the points represent the difference between the calculated SF fraction adopting two different different reference pairs (triangle and square symbols) to highlight the systematic uncertainties introduced by a certain choice. Right panel: MCMC fitting results of the star formation fraction adopting a large cloud of potential reference points for the AGN (red points) and H II region (magenta points) which are explored during the MCMC sampling.

In the text
thumbnail Fig. 4.

Schematic sketch of the BPT morphology populations. Upper row: on the right there is a representation of a galaxy with a central region in green, the rest of the galactic body in blue, and an AGN ionization cone in yellow. These regions populate specific areas on the BPT diagram: left diagram: star-forming cloud with two different sizes; middle diagram: AGN-ionized cloud with two sizes; right diagram: mixing sequence elongated toward the AGN area and a mixing sequence elongated toward the LINER area. Middle row: BPT diagram schematic examples for the SF cloud, AGN cloud, mixing sequence and LINER mixing sequence dominated galaxies. Lower row: BPT diagram and spatial map examples for SF cloud and AGN cloud combined with a mixing sequence in one galaxy.

In the text
thumbnail Fig. 5.

Example of the 2D and radial metallicity distribution for HE 2211−3903. The oxygen abundance was measured using the N2S2 index which is almost insensitive to the excitation. On the left panel the full 2D distribution is shown and on the right panel a radial projection in units of the effective radius taken from the single Sersic model. A linear fit is performed and shown as the red line.

In the text
thumbnail Fig. 6.

Comparison of the sSFR against stellar mass for Hα- and IR-based SFR tracers. The CARS objects are shown as red points with error bars, while xCOLD GASS (xCG) are in small gray dots, KINGFISH are in blue and Herschel Reference Survey (HRS) in green. Specific SFR is calculated from Hα luminosity (left panel) derived after BPT morphology analysis and RAINBOW fitting using the SFR formula SFRHα = 5.3 × 10−42 × LHα (Calzetti et al. 2007). IR specific SFR (right panel) derived from the AGNFITTER modeling using the Murphy et al. (2011) SFR formula. The stellar mass is taken from the AGNFITTER modeling parameters. The background on the left plot is the density map of SDSS galaxies (Brinchmann et al. 2004). Linear models to the star-forming main sequence are shown as black solid lines from Renzini & Peng (2015) and Shimizu et al. (2015), respectively. The best-fit linear relation from the xCG sample is shown as the blue dashed line in both cases. Lower panel: residuals around the adopted star-forming main sequence where the colored bands (red, blue, and purple) emphasize three bins used further in Fig. 8.

In the text
thumbnail Fig. 7.

Difference between measured and expected Hα (upper panels) and FIR luminosity (lower panels) as a function of stellar mass, CO(1−0) luminosity and 12 + log(O/H) ionized gas metallicity. Here we used our more complex model (Eqs. (1) and (2)) to predict the SFR-sensitive luminosities. The black dots are the data from our xCG training sample and the red symbols are the CARS data with triangles indicating upper limits. The colored bands (red, blue and purple) define three bins in residual luminosity as discussed in the text and used for Fig. 8.

In the text
thumbnail Fig. 8.

Comparison of SFR residuals from the SFMS and SFMS+gas model against various AGN parameters. Upper left panel: histograms of the residuals for SFMS and SFMS+gas for both Hα and IR models, with the red, blue, and purple colors corresponding to the previous plots (Figs. 6 and 7, background colors). The other panels show the comparison of the AGN parameters taken from Husemann et al. (2022) and Singha et al. (2022): logarithm of BH mass (upper center), logarithm of AGN bolometric luminosity (upper right), number of the objects, containing extended outflows (lower left), logarithm of Eddington ratio (lower center), and b/a morphological parameter (lower right). Each of the four models is represented with three bins (from left to right: negative red, central blue, positive purple) and the forth bin with the upper limits (gray). The error bars of the bins also represent the distribution of the parameters within the bins. The black dots and lines highlight the trend of the mean values that are written over the bins.

In the text
thumbnail Fig. 9.

8−1000 μm IR luminosity LIR (left panel) and the CO (1−0) luminosity LCO (right panel) plotted against Hα luminosity. The training sample from the xCG (gray dots) is compared to the CARS objects (red circles and triangles for the 5σ upper limits). Galaxy mergers in the xCG sample as classified by GalaxyZoo2 are additionally marked (black crosses) on the left panel. The blue dashed line corresponds to Eqs. (12) and (13) are derived by the combination of literature calibrations as described in the main text. The black lines represent a linear relation fitted to the XCOLD GASS reference sample, see Eqs. (14) and (15). The red dotted line is fitted to the CARS objects with the same slope as the black line, therefore highlight the shift in the normalization.

In the text
thumbnail Fig. 10.

Comparison of different SFR measurements for individual objects. Left column: Hα-, IR-, and CO(1−0)-based SFRs for the CARS objects sorted by the SFRIR − SFRHα difference. Central column: near infrared variability; the objects with no AGN or known changing look AGN are also marked here. Right column: morphology of the host galaxies. Lower panel: contains the objects with either Hα and/or IR SFR uncertainties (upper limits or missing data).

In the text
thumbnail Fig. A.1.

Overview of the AGNfitter SED modeling results for the entire CARS sample. Black data points represent the measured panchromatic photometry collected for each object. Upper limits are highlighted as black arrows. The red dots are the predicted photometry for XX MCMC relations from the superposition of individual SED components. The red, yellow, blue, purple and green lines correspond to the total SED, the stellar component, the AGN component, the hot dust component and the cold dust component, respectively.

In the text
thumbnail Fig. B.1.

BPT morphology and RAINBOW results for the sample. The BPT and a spatial map for each galaxy show SF fraction. The red dotted line shows the threshold for the AGN basis spaxels. The contours on the maps represent MUSE whitelight contours. The gray pixels on the maps represent those pixels, where Hα S/N is more than 3σ, but the other lines are weaker so the datapoint cannot be shown on a BPT.

In the text
thumbnail Fig. C.1.

Same as Fig. 5 for the entire sample.

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.