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.

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.

The cosmic X-ray background (CXB), which peaks at an energy of ~30 keV, is produced primarily by emission from accreting supermassive black holes (SMBHs). The CXB therefore serves as a constraint on the integrated SMBH growth in the Universe and the accretion physics and obscuration in active galactic nuclei (AGNs). This paper gives an overview of recent progress in understanding the high-energy (>~10 keV) X-ray emission from AGNs and the synthesis of the CXB, with an emphasis on results from NASA's NuSTAR hard X-ray mission. We then discuss remaining challenges and open questions regarding the nature of AGN obscuration and AGN physics. Finally, we highlight the exciting opportunities for a next-generation, high-resolution hard X-ray mission to achieve the long-standing goal of resolving and characterizing the vast majority of the accreting SMBHs that produce the CXB.

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.

We report on the results of active galactic nuclei (AGNs) detection by NuSTAR performed in three extragalactic survey fields (COSMOS, UDS, ECDFS) in three hard bands, namely H1 (8-16 keV), H2 (16-24 keV) and VH (35-55 keV). The aggregated area of the surveys is $\sim 2.7$ deg$^2$. While a large number of sources is detected in the H1 band (72 at the $97\%$ level of reliability), the H2 band directly probing close to the peak of the Cosmic X-ray Background (CXB) returns four significant detections, and two tentative, although not significant, detections are found in the VH band. All the sources detected above 16 keV are also detected at lower energies. We compute the integral number counts for sources in such bands, which show broad consistency with population synthesis models of the CXB. We furthermore identify two Compton-thick AGNs, one in the COSMOS field, associated with a hard and faint Chandra source, and one in the UDS field, never detected in the X-ray band before. Both sources are at the same redshift $z \sim 1.25$, which shifts their Compton-hump into the H1 band, and were previously missed in the usually employed NuSTAR bands, confirming the potential of using the H1 band to discover obscured AGNs at $z > 1$ in deep surveys.

We discuss the spectral analysis of a sample of 63 Active Galactic Nuclei (AGN) detected above a limiting flux of $S(8-24$ keV$)=7\times10^{-14}$ erg/s/cm$^2$ in the multi-tiered NuSTAR Extragalactic Survey program. The sources span a redshift range z=0-2.1 (median $\langle$z$\rangle=$0.58). The spectral analysis is performed over the broad 0.5-24 keV energy range, combining NuSTAR with Chandra and/or XMM-Newton data and employing empirical and physically motivated models. This constitutes the largest sample of AGN selected at $ >10$ keV to be homogeneously spectrally analyzed at these flux levels. We study the distribution of spectral parameters such as photon index, column density ($N_{\rm H}$), reflection parameter ($R$) and 10-40 keV luminosity ($L_{X}$). Heavily obscured ($\log [N_{\rm H}/cm^{-2}]\ge23$) and Compton Thick (CT; $\log [N_{\rm H}/cm^{-2}]\ge24$) AGN constitute $\sim$25% (15-17 sources) and $\sim$2-3% ( 1-2 sources) of the sample, respectively. The observed $N_{\rm H}$ distribution fairly agrees with predictions of Cosmic X-ray Background population synthesis models (CXBPSM). We estimate the intrinsic fraction of AGN as a function of $N_{H}$, accounting for the bias against obscured AGN in a flux-selected sample. The fraction of CT AGN relative to $\log [N_{\rm H}/cm^{-2}]=20-24$ AGN is poorly constrainted, formally in the range 2-56% (90% upper limit of 66%). We derived a fraction ($f_{abs}$) of obscured AGN ($\log [N_{\rm H}/cm^{-2}]=22-24$) as a function of $L_{X}$ in agreement with CXBPSM and previous z < 1 X-ray determinations. Furthermore $f_{abs}$ at z=0.1-0.5 and $\log (L_{X}/erg$ $s^{-1})\approx43.6-44.3$ agrees with observational measurements/trends obtained over larger redshift intervals. We report a significant anti-correlation of $R$ with $L_{X}$ (confirmed by our companion paper on stacked spectra) with considerable scatter around the median $R$ values.

We present the results and the source catalog of the NuSTAR survey in the UKIDSS Ultra Deep Survey (UDS) field, bridging the gap in depth and area between NuSTAR's ECDFS and COSMOS surveys. The survey covers a $\sim 0.6$ deg$^2$ area of the field for a total observing time of $\sim$ 1.75 Ms, to a half-area depth of $\sim$ 155 ks corrected for vignetting at $3-24$ keV, and reaching sensitivity limits at half-area in the full ($3-24$ keV), soft ($3-8$ keV) and hard ($8-24$ keV) bands of $2.2 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, $1.0 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, and $2.7 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, respectively. A total of 67 sources are detected in at least one of the three bands, 56 of which have a robust optical redshift with a median of $\langle z\rangle \sim 1.1$. Through a broadband ($0.5-24$ keV) spectral analysis of the whole sample combined with the NuSTAR hardness ratios, we compute the observed Compton-thick (CT; $N_{\rm H} > 10^{24}$ cm$^{-2}$) fraction. Taking into account the uncertainties on each $N_{\rm H}$ measurement, the final number of CT sources is $6.8\pm1.2$. This corresponds to an observed CT fraction of $(11.5\pm2.0)\%$, providing a robust lower limit to the intrinsic fraction of CT AGN and placing constraints on cosmic X-ray background synthesis models.

Disk megamasers are a unique tool to study active galactic nuclei (AGN) sub-pc environment, and precisely measure some of their fundamental parameters. While the majority of disk megamasers are hosted in heavily obscured (i.e., Seyfert 2, Sy2) AGN, the converse is not true, and disk megamasers are very rarely found even in obscured AGN. The very low detection rate of such systems in Sy2 AGN could be due to the geometry of the maser beaming, which requires a strict edge-on condition. We explore some other fundamental factors which could play a role in a volume-limited survey of disk megamasers in Sy2 galaxies, most importantly the radio luminosity.

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.

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.

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.

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.

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.

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.

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.

Due to their heavily obscured central engines, the growth rate of Compton-thick (CT) active galactic nuclei (AGN) is difficult to measure. A statistically significant correlation between the Eddington ratio, {\lambda}$_{Edd}$, and the X-ray power-law index, {\Gamma}, observed in unobscured AGN offers an estimate of their growth rate from X-ray spectroscopy (albeit with large scatter). However, since X-rays undergo reprocessing by Compton scattering and photoelectric absorption when the line-of-sight to the central engine is heavily obscured, the recovery of the intrinsic {\Gamma} is challenging. Here we study a sample of local, predominantly Compton-thick megamaser AGN, where the black hole mass, and thus Eddington luminosity, are well known. We compile results on X-ray spectral fitting of these sources with sensitive high-energy (E> 10 keV) NuSTAR data, where X-ray torus models which take into account the reprocessing effects have been used to recover the intrinsic {\Gamma} values and X-ray luminosities, L$_X$. With a simple bolometric correction to L$_X$ to calculate {\lambda}$_{Edd}$, we find a statistically significant correlation between {\Gamma} and {\lambda}$_{Edd}$ (p = 0.007). A linear fit to the data yields {\Gamma} = (0.41$\pm$0.18)log$_{10}${\lambda}$_{Edd}$+(2.38$\pm$ 0.20), which is statistically consistent with results for unobscured AGN. This result implies that torus modeling successfully recovers the intrinsic AGN parameters. Since the megamasers have low-mass black holes (M$_{BH}\approx10^6-10^7$ M$_{sol}$) and are highly inclined, our results extend the {\Gamma}-{\lambda}$_{Edd}$ relationship to lower masses and argue against strong orientation effects in the corona, in support of AGN unification. Finally this result supports the use of {\Gamma} as a growth-rate indicator for accreting black holes, even for Compton-thick AGN.

PG1247+267 is one of the most luminous known quasars at $z\sim2$ and is a strongly super-Eddington accreting SMBH candidate. We obtained NuSTAR data of this intriguing source in December 2014 with the aim of studying its high-energy emission, leveraging the broad band covered by the new NuSTAR and the archival XMM-Newton data. Several measurements are in agreement with the super-Eddington scenario for PG1247+267: the soft power law ($\Gamma=2.3\pm0.1$); the weak ionized Fe emission line and a hint of the presence of outflowing ionized gas surrounding the SMBH. The presence of an extreme reflection component is instead at odds with the high accretion rate proposed for this quasar. This can be explained with three different scenarios; all of them are in good agreement with the existing data, but imply very different conclusions: i) a variable primary power law observed in a low state, superimposed on a reflection component echoing a past, higher flux state; ii) a power law continuum obscured by an ionized, Compton thick, partial covering absorber; and iii) a relativistic disk reflector in a lamp-post geometry, with low coronal height and high BH spin. The first model is able to explain the high reflection component in terms of variability. The second does not require any reflection to reproduce the hard emission, while a rather low high-energy cutoff of $\sim100$ keV is detected for the first time in such a high redshift source. The third model require a face-on geometry, which may affect the SMBH mass and Eddington ratio measurements. Deeper X-ray broad-band data are required in order to distinguish between these possibilities.

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.