arXiv:2107.04213v1 [astro-ph.GA] 9 Jul 2021 in hc,i un a rdc h togcorrelations forma- strong star the produce and can growth turn, in BH which, regulating nuclei tion, im- in an galactic play role SMBH active portant accreting from studies, originate that (AGN) theoretical to cording errs ert 2000 Merritt & Ferrarese n t otglx ntrso hi omto and formation their of SMBH terms host a in the between galaxy ( of evolution connection host component its causal spheroidal and a stel- the imply and of galaxy dispersion, mass) prop- velocity physical lar luminosity, the cor- (e.g., and strong mass erties The SMBH ubiquitous. between are relations (SMBHs) holes black 1 eateto srnm n topei cecs Kyungpo Sciences, Atmospheric and Astronomy of Department ∗ yee sn L using 2021 Typeset 12, July version Draft nmsieglxe htehbtble,supermassive bulges, exhibit that galaxies massive In nAtooy n. ne AAcnrc A5255 These #15444. program NAS5-26555. with contract associated NASA Res are Institu for under observations Universities Science Inc., of Astronomy, Telescope Association the in Space by operated the is Space which at Hubble NASA/ESA obtained the with Telescope, made observations on Based 2 eateto hsc n srnm,42 rdrc Reines Frederick 4129 Astronomy, and Physics of Department ubeSaeTlsoeIaigSre fLwRdhf Swi Low-Redshift of Survey Imaging Telescope Space Hubble A Keywords: t of stages those evolutionary with Howe and compared mechanism luminosity triggering fraction. identical. bolometric the merging lower that and a imply at distribution may A peaks type 2 AGN type Hubble 2 and their type 1 Type s of ratio, ratio. terms Eddington Eddington in enhanced the AG an with to luminous correlated lead triggering weakly necessarily funct in is role a fraction important merging as an the strongly plays increases merging that fraction revealing merging The t (18% respe because features galaxies galaxies galaxies, disturbance). host merging peculiar of late-type as fraction and and classified minor are type a early-type sample by the Hubble of hosted the 8% is remaining as th sample such investigate the to galaxies, of us host enable optic the AGN to target of compared 154 obscuration of of images terms high-resolution in at unbiased (AGN) relatively nuclei sample galactic active Swift-BAT of epeetiiilrslsfo ubeSaeTlsoesaso im snapshot Telescope Space Hubble a from results initial present We 3 aora ta.1998 al. et Magorrian al nttt o srnm n srpyis eigUni Peking Astrophysics, and Astronomy for Institute Kavli A T 1. E X aais cie—glxe:ble aais udmna paramet fundamental general galaxies: quasars: — — bulges tometry galaxies: — active galaxies: INTRODUCTION 4 twocolumn eateto srnm,Sho fPyis eigUnivers Peking Physics, of School Astronomy, of Department ; ijnKim, Minjin omny&H 2013 Ho & Kormendy tl nAASTeX63 in style ; ehrte l 2000 al. et Gebhardt 1 ao .Barth, J. Aaron − 5)ehbtmrigfaue eg,tdltis hls rmajor or shells, tails, tidal (e.g., features merging exhibit 25%) .Ac- ). ABSTRACT earch ; < z kNtoa nvriy ag 0-0,Kra mkim.astro@ Korea; 702-701, Daegu University, National ok te, 2 usC Ho, C. Luis 0 al nvriyo aiona rie A96747,USA 92697-4575, CA Irvine, California, of University Hall, . .Tehr -a eeto ae hssample this makes selection X-ray hard The 1. lyosue ytesbpre ut ou ntp 2 type in visu- torus are dusty sub-parsec structures ( the AGN central by those obscured in while ally observed AGN, directly are 1 region type broad-line and disk tion galaxy the understand to nection. vital is hosts AGN with interactions or mergers ( can galaxies gas-rich which companion by phase, AGN triggered the be during primarily occur tion ewe MH n h rpriso hi otgalax- host their of ( properties ies the and SMBHs between e h urudn a n uts htteANis AGN the that so dust ex- and can gas outflow surrounding AGN-driven early the the the pel which in during occurs, enhanced ( central then merger be the the can gas-rich in of host by content stages the dust triggered of and mainly gas regions is the then AGN merging, an if contrary, 2008 2005 est,Biig107,Cia [email protected] China; 100871, Beijing versity, codn oteANuicto oe,teaccre- the model, unification AGN the to According ain1999 Fabian .Teeoe dniyn h egn etrsof features merging the identifying Therefore, ). forma- star and growth SMBH model, this In ). nouc 1993 Antonucci ,4 3, n uenSon Suyeon and t,Biig107,China 100871, Beijing ity, o fblmti G luminosity, AGN bolometric of ion tvl.Tehs aaiso the of galaxies host The ctively. ealdpooercstructure photometric detailed e getn htmrigde not does merging that uggesting lANslcinmtos The methods. selection AGN al Naeams indistinguishable almost are GN gn uvyo h otgalaxies host the of survey aging e r evl itre.Only disturbed. heavily are hey efidta 8 n 44% and 48% that find We . p n ye2ANaenot are AGN 2 type and 1 ype ; iMte ta.2005 al. et Matteo Di e,temrigfato of fraction merging the ver, tBTAtv Galaxies Active ft-BAT nti ape However, sample. this in N ftp G.Ti result This AGN. 1 type of adr ta.1988 al. et Sanders ; lcae l 2018 al. et Blecha ry&Pdvn 1995 Padovani & Urry r aais pho- galaxies: — ers 1 ; ; .Alt stage late A ). okn tal. et Hopkins al. et Hopkins − gmail.com .O the On ). otcon- host ∗ 2 KIM et al. eventually observed to be type 1 (Sanders et al. 1988; Ishibashi & Fabian 2016). Therefore, the nucleus can be 30 naturally obscured due to the dense material in the host galaxy on spatial scales larger than the dusty torus. Ac- N 15 cording to this evolutionary sequence, the host galaxies of type 2 AGN are more likely to be disturbed than those 46 of type 1 AGN, or might show an association with earlier merger stages. On the other hand, Elitzur et al. (2014) ) argued that type 1 AGN can evolve to intermediate- 1 45 − type (1.8-1.9), and finally to type 2 AGN as luminos- ity decreases. In this light, investigating the physical / erg s erg / properties of host galaxies is crucial not only to deter- bol 44 mine the AGN triggering mechanism but also to test L
AGN unification models. Host galaxy properties have ( log Type 1 Type 2 been used to test unified models for many AGN samples 43 Swift-BAT sample (e.g., Simcoe et al. 1997). Several studies have found that the merging fraction for obscured AGN is higher 0.00 0.05 0.10 15 30 than that for unobscured AGN (e.g., Urrutia et al. 2008; redshift N Glikman et al. 2015; Donley et al. 2018). Additionally, Compton-thick AGN appear to be preferentially hosted Figure 1. Redshift and bolometric luminosity distribu- by merging galaxies (Kocevski et al. 2015; Ricci et al. tions of the parent sample of Swift-BAT AGN (small gray dots). The Type 1 AGN observed for this HST program are 2017b; Koss et al. 2018). Koss et al. (2011), however, denoted by blue stars and blue histograms, while Type 2 ob- argued that the properties of the host galaxies of type 1 jects are given by red circles and red histograms. Bolometric and type 2 AGN are indistinguishable. luminosity is estimated from the intrinsic X-ray luminosity Numerous studies have investigated the properties of (Ricci et al. 2017) assuming a bolometric correction of 8. host galaxies of nearby and distant AGN. One of the primary goals of those studies is to identify the merg- tween SMBHs and host galaxies. The stellar populations ing features in the AGN hosts as a possible origin of of AGN host galaxies can provide a useful constraint AGN triggering. Several studies have claimed that on the SMBH-galaxy connection (e.g., S´anchez et al. nearby AGN are preferentially hosted by galaxies that 2004; Canalizo & Stockton 2013; Matsuoka et al. 2014; exhibit merging features compared with normal galax- Kim & Ho 2019; Zhao et al. 2019). Some studies ies (Koss et al. 2011; Hong et al. 2015; Goulding et al. have shown that host galaxies of type 1 AGN tend 2018; Ellison et al. 2019; Marian et al. 2020). In addi- to be bluer or overluminous compared with inactive tion, luminous AGN are more likely to be associated galaxies due to the former’s young stellar population with merging or interaction than less luminous AGN (e.g., S´anchez et al. 2004; Kim & Ho 2019). However, (e.g., Treister et al. 2012; Kim et al. 2017). However, other investigations have found no sign of young stel- some observational studies have found that there is lar populations in AGN hosts (e.g., Nolan et al. 2001; no evidence of an excess of merging fraction in vari- Bettoni et al. 2015). Additionally, the spheroid lumi- ous types of luminous AGN (e.g., Dunlop et al. 2003; nosity of type 1 AGN hosts has also been used to inves- Cisternas et al. 2011; B¨ohm et al. 2013; Villforth et al. tigate the relation between BH mass (MBH) and bulge 2014, 2017). For example, Zhao et al. (2019) argued luminosity (Lbul) (e.g., Kim et al. 2008; Greene et al. that only a minor fraction Type 2 quasar hosts (34%) 2010; Park et al. 2015; Kim & Ho 2019; Li et al. 2021). are mergers, suggesting that major merging plays only a The aforementioned studies have employed widely dif- limited role in triggering AGN. For narrow-line Seyfert ferent sample selection methods and datasets with di- 1 (NLS1) galaxies, the fraction of barred hosts is signif- verse properties in terms of wavelength or filter bands, icantly higher than the bar fraction of ordinary Type 1 spatial resolution, and depth, and these differences may AGN, implying that this subclass is predominantly trig- be largely responsible for the sometimes divergent con- gered by internal secular evolution rather than by merg- clusions reached in different investigations. Thus, de- ing (Crenshaw et al. 2003; Deo et al. 2006; Sani et al. spite the large observational effort invested into AGN 2010; Orban de Xivry et al. 2011; Kim et al. 2017). host galaxy studies over many years, further work based Finally, the photometric properties of host galaxies on large, uniformly selected samples with high-quality have been widely used to explore the coevolution be- imaging data can still provide new insights. Indeed, Host Galaxies of Swift-BAT AGN 3 there have been several attempts to conduct imaging surveys of nearby Seyferts using HST (e.g., Nelson et al. 1996; Malkan et al. 1998; Schade et al. 2000). Those surveys used early HST images with a relatively nar- N row field of view that did not always cover the full host galaxy (the Wide Field Planetary Camera 1 and 2), and early snapshot programs with these cameras obtained relatively shallow depth. Images of Type 1 AGN are of- ten saturated due to the bright nuclei, which makes it difficult to obtain a robust decomposition to determine Edd L bulge properties accurately. Furthermore, some of these bol
studies (e.g., Malkan et al. 1998) observed samples se- L lected from highly heterogeneous criteria based on AGN catalogs available at the time. To overcome these limitations, we have used HST to conduct an imaging survey of nearby AGN at z < 0.1 selected from the Swift-BAT AGN catalog. Swift-BAT 6 AGN are identified from hard X-ray data that, compared MBH/M ⊙ N with optical and infrared data, is less affected by extinc- tion and contamination from star formation. Therefore, Figure 2. The distributions of BH mass and Eddington ra- hard X-ray data provides a relatively uniform and homo- tio of the sample. MBH is measured using the virial method, the MBH − σ⋆ relation, or the MBH − M⋆ relation. geneous sample in terms of obscuration, with the excep- tion of very highly obscured (Compton-thick) AGN. The deep and coherent high-resolution imaging dataset from 2. SAMPLE AND DATA HST allows us to robustly investigate the structural 2.1. Sample Selection properties of the host galaxies for a large and uniformly In 2017, the Space Telescope Science Institute an- selected sample. Moreover, the Swift-BAT AGN sample nounced a special mid-cycle call proposals for “gap- has the benefit of extensive multiwavelength measure- filler” snapshot surveys to be conducted using the Ad- ments (e.g., BH mass, bolometric luminosity, Edding- vanced Camera for Surveys (ACS) Wide-Field Camera ton ratio, hydrogen column density, and optical classi- (WFC). The programmatic goal for this call was to pro- fication) that have been carried out in numerous stud- vide large, full-sky catalogs of snapshot targets that ies (Koss et al. 2017; Oh et al. 2017; Ricci et al. 2017a). would be used to fill HST scheduling gaps, at lower pri- Hence, our imaging survey is an excellent complement to ority than standard general observer and snapshot pro- the existing body of data for this benchmark AGN sam- grams. This survey (HST program 15444) was selected ple, providing an excellent dataset to study the detailed as one of three gap-filler programs in the 2017 call. The properties of nearby AGN hosts and to investigate the proposal call dictated that the total duration of snap- relationships between AGN and host properties while shot visits should be . 26 minutes, which is sufficient to minimizing sample selection biases. allow for two full-frame ACS/WFC exposures with a to- This paper presents a description of the morphological tal on-source exposure time of ∼ 700 s after acquisition properties of the host galaxies of these low-z Swift-BAT and other overheads. AGN using newly acquired HST I-band images for 154 Our sample is drawn from the 70-month Swift-BAT objects. Our primary goal is to explore the connection X-ray source catalog (Baumgartner et al. 2013). Since between the physical properties of AGN and those of our primary science goals involve morphological studies their hosts. Section 2 describes the sample selection, of AGN host galaxies, we focus on low-redshift objects, our imaging data, and the physical properties of the and we selected Swift-BAT AGN at z < 0.1. After cross- sample. Section 3 presents visual morphological classi- matching this catalog with the HST archive, we excluded fications. We discuss the correlation between the mor- objects already having images in an I-band equivalent phological structures of the hosts and the properties of filter (in most cases F814W) with the WF/PC2, ACS, the AGN in Section 4. A summary is given in Section 5. or WFC3 cameras, to avoid duplications. The remain- This work adopts the following cosmological parame- ing set of 543 objects provides an ideal gap-filler sample −1 −1 ters: H0 = 100h = 67.8 km s Mpc , Ωm = 0.308, with a large number of targets widely distributed over and ΩΛ =0.692 (Planck Collaboration et al. 2016). the full sky. No cuts to the sample were made based on 4 KIM et al.
Galactic latitude or foreground extinction, because HST of ∼0.′′05 and a field-of-view of ∼ 202′′×202′′. We mea- I-band imaging could in principle provide new morpho- sured the variance of the sky values after masking all logical classifications and host galaxy structural infor- objects to determine the depth of the images, which on mation even for targets in fields with moderate to high average we estimate to be ∼ 23.6 mag arcsec−1 (1 σ). extinction that otherwise might not be observed in other The short-exposure sub-array images for the type 1 surveys. Figure 1 shows the distribution of redshifts and objects exhibit horizontal bands introduced by bias bolometric luminosities for the parent Swift-BAT sam- striping noise (Grogin et al. 2010, 2011), which were re- ple and for the set of 543 targets making up the sample moved separately with the acs destripe plus task. for program 15444. The saturated pixels in the combined long-exposure im- age were replaced with the same pixels in the short, 2.2. Observations unsaturated image. The short and long exposures were Observations for Program 15444 began in 2018 March aligned using the position of the nucleus. The central and are ongoing (at a low rate) during the first half position of the nucleus for the long exposures was deter- of 2021. This paper includes data taken through 2021 mined after masking the saturated pixels. June. We used ACS/WFC with the F814W filter, which is approximately equivalent to the I band, to maximize 2.4. Physical Properties of AGN the brightness contrast between the host galaxy and the To investigate the connection between AGN and active nucleus. Suffering less dust attenuation than fil- their host galaxies, we assembled the optical spectral ters at shorter wavelengths, F814W also allows us to properties of the sample from the literature to esti- better probe the photometric substructures of the hosts. mate physical properties such as BH mass, accretion We obtained two dithered exposures per target with in- rate, and bolometric luminosity. For type 1 AGN, tegration times of 337 s, the minimum time required 2 we use the virial method (MBH ∝ Rv /G) to esti- to avoid buffer dump overheads. The AGN was cen- mate BH mass by combining the size (RBLR) and ve- tered on the WFC1 detector. For Type 1 (unobscured) locity dispersion (v) of the broad-line region (BLR). AGN, a nuclear point source is generally expected to The BLR size can be estimated through reverbera- be present at optical wavelengths (Nelson et al. 1996). tion mapping experiments (Blandford & McKee 1982) Accordingly, we took an additional 5 s exposure with and from single-epoch optical spectroscopy using the the WFC1-B512 subarray to obtain an unsaturated empirical relation between RBLR and AGN luminos- image of the nuclear point source and its immediate ity (Kaspi et al. 2000; Bentz et al. 2013). The intrin- surroundings for objects listed as Type 1 AGN by sic scatter of the BLR size-luminosity relation (Du et al. Baumgartner et al. (2013) or Koss et al. (2017). Of the 2018; Fonseca Alvarez et al. 2020) introduces significant 157 AGN targeted for observations, useful data were systematic uncertainty into the BH masses derived from obtained for 154 sources (Table 1), while three obser- single-epoch spectra. For single-epoch mass estimates, vations were lost due to guide star acquisition failures. the assumption of a single virial scaling factor con- For other targets, delayed guide star acquisition led to tributes additional error to the mass estimates, since slightly shorter exposure times than the planned 2×337 the true virial factor for a given AGN will depend on the s (SWIFT J0533.9−1318, SWIFT J0747.5+6057, and geometry and kinematics of the BLR. While the mean SWIFT J0753.1+4559). The set of 154 successfully ob- virial factor has been determined by assuming that AGN served targets includes 75 Type 1 and 79 Type 2 AGN. follow the same MBH −σ⋆ relation as the inactive galax- ies, its calibration appears to be sensitive morphologies 2.3. Data Reduction of galaxies and types of AGNs (e.g., Ho & Kim 2014). For the long exposures, we adope the standard re- Here, we recalculate the BH masses rather than using duction provided by the pipeline, including bias and the values from Koss et al. (2017). We adopt the virial dark subtraction, flat-field correction, and correction mass estimator based on Hβ and 5100 A˚ AGN con- for charge-transfer efficiency. The cosmic ray rejection tinuuum luminosity as calibrated by Ho & Kim (2015) was performed when combining the two dithered im- for all bulge types, using spectral measurements from ages is imperfect. Therefore, we additionally applied the Koss et al. (2017) of 56 type 1 AGN. We use the scaling LACosmic algorithm (van Dokkum 2001) to the individ- relations of Greene & Ho (2005) to convert Hα and Hβ ual images to improve the cosmic ray rejection. We used luminosities to 5100 A˚ AGN luminosity, and we scale the PyRAF task TweakReg to determine a robust shift broad Hα line widths and broad Hβ line widths. between the dithered images. The final combined im- For the Type 1 AGN lacking useful measurements of ages, generated with AstroDrizzle, have a pixel scale broad emission lines but that do have available stellar Host Galaxies of Swift-BAT AGN 5