sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Koss, M.

Normalized to: Koss, M.

122 article(s) in total. 670 co-authors, from 1 to 57 common article(s). Median position in authors list is 6,0.

[1]  [pdf] - 2127047
A hard look at local, optically-selected, obscured Seyfert galaxies
Comments: Submitted to ApJ. Comments are welcome
Submitted: 2020-07-06
We study the X-ray spectra of a sample of 19 obscured, optically-selected Seyfert galaxies (Sy 1.8, 1.9 and 2) in the local universe ($d \leq 175$~Mpc), drawn from the CfA Seyfert sample. Our analysis is driven by the high sensitivity of NuSTAR in the hard X-rays, coupled with soft X-ray spectra using XMM-Newton, Chandra, Suzaku, and Swift/XRT. We also analyze the optical spectra of these sources in order to obtain accurate mass estimates and Eddington fractions. We employ four different models to analyze the X-ray spectra of these sources, which all result in consistent results. We find that 79-90 % of the sources are heavily obscured with line-of-sight column density $N_{\rm H} > 10^{23}~\rm cm^{-2}$. We also find a Compton-thick ($N_{\rm H} > 10^{24}~\rm cm^{-2}$) fraction of $37-53$ %. These results are consistent with previous estimates based on multi-wavelength analyses. We find that the fraction of reprocessed to intrinsic emission is positively correlated with $N_{\rm H}$ and negatively correlated with the intrinsic, unabsorbed, X-ray luminosity (in agreement with the Iwasawa-Taniguchi effect). Our results support the hypothesis that radiation pressure regulates the distribution of the circumnuclear material.
[2]  [pdf] - 2121557
$NuSTAR$ Observations of Four Nearby X-ray Faint AGN: Low Luminosity or Heavy Obscuration?
Comments: 18 pages, 11 figures, 3 tables, accepted for publication in MNRAS
Submitted: 2020-06-24
We present $NuSTAR$ observations of four active galactic nuclei (AGN) located within 15 Mpc. These AGN, namely ESO 121-G6, NGC 660, NGC 3486 and NGC 5195, have observed X-ray luminosities of $L_{\rm 2-10\ keV, obs} \lesssim$ 10$^{39}$ erg s$^{-1}$, classifying them as low luminosity AGN (LLAGN). We perform broadband X-ray spectral analysis for the AGN by combining our $NuSTAR$ data with $Chandra$ or $XMM-Newton$ observations to directly measure their column densities ($N_{\rm H}$) and infer their intrinsic power. We complement our X-ray data with archival and new high angular resolution mid-infrared (mid-IR) data for all objects, except NGC 5195. Based on our X-ray spectral analysis, we found that both ESO 121-G6 and NGC 660 are heavily obscured ($N_{\rm H}$ > 10$^{23}$ cm$^{-2}$; $L_{\rm 2-10\ keV,\ int} \sim$ 10$^{41}$ erg s$^{-1}$), and NGC 660 may be Compton-thick. We also note that the X-ray flux and spectral slope for ESO 121-G6 have significantly changed over the last decade, indicating significant changes in the obscuration and potentially accretion rate. On the other hand, NGC 3486 and NGC 5195 appear to be unobscured and just mildly obscured, respectively, with $L_{\rm 2-10\ keV,\ int} <$ 10$^{39}$ erg s$^{-1}$; i.e., genuine LLAGN. Both of the heavily obscured AGN have $L_{\rm bol} >$ 10$^{41}$ erg s$^{-1}$ and $\lambda_{\rm Edd} \gtrsim$ 10$^{-3}$, and are detected in high angular resolution mid-IR imaging, indicating the presence of obscuring dust on nuclear scale. NGC 3486 however, is undetected in high-resolution mid-IR imaging, and the current data do not provide stringent constraints on the presence or absence of obscuring nuclear dust in the AGN.
[3]  [pdf] - 2124578
The BAT AGN Spectroscopic Survey -- XVIII. Searching for Supermassive Black Hole Binaries in the X-rays
Comments: 17 pages, including 8 figures and 4 tables. Accepted for publication in ApJ
Submitted: 2019-12-05, last modified: 2020-05-22
Theory predicts that a supermassive black hole binary (SMBHB) could be observed as a luminous active galactic nucleus (AGN) that periodically varies on the order of its orbital timescale. In X-rays, periodic variations could be caused by mechanisms including relativistic Doppler boosting and shocks. Here we present the first systematic search for periodic AGNs using $941$ hard X-ray light curves (14-195 keV) from the first 105 months of the Swift Burst Alert Telescope (BAT) survey (2004-2013). We do not find evidence for periodic AGNs in Swift-BAT, including the previously reported SMBHB candidate MCG+11$-$11$-$032. We find that the null detection is consistent with the combination of the upper-limit binary population in AGNs in our adopted model, their expected periodic variability amplitudes, and the BAT survey characteristics. We have also investigated the detectability of SMBHBs against normal AGN X-ray variability in the context of the eROSITA survey. Under our assumptions of a binary population and the periodic signals they produce which have long periods of hundreds of days, up to $13$% true periodic binaries can be robustly distinguished from normal variable AGNs with the ideal uniform sampling. However, we demonstrate that realistic eROSITA sampling is likely to be insensitive to long-period binaries because longer observing gaps reduce their detectability. In contrast, large observing gaps do not diminish the prospect of detecting binaries of short, few-day periods, as 19% can be successfully recovered, the vast majority of which can be identified by the first half of the survey.
[4]  [pdf] - 2119843
Characterization of the Nucleus, Morphology and Activity of Interstellar Comet 2I/Borisov by Optical and Near-Infrared GROWTH, Apache Point, IRTF, ZTF and Keck Observations
Comments: 22 pages, 7 figures, 1 table. Accepted for publication in AJ on 12 May 2020
Submitted: 2019-10-30, last modified: 2020-05-12
We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.
[5]  [pdf] - 2063856
Ionized outflows in local luminous AGN: what are the real densities and outflow rates?
Comments: Submitted to MNRAS; 25 pages (the last 5 are extra figures in an appendix)
Submitted: 2020-03-13
We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates.
[6]  [pdf] - 2030722
BAT AGN Spectroscopic Survey -- XV: The High Frequency Radio Cores of Ultra-hard X-ray Selected AGN
Comments: 22 pages, 10 figures, 2 tables. Accepted for publication in MNRAS
Submitted: 2020-01-03, last modified: 2020-01-08
We have conducted 22 GHz radio imaging at 1" resolution of 100 low-redshift AGN selected at 14-195 keV by the Swift-BAT. We find a radio core detection fraction of 96%, much higher than lower-frequency radio surveys. Of the 96 radio-detected AGN, 55 have compact morphologies, 30 have morphologies consistent with nuclear star formation, and 11 have sub-kpc to kpc-scale jets. We find that the total radio power does not distinguish between nuclear star formation and jets as the origin of the radio emission. For 87 objects, we use optical spectroscopy to test whether AGN physical parameters are distinct between radio morphological types. We find that X-ray luminosities tend to be higher if the 22 GHz morphology is jet-like, but find no significant difference in other physical parameters. We find that the relationship between the X-ray and core radio luminosities is consistent with the $L_R/L_X \sim 10^{-5}$ of coronally active stars. We further find that the canonical fundamental planes of black hole activity systematically over-predict our radio luminosities, particularly for objects with star formation morphologies.
[7]  [pdf] - 2057842
The Molecular Gas in the NGC 6240 Merging Galaxy System at the Highest Spatial Resolution
Comments: 27 pages, 13 figures, accepted for publication by The Astrophysical Journal
Submitted: 2020-01-02
We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.
[8]  [pdf] - 2050306
LLAMA: The $M_{BH}$ - $\sigma_{\star}$ Relation of the most luminous local AGNs
Comments: Accepted by A&A, 30 pages, 19 figures
Submitted: 2019-12-16
The $M_{BH}$ - $\sigma_{\star}$ relation is considered a result of co-evolution between the host galaxies and their super-massive black holes. For elliptical-bulge hosting inactive galaxies, this relation is well established, but there is still a debate whether active galaxies follow the same relation. In this paper, we estimate black hole masses for a sample of 19 local luminous AGNs (LLAMA) in order to test their location on the $M_{BH}$ - $\sigma_{\star}$ relation. Super-massive black hole masses ($M_{BH}$) were derived from the broad-line based relations for H$\alpha$, H$\beta$ and Pa$\beta$ emission line profiles for the Type 1 AGNs. We compare the bulge stellar velocity dispersion ($\sigma_{\star}$) as determined from the Ca II triplet (CaT) with the dispersion measured from the near-infrared CO (2-0) absorption features for each AGN and find them to be consistent with each other. We apply an extinction correction to the observed broad line fluxes and we correct the stellar velocity dispersion by an average rotation contribution as determined from spatially resolved stellar kinematic maps. The H$\alpha$-based black hole masses of our sample of AGNs were estimated in the range 6.34 $\leq$ $\log{M_{BH}}$ $\leq$ 7.75 M$_\odot$ and the $\sigma_{\star CaT}$ estimates range between 73 $\leq$ $\sigma_{\star CaT}$ $\leq$ 227 km s$^{-1}$. From the so-constructed $M_{BH}$ - $\sigma_{\star}$ relation for our Type 1 AGNs, we estimate the black hole masses for the Type 2 AGNs and the inactive galaxies in our sample. In conclusion, we find that our sample of local luminous AGNs is consistent with the $M_{BH}$ - $\sigma_{\star}$ relation of lower luminosity AGNs and inactive galaxies, after correcting for dust extinction and the rotational contribution to the stellar velocity dispersion.
[9]  [pdf] - 2046354
The CO(3-2)/CO(1-0) luminosity line ratio in nearby star-forming galaxies and AGN from xCOLD GASS, BASS and SLUGS
Submitted: 2019-12-02
We study the r31=L'CO(3-2)/L'CO(1-0) luminosity line ratio in a sample of nearby (z < 0.05) galaxies: 25 star-forming galaxies (SFGs) from the xCOLD GASS survey, 36 hard X-ray selected AGN host galaxies from BASS and 37 infrared luminous galaxies from SLUGS. We find a trend for r31 to increase with star-formation efficiency (SFE). We model r31 using the UCL-PDR code and find that the gas density is the main parameter responsible for variation of r31, while the interstellar radiation field and cosmic ray ionization rate play only a minor role. We interpret these results to indicate a relation between SFE and gas density. We do not find a difference in the r31 value of SFGs and AGN host galaxies, when the galaxies are matched in SSFR (<r31>= 0.52 +/- 0.04 for SFGs and <r31> = 0.53 +/- 0.06 for AGN hosts). According to the results of UCL-PDR models, the X-rays can contribute to the enhancement of the CO line ratio, but only for strong X-ray fluxes and for high gas density (nH > 10$^4$ cm-3). We find a mild tightening of the Kennicutt-Schmidt relation when we use the molecular gas mass surface density traced by CO(3-2) (Pearson correlation coefficient R=0.83), instead of the molecular gas mass surface density traced by CO(1-0) (R=0.78), but the increase in correlation is not statistically significant (p-value=0.06). This suggests that the CO(3-2) line can be reliably used to study the relation between SFR and molecular gas for normal SFGs at high redshift, and to compare it with studies of low-redshift galaxies, as is common practice.
[10]  [pdf] - 2026453
BAT AGN Spectroscopic Survey -- XIX: Type 1 versus Type 2 AGN dichotomy from the point of view of ionized outflows
Comments: Accepted for publication in MNRAS, 15 pages, 19 figures, 1 table
Submitted: 2019-11-27
We present a detailed study of ionized outflows in a large sample of ~650 hard X-ray detected AGN. Using optical spectroscopy from the BAT AGN Spectroscopic Survey (BASS) we are able to reveal the faint wings of the [OIII] emission lines associated with outflows covering, for the first time, an unexplored range of low AGN bolometric luminosity at low redshift (z~0.05). We test if and how the incidence and velocity of ionized outflow is related to AGN physical parameters: black hole mass, gas column density, Eddington Ratio, [OIII], X-ray, and bolometric luminosities. We find a higher occurrence of ionized outflows in type 1.9 (55%) and type 1 AGN (46%) with respect to type 2 AGN (24%). While outflows in type 2 AGN are evenly balanced between blue and red velocity offsets with respect to the [OIII] narrow component, they are almost exclusively blueshifted in type 1 and type 1.9 AGN. We observe a significant dependence between the outflow occurrence and accretion rate, which becomes relevant at high Eddington ratios (> -1.7). We interpret such behaviour in the framework of covering factor-Eddington ratio dependence. We don't find strong trends of the outflow maximum velocity with AGN physical parameters, as an increase with bolometric luminosity can be only identified when including samples of AGN at high luminosity and high redshift taken from literature.
[11]  [pdf] - 2026347
The Broadband X-ray Spectrum of the X-ray Obscured Type 1 AGN 2MASX J193013.80+341049.5
Comments: 13 pages, 9 figures. Accepted for publication in ApJ
Submitted: 2019-11-13
We present results from modeling the broadband X-ray spectrum of the Type 1 AGN 2MASX J193013.80+341049.5 using NuSTAR, Swift and archival XMM-Newton observations. We find this source to be highly X-ray obscured, with column densities exceeding 10$^{23}$ cm$^{-2}$ across all epochs of X-ray observations, spanning an 8 year period. However, the source exhibits prominent broad optical emission lines, consistent with an unobscured Type 1 AGN classification. We fit the X-ray spectra with both phenomenological reflection models and physically-motivated torus models to model the X-ray absorption. We examine the spectral energy distribution of this source and investigate some possible scenarios to explain the mismatch between X-ray and optical classifications. We compare the ratio of reddening to X-ray absorbing column density ($E_{B-V}/N_{H}$) and find that 2MASX J193013.80+341049.5 likely has a much lower dust-to-gas ratio relative to the Galactic ISM, suggesting that the Broad Line Region (BLR) itself could provide the source of extra X-ray obscuration, being composed of low-ionization, dust-free gas.
[12]  [pdf] - 2026345
A Dissection of Spatially Resolved AGN Feedback across the Electromagnetic Spectrum
Comments: Accepted for publication in ApJ. 20 pages, 20 figures
Submitted: 2019-11-13
We present optical SuperNova Integral Field Spectrograph (SNIFS) integral field spectroscopy, Hubble Space Telescope optical imaging, Chandra X-ray imaging, and Very Large Array radio interferometry of the merging galaxy 2MASX J04234080+0408017, which hosts a Seyfert 2 active galactic nucleus (AGN) at z = 0.046. Our observations reveal that radiatively driven, ionized gas outflows are successful to distances > 10 kpc due to the low mass of the host system, encompassing the entirety of the observed optical emission. We also find that at large radii, where observed velocities cannot be reproduced by radiative driving models, high velocity kinematics are likely due to mechanical driving from AGN winds impacting high density host material. This impacting deposits sufficient energy to shock the host material, producing thermal X-ray emission and cosmic rays, which in turn promote the formation of in situ radio structure in a pseudo-jet morphology along the high density lanes.
[13]  [pdf] - 2026132
The nature of the broadband X-ray variability in the dwarf Seyfert galaxy NGC 4395
Comments: Accepted for publication in ApJ
Submitted: 2019-10-24, last modified: 2019-10-30
We present a flux-resolved X-ray analysis of the dwarf Seyfert 1.8 galaxy NGC 4395, based on three archival $XMM-Newton$ and one archival $NuSTAR$ observations. The source is known to harbor a low mass black hole ($\sim 10^4- {\rm a~ few~}\times 10^{5}~\rm M_\odot$) and shows strong variability in the full X-ray range during these observations. We model the flux-resolved spectra of the source assuming three absorbing layers: neutral, mildly ionized, and highly ionized ($N_{\rm H} \sim 1.6\times 10^{22}-3.4 \times 10^{23}~\rm cm^{-2}$, $\sim 0.8-7.8 \times 10^{22}~\rm cm^{-2}$, and $ 3.8 \times 10^{22}~\rm cm^{-2}$, respectively. The source also shows intrinsic variability by a factor of $\sim 3$, on short timescales, due to changes in the nuclear flux, assumed to be a power law ($\Gamma = 1.6-1.67$). Our results show a positive correlation between the intrinsic flux and the absorbers' ionization parameter. The covering fraction of the neutral absorber varies during the first $XMM-Newton$ observation, which could explain the pronounced soft X-ray variability. However, the source remains fully covered by this layer during the other two observations, largely suppressing the soft X-ray variability. This suggests an inhomogeneous and layered structure in the broad line region. We also find a difference in the characteristic timescale of the power spectra between different energy ranges and observations. We finally show simulated spectra with $XRISM$, $Athena$, and $eXTP$, which will allow us to characterize the different absorbers, study their dynamics, and will help us identify their locations and sizes.
[14]  [pdf] - 1983882
The multiphase gas structure and kinematics in the circumnuclear region of NGC 5728
Comments: accepted for publication in MNRAS; 22 pages + Appendix, 16 figures
Submitted: 2019-07-08, last modified: 2019-10-02
We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar bar which is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with HII regions. We model the kinematics of the ring using the velocity field of the CO (2--1) emission and stars and find it is consistent with a rotating disk. The outer regions of the disk, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$_{\sun}$ yr$^{-1}$. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disk, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$_{\sun}$ yr$^{-1}$ consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region.
[15]  [pdf] - 1953748
BAT AGN Spectroscopic Survey -- XIII. The nature of the most luminous obscured AGN in the low-redshift universe
Comments: Accepted for publication in MNRAS (ID: stz2309)
Submitted: 2019-08-20
We present a multi wavelength analysis of 28 of the most luminous low-redshift narrow-line, ultra-hard X-ray selected active galactic nuclei (AGN) drawn from the 70 month Swift/BAT all-sky survey, with bolometric luminosities of log(L_bol/erg/s) > 45.25. The broad goal of our study is to determine whether these objects have any distinctive properties, potentially setting them aside from lower-luminosity obscured AGN in the local Universe. Our analysis relies on the first data release of the BAT AGN Spectroscopic Survey (BASS/DR1) and on dedicated observations with the VLT, Palomar, and Keck observatories. We find that the vast majority of our sources agree with commonly used AGN selection criteria which are based on emission line ratios and on mid-infrared colours. Our AGN are predominantly hosted in massive galaxies (9.8 < log(M_*/M_sun) < 11.7); based on visual inspection of archival optical images, they appear to be mostly ellipticals. Otherwise, they do not have distinctive properties. Their radio luminosities, determined from publicly available survey data, show a large spread of almost 4 orders of magnitude - much broader than what is found for lower X-ray luminosity obscured AGN in BASS. Moreover, our sample shows no preferred combination of black hole masses (M_BH) and/or Eddington ratio (lambda_Edd), covering 7.5 < log(M_BH/M_sun) < 10.3 and 0.01 < lambda_Edd < 1. Based on the distribution of our sources in the lambda_Edd-N_H plane, we conclude that our sample is consistent with a scenario where the amount of obscuring material along the line of sight is determined by radiation pressure exerted by the AGN on the dusty circumnuclear gas.
[16]  [pdf] - 1920174
BAT AGN Spectroscopic Survey -- XVII: The Parsec-scale Jet Properties of the Ultra Hard X-ray Selected Local AGN
Comments: 11+2 pages, 6 figures, Accepted for publication in MNRAS, Comments are welcome
Submitted: 2019-07-19
We have performed a very long baseline interferometry (VLBI) survey of local (z < 0.05) ultra hard X-ray (14-195 keV) selected active galactic nuclei (AGN) from the Swift Burst Alert Telescope (BAT) using KVN, KaVA, and VLBA. We first executed fringe surveys of 142 BAT-detected AGN at 15 or 22 GHz. Based on the fringe surveys and archival data, we find 10/279 nearby AGN (~4%) VLBI have 22 GHz flux above 30 mJy. This implies that the X-ray AGN with a bright nuclear jet are not common. Among these 10 radio-bright AGN, we obtained 22 GHz VLBI imaging data of our own for four targets and reprocessed archival data for six targets. We find that, although our 10 AGN observed with VLBI span a wide range of pc-scale morphological types, they lie on a tight linear relation between accretion luminosity and nuclear jet luminosity. Our result suggests that a powerful nuclear radio jet correlates with the accretion disc luminosity. We also probed the fundamental plane of black hole activity at VLBI scales (e.g., few milli-arcsecond). The jet luminosity and size distribution among our sample roughly fit into the proposed AGN evolutionary scenario, finding powerful jets after the blow-out phase based on the Eddington ratio (\lambda_{Edd})-hydrogen column density (N_{H}) relation. In addition, we find some hints of gas inflow or galaxy-galaxy merger in the majority of our sample. This implies that gas supply via tidal interactions in galactic scale may help the central AGN to launch a powerful parsec-scale jet.
[17]  [pdf] - 1958191
Measuring the Obscuring Column of a Disk Megamaser AGN in a Nearby Merger
Comments: 10 pages, 8 figures. Accepted for publication in ApJ
Submitted: 2019-07-11
Active Galactic Nuclei (AGNs) hosting disk water megamasers are well known to be obscured by large amounts of gas, likely due to the presence along the line of sight of an almost edge-on disky structure orbiting the supermassive black hole. Correcting for the high obscuration is crucial to infer parameters intrinsic to the source, like its luminosity. We present a broadband X-ray spectral analysis of a water megamaser AGN in an early merger (NGC 5765B), combining Chandra and NuSTAR data. NGC 5765B is highly Compton-thick and reflection-dominated, following the general trend among disk megamasers. Combining the exquisite black hole mass from masers with our X-ray spectroscopy, the Eddington ratio of the megamaser is estimated to be in the $2-14\%$ range, and its robustness is confirmed through SED fitting.
[18]  [pdf] - 1948804
BAT AGN Spectroscopic Survey: XVI. General Physical Characteristics of BAT Blazars
Comments: 16 pages, 7 figures, 6 tables, to appear in the Astrophysical Journal. Full data tables and SED plots are included in the arXiv uploaded files
Submitted: 2019-07-01
The recently released 105-month {\it Swift}-Burst Alert Telescope (BAT) all-sky hard X-ray survey catalog presents an opportunity to study astrophysical objects detected in the deepest look at the entire hard X-ray (14$-$195 keV) sky. Here we report the results of a multifrequency study of 146 blazars from this catalog, quadrupling the number compared to past studies, by utilizing recent data from the {\it Fermi}-Large Area Telescope (LAT), Swift-BAT, and archival measurements. In our $\gamma$-ray analysis of $\sim$10 years of the LAT data, 101 are found as $\gamma$-ray emitters, whereas, 45 remains LAT undetected. We model the broadband spectral energy distributions with a synchrotron-inverse Compton radiative model. On average, BAT detected sources host massive black holes ($M_{\rm bh}\sim10^9$ M$_{\odot}$) and luminous accretion disks ($L_{\rm d}\sim10^{46}$ erg s$^{-1}$). At high-redshifts ($z>2$), BAT blazars host more powerful jets with luminous accretion disks compared to those detected only with the {\it Fermi}-LAT. We find good agreement in the black hole masses derived from the single-epoch optical spectroscopic measurements and standard accretion disk modeling approaches. Other physical properties of BAT blazars are similar to those known for {\it Fermi}-LAT detected objects.
[19]  [pdf] - 1929825
An observational link between AGN Eddington ratio and [NII]{\lambda}6583/Ha at 0.6<z<1.7
Comments: 15 pages, 7 figures, and 2 tables; Accepted for publication in ApJ
Submitted: 2019-06-10
We present an observed relationship between Eddington ratio ({\lambda}Edd) and optical narrow-emission-line ratio ([N II]{\lambda}6583/H{\alpha}) of X-ray-selected broad-line active galactic nuclei (AGN) at 0.6 < z < 1.7. We use 27 near-infrared spectra from the Fiber Multi-Object Spectrograph along with 26 sources from the literature. We show that the {\lambda}Edd and [N II]{\lambda}6583/H{\alpha} ratio at 0.6 < z < 1.7 exhibits a similar anti-correlation distribution of {\lambda}Edd-[N II]{\lambda}6583/H{\alpha} as has been found for local (<z> = 0.036), ultra-hard X-ray selected AGN. The observed distribution implies that there is a consistent relationship from local to z \sim 1.7 which corresponds from the present time to 4 Gyr old. Further study of high redshift low Eddington ratio AGN (log{\lambda}Edd < -2) is necessary to determine fully whether the {\lambda}Edd-[N II]{\lambda}6583/H{\alpha} anti-correlation still holds in high-redshift AGN at low Eddington ratios.
[20]  [pdf] - 1894331
A hard look at NGC 5347: revealing a nearby Compton-thick AGN
Comments: Accepted for publication in ApJ
Submitted: 2019-04-24
Current measurements show that the observed fraction of Compton-thick (CT) AGN is smaller than the expected values needed to explain the cosmic X-ray background. Prior fits to the X-ray spectrum of the nearby Seyfert-2 galaxy NGC 5347 ($z=0.00792,\, D =35.5 \rm ~Mpc $) have alternately suggested a CT and Compton-thin source. Combining archival data from $Suzaku$, $Chandra$, and - most importantly - new data from $NuSTAR$, and using three distinct families of models, we show that NGC 5347 is an obscured CTAGN ($N_{\rm H} > 2.23\times 10^{24}~\rm cm^{-2}$). Its 2-30~keV spectrum is dominated by reprocessed emission from distant material, characterized by a strong Fe K$\alpha$ line and a Compton hump. We found a large equivalent width of the Fe K$\alpha$ line ($\rm EW = 2.3 \pm 0.3$ keV) and a high intrinsic-to-observed flux ratio ($\sim 100$). All of these observations are typical for bona fide CTAGN. We estimate a bolometric luminosity of $L_{\rm bol} \simeq 0.014 \pm 0.005~L_{\rm Edd.}$. The $Chandra$ image of NGC 5347 reveals the presence of extended emission dominating the soft X-ray spectrum ($E < 2\,\rm keV$), which coincides with the [O III] emission detected in the $Hubble ~Space~ Telescope$ images. Comparison to other CTAGN suggests that NGC 5347 is broadly consistent with the average properties of this source class. We simulated $XRISM$ and $Athena$/X-IFU spectra of the source, showing the potential of these future missions in identifying CTAGN in the soft X-rays.
[21]  [pdf] - 1863882
Investigating the Covering Fraction Distribution of Swift/BAT AGN with X-ray and IR Observations
Comments: Published in the Astrophysical Journal: Lanz, L. et al. 2019, ApJ, 870, 26; 19 pages; 9 figures
Submitted: 2018-11-06, last modified: 2019-04-08
We present an analysis of a sample of 69 local obscured Swift/Burst Alert Telescope active galactic nuclei (AGNs) with X-ray spectra from NuSTAR and infrared (IR) spectral energy distributions from Herschel and WISE. We combine this X-ray and IR phenomenological modeling and find a significant correlation between reflected hard X-ray emission and IR AGN emission, with suggestive indications that this correlation may be stronger than the one between intrinsic hard X-ray and IR emissions. This relation between the IR and reflected X-ray emission suggests that both are the result of the processing of intrinsic emission from the corona and accretion disk by the same structure. We explore the resulting implications on the underlying distribution of covering fraction for all AGNs, by generating mock observables for the reflection parameter and IR luminosity ratio using empirical relations found for the covering fraction with each quantity. We find that the observed distributions of the reflection parameter and IR-to-X-ray ratio are reproduced with broad distributions centered around covering fractions of at least ~40%-50%, whereas narrower distributions match our observations only when centered around covering fractions of ~70%-80%. Our results are consistent with both independent estimates of the covering fractions of individual objects and the typical covering fraction obtained on the basis of obscured fractions for samples of AGNs. These results suggest that the level of reprocessing in AGNs, including X-ray reflection, is related in a relatively straightforward way to the geometry of the obscuring material.
[22]  [pdf] - 1850842
Astro2020 Science White Paper: Black Hole Growth in Mergers and Dual AGN
Comments: Submission to the 2020-2030 Astronomy and Astrophysics Decadal Survey (Astro2020), 10 pages, 3 figures
Submitted: 2019-03-15
Hierarchical models of galaxy formation predict that galaxy mergers represent a significant transitional stage of rapid supermassive black hole (SMBH) growth. Yet, the connection between the merging process and enhanced active galactic nuclei (AGN) activity as well as the timescale of SMBH mergers remains highly uncertain. The breakthrough in reconciling the importance of galaxy mergers with black hole growth lies in a thoroughly-studied census of dual AGN across cosmic history, which will be enabled by next-generation observational capabilities, theoretical advances, and simulations. This white paper outlines the key questions in galaxy mergers, dual and offset AGN, and proposes multiwavelength solutions using future high-resolution observatories in the X-rays (AXIS, Lynx), near and mid-infrared (30 meter class telescopes, JWST), and submillimeter (ALMA).
[23]  [pdf] - 1851763
The Advanced X-ray Imaging Satellite
Comments: A Probe-class mission study commissioned by NASA for the NAS Astro2020 Decadal Survey. Cost section redacted. 66 pages, 41 figures. v2: minor fixes
Submitted: 2019-03-10, last modified: 2019-03-14
Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline.
[24]  [pdf] - 1845076
STROBE-X: X-ray Timing and Spectroscopy on Dynamical Timescales from Microseconds to Years
Ray, Paul S.; Arzoumanian, Zaven; Ballantyne, David; Bozzo, Enrico; Brandt, Soren; Brenneman, Laura; Chakrabarty, Deepto; Christophersen, Marc; DeRosa, Alessandra; Feroci, Marco; Gendreau, Keith; Goldstein, Adam; Hartmann, Dieter; Hernanz, Margarita; Jenke, Peter; Kara, Erin; Maccarone, Tom; McDonald, Michael; Nowak, Michael; Phlips, Bernard; Remillard, Ron; Stevens, Abigail; Tomsick, John; Watts, Anna; Wilson-Hodge, Colleen; Wood, Kent; Zane, Silvia; Ajello, Marco; Alston, Will; Altamirano, Diego; Antoniou, Vallia; Arur, Kavitha; Ashton, Dominic; Auchettl, Katie; Ayres, Tom; Bachetti, Matteo; Balokovic, Mislav; Baring, Matthew; Baykal, Altan; Begelman, Mitch; Bhat, Narayana; Bogdanov, Slavko; Briggs, Michael; Bulbul, Esra; Bult, Petrus; Burns, Eric; Cackett, Ed; Campana, Riccardo; Caspi, Amir; Cavecchi, Yuri; Chenevez, Jerome; Cherry, Mike; Corbet, Robin; Corcoran, Michael; Corsi, Alessandra; Degenaar, Nathalie; Drake, Jeremy; Eikenberry, Steve; Enoto, Teruaki; Fragile, Chris; Fuerst, Felix; Gandhi, Poshak; Garcia, Javier; Goldstein, Adam; Gonzalez, Anthony; Grefenstette, Brian; Grinberg, Victoria; Grossan, Bruce; Guillot, Sebastien; Guver, Tolga; Haggard, Daryl; Heinke, Craig; Heinz, Sebastian; Hemphill, Paul; Homan, Jeroen; Hui, Michelle; Huppenkothen, Daniela; Ingram, Adam; Irwin, Jimmy; Jaisawal, Gaurava; Jaodand, Amruta; Kalemci, Emrah; Kaplan, David; Keek, Laurens; Kennea, Jamie; Kerr, Matthew; van der Klis, Michiel; Kocevski, Daniel; Koss, Mike; Kowalski, Adam; Lai, Dong; Lamb, Fred; Laycock, Silas; Lazio, Joseph; Lazzati, Davide; Longcope, Dana; Loewenstein, Michael; Maitra, Dipankair; Majid, Walid; Maksym, W. Peter; Malacaria, Christian; Margutti, Raffaella; Martindale, Adrian; McHardy, Ian; Meyer, Manuel; Middleton, Matt; Miller, Jon; Miller, Cole; Motta, Sara; Neilsen, Joey; Nelson, Tommy; Noble, Scott; O'Brien, Paul; Osborne, Julian; Osten, Rachel; Ozel, Feryal; Palliyaguru, Nipuni; Pasham, Dheeraj; Patruno, Alessandro; Pelassa, Vero; Petropoulou, Maria; Pilia, Maura; Pohl, Martin; Pooley, David; Prescod-Weinstein, Chanda; Psaltis, Dimitrios; Raaijmakers, Geert; Reynolds, Chris; Riley, Thomas E.; Salvesen, Greg; Santangelo, Andrea; Scaringi, Simone; Schanne, Stephane; Schnittman, Jeremy; Smith, David; Smith, Krista Lynne; Snios, Bradford; Steiner, Andrew; Steiner, Jack; Stella, Luigi; Strohmayer, Tod; Sun, Ming; Tauris, Thomas; Taylor, Corbin; Tohuvavohu, Aaron; Vacchi, Andrea; Vasilopoulos, Georgios; Veledina, Alexandra; Walsh, Jonelle; Weinberg, Nevin; Wilkins, Dan; Willingale, Richard; Wilms, Joern; Winter, Lisa; Wolff, Michael; Zand, Jean in 't; Zezas, Andreas; Zhang, Bing; Zoghbi, Abdu
Comments: 50 pages, Probe class mission concept study report submitted to NASA for Astro2020 Decadal Survey
Submitted: 2019-03-07, last modified: 2019-03-08
We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.
[25]  [pdf] - 1838314
A deep X-ray view of the bare AGN Ark120. V. Spin determination from disc-Comptonisation efficiency method
Comments: 15 pages, 7 figures. Accepted for publication in A&A
Submitted: 2019-01-07
[Abridged] In our previous work on Ark 120, we found that its 2014 X-ray spectrum is dominated by Comptonisation, while the relativistic reflection emission only originates at tens of $R_{\rm g}$ from the SMBH. As a result, we could not constrain the SMBH spin from disc reflection alone. Our aim is to determine its SMBH spin from an alternative technique based on the global energetics of the disc-corona system. The spectral analysis uses simultaneous XMM-Newton (OM and pn) and NuSTAR observations on 2014 March 22 and 2013 February 18. We applied the optxconv model (based on optxagnf) to self consistently reproduce the emission from the inner corona (warm and hot thermal Comptonisation) and the outer disc (colour temperature corrected black body), taking into account both the disc inclination angle and relativistic effects. We modelled the mild relativistic reflection of the incident Comptonisation components using the xilconv model. We infer a SMBH spin of 0.83$^{+0.05}_{-0.03}$, adopting the SMBH reverberation mass of 1.50$\times$10$^{8}$ M$_{\odot}$. In addition, we find that the coronal radius decreases with increasing flux (by about a factor of two), from 85$^{+13}_{-10}$ $R_{\rm g}$ in 2013 to 14$\pm$3 $R_{\rm g}$ in 2014. This is the first time that such a constraint is obtained for a SMBH spin from this technique, thanks to the bare properties of Ark 120, its well determined SMBH mass, and the presence of a mild relativistic reflection component in 2014 which allows us to constrain the disc inclination angle. We caution that these results depend on the detailed disc-corona structure, which is not yet fully established. However, the realistic parameter values found suggest that this is a promising method to determine spin in moderate accretion rate AGN.
[26]  [pdf] - 1782223
A population of luminous accreting black holes with hidden mergers
Comments: To appear in the 8 November 2018 issue of Nature. This is the authors' version of the work
Submitted: 2018-11-08
Major galaxy mergers are thought to play an important part in fuelling the growth of supermassive black holes. However, observational support for this hypothesis is mixed, with some studies showing a correlation between merging galaxies and luminous quasars and others showing no such association. Recent observations have shown that a black hole is likely to become heavily obscured behind merger-driven gas and dust, even in the early stages of the merger, when the galaxies are well separated (5 to 40 kiloparsecs). Merger simulations further suggest that such obscuration and black-hole accretion peaks in the final merger stage, when the two galactic nuclei are closely separated (less than 3 kiloparsecs). Resolving this final stage requires a combination of high-spatial-resolution infrared imaging and high-sensitivity hard-X-ray observations to detect highly obscured sources. However, large numbers of obscured luminous accreting supermassive black holes have been recently detected nearby (distances below 250 megaparsecs) in X-ray observations. Here we report high-resolution infrared observations of hard-X-ray-selected black holes and the discovery of obscured nuclear mergers, the parent populations of supermassive-black-hole mergers. We find that obscured luminous black holes (bolometric luminosity higher than 2x10^44 ergs per second) show a significant (P<0.001) excess of late-stage nuclear mergers (17.6 per cent) compared to a sample of inactive galaxies with matching stellar masses and star formation rates (1.1 per cent), in agreement with theoretical predictions. Using hydrodynamic simulations, we confirm that the excess of nuclear mergers is indeed strongest for gas-rich major-merger hosts of obscured luminous black holes in this final stage.
[27]  [pdf] - 1811184
NuSTAR and Keck Observations of Heavily Obscured Quasars Selected by WISE
Comments: Accepted for publication in the Astrophysical Journal, 11 pages, 10 figures
Submitted: 2018-11-06
A primary aim of the Nuclear Spectroscopic Telescope Array (NuSTAR) mission is to find and characterize heavily obscured Active Galactic Nuclei (AGNs). Based on mid-infrared photometry from the Wide-Field Infrared Survey Explorer (WISE) and optical photometry from the Sloan Digital Sky Surveys, we have selected a large population of luminous obscured AGN (i.e., obscured quasars). Here we report NuSTAR observations of four WISE-selected heavily obscured quasars for which we have optical spectroscopy from the Southern African Large Telescope and W. M. Keck Observatory. Optical diagnostics confirm that all four targets are AGNs. With NuSTAR hard X-ray observations, three of the four objects are undetected, while the fourth has a marginal detection. We confirm that these objects have observed hard X-ray (10-40 keV) luminosities at or below ~ 10^43 erg s^-1. We compare X-ray and IR luminosities to obtain estimates of the hydrogen column densities (N_H) based on the suppression of the hard X-ray emission. We estimate N_H of these quasars to be at or larger than 10^25 cm^-2, confirming that WISE and optical selection can identify very heavily obscured quasars that may be missed in X-ray surveys, and do not contribute significantly to the cosmic X-ray background. From the optical Balmer decrements, we found that our three extreme obscured targets lie in highly reddened host environments. This galactic extinction is not adequate to explain the more obscured AGN, but it may imply a different scale of obscuration in the galaxy.
[28]  [pdf] - 1811183
BAT AGN Spectroscopic Survey - XI. The Covering Factor of Dust and Gas in Swift/BAT Active Galactic Nuclei
Comments: 21 pages, 15 figures, accepted for publication in ApJ. The full list of Table 1 is available at
Submitted: 2018-11-06
We quantify the luminosity contribution of active galactic nuclei (AGN) to the 12 $\mu$m, mid-infrared (MIR; 5-38 $\mu$m), and the total IR (5-1000 $\mu$m) emission in the local AGN detected in the all-sky 70-month Swift/Burst Alert Telescope (BAT) ultra hard X-ray survey. We decompose the IR spectral energy distributions (SEDs) of 587 objects into AGN and starburst components using AGN torus and star-forming galaxy templates. This enables us to recover the AGN torus emission also for low-luminosity end, down to $\log (L_{14-150}/{\rm erg}~{\rm s}^{-1}) \simeq 41$, which typically have significant host galaxy contamination. We find that the luminosity contribution of the AGN to the 12 $\mu$m, the MIR, and the total IR band is an increasing function of the 14-150 keV luminosity. We also find that for the most extreme cases, the IR pure-AGN emission from the torus can extend up to 90 $\mu$m. The obtained total IR AGN luminosity through the IR SED decomposition enables us to estimate the fraction of the sky obscured by dust, i.e., the dust covering factor. We demonstrate that the median of the dust covering factor is always smaller than that of the X-ray obscuration fraction above the AGN bolometric luminosity of $\log (L_{\rm bol}/{\rm erg}~{\rm s}^{-1}) \simeq 42.5$. Considering that X-ray obscuration fraction is equivalent to the covering factor coming from both the dust and gas, it indicates that an additional neutral gas component, along with the dusty torus, is responsible for the absorption of X-ray emission.
[29]  [pdf] - 1783780
A long hard-X-ray look at the dual active galactic nuclei of M51 with NuSTAR
Comments: Accepted for publication in ApJ
Submitted: 2018-05-30, last modified: 2018-09-13
We present a broadband X-ray spectral analysis of the M51 system, including the dual active galactic nuclei (AGN) and several off-nuclear point sources. Using a deep observation by NuSTAR, new high-resolution coverage of M51b by Chandra, and the latest X-ray torus models, we measure the intrinsic X-ray luminosities of the AGN in these galaxies. The AGN of M51a is found to be Compton thick, and both AGN have very low accretion rates ($\lambda_{\rm Edd} <10^{-4}$). The latter is surprising considering that the galaxies of M51 are in the process of merging, which is generally predicted to enhance nuclear activity. We find that the covering factor of the obscuring material in M51a is $0.26 \pm 0.03$, consistent with the local AGN obscured fraction at $L_{\rm X}\sim 10^{40}$ erg s$^{-1}$. The substantial obscuring column does not support theories that the torus, presumed responsible for the obscuration, disappears at these low accretion luminosities. However, the obscuration may have resulted from the gas infall driven by the merger rather than the accretion process. We report on several extra-nuclear sources with $L_{\rm X}>10^{39}$ erg s$^{-1}$ and find that a spectral turnover is present below 10 keV in most such sources, in line with recent results on ultraluminous X-ray sources.
[30]  [pdf] - 1752012
BAT AGN Spectroscopic Survey -- XII. The relation between coronal properties of Active Galactic Nuclei and the Eddington ratio
Comments: Accepted for publication in MNRAS
Submitted: 2018-09-11
The bulk of the X-ray emission in Active Galactic Nuclei (AGN) is produced very close to the accreting supermassive black hole (SMBH), in a corona of hot electrons which up scatters optical and ultraviolet photons from the accretion flow. The cutoff energy ($E_{\rm C}$) of the primary X-ray continuum emission carries important information on the physical characteristics of the X-ray emitting plasma, but little is currently known about its potential relation with the properties of accreting SMBHs. Using the largest broad-band (0.3-150 keV) X-ray spectroscopic study available to date, we investigate how the corona is related to the AGN luminosity, black hole mass and Eddington ratio ($\lambda_{\rm Edd}$). Assuming a slab corona the median values of the temperature and optical depth of the Comptonizing plasma are $kT_{\rm e}=105 \pm 18$ keV and $\tau=0.25\pm0.06$, respectively. When we properly account for the large number of $E_{\rm C}$ lower limits, we find a statistically significant dependence of the cutoff energy on the Eddington ratio. In particular, objects with $ \lambda_{\rm Edd}>0.1$ have a significantly lower median cutoff energy ($E_{\rm C}=160\pm41$ keV) than those with $\lambda_{\rm Edd}\leq 0.1$ ($E_{\rm C}=370\pm51$ keV). This is consistent with the idea that radiatively compact coronae are also cooler, because they tend to avoid the region in the temperature-compactness parameter space where runaway pair production would dominate. We show that this behaviour could also straightforwardly explain the suggested positive correlation between the photon index ($\Gamma$) and the Eddington ratio, being able to reproduce the observed slope of the $\Gamma-\lambda_{\rm Edd}$ trend.
[31]  [pdf] - 1682626
The Swift/BAT AGN Spectroscopic Survey -- IX. The Clustering Environments of an Unbiased Sample of Local AGN
Comments: Accepted for publication in ApJ
Submitted: 2018-03-20, last modified: 2018-04-26
We characterize the environments of local accreting supermassive black holes by measuring the clustering of AGN in the Swift/BAT Spectroscopic Survey (BASS). With 548 AGN in the redshift range 0.01<z<0.1 over the full sky from the DR1 catalog, BASS provides the largest, least biased sample of local AGN to date due to its hard X-ray selection (14-195 keV) and rich multiwavelength/ancillary data. By measuring the projected cross-correlation function between the AGN and 2MASS galaxies, and interpreting it via halo occupation distribution (HOD) and subhalo-based models, we constrain the occupation statistics of the full sample, as well as in bins of absorbing column density and black hole mass. We find that AGN tend to reside in galaxy group environments, in agreement with previous studies of AGN throughout a large range of luminosity and redshift, and that on average they occupy their dark matter halos similar to inactive galaxies of comparable stellar mass. We also find evidence that obscured AGN tend to reside in denser environments than unobscured AGN, even when samples were matched in luminosity, redshift, stellar mass, and Eddington ratio. We show that this can be explained either by significantly different halo occupation distributions or statistically different host halo assembly histories. Lastly, we see that massive black holes are slightly more likely to reside in central galaxies than black holes of smaller mass.
[32]  [pdf] - 1656404
PSFGAN: a generative adversarial network system for separating quasar point sources and host galaxy light
Comments: 17 pages, 18 figures, accepted for publication in MNRAS
Submitted: 2018-03-23
The study of unobscured active galactic nuclei (AGN) and quasars depends on the reliable decomposition of the light from the AGN point source and the extended host galaxy light. The problem is typically approached using parametric fitting routines using separate models for the host galaxy and the point spread function (PSF). We present a new approach using a Generative Adversarial Network (GAN) trained on galaxy images. We test the method using Sloan Digital Sky Survey (SDSS) r-band images with artificial AGN point sources added which are then removed using the GAN and with parametric methods using GALFIT. When the AGN point source PS is more than twice as bright as the host galaxy, we find that our method, PSFGAN, can recover PS and host galaxy magnitudes with smaller systematic error and a lower average scatter ($49\%$). PSFGAN is more tolerant to poor knowledge of the PSF than parametric methods. Our tests show that PSFGAN is robust against a broadening in the PSF width of $\pm 50\%$ if it is trained on multiple PSF's. We demonstrate that while a matched training set does improve performance, we can still subtract point sources using a PSFGAN trained on non-astronomical images. While initial training is computationally expensive, evaluating PSFGAN on data is more than $40$ times faster than GALFIT fitting two components. Finally, PSFGAN it is more robust and easy to use than parametric methods as it requires no input parameters.
[33]  [pdf] - 1648895
A model for AGN variability on multiple timescales
Comments: 5 pages, 2 figures, letter accepted for publication in MNRAS
Submitted: 2018-02-15, last modified: 2018-03-13
We present a framework to link and describe AGN variability on a wide range of timescales, from days to billions of years. In particular, we concentrate on the AGN variability features related to changes in black hole fuelling and accretion rate. In our framework, the variability features observed in different AGN at different timescales may be explained as realisations of the same underlying statistical properties. In this context, we propose a model to simulate the evolution of AGN light curves with time based on the probability density function (PDF) and power spectral density (PSD) of the Eddington ratio ($L/L_{\rm Edd}$) distribution. Motivated by general galaxy population properties, we propose that the PDF may be inspired by the $L/L_{\rm Edd}$ distribution function (ERDF), and that a single (or limited number of) ERDF+PSD set may explain all observed variability features. After outlining the framework and the model, we compile a set of variability measurements in terms of structure function (SF) and magnitude difference. We then combine the variability measurements on a SF plot ranging from days to Gyr. The proposed framework enables constraints on the underlying PSD and the ability to link AGN variability on different timescales, therefore providing new insights into AGN variability and black hole growth phenomena.
[34]  [pdf] - 1641442
Optical, near-IR and sub-mm IFU Observations of the nearby dual AGN Mrk 463
Comments: Accepted by The Astrophysical Journal, 23 pages, 19 figures
Submitted: 2018-01-18
We present optical and near-IR Integral Field Unit (IFU) and ALMA band 6 observations of the nearby dual Active Galactic Nuclei (AGN) Mrk 463. At a distance of 210 Mpc, and a nuclear separation of $\sim$4 kpc, Mrk 463 is an excellent laboratory to study the gas dynamics, star formation processes and supermassive black hole (SMBH) accretion in a late-stage gas-rich major galaxy merger. The IFU observations reveal a complex morphology, including tidal tails, star-forming clumps, and emission line regions. The optical data, which map the full extent of the merger, show evidence for a biconical outflow and material outflowing at $>$600 km s$^{-1}$, both associated with the Mrk 463E nucleus, together with large scale gradients likely related to the ongoing galaxy merger. We further find an emission line region $\sim$11 kpc south of Mrk 463E that is consistent with being photoionized by an AGN. Compared to the current AGN luminosity, the energy budget of the cloud implies a luminosity drop in Mrk 463E by a factor 3-20 over the last 40,000 years. The ALMA observations of $^{12}$CO(2-1) and adjacent 1mm continuum reveal the presence of $\sim$10$^{9}$M$_\odot$ in molecular gas in the system. The molecular gas shows velocity gradients of $\sim$800 km/s and $\sim$400 km/s around the Mrk 463E and 463W nuclei, respectively. We conclude that in this system the infall of $\sim$100s $M_\odot$/yr of molecular gas is in rough balance with the removal of ionized gas by a biconical outflow being fueled by a relatively small, $<$0.01% of accretion onto each SMBH.
[35]  [pdf] - 1648783
The 105 month Swift-BAT all-sky hard X-ray survey
Comments: Accepted for publication in ApJS. The Swift-BAT 105-month Survey public website can be found at this URL:
Submitted: 2018-01-05
We present a catalog of hard X-ray sources detected in the first 105 months of observations with the Burst Alert Telescope (BAT) coded mask imager on board the Swift observatory. The 105 month Swift-BAT survey is a uniform hard X-ray all-sky survey with a sensitivity of $8.40\times 10^{-12}\ {\rm erg\ s^{-1}\ cm^{-2}}$ over 90% of the sky and $7.24\times 10^{-12}\ {\rm erg\ s^{-1}\ cm^{-2}}$ over 50% of the sky in the 14-195 keV band. The Swift-BAT 105 month catalog provides 1632 (422 new detections) hard X-ray sources in the 14-195 keV band above the 4.8{\sigma} significance level. Adding to the previously known hard X-ray sources, 34% (144/422) of the new detections are identified as Seyfert AGN in nearby galaxies (z<0.2). The majority of the remaining identified sources are X-ray binaries (7%, 31) and blazars/BL Lac objects (10%, 43). As part of this new edition of the Swift-BAT catalog, we release eight-channel spectra and monthly sampled light curves for each object in the online journal and at the Swift-BAT 105 month Web site.
[36]  [pdf] - 1604909
BAT AGN Spectroscopic Survey - V. X-ray properties of the Swift/BAT 70-month AGN catalog
Comments: ApJS in press. The full tables will be available on ApJS and the on the BASS website [] once the paper is published
Submitted: 2017-09-12, last modified: 2017-11-28
Hard X-ray ($\geq 10$ keV) observations of Active Galactic Nuclei (AGN) can shed light on some of the most obscured episodes of accretion onto supermassive black holes. The 70-month Swift/BAT all-sky survey, which probes the 14-195 keV energy range, has currently detected 838 AGN. We report here on the broad-band X-ray (0.3-150 keV) characteristics of these AGN, obtained by combining XMM-Newton, Swift/XRT, ASCA, Chandra, and Suzaku observations in the soft X-ray band ($\leq 10$ keV) with 70-month averaged Swift/BAT data. The non-blazar AGN of our sample are almost equally divided into unobscured ($N_{\rm H}< 10^{22}\rm cm^{-2}$) and obscured ($N_{\rm H}\geq 10^{22}\rm cm^{-2}$) AGN, and their Swift/BAT continuum is systematically steeper than the 0.3-10 keV emission, which suggests that the presence of a high-energy cutoff is almost ubiquitous. We discuss the main X-ray spectral parameters obtained, such as the photon index, the reflection parameter, the energy of the cutoff, neutral and ionized absorbers, and the soft excess for both obscured and unobscured AGN.
[37]  [pdf] - 1591353
Joint NuSTAR and Chandra analysis of the obscured quasar in IC 2497 - Hanny's Voorwerp system
Comments: 9 pages, 3 figures. Accepted for publication in MNRAS
Submitted: 2017-11-16, last modified: 2017-11-22
We present new Nuclear Spectroscopic Telescope Array (NuSTAR) observations of the core of IC 2497, the galaxy associated with Hanny's Voorwerp. The combined fits of the Chandra (0.5-8 keV) and NuSTAR (3-24 keV) X-ray spectra, together with WISE mid-IR photometry, optical longslit spectroscopy and optical narrow-band imaging, suggest that the galaxy hosts a Compton-thick AGN ($N_{\rm H} \sim 2 \times 10^{24}$ cm$^{-2}$, current intrinsic luminosity $L_{\rm bol} \sim 2-5 \times 10^{44}$ erg s$^{-1}$) whose luminosity dropped by a factor of $\sim$50 within the last $\sim 100$ kyr. This corresponds to a change in Eddington ratio from $\rm \lambda_{Edd} \sim$ 0.35 to $\rm \lambda_{Edd} \sim$ 0.007. We argue that the AGN in IC 2497 should not be classified as a changing-look AGN, but rather we favour the interpretation where the AGN is undergoing a change in accretion state (from radiatively efficient to radiatively inefficient). In this scenario the observed drop in luminosity and Eddington ratio corresponds to the final stage of an AGN accretion phase. Our results are consistent with previous studies in the optical, X-ray and radio although the magnitude of the drop is lower than previously suggested. In addition, we discuss a possible analogy between X-ray binaries and an AGN.
[38]  [pdf] - 1666761
BAT AGN Spectroscopic Survey--VIII. Type 1 AGN With Massive Absorbing Columns
Comments: 14 pages, 10 figures, 1 Appendix, submitted to ApJ
Submitted: 2017-10-25
We explore the relationship between X-ray absorption and optical obscuration within the BAT AGN Spectroscopic Survey (BASS) which has been collecting and analyzing the optical and X-ray spectra for 641 hard X-ray selected ($E>14$ keV) active galactic nuclei (AGN). We use the deviation from a linear broad H$\alpha$-to-X-ray relationship as an estimate of the maximum optical obscuration towards the broad line region and compare the $A_{\rm V}$ to the hydrogen column densities ($N_{\rm H}$) found through systematic modeling of their X-ray spectra. We find that the inferred columns implied by $A_{\rm V}$ towards the broad line region (BLR) are often orders of magnitude less than the columns measured towards the X-ray emitting region indicating a small scale origin for the X-ray absorbing gas. After removing 30\% of Sy 1.9s that potentially have been misclassified due to outflows, we find that 86\% (164/190) of the Type 1 population (Sy 1--1.9) are X-ray unabsorbed as expected based on a single obscuring structure. However, 14\% (26/190), of which 70\% (18/26) are classified as Sy 1.9, are X-ray absorbed, suggesting the broad line region itself is providing extra obscuration towards the X-ray corona. The fraction of X-ray absorbed Type 1 AGN remains relatively constant with AGN luminosity and Eddington ratio, indicating a stable broad line region covering fraction.
[39]  [pdf] - 1589584
LLAMA: Nuclear stellar properties of Swift BAT AGN and matched inactive galaxies
Comments: 34 pages, 25 figures. Accepted for publication in MNRAS
Submitted: 2017-10-11
In a complete sample of local 14-195 keV selected AGNs and inactive galaxies, matched by their host galaxy properties, we study the spatially resolved stellar kinematics and luminosity distributions at near-infrared wavelengths on scales of 10-150 pc, using SINFONI on the VLT. In this paper, we present the first half of the sample, which comprises 13 galaxies, 8 AGNs and 5 inactive galaxies. The stellar velocity fields show a disk-like rotating pattern, for which the kinematic position angle is in agreement with the photometric position angle obtained from large scale images. For this set of galaxies, the stellar surface brightness of the inactive galaxy sample is generally comparable to the matched sample of AGN but extends to lower surface brightness. After removal of the bulge contribution, we find a nuclear stellar light excess with an extended nuclear disk structure, and which exhibits a size-luminosity relation. While we expect the excess luminosity to be associated with a dynamically cooler young stellar population, we do not typically see a matching drop in dispersion. This may be because these galaxies have pseudo-bulges in which the intrinsic dispersion increases towards the centre. And although the young stars may have an impact in the observed kinematics, their fraction is too small to dominate over the bulge and compensate the increase in dispersion at small radii, so no dispersion drop is seen. Finally, we find no evidence for a difference in the stellar kinematics and nuclear stellar luminosity excess between these active and inactive galaxies.
[40]  [pdf] - 1604978
LLAMA: Normal star formation efficiencies of molecular gas in the centres of luminous Seyfert galaxies
Comments: 24 pages, 18 figures, 2 part Appendix; accepted for publication in MNRAS
Submitted: 2017-10-11
Using new APEX and JCMT spectroscopy of the CO 2-1 line, we undertake a controlled study of cold molecular gas in moderately luminous Active Galactic Nuclei (AGN) and inactive galaxies from the Luminous Local AGN with Matched Analogs (LLAMA) survey. We use spatially resolved infrared photometry of the LLAMA galaxies from 2MASS, WISE, IRAS & Herschel, corrected for nuclear emission using multi-component spectral energy distribution (SED) fits, to examine the dust-reprocessed star-formation rates (SFRs), molecular gas fractions and star formation efficiencies (SFEs) over their central 1 - 3 kpc. We find that the gas fractions and central SFEs of both active and inactive galaxies are similar when controlling for host stellar mass and morphology (Hubble type). The equivalent central molecular gas depletion times are consistent with the discs of normal spiral galaxies in the local Universe. Despite energetic arguments that the AGN in LLAMA should be capable of disrupting the observable cold molecular gas in their central environments, our results indicate that nuclear radiation only couples weakly with this phase. We find a mild preference for obscured AGN to contain higher amounts of central molecular gas, which suggests a connection between AGN obscuration and the gaseous environment of the nucleus. Systems with depressed SFEs are not found among the LLAMA AGN. We speculate that the processes that sustain the collapse of molecular gas into dense pre-stellar cores may also be a prerequisite for the inflow of material on to AGN accretion disks.
[41]  [pdf] - 1602676
Investigating the evolution of the dual AGN system ESO~509-IG066
Comments: Accepted for publication in ApJ
Submitted: 2017-10-09
We analyze the evolution of the dual AGN in ESO 509-IG066, a galaxy pair located at z=0.034 whose nuclei are separated by 11 kpc. Previous observations with XMM-Newton on this dual AGN found evidence for two moderately obscured ($N_H\sim10^{22}$ cm$^{-2}$) X-ray luminous ($L_X\sim10^{43}$ erg/s) nuclear sources. We present an analysis of subsequent Chandra, NuSTAR and Swift/XRT observations that show one source has dropped in flux by a factor of 10 between 2004 and 2011, which could be explained by either an increase in the absorbing column or an intrinsic fading of the central engine possibly due to a decrease in mass accretion. Both of these scenarios are predicted by galaxy merger simulations. The source which has dropped in flux is not detected by NuSTAR, which argues against absorption, unless it is extreme. However, new Keck/LRIS optical spectroscopy reveals a previously unreported broad H-alpha line which is highly unlikely to be visible under the extreme absorption scenario. We therefore conclude that the black hole in this nucleus has undergone a dramatic drop in accretion rate. From AO-assisted near-infrared integral-field spectroscopy of the other nucleus, we find evidence that the galaxy merger is having a direct effect on the kinematics of the gas close to the nucleus of the galaxy, providing a direct observational link between the galaxy merger and the mass accretion rate on to the black hole.
[42]  [pdf] - 1588983
The close environments of accreting massive black holes are shaped by radiative feedback
Comments: To appear in the 28 September 2017 issue of Nature. This is the authors' version of the work
Submitted: 2017-09-27
The large majority of the accreting supermassive black holes in the Universe are obscured by large columns of gas and dust. The location and evolution of this obscuring material have been the subject of intense research in the past decades, and are still highly debated. A decrease in the covering factor of the circumnuclear material with increasing accretion rates has been found by studies carried out across the electromagnetic spectrum. The origin of this trend has been suggested to be driven either by the increase in the inner radius of the obscuring material with incident luminosity due to the sublimation of dust; by the gravitational potential of the black hole; by radiative feedback; or by the interplay between outflows and inflows. However, the lack of a large, unbiased and complete sample of accreting black holes, with reliable information on gas column density, luminosity and mass, has left the main physical mechanism regulating obscuration unclear. Using a systematic multi-wavelength survey of hard X-ray-selected black holes, here we show that radiation pressure on dusty gas is indeed the main physical mechanism regulating the distribution of the circumnuclear material. Our results imply that the bulk of the obscuring dust and gas in these objects is located within the sphere of influence of the black hole (i.e., a few to tens of parsecs), and that it can be swept away even at low radiative output rates. The main physical driver of the differences between obscured and unobscured accreting black holes is therefore their mass-normalized accretion rate.
[43]  [pdf] - 1587362
NuSTAR Resolves the First Dual AGN above 10 keV in SWIFT J2028.5+2543
Comments: Published in ApJ, 7 pages, 4 figures
Submitted: 2017-08-22
We have discovered heavy obscuration in the dual active galactic nucleus (AGN) in the Swift/Burst Alert Telescope (BAT) source SWIFT J2028.5+2543 using Nuclear Spectroscopic Telescope Array (NuSTAR). While an early XMM-Newton study suggested the emission was mainly from NGC 6921, the superior spatial resolution of NuSTAR above 10 keV resolves the Swift/BAT emission into two sources associated with the nearby galaxies MCG +04-48-002 and NGC 6921 (z = 0.014) with a projected separation of 25.3 kpc (91"). NuSTAR's sensitivity above 10 keV finds both are heavily obscured to Compton-thick levels (N H=(1-2)x10^24 cm-2) and contribute equally to the BAT detection ({L}10-50 {keV}{{int}}= 6x10^42 erg s-1). The observed luminosity of both sources is severely diminished in the 2-10 keV band, illustrating the importance of >10 keV surveys like those with NuSTAR and Swift/BAT. Compared to archival X-ray data, MCG +04-48-002 shows significant variability (>3) between observations. Despite being bright X-ray AGNs, they are difficult to detect using optical emission-line diagnostics because MCG +04-48-002 is identified as a starburst/composite because of the high rates of star formation from a luminous infrared galaxy while NGC 6921 is only classified as a LINER using line detection limits. SWIFT J2028.5+2543 is the first dual AGN resolved above 10 keV and is the second most heavily obscured dual AGN discovered to date in the X-rays other than NGC 6240.
[44]  [pdf] - 1586197
The NuSTAR Serendipitous Survey: Hunting for The Most Extreme Obscured AGN at >10 keV
Comments: 19 pages, 15 figures, 5 tables; Accepted for publication in ApJ; Author list updated
Submitted: 2017-07-20, last modified: 2017-08-16
We identify sources with extremely hard X-ray spectra (i.e., with photon indices of Gamma<0.6 in the 13 sq. deg. NuSTAR serendipitous survey, to search for the most highly obscured AGNs detected at >10 keV. Eight extreme NuSTAR sources are identified, and we use the NuSTAR data in combination with lower energy X-ray observations (from Chandra, Swift XRT, and XMM-Newton) to characterize the broad-band (0.5-24 keV) X-ray spectra. We find that all of the extreme sources are highly obscured AGNs, including three robust Compton-thick (CT; N_H > 1.5e24 cm^-2) AGNs at low redshift (z<0.1), and a likely-CT AGN at higher redshift (z=0.16). Most of the extreme sources would not have been identified as highly obscured based on the low energy (<10 keV) X-ray coverage alone. The multiwavelength properties (e.g., optical spectra and X-ray/MIR luminosity ratios) provide further support for the eight sources being significantly obscured. Correcting for absorption, the intrinsic rest-frame 10-40 keV luminosities of the extreme sources cover a broad range, from ~ 5 x 10^42 to 10^45 erg s^-1. The estimated number counts of CT AGNs in the NuSTAR serendipitous survey are in broad agreement with model expectations based on previous X-ray surveys, except for the lowest redshifts (z<0.07) where we measure a high CT fraction of f_CT^obs = 30 (+16 -12) %. For the small sample of CT AGNs, we find a high fraction of galaxy major mergers (50 +/- 33%) compared to control samples of "normal" AGNs.
[45]  [pdf] - 1593566
BAT AGN Spectroscopic Survey I: Spectral Measurements, Derived Quantities, and AGN Demographics
Comments: Accepted ApJ, see for data
Submitted: 2017-07-25
We present the first catalog and data release of the Swift-BAT AGN Spectroscopic Survey (BASS). We analyze optical spectra of the majority of AGN (77%, 641/836) detected based on their 14-195 keV emission in the 70-month Swift BAT all-sky catalog. This includes redshift determination, absorption and emission line measurements, and black hole mass and accretion rate estimates for the majority of obscured and un-obscured AGN (74%, 473/641) with 340 measured for the first time. With ~90% of sources at z<0.2, the survey represents a significant census of hard-X-ray selected AGN in the local universe. In this first catalog paper, we describe the spectroscopic observations and datasets, and our initial spectral analysis. The FWHM of the emission lines show broad agreement with the X-ray obscuration (~94%), such that Sy 1-1.8 have NH<10^21.9 cm^-2, and Seyfert 2, have NH>10^21.9 cm^-2. Seyfert 1.9 show a range of column densities. Compared to narrow line AGN in the SDSS, the X-ray selected AGN have a larger fraction of dusty host galaxies suggesting these types of AGN are missed in optical surveys. Using the most sensitive [OIII]/Hbeta and [NII]/Halpha emission line diagnostic, about half of the sources are classified as Seyferts, ~15% reside in dusty galaxies that lack an Hbeta detection, but for which the line upper limits imply either a Seyfert or LINER, ~15% are in galaxies with weak or no emission lines despite high quality spectra, and a few percent each are LINERS, composite galaxies, HII regions, or in known beamed AGN.
[46]  [pdf] - 1584942
Inferring Compton-thick AGN candidates at z>2 with Chandra using the >8 keV restframe spectral curvature
Comments: 10 pages, 8 figures
Submitted: 2017-06-20
To fully understand cosmic black hole growth we need to constrain the population of heavily obscured active galactic nuclei (AGN) at the peak of cosmic black hole growth ($z\sim$1-3). Sources with obscuring column densities higher than $\mathrm{10^{24}}$ atoms $\mathrm{cm^{-2}}$, called Compton-thick (CT) AGN, can be identified by excess X-ray emission at $\sim$20-30 keV, called the "Compton hump". We apply the recently developed Spectral Curvature (SC) method to high-redshift AGN (2<z<5) detected with Chandra. This method parametrizes the characteristic "Compton hump" feature cosmologically redshifted into the X-ray band at observed energies <10 keV. We find good agreement in CT AGN found using the SC method and bright sources fit using their full spectrum with X-ray spectroscopy. In the Chandra deep field south, we measure a CT fraction of $\mathrm{17^{+19}_{-11}\%}$ (3/17) for sources with observed luminosity $\mathrm{>5\times 10^{43}}$ erg $\mathrm{s^{-1}}$. In the Cosmological evolution survey (COSMOS), we find an observed CT fraction of $\mathrm{15^{+4}_{-3}\%}$ (40/272) or $\mathrm{32\pm11 \%}$ when corrected for the survey sensitivity. When comparing to low redshift AGN with similar X-ray luminosities, our results imply the CT AGN fraction is consistent with having no redshift evolution. Finally, we provide SC equations that can be used to find high-redshift CT AGN (z>1) for current (XMM-Newton) and future (eROSITA and ATHENA) X-ray missions.
[47]  [pdf] - 1579951
The weak neutral Fe fluorescence line and long-term X-ray evolution of the Compton-thick AGN in NGC 7674
Comments: MNRAS accepted version
Submitted: 2016-05-25, last modified: 2017-06-15
We present $NuSTAR$ X-ray observations of the active galactic nucleus (AGN) in NGC 7674. The source shows a flat X-ray spectrum, suggesting that it is obscured by Compton-thick gas columns. Based upon long-term flux dimming, previous work suggested the alternate possibility that the source is a recently switched-off AGN with the observed X-rays being the lagged echo from the torus. Our high-quality data show the source to be reflection-dominated in hard X-rays, but with a relatively weak neutral Fe K$\alpha$ emission line (equivalent width [EW] of $\approx$ 0.4 keV) and a strong Fe XXVI ionised line (EW $\approx$ 0.2 keV). We construct an updated long-term X-ray light curve of NGC 7674 and find that the observed 2-10 keV flux has remained constant for the past $\approx$ 20 years, following a high flux state probed by $Ginga$. Light travel time arguments constrain the minimum radius of the reflector to be $\sim$ 3.2 pc under the switched-off AGN scenario, $\approx$ 30 times larger than the expected dust sublimation radius, rendering this possibility unlikely. A patchy Compton-thick AGN (CTAGN) solution is plausible, requiring a minimum line-of-sight column density ($N_{\rm H}$) of 3 $\times$ 10$^{24}$ cm$^{-2}$ at present, and yields an intrinsic 2-10 keV luminosity of (3-5) $\times$ 10$^{43}$ erg s$^{-1}$. Realistic uncertainties span the range of $\approx$ (1-13) $\times$ 10$^{43}$ erg s$^{-1}$. The source has one of the weakest fluorescence lines amongst {\em bona fide} CTAGN, and is potentially a local analogue of bolometrically luminous systems showing complex neutral and ionised Fe emission. It exemplifies the difficulty of identification and proper characterisation of distant CTAGN based on the strength of the neutral Fe K$\alpha$ line.
[48]  [pdf] - 1584057
X-ray bolometric corrections for Compton-thick active galactic nuclei
Comments: Accepted for publication in ApJ
Submitted: 2017-05-30
We present X-ray bolometric correction factors, $\kappa_{Bol}$ ($\equiv L_{Bol}/L_X$), for Compton-thick (CT) active galactic nuclei (AGN) with the aim of testing AGN torus models, probing orientation effects, and estimating the bolometric output of the most obscured AGN. We adopt bolometric luminosities, $L_{Bol}$, from literature infrared (IR) torus modeling and compile published intrinsic 2--10 keV X-ray luminosities, $L_{X}$, from X-ray torus modeling of NuSTAR data. Our sample consists of 10 local CT AGN where both of these estimates are available. We test for systematic differences in $\kappa_{Bol}$ values produced when using two widely used IR torus models and two widely used X-ray torus models, finding consistency within the uncertainties. We find that the mean $\kappa_{Bol}$ of our sample in the range $L_{Bol}\approx10^{42}-10^{45}$ erg/s is log$_{10}\kappa_{Bol}=1.44\pm0.12$ with an intrinsic scatter of $\sim0.2$ dex, and that our derived $\kappa_{Bol}$ values are consistent with previously established relationships between $\kappa_{Bol}$ and $L_{Bol}$ and $\kappa_{Bol}$ and Eddington ratio. We investigate if $\kappa_{Bol}$ is dependent on $N_H$ by comparing our results on CT AGN to published results on less-obscured AGN, finding no significant dependence. Since many of our sample are megamaser AGN, known to be viewed edge-on, and furthermore under the assumptions of AGN unification whereby unobscured AGN are viewed face-on, our result implies that the X-ray emitting corona is not strongly anisotropic. Finally, we present $\kappa_{Bol}$ values for CT AGN identified in X-ray surveys as a function of their observed $L_X$, where an estimate of their intrinsic $L_{X}$ is not available, and redshift, useful for estimating the bolometric output of the most obscured AGN across cosmic time.
[49]  [pdf] - 1583083
The Swift/BAT AGN Spectroscopic Survey (BASS) -- VI. The Gamma_X - L/L_Edd relation
Comments: Resubmitted to MNRAS after addressing referee's comments
Submitted: 2017-05-03
We study the observed relation between accretion rate (in terms of L/L_Edd) and shape of the hard X-ray spectral energy distribution (namely the photon index Gamma_X) for a large sample of 228 hard X-ray selected, low-redshift active galactic nuclei (AGN), drawn from the Swift/BAT AGN Spectroscopic Survey (BASS). This includes 30 AGN for which black hole mass (and therefore L/L_Edd) is measured directly through masers, spatially resolved gas or stellar dynamics, or reverberation mapping. The high quality and broad energy coverage of the data provided through BASS allow us to examine several alternative determinations of both Gamma_X and L/L_Edd. For the BASS sample as a whole, we find a statistically significant, albeit very weak correlation between Gamma_X and L/L_Edd. The best-fitting relations we find, Gamma_X=0.15 log(L/L_Edd)+const., are considerably shallower than those reported in previous studies. Moreover, we find no corresponding correlations among the subsets of AGN with different M_BH determination methodology. In particular, we find no robust evidence for a correlation when considering only those AGN with direct or single-epoch M_BH estimates. This latter finding is in contrast to several previous studies which focused on z>0.5 broad-line AGN. We discuss this tension and conclude that it can be partially accounted for if one adopts a simplified, power-law X-ray spectral model, combined with L/L_Edd estimates that are based on the continuum emission and on single-epoch broad line spectroscopy in the optical regime. We finally highlight the limitations on using Gamma_X as a probe of supermassive black hole evolution in deep extragalactic X-ray surveys.
[50]  [pdf] - 1542942
Hard X-ray selected AGNs in low-mass galaxies from the NuSTAR serendipitous survey
Comments: Accepted for publication in ApJ. 27 pages, 6 tables, 10 figures, including one online figure set available at (v2: updated acknowledgement.)
Submitted: 2017-01-30, last modified: 2017-02-22
We present a sample of 10 low-mass active galactic nuclei (AGNs) selected from the 40-month NuSTAR serendipitous survey. The sample is selected to have robust NuSTAR detections at $3 - 24$~keV, to be at $z < 0.3$, and to have optical r-band magnitudes at least 0.5~mag fainter than an $L_\star$ galaxy at its redshift. The median values of absolute magnitude, stellar mass and 2--10 X-ray luminosity of our sample are $\langle M_r\rangle = -20.03$, $\langle M_\star\rangle = 4.6\times10^{9}M_\odot$, and $\langle L_{2-10\mathrm{keV}}\rangle = 3.1\times10^{42}$ erg s$^{-1}$, respectively. Five objects have detectable broad H$\alpha$ emission in their optical spectra, indicating black-hole masses of $(1.1-10.4)\times 10^6 M_\odot$. We find that $30^{+17}_{-10}\%$ of the galaxies in our sample do not show AGN-like optical narrow emission lines, and one of the ten galaxies in our sample, J115851+4243.2, shows evidence for heavy X-ray absorption. This result implies that a non-negligible fraction of low-mass galaxies might harbor accreting massive black holes that are missed by optical spectroscopic surveys and $<10$ keV X-ray surveys. The mid-IR colors of our sample also indicate these optically normal low-mass AGNs cannot be efficiently identified with typical AGN selection criteria based on WISE colors. While the hard ($>10$ keV) X-ray selected low-mass AGN sample size is still limited, our results show that sensitive NuSTAR observations are capable of probing faint hard X-ray emission originating from the nuclei of low-mass galaxies out to moderate redshift ($z<0.3$), thus providing a critical step in understanding AGN demographics in low-mass galaxies.
[51]  [pdf] - 1542949
Evidence for Relativistic Disk Reflection in the Seyfert 1h Galaxy/ULIRG IRAS 05189-2524 Observed by NuSTAR and XMM-Newton
Comments: 8 pages, 5 figures. Accepted for publication in ApJ
Submitted: 2017-01-31
We present a spectral analysis of the NuSTAR and XMM-Newton observations of the Seyfert 1h galaxy/ULIRG IRAS 05189-2524 taken in 2013. We find evidence for relativistic disk reflection in the broadband X-ray spectrum: a highly asymmetric broad Fe K$\alpha$ emission line extending down to 3 keV and a Compton scattering component above 10 keV. Physical modeling with a self-consistent disk reflection model suggests the accretion disk is viewed at an intermediate angle with a super-solar iron abundance, and a mild constraint can be put on the high-energy cutoff of the power-law continuum. We test the disk reflection modeling under different absorption scenarios. A rapid black hole (BH) spin is favored, however we cannot place a model-independent tight constraint on the value. The high reflection fraction ($R_{\rm ref} \simeq$ 2.0-3.2) suggests the coronal illuminating source is compact and close to the BH (lying within 8.7 $R_{\rm g}$ above the central BH), where light-bending effects are important.
[52]  [pdf] - 1534394
BAT AGN Spectroscopic Survey - IV: Near-Infrared Coronal Lines, Hidden Broad Lines, and Correlation with Hard X-ray Emission
Comments: 38 pages, 19 figures. Accepted for publication in MNRAS
Submitted: 2017-01-10
We provide a comprehensive census of the near-Infrared (NIR, 0.8-2.4 $\mu$m) spectroscopic properties of 102 nearby (z < 0.075) active galactic nuclei (AGN), selected in the hard X-ray band (14-195 keV) from the Swift-Burst Alert Telescope (BAT) survey. With the launch of the James Webb Space Telescope this regime is of increasing importance for dusty and obscured AGN surveys. We measure black hole masses in 68% (69/102) of the sample using broad emission lines (34/102) and/or the velocity dispersion of the Ca II triplet or the CO band-heads (46/102). We find that emission line diagnostics in the NIR are ineffective at identifying bright, nearby AGN galaxies because ([Fe II] 1.257$\mu$m/Pa$\beta$ and H$_2$ 2.12$\mu$m/Br$\gamma$) identify only 25% (25/102) as AGN with significant overlap with star forming galaxies and only 20% of Seyfert 2 have detected coronal lines (6/30). We measure the coronal line emission in Seyfert 2 to be weaker than in Seyfert 1 of the same bolometric luminosity suggesting obscuration by the nuclear torus. We find that the correlation between the hard X-ray and the [Si VI] coronal line luminosity is significantly better than with the [O III] luminosity. Finally, we find 3/29 galaxies (10%) that are optically classified as Seyfert 2 show broad emission lines in the NIR. These AGN have the lowest levels of obscuration among the Seyfert 2s in our sample ($\log N_{\rm H} < 22.43$ cm$^{-2}$), and all show signs of galaxy-scale interactions or mergers suggesting that the optical broad emission lines are obscured by host galaxy dust.
[53]  [pdf] - 1534312
A Long Look At MCG-5-23-16 With NuSTAR: I- Relativistic Reflection And Coronal Properties
Comments: Accepted for publication in ApJ
Submitted: 2017-01-09
MCG-5-23-16 was targeted in early 2015 with a half mega-seconds observing campaign using NuSTAR. Here we present the spectral analysis of these datasets along with an earlier observation and study the relativistic reflection and the primary coronal source. The data show strong reflection features in the form of both narrow and broad iron lines plus a Compton reflection hump. A cutoff energy is significantly detected in all exposures. The shape of the reflection spectrum does not change in the two years spanned by the observations, suggesting a stable geometry. A strong positive correlation is found between the cutoff energy and both the hard X-ray flux and spectral index. The measurements imply that the coronal plasma is not at the runaway electron-positron pair limit, and instead contains mostly electrons. The observed variability in the coronal properties is driven by a variable optical depth. A constant heating to cooling ratio is measured implying that there is a feedback mechanism in which a significant fraction of the photons cooling the corona are due to reprocessed hard X-rays.
[54]  [pdf] - 1532342
The Role of Host Galaxy for the Environmental Dependence of Active Nuclei in Local Galaxies
Comments: accepted for publication in MNRAS. 12 pages with 10 figures
Submitted: 2016-10-31, last modified: 2017-01-06
We discuss the environment of local hard X-ray selected active galaxies, with reference to two independent group catalogues. We find that the fraction of these AGN in S0 host galaxies decreases strongly as a function of galaxy group size (halo mass) - which contrasts with the increasing fraction of galaxies of S0 type in denser environments. However, there is no evidence for an environmental dependence of AGN in spiral galaxies. Because most AGN are found in spiral galaxies, this dilutes the signature of environmental dependence for the population as a whole. We argue that the differing results for AGN in disk-dominated and bulge-dominated galaxies is related to the source of the gas fuelling the AGN, and so may also impact the luminosity function, duty cycle, and obscuration. We find that there is a significant difference in the luminosity function for AGN in spiral and S0 galaxies, and tentative evidence for some difference in the fraction of obscured AGN.
[55]  [pdf] - 1534049
A New Compton-thick AGN in our Cosmic Backyard: Unveiling the Buried Nucleus in NGC 1448 with NuSTAR
Comments: Accepted for publication in ApJ on Dec. 27th, 2016; 14 pages; 9 figures and 2 tables
Submitted: 2017-01-02
NGC 1448 is one of the nearest luminous galaxies ($L_{8-1000\mu m} >$ 10$^{9} L_{\odot}$) to ours ($z$ $=$ 0.00390), and yet the active galactic nucleus (AGN) it hosts was only recently discovered, in 2009. In this paper, we present an analysis of the nuclear source across three wavebands: mid-infrared (MIR) continuum, optical, and X-rays. We observed the source with the Nuclear Spectroscopic Telescope Array (NuSTAR), and combined this data with archival Chandra data to perform broadband X-ray spectral fitting ($\approx$0.5-40 keV) of the AGN for the first time. Our X-ray spectral analysis reveals that the AGN is buried under a Compton-thick (CT) column of obscuring gas along our line-of-sight, with a column density of $N_{\rm H}$(los) $\gtrsim$ 2.5 $\times$ 10$^{24}$ cm$^{-2}$. The best-fitting torus models measured an intrinsic 2-10 keV luminosity of $L_{2-10\rm{,int}}$ $=$ (3.5-7.6) $\times$ 10$^{40}$ erg s$^{-1}$, making NGC 1448 one of the lowest luminosity CTAGNs known. In addition to the NuSTAR observation, we also performed optical spectroscopy for the nucleus in this edge-on galaxy using the European Southern Observatory New Technology Telescope. We re-classify the optical nuclear spectrum as a Seyfert on the basis of the Baldwin-Philips-Terlevich diagnostic diagrams, thus identifying the AGN at optical wavelengths for the first time. We also present high spatial resolution MIR observations of NGC 1448 with Gemini/T-ReCS, in which a compact nucleus is clearly detected. The absorption-corrected 2-10 keV luminosity measured from our X-ray spectral analysis agrees with that predicted from the optical [OIII]$\lambda$5007\AA\ emission line and the MIR 12$\mu$m continuum, further supporting the CT nature of the AGN.
[56]  [pdf] - 1533626
The NuSTAR Serendipitous Survey: The 40 month Catalog and the Properties of the Distant High Energy X-ray Source Population
Comments: 49 pages, 23 figures, 8 tables. Accepted for publication in ApJ
Submitted: 2016-12-19
We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of <z> = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy (<10 keV) X-ray counterparts from XMM-Newton, Chandra, and Swift/XRT. The mid-infrared (MIR) analysis, using WISE all-sky survey data, shows that MIR AGN color selections miss a large fraction of the NuSTAR-selected AGN population, from ~15% at the highest luminosities (Lx > 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by <10 keV X-ray missions.
[57]  [pdf] - 1533569
Broadband X-ray spectral analysis of the Seyfert 1 galaxy GRS 1734-292
Submitted: 2016-12-18
We discuss the broadband X-ray spectrum of GRS 1734-292 obtained from non-simultaneous XMM-Newton and NuSTAR observations, performed in 2009 and 2014, respectively. GRS1734-292 is a Seyfert 1 galaxy, located near the Galactic plane at $z=0.0214$. The NuSTAR spectrum ($3-80$ keV) is dominated by a primary power-law continuum with $\Gamma=1.65 \pm 0.05$ and a high-energy cutoff $E_c=53^{+11}_{-8}$ keV, one of the lowest measured by NuSTAR in a Seyfert galaxy. Comptonization models show a temperature of the coronal plasma of $kT_e=11.9^{+1.2}_{-0.9}$ keV and an optical depth, assuming a slab geometry, $\tau=2.98^{+0.16}_{-0.19}$ or a similar temperature and $\tau=6.7^{+0.3}_{-0.4}$ assuming a spherical geometry. The 2009 XMM-Newton spectrum is well described by a flatter intrinsic continuum ($\Gamma=1.47^{+0.07}_{-0.03}$) and one absorption line due to Fe\textsc{XXV} K$\alpha$ produced by a warm absorber. Both data sets show a modest iron K$\alpha$ emission line at $6.4$ keV and the associated Compton reflection, due to reprocessing from neutral circumnuclear material.
[58]  [pdf] - 1533415
Active Galactic Nuclei from HeII: a more complete census of AGN in SDSS galaxies yields a new population of low-luminosity AGN in highly star-forming galaxies
Comments: Accepted for publication in MNRAS
Submitted: 2016-12-14
In order to perform a more complete census of active galactic nuclei (AGN) in the local Universe, we investigate the use of the HeII~$\lambda 4685$ emission line diagnostic diagram by Shirazi & Brinchmann (2012) in addition to the standard methods based on other optical emission lines. The He II based diagnostics is more sensitive to AGN ionization in the presence of strong star formation than conventional line diagnostics. We survey a magnitude-limited sample of 63,915 galaxies from the Sloan Digital Sky Survey Data Release 7 at $0.02<z<0.05$ and use both the conventional BPT emission line diagnostic diagrams, as well as the HeII diagram to identify AGN. In this sample, 1,075 galaxies are selected as AGN using the BPT diagram, while an additional 234 galaxies are identified as AGN using the HeII diagnostic diagram, representing a 22% increase of AGN in the parent galaxy sample. We explore the host galaxy properties of these new HeII selected AGN candidates and find that they are most common in star-forming galaxies on the blue cloud and on the main sequence where ionization from star-formation is most likely to mask AGN emission in the BPT lines. We note in particular a high HeII AGN fraction in galaxies above the high-mass end of the main sequence where quenching is expected to occur. We use archival Chandra observations to confirm the AGN nature of candidates selected through HeII based diagnostic. Finally, we discuss how this technique can help inform galaxy/black hole co-evolution scenarios.
[59]  [pdf] - 1533302
Herschel far-infrared photometry of the Swift Burst Alert Telescope active galactic nuclei sample of the local universe - III. Global star-forming properties and the lack of a connection to nuclear activity
Comments: 26 pages, 17 figures, 3 Tables, accepted for publication in MNRAS
Submitted: 2016-12-12
We combine the Herschel Space Observatory PACS and SPIRE photometry with archival WISE photometry to construct the spectral energy distributions (SED) for over 300 local ($z < 0.05$), ultra-hard X-ray (14 - 195 keV) selected active galactic nuclei (AGN) from the Swift Burst Alert Telescope (BAT) 58 month catalogue. Using a simple analytical model that combines an exponentially cut-off powerlaw with a single temperature modified blackbody, we decompose the SEDs into a host-galaxy and AGN component. We calculate dust masses, dust temperatures, and star-formation rates (SFR) for our entire sample and compare them to a stellar mass-matched sample of local non-AGN galaxies. We find AGN host galaxies have systematically higher dust masses, dust temperatures, and SFRs due to the higher prevalence of late-type galaxies to host an AGN, in agreement with previous studies of the Swift/BAT AGN. We provide a scaling to convert X-ray luminosities into 8 - 1000 $\mu$m AGN luminosities, as well as determine the best mid-to-far IR colors for identifying AGN dominated galaxies in the IR regime. We find that for nearly 30 per cent of our sample, the 70 $\mu$m emission contains a significant contribution from the AGN ($> 0.5$), especially at higher luminosities ($L_{14-195\,\rm{keV}} > 10^{42.5}$ ergs s$^{-1}$). Finally, we measure the local SFR-AGN luminosity relationship, finding a slope of 0.18, large scatter (0.37 dex), and no evidence for an upturn at high AGN luminosity. We conclude with a discussion on the implications of our results within the context of galaxy evolution with and without AGN feedback.
[60]  [pdf] - 1532891
The Complete Infrared View of Active Galactic Nuclei from the 70-month Swift/BAT Catalog
Comments: 24 pages, 11 figures, accepted for publication in ApJ. The full list of Table 1 is available at
Submitted: 2016-11-29
We systematically investigate the near- (NIR) to far-infrared (FIR) photometric properties of a nearly complete sample of local active galactic nuclei (AGN) detected in the Swift/Burst Alert Telescope (BAT) all-sky ultra hard X-ray (14-195 keV) survey. Out of 606 non-blazar AGN in the Swift/BAT 70-month catalog at high galactic latitude of $|b|>10^{\circ}$, we obtain IR photometric data of 604 objects by cross-matching the AGN positions with catalogs from the WISE, AKARI, IRAS, and Herschel infrared observatories. We find a good correlation between the ultra-hard X-ray and mid-IR (MIR) luminosities over five orders of magnitude ($41 < \log (L_{14-195}/{\rm erg}~{\rm s}^{-1})< 46$). Informed by previous measures of the intrinsic spectral energy distribution of AGN, we find FIR pure-AGN candidates whose FIR emission is thought to be AGN-dominated with low starformation activity. We demonstrate that the dust covering factor decreases with the bolometric AGN luminosity, confirming the luminosity-dependent unified scheme. We also show that the completeness of the WISE color-color cut in selecting Swift/BAT AGN increases strongly with 14-195 keV luminosity.
[61]  [pdf] - 1531610
NuSTAR observations of WISE J1036+0449, a Galaxy at z$\sim1$ obscured by hot dust
Comments: Accepted for publication in ApJ, 16 pages
Submitted: 2016-09-15, last modified: 2016-11-21
Hot, Dust-Obscured Galaxies (Hot DOGs), selected from the WISE all sky infrared survey, host some of the most powerful Active Galactic Nuclei (AGN) known, and might represent an important stage in the evolution of galaxies. Most known Hot DOGs are at $z> 1.5$, due in part to a strong bias against identifying them at lower redshift related to the selection criteria. We present a new selection method that identifies 153 Hot DOG candidates at $z\sim 1$, where they are significantly brighter and easier to study. We validate this approach by measuring a redshift $z=1.009$, and an SED similar to higher redshift Hot DOGs for one of these objects, WISE J1036+0449 ($L_{\rm\,Bol}\simeq 8\times 10^{46}\rm\,erg\,s^{-1}$), using data from Keck/LRIS and NIRSPEC, SDSS, and CSO. We find evidence of a broadened component in MgII, which, if due to the gravitational potential of the supermassive black hole, would imply a black hole mass of $M_{\rm\,BH}\simeq 2 \times 10^8 M_{\odot}$, and an Eddington ratio of $\lambda_{\rm\,Edd}\simeq 2.7$. WISE J1036+0449 is the first Hot DOG detected by NuSTAR, and the observations show that the source is heavily obscured, with a column density of $N_{\rm\,H}\simeq(2-15)\times10^{23}\rm\,cm^{-2}$. The source has an intrinsic 2-10 keV luminosity of $\sim 6\times 10^{44}\rm\,erg\,s^{-1}$, a value significantly lower than that expected from the mid-infrared/X-ray correlation. We also find that the other Hot DOGs observed by X-ray facilities show a similar deficiency of X-ray flux. We discuss the origin of the X-ray weakness and the absorption properties of Hot DOGs. Hot DOGs at $z\lesssim1$ could be excellent laboratories to probe the characteristics of the accretion flow and of the X-ray emitting plasma at extreme values of the Eddington ratio.
[62]  [pdf] - 1532300
IC 3639 - A new bona fide Compton thick AGN unveiled by NuSTAR
Comments: Accepted for publication in the The Astrophysical Journal, 26/10/2016
Submitted: 2016-10-27
We analyse high-quality NuSTAR observations of the local (z = 0.011) Seyfert 2 active galactic nucleus (AGN) IC 3639, in conjunction with archival Suzaku and Chandra data. This provides the first broadband X-ray spectral analysis of the source, spanning nearly two decades in energy (0.5-30 keV). Previous X-ray observations of the source below 10 keV indicated strong reflection/obscuration on the basis of a pronounced iron fluorescence line at 6.4 keV. The hard X-ray energy coverage of NuSTAR, together with self-consistent toroidal reprocessing models, enables direct broadband constraints on the obscuring column density of the source. We find the source to be heavily Compton-thick (CTK) with an obscuring column in excess of $3.6\times10^{24}$ cm$^{-2}$, unconstrained at the upper end. We further find an intrinsic 2-10 keV luminosity of $\textrm{log}_{10}(L_{\textrm{2-10 keV}} \textrm{[erg s}^{-1}]) = 43.4^{+0.6}_{-1.1}$ to 90% confidence, almost 400 times the observed flux, and consistent with various multi-wavelength diagnostics. Such a high intrinsic to observed flux ratio in addition to an Fe-K$\alpha$ fluorescence line equivalent width exceeding 2 keV is extreme amongst known bona fide CTK AGN, which we suggest are both due to the high level of obscuration present around IC 3639. Our study demonstrates that broadband spectroscopic modelling with NuSTAR enables large corrections for obscuration to be carried out robustly, and emphasises the need for improved modelling of AGN tori showing intense iron fluorescence.
[63]  [pdf] - 1513717
BAT AGN Spectroscopic Survey-III. An observed link between AGN Eddington ratio and narrow emission line ratios
Comments: Accepted for publication in MNRAS, 9 pages, 5 figures, 1 table
Submitted: 2016-09-27
We investigate the observed relationship between black hole mass ($M_{\rm BH}$), bolometric luminosity ($L_{\rm bol}$), and Eddington ratio (${\lambda}_{\rm Edd}$) with optical emission line ratios ([NII] {\lambda}6583/H{\alpha}, [SII] {\lambda}{\lambda}6716,6731/H{\alpha}, [OI] {\lambda}6300/H{\alpha}, [OIII] {\lambda}5007/H{\beta}, [NeIII] {\lambda}3869/H{\beta}, and HeII {\lambda}4686/H{\beta}) of hard X-ray-selected AGN from the BAT AGN Spectroscopic Survey (BASS). We show that the [NII] {\lambda}6583/H{\alpha} ratio exhibits a significant correlation with ${\lambda}_{\rm Edd}$ ($R_{\rm Pear}$ = -0.44, $p$-value=$3\times10^{-13}$, {\sigma} = 0.28 dex), and the correlation is not solely driven by $M_{\rm BH}$ or $L_{\rm bol}$. The observed correlation between [NII] {\lambda}6583/H{\alpha} ratio and $M_{\rm BH}$ is stronger than the correlation with $L_{\rm bol}$, but both are weaker than the ${\lambda}_{\rm Edd}$ correlation. This implies that the large-scale narrow lines of AGN host galaxies carry information about the accretion state of the AGN central engine. We propose that the [NII] {\lambda}6583/H{\alpha} is a useful indicator of Eddington ratio with 0.6 dex of rms scatter, and that it can be used to measure ${\lambda}_{\rm Edd}$ and thus $M_{\rm BH}$ from the measured $L_{\rm bol}$, even for high redshift obscured AGN. We briefly discuss possible physical mechanisms behind this correlation, such as the mass-metallicity relation, X-ray heating, and radiatively driven outflows.
[64]  [pdf] - 1531436
The Phoenix galaxy as seen by NuSTAR
Comments: 10 pages, 11 figures, accepted for publication in A&A
Submitted: 2016-09-01
Aims. We study the long-term variability of the well-known Seyfert 2 galaxy Mrk 1210 (a.k.a. UGC 4203, or the Phoenix galaxy). Methods. The source was observed by many X-ray facilities in the last 20 years. Here we present a NuSTAR observation and put the results in context of previously published observations. Results. NuSTAR observed Mrk 1210 in 2012 for 15.4 ks. The source showed Compton-thin obscuration similar to that observed by Chandra, Suzaku, BeppoSAX and XMM-Newton over the past two decades, but different from the first observation by ASCA in 1995, in which the active nucleus was caught in a low flux state - or obscured by Compton-thick matter, with a reflection-dominated spectrum. Thanks to the high-quality hard X-ray spectrum obtained with NuSTAR and exploiting the long-term spectral coverage spanning 16.9 years, we can precisely disentangle the transmission and reflection components and put constraints on both the intrinsic long-term variability and hidden nucleus scenarios. In the former case, the distance between the reflector and the source must be at least ~ 2 pc, while in the latter one the eclipsing cloud may be identified with a water maser-emitting clump.
[65]  [pdf] - 1475452
Constraints on the Broad Line Region Properties and Extinction in Local Seyferts
Comments: 19 pages, 14 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2016-07-25
We use high spectral resolution (R > 8000) data covering 3800-13000\r{A} to study the physical conditions of the broad line region (BLR) of nine nearby Seyfert 1 galaxies. Up to six broad HI lines are present in each spectrum. A comparison - for the first time using simultaneous optical to near-infrared observations - to photoionisation calculations with our devised simple scheme yields the extinction to the BLR at the same time as determining the density and photon flux, and hence distance from the nucleus, of the emitting gas. This points to a typical density for the HI emitting gas of 10$^{11}$cm$^{-3}$ and shows that a significant amount of this gas lies at regions near the dust sublimation radius, consistent with theoretical predictions. We also confirm that in many objects the line ratios are far from case B, the best-fit intrinsic broad-line H$\alpha$/H$\beta$ ratios being in the range 2.5-6.6 as derived with our photoionization modeling scheme. The extinction to the BLR, based on independent estimates from HI and HeII lines, is A$_V$ $\le$ 3 for Seyfert 1-1.5s, while Seyfert 1.8-1.9s have A$_V$ in the range 4-8. A comparison of the extinction towards the BLR and narrow line region (NLR) indicates that the structure obscuring the BLR exists on scales smaller than the NLR. This could be the dusty torus, but dusty nuclear spirals or filaments could also be responsible. The ratios between the X-ray absorbing column N$_H$ and the extinction to the BLR are consistent with the Galactic gas-to-dust ratio if N$_H$ variations are considered.
[66]  [pdf] - 1530859
The Geometry of the Infrared and X-ray Obscurer in a Dusty Hyperluminous Quasar
Comments: Accepted to ApJ
Submitted: 2016-06-17
We study the geometry of the AGN obscurer in IRAS 09104+4109, an IR-luminous, radio-intermediate FR-I source at $z=0.442$, using infrared data from Spitzer and Herschel, X-ray data from Nustar, Swift, Suzaku, and Chandra, and an optical spectrum from Palomar. The infrared data imply a total rest-frame 1-1000$\mu$m luminosity of $5.5\times10^{46}$erg s$^{-1}$ and require both an AGN torus and starburst model. The AGN torus has an anisotropy-corrected IR luminosity of $4.9\times10^{46}$erg s$^{-1}$, and a viewing angle and half opening angle both of approximately $36$ degrees from pole-on. The starburst has a star formation rate of $(110\pm34)$M$_{\odot}$ yr$^{-1}$ and an age of $<50$Myr. These results are consistent with two epochs of luminous activity in 09104: one approximately $150$Myr ago, and one ongoing. The X-ray data suggest a photon index of $\Gamma \simeq 1.8$ and a line-of-sight column of $N_{\rm H} \simeq 5\times10^{23}$cm$^{-2}$. This argues against a reflection-dominated hard X-ray spectrum, which would have implied a much higher $N_{\rm H}$ and luminosity. The X-ray and infrared data are consistent with a bolometric AGN luminosity of $L_{\rm bol}\sim(0.5-2.5)\times10^{47}$erg s$^{-1}$. The X-ray and infrared data are further consistent with coaligned AGN obscurers in which the line of sight "skims" the torus. This is also consistent with the optical spectra, which show both coronal iron lines and broad lines in polarized but not direct light. Combining constraints from the X-ray, optical, and infrared data suggests that the AGN obscurer is within a vertical height of $20$pc, and a radius of $125$pc, of the nucleus.
[67]  [pdf] - 1443937
The nature of the torus in the heavily obscured AGN Markarian 3: an X-ray study
Comments: Accepted for publication in MNRAS, 17 pages, 11 figures, 5 tables
Submitted: 2016-05-09
In this paper we report the results of an X-ray monitoring campaign on the heavily obscured Seyfert galaxy Markarian 3 carried out between the fall of 2014 and the spring of 2015 with NuSTAR, Suzaku and XMM-Newton. The hard X-ray spectrum of Markarian 3 is variable on all the time scales probed by our campaign, down to a few days. The observed continuum variability is due to an intrinsically variable primary continuum seen in transmission through a large, but still Compton-thin column density (N_H~0.8-1.1$\times$10$^{24}$ cm$^{-2}$). If arranged in a spherical-toroidal geometry, the Compton scattering matter has an opening angle ~66 degrees and is seen at a grazing angle through its upper rim (inclination angle ~70 degrees). We report a possible occultation event during the 2014 campaign. If the torus is constituted by a system of clouds sharing the same column density, this event allows us to constrain their number (17$\pm$5) and individual column density, [~(4.9$\pm$1.5)$\times$10$^{22}$ cm$^{-2}$]. The comparison of IR and X-ray spectroscopic results with state-of-the art "torus" models suggests that at least two thirds of the X-ray obscuring gas volume might be located within the dust sublimation radius. We report also the discovery of an ionized absorber, characterised by variable resonant absorption lines due to He- and H-like iron. This discovery lends support to the idea that moderate column density absorbers could be due to clouds evaporated at the outer surface of the torus, possibly accelerated by the radiation pressure due to the central AGN emission leaking through the patchy absorber.
[68]  [pdf] - 1436358
A New Population of Compton-Thick AGN Identified Using the Spectral Curvature Above 10 keV
Comments: Accepted to ApJ, 19 pages, 12 figures
Submitted: 2016-04-26
We present a new metric that uses the spectral curvature (SC) above 10 keV to identify Compton-thick AGN in low-quality Swift BAT X-ray data. Using NuSTAR, we observe nine high SC-selected AGN. We find that high-sensitivity spectra show the majority are Compton-thick (78% or 7/9) and the remaining two are nearly Compton-thick (NH~5-8x10^23 cm^-2). We find the SC_bat and SC_nustar measurements are consistent, suggesting this technique can be applied to future telescopes. We tested the SC method on well-known Compton-thick AGN and find it is much more effective than broad band ratios (e.g. 100% using SC vs. 20% using 8-24/3-8 keV). Our results suggest that using the >10 keV emission may be the only way to identify this population since only two sources show Compton-thick levels of excess in the OIII to X-ray emission ratio (F_OIII/F_2-10 keV>1) and WISE colors do not identify most of them as AGN. Based on this small sample, we find that a higher fraction of these AGN are in the final merger stage than typical BAT AGN. Additionally, these nine obscured AGN have, on average, four times higher accretion rates than other BAT-detected AGN (Edd ratio=0.068 compared to 0.016). The robustness of SC at identifying Compton-thick AGN implies a higher fraction of nearby AGN may be Compton-thick (~22%) and the sum of black hole growth in Compton-thick AGN (Eddington ratio times population percentage), is nearly as large as mildly obscured and unobscured AGN.
[69]  [pdf] - 1400563
Determining the radio AGN contribution to the radio-FIR correlation using the black hole fundamental plane relation
Comments: (11 pages, 5 figures, accepted for publication in MNRAS)
Submitted: 2016-04-20
We investigate the 1.4 GHz radio properties of 92 nearby (z<0.05) ultra hard X-ray selected Active Galactic Nuclei (AGN) from the Swift Burst Alert Telescope (BAT) sample. Through the ultra hard X-ray selection we minimise the biases against obscured or Compton-thick AGN as well as confusion with emission derived from star formation that typically affect AGN samples selected from the UV, optical and infrared wavelengths. We find that all the objects in our sample of nearby, ultra-hard X-ray selected AGN are radio quiet; 83\% of the objects are classed as high-excitation galaxies (HEGs) and 17\% as low-excitation galaxies (LEGs). While these low-z BAT sources follow the radio--far-infrared correlation in a similar fashion to star forming galaxies, our analysis finds that there is still significant AGN contribution in the observed radio emission from these radio quiet AGN. In fact, the majority of our BAT sample occupy the same X-ray--radio fundamental plane as have been observed in other samples, which include radio loud AGN --- evidence that the observed radio emission (albeit weak) is connected to the AGN accretion mechanism, rather than star formation.
[70]  [pdf] - 1553829
Compton-thick Accretion in the local Universe
Comments: Published in ApJL
Submitted: 2016-03-15, last modified: 2016-03-22
Heavily obscured accretion is believed to represent an important stage in the growth of supermassive black holes, and to play an important role in shaping the observed spectrum of the Cosmic X-ray Background (CXB). Hard X-ray (E$>$10 keV) selected samples are less affected by absorption than samples selected at lower energies, and are therefore one of the best ways to detect and identify Compton-thick (CT, $\log N_{\rm\,H}\geq 24$) Active Galactic Nuclei (AGN). In this letter we present the first results of the largest broad-band (0.3-150 keV) X-ray spectral study of hard X-ray selected AGN to date, focusing on the properties of heavily obscured sources. Our sample includes the 834 AGN (728 non-blazar, average redshift $z\simeq 0.055$) reported in the 70-months catalog of the all-sky hard X-ray Swift/BAT survey. We find 55 CT AGN, which represent $7.6^{+1.1}_{-2.1}\%$ of our non-blazar sample. Of these, 26 are reported as candidate CT AGN for the first time. We correct for selection bias and derive the intrinsic column density distribution of AGN in the local Universe in two different luminosity ranges. We find a significant decrease in the fraction of obscured Compton-thin AGN for increasing luminosity, from $46\pm3\%$ (for $\log L_{\rm\,14-195} = 40-43.7$) to $39\pm3\%$ (for $\log L_{\rm\,14-195} = 43.7-46$). A similar trend is also found for CT AGN. The intrinsic fraction of CT AGN with $\log N_{\rm\,H}=24-25$ normalised to unity in the $\log N_{\rm H} = 20-25$ range is $27\pm4\%$, and is consistent with the observed value obtained for AGN located within 20 Mpc.
[71]  [pdf] - 1385487
Peering through the Dust: NuSTAR Observations of Two FIRST-2MASS Red Quasars
Comments: accepted for publication to ApJ; 34 pages, 8 figures, 5 tables
Submitted: 2016-02-10
Some reddened quasars appear to be transitional objects in the merger-induced black hole growth/galaxy evolution paradigm, where a heavily obscured nucleus starts to be unveiled by powerful quasar winds evacuating the surrounding cocoon of dust and gas. Hard X-ray observations are able to peer through this gas and dust, revealing the properties of circumnuclear obscuration. Here, we present NuSTAR and XMM-Newton/Chandra observations of FIRST-2MASS selected red quasars F2M 0830+3759 and F2M 1227+3214. We find that though F2M 0830+3759 is moderately obscured ($N_{\rm H,Z} = 2.1\pm0.2 \times10^{22}$ cm$^{-2}$) and F2M 1227+3214 is mildly absorbed ($N_{\rm H,Z} = 3.4^{+0.8}_{-0.7}\times10^{21}$ cm$^{-2}$) along the line-of-sight, heavier global obscuration may be present in both sources, with $N_{\rm H,S} = 3.7^{+4.1}_{-2.6} \times 10^{23}$ cm$^{-2}$ and $< 5.5\times10^{23}$ cm$^{-2}$, for F2M 0830+3759 and F2M 1227+3214, respectively. F2M 0830+3759 also has an excess of soft X-ray emission below 1 keV which is well accommodated by a model where 7% of the intrinsic AGN X-ray emission is scattered into the line-of-sight. While F2M 1227+3214 has a dust-to-gas ratio ($E(B-V)$/$N_{\rm H}$) consistent with the Galactic value, the $E(B-V)$/$N_{\rm H}$ value for F2M 0830+3759 is lower than the Galactic standard, consistent with the paradigm that the dust resides on galactic scales while the X-ray reprocessing gas originates within the dust-sublimation zone of the broad-line-region. The X-ray and 6.1$\mu$m luminosities of these red quasars are consistent with the empirical relations derived for high-luminosity, unobscured quasars, extending the parameter space of obscured AGN previously observed by NuSTAR to higher luminosities.
[72]  [pdf] - 1388926
NuSTAR observations of water megamaser AGN
Comments: 14 pages, 8 figures, accepted for publication by A&A
Submitted: 2016-02-09
Aims. Study the connection between the masing disk and obscuring torus in Seyfert 2 galaxies. Methods. We present a uniform X-ray spectral analysis of the high energy properties of 14 nearby megamaser Active Galactic Nuclei observed by NuSTAR. We use a simple analytical model to localize the maser disk and understand its connection with the torus by combining NuSTAR spectral parameters with available physical quantities from VLBI mapping. Results. Most of the sources analyzed are heavily obscured, showing a column density in excess of $\sim 10^{23}$ cm$^{-2}$. In particular, $79\%$ are Compton-thick ($N_{\rm H} > 1.5 \times 10^{24}$ cm$^{-2}$). Using column densities measured by NuSTAR, with the assumption that the torus is the extension of the maser disk, and further assuming a reasonable density profile, the torus dimensions can be predicted. They are found to be consistent with mid-IR interferometry parsec-scale observations of Circinus and NGC 1068. In this picture, the maser disk is intimately connected to the inner part of the torus. It is probably made of a large number of molecular clouds connecting the torus and the outer part of the accretion disk, giving rise to a thin disk rotating in most cases in Keplerian or sub-Keplerian motion. This toy model explains the established close connection between water megamaser emission and nuclear obscuration as a geometric effect.
[73]  [pdf] - 1378872
IC 751: a new changing-look AGN discovered by NuSTAR
Comments: Accepted for publication in ApJ, 11 pages, 6 figures
Submitted: 2016-02-01
We present the results of five NuSTAR observations of the type 2 active galactic nucleus (AGN) in IC 751, three of which were performed simultaneously with XMM-Newton or Swift/XRT. We find that the nuclear X-ray source underwent a clear transition from a Compton-thick ($N_{\rm\,H}\simeq 2\times 10^{24}\rm\,cm^{-2}$) to a Compton-thin ($N_{\rm\,H}\simeq 4\times 10^{23}\rm\,cm^{-2}$) state on timescales of $\lesssim 3$ months, which makes IC 751 the first changing-look AGN discovered by NuSTAR. Changes of the line-of-sight column density at a $\sim2\sigma$ level are also found on a time-scale of $\sim 48$ hours ($\Delta N_{\rm\,H}\sim 10^{23}\rm\,cm^{-2}$). From the lack of spectral variability on timescales of $\sim 100$ ks we infer that the varying absorber is located beyond the emission-weighted average radius of the broad-line region, and could therefore be related either to the external part of the broad-line region or a clumpy molecular torus. By adopting a physical torus X-ray spectral model, we are able to disentangle the column density of the non-varying absorber ($N_{\rm\,H}\sim 3.8\times 10^{23}\rm\,cm^{-2}$) from that of the varying clouds [$N_{\rm\,H}\sim(1-150)\times10^{22}\rm\,cm^{-2}$], and to constrain that of the material responsible for the reprocessed X-ray radiation ($N_{\rm\,H} \sim 6 \times 10^{24}\rm\,cm^{-2}$). We find evidence of significant intrinsic X-ray variability, with the flux varying by a factor of five on timescales of a few months in the 2-10 and 10-50 keV band.
[74]  [pdf] - 1351725
Extended X-ray emission in the IC 2497 - Hanny's Voorwerp system: energy injection in the gas around a fading AGN
Comments: 9 pages, 5 figures, accepted for publication in MNRAS
Submitted: 2016-01-27
We present deep Chandra X-ray observations of the core of IC 2497, the galaxy associated with Hanny's Voorwerp and hosting a fading AGN. We find extended soft X-ray emission from hot gas around the low intrinsic luminosity (unobscured) AGN ($L_{\rm bol} \sim 10^{42}-10^{44}$ erg s$^{-1}$). The temperature structure in the hot gas suggests the presence of a bubble or cavity around the fading AGN ($\mbox{E$_{\rm bub}$} \sim 10^{54} - 10^{55}$ erg). A possible scenario is that this bubble is inflated by the fading AGN, which after changing accretion state is now in a kinetic mode. Other possibilities are that the bubble has been inflated by the past luminous quasar ($L_{\rm bol} \sim 10^{46}$ erg s$^{-1}$), or that the temperature gradient is an indication of a shock front from a superwind driven by the AGN. We discuss the possible scenarios and the implications for the AGN-host galaxy interaction, as well as an analogy between AGN and X-ray binaries lifecycles. We conclude that the AGN could inject mechanical energy into the host galaxy at the end of its lifecycle, and thus provide a source for mechanical feedback, in a similar way as observed for X-ray binaries.
[75]  [pdf] - 1347639
On the relation of optical obscuration and X-ray absorption in Seyfert galaxies
Comments: 7 pages, 3 figures, accepted for publication by A&A; updated PDF to include abstract
Submitted: 2015-11-17, last modified: 2016-01-12
The optical classification of a Seyfert galaxy and whether it is considered X-ray absorbed are often used interchangeably. But there are many borderline cases and also numerous examples where the optical and X-ray classifications appear to be in conflict. In this article we re-visit the relation between optical obscuration and X-ray absorption in AGNs. We make use of our "dust color" method (Burtscher et al. 2015) to derive the optical obscuration A_V and consistently estimated X-ray absorbing columns using 0.3--150 keV spectral energy distributions. We also take into account the variable nature of the neutral gas column N_H and derive the Seyfert sub-classes of all our objects in a consistent way. We show in a sample of 25 local, hard-X-ray detected Seyfert galaxies (log L_X / (erg/s) ~ 41.5 - 43.5) that there can actually be a good agreement between optical and X-ray classification. If Seyfert types 1.8 and 1.9 are considered unobscured, the threshold between X-ray unabsorbed and absorbed should be chosen at a column N_H = 10^22.3 / cm^2 to be consistent with the optical classification. We find that N_H is related to A_V and that the N_H/A_V ratio is approximately Galactic or higher in all sources, as indicated previously. But in several objects we also see that deviations from the Galactic ratio are only due to a variable X-ray column, showing that (1) deviations from the Galactic N_H/A_V can simply be explained by dust-free neutral gas within the broad line region in some sources, that (2) the dust properties in AGNs can be similar to Galactic dust and that (3) the dust color method is a robust way to estimate the optical extinction towards the sublimation radius in all but the most obscured AGNs.
[76]  [pdf] - 1342028
Herschel Far-Infrared Photometry of the Swift Burst Alert Telescope Active Galactic Nuclei Sample of the Local Universe. II. SPIRE Observations
Comments: 22 pages, 11 figures, accepted for publication in MNRAS. We have created a Dataverse that contains all of the FITS images for the Herschel-BAT AGN as well as the photometry at . Full versions of Tables 1 and 2 and Figure 1 are available at
Submitted: 2015-12-08
We present far-infrared (FIR) and submillimeter photometry from the Herschel Space Observatory's Spectral and Photometric Imaging Receiver (SPIRE) for 313 nearby z<0.05 active galactic nuclei (AGN). We selected AGN from the 58 month Swift Burst Alert Telescope (BAT) catalog, the result of an all-sky survey in the 14-195 keV energy band, allowing for a reduction in AGN selection effects due to obscuration and host galaxy contamination. We find 46% (143/313) of our sample is detected at all three wavebands and combined with our PACS observations represents the most complete FIR spectral energy distributions of local, moderate luminosity AGN. We find no correlation between the 250, 350, and 500 micron luminosities with 14-195 keV luminosity, indicating the bulk of the FIR emission is not related to the AGN. However, Seyfert 1s do show a very weak correlation with X-ray luminosity compared to Seyfert 2s and we discuss possible explanations. We compare the SPIRE colors (F250/F350 and F350/F500) to a sample of normal star-forming galaxies, finding the two samples are statistically similar, especially after matching in stellar mass. But a color-color plot reveals a fraction of the Herschel-BAT AGN are displaced from the normal star-forming galaxies due to excess 500 micron emission E500). Our analysis shows E500 is strongly correlated with the 14-195 keV luminosity and 3.4/4.6 micron flux ratio, evidence the excess is related to the AGN. We speculate these sources are experiencing millimeter excess emission originating in the corona of the accretion disk.
[77]  [pdf] - 1312145
The NuSTAR Extragalactic Surveys: Initial Results and Catalog from the Extended Chandra Deep Field South
Comments: 20 pages, published in the Astrophysical Journal
Submitted: 2015-11-13
We present initial results and the source catalog from the NuSTAR survey of the Extended Chandra Deep Field South (hereafter, ECDFS) - currently the deepest contiguous component of the NuSTAR extragalactic survey program. The survey covers the full ~30 arcmin x 30 arcmin area of this field to a maximum depth of ~360 ks (~220 ks when corrected for vignetting at 3-24 keV), reaching sensitivity limits of ~1.3 x 10^-14 erg/cm2/s (3-8 keV), ~3.4 x 10^-14 erg/cm2/s (8-24 keV) and ~3.0 x 10^-14 erg/cm2/s (3-24 keV). Fifty four (54) sources are detected over the full field, although five of these are found to lie below our significance threshold once contaminating flux from neighboring (i.e., blended) sources is taken into account. Of the remaining 49 that are significant, 19 are detected in the 8-24 keV band. The 8-24 keV to 3-8 keV band ratios of the twelve sources that are detected in both bands span the range 0.39-1.7, corresponding to a photon index range of Gamma ~ 0.5-2.3, with a median photon index of 1.70 +/- 0.52. The redshifts of the 49 sources in our main sample span the range z = 0.21-2.7, and their rest-frame 10-40 keV luminosities (derived from the observed 8-24 keV fluxes) span the range L(10-40 keV) ~ (0.7-300) x 10^43 erg/s, sampling below the "knee" of the X-ray luminosity function out to z ~ 0.8-1. Finally, we identify one NuSTAR source that has neither a Chandra nor an XMM-Newton counterpart, but that shows evidence of nuclear activity at infrared wavelengths, and thus may represent a genuine, new X-ray source detected by NuSTAR in the ECDFS.
[78]  [pdf] - 1327579
The NuSTAR Extragalactic Surveys: First Direct Measurements of the >10 keV X-Ray Luminosity Function for Active Galactic Nuclei at z>0.1
Comments: 14 pages, accepted for publication in the Astrophysical Journal
Submitted: 2015-11-13
We present the first direct measurements of the rest-frame 10-40 keV X-ray luminosity function (XLF) of Active Galactic Nuclei (AGNs) based on a sample of 94 sources at 0.1 < z <3, selected at 8-24 keV energies from sources in the NuSTAR extragalactic survey program. Our results are consistent with the strong evolution of the AGN population seen in prior, lower-energy studies of the XLF. However, different models of the intrinsic distribution of absorption, which are used to correct for selection biases, give significantly different predictions for the total number of sources in our sample, leading to small, systematic differences in our binned estimates of the XLF. Adopting a model with a lower intrinsic fraction of Compton-thick sources and a larger population of sources with column densities N_H ~ 10^{23-24} /cm2 or a model with stronger Compton reflection component (with a relative normalization of R ~ 2 at all luminosities) can bring extrapolations of the XLF from 2-10 keV into agreement with our NuSTAR sample. Ultimately, X-ray spectral analysis of the NuSTAR sources is required to break this degeneracy between the distribution of absorbing column densities and the strength of the Compton reflection component and thus refine our measurements of the XLF. Furthermore, the models that successfully describe the high-redshift population seen by NuSTAR tend to over-predict previous, high-energy measurements of the local XLF, indicating that there is evolution of the AGN population that is not fully captured by the current models.
[79]  [pdf] - 1331118
NuSTAR catches the unveiling nucleus of NGC 1068
Comments: 5 pages, 3 figures, accepted for publication on MNRAS letters
Submitted: 2015-11-11
We present a NuSTAR and XMM-Newton monitoring campaign in 2014/2015 of the Compton-thick Seyfert 2 galaxy, NGC 1068. During the August 2014 observation, we detect with NuSTAR a flux excess above 20 keV ($32\pm6 \%$) with respect to the December 2012 observation and to a later observation performed in February 2015. We do not detect any spectral variation below 10 keV in the XMM-Newton data. The transient excess can be explained by a temporary decrease of the column density of the obscuring material along the line of sight (from N$_{\rm H}\simeq10^{25}$ cm$^{-2}$ to N$_{\rm H}=6.7\pm1.0\times10^{24}$ cm$^{-2}$), which allows us for the first time to unveil the direct nuclear radiation of the buried AGN in NGC 1068 and to infer an intrinsic 2-10 keV luminosity L$_{\rm X}=7^{+7}_{-4} \times 10^{43}$ erg s$^{-1}$.
[80]  [pdf] - 1327560
The NuSTAR view of reflection and absorption in NGC 7582
Comments: Accepted to ApJ, Nov 2015
Submitted: 2015-11-05
NGC 7582 is a well-studied X-ray bright Seyfert 2 with moderately heavy ($N_{\text{H}}\sim10^{23}-10^{24}$~cm$^{-2}$), highly variable absorption and strong reflection spectral features. The spectral shape changed around the year 2000, dropping in observed flux and becoming much more highly absorbed. Two scenarios have been put forth to explain this spectral change: 1) the central X-ray source partially ``shut off'' around this time, decreasing in intrinsic luminosity, with a delayed decrease in reflection features due to the light-crossing time of the Compton-thick material or 2) the source became more heavily obscured, with only a portion of the power law continuum leaking through. NuSTAR observed NGC~7582 twice in 2012, two weeks apart, in order to quantify the reflection using high-quality data above 10 keV. We find that the most plausible scenario is that NGC 7582 has recently become more heavily absorbed by a patchy torus with a covering fraction of $\sim\,80-90\%$ and an equatorial column density of $\sim 3 \times10^{24}$ cm$^{-2}$. We find the need for an additional highly variable full-covering absorber with $N_{\text{H}}= 4-6 \times10^{23}$ cm$^{-2}$ in the line of sight, possibly associated with a hidden broad line region.
[81]  [pdf] - 1323501
[O III]$\lambda 5007$ and X-ray Properties of a Complete Sample of Hard X-ray Selected AGNs in the Local Universe
Comments: 14 pages including 11 figures and 3 tables, accepted for publication in ApJ. In this manuscript, the observed 14-195 keV luminosities in Table 1 have been corrected to be exactly the same as in the original Swift/BAT 9-month catalog. Accordingly, Figures 2(a) and 3(a) and a part of Tables 2 and 3 have been updated. The changes from the previous version are small and do not affect the text
Submitted: 2015-10-12, last modified: 2015-11-05
We study the correlation between the [O III]$\lambda 5007$ and X-ray luminosities of local Active Galactic Nuclei (AGNs), using a complete, hard X-ray ($>10$ keV) selected sample in the Swift/BAT 9-month catalog. From our optical spectroscopic observations at the South African Astronomical Observatory and the literature, a catalog of [O III]$\lambda 5007$ line flux for all 103 AGNs at Galactic latitudes of $|b|>15^\circ$ is complied. Significant correlations with intrinsic X-ray luminosity ($L_{\rm X}$) are found both for observed ($L_{\rm [O~III]}$) and extinction-corrected ($L_{\rm [O~III]}^{\rm cor}$) luminosities, separately for X-ray unabsorbed and absorbed AGNs. We obtain the regression form of $L_{\rm [O~III]}$ $\propto L_{\rm 2-10\; keV}^{1.18\pm0.07}$ and $L_{\rm [O~III]}^{\rm cor}$ $\propto L_{\rm 2-10\; keV}^{1.16\pm0.09}$ from the whole sample. The absorbed AGNs with low ($<$0.5\%) scattering fractions in soft X-rays show on average smaller $L_{\rm [O~III]}/L_{\rm X}$ and $L_{\rm [O~III]}^{\rm cor}/L_{\rm X}$ ratios than the other absorbed AGNs, while those in edge-on host galaxies do not. These results suggest that a significant fraction of this population are buried in tori with small opening angles. By using these $L_{\rm [O~III]}$ vs. $L_{\rm X}$ correlations, the X-ray luminosity function of local AGNs (including Compton thick AGNs) in a standard population synthesis model gives much better agreement with the [O III]$\lambda 5007$ luminosity function derived from the Sloan Digital Sky Survey than previously reported. This confirms that hard X-ray observations are a very powerful tool to find AGNs with high completeness.
[82]  [pdf] - 1294142
A NuSTAR Survey of Nearby Ultraluminous Infrared Galaxies
Comments: 17 pages, 9 figures, 4 tables; Accepted for publication in the Astrophysical Journal
Submitted: 2015-10-15
We present a NuSTAR, Chandra, and XMM--Newton survey of nine of the nearest ultraluminous infrared galaxies (ULIRGs). The unprecedented sensitivity of NuSTAR at energies above 10 keV enables spectral modeling with far better precision than was previously possible. Six of the nine sources observed were detected sufficiently well by NuSTAR to model in detail their broadband X-ray spectra, and recover the levels of obscuration and intrinsic X-ray luminosities. Only one source (IRAS 13120--5453) has a spectrum consistent with a Compton--thick AGN, but we cannot rule out that a second source (Arp 220) harbors an extremely highly obscured AGN as well. Variability in column density (reduction by a factor of a few compared to older observations) is seen in IRAS 05189--2524 and Mrk 273, altering the classification of these border-line sources from Compton-thick to Compton-thin. The ULIRGs in our sample have surprisingly low observed fluxes in high energy (>10 keV) X-rays, especially compared to their bolometric luminosities. They have lower ratios of unabsorbed 2--10 keV to bolometric luminosity, and unabsorbed 2--10 keV to mid-IR [O IV] line luminosity than do Seyfert 1 galaxies. We identify IRAS 08572+3915 as another candidate intrinsically X-ray weak source, similar to Mrk 231. We speculate that the X-ray weakness of IRAS 08572+3915 is related to its powerful outflow observed at other wavelengths.
[83]  [pdf] - 1342872
The hard X-ray emission of the luminous infrared galaxy NGC 6240 as observed by NuSTAR
Comments: 14 pages, 10 figures, accepted for publication in A&A
Submitted: 2015-10-15
We present a broad-band (~0.3-70 keV) spectral and temporal analysis of NuSTAR observations of the luminous infrared galaxy NGC 6240, combined with archival Chandra, XMM-Newton and BeppoSAX data. NGC 6240 is a galaxy in a relatively early merger state with two distinct nuclei separated by ~1."5. Previous Chandra observations have resolved the two nuclei, showing that they are both active and obscured by Compton-thick material. Although they cannot be resolved by NuSTAR, thanks to the unprecedented quality of the NuSTAR data at energies >10 keV, we clearly detect, for the first time, both the primary and the reflection continuum components. The NuSTAR hard X-ray spectrum is dominated by the primary continuum piercing through an absorbing column density which is mildly optically thick to Compton scattering (tau ~ 1.2, N_H ~ 1.5 x 10^(24) cm^-2). We detect moderate hard X-ray (> 10 keV) flux variability up to 20% on short (15-20 ksec) timescales. The amplitude of the variability is maximum at ~30 keV and is likely to originate from the primary continuum of the southern nucleus. Nevertheless, the mean hard X-ray flux on longer timescales (years) is relatively constant. Moreover, the two nuclei remain Compton-thick, although we find evidence of variability of the material along the line of sight with column densities N_H <~ 2 x 10^(23) cm-2 over long (~3-15 years) timescales. The observed X-ray emission in the NuSTAR energy range is fully consistent with the sum of the best-fit models of the spatially resolved Chandra spectra of the two nuclei.
[84]  [pdf] - 1300363
The search for active black holes in nearby low-mass galaxies using optical and mid-IR data
Comments: 24 pages, 14 figures, accepted for publication in MNRAS
Submitted: 2015-09-28
We investigated AGN activity in low-mass galaxies, an important regime that can shed light onto BH formation and evolution, and their interaction with their host galaxies. We identified 336 AGN candidates from a parent sample of $\sim 48,000$ nearby low-mass galaxies ($M_{\rm \star} \leq 10^{9.5}M_\odot$, $z < 0.1$) in the SDSS. We selected the AGN using the classical BPT diagram, a similar optical emission line diagnostic based on the HeII$\lambda$4686 line, and mid-IR color cuts. Different criteria select host galaxies with different physical properties such as stellar mass and optical color, and only 3 out of 336 sources fulfill all three criteria. This could be in part due to selection biases. The resulting AGN fraction of $\sim 0.7 \%$ is at least one order of magnitude below the one estimated for more massive galaxies. At optical wavelengths, the HeII-based AGN selection appears to be more sensitive to AGN hosted in star-forming galaxies than the classical BPT diagram, at least in the low-mass regime. The archival X-ray and radio data available for some of the optically selected AGN candidates seem to confirm their AGN nature, but follow-up observations are needed to confirm the AGN nature of the rest of the sample, especially in the case of mid-IR selection. Our sample will be important for future follow-up studies aiming to understand the relation between BHs and host galaxies in the low-mass regime.
[85]  [pdf] - 1300325
BAT AGN Spectroscopic Survey II: X-ray Emission and High Ionization Optical Emission Lines
Comments: About to be accepted in MNRAS upon minor grammar changes
Submitted: 2015-09-17
We investigate the relationship between X-ray and optical line emission in 340 nearby AGN selected above 10 keV using Swift BAT. We find a weak correlation between the extinction corrected [O III] and hard X-ray luminosity (14-195 keV) with a [OIII] large scatter (R_Pear = 0.64, sigma = 0.62 dex) and a similarly large scatter with the intrinsic 2-10 keV to [O III] luminosities (RPear=0.63, sigma = 0.63 dex). Correlations of the hard X-ray fluxes with the fluxes of high-ionization narrow lines ([O III], He II, [Ne III] and [Ne V]) are not significantly better than with the low ionization lines (Halpha, [SII]). Factors like obscuration or physical slit size are not found to be a significant part of the large scatter. In contrast, the optical emission lines show much better correlations with each other (sigma = 0.3 dex) than with the X-ray flux. The inherent large scatter questions the common usage of narrow emission lines as AGN bolometric luminosity indicators and suggests that other issues such as geometrical differences in the scattering of the ionized gas or long term AGN variability are important.
[86]  [pdf] - 1251965
Decreased Specific Star Formation Rates in AGN Host Galaxies
Comments: 23 pages, 14 figures, Accepted for publication in MNRAS 2015 June 23. In original form 2015 January 22
Submitted: 2015-06-23
We investigate the location of an ultra-hard X-ray selected sample of AGN from the Swift Burst Alert Telescope (BAT) catalog with respect to the main sequence (MS) of star-forming galaxies using Herschel-based measurements of the star formation rate (SFR) and stellar mass (\mstar) from Sloan Digital Sky Survey (SDSS) photometry where the AGN contribution has been carefully removed. We construct the MS with galaxies from the Herschel Reference Survey and Herschel Stripe 82 Survey using the exact same methods to measure the SFR and \mstar{} as the Swift/BAT AGN. We find a large fraction of the Swift/BAT AGN lie below the MS indicating decreased specific SFR (sSFR) compared to non-AGN galaxies. The Swift/BAT AGN are then compared to a high-mass galaxy sample (COLD GASS), where we find a similarity between the AGN in COLD GASS and the Swift/BAT AGN. Both samples of AGN lie firmly between star-forming galaxies on the MS and quiescent galaxies far below the MS. However, we find no relationship between the X-ray luminosity and distance from the MS. While the morphological distribution of the BAT AGN is more similar to star-forming galaxies, the sSFR of each morphology is more similar to the COLD GASS AGN. The merger fraction in the BAT AGN is much higher than the COLD GASS AGN and star-forming galaxies and is related to distance from the MS. These results support a model in which bright AGN tend to be in high mass star-forming galaxies in the process of quenching which eventually starves the supermassive black hole itself.
[87]  [pdf] - 1264053
NuSTAR Reveals Extreme Absorption in z < 0.5 Type 2 Quasars
Comments: 19 pages, 13 figures, 5 tables. Accepted for publication in ApJ
Submitted: 2015-06-16
The intrinsic column density (NH) distribution of quasars is poorly known. At the high obscuration end of the quasar population and for redshifts z<1, the X-ray spectra can only be reliably characterized using broad-band measurements which extend to energies above 10 keV. Using the hard X-ray observatory NuSTAR, along with archival Chandra and XMM-Newton data, we study the broad-band X-ray spectra of nine optically selected (from the SDSS), candidate Compton-thick (NH > 1.5e24 cm^-2) type 2 quasars (CTQSO2s); five new NuSTAR observations are reported herein, and four have been previously published. The candidate CTQSO2s lie at z<0.5, have observed [OIII] luminosities in the range 8.4 < log (L_[OIII]/L_solar) < 9.6, and show evidence for extreme, Compton-thick absorption when indirect absorption diagnostics are considered. Amongst the nine candidate CTQSO2s, five are detected by NuSTAR in the high energy (8-24 keV) band: two are weakly detected at the ~ 3 sigma confidence level and three are strongly detected with sufficient counts for spectral modeling (>~ 90 net source counts at 8-24 keV). For these NuSTAR-detected sources direct (i.e., X-ray spectral) constraints on the intrinsic AGN properties are feasible, and we measure column densities ~2.5-1600 times higher and intrinsic (unabsorbed) X-ray luminosities ~10-70 times higher than pre-NuSTAR constraints from Chandra and XMM-Newton. Assuming the NuSTAR-detected type 2 quasars are representative of other Compton-thick candidates, we make a correction to the NH distribution for optically selected type 2 quasars as measured by Chandra and XMM-Newton for 39 objects. With this approach, we predict a Compton-thick fraction of f_CT = 36^{+14}_{-12} %, although higher fractions (up to 76%) are possible if indirect absorption diagnostics are assumed to be reliable.
[88]  [pdf] - 1056019
Active galactic nuclei flicker: an observational estimate of the duration of black hole growth phases of ~1e5 years
Comments: Accepted for publication in MNRAS
Submitted: 2015-05-25
We present an observational constraint for the typical active galactic nucleus (AGN) phase lifetime. The argument is based on the time lag between an AGN central engine switching on and becoming visible in X-rays, and the time the AGN then requires to photoionize a large fraction of the host galaxy. Based on the typical light travel time across massive galaxies, and the observed fraction of X-ray selected AGN without AGN-photoionized narrow lines, we estimate that the AGN phase typically lasts $\sim10^{5}$ years. This lifetime is short compared to the total growth time of $10^{7}-10^{9}$ years estimated from e.g. the Soltan argument and implies that black holes grow via many such short bursts and that AGN therefore "flicker" on and off. We discuss some consequences of this flickering behavior for AGN feedback and the analogy of X-ray binaries and AGN lifecycles.
[89]  [pdf] - 1245778
Broadband Observations of the Compton-thick Nucleus of NGC 3393
Comments: Accepted for publication in ApJ. 15 Figures and 4 tables
Submitted: 2015-05-13
We present new NuSTAR and Chandra observations of NGC 3393, a galaxy reported to host the smallest separation dual AGN resolved in the X-rays. While past results suggested a 150 pc separation dual AGN, three times deeper Chandra imaging, combined with adaptive optics and radio imaging suggest a single, heavily obscured, radio-bright AGN. Using VLA and VLBA data, we find an AGN with a two-sided jet rather than a dual AGN and that the hard X-ray, UV, optical, NIR, and radio emission are all from a single point source with a radius <0.2". We find that the previously reported dual AGN is most likely a spurious detection resulting from the low number of X-ray counts (<160) at 6-7 keV and Gaussian smoothing of the data on scales much smaller than the PSF (0.25" vs. 0.80" FWHM). We show that statistical noise in a single Chandra PSF generates spurious dual peaks of the same separation (0.55$\pm$0.07" vs. 0.6") and flux ratio (39$\pm$9% vs. 32% of counts) as the purported dual AGN. With NuSTAR, we measure a Compton-thick source (NH=$2.2\pm0.4\times10^{24}$ cm$^{-2}$) with a large torus half-opening angle, {\theta}=79 which we postulate results from feedback from strong radio jets. This AGN shows a 2-10 keV intrinsic to observed flux ratio of 150. Using simulations, we find that even the deepest Chandra observations would severely underestimate the intrinsic luminosity of NGC 3393 above z>0.2, but would detect an unobscured AGN of this luminosity out to high redshift (z=5).
[90]  [pdf] - 987922
Insights on the Dusty Torus and Neutral Torus from Optical and X-ray Obscuration in a Complete Volume Limited Hard X-ray AGN Sample
Comments: Accepted for ApJ; 15 pages with 10 figures
Submitted: 2015-05-04
We describe a complete volume limited sample of nearby active galaxies selected by their 14-195keV luminosity, and outline its rationale for studying the mechanisms regulating gas inflow and outflow. We describe also a complementary sample of inactive galaxies, selected to match the AGN host galaxy properties. The active sample appears to have no bias in terms of AGN type, the only difference being the neutral absorbing column which is two orders of magnitude greater for the Seyfert 2s. In the luminosity range spanned by the sample, log L_{14-195keV} [erg/s] = 42.4-43.7, the optically obscured and X-ray absorbed fractions are 50-65%. The similarity of these fractions to more distant spectroscopic AGN samples, although over a limited luminosity range, suggests that the torus does not strongly evolve with redshift. Our sample confirms that X-ray unabsorbed Seyfert 2s are rare, comprising not more than a few percent of the Seyfert 2 population. At higher luminosities, the optically obscured fraction decreases (as expected for the increasing dust sublimation radius), but the X-ray absorbed fraction changes little. We argue that the cold X-ray absorption in these Seyfert 1s can be accounted for by neutral gas in clouds that also contribute to the broad line region (BLR) emission; and suggest that a geometrically thick neutral gas torus co-exists with the BLR and bridges the gap to the dusty torus.
[91]  [pdf] - 1245756
A New Catalogue of Type 1 AGN and its Implication on the AGN Unified Model
Comments: Accepted for publication in ApJS. 17 pages including 17 figures and 9 tables. 2 machine-readable tables to appear online
Submitted: 2015-04-27
We have newly identified a substantial number of type 1 active galactic nuclei (AGN) featuring weak broad-line regions (BLRs) at z < 0.2 from detailed analysis of galaxy spectra in the Sloan Digital Sky Survey Data Release 7. These objects predominantly show a stellar continuum but also a broad H-alpha emission line, indicating the presence of a low-luminosity AGN oriented so that we are viewing the central engine directly without significant obscuration. These accreting black holes have previously eluded detection due to their weak nature. The new BLR AGNs we found increased the number of known type 1 AGNs by 49%. Some of these new BLR AGNs were detected at the Chandra X-ray Observatory, and their X-ray properties confirm that they are indeed type 1 AGN. Based on our new and more complete catalogue of type 1 AGNs, we derived the type 1 fraction of AGNs as a function of [OIII] 5007 emission luminosity and explored the possible dilution effect on the obscured AGN due to star-formation. The new type 1 AGN fraction shows much more complex behavior with respect to black hole mass and bolometric luminosity than suggested by the existing receding torus model. The type 1 AGN fraction is sensitive to both of these factors, and there seems to be a sweet spot (ridge) in the diagram of black hole mass and bolometric luminosity. Furthermore, we present a hint that the Eddington ratio plays a role in determining the opening angles.
[92]  [pdf] - 1055806
Determining the covering factor of Compton-thick active galactic nuclei with NuSTAR
Comments: Accepted for publication in ApJ
Submitted: 2015-02-25
The covering factor of Compton-thick obscuring material associated with the torus in active galactic nuclei (AGN) is at present best understood through the fraction of sources exhibiting Compton-thick absorption along the line of sight ($N_{H}>1.5\times10^{24}$ cm$^{-2}$) in the X-ray band, which reveals the average covering factor. Determining this Compton-thick fraction is difficult however, due to the extreme obscuration. With its spectral coverage at hard X-rays ($>$10 keV), NuSTAR is sensitive to the AGN covering factor since Compton scattering of X-rays off optically thick material dominates at these energies. We present a spectral analysis of 10 AGN observed with NuSTAR where the obscuring medium is optically thick to Compton scattering, so called Compton-thick (CT) AGN. We use the torus models of Brightman & Nandra which predict the X-ray spectrum from reprocessing in a torus and include the torus opening angle as a free parameter and aim to determine the covering factor of the Compton-thick gas in these sources individually. Across the sample we find mild to heavy Compton-thick columns, with $N_{H}$ measured from $10^{24}-10^{26}$ cm$^{-2}$, and a wide range of covering factors, where individual measurements range from 0.2-0.9. We find that the covering factor, $f_{c}$, is a strongly decreasing function of the intrinsic 2-10 keV luminosity, $L_{X}$, where $f_{c}=(-0.41\pm0.13)$log$_{10}$($L_{X}$/erg s$^{-1}$)$+18.31\pm5.33$, across more than two orders of magnitude in $L_{X}$ (10$^{41.5}-10^{44}$ erg s$^{-1}$). The covering factors measured here agree well with the obscured fraction as a function of $L_{X}$ as determined by studies of local AGN with $L_{X}>10^{42.5}$ erg s$^{-1}$.
[93]  [pdf] - 1224188
The systematic search for $z\gtrsim5$ active galactic nuclei in the $Chandra$ Deep Field South
Comments: Accepted for publication in MNRAS, 22 pages, 15 figures, for an overview of our analysis see figure 2
Submitted: 2015-01-26
We investigate early black hole (BH) growth through the methodical search for $z\gtrsim5$ AGN in the $Chandra$ Deep Field South. We base our search on the $Chandra$ 4-Ms data with flux limits of $9.1\times\ 10^{-18}$ (soft, 0.5 - 2 keV) and $5.5\times\ 10^{-17}\ \mathrm{erg}\ \mathrm{s}^{-1}\ \mathrm{cm}^{-2}$ (hard, 2 - 8 keV). At $z\sim5$ this corresponds to luminosities as low as $\sim10^{42}$ ($\sim10^{43}$) $\mathrm{erg}\ \mathrm{s}^{-1}$ in the soft (hard) band and should allow us to detect Compton-thin AGN with $M_\mathrm{BH}>10^7 M_{\odot}$ and Eddington ratios > 0.1. Our field ($0.03~\mathrm{deg}^2$) contains over 600 $z\sim5$ Lyman Break Galaxies. Based on lower redshift relations we would expect $\sim20$ of them to host AGN. After combining the $Chandra$ data with GOODS/ACS, CANDELS/WFC3 and $Spitzer$/IRAC data, the sample consists of 58 high-redshift candidates. We run a photometric redshift code, stack the GOODS/ACS data, apply colour criteria and the Lyman Break Technique and use the X-ray Hardness Ratio. We combine our tests and using additional data find that all sources are most likely at low redshift. We also find five X-ray sources without a counterpart in the optical or infrared which might be spurious detections. We conclude that our field does not contain any convincing $z\gtrsim5$ AGN. Explanations for this result include a low BH occupation fraction, a low AGN fraction, short, super-Eddington growth modes, BH growth through BH-BH mergers or in optically faint galaxies. By searching for $z\gtrsim5$ AGN we are setting the foundation for constraining early BH growth and seed formation scenarios.
[94]  [pdf] - 1202997
SDSS1133: An Unusually Persistent Transient in a Nearby Dwarf Galaxy
Comments: Accepted for publication in MNRAS
Submitted: 2014-01-27, last modified: 2014-11-18
While performing a survey to detect recoiling supermassive black holes, we have identified an unusual source having a projected offset of 800 pc from a nearby dwarf galaxy. The object, SDSS J113323.97+550415.8, exhibits broad emission lines and strong variability. While originally classified as a supernova (SN) because of its nondetection in 2005, we detect it in recent and past observations over 63 yr and find over a magnitude of rebrightening in the last 2 years. Using high-resolution adaptive optics observations, we constrain the source emission region to be <12 pc and find a disturbed host-galaxy morphology indicative of recent merger activity. Observations taken over more than a decade show narrow [O III] lines, constant ultraviolet emission, broad Balmer lines, a constant putative black hole mass over a decade of observations despite changes in the continuum, and optical emission-line diagnostics consistent with an active galactic nucleus (AGN). However, the optical spectra exhibit blueshifted absorption, and eventually narrow Fe II and [Ca II] emission, each of which is rarely found in AGN spectra. While this peculiar source displays many of the observational properties expected of a potential black hole recoil candidate, some of the properties could also be explained by a luminous blue variable star (LBV) erupting for decades since 1950, followed by a Type IIn SN in 2001. Interpreted as an LBV followed by a SN analogous to SN 2009ip, the multi-decade LBV eruptions would be the longest ever observed, and the broad Halpha emission would be the most luminous ever observed at late times (>10 yr), larger than that of unusually luminous supernovae such as SN 1988Z, suggesting one of the most extreme episodes of pre-SN mass loss ever discovered.
[95]  [pdf] - 1296083
NuSTAR Spectroscopy of Multi-Component X-ray Reflection from NGC 1068
Comments: Submitted to ApJ; 23 pages (ApJ format); 11 figures and 3 tables; Comments welcomed!
Submitted: 2014-11-03
We report on observations of NGC1068 with NuSTAR, which provide the best constraints to date on its $>10$~keV spectral shape. We find no strong variability over the past two decades, consistent with its Compton-thick AGN classification. The combined NuSTAR, Chandra, XMM-Newton, and Swift-BAT spectral dataset offers new insights into the complex reflected emission. The critical combination of the high signal-to-noise NuSTAR data and a spatial decomposition with Chandra allow us to break several model degeneracies and greatly aid physical interpretation. When modeled as a monolithic (i.e., a single N_H) reflector, none of the common Compton-reflection models are able to match the neutral fluorescence lines and broad spectral shape of the Compton reflection. A multi-component reflector with three distinct column densities (e.g., N_H~1.5e23, 5e24, and 1e25 cm^{-2}) provides a more reasonable fit to the spectral lines and Compton hump, with near-solar Fe abundances. In this model, the higher N_H components provide the bulk of the Compton hump flux while the lower N_H component produces much of the line emission, effectively decoupling two key features of Compton reflection. We note that ~30% of the neutral Fe Kalpha line flux arises from >2" (~140 pc), implying that a significant fraction of the <10 keV reflected component arises from regions well outside of a parsec-scale torus. These results likely have ramifications for the interpretation of poorer signal-to-noise observations and/or more distant objects [Abridged].
[96]  [pdf] - 1216653
An Infrared and Optical Analysis of a Sample of XBONGs and Optically Elusive AGN
Comments: 38 pages including 3 tables and 11 figures, Accepted for publication in ApJ
Submitted: 2014-09-02
We present near-infrared (NIR) spectra of four optically-elusive AGN and four X- ray bright, optically normal galaxies (XBONGs) from the Swift-BAT survey. With archival observations from SDSS, 2MASS, Spitzer and WISE, we test a number of AGN indicators in the NIR and MIR; namely, NIR emission line diagnostic ratios, the presence of coronal high-ionization lines, and infrared photometry. Of our eight hard X-ray selected AGN, we find that optical normalcy has a variety of causes from object to object, and no one explanation applies. Our objects have normal Eddington ratios, and so are unlikely to host radiatively-inefficient accretion flows (RIAFs). It is unlikely that star formation in the host or starlight dilution is contributing to their failure of optical diagnostics, except perhaps in two cases. The NIR continua are well-fit by two blackbodies: one at the stellar temperature, and a hot dust component near the dust sublimation temperature. XBONGs are more likely to have significant hot dust components, while these components are small relative to starlight in the optically-elusive AGN. Some of our sample have NIR line ratios typical of AGN, but NIR diagnostics are unsuccessful in distinguishing HII regions from AGN in general. In one object, we discover a hidden broad line region in the NIR. These results have strong relevance to the origin of optically normal AGN in deep X-ray surveys.
[97]  [pdf] - 910556
The NuSTAR View of Nearby Compton-thick AGN: The Cases of NGC 424, NGC 1320 and IC 2560
Comments: Accepted for publication in ApJ
Submitted: 2014-08-22
We present X-ray spectral analyses for three Seyfert 2 active galactic nuclei, NGC 424, NGC 1320, and IC 2560, observed by NuSTAR in the 3-79 keV band. The high quality hard X-ray spectra allow detailed modeling of the Compton reflection component for the first time in these sources. Using quasi-simultaneous NuSTAR and Swift/XRT data, as well as archival XMM-Newton data, we find that all three nuclei are obscured by Compton-thick material with column densities in excess of ~5 x $10^{24}$ cm$^{-2}$, and that their X-ray spectra above 3 keV are dominated by reflection of the intrinsic continuum on Compton-thick material. Due to the very high obscuration, absorbed intrinsic continuum components are not formally required by the data in any of the sources. We constrain the intrinsic photon indices and the column density of the reflecting medium through the shape of the reflection spectra. Using archival multi-wavelength data we recover the intrinsic X-ray luminosities consistent with the broadband spectral energy distributions. Our results are consistent with the reflecting medium being an edge-on clumpy torus with a relatively large global covering factor and overall reflection efficiency of the order of 1%. Given the unambiguous confirmation of the Compton-thick nature of the sources, we investigate whether similar sources are likely to be missed by commonly used selection criteria for Compton-thick AGN, and explore the possibility of finding their high-redshift counterparts.
[98]  [pdf] - 1216329
Weak Hard X-ray Emission from Broad Absorption Line Quasars: Evidence for Intrinsic X-ray Weakness
Comments: 11 pages, 4 figures. ApJ in press
Submitted: 2014-08-15
We report NuSTAR observations of a sample of six X-ray weak broad absorption line (BAL) quasars. These targets, at z=0.148-1.223, are among the optically brightest and most luminous BAL quasars known at z<1.3. However, their rest-frame 2 keV luminosities are 14 to >330 times weaker than expected for typical quasars. Our results from a pilot NuSTAR study of two low-redshift BAL quasars, a Chandra stacking analysis of a sample of high-redshift BAL quasars, and a NuSTAR spectral analysis of the local BAL quasar Mrk 231 have already suggested the existence of intrinsically X-ray weak BAL quasars, i.e., quasars not emitting X-rays at the level expected from their optical/UV emission. The aim of the current program is to extend the search for such extraordinary objects. Three of the six new targets are weakly detected by NuSTAR with <45 counts in the 3-24 keV band, and the other three are not detected. The hard X-ray (8-24 keV) weakness observed by NuSTAR requires Compton-thick absorption if these objects have nominal underlying X-ray emission. However, a soft stacked effective photon index ({\Gamma}~1.8) for this sample disfavors Compton-thick absorption in general. The uniform hard X-ray weakness observed by NuSTAR for this and the pilot samples selected with <10 keV weakness also suggests that the X-ray weakness is intrinsic in at least some of the targets. We conclude that the NuSTAR observations have likely discovered a significant population (>33%) of intrinsically X-ray weak objects among the BAL quasars with significantly weak <10 keV emission. We suggest that intrinsically X-ray weak quasars might be preferentially observed as BAL quasars.
[99]  [pdf] - 1215677
The variable hard X-ray emission of NGC4945 as observed by NuSTAR
Comments: Accepted for publication in ApJ
Submitted: 2014-07-15
We present a broadband (~0.5-79 keV) spectral and temporal analysis of multiple NuSTAR observations combined with archival Suzaku and Chandra data of NGC4945, the brightest extragalactic source at 100 keV. We observe hard X-ray (> 10 keV) flux and spectral variability, with flux variations of a factor 2 on timescales of 20 ksec. A variable primary continuum dominates the high energy spectrum (>10 keV) in all the states, while the reflected/scattered flux which dominates at E< 10 keV stays approximately constant. From modelling the complex reflection/transmission spectrum we derive a Compton depth along the line of sight of tau_Thomson ~ 2.9, and a global covering factor for the circumnuclear gas of ~ 0.15. This agrees with the constraints derived from the high energy variability, which implies that most of the high energy flux is transmitted, rather that Compton-scattered. This demonstrates the effectiveness of spectral analysis in constraining the geometric properties of the circumnuclear gas, and validates similar methods used for analyzing the spectra of other bright, Compton-thick AGN. The lower limits on the e-folding energy are between 200-300 keV, consistent with previous BeppoSAX, Suzaku and Swift BAT observations. The accretion rate, estimated from the X-ray luminosity and assuming a bolometric correction typical of type 2 AGN, is in the range ~0.1-0.3 lambda_Edd depending on the flux state. The substantial observed X-ray luminosity variability of NGC4945 implies that large errors can arise from using single-epoch X-ray data to derive L/L_Edd values for obscured AGNs.
[100]  [pdf] - 1215466
NuSTAR unveils a Compton-thick Type 2 quasar in Mrk 34
Comments: ApJ, in press
Submitted: 2014-07-07
We present Nustar 3-40 keV observations of the optically selected Type 2 quasar (QSO2) SDSS J1034+6001 or Mrk 34. The high-quality hard X-ray spectrum and archival XMM-Newton data can be fitted self-consistently with a reflection-dominated continuum and strong Fe Kalpha fluorescence line with equivalent-width >1 keV. Prior X-ray spectral fitting below 10 keV showed the source to be consistent with being obscured by Compton-thin column densities of gas along the line-of-sight, despite evidence for much higher columns from multiwavelength data. NuSTAR now enables a direct measurement of this column, and shows that Nh lies in the Compton-thick (CT) regime. The new data also show a high intrinsic 2-10 keV luminosity of L_{2-10}~10^{44} erg/s, in contrast to previous low-energy X-ray measurements for which L_{2-10}<~10^{43} erg/s (i.e. X-ray selection below 10 keV does not pick up this source as an intrinsically luminous obscured quasar). Both the obscuring column and the intrinsic power are about an order of magnitude (or more) larger than inferred from pre-NuSTAR X-ray spectral fitting. Mrk34 is thus a 'gold standard' CT QSO2 and is the nearest non-merging system in this class, in contrast to the other local CT quasar NGC6240 which is currently undergoing a major merger coupled with strong star-formation. For typical X-ray bolometric correction factors, the accretion luminosity of Mrk34 is high enough to potentially power the total infrared luminosity. X-ray spectral fitting also shows that thermal emission related to star-formation is unlikely to drive the observed bright soft component below ~3 keV, favoring photionization instead.
[101]  [pdf] - 1179799
The Man Behind the Curtain: X-rays Drive the UV through NIR Variability in the 2013 AGN Outburst in NGC 2617
Comments: 36 pages, 11 figures, 3 Tables. Accepted for publication in ApJ. Spectroscopic and photometric data presented in this submission are included as ancillary files. To see a video of the Swift UV/optical data see For a brief video explaining this paper, see
Submitted: 2013-10-08, last modified: 2014-06-26
After the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered a significant brightening of the inner region of NGC 2617, we began a ~70 day photometric and spectroscopic monitoring campaign from the X-ray through near-infrared (NIR) wavelengths. We report that NGC 2617 went through a dramatic outburst, during which its X-ray flux increased by over an order of magnitude followed by an increase of its optical/ultraviolet (UV) continuum flux by almost an order of magnitude. NGC 2617, classified as a Seyfert 1.8 galaxy in 2003, is now a Seyfert 1 due to the appearance of broad optical emission lines and a continuum blue bump. Such "changing look Active Galactic Nuclei (AGN)" are rare and provide us with important insights about AGN physics. Based on the Hbeta line width and the radius-luminosity relation, we estimate the mass of central black hole to be (4 +/- 1) x 10^7 M_sun. When we cross-correlate the light curves, we find that the disk emission lags the X-rays, with the lag becoming longer as we move from the UV (2-3 days) to the NIR (6-9 days). Also, the NIR is more heavily temporally smoothed than the UV. This can largely be explained by a simple model of a thermally emitting thin disk around a black hole of the estimated mass that is illuminated by the observed, variable X-ray fluxes.
[102]  [pdf] - 1215039
A comparison of the morphological properties between local and z~1 infrared luminous galaxies. Are local and high-z (U)LIRGs different?
Comments: 10 pages, 6 figures, ApJ accepted
Submitted: 2014-06-17
Ultraluminous and luminous infrared galaxies (ULIRGs and LIRGs) are the most extreme star-forming galaxies in the universe, and dominate the total star formation rate density at z>1. In the local universe (z<0.3), the majority of ULIRGs and a significant portion of LIRGs are triggered by interactions between gas-rich spiral galaxies, yet it is unclear if this is still the case at high-z. To investigate the relative importance of galaxy interactions in infrared luminous galaxies, we carry out a comparison of optical morphological properties between local (U)LIRGs and (U)LIRGs at z=0.5-1.5 based on the same sample selection, morphology classification scheme, and optical morphology at similar rest-frame wavelengths. In addition, we quantify the systematics in comparing local and high-z datasets by constructing a redshifted dataset from local (U)LIRGs, in which its data quality mimics the high-z dataset. Based on the Gini-M20 classification scheme, we find that the fraction of interacting systems decreases by ~8% from local to z<~1, and it is consistent with the reduction between local and redshifted datasets (6(+14-6)%). Based on visual classifications, the merger fraction of local ULIRGs is found to be ~20% lower compared to published results, and the reduction due to redshifiting is 15(+10-8)%. Consequently, the differences of merger fractions between local and z<~1 (U)LIRGs is only ~17%. These results demonstrate that there is no strong evolution in the fraction of (U)LIRGs classified as mergers at least out to z~1. At z>1, the morphology types of ~30% of (U)LIRGs can not be determined due to their faintness in the F814W-band, and thus the merger fraction measured at z>1 suffers from large uncertainties.
[103]  [pdf] - 1210126
The 2-79 keV X-ray Spectrum of the Circinus Galaxy with NuSTAR, XMM-Newton and Chandra: a Fully Compton-Thick AGN
Comments: Accepted for publication in ApJ
Submitted: 2014-06-12
The Circinus galaxy is one of the nearest obscured AGN, making it an ideal target for detailed study. Combining archival Chandra and XMM-Newton data with new NuSTAR observations, we model the 2-79 keV spectrum to constrain the primary AGN continuum and to derive physical parameters for the obscuring material. Chandra's high angular resolution allows a separation of nuclear and off-nuclear galactic emission. In the off-nuclear diffuse emission we find signatures of strong cold reflection, including high equivalent-width neutral Fe lines. This Compton-scattered off-nuclear emission amounts to 18% of the nuclear flux in the Fe line region, but becomes comparable to the nuclear emission above 30 keV. The new analysis no longer supports a prominent transmitted AGN component in the observed band. We find that the nuclear spectrum is consistent with Compton-scattering by an optically-thick torus, where the intrinsic spectrum is a powerlaw of photon index $\Gamma = 2.2-2.4$, the torus has an equatorial column density of $N_{\rm H} = (6-10)\times10^{24}$cm$^{-2}$ and the intrinsic AGN $2-10$ keV luminosity is $(2.3-5.1)\times 10^{42}$ erg/s. These values place Circinus along the same relations as unobscured AGN in accretion rate-vs-$\Gamma$ and $L_X$-vs-$L_{IR}$ phase space. NuSTAR's high sensitivity and low background allow us to study the short time-scale variability of Circinus at X-ray energies above 10 keV for the first time. The lack of detected variability favors a Compton-thick absorber, in line with the the spectral fitting results.
[104]  [pdf] - 1208347
NuSTAR and XMM-Newton Observations of Luminous, Heavily Obscured, WISE-Selected Quasars at z ~ 2
Comments: 7 pages; submitted to ApJ
Submitted: 2014-03-12
We report on a NuSTAR and XMM-Newton program that has observed a sample of three extremely luminous, heavily obscured WISE-selected AGN at z~2 in a broad X-ray band (0.1 - 79 keV). The parent sample, selected to be faint or undetected in the WISE 3.4um (W1) and 4.6um (W2) bands but bright at 12um (W3) and 22um (W4), are extremely rare, with only ~1000 so-called W1W2-dropouts across the extragalactic sky. Optical spectroscopy reveals typical redshifts of z~2 for this population, implying rest-frame mid-IR luminosities of L(6um)~6e46 erg/s and bolometric luminosities that can exceed L(bol)~1e14 L(sun). The corresponding intrinsic, unobscured hard X-ray luminosities are L(2-10)~4e45 erg/s for typical quasar templates. These are amongst the most luminous AGN known, though the optical spectra rarely show evidence of a broad-line region and the selection criteria imply heavy obscuration even at rest-frame 1.5um. We designed our X-ray observations to obtain robust detections for gas column densities N(H)<1e24 /cm2. In fact, the sources prove to be fainter than these predictions. Two of the sources were observed by both NuSTAR and XMM-Newton, with neither being detected by NuSTAR and one being faintly detected by XMM-Newton. A third source was observed only with XMM-Newton, yielding a faint detection. The X-ray data require gas column densities N(H)>1e24 /cm2, implying the sources are extremely obscured, consistent with Compton-thick, luminous quasars. The discovery of a significant population of heavily obscured, extremely luminous AGN does not conform to the standard paradigm of a receding torus, in which more luminous quasars are less likely to be obscured. If a larger sample conforms with this finding, then this suggests an additional source of obscuration for these extreme sources.
[105]  [pdf] - 1203305
NuSTAR Observations of Heavily Obscured Quasars at z ~ 0.5
Comments: 10 pages, 6 figures, 3 tables. Accepted for publication in ApJ
Submitted: 2014-02-11
We present NuSTAR hard X-ray (3-79 keV) observations of three Type 2 quasars at z ~ 0.4-0.5, optically selected from the Sloan Digital Sky Survey (SDSS). Although the quasars show evidence for being heavily obscured Compton-thick systems on the basis of the 2-10 keV to [OIII] luminosity ratio and multiwavelength diagnostics, their X-ray absorbing column densities (N_H) are poorly known. In this analysis: (1) we study X-ray emission at >10 keV, where X-rays from the central black hole are relatively unabsorbed, in order to better constrain N_H; (2) we further characterize the physical properties of the sources through broad-band near-UV to mid-IR spectral energy distribution (SED) analyses. One of the quasars is detected with NuSTAR at >8 keV with a no-source probability of <0.1%, and its X-ray band ratio suggests near Compton-thick absorption with N_H \gtrsim 5 x 10^23 cm^-2. The other two quasars are undetected, and have low X-ray to mid-IR luminosity ratios in both the low energy (2-10 keV) and high energy (10-40 keV) X-ray regimes that are consistent with extreme, Compton-thick absorption (N_H \gtrsim 10^24 cm^-2). We find that for quasars at z ~ 0.5, NuSTAR provides a significant improvement compared to lower energy (<10 keV) Chandra and XMM-Newton observations alone, as higher column densities can now be directly constrained.
[106]  [pdf] - 1322324
Do Most Active Galactic Nuclei Live in High Star Formation Nuclear Cusps?
Comments: 5 pages, 5 figures, accepted for publication in ApJ Letters
Submitted: 2013-12-30
We present early results of the Herschel PACS (70 and 160 \micron{}) and SPIRE (250, 350, and 500 \micron{}) survey of 313 low redshift ($\rm{z} < 0.05$), ultra-hard X-ray (14--195 keV) selected AGN from the 58 month Swift/BAT catalog. Selection of AGN from ultra-hard X-rays avoids bias from obscuration providing a complete sample of AGN to study the connection between nuclear activity and star formation in host galaxies. With the high angular resolution of PACS, we find that $>$35\%$ and $>$20\%$ of the sources are "point-like" at 70 and 160 \micron{} respectively and many more that have their flux dominated by a point source located at the nucleus. The inferred star formation rates (SFR) of 0.1 - 100 M$_{\sun}$ yr$^{-1}$ using the 70 and 160 \micron{} flux densities as SFR indicators are consistent with those inferred from Spitzer NeII fluxes, but we find that 11.25 \micron{} PAH data give $\sim$3x lower SFR. Using GALFIT to measure the size of the FIR emitting regions, we determined the SFR surface density [M$_{\sun}$ yr$^{-1}$ kpc$^{-2}$] for our sample, finding a significant fraction of these sources exceed the threshold for star formation driven winds (0.1 M$_{\sun}$ yr$^{-1}$ kpc$^{-2}$).
[107]  [pdf] - 1179314
The role of galaxy interaction in the SFR-M relation: characterizing morphological properties of Herschel-selected galaxies at 0.2<z<1.5
Comments: 14 pages, 11 figures, ApJ accepted
Submitted: 2013-09-17
Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z<0.25). Recent studies based on the relation between galaxies' star formation rates and stellar mass (the SFR-M relation or the galaxy main sequence (MS)) have suggested that galaxy interaction/mergers may only become significant when galaxies fall well above the galaxy MS. Since the typical SFR at given M increases with redshift, the existence of galaxy MS implies that massive, IR-luminous galaxies at high-z may not necessarily be driven by galaxy interactions. We examine the role of galaxy interactions in the SFR-M relation by carrying out a morphological analysis of 2084 Herschel-selected galaxies at 0.2 < z < 1.5 in the COSMOS field. Herschel-PACS and -SPIRE observations covering the full 2-deg^2 COSMOS field provide one of the largest far-IR selected samples of high-redshift galaxies with well-determined redshifts to date, with sufficient sensitivity at z ~ 1, to sample objects lying on and above the galaxy MS. Using a detailed visual classification scheme, we show that the fraction of "disk galaxies" decreases and the fraction of "irregular" galaxies increases systematically with increasing LIR out to z ~ 1.5 and z ~ 1.0, respectively. At log(LIR) > 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model.
[108]  [pdf] - 1164540
Characterisation of SCUBA-2 450um and 850um-selected Galaxies in the COSMOS Field
Comments: 37 pages (21 of text),14 figures, updated to MNRAS accepted version
Submitted: 2013-02-11, last modified: 2013-09-04
We present deep 450um and 850um observations of a large, uniformly covered 394arcmin^2 area in the COSMOS field obtained with the SCUBA-2 instrument on the James Clerk Maxwell Telescope (JCMT). We achieve root-mean-square noise values of 4.13mJy at 450um and 0.80mJy at 850um. The differential and cumulative number counts are presented and compared to similar previous works. Individual point sources are identified at >3.6sigma significance, a threshold corresponding to a 3-5% sample contamination rate. We identify 78 sources at 450um and 99 at 850um, with flux densities S450=13-37mJy and S850=2-16mJy. Only 62-76% of 450um sources are 850um detected and 61-81% of 850um sources are 450um detected. The positional uncertainties at 450um are small (1-2.5") and therefore allow a precise identification of multiwavelength counterparts without reliance on detection at 24um or radio wavelengths; we find that only 44% of 450um-selected galaxies and 60% of 850um-sources have 24um or radio counterparts. 450um-selected galaxies peak at <z>=1.95+-0.19 and 850um=selected galaxies peak at <z>=2.16+-0.11. The two samples occupy similar parameter space in redshift and luminosity, while their median SED peak wavelengths differ by ~10-50um (translating to deltaTdust =8-12K, where 450um-selected galaxies are warmer). The similarities of the 450um and 850um populations, yet lack of direct overlap between them, suggests that submillimeter surveys conducted at any single far-infrared wavelength will be significantly incomplete (~>30%) at censusing infrared-luminous star formation at high-z.
[109]  [pdf] - 1173015
A SCUBA-2 850 micron survey of protoplanetary discs in the sigma Orionis cluster
Comments: 10 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2013-07-26
We present the results from a large 850 micron survey of the sigma Orionis cluster using the SCUBA-2 camera on the James Clerk Maxwell Telescope. The 0.5-degree diameter circular region we surveyed contains 297 young stellar objects with an age estimated at about 3Myr. We detect 9 of these objects, 8 of which have infrared excesses from an inner disc. We also serendipitously detect 3 non-stellar sources at > 5sigma that are likely background submillimetre galaxies. The 9 detected stars have inferred disc masses ranging from 5 to about 17MJup, assuming similar dust properties as Taurus discs and an ISM gas-to-dust ratio of 100. There is a net positive signal toward the positions of the individually undetected infrared excess sources indicating a mean disc mass of 0.5 MJup . Stacking the emission toward those stars without infrared excesses constrains their mean disc mass to less than 0.3MJup, or an equivalent Earth mass in dust. The submillimetre luminosity distribution is significantly different from that in the younger Taurus region, indicating disc mass evolution as star forming regions age and the infrared excess fraction decreases. Submillimeter Array observations reveal CO emission toward 4 sources demonstrating that some, but probably not much, molecular gas remains in these relatively evolved discs. These observations provide new constraints on the dust and gas mass of protoplanetary discs during the giant planet building phase and provide a reference level for future studies of disc evolution.
[110]  [pdf] - 1166671
X-ray view of four high-luminosity Swift/BAT AGN: Unveiling obscuration and reflection with Suzaku
Comments: A&A paper in press, 17 pages
Submitted: 2013-05-15
The Swift/BAT nine-month survey observed 153 AGN, all with ultra-hard X-ray BAT fluxes in excess of 10^-11 erg cm^-2 s^-1 and an average redshift of 0.03. Among them, four of the most luminous BAT AGN (44.73 < Log L(BAT) < 45.31) were selected as targets of Suzaku follow-up observations: J2246.0+3941 (3C 452), J0407.4+0339 (3C 105), J0318.7+6828, and J0918.5+0425. The column density, scattered/reflected emission, the properties of the Fe K line, and a possible variability are fully analyzed. For the latter, the spectral properties from Chandra, XMM-Newton and Swift/XRT public observations were compared with the present Suzaku analysis. Of our sample, 3C 452 is the only certain Compton-thick AGN candidate because of i) the high absorption and strong Compton reflection; ii) the lack of variability; iii) the "buried" nature, i.e. the low scattering fraction (<0.5%) and the extremely low relative [OIII] luminosity. In contrast 3C 105 is not reflection-dominated, despite the comparable column density, X-ray luminosity and radio morphology, but shows a strong long-term variability in flux and scattering fraction, consistent with the soft emission being scattered from a distant region (e.g., the narrow emission line region). The sample presents high (>100) X-to-[OIII] luminosity ratios, confirming the [OIII] luminosity to be affected by residual extinction in presence of mild absorption, especially for "buried" AGN such as 3C 452. Three of our targets are powerful FRII radio galaxies, making them the most luminous and absorbed AGN of the BAT Seyfert survey despite the inversely proportional N_H - L_X relation.
[111]  [pdf] - 1159492
Studying Faint Ultra Hard X-ray Emission from AGN in GOALS LIRGs with Swift BAT
Comments: Accepted for publication in ApJL
Submitted: 2013-02-04
We present the first analysis of the all-sky Swift BAT ultra hard X-ray (14-195 keV) data for a targeted list of objects. We find the BAT data can be studied at 3x fainter limits than in previous blind detection catalogs based on prior knowledge of source positions and using smaller energy ranges for source detection. We determine the AGN fraction in 134 nearby (z<0.05) luminous infrared galaxies (LIRGS) from the GOALS sample. We find that LIRGs have a higher detection frequency than galaxies matched in stellar mass and redshift at 14-195 keV and 24-35 keV. In agreement with work at other wavelengths, the AGN detection fraction increases strongly at high IR luminosity with half of high luminosity LIRGs (50%, 6/12, log L_IR/L_sun>11.8) detected. The BAT AGN classification shows 97% (37/38) agreement with Chandra and XMM AGN classification using hardness ratios or detection of a iron K-alpha line. This confirms our statistical analysis and supports the use of the Swift BAT all-sky survey to study fainter populations of any category of sources in the ultra hard X-ray band. BAT AGN in LIRGs tend to show higher column densities with 40\pm9% showing 14-195 keV/2-10 keV hardness flux ratios suggestive of high or Compton-thick column densities (log N_H>24 cm^-2), compared to only 12\pm5% of non-LIRG BAT AGN. We also find that using specific energy ranges of the BAT detector can yield additional sources over total band detections with 24% (5/21) of detections in LIRGs at 24-35 keV not detected at 14-195 keV.
[112]  [pdf] - 1092962
Understanding Dual AGN Activation in the Nearby Universe
Comments: Accepted for publication in ApJL
Submitted: 2012-01-13
We study the fraction of dual AGN in a sample of 167 nearby (z<0.05), moderate luminosity, ultra hard X-ray selected AGN from the all-sky Swift BAT survey. Combining new Chandra and Gemini observations together with optical and X-ray observations, we find that the dual AGN frequency at scales <100 kpc is 10% (16/167). Of the 16 dual AGN, 3 (19%) were detected using X-ray spectroscopy and were not detected using emission line diagnostics. Close dual AGN (<30 kpc) tend to be more common among the most X-ray luminous systems. In dual AGN, the X-ray luminosity of both AGN increases strongly with decreasing galaxy separation, suggesting that the merging event is key in powering both AGN. 50% of the AGN with a very close companion (<15 kpc), are dual AGN. We also find that dual AGN are more likely to occur in major mergers and tend to avoid absorption line galaxies with elliptical morphologies. Finally, we find SDSS Seyferts are much less likely than BAT AGN (0.25% vs. 7.8%) to be found in dual AGN at scales <30 kpc because of a smaller number of companions galaxies, fiber collision limits, a tendency for AGN at small separations to be detected only in X-rays, and a higher fraction of dual AGN companions with increasing AGN luminosity.
[113]  [pdf] - 413602
Host Galaxy Properties of the Swift BAT Ultra Hard X-ray Selected AGN
Comments: 30 pages, 30 figures, accepted for publication in ApJ
Submitted: 2011-07-06
We have assembled the largest sample of ultra hard X-ray selected (14-195 keV) AGN with host galaxy optical data to date, with 185 nearby (z<0.05), moderate luminosity AGN from the Swift BAT sample. The BAT AGN host galaxies have intermediate optical colors (u-r and g-r) that are bluer than a comparison sample of inactive galaxies and optically selected AGN from the Sloan Digital Sky Survey (SDSS) which are chosen to have the same stellar mass. Based on morphological classifications from the RC3 and the Galaxy Zoo, the bluer colors of BAT AGN are mainly due to a higher fraction of mergers and massive spirals than in the comparison samples. BAT AGN in massive galaxies (log M_*>10.5) have a 5 to 10 times higher rate of spiral morphologies than in SDSS AGN or inactive galaxies. We also see enhanced far-IR emission in BAT AGN suggestive of higher levels of star formation compared to the comparison samples. BAT AGN are preferentially found in the most massive host galaxies with high concentration indexes indicative of large bulge-to-disk ratios and large supermassive black holes. The narrow-line BAT AGN have similar intrinsic luminosities as the SDSS NL Seyferts based on measurements of [O III]. There is also a correlation between the stellar mass and X-ray emission. The BAT AGN in mergers have bluer colors and greater ultra hard X-ray emission compared to the BAT sample as whole. In agreement with the Unified Model of AGN, and the relatively unbiased nature of the BAT sources, the host galaxy colors and morphologies are independent of measures of obscuration such as X-ray column density or Seyfert type. The high fraction of massive spiral galaxies and galaxy mergers in BAT AGN suggest that host galaxy morphology is related to the activation and fueling of local AGN.
[114]  [pdf] - 382619
Chandra Discovery of a Binary Active Galactic Nucleus in Mrk 739
Comments: 7 pages, 5 figures, Published in Astrophysical Journal Letters. Reformatted and fixed typos in v1
Submitted: 2011-06-10, last modified: 2011-07-06
We have discovered a binary AGN in the galaxy Mrk 739 using Chandra and Swift BAT. We find two luminous L_2-10 keV=1.1*10^43 and 1.0*10^42 erg/s, unresolved nuclei with a projected separation of 3.4 kpc (5.8\pm0.1) coincident with two bulge components in the optical image. The western X-ray source (Mrk 739W) is highly variable (2.5x) during the 4-hour Chandra observation and has a very hard spectrum consistent with an AGN. While the eastern component was already known to be an AGN based on the presence of broad optical recombination lines, Mrk 739W shows no evidence of being an AGN in optical, UV, and radio observations, suggesting the critical importance of high spatial resolution hard X-ray observations (>2 keV) in finding these binary AGN. A high level of star formation combined with a very low L_[O III]/L_2-10 keV ratio cause the AGN to be missed in optical observations. CO observations of the (3-2) and (2-1) lines indicate large amounts of molecular gas in the system that could be driven towards the black holes during the violent galaxy collision and be key to fueling the binary AGN. Mrk 739E has a high Eddington ratio of 0.71 and a small black hole (log M_BH=7.05\pm0.3) consistent with an efficiently accreting AGN. Other than NGC 6240, this stands as the nearest case of a binary AGN discovered to date.
[115]  [pdf] - 1076411
XMM Follow-Up Observations of Three Swift BAT-Selected Active Galactic Nuclei
Comments: 18 pages including 4 figures, accepted by ApJ
Submitted: 2011-05-04
We present XMM-Newton observations of three AGN taken as part of a hunt to find very heavily obscured Compton-thick AGN. For obscuring columns greater than 10^25 cm^-2, AGN are only visible at energies below 10 keV via reflected/scattered radiation, characterized by a flat power-law. We therefore selected three objects (ESO 417-G006, IRAS 05218-1212, and MCG -01-05-047) from the Swift BAT hard X-ray survey catalog with Swift X-ray Telescope XRT 0.5-10 keV spectra with flat power-law indices as candidate Compton-thick sources for follow-up observations with the more sensitive instruments on XMM-Newton. The XMM spectra, however, rule out reflection-dominated models based on the weakness of the observed Fe K-alpha lines. Instead, the spectra are well-fit by a model of a power-law continuum obscured by a Compton-thin absorber, plus a soft excess. This result is consistent with previous follow-up observations of two other flat-spectrum BAT-detected AGN. Thus, out of the six AGN in the 22-month BAT catalog with apparently flat Swift XRT spectra, all five that have had follow-up observations are not likely Compton-thick. We also present new optical spectra of two of these objects, IRAS 05218-1212 and MCG -01-05-047. Interestingly, though both these AGN have similar X-ray spectra, their optical spectra are completely different, adding evidence against the simplest form of the geometric unified model of AGN. IRAS 05218-1212 appears in the optical as a Seyfert 1, despite the ~8.5x10^22 cm^-2 line-of-sight absorbing column indicated by its X-ray spectrum. MCG -01-05-047's optical spectrum shows no sign of AGN activity; it appears as a normal galaxy.
[116]  [pdf] - 176517
Merging and Clustering of the Swift BAT AGN Sample
Comments: 15 pages, 4 figures, published in the Astrophysical Journal Letters
Submitted: 2010-06-01
We discuss the merger rate, close galaxy environment, and clustering on scales up to a Mpc of the SWIFT BAT hard X-ray sample of nearby (z<0.05), moderate-luminosity active galactic nuclei (AGN). We find a higher incidence of galaxies with signs of disruption compared to a matched control sample (18% versus 1%) and of close pairs within 30 kpc (24% versus 1%). We also find a larger fraction with companions compared to normal galaxies and optical emission line selected AGN at scales up to 250 kpc. We hypothesize that these merging AGN may not be identified using optical emission line diagnostics because of optical extinction and dilution by star formation. In support of this hypothesis, in merging systems we find a higher hard X-ray to [OIII] flux ratio, as well as emission line diagnostics characteristic of composite or star-forming galaxies, and a larger IRAS 60 um to stellar mass ratio.
[117]  [pdf] - 1024995
Scattered X-rays in Obscured Active Galactic Nuclei and their Implications for Geometrical Structure and Evolution
Comments: 13 pages, 10 figures, accepted for publication in ApJ
Submitted: 2010-02-04
We construct a new sample of 32 obscured active galactic nuclei (AGNs) selected from the Second XMM-Newton Serendipitous Source Catalogue to investigate their multiwavelength properties in relation to the "scattering fraction", the ratio of the soft X-ray flux to the absorption-corrected direct emission. The sample covers a broad range of the scattering fraction (0.1%-10%). A quarter of the 32 AGNs have a very low scattering fraction (smaller than 0.5%), which suggests that they are buried in a geometrically thick torus with a very small opening angle. We investigate correlations between the scattering fraction and multiwavelength properties. We find that AGNs with a small scattering fraction tend to have low [O III]lambda5007/X-ray luminosity ratios. This result agrees with the expectation that the extent of the narrow-line region is small because of the small opening angle of the torus. There is no significant correlation between scattering fraction and far-infrared luminosity. This implies that a scale height of the torus is not primarily determined by starburst activity. We also compare scattering fraction with black hole mass or Eddington ratio and find a weak anti-correlation between the Eddington ratio and scattering fraction. This implies that more rapidly growing supermassive black holes tend to have thicker tori.
[118]  [pdf] - 1018940
Optical Spectral Properties of Swift BAT Hard X-ray Selected Active Galactic Nuclei Sources
Comments: 84 pages, 20 figures, 17 tables, accepted in ApJ
Submitted: 2009-12-30
The Swift Burst Alert Telescope (BAT) survey of Active Galactic Nuclei (AGN) is providing an unprecedented view of local AGNs (<z> = 0.03) and their host galaxy properties. In this paper, we present an analysis of the optical spectra of a sample of 64 AGNs from the 9-month survey, detected solely based on their 14-195 keV flux. Our analysis includes both archived spectra from the Sloan Digital Sky Survey and our own observations from the 2.1-m Kitt Peak National Observatory telescope. Among our results, we include line ratio classifications utilizing standard emission line diagnostic plots, [O III] 5007 A luminosities, and H-beta derived black hole masses. As in our X-ray study, we find the type 2 sources to be less luminous (in [O III] 5007 A and 14-195 keV luminosities) with lower accretion rates than the type 1 sources. We find that the optically classified LINERs, H II/composite galaxies, and ambiguous sources have the lowest luminosities, while both broad line and narrow line Seyferts have similar luminosities. From a comparison of the hard X-ray (14-195 keV) and [O III] luminosities, we find that both the observed and extinction-corrected [O III] luminosities are weakly correlated with X-ray luminosity. In a study of the host galaxy properties from both continuum fits and measurements of the stellar absorption indices, we find that the hosts of the narrow line sources have properties consistent with late type galaxies.
[119]  [pdf] - 1018633
Optical and X-ray Properties of the Swift BAT-detected AGN
Comments: 4 pages, 2 figures, to appear in proceedings for 'X-ray Astronomy 2009', Bologna 09/2009, AIP Conference Series, Eds. A. Comastri, M. Cappi, L. Angelini
Submitted: 2009-12-04
The Swift Gamma-Ray Burst satellite has detected a largely unbiased towards absorption sample of local ($<z> \approx 0.03$) AGN, based solely on their 14--195 keV flux. In the first 9 months of the survey, 153 AGN sources were detected. The X-ray properties in the 0.3--10 keV band have been compiled and presented based on analyses with XMM-Newton, Chandra, Suzaku, and the Swift XRT (Winter et al. 2009). Additionally, we have compiled a sub-sample of sources with medium resolution optical ground-based spectra from the SDSS or our own observations at KPNO. In this sample of 60 sources, we have classified the sources using standard emission line diagnostic plots, obtained masses for the broad line sources through measurement of the broad H$\beta$ emission line, and measured the [OIII] 5007\AA luminosity of this sample. Based on continuum fits to the intrinsic absorption features, we have obtained clues about the stellar populations of the host galaxies. We now present the highlights of our X-ray and optical studies of this unique sample of local AGNs, including a comparison of the 2--10 keV and 14--195 keV X-ray luminosities with the [OIII] 5007\AA luminosity and the implications of our results towards measurements of bolometric luminosities.
[120]  [pdf] - 133450
The 22-Month Swift-BAT All-Sky Hard X-ray Survey
Comments: Accepted by the Astrophysical Journal Supplement Series. 27 pages. This version was accepted by the journal and includes changes to the text and figures in response to the referee's comments. The main data table remains substantially the same as the previous version
Submitted: 2009-03-18, last modified: 2009-09-17
We present the catalog of sources detected in the first 22 months of data from the hard X-ray survey (14--195 keV) conducted with the BAT coded mask imager on the \swift satellite. The catalog contains 461 sources detected above the 4.8 sigma level with BAT. High angular resolution X-ray data for every source from Swift XRT or archival data have allowed associations to be made with known counterparts in other wavelength bands for over 97% of the detections, including the discovery of ~30 galaxies previously unknown as AGN and several new Galactic sources. A total of 266 of the sources are associated with Seyfert galaxies (median redshift z ~ 0.03) or blazars, with the majority of the remaining sources associated with X-ray binaries in our Galaxy. This ongoing survey is the first uniform all sky hard X-ray survey since HEAO-1 in 1977. Since the publication of the 9-month BAT survey we have increased the number of energy channels from 4 to 8 and have substantially increased the number of sources with accurate average spectra. The BAT 22-month catalog is the product of the most sensitive all-sky survey in the hard X-ray band, with a detection sensitivity (4.8 sigma) of 2.2e-11 erg/cm2/s (1 mCrab) over most of the sky in the 14--195 keV band.
[121]  [pdf] - 8338
Swift and Chandra Detections of Supernova 2006jc: Evidence for Interaction of the Supernova Shock with a Circumstellar Shell
Comments: 4 pages, 4 figures, submitted to ApJL (v2)
Submitted: 2007-12-19, last modified: 2008-01-08
The peculiar Type Ib supernova (SN) 2006jc has been observed with the UV/Optical Telescope (UVOT) and X-Ray Telescope (XRT) on board the Swift observatory over a period of 19 to 183 days after the explosion. Signatures of interaction of the outgoing SN shock with dense circumstellar material (CSM) are detected, such as strong X-ray emission (L_{0.2-10} > E39 erg/s) and the presence of MgII 2800A line emission visible in the UV spectra. In combination with a Chandra observation obtained on day 40 after the explosion, the X-ray light curve is constructed, which shows a unique rise of the X-ray emission by a factor of ~5 over a period of ~4 months, followed by a rapid decline. We interpret the unique X-ray and UV properties as a result of the SN shock interacting with a shell of material that was deposited by an outburst of the SN progenitor two years prior to the explosion. Our results are consistent with the explosion of a Wolf-Rayet star that underwent an episodic mass ejection qualitatively similar to those of luminous blue variable stars prior to its explosion. This led to the formation of a dense (>E7 cm**-3) shell at a distance of ~E16 cm from the site of the explosion, which expands with the WR wind at a velocity of (1300+/-300) km/s.
[122]  [pdf] - 133786
X-Ray, UV, and Optical Observations of Supernova 2006bp with Swift: Detection of Early X-Ray Emission
Comments: 8 pages, 5 figures, accepted for publication in ApJ
Submitted: 2007-03-29
We present results on the X-ray and optical/UV emission from the type IIP SN 2006bp and the interaction of the SN shock with its environment, obtained with the X-Ray Telescope (XRT) and UV/Optical Telescope (UVOT) on-board the Swift observatory. SN 2006bp is detected in X-rays at a 4.5 sigma level of significance in the merged XRT data from days 1 to 12 after the explosion. If the X-ray luminosity of (1.8+/-0.4)E39 ergs/s is caused by interaction of the SN shock with circumstellar material (CSM), deposited by a stellar wind from the progenitor's companion star, a mass-loss rate of ~E-05 M_sun/yr is inferred. The mass-loss rate is consistent with the non-detection in the radio with the VLA on days 2, 9, and 11 after the explosion and characteristic of a red supergiant progenitor with a mass around 12-15 M_sun prior to the explosion. In combination with a follow-up XMM-Newton observation obtained on day 21 after the explosion, an X-ray rate of decline with index 1.2+/-0.6 is inferred. Since no other SN has been detected in X-rays prior to the optical peak and since type IIP SNe have an extended 'plateau' phase in the optical, we discuss the scenario that the X-rays might be due to inverse Compton scattering of photospheric optical photons off relativistic electrons produced in circumstellar shocks. However, due to the high required value of the Lorentz factor (~10-100) we conclude that Inverse Compton scattering is an unlikely explanation for the observed X-ray emission. The fast evolution of the optical/ultraviolet spectral energy distribution and the spectral changes observed with Swift reveal the onset of metal line-blanketing and cooling of the expanding photosphere during the first few weeks after the outburst.