Search

The first very long baseline interferometry (VLBI) detections at 870$\mu$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$\lambda$ corresponding to an angular resolution, or fringe spacing, of 19$\mu$as. The Allan deviation of the visibility phase at 870$\mu$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$\mu$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time., Comment: Corresponding author: S. Doeleman

Hard X-ray-selected samples of Active Galactic Nuclei (AGN) provide one of the cleanest views of supermassive black hole accretion, but are biased against objects obscured by Compton-thick gas column densities of $N_{\rm H}$ $>$ 10$^{24}$ cm$^{-2}$. To tackle this issue, we present the NuSTAR Local AGN $N_{\rm H}$ Distribution Survey (NuLANDS)$-$a legacy sample of 122 nearby ($z$ $<$ 0.044) AGN primarily selected to have warm infrared colors from IRAS between 25$-$60 $\mu$m. We show that optically classified type 1 and 2 AGN in NuLANDS are indistinguishable in terms of optical [OIII] line flux and mid-to-far infrared AGN continuum bolometric indicators, as expected from an isotropically selected AGN sample, while type 2 AGN are deficient in terms of their observed hard X-ray flux. By testing many X-ray spectroscopic models, we show the measured line-of-sight column density varies on average by $\sim$ 1.4 orders of magnitude depending on the obscurer geometry. To circumvent such issues we propagate the uncertainties per source into the parent column density distribution, finding a directly measured Compton-thick fraction of 35 $\pm$ 9%. By construction, our sample will miss sources affected by severe narrow-line reddening, and thus segregates sources dominated by small-scale nuclear obscuration from large-scale host-galaxy obscuration. This bias implies an even higher intrinsic obscured AGN fraction may be possible, although tests for additional biases arising from our infrared selection find no strong effects on the measured column-density distribution. NuLANDS thus holds potential as an optimized sample for future follow-up with current and next-generation instruments aiming to study the local AGN population in an isotropic manner., Comment: Accepted for publication in ApJ. 50 pages (78 including appendix and bibliography), 21 figures

Since the discovery of the cosmic X-ray background (CXB), astronomers have strived to understand the accreting supermassive black holes (SMBHs) contributing to its peak in the 10-40 keV band. Existing soft X-ray telescopes could study this population up to only 10 keV, and, while NuSTAR (focusing on 3--24 keV) made great progress, it also left significant uncertainties in characterizing the hard X-ray population, crucial for calibrating current population synthesis models. This paper presents an in-depth analysis of simulations of two extragalactic surveys (deep and wide) with the High-Energy X-ray Probe (HEX-P), each observed for 2 Ms. Applying established source detection techniques, we show that HEX-P surveys will reach a flux of $\sim$10$^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the 10-40 keV band, an order of magnitude fainter than current NuSTAR surveys. With the large sample of new hard X-ray detected sources ($\sim2000$), we showcase HEX-P's ability to resolve more than 80% of the CXB up to 40 keV into individual sources. The expected precision of HEX-P's resolved background measurement will allow us to distinguish between population synthesis models of SMBH growth. HEX-P leverages accurate broadband (0.5-40 keV) spectral analysis and the combination of soft and hard X-ray colors to provide obscuration constraints even for the fainter sources, with the overall objective of measuring the Compton-thick fraction. With unprecedented sensitivity in the 10--40 keV band, HEX-P will explore the hard X-ray emission from AGN to flux limits never reached before, thus expanding the parameter space for serendipitous discoveries. Consequently, it is plausible that new models will be needed to capture the population HEX-P will unveil., Comment: 28 pages, 13 Figures, 2 tables; submitted to Frontiers in Astronomy and Space Sciences; minor edits in the text, abstract and figure labels

We present a catalog of the millimeter-wave (mm-wave) continuum properties of 98 nearby ($z <$ 0.05) active galactic nuclei (AGNs) selected from the 70-month Swift/BAT hard X-ray catalog that have precisely determined X-ray spectral properties and subarcsec-resolution ALMA Band-6 (211--275 GHz) observations as of 2021 April. Due to the hard-X-ray ($>$ 10 keV) selection, the sample is nearly unbiased for obscured systems at least up to Compton-thick-level obscuration, and provides the largest number of AGNs with high physical resolution mm-wave data ($\lesssim$ 100--200 pc). Our catalog reports emission peak coordinates, spectral indices, and peak fluxes and luminosities at 1.3 mm (230 GHz). Additionally, high-resolution mm-wave images are provided. Using the images and creating radial surface brightness profiles of mm-wave emission, we identify emission extending from the central source and isolated blob-like emission. Flags indicating the presence of these emission features are tabulated. Among 90 AGNs with significant detections of nuclear emission, 37 AGNs ($\approx$ 41%) appear to have both or one of extended or blob-like components. We, in particular, investigate AGNs that show well-resolved mm-wave components and find that these seem to have a variety of origins (i.e., a jet, radio lobes, a secondary AGN, stellar clusters, a narrow line region, galaxy disk, active star-formation regions, and AGN-driven outflows), and some components have currently unclear origins., Comment: 49 pages, 7 figures, 3 tables, accepted for publication in ApJS

The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($\lambda$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region., Comment: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

Spectral energy distributions (SEDs) from X-ray to far-infrared (FIR) wavelengths are presented for a sample of 1246 X-ray luminous active galactic nuclei (AGN; $L_{0.5-10\rm{keV}}>10^{43}$ erg s$^{-1}$), with $z_{\rm{spec}}<1.2$, selected from Stripe 82X, COSMOS, and GOODS-N/S. The rest-frame SEDs show a wide spread ($\sim2.5$ dex) in the relative strengths of broad continuum features at X-ray, ultraviolet (UV), mid-infrared (MIR), and FIR wavelengths. A linear correlation (log-log slope of 0.7$\pm0.04$) is found between $L_{\rm{MIR}}$ and $L_{\rm{X}}$. There is significant scatter in the relation between the $L_{\rm{UV}}$ and $L_{\rm{X}}$ due to heavy obscuration, however the most luminous and unobscured AGN show a linear correlation (log-log slope of 0.8$\pm0.06$) in the relation above this scatter. The relation between $L_{\rm{FIR}}$ and $L_{\rm{X}}$ is predominantly flat, but with decreasing dispersion at $L_{\rm{X}}>10^{44}$ erg s$^{-1}$. The ratio between the "galaxy subtracted" bolometric luminosity and the intrinsic $L_{\rm{X}}$ increases from a factor of $\sim$$10-70$ from log $L_{\rm{bol}}/{\rm(erg\; s}^{-1})=44.5-46.5$. Characteristic SED shapes have been determined by grouping AGN based on relative strengths of the UV and MIR emission. The average $L_{1\mu\rm{m}}$ is constant for the majority of these SED shapes, while AGN with the strongest UV and MIR emission have elevated $L_{1\mu\rm{m}}$, consistent with the AGN emission dominating their SEDs at optical and NIR wavelengths. A strong correlation is found between the SED shape and both the $L_{\rm{X}}$ and $L_{\rm{bol}}$, such that $L_{\rm{bol}}/L_{\rm{X}}=20.4\pm1.8$, independent of the SED shape. This is consistent with an evolutionary scenario of increasing $L_{\rm{bol}}$ with decreasing obscuration as the AGN blows away circumnuclear gas., Comment: 24 pages, 20 figures, 4 appendix tables, Accepted to ApJ

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties., Comment: Accepted for publication in ApJ

Kilo-parsec scale hard ($>$ 3 keV) X-ray continuum and fluorescent Fe K$\alpha$ line emission has been recently discovered in nearby Compton-thick (CT) active galactic nuclei (AGN), which opens new opportunities to improve AGN torus modeling and investigate how the central supermassive black hole interacts with and impacts the host galaxy. Following a pilot Chandra survey of nearby CT AGN, we present in this paper the Chandra spatial analysis results of five uniformly selected non-CT but still heavily obscured AGN to investigate the extended hard X-ray emission by measuring the excess emission counts, excess fractions, and physical scales. Three of them show extended emission in the 3.0-7.0 keV band detected at $>$ 3$\sigma$ above the Chandra PSF with total excess fractions ranging from $\sim$8% - 20%. The extent of the hard emission ranges from at least $\sim$250 pc to 1.1 kpc in radius. We compare these new sources with CT AGN and find that CT AGN appear to be more extended in the hard band than the non-CT AGN. Similar to CT AGN, the amounts of extended hard X-ray emission relative to the total emission of these obscured AGN are not negligible. Together with other extended hard X-ray detected AGN in the literature, we further explore potential correlations between the extended hard X-ray component and AGN parameters. We also discuss the implications for torus modeling and AGN feedback. Considering potential contributions from X-ray binaries (XRBs) to the extended emission, we do not see strong XRB contamination in the overall sample., Comment: 14 pages, 12 figures; Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:2008.02175

We constrain the intrinsic Eddington ratio (\lamEdd ) distribution function for local AGN in bins of low and high obscuration (log NH <= 22 and 22 < log NH < 25), using the Swift-BAT 70-month/BASS DR2 survey. We interpret the fraction of obscured AGN in terms of circum-nuclear geometry and temporal evolution. Specifically, at low Eddington ratios (log lamEdd < -2), obscured AGN outnumber unobscured ones by a factor of ~4, reflecting the covering factor of the circum-nuclear material (0.8, or a torus opening angle of ~ 34 degrees). At high Eddington ratios (\log lamEdd > -1), the trend is reversed, with < 30% of AGN having log NH > 22, which we suggest is mainly due to the small fraction of time spent in a highly obscured state. Considering the Eddington ratio distribution function of narrow-line and broad-line AGN from our prior work, we see a qualitatively similar picture. To disentangle temporal and geometric effects at high lamEdd, we explore plausible clearing scenarios such that the time-weighted covering factors agree with the observed population ratio. We find that the low fraction of obscured AGN at high lamEdd is primarily due to the fact that the covering factor drops very rapidly, with more than half the time is spent with < 10% covering factor. We also find that nearly all obscured AGN at high-lamEdd exhibit some broad-lines. We suggest that this is because the height of the depleted torus falls below the height of the broad-line region, making the latter visible from all lines of sight., Comment: Accepted by ApJL

We present X-ray spectral analysis of XMM and Chandra observations in the 31.3 deg$^2$ Stripe-82X (S82X) field. Of the 6181 X-ray sources in this field, we analyze a sample of 2937 active galactic nuclei (AGN) with solid redshifts and sufficient counts determined by simulations. Our results show a population with median values of spectral index $\Gamma=1.94_{-0.39}^{+0.31}$, column density log$\,N_{\mathrm{H}}/\mathrm{cm}^{-2}=20.7_{-0.5}^{+1.2}$ and intrinsic, de-absorbed, 2-10 keV luminosity log$\,L_{\mathrm{X}}/\mathrm{erg\,s}^{-1}=44.0_{-1.0}^{+0.7}$, in the redshift range 0-4. We derive the intrinsic fraction of AGN that are obscured ($22\leq\mathrm{log}\,N_{\mathrm{H}}/\mathrm{cm}^{-2}<24$), finding a significant increase in the obscured AGN fraction with redshift and a decline with increasing luminosity. The average obscured AGN fraction is $57\pm4\%$ for log$\,L_{\mathrm{X}}/\mathrm{erg\,s}^{-1}>43$. This work constrains the AGN obscuration and spectral shape of the still uncertain high-luminosity and high-redshift regimes (log$\,L_{\mathrm{X}}/\mathrm{erg\,s}^{-1}>45.5$, $z>3$), where the obscured AGN fraction rises to $64\pm12\%$. We report a luminosity and density evolution of the X-ray luminosity function, with obscured AGN dominating at all luminosities at $z>2$ and unobscured sources prevailing at log$\,L_{\mathrm{X}}/\mathrm{erg\,s}^{-1}>45$ at lower redshifts. Our results agree with evolutionary models in which the bulk of AGN activity is triggered by gas-rich environments and in a downsizing scenario. Also, the black hole accretion density (BHAD) is found to evolve similarly to the star formation rate density, confirming the co-evolution between AGN and host-galaxy, but suggesting different time scales in their growing history. The derived BHAD evolution shows that Compton-thick AGN contribute to the accretion history of AGN as much as all other AGN populations combined., Comment: 31 pages, 35 figures, Accepted for publication in ApJ

To understand the origin of nuclear ($\lesssim$ 100 pc) millimeter-wave (mm-wave) continuum emission in active galactic nuclei (AGNs), we systematically analyzed sub-arcsec resolution Band-6 (211-275 GHz) ALMA data of 98 nearby AGNs ($z <$ 0.05) from the 70-month Swift/BAT catalog. The sample, almost unbiased for obscured systems, provides the largest number of AGNs to date with high mm-wave spatial resolution sampling ($\sim$ 1-200 pc), and spans broad ranges of 14-150 keV luminosity {$40 < \log[L_{\rm 14-150}/({\rm erg\,s^{-1}})] < 45$}, black hole mass [$5 < \log(M_{\rm BH}/M_\odot) < 10$], and Eddington ratio ($-4 < \log \lambda_{\rm Edd} < 2$). We find a significant correlation between 1.3 mm (230 GHz) and 14-150 keV luminosities. Its scatter is $\approx$ 0.36 dex, and the mm-wave emission may serve as a good proxy of the AGN luminosity, free of dust extinction up to $N_{\rm H} \sim 10^{26}$ cm$^{-2}$. While the mm-wave emission could be self-absorbed synchrotron radiation around the X-ray corona according to past works, we also discuss different possible origins of the mm-wave emission; AGN-related dust emission, outflow-driven shocks, and a small-scale ($<$ 200 pc) jet. The dust emission is unlikely to be dominant, as the mm-wave slope is generally flatter than expected. Also, due to no increase in the mm-wave luminosity with the Eddington ratio, a radiation-driven outflow model is possibly not the common mechanism. Furthermore, we find independence of the mm-wave luminosity on indicators of the inclination angle from the polar axis of the nuclear structure, which is inconsistent with a jet model whose luminosity depends only on the angle., Comment: 44 pages, 35 figures, 3 tables, accepted for publication in ApJ; modified format, added supplementary figure (Fig. 32)

We present the AGN catalog and optical spectroscopy for the second data release of the Swift BAT AGN Spectroscopic Survey (BASS DR2). With this DR2 release we provide 1425 optical spectra, of which 1181 are released for the first time, for the 858 hard X-ray selected AGN in the Swift BAT 70-month sample. The majority of the spectra (813/1425, 57%) are newly obtained from VLT/Xshooter or Palomar/Doublespec. Many of the spectra have both higher resolution (R>2500, N~450) and/or very wide wavelength coverage (3200-10000 A, N~600) that are important for a variety of AGN and host galaxy studies. We include newly revised AGN counterparts for the full sample and review important issues for population studies, with 44 AGN redshifts determined for the first time and 780 black hole mass and accretion rate estimates. This release is spectroscopically complete for all AGN (100%, 858/858) with 99.8% having redshift measurements (857/858) and 96% completion in black hole mass estimates of unbeamed AGN (outside the Galactic plane). This AGN sample represents a unique census of the brightest hard X-ray selected AGN in the sky, spanning many orders of magnitude in Eddington ratio (Ledd=10^-5-100), black hole mass (MBH=10^5-10^10 Msun), and AGN bolometric luminosity (Lbol=10^40-10^47 ergs/s)., Comment: 34 pages, 11 figures, 15 tables, published in ApJS as part of BASS DR2 special issue

The BAT AGN Spectroscopic Survey (BASS) is designed to provide a highly complete census of the key physical parameters of supermassive black holes (SMBHs) that power local active galactic nuclei (AGN) (z<0.3), including their bolometric luminosity, black hole mass, accretion rates, and line-of-sight gas obscuration, and the distinctive properties of their host galaxies (e.g., star formation rates, masses, and gas fractions). We present an overview of the BASS data release 2 (DR2), an unprecedented spectroscopic survey in spectral range, resolution, and sensitivity, including 1449 optical (3200-10000 A) and 233 NIR (1-2.5 um) spectra for the brightest 858 ultra-hard X-ray (14-195 keV) selected AGN across the entire sky and essentially all levels of obscuration. This release provides a highly complete set of key measurements (emission line measurements and central velocity dispersions), with 99.9% measured redshifts and 98% black hole masses estimated (for unbeamed AGN outside the Galactic plane). The BASS DR2 AGN sample represents a unique census of nearby powerful AGN, spanning over 5 orders of magnitude in AGN bolometric luminosity, black hole mass, Eddington ratio, and obscuration. The public BASS DR2 sample and measurements can thus be used to answer fundamental questions about SMBH growth and its links to host galaxy evolution and feedback in the local universe, as well as open questions concerning SMBH physics. Here we provide a brief overview of the survey strategy, the key BASS DR2 measurements, data sets and catalogs, and scientific highlights from a series of DR2-based works., Comment: 21 pages, 6 figures, 5 tables, published in ApJS as part of BASS DR2 special issue

The corona is an integral component of Active Galactic Nuclei (AGN) which produces the bulk of the X-ray emission above 1--2 keV. However, many of its physical properties and the mechanisms powering this emission remain a mystery. In particular, the temperature of the coronal plasma has been difficult to constrain for large samples of AGN, as constraints require high quality broadband X-ray spectral coverage extending above 10 keV in order to measure the high energy cutoff, which provides constraints on the combination of coronal optical depth and temperature. We present constraints on the coronal temperature for a large sample of Seyfert 1 AGN selected from the Swift/BAT survey using high quality hard X-ray data from the NuSTAR observatory combined with simultaneous soft X-ray data from Swift/XRT or XMM-Newton. When applying a physically-motivated, non-relativistic disk reflection model to the X-ray spectra, we find a mean coronal temperature kT $=$ 84$\pm$9 keV. We find no significant correlation between the coronal cutoff energy and accretion parameters such as the Eddington ratio and black hole mass. We also do not find a statistically significant correlation between the X-ray photon index, $\Gamma$, and Eddington ratio. This calls into question the use of such relations to infer properties of supermassive black hole systems., Comment: 16 pages, 10 figures. Accepted November 2021 for publication in ApJ; passed to production January 2022

Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below ${\lambda}1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at ${\lambda}1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses., Comment: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas., Comment: Accepted for Publication in ApJ

We present the joint Chandra, XMM-Newton and NuSTAR analysis of two nearby Seyfert galaxies, NGC 3081 and ESO 565-G019. These are the only two having Chandra data in a larger sample of ten low redshift ($z \le 0.05$), candidates Compton-thick Active Galactic Nuclei (AGN) selected in the 15-150 keV band with Swift-BAT that were still lacking NuSTAR data. Our spectral analysis, performed using physically-motivated models, provides an estimate of both the line-of-sight (l.o.s.) and average (N$_{H,S}$) column densities of the two torii. NGC 3081 has a Compton-thin l.o.s. column density N$_{H,z}$=[0.58-0.62] $\times 10^{24}$cm$^{-2}$, but the N$_{H,S}$, beyond the Compton-thick threshold (N$_{H,S}$=[1.41-1.78] $\times 10^{24}$cm$^{-2}$), suggests a "patchy" scenario for the distribution of the circumnuclear matter. ESO 565-G019 has both Compton-thick l.o.s. and N$_{H,S}$ column densities (N$_{H,z}>$2.31 $\times 10^{24}$cm$^{-2}$ and N$_{H,S} >$2.57 $\times 10^{24}$cm$^{-2}$, respectively). The use of physically-motivated models, coupled with the broad energy range covered by the data (0.6-70 keV and 0.6-40 keV, for NGC 3081 and ESO 565-G019, respectively) allows us to constrain the covering factor of the obscuring material, which is C$_{TOR}$=[0.63-0.82] for NGC 3081, and C$_{TOR}$=[0.39-0.65] for ESO 565-G019., Comment: 15 pages, 22 figures, Accepted for publication on ApJ

The nuclear obscurer of Active Galactic Nuclei (AGN) is poorly understood in terms of its origin, geometry and dynamics. We investigate whether physically motivated geometries emerging from hydro-radiative simulations can be differentiated with X-ray reflection spectroscopy. For two new geometries, the radiative fountain model of Wada (2012) and a warped disk, we release spectral models produced with the ray tracing code XARS. We contrast these models with spectra of three nearby AGN taken by NuSTAR and Swift/BAT. Along heavily obscured sight-lines, the models present different 4-20keV continuum spectra. These can be differentiated by current observations. Spectral fits of the Circinus Galaxy favor the warped disk model over the radiative fountain, and clumpy or smooth torus models. The necessary reflector (NH>10^25/cm^2) suggests a hidden population of heavily Compton-thick AGN amongst local galaxies. X-ray reflection spectroscopy is a promising pathway to understand the nuclear obscurer in AGN., Comment: Accepted in A&A. X-ray spectral models can be downloaded from https://github.com/JohannesBuchner/xars/blob/master/doc/warpeddisk.rst (warped disk) and https://github.com/JohannesBuchner/xars/blob/master/doc/wada.rst (radiative fountain model). The XARS code is at https://github.com/JohannesBuchner/xars/

Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes., Comment: 15 pages, 4 figures, published in PRD on May 19

We present torus modelling for the X-ray spectra of a nearby narrow-line Seyfert 1 galaxy Mrk 1239 ($z=0.0199$), based on archival Suzaku, NuSTAR and Swift observations. Our model suggests very soft intrinsic power-law continuum emission of $\Gamma\approx2.57$ in 2019 and $\Gamma\approx2.98$ in 2007. By applying a correction factor to the unabsorbed X-ray luminosity, we find that Mrk 1239 is accreting near or around the Eddington limit. Our best-fit spectral model also suggests a torus with a column density of $\log(N_{\rm H, ave}/$cm$^{-2})=25.0\pm0.2$ and a high covering factor of $0.90$ in Mrk 1239, indicating that this source is most likely to be viewed almost face-on with $i\approx26^{\circ}$. Our line of sight might cross the edge of the torus with $N_{\rm H, los}=2-5\times10^{23}$cm$^{-2}$. The high Eddington ratio and the high line-of-sight column density makes Mrk 1239 one of the AGNs that are close to the limit where wind may form near the edge of the torus due to high radiation pressure., Comment: 11 pages, 8 figures (accepted by MNRAS)

Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images., Comment: 29 pages, 11 figures, published in ApJ on May 3

The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources., Comment: Physical Review Letters

The recent Chandra discovery of extended $\sim$kpc-scale hard ($>$ 3 keV) X-ray emission in nearby Compton-thick (CT) active galactic nuclei (AGN) opens a new window to improving AGN torus modeling and investigating how the central super massive black hole interacts with and impacts the host galaxy. Since there are only a handful of detections so far, we need to establish a statistical sample to determine the ubiquity of the extended hard X-ray emission in CT AGN, and quantify the amount and extent of this component. In this paper, we present the spatial analysis results of a pilot Chandra imaging survey of 7 nearby ($0.006 < z < 0.013$) CT AGN selected from the Swift-BAT spectroscopic AGN survey. We find that five out of the seven CT AGN show extended emission in the 3-7 keV band detected at $>$ 3$\sigma$ above the Chandra PSF with $\sim$12% to 22% of the total emission in the extended components. ESO 137-G034 and NGC 3281 display biconical ionization structures with extended hard X-ray emission reaching kpc-scales ($\sim$ 1.9 kpc and 3.5 kpc in diameter). The other three show extended hard X-ray emission above the PSF out to at least $\sim$360 pc in radius. We find a trend that a minimum 3-7 keV count rate of 0.01 cts/s and total excess fraction $>$20% is required to detect a prominent extended hard X-ray component. Given that this extended hard X-ray component appears to be relatively common in this uniformly selected CT AGN sample, we further discuss the implications for torus modeling and AGN feedback., Comment: 20 pages, 20 figures, accepted for publication in ApJ

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., Comment: 13 pages, 9 figures. Accepted for publication in ApJ

We present a multi-frequency study of PKS J1222$+$0413 (4C$+$04.42), currently the highest redshift $\gamma$-ray emitting narrow-line Seyfert 1 ($\gamma$-NLS1). We assemble a broad spectral energy distribution (SED) including previously unpublished datasets: X-ray data obtained with the NuSTAR and Neil Gehrels Swift observatories; near-infrared, optical and UV spectroscopy obtained with VLT X-shooter; and multiband radio data from the Effelsberg telescope. These new observations are supplemented by archival data from the literature. We apply physical models to the broadband SED, parameterising the accretion flow and jet emission to investigate the disc-jet connection. PKS J1222$+$0413 has a much greater black hole mass than most other NLS1s, $M_\mathrm{BH}\approx2\times10^{8}$ M$_\odot$, similar to those found in flat spectrum radio quasars (FSRQs). Therefore this source provides insight into how the jets of $\gamma$-NLS1s relate to those of FSRQs., Comment: 18 pages, 10 figures and 6 tables. Accepted for publication in MNRAS

The corona is a key component of most luminous accreting black holes, carrying 5 - 30 % of the power and in non-jetted Active Galactic Nuclei (AGN), creating all the X-ray emission above $\simeq 1-2$ keV. Its emission illuminates the inner accretion disc, creating the atomic line-rich reflection spectrum used to diagnose and map the accretion flow and measure black hole spin. The corona is likely powered magnetically by the strong differential rotation of the accretion disc and it may be intimately related to relativistic jets. Recent work shows that many black hole coronae may be dominated by electron-positron pairs produced by photon-photon collisions in the compact coronal environment. Despite the corona being an integral component of AGN and black hole binary systems, much is still unknown about the nature of the corona, such as its geometry, location, and the physical mechanisms powering the emission. In this white paper we explore our current understanding of coronal properties, such as its temperature, obtained from measurements with existing hard X-ray telescopes such as NuSTAR, and discuss important questions to be addressed in the coming decade surrounding the nature of the corona. Hard X-ray observations will continue to dispel the mystery of coronae and open up this part of the quasar engine to full understanding., Comment: White Paper submitted to the Astro2020 Decadal Survey

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., Comment: 50 pages, Probe class mission concept study report submitted to NASA for Astro2020 Decadal Survey

Accreting supermassive black holes (SMBHs) can exhibit variable emission across the electromagnetic spectrum and over a broad range of timescales. The variability of active galactic nuclei (AGNs) in the ultraviolet and optical is usually at the few tens of per cent level over timescales of hours to weeks. Recently, rare, more dramatic changes to the emission from accreting SMBHs have been observed, including tidal disruption events, 'changing look' AGNs and other extreme variability objects. The physics behind the 're-ignition', enhancement and 'shut-down' of accretion onto SMBHs is not entirely understood. Here we present a rapid increase in ultraviolet-optical emission in the centre of a nearby galaxy, marking the onset of sudden increased accretion onto a SMBH. The optical spectrum of this flare, dubbed AT 2017bgt, exhibits a mix of emission features. Some are typical of luminous, unobscured AGNs, but others are likely driven by Bowen fluorescence - robustly linked here with high-velocity gas in the vicinity of the accreting SMBH. The spectral features and increased ultraviolet flux show little evolution over a period of at least 14 months. This disfavours the tidal disruption of a star as their origin, and instead suggests a longer-term event of intensified accretion. Together with two other recently reported events with similar properties, we define a new class of SMBH-related flares. This has important implications for the classification of different types of enhanced accretion onto SMBHs., Comment: Published in Nature Astronomy

We present the analysis of a sample of 35 candidate Compton thick (CT-) active galactic nuclei (AGNs) selected in the nearby Universe (average redshift ~0.03) with the Swift-BAT 100-month survey. All sources have available NuSTAR data, thus allowing us to constrain with unprecedented quality important spectral parameters such as the obscuring torus line-of-sight column density (N_{H, z}), the average torus column density (N_{H, tor}) and the torus covering factor (f_c). We compare the best-fit results obtained with the widely used MyTorus (Murphy et al. 2009) model with those of the recently published borus02 model (Balokovic et al. 2018) used in the same geometrical configuration of MyTorus (i.e., with f_c=0.5). We find a remarkable agreement between the two, although with increasing dispersion in N_{H, z} moving towards higher column densities. We then use borus02 to measure f_c. High-f_c sources have, on average, smaller offset between N_{H, z} and N_{H, tor} than low-f_c ones. Therefore, low f_c values can be linked to a "patchy torus" scenario, where the AGN is seen through an over-dense region in the torus, while high-f_c objects are more likely to be obscured by a more uniform gas distribution. Finally, we find potential evidence of an inverse trend between f_c and the AGN 2-10 keV luminosity, i.e., sources with higher f_c values have on average lower luminosities., Comment: 35 Pages, 23 Figures. Accepted for publication in ApJ

We present the analysis of the first NuSTAR observations ($\sim 220$ ks), simultaneous with the last SUZAKU observations ($\sim 50$ ks), of the active galactic nucleus of the bright Seyfert 1 galaxy Mrk 509. The time-averaged spectrum in the $1-79$ keV X-ray band is dominated by a power-law continuum ($\Gamma\sim 1.8-1.9$), a strong soft excess around 1 keV, and signatures of X-ray reflection in the form of Fe K emission ($\sim 6.4$ keV), an Fe K absorption edge ($\sim 7.1$ keV), and a Compton hump due to electron scattering ($\sim 20-30$ keV). We show that these data can be described by two very different prescriptions for the soft excess: a warm ($kT\sim 0.5-1$ keV) and optically thick ($\tau\sim10-20$) Comptonizing corona, or a relativistically blurred ionized reflection spectrum from the inner regions of the accretion disk. While these two scenarios cannot be distinguished based on their fit statistics, we argue that the parameters required by the warm corona model are physically incompatible with the conditions of standard coronae. Detailed photoionization calculations show that even in the most favorable conditions, the warm corona should produce strong absorption in the observed spectrum. On the other hand, while the relativistic reflection model provides a satisfactory description of the data, it also requires extreme parameters, such as maximum black hole spin, a very low and compact hot corona, and a very high density for the inner accretion disk. Deeper observations of this source are thus necessary to confirm the presence of relativistic reflection, and to further understand the nature of its soft excess., Comment: Accepted for publication in ApJ, 18 pages, 7 figures

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., Comment: Published in the Astrophysical Journal: Lanz, L. et al. 2019, ApJ, 870, 26; 19 pages; 9 figures

The NuSTAR observatory, with its high sensitivity in hard X-rays, has enabled detailed broadband modeling of the X-ray spectra of Active Galactic Nuclei (AGN), thereby allowing constraints to be placed on the high-energy cutoff of the X-ray coronal continuum. We investigate the spectral properties of a sample of 46 NuSTAR-observed Seyfert 1 AGN selected from the Swift/BAT 70-month hard X-ray survey. Our measurements of the high-energy cutoff of the continuum from modeling the NuSTAR X-ray spectra are used to map out the temperature - compactness ($\theta-l$) plane for AGN coronae. We find that most of the coronae lie clustered near the boundary for runaway pair production, suggesting that annihilation and pair production act to regulate the temperature of the corona. We discuss the implications of coronae whose high-energy cutoff may indicate a low coronal temperature on the heating and thermalization mechanisms in the corona., Comment: 10 pages, 8 figures, 3 tables, accepted for publication in the Astrophysical Journal

We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 \mu$m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of fifteen decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 which, when compared with the 2006 archival SDSS spectrum, reveals that the broad H$\beta$ emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes in the quasar luminosity over this timescale, and favor changes in the innermost regions of the accretion disk that occur on the thermal and heating/cooling front timescales. We discuss possible physical triggers that could cause these changes, and predict the multiwavelength signatures that could distinguish these physical scenarios., Comment: 9 pages, 4 figures, submitted to MNRAS

We present the first study of an Iwasawa-Taniguchi/X-ray Baldwin effect for Compton-thick active galactic nuclei (AGN). We report a statistically significant anti-correlation between the rest-frame equivalent width (EW) of the narrow core of the neutral Fe K$\alpha$ fluorescence emission line, ubiquitously observed in the reflection spectra of obscured AGN, and the mid-infrared 12$\,\mu$m continuum luminosity (taken as a proxy for the bolometric AGN luminosity). Our sample consists of 72 Compton-thick AGN selected from pointed and deep-field observations covering a redshift range of $z\sim0.0014-3.7$. We employ a Monte Carlo-based fitting method, which returns a Spearman's Rank correlation coefficient of $\rho=-0.28\pm0.12$, significant to 98.7% confidence. The best fit found is ${\rm log}({\rm EW}_{{\rm Fe\,K}\alpha})\,\propto\,-0.08\pm0.04\,{\rm log}(L_{12\,\mu{\rm m}})$, which is consistent with multiple studies of the X-ray Baldwin effect for unobscured and mildly obscured AGN. This is an unexpected result, as the Fe K$\alpha$ line is conventionally thought to originate from the same region as the underlying reflection continuum, which together constitute the reflection spectrum. We discuss the implications this could have if confirmed on larger samples, including a systematic underestimation of the line of sight X-ray obscuring column density and hence the intrinsic luminosities and growth rates for the most luminous AGN., Comment: Accepted for publication in MNRAS. 58 pages (22 page manuscript; 13 figures + 26 page appendix; 66 figures). This version update combines the manuscript and appendix

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., Comment: 9 pages, 3 figures. Accepted for publication in MNRAS

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., Comment: ApJS in press. The full tables will be available on ApJS and the on the BASS website [www.bass-survey.com] once the paper is published

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., Comment: Published in ApJ, 7 pages, 4 figures

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., Comment: Accepted ApJ, see www.bass-survey.com for data

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., Comment: Resubmitted to MNRAS after addressing referee's comments

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., Comment: 8 pages, 5 figures. Accepted for publication in ApJ

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., Comment: Accepted to ApJ

We present the results from a joint Suzaku/NuSTAR broad-band spectral analysis of 3C 390.3. The high quality data enables us to clearly separate the primary continuum from the reprocessed components allowing us to detect a high energy spectral cut-off ($E_\text{cut}=117_{-14}^{+18}$ keV), and to place constraints on the Comptonization parameters of the primary continuum for the first time. The hard over soft compactness is 69$_{-24}^{+124}$ and the optical depth 4.1$_{-3.6}^{+0.5}$, this leads to an electron temperature of $30_{-8}^{+32}$ keV. Expanding our study of the Comptonization spectrum to the optical/UV by studying the simultaneous Swift-UVOT data, we find indications that the compactness of the corona allows only a small fraction of the total UV/optical flux to be Comptonized. Our analysis of the reprocessed emission show that 3C 390.3 only has a small amount of reflection (R~0.3), and of that the vast majority is from distant neutral matter. However we also discover a soft X-ray excess in the source, which can be described by a weak ionized reflection component from the inner parts of the accretion disk. In addition to the backscattered emission, we also detect the highly ionized iron emission lines Fe XXV and Fe XXVI., Comment: accepted for publication in ApJ

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., Comment: About to be accepted in MNRAS upon minor grammar changes

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]., Comment: Submitted to ApJ; 23 pages (ApJ format); 11 figures and 3 tables; Comments welcomed!

The AllWISE processing pipeline has measured motions for all objects detected on WISE images taken between 2010 January and 2011 February. In this paper, we discuss new capabilities made to the software pipeline in order to make motion measurements possible, and we characterize the resulting data products for use by future researchers. Using a stringent set of selection criteria, we find 22,445 objects that have significant AllWISE motions, of which 3,525 have motions that can be independently confirmed from earlier 2MASS images yet lack any published motions in SIMBAD. Another 58 sources lack 2MASS counterparts and are presented as motion candidates only. Limited spectroscopic follow-up of this list has already revealed eight new L subdwarfs. These may provide the first hints of a "subdwarf gap" at mid-L types that would indicate the break between the stellar and substellar populations at low metallicities (i.e., old ages). Another object in the motion list -- WISEA J154045.67-510139.3 -- is a bright (J ~ 9 mag) object of type M6; both the spectrophotometric distance and a crude preliminary parallax place it ~6 pc from the Sun. We also compare our list of motion objects to the recently published list of 762 WISE motion objects from Luhman (2014). While these first large motion studies with WISE data have been very successful in revealing previously overlooked nearby dwarfs, both studies missed objects that the other found, demonstrating that many other nearby objects likely await discovery in the AllWISE data products., Comment: Accepted for publication in The Astrophysical Journal. Table 3 is a machine-readable ASCII file

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014., Comment: 20 pages, 12 figures, submitted to the Astrophysical Journal

We investigate the dichotomy in the radio loudness distribution of quasars by modelling their radio emission and various selection effects using a Monte Carlo approach. The existence of two physically distinct quasar populations, the radio-loud and radio-quiet quasars, is controversial and over the last decade a bimodal distribution of radio loudness of quasars has been both affirmed and disputed. We model the quasar radio luminosity distribution with simple unimodal and bimodal distribution functions. The resulting simulated samples are compared to a fiducial sample of 8,300 quasars drawn from the SDSS DR7 Quasar Catalog and combined with radio observations from the FIRST survey. Our results indicate that the SDSS-FIRST sample is best described by a radio loudness distribution which consists of two components, with 12+/-1 % of sources in the radio-loud component. On the other hand, the evidence for a local minimum in the loudness distribution (bimodality) is not strong and we find that previous claims for its existence were probably affected by the incompleteness of the FIRST survey close to its faint limit. We also investigate the redshift and luminosity dependence of the radio loudness distribution and find tentative evidence that at high redshift radio-loud quasars were rarer, on average "louder", and exhibited a smaller range in radio loudness. In agreement with other recent work, we conclude that the SDSS-FIRST sample strongly suggests that the radio loudness distribution of quasars is not a universal function, and that more complex models than presented here are needed to fully explain available observations., Comment: Accepted for publication in ApJ; 13 pages, 10 figures