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.

Hot Dust-Obscured Galaxies (Hot DOGs) are among the most luminous galaxies in the Universe. Powered by highly obscured, possibly Compton-thick, active galactic nuclei (AGNs), Hot DOGs are characterized by SEDs that are very red in the mid-IR yet dominated by the host galaxy stellar emission in the UV and optical. An earlier study identified a sub-sample of Hot DOGs with significantly enhanced UV emission. One target, W0204-0506, was studied in detail and, based on Chandra observations, it was concluded that the enhanced emission was most likely due to either extreme unobscured star-formation (${\rm SFR}>1000~M_{\odot}~\rm yr^{-1}$) or to light from the highly obscured AGN scattered by gas or dust into our line of sight. Here, we present a follow-up study of W0204-0506 as well as two more Hot DOGs with excess UV emission. For the two new objects we obtained Chandra/ACIS-S observations, and for all three targets we obtained HST/WFC3 F555W and F160W imaging. We conclude that the excess UV emission is primarily dominated by light from the central highly obscured, hyper-luminous AGN that has been scattered into our line of sight. We cannot rule out, however, that star-formation may significantly contribute to the UV excess of W0204-0506.

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generally observed at an X-ray luminosity between $10^{39}$ and $10^{40}$ erg s$^{-1}$, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2-d binary with a projected semi-major axis $a_\mathrm{X} \sin i \simeq$ 28 lt-s. For a neutron star (NS) of 1.4 $M_{\odot}$, this implies a lower limit on the companion mass of 8 $M_{\odot}$, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by $\simeq$0.4 ms in the 31 d spanned by our May -- June 2018 observations, corresponding to a spin-up rate $\dot{P} \simeq -1.5\times10^{-10}\text{s s}^{-1}$. In an archival 2005 XMM-Newton exposure, we measured a spin period of $\sim$3.3 s, indicating a secular spin-up of $\dot{P}_{\mathrm{sec}}\simeq -10^{-9}\text{ s s}^{-1}$, a value in the range of other known PULXs. Our findings suggest that the system consists of an OB giant and a moderately magnetic (dipole field component in the range $10^{12}$ G $\lesssim B_{\mathrm{dip}}\lesssim 10^{13}$G) accreting NS with weakly beamed emission ($1/12\lesssim b\lesssim1/4$).

A rip in the MLI at the exit aperture of OMA, the NuSTAR optic aligned with detector focal plane module FPMA, has resulted in an increased photon flux through OMA that has manifested itself as a low energy excess. Overall, the MLI coverage has decreased by 10%, but there is an additional time-varying component, which occasionally causes the opening to increase by up to 20%. We address the problem with a calibration update, and in this paper, we describe the attributes of the problem, the implications it has on data analysis, and the solution.

We present results from the major coordinated X-ray observing program on the ULX NGC 1313 X-1 performed in 2017, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$, focusing on the evolution of the broadband ($\sim$0.3-30.0 keV) continuum emission. Clear and unusual spectral variability is observed, but this is markedly suppressed above $\sim$10-15 keV, qualitatively similar to the ULX Holmberg IX X-1. We model the multi-epoch data with two-component accretion disc models designed to approximate super-Eddington accretion, allowing for both a black hole and a neutron star accretor. With regards to the hotter disc component, the data trace out two distinct tracks in the luminosity-temperature plane, with larger emitting radii and lower temperatures seen at higher observed fluxes. Despite this apparent anti-correlation, each of these tracks individually shows a positive luminosity-temperature relation. Both are broadly consistent with $L\propto{T}^{4}$, as expected for blackbody emission with a constant area, and also with $L\propto{T}^{2}$, as may be expected for an advection-dominated disc around a black hole. We consider a variety of possibilities for this unusual behaviour. Scenarios in which the innermost flow is suddenly blocked from view by outer regions of the super-Eddington disc/wind can explain the luminosity-temperature behaviour, but are difficult to reconcile with the lack of strong variability at higher energies, assuming this emission arises from the most compact regions. Instead, we may be seeing evidence for further radial stratification of the accretion flow than is included in the simple models considered, with a combination of winds and advection resulting in the suppressed high-energy variability.

We report the discovery of a third ULX in NGC 925 (ULX-3), detected in November 2017 by Chandra at a luminosity of $L_{\rm X} = (7.8\pm0.8)\times10^{39}$ erg s$^{-1}$. Examination of archival data for NGC 925 reveals that ULX-3 was detected by Swift at a similarly high luminosity in 2011, as well as by XMM-Newton in January 2017 at a much lower luminosity of $L_{\rm X} = (3.8\pm0.5)\times10^{38}$ erg s$^{-1}$. With an additional Chandra non-detection in 2005, this object demonstrates a high dynamic range of flux of factor >26. In its high-luminosity detections, ULX-3 exhibits a hard power-law spectrum with $\Gamma=1.6\pm0.1$, whereas the XMM-Newton detection is slightly softer, with $\Gamma=1.8^{+0.2}_{-0.1}$ and also well-fitted with a broadened disc model. The long-term light curve is sparsely covered and could be consistent either with the propeller effect or with a large-amplitude superorbital period, both of which are seen in ULXs, in particular those with neutron star accretors. Further systematic monitoring of ULX-3 will allow us to determine the mechanism by which ULX-3 undergoes its extreme variability and to better understand the accretion processes of ULXs.

Ultraluminous X-ray sources (ULXs) are a class of accreting compact objects with X-ray luminosities above 1e39 erg/s. The ULX population counts several hundreds objects but only a minor fraction is well studied. Here we present a detailed analysis of all ULXs hosted in the galaxy NGC 7456. It was observed in X-rays only once in the past (in 2005) by XMM-Newton, but the observation was short and strongly affected by high background. In 2018, we obtained a new, deeper (~90 ks) XMM-Newton observation that allowed us to perform a detailed characterization of the ULXs hosted in the galaxy. ULX-1 and ULX-2, the two brightest objects (Lx~(6-10)e39 erg/s), have spectra that can be described by a two-thermal component model as often found in ULXs. ULX-1 shows also one order of magnitude in flux variability on short-term timescales (hundreds to thousand ks). The other sources (ULX-3 and ULX-4) show flux changes of at least an order of magnitude, and these objects may be candidate transient ULXs although longer X-ray monitoring or further studies are required to ascribe them to the ULX population. In addition, we found a previously undetected source that might be a new candidate ULX (labelled as ULX-5) with a luminosity of ~1e39 erg/s and hard power-law spectral shape, whose nature is still unclear and for which a background Active Galactic Nucleus cannot be excluded. We discuss the properties of all the ULXs in NGC 7456 within the framework of super-Eddington accretion onto stellar mass compact objects. Although no pulsations were detected, we cannot exclude that the sources host neutron stars.

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$\alpha$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $\Gamma$, and the X-ray-to-optical power slopes, $\alpha_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

M82 hosts two well-known ultraluminous X-ray sources (ULXs). X-1, an intermediate-mass black hole (IMBH) candidate, and X-2, an ultraluminous X-ray pulsar (ULXP). Here we present a broadband X-ray spectral analysis of both sources based on ten observations made simultaneously with Chandra and NuSTAR. Chandra provides the high spatial resolution to resolve the crowded field in the 0.5--8 keV band, and NuSTAR provides the sensitive hard X-ray spectral data, extending the bandpass of our study above 10 keV. The observations, taken in the period 2015--2016, cover a period of flaring from X-1, allowing us to study the spectral evolution of this source with luminosity. During four of these observations, X-2 was found to be at a low flux level, allowing an unambiguous view of the emission from X-1. We find that the broadband X-ray emission from X-1 is consistent with that seen in other ULXs observed in detail with NuSTAR, with a spectrum that includes a broadened disk-like component and a high-energy tail. We find that the luminosity of the disk scales with inner disk temperature as L~T^-3/2 contrary to expectations of a standard accretion disk and previous results. These findings rule out a thermal state for sub-Eddington accretion and therefore do not support M82 X-1 as an IMBH candidate. We also find evidence that the neutral column density of the material in the line of sight increases with L$_X$, perhaps due to an increased mass outflow with accretion rate. For X-2, we do not find any significant spectral evolution, but we find the spectral parameters of the phase-averaged broadband emission are consistent with the pulsed emission at the highest X-ray luminosities.

We present the results from a monitoring campaign made with the Neil Gehrels Swift Observatory of the M51 galaxies, which contain several variable ultraluminous X-ray sources (ULXs). The ongoing campaign started in May 2018, and we report here on $\sim1.5$ years of observations. The campaign, which consists of 105 observations, has a typical cadence of 3--6 days, and has the goal of determining the long-term X-ray variability of the ULXs. Two of the most variable sources were ULX7 and ULX8, both of which are known to be powered by neutron stars that are exceeding their isotropic Eddington luminosities by factors of up to 100. This is further evidence that neutron star powered ULXs are the most variable. Our two main results are, first, that ULX7 exhibits a periodic flux modulation with a period of 38 days varying over a magnitude and a half in flux from peak to trough. Since the orbital period of the system is known to be 2 days, the modulation is super-orbital, which is a near-ubiquitous property of ULX pulsars. Secondly we identify a new transient ULX, M51 XT-1, the onset of which occurred during our campaign, reaching a peak luminosity of $\sim10^{40}$ erg s$^{-1}$, before gradually fading over the next $\sim200$ days until it slipped below the detection limit of our observations. Combined with the high-quality Swift/XRT lightcurve of the transient, serendipitous observations made with Chandra and XMM-Newton provide insights into the onset and evolution of a likely super-Eddington event.

Ultraluminous X-ray Sources (ULXs) provide a unique opportunities to probe the geometry and energetics of super-Eddington accretion. The radiative processes involved in super-Eddington accretion are not well understood, and so studying correlated variability between different energy bands can provide insights into the causal connection between different emitting regions. We present a spectral-timing analysis of NGC 1313 X-1 from a recent XMM-Newton campaign. The spectra can be decomposed into two thermal-like components, the hotter of which may originate from the inner accretion disc, and the cooler from an optically thick outflow. We find correlated variability between hard (2-10 keV) and soft (0.3-2 keV) bands on kilosecond timescales, and find a soft lag of ~150 seconds. The covariance spectrum suggests that emission contributing to the lags is largely associated with the hotter of the two thermal-like components, likely originating from the inner accretion flow. This is only the third ULX to exhibit soft lags. The lags range over three orders of magnitude in amplitude, but all three are ~5 to ~20 percent of the corresponding characteristic variability timescales. If these soft lags can be understood in the context of a unified picture of ULXs, then lag timescales may provide constraints on the density and extent of radiatively-driven outflows.

M82 X-2 is the first pulsating ultraluminous X-ray source (PULX)to be identified. Since the discovery in 2014, NuSTAR has observed the M82 field 15 times throughout 2015 and 2016. In this paper, we report the results of pulsation searches in all these datasets, and find only one new detection. This new detection allows us to refine the orbital period of the source and measure an average spin down rate between 2014 and 2016 of 5x10^-11 Hz/s, which is in contrast to the strong spin up seen during the 2014 observations and represent the first detection of spin down in a PULX system. Thanks to the improved orbital solution allowed by this new detection, we are also able to detect pulsations in additional segments of the original 2014 dataset. We find a glitch superimposed on the very strong and variable spin-up already reported, the first positive glitch identified in a PULX system. We discuss the new findings in the context of current leading models for PULXs.

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.

SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS stellar atmosphere with $T_{\rm eff} = 3650 - 3900$ K and $\log(L_{\rm bol}/L_{\odot}) = 4.25\pm0.10$, which yields a stellar radius $R = 310 \pm 70 R_{\odot}$. To fit the full spectrum, two additional components are required: a blue excess that can be fitted either by a hot blackbody (T $\gtrsim 20,000$ K) or a power law (spectral index $\alpha \approx 4$) and is likely due to X-ray emission reprocessed in the outer accretion disk or the donor star; and a red excess that is well fitted by a blackbody with a temperature of $\sim 1100$ K, and is likely due to warm dust in the vicinity of SN2010da. The presence of a red supergiant in this system implies an orbital period of at least 0.8-2.1 years, assuming Roche lobe overflow. Given the large donor-to-compact object mass ratio, orbital modulations of the radial velocity of the red supergiant are likely undetectable. However, the radial velocity amplitude of the neutron star is large enough (up to 40-60 km s$^{-1}$) to potentially be measured in the future, unless the system is viewed at a very unfavorable inclination.

We present a unification model for a clumpy obscurer in active galactic nuclei (AGN) and investigate the properties of the resulting X-ray spectrum. Our model is constructed to reproduce the column density distribution of the AGN population and cloud eclipse events in terms of their angular sizes and frequency. We developed and release a generalised Monte Carlo X-ray radiative transfer code, XARS, to compute X-ray spectra of obscurer models. The geometry results in strong Compton scattering, causing soft photons to escape also along Compton-thick sight lines. This makes our model spectra very similar to the Brightman & Nandra TORUS model. However, only if we introduce an additional Compton-thick reflector near the corona, we achieve good fits to NuSTAR spectra. This additional component in our model can be interpreted as part of the dust-free broad-line region, an inner wall or rim, or a warped disk. It cannot be attributed to a simple disk because the reflector must simultaneously block the line of sight to the corona and reflect its radiation. We release our model as an Xspec table model and present corresponding CLUMPY infrared spectra, paving the way for self-consistent multi-wavelength analyses.

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.

Two recent observations of the nearby galaxy NGC 6946 with NuSTAR, one simultaneous with an XMM-Newton observation, provide an opportunity to examine its population of bright accreting sources from a broadband perspective. We study the three known ultraluminous X-ray sources (ULXs) in the galaxy, and find that ULX-1 and ULX-2 have very steep power-law spectra with $\Gamma=3.6^{+0.4}_{-0.3}$ in both cases. Their properties are consistent with being super-Eddington accreting sources with the majority of their hard emission obscured and down-scattered. ULX-3 (NGC 6946 X-1) is significantly detected by both XMM-Newton and NuSTAR at $L_{\rm X}=(6.5\pm0.1)\times10^{39}$ erg s$^{-1}$, and has a power-law spectrum with $\Gamma=2.51\pm0.05$. We are unable to identify a high-energy break in its spectrum like that found in other ULXs, but the soft spectrum likely hinders our ability to detect one. We also characterise the new source, ULX-4, which is only detected in the joint XMM-Newton and NuSTAR observation, at $L_{\rm X}=(2.27\pm0.07)\times10^{39}$ erg s$^{-1}$, and is absent in a Chandra observation ten days later. It has a very hard cut-off power-law spectrum with $\Gamma=0.7\pm0.1$ and $E_{\rm cut}=11^{+9}_{-4}$ keV. We do not detect pulsations from ULX-4, but its transient nature can be explained either as a neutron star ULX briefly leaving the propeller regime or as a micro-tidal disruption event induced by a stellar-mass compact object.

SS433 is the only Galactic binary system known to accrete at highly super-critical rates, analogous to tidal disruption events, and needed to explain the mass of some high redshift quasars. Probing the inner regions of SS433 in the X-rays is crucial to understanding this system, and super-critical accretion in general, but has not yet been possible due to obscuration. NuSTAR observed SS433 in the hard X-ray band across multiple phases of its super-orbital precession period. Spectral-timing tools have allowed us to confirm that the hard X-ray emission from the inner regions is scattered towards us by the walls of the wind-cone. By comparing to numerical models, we determine an intrinsic X-ray luminosity of >= 3x10^37 erg/s and that, if viewed face on, the apparent luminosity would be > 1x10^39 erg/s, confirming its long-suspected nature as an ultraluminous X-ray source (ULX). A lag due to absorption by Fe XXV/XXVI in outflowing material travelling at least 0.14-0.29c matches absorption lines seen in ULXs and - in future - will allow us to map a super-critical outflow for the first time.

Over the last decades X-ray spectroscopy has proven to be a powerful tool for the estimation of black hole spin and several other key parameters in dozens of AGN and black hole X-ray binaries. In this White Paper, we discuss the observational and theoretical challenges expected in the exploration, discovery, and study of astrophysical black holes in the next decade. We focus on the case of accreting black holes and their electromagnetic signatures, with particular emphasis on the measurement of the relativistic reflection component in their X-ray spectra.

With the recent discoveries of massive and highly luminous quasars at high redshifts ($z\sim7$; e.g. Mortlock et al. 2011), the question of how black holes (BHs) grow in the early Universe has been cast in a new light. In order to grow BHs with $M_{\rm BH} > 10^9$ M$_{\odot}$ by less than a billion years after the Big Bang, mass accretion onto the low-mass seed BHs needs to have been very rapid (Volonteri & Rees, 2005). Indeed, for any stellar remnant seed, the rate required would need to exceed the Eddington limit. This is the point at which the outward force produced by radiation pressure is equal to the gravitational attraction experienced by the in-falling matter. In principle, this implies that there is a maximum luminosity an object of mass $M$ can emit; assuming spherical accretion and that the opacity is dominated by Thompson scattering, this Eddington luminosity is $L_{\rm{E}} = 1.38 \times 10^{38} (M/M_{\odot})$ erg s$^{-1}$. In reality, it is known that this limit can be violated, due to non-spherical geometry or various kinds of instabilities. Nevertheless, the Eddington limit remains an important reference point, and many of the details of how accretion proceeds above this limit remain unclear. Understanding how this so-called super-Eddington accretion occurs is of clear cosmological importance, since it potentially governs the growth of the first supermassive black holes (SMBHs) and the impact this growth would have had on their host galaxies (`feedback') and the epoch of reionization, as well as improving our understanding of accretion physics more generally.

The reported discovery of a cyclotron resonance scattering feature (CRSF) in the spectrum of M51 ULX-8 may provide an important clue as to the nature of the magnetic field in those ultraluminous X-ray sources hosting neutron stars. In this paper we present the covariance (linearly correlated variability) spectrum of M51 ULX-8 on long (> 2000s) timescales. This allows us to unambiguously decompose the spectrum which requires multiple components in order to explain the broad-band emission. Having a clearer picture of the spectral decomposition leads to various tests for the dipole field strength of the neutron star which can be extended to other ULXs when certain criteria are met. In the case of M51 ULX-8, we rule out a very strong (~10^15 G) dipole solution with either a sub- or super-critical disc. Instead, our tests indicate an upper limit on the dipole field of ~10^12 G, and a classical super-critical inflow, similar to that inferred in other ULXs found to harbour neutron stars, although we do not rule out the presence of an additional, strong (~10^15 G) multipole field falling off steeply with distance from the neutron star.

Ultraluminous X-ray (ULX) pulsars are a new class of object powered by apparent super-critical accretion onto magnetized neutron stars. Three sources in this class identified so far; M82 X-2, NGC 5907 ULX-1 and NGC 7793 P13, have been found to have two properties in common; $\sim1$-s spin periods, and for NGC 5907 ULX-1 and NGC 7793 P13 periodic X-ray flux modulations on timescales of $\sim60-80$ days. M82 X-2 resides in a crowded field that includes the ULX M82 X-1 separated from X-2 by 5", as well as other bright point sources. A 60-day modulation has been observed from the region but the origin has been difficult to identify; both M82 X-1 and X-2 have been suggested as the source. In this paper we present the analysis of a systematic monitoring campaign by Chandra, the only X-ray telescope capable of resolving the crowded field. From a simple Lomb-Scargle periodogram analysis and a more sophisticated Gaussian Process analysis we find that only X-2 exhibits a periodic signal around 60 days supporting previous claims that it is the origin. We also construct a phase-averaged flux profile of the modulations from higher cadence Swift/XRT data and find that the flux variations in the Chandra data are fully consistent with the flux profile. Since the orbit of the neutron star and its companion is known to be 2.5 days, the $\sim60$-day period must be super-orbital in origin. The flux of the modulations varies by a factor of $\sim$100 from minimum to maximum, with no evidence for spectral variations, making the origin difficult to explain.

We report on our search for the optical counterparts of two ultraluminous X-ray pulsars with known orbital periods, M82 X-2 and NGC 5907 X-1, in new and archival HST observations, in an effort to characterize the donor stars in these systems. We detect five near-infrared sources consistent with the position of M82 X-2 that are too bright to be single stars. We also detect seven sources in the WFC3/UVIS F336W image whose photometry matches that of 10-15 M$_\odot$ stars turning off the main sequence. Such stars have densities consistent with the properties of the donor star of M82 X-2 as inferred from X-ray timing analysis, although it is also possible that the donor is a lower mass star below our detection limit or that there is a significant contribution from the accretion disc to the optical emission. We detect three candidate counterparts to NGC 5907 X-1 in the near-infrared. All of these are too bright to be the donor star of the ULX, which based on its orbital period is a red giant. The high background at the location of NGC 5907 X-1 precludes us from detecting this expected donor star. The recently discovered NGC 5907 ULX-2 also falls within the field of view of the near-infrared imaging; we detect four sources in the error circle, with photometry that matches AGB stars. The star suggested to be the counterpart of NGC 5907 ULX-2 by Pintore et al. (2018) falls outside our 2-$\sigma$ error circle.

We present the analysis of a sample of 35 candidate Compton thick (CT-) active galactic nuclei (AGNs) selected in the nearby Universe (average redshift <z>~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.

The existence of a large population of Compton thick (CT, $N_{H}>10^{24} cm^{-2}$) AGN is a key ingredient of most Cosmic X-ray background synthesis models. However, direct identification of these sources, especially at high redshift, is difficult due to the flux suppression and complex spectral shape produced by CT obscuration. We explored the Chandra COSMOS Legacy point source catalog, comprising 1855 sources, to select via X-ray spectroscopy, a large sample of CT candidates at high redshift. Adopting a physical model to reproduce the toroidal absorber, and a Monte-Carlo sampling method, we selected 67 individual sources with >5% probability of being CT, in the redshift range $0.04<z<3.5$. The sum of the probabilities above $N_{H}>10^{24} cm^{-2}$, gives a total of 41.9 effective CT, corrected for classification bias. We derive number counts in the 2-10 keV band in three redshift bins. The observed logN-logS is consistent with an increase of the intrinsic CT fraction ($f_{CT}$) from $\sim0.30$ to $\sim0.55$ from low to high redshift. When rescaled to a common luminosity (log(L$_{\rm X}$/erg/s)$=44.5$) we find an increase from $f_{CT}=0.19_{-0.06}^{+0.07}$ to $f_{CT}=0.30_{-0.08}^{+0.10}$ and $f_{CT}=0.49_{-0.11}^{+0.12}$ from low to high z. This evolution can be parametrized as $f_{CT}=0.11_{-0.04}^{+0.05}(1+z)^{1.11\pm0.13}$. Thanks to HST-ACS deep imaging, we find that the fraction of CT AGN in mergers/interacting systems increases with luminosity and redshift and is significantly higher than for non-CT AGN hosts.

We present the X-UDS survey, a set of wide and deep Chandra observations of the Subaru-XMM Deep/UKIDSS Ultra Deep Survey (SXDS/UDS) field. The survey consists of 25 observations that cover a total area of 0.33 deg$^{2}$. The observations are combined to provide a nominal depth of ~600 ksec in the central 100 arcmin$^{2}$ region of the field that has been imaged with Hubble/WFC3 by the CANDELS survey and $\sim$200 ksec in the remainder of the field. In this paper, we outline the survey's scientific goals, describe our observing strategy, and detail our data reduction and point source detection algorithms. Our analysis has resulted in a total of 868 band-merged point sources detected with a false-positive Poisson probability of $<1\times10^{-4}$. In addition, we present the results of an X-ray spectral analysis and provide best-fitting neutral hydrogen column densities, $N_{\rm H}$, as well as a sample of 51 Compton-thick active galactic nucleus candidates. Using this sample, we find the intrinsic Compton-thick fraction to be 30-35% over a wide range in redshift ($z=0.1-3$), suggesting the obscured fraction does not evolve very strongly with epoch. However, if we assume that the Compton-thick fraction is dependent on luminosity, as is seen for Compton-thin sources, then our results are consistent with a rise in the obscured fraction out to $z\sim3$. Finally, an examination of the host morphologies of our Compton-thick candidates shows a high fraction of morphological disturbances, in agreement with our previous results. All data products described in this paper are made available via a public website.

We present a broadband X-ray spectral analysis of the M51 system, including the dual active galactic nuclei (AGN) and several off-nuclear point sources. Using a deep observation by NuSTAR, new high-resolution coverage of M51b by Chandra, and the latest X-ray torus models, we measure the intrinsic X-ray luminosities of the AGN in these galaxies. The AGN of M51a is found to be Compton thick, and both AGN have very low accretion rates ($\lambda_{\rm Edd} <10^{-4}$). The latter is surprising considering that the galaxies of M51 are in the process of merging, which is generally predicted to enhance nuclear activity. We find that the covering factor of the obscuring material in M51a is $0.26 \pm 0.03$, consistent with the local AGN obscured fraction at $L_{\rm X}\sim 10^{40}$ erg s$^{-1}$. The substantial obscuring column does not support theories that the torus, presumed responsible for the obscuration, disappears at these low accretion luminosities. However, the obscuration may have resulted from the gas infall driven by the merger rather than the accretion process. We report on several extra-nuclear sources with $L_{\rm X}>10^{39}$ erg s$^{-1}$ and find that a spectral turnover is present below 10 keV in most such sources, in line with recent results on ultraluminous X-ray sources.

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.

Ultraluminous pulsars are a definite proof that persistent super-Eddington accretion occurs in nature. They support the scenario according to which most Ultraluminous X-ray Sources (ULXs) are super-Eddington accretors of stellar mass rather than sub-Eddington intermediate mass black holes. An important prediction of theories of supercritical accretion is the existence of powerful outflows of moderately ionized gas at mildly relativistic speeds. In practice, the spectral resolution of X-ray gratings such as RGS onboard XMM-Newton is required to resolve their observational signatures in ULXs. Using RGS, outflows have been discovered in the spectra of 3 ULXs (none of which are currently known to be pulsars). Most recently, the fourth ultraluminous pulsar was discovered in NGC 300. Here we report detection of an ultrafast outflow (UFO) in the X-ray spectrum of the object, with a significance of more than 3{\sigma}, during one of the two simultaneous observations of the source by XMM-Newton and NuSTAR in December 2016. The outflow has a projected velocity of 65000 km/s (0.22c) and a high ionisation factor with a log value of 3.9. This is the first direct evidence for a UFO in a neutron star ULX and also the first time that this its evidence in a ULX spectrum is seen in both soft and hard X-ray data simultaneously. We find no evidence of the UFO during the other observation of the object, which could be explained by either clumpy nature of the absorber or a slight change in our viewing angle of the accretion flow.

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.

We present a timing analysis of multiple XMM-Newton and NuSTAR observations of the ultra-luminous pulsar NGC 7793 P13 spread over its 65d variability period. We use the measured pulse periods to determine the orbital ephemeris, confirm a long orbital period with P_orb = 63.9 (+0.5,-0.6) d, and find an eccentricity of e <= 0.15. The orbital signature is imprinted on top of a secular spin-up, which seems to get faster as the source becomes brighter. We also analyse data from dense monitoring of the source with Swift and find an optical photometric period of 63.9 +/- 0.5 d and an X-ray flux period of 66.8 +/- 0.4 d. The optical period is consistent with the orbital period, while the X-ray flux period is significantly longer. We discuss possible reasons for this discrepancy, which could be due to a super-orbital period caused by a precessing accretion disk or an orbital resonance. We put the orbital period of P13 into context with the orbital periods implied for two other ultra-luminous pulsars, M82 X-2 and NGC 5907 ULX and discuss possible implications for the system parameters.

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.

Based on phase-resolved broadband spectroscopy using $XMM$-$Newton$ and $NuSTAR$, we report on a potential cyclotron resonant scattering feature at $E \sim 13$ keV in the pulsed spectrum of the recently discoverd ULX pulsar NGC 300 ULX1. If this interpretation is correct, the implied magnetic field of the central neutron star is $B \sim 10^{12}$ G (assuming scattering off electrons), similar to that estimated from the observed spin-up of the star, and also similar to known Galactic X-ray pulsars. We discuss the implications of this result for the connection between NGC 300 ULX1 and the other known ULX pulsars, particularly in light of the recent discovery of a likely proton Cyclotron line in another ULX, M51 ULX-8.

We present broadband X-ray analyses of a sample of bright ultraluminous X-ray sources with the goal of investigating the spectral similarity of this population to the known ULX pulsars, M82 X-2, NGC7793 P13 and NGC5907 ULX. We perform a phase-resolved analysis of the broadband XMM-Newton+NuSTAR dataset of NGC5907 ULX, finding that the pulsed emission from the accretion column in this source exhibits a similar spectral shape to that seen in both M82 X-2 and NGC7793 P13, and that this is responsible for the excess emission observed at the highest energies when the spectra are fit with accretion disk models. We then demonstrate that similar 'hard' excesses are seen in all the ULXs in the broadband sample. Finally, for the ULXs where the nature of the accretor is currently unknown, we test whether the hard excesses are all consistent with being produced by an accretion column similar to those present in M82 X-2, NGC7793 P13 and NGC5907 ULX. Based on the average shape of the pulsed emission, we find that in all cases a similar accretion column can successfully reproduce the observed data, consistent with the hypothesis that this ULX sample may be dominated by neutron star accretors. Compared to the known pulsar ULXs, our spectral fits for the remaining ULXs suggest that the non-pulsed emission from the accretion flow beyond the magnetosphere makes a stronger relative contribution than the component associated with the accretion column. If these sources do also contain neutron star accretors, this may help to explain the lack of detected pulsations.

Ultraluminous X-ray sources (ULXs) are bright X-ray sources in nearby galaxies not associated with the central supermassive black hole. Their luminosities imply they are powered by either an extreme accretion rate onto a compact stellar remnant, or an intermediate mass ($\sim100-10^5$ M$_{\odot}$) black hole. The recent detection of coherent pulsations coming from three bright ULXs demonstrates that some of these sources are powered by accretion onto a neutron star, implying accretion rates significantly in excess of the Eddington limit, a high degree of geometric beaming, or both. The physical challenges associated with the high implied accretion rates can be mitigated if the neutron star surface field is high - in the magnetar regime ($10^{14}$ G), since this suppresses the electron scattering cross section, reducing the radiation pressure that chokes off accretion for high luminosities. One of the few ways to determine surface magnetic fields is through the detection of cyclotron resonance scattering features (CRSFs) produced by the transition of charged particles between quantized Landau levels. To date, CRSFs have only been detected in Galactic accreting pulsars. Here we present the detection at 3.8-{\sigma} significance of a strong absorption line at a rest-frame energy of 4.5 keV in the Chandra spectrum of a ULX in M51. We find that this feature is likely to be a CRSF produced by the strong magnetic field of a neutron star. Assuming scattering off electrons, the magnetic field strength is implied to be $\sim10^{11}$ G, however the line is narrower than any electron CRSFs previously observed, and assuming thermal broadening, the implied temperature is significantly cooler than the accretion disk or column. The line shape is, however, consistent with a proton resonance scattering feature, implying that the neutron star has a magnetic field near the surface of B$\sim10^{15}$ G.

We studied the nearby, interacting galaxy NGC 5195 (M51b) in the radio, optical and X-ray bands. We mapped the extended, low-surface-brightness features of its radio-continuum emission; determined the energy content of its complex structure of shock-ionized gas; constrained the current activity level of its supermassive nuclear black hole. In particular, we combined data from the European Very Long Baseline Interferometry Network (~1-pc scale), from our new e-MERLIN observations (~10-pc scale), and from the Very Large Array (~100-1000-pc scale), to obtain a global picture of energy injection in this galaxy. We put an upper limit to the luminosity of the (undetected) flat-spectrum radio core. We find steep-spectrum, extended emission within 10 pc of the nuclear position, consistent with optically-thin synchrotron emission from nuclear star formation or from an outflow powered by an active galactic nucleus (AGN). A linear spur of radio emission juts out of the nuclear source towards the kpc-scale arcs (detected in radio, Halpha and X-ray bands). From the size, shock velocity, and Balmer line luminosity of the kpc-scale bubble, we estimate that it was inflated by a long-term-average mechanical power ~3-6 x 10^{41} erg/s over the last 3-6 Myr. This is an order of magnitude more power than can be provided by the current level of star formation, and by the current accretion power of the supermassive black hole. We argue that a jet-inflated bubble scenario associated with previous episodes of AGN activity is the most likely explanation for the kpc-scale structures.

The basic unified model of active galactic nuclei (AGN) invokes an anisotropic obscuring structure, usually referred to as a torus, to explain AGN obscuration as an angle-dependent effect. We present a new grid of X-ray spectral templates based on radiative transfer calculations in neutral gas in an approximately toroidal geometry, appropriate for CCD-resolution X-ray spectra (FWHM > 130 eV). Fitting the templates to broadband X-ray spectra of AGN provides constraints on two important geometrical parameters of the gas distribution around the supermassive black hole: the average column density and the covering factor. Compared to the currently available spectral templates, our model is more flexible, and capable of providing constraints on the main torus parameters in a wider range of AGN. We demonstrate the application of this model using hard X-ray spectra from NuSTAR (3-79 keV) for four AGN covering a variety of classifications: 3C 390.3, NGC 2110, IC 5063, and NGC 7582. This small set of examples was chosen to illustrate the range of possible torus configurations, from disk-like to sphere-like geometries with column densities below, as well as above, the Compton-thick threshold. This diversity of torus properties challenges the simple assumption of a standard geometrically and optically thick toroidal structure commonly invoked in the basic form of the unified model of AGN. Finding broad consistency between the our constraints and those from infrared modeling, we discuss how the approach from the X-ray band complements similar measurements of AGN structures at other wavelengths.

Hot Dust-Obscured Galaxies (Hot DOGs) are hyperluminous ($L_{\mathrm{8-1000\,\mu m}}>10^{13}\,\mathrm{L_\odot}$) infrared galaxies with extremely high (up to hundreds of K) dust temperatures. The sources powering both their extremely high luminosities and dust temperatures are thought to be deeply buried and rapidly accreting supermassive black holes (SMBHs). Hot DOGs could therefore represent a key evolutionary phase in which the SMBH growth peaks. X-ray observations can be used to study their obscuration levels and luminosities. In this work, we present the X-ray properties of the 20 most-luminous ($L_{\mathrm{bol}}\gtrsim10^{14}\, L_\odot$) known Hot DOGs at $z=2-4.6$. Five of them are covered by long-exposure ($10-70$ ks) Chandra and XMM-Newton observations, with three being X-ray detected, and we study their individual properties. One of these sources (W0116$-$0505) is a Compton-thick candidate, with column density $N_H=(1.0-1.5)\times10^{24}\,\mathrm{cm^{-2}}$ derived from X-ray spectral fitting. The remaining 15 Hot DOGs have been targeted by a Chandra snapshot (3.1 ks) survey. None of these 15 is individually detected; therefore we applied a stacking analysis to investigate their average emission. From hardness-ratio analysis, we constrained the average obscuring column density and intrinsic luminosity to be log$N_H\,\mathrm{[cm^{-2}]}>23.5$ and $L_X\gtrsim10^{44}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, which are consistent with results for individually detected sources. We also investigated the $L_X-L_{6\mu\mathrm{m}}$ and $L_X-L_{bol}$ relations, finding hints that Hot DOGs are typically X-ray weaker than expected, although larger samples of luminous obscured QSOs are needed to derive solid conclusions.

The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in the disc (implying surface dipole fields that are sub-magnetar in strength), then the large scale-height flow can precess through the Lense-Thirring effect which could provide an explanation for the observed super-orbital periods. By connecting the details of the Lense-Thirring effect with the observed pulsar spin period, we are able to infer the moment-of-inertia and therefore equation-of-state of the neutron star without relying on the inclination of, or distance to the system. We apply our technique to the case of NGC 7793 P13 and demonstrate that stronger magnetic fields imply stiffer equations of state. We discuss the caveats and uncertainties, many of which can be addressed through forthcoming radiative magnetohydrodynamic (RMHD) simulations and their connection to observation.

We present results from a coordinated $XMM$-$Newton$+$NuSTAR$ observation of the type 1.8 Seyfert galaxy IRAS 13197-1627. This is a highly complex source, with strong contributions from relativistic reflection from the inner accretion disk, neutral absorption and further reprocessing by more distant material, and ionised absorption from an outflow. We undertake a detailed spectral analysis combining the broadband coverage provided by $XMM$-$Newton$+$NuSTAR$ with a multi-epoch approach incorporating archival observations performed by $XMM$-$Newton$ and $Suzaku$. Our focus is on characterising the reflection from the inner accretion disk, which previous works have suggested may dominate the AGN emission, and constraining the black hole spin. Using lamppost disk reflection models, we find that the results for the inner disk are largely insensitive to assumptions regarding the geometry of the distant reprocessor and the precise form of the illuminating X-ray continuum. However, these results do depend on the treatment of the iron abundance of the distant absorber/reprocessor. The multi-epoch data favour a scenario in which the AGN is chemically homogeneous, and we find that a rapidly rotating black hole is preferred, with $a^* \geq 0.7$, but a slowly-rotating black hole is not strongly excluded. In addition to the results for the inner disk, we also find that both the neutral and ionised absorbers vary from epoch to epoch, implying that both have some degree of inhomogeneity in their structure.

We present a detailed, broadband X-ray spectral analysis of the ULX pulsar NGC 7793 P13, a known super-Eddington source, utilizing data from the $XMM$-$Newton$, $NuSTAR$ and $Chandra$ observatories. The broadband $XMM$-$Newton+NuSTAR$ spectrum of P13 is qualitatively similar to the rest of the ULX sample with broadband coverage, suggesting that additional ULXs in the known population may host neutron star accretors. Through time-averaged, phase-resolved and multi-epoch studies, we find that two non-pulsed thermal blackbody components with temperatures $\sim$0.5 and $\sim$1.5 keV are required to fit the data below 10 keV, in addition to a third continuum component which extends to higher energies and is associated with the pulsed emission from the accretion column. The characteristic radii of the thermal components appear to be similar, and are too large to be associated with the neutron star itself, so the need for two components likely indicates the accretion flow outside the magnetosphere is complex. We suggest a scenario in which the thick inner disc expected for super-Eddington accretion begins to form, but is terminated by the neutron star's magnetic field soon after its onset, implying a limit of $B \lesssim 6 \times 10^{12}$ G for the dipolar component of the central neutron star's magnetic field. Evidence of similar termination of the disc in other sources may offer a further means of identifying additional neutron star ULXs. Finally, we examine the spectrum exhibited by P13 during one of its unusual 'off' states. These data require both a hard powerlaw component, suggesting residual accretion onto the neutron star, and emission from a thermal plasma, which we argue is likely associated with the P13 system.

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 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 present analysis of NuSTAR X-ray observations in the 3-79 keV energy band of the Seyfert 2 galaxy NGC 4388, taken in 2013. The broadband sensitivity of NuSTAR, covering the Fe K$\alpha$ line and Compton reflection hump, enables tight constraints to be placed on reflection features in AGN X-ray spectra, thereby providing insight into the geometry of the circumnuclear material. In this observation, we found the X-ray spectrum of NGC 4388 to be well described by a moderately absorbed power law with non-relativistic reflection. We fit the spectrum with phenomenological reflection models and a physical torus model, and find the source to be absorbed by Compton-thin material (N$_{H} = (6.5\pm0.8)\times10^{23}$ cm$^{-2}$) with a very weak Compton reflection hump (R $<$ 0.09) and an exceptionally large Fe K$\alpha$ line (EW $= 368^{+56}_{-53}$ eV) for a source with weak or no reflection. Calculations using a thin-shell approximation for the expected Fe K$\alpha$ EW indicate that an Fe K$\alpha$ line originating from Compton-thin material presents a possible explanation.

We present first results from a series of $NuSTAR$ observations of the black hole X-ray binary V404 Cyg obtained during its summer 2015 outburst, primarily focusing on observations during the height of this outburst activity. The $NuSTAR$ data show extreme variability in both the flux and spectral properties of the source. This is partly driven by strong and variable line-of-sight absorption, similar to previous outbursts. The latter stages of this observation are dominated by strong flares, reaching luminosities close to Eddington. During these flares, the central source appears to be relatively unobscured and the data show clear evidence for a strong contribution from relativistic reflection, providing a means to probe the geometry of the innermost accretion flow. Based on the flare properties, analogy with other Galactic black hole binaries, and also the simultaneous onset of radio activity, we argue that this intense X-ray flaring is related to transient jet activity during which the ejected plasma is the primary source of illumination for the accretion disk. If this is the case, then our reflection modelling implies that these jets are launched in close proximity to the black hole (as close as a few gravitational radii), consistent with expectations for jet launching models that tap either the spin of the central black hole, or the very innermost accretion disk. Our analysis also allows us to place the first constraints on the black hole spin for this source, which we find to be $a^* > 0.92$ (99% statistical uncertainty, based on an idealized lamppost geometry).

NGC 1448 is one of the nearest luminous galaxies ($L_{8-1000\mu m} >$ 10$^{9} L_{\odot}$) to ours ($z$ $=$ 0.00390), and yet the active galactic nucleus (AGN) it hosts was only recently discovered, in 2009. In this paper, we present an analysis of the nuclear source across three wavebands: mid-infrared (MIR) continuum, optical, and X-rays. We observed the source with the Nuclear Spectroscopic Telescope Array (NuSTAR), and combined this data with archival Chandra data to perform broadband X-ray spectral fitting ($\approx$0.5-40 keV) of the AGN for the first time. Our X-ray spectral analysis reveals that the AGN is buried under a Compton-thick (CT) column of obscuring gas along our line-of-sight, with a column density of $N_{\rm H}$(los) $\gtrsim$ 2.5 $\times$ 10$^{24}$ cm$^{-2}$. The best-fitting torus models measured an intrinsic 2-10 keV luminosity of $L_{2-10\rm{,int}}$ $=$ (3.5-7.6) $\times$ 10$^{40}$ erg s$^{-1}$, making NGC 1448 one of the lowest luminosity CTAGNs known. In addition to the NuSTAR observation, we also performed optical spectroscopy for the nucleus in this edge-on galaxy using the European Southern Observatory New Technology Telescope. We re-classify the optical nuclear spectrum as a Seyfert on the basis of the Baldwin-Philips-Terlevich diagnostic diagrams, thus identifying the AGN at optical wavelengths for the first time. We also present high spatial resolution MIR observations of NGC 1448 with Gemini/T-ReCS, in which a compact nucleus is clearly detected. The absorption-corrected 2-10 keV luminosity measured from our X-ray spectral analysis agrees with that predicted from the optical [OIII]$\lambda$5007\AA\ emission line and the MIR 12$\mu$m continuum, further supporting the CT nature of the AGN.

We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of <z> = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy (<10 keV) X-ray counterparts from XMM-Newton, Chandra, and Swift/XRT. The mid-infrared (MIR) analysis, using WISE all-sky survey data, shows that MIR AGN color selections miss a large fraction of the NuSTAR-selected AGN population, from ~15% at the highest luminosities (Lx > 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by <10 keV X-ray missions.

Hot, Dust-Obscured Galaxies (Hot DOGs), selected from the WISE all sky infrared survey, host some of the most powerful Active Galactic Nuclei (AGN) known, and might represent an important stage in the evolution of galaxies. Most known Hot DOGs are at $z> 1.5$, due in part to a strong bias against identifying them at lower redshift related to the selection criteria. We present a new selection method that identifies 153 Hot DOG candidates at $z\sim 1$, where they are significantly brighter and easier to study. We validate this approach by measuring a redshift $z=1.009$, and an SED similar to higher redshift Hot DOGs for one of these objects, WISE J1036+0449 ($L_{\rm\,Bol}\simeq 8\times 10^{46}\rm\,erg\,s^{-1}$), using data from Keck/LRIS and NIRSPEC, SDSS, and CSO. We find evidence of a broadened component in MgII, which, if due to the gravitational potential of the supermassive black hole, would imply a black hole mass of $M_{\rm\,BH}\simeq 2 \times 10^8 M_{\odot}$, and an Eddington ratio of $\lambda_{\rm\,Edd}\simeq 2.7$. WISE J1036+0449 is the first Hot DOG detected by NuSTAR, and the observations show that the source is heavily obscured, with a column density of $N_{\rm\,H}\simeq(2-15)\times10^{23}\rm\,cm^{-2}$. The source has an intrinsic 2-10 keV luminosity of $\sim 6\times 10^{44}\rm\,erg\,s^{-1}$, a value significantly lower than that expected from the mid-infrared/X-ray correlation. We also find that the other Hot DOGs observed by X-ray facilities show a similar deficiency of X-ray flux. We discuss the origin of the X-ray weakness and the absorption properties of Hot DOGs. Hot DOGs at $z\lesssim1$ could be excellent laboratories to probe the characteristics of the accretion flow and of the X-ray emitting plasma at extreme values of the Eddington ratio.

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.

We report the detection of coherent pulsations from the ultraluminous X-ray source NGC 7793 P13. The ~0.42s nearly sinusoidal pulsations were initially discovered in broadband X-ray observations using XMM-Newton and NuSTAR taken in 2016. We subsequently also found pulsations in archival XMM-Newton data taken in 2013 and 2014. The significant (>>5 sigma) detection of coherent pulsations demonstrates that the compact object in P13 is a neutron star with an observed peak luminosity of ~1e40 erg/s (assuming isotropy), well above the Eddington limit for a 1.4 M_sun accretor. This makes P13 the second ultraluminous X-ray source known to be powered by an accreting neutron star. The pulse period varies between epochs, with a slow but persistent spin up over the 2013-2016 period. This spin-up indicates a magnetic field of B ~ 1.5e12 G, typical of many accreting pulsars. The most likely explanation for the extreme luminosity is a high degree of beaming, however this is difficult to reconcile with the sinusoidal pulse profile.

We present broad-band, multi-epoch X-ray spectroscopy of the pulsating ultra-luminous X-ray source (ULX) in NGC 5907. Simultaneous XMM-Newton and NuSTAR data from 2014 are best described by a multi-color black-body model with a temperature gradient as a function of accretion disk radius significantly flatter than expected for a standard thin accretion disk (T(r) ~ r^{-p}, with p=0.608^{+0.014}_{-0.012}). Additionally, we detect a hard power-law tail at energies above 10 keV, which we interpret as being due to Comptonization. We compare this observation to archival XMM-Newton, Chandra, and NuSTAR data from 2003, 2012, and 2013, and investigate possible spectral changes as a function of phase over the 78d super-orbital period of this source. We find that observations taken around phases 0.3-0.4 show very similar temperature profiles, even though the observed flux varies significantly, while one observation taken around phase 0 has a significantly steeper profile. We discuss these findings in light of the recent discovery that the compact object is a neutron star and show that precession of the accretion disk or the neutron star can self-consistently explain most observed phenomena.

We present the 3-8 keV and 8-24 keV number counts of active galactic nuclei (AGN) identified in the NuSTAR extragalactic surveys. NuSTAR has now resolved 33-39% of the X-ray background in the 8-24 keV band, directly identifying AGN with obscuring columns up to approximately 1e25 /cm2. In the softer 3-8 keV band the number counts are in general agreement with those measured by XMM-Newton and Chandra over the flux range 5e-15 < S(3 - 8 keV)/(erg/cm2/s) < 1e-12 probed by NuSTAR. In the hard 8-24 keV band NuSTAR probes fluxes over the range 2e-14 < S(8-24 keV)/(erg/cm2/s) < 1e-12, a factor of approximately 100 fainter than previous measurements. The 8-24 keV number counts match predictions from AGN population synthesis models, directly confirming the existence of a population of obscured and/or hard X-ray sources inferred from the shape of the integrated cosmic X-ray background. The measured NuSTAR counts lie significantly above simple extrapolation with a Euclidian slope to low flux of the Swift/BAT 15-55 keV number counts measured at higher fluxes S(15-55 keV) > 1e-11 erg/cm2/s, reflecting the evolution of the AGN population between the Swift/BAT local (z<0.1) sample and NuSTAR's z~1 sample. CXB synthesis models, which account for AGN evolution, lie above the Swift/BAT measurements, suggesting that they do not fully capture the evolution of obscured AGN at low redshifts.

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.

We report the detection of a $78.1\pm0.5$ day period in the X-ray lightcurve of the extreme ultraluminous X-ray source NGC 5907 ULX1 ($L_{\rm{X,peak}}\sim5\times10^{40}$ erg s$^{-1}$), discovered during an extensive monitoring program with Swift. These periodic variations are strong, with the observed flux changing by a factor of $\sim$3-4 between the peaks and the troughs of the cycle; our simulations suggest that the observed periodicity is detected comfortably in excess of 3$\sigma$ significance. We discuss possible origins for this X-ray period, but conclude that at the current time we cannot robustly distinguish between orbital and super-orbital variations.

We present an analysis of NuSTAR observations of a hard intermediate state of the transient black hole GX 339-4 taken in January 2015. As the source softened significantly over the course of the 1.3 d-long observation we split the data into 21 sub-sets and find that the spectrum of all of them can be well described by a power-law continuum with an additional relativistically blurred reflection component. The photon index increases from ~1.69 to ~1.77 over the course of the observation. The accretion disk is truncated at around 9 gravitational radii in all spectra. We also perform timing analysis on the same 21 individual data sets, and find a strong type-C quasi-periodic oscillation (QPO), which increase in frequency from ~0.68 to ~1.05 Hz with time. The frequency change is well correlated with the softening of the spectrum. We discuss possible scenarios for the production of the QPO and calculate predicted inner radii in the relativistic precession model as well as the global disk mode oscillations model. We find discrepancies with respect to the observed values in both models unless we allow for a black hole mass of ~100 M_sun , which is highly unlikely. We discuss possible systematic uncertainties, in particular with the measurement of the inner accretion disk radius in the relativistic reflection model. We conclude that the combination of observed QPO frequencies and inner accretion disk radii, as obtained from spectral fitting, is difficult to reconcile with current models.

M82 X-1 is one of the brightest ultraluminous X-ray sources (ULXs) known, which, assuming Eddington-limited accretion and other considerations, makes it one of the best intermediate-mass black hole (IMBH) candidates. However, the ULX may still be explained by super-Eddington accretion onto a stellar-remnant black hole. We present simultaneous NuSTAR, Chandra and Swift/XRT observations during the peak of a flaring episode with the aim of modeling the emission of M82 X-1 and yielding insights into its nature. We find that thin-accretion disk models all require accretion rates at or above the Eddington limit in order to reproduce the spectral shape, given a range of black hole masses and spins. Since at these high Eddington ratios the thin-disk model breaks down due to radial advection in the disk, we discard the results of the thin-disk models as unphysical. We find that the temperature profile as a function of disk radius ($T(r)\propto r^{-p}$) is significantly flatter ($p=0.55^{+ 0.07}_{- 0.04}$) than expected for a standard thin disk ($p=0.75$). A flatter profile is instead characteristic of a slim disk which is highly suggestive of super-Eddington accretion. Furthermore, radiation hydrodynamical simulations of super-Eddington accretion have shown that the predicted spectra of these systems are very similar to what we observe for M82 X-1. We therefore conclude that M82 X-1 is a super-Eddington accretor. Our mass estimates inferred from the inner disk radius imply a stellar-remnant black hole ($M_{\rm BH}=26^{+9}_{-6} M_{\odot}$) when assuming zero spin, or an IMBH ($M_{\rm BH}=125^{+45}_{-30} M_{\odot}$) when assuming maximal spin.

We present the detection of an absorpton feature at $E=8.77^{+0.05}_{-0.06}$ keV in the combined X-ray spectrum of the ultraluminous X-ray source NGC 1313 X-1 observed with XMM-Newton and NuSTAR, significant at the 3$\sigma$ level. If associated with blueshifted ionized iron, the implied outflow velocity is ~0.2$c$ for Fe XXVI, or ~0.25$c$ for Fe XXV. These velocities are similar to the ultrafast outflow seen in absorption recently discovered in this source at lower energies by XMM-Newton, and we therefore conclude that this is an iron component to the same outflow. Photoionization modeling marginally prefers the Fe XXV solution, but in either case the outflow properties appear to be extreme, potentially supporting a super-Eddington hypothesis for NGC 1313 X-1.

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.

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.

In this paper we report the results of an X-ray monitoring campaign on the heavily obscured Seyfert galaxy Markarian 3 carried out between the fall of 2014 and the spring of 2015 with NuSTAR, Suzaku and XMM-Newton. The hard X-ray spectrum of Markarian 3 is variable on all the time scales probed by our campaign, down to a few days. The observed continuum variability is due to an intrinsically variable primary continuum seen in transmission through a large, but still Compton-thin column density (N_H~0.8-1.1$\times$10$^{24}$ cm$^{-2}$). If arranged in a spherical-toroidal geometry, the Compton scattering matter has an opening angle ~66 degrees and is seen at a grazing angle through its upper rim (inclination angle ~70 degrees). We report a possible occultation event during the 2014 campaign. If the torus is constituted by a system of clouds sharing the same column density, this event allows us to constrain their number (17$\pm$5) and individual column density, [~(4.9$\pm$1.5)$\times$10$^{22}$ cm$^{-2}$]. The comparison of IR and X-ray spectroscopic results with state-of-the art "torus" models suggests that at least two thirds of the X-ray obscuring gas volume might be located within the dust sublimation radius. We report also the discovery of an ionized absorber, characterised by variable resonant absorption lines due to He- and H-like iron. This discovery lends support to the idea that moderate column density absorbers could be due to clouds evaporated at the outer surface of the torus, possibly accelerated by the radiation pressure due to the central AGN emission leaking through the patchy absorber.

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.

Hot Dust-Obscured Galaxies (Hot DOGs) are a population of hyper-luminous infrared galaxies identified by the WISE mission from their very red mid-IR colors, and characterized by hot dust temperatures ($T>60~\rm K$). Several studies have shown clear evidence that the IR emission in these objects is powered by a highly dust-obscured AGN that shows close to Compton-thick absorption at X-ray wavelengths. Thanks to the high AGN obscuration, the host galaxy is easily observable, and has UV/optical colors usually consistent with those of a normal galaxy. Here we discuss a sub-population of 8 Hot DOGs that show enhanced rest-frame UV/optical emission. We discuss three scenarios that might explain the excess UV emission: (i) unobscured light leaked from the AGN by reflection over the dust or by partial coverage of the accretion disk; (ii) a second unobscured AGN in the system; or (iii) a luminous young starburst. X-ray observations can help discriminate between these scenarios. We study in detail the blue excess Hot DOG WISE J020446.13-050640.8, which was serendipitously observed by Chandra/ACIS-I for 174.5 ks. The X-ray spectrum is consistent with a single, hyper-luminous, highly absorbed AGN, and is strongly inconsistent with the presence of a secondary unobscured AGN. Based on this, we argue that the excess blue emission in this object is most likely either due to reflection or a co-eval starburst. We favor the reflection scenario as the unobscured star-formation rate needed to power the UV/optical emission would be $\gtrsim 1000~\rm M_{\odot}~\rm yr^{-1}$. Deep polarimetry observations could confirm the reflection hypothesis.

We present simultaneous XMM-Newton and NuSTAR observations spanning 3-78 keV of the nearest radio galaxy, Centaurus A (Cen A). The accretion geometry around the central engine in Cen A is still debated, and we investigate possible configurations using detailed X-ray spectral modeling. NuSTAR imaged the central region of Cen A with sub-arcminute resolution at X-ray energies above 10 keV for the first time, but finds no evidence for an extended source or other off-nuclear point-sources. The XMM-Newton and NuSTAR spectra agree well and can be described with an absorbed power-law with a photon index {\Gamma} = 1.815 +/- 0.005 and a fluorescent Fe K{\alpha} line in good agreement with literature values. The spectrum does not require a high-energy exponential rollover, with a constraint of E_fold > 1 MeV. A thermal Comptonization continuum describes the data well, with parameters that agree with values measured by INTEGRAL, in particular an electron temperature kTe between ~100-300 keV, seed photon input temperatures between 5-50 eV. We do not find evidence for reflection or a broad iron line and put stringent upper limits of R < 0.01 on the reflection fraction and accretion disk illumination. We use archival Chandra data to estimate the contribution from diffuse emission, extra-nuclear point-sources, and the outer X-ray jet to the observed NuSTAR and XMM-Newton X-ray spectra and find the contribution to be negligible. We discuss different scenarios for the physical origin of the observed hard X-ray spectrum, and conclude that the inner disk is replaced by an advection-dominated accretion flow or that the X-rays are dominated by synchrotron self-Compton emission from the inner regions of the radio jet or a combination thereof.

An enigmatic group of objects, unabsorbed Seyfert 2s may have intrinsically weak broad line regions, obscuration in the line of sight to the BLR but not to the X-ray corona, or so much obscuration that the X-ray continuum is completely suppressed and the observed spectrum is actually scattered into the line of sight from nearby material. NGC 3660 has been shown to have weak broad optical/near infrared lines, no obscuration in the soft X-ray band, and no indication of "changing look" behavior. The only previous hard X-ray detection of this source by Beppo-SAX seemed to indicate that the source might harbor a heavily obscured nucleus. However, our analysis of a long-look Suzaku observation of this source shows that this is not the case, and that this source has a typical power law X-ray continuum with normal reflection and no obscuration. We conclude that NGC 3660 is confirmed to have no unidentified obscuration and that the anomolously high Beppo-SAX measurement must be due to source confusion or similar, being inconsistent with our Suzaku measurements as well as non-detections from Swift-BAT and RXTE.

The recent discovery by Bachetti et al. (2014) of a pulsar in M82 that can reach luminosities of up to 10^40 ergs s^-1, a factor of ~100 the Eddington luminosity for a 1.4 Msol compact object, poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in the light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.5-8 keV X-ray emission from this source from 15 years of Chandra observations. We fit the Chandra spectra of the pulsar with a power-law model and a disk black body model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where 4/19 observations have a pile-up fraction >10%, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (L_X>10^39 ergs s^-1) is Gamma=1.33+/-0.15. For the disk black body model, the average temperature is T=3.24+/-0.65 keV, consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cut-off, where Gamma=0.6+/-0.3 and E_C=14^{+5}_{-3} keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10^39 ergs s^-1, greater than 10 times its Eddington limit.

We present the first direct measurements of the rest-frame 10-40 keV X-ray luminosity function (XLF) of Active Galactic Nuclei (AGNs) based on a sample of 94 sources at 0.1 < z <3, selected at 8-24 keV energies from sources in the NuSTAR extragalactic survey program. Our results are consistent with the strong evolution of the AGN population seen in prior, lower-energy studies of the XLF. However, different models of the intrinsic distribution of absorption, which are used to correct for selection biases, give significantly different predictions for the total number of sources in our sample, leading to small, systematic differences in our binned estimates of the XLF. Adopting a model with a lower intrinsic fraction of Compton-thick sources and a larger population of sources with column densities N_H ~ 10^{23-24} /cm2 or a model with stronger Compton reflection component (with a relative normalization of R ~ 2 at all luminosities) can bring extrapolations of the XLF from 2-10 keV into agreement with our NuSTAR sample. Ultimately, X-ray spectral analysis of the NuSTAR sources is required to break this degeneracy between the distribution of absorbing column densities and the strength of the Compton reflection component and thus refine our measurements of the XLF. Furthermore, the models that successfully describe the high-redshift population seen by NuSTAR tend to over-predict previous, high-energy measurements of the local XLF, indicating that there is evolution of the AGN population that is not fully captured by the current models.

We present initial results and the source catalog from the NuSTAR survey of the Extended Chandra Deep Field South (hereafter, ECDFS) - currently the deepest contiguous component of the NuSTAR extragalactic survey program. The survey covers the full ~30 arcmin x 30 arcmin area of this field to a maximum depth of ~360 ks (~220 ks when corrected for vignetting at 3-24 keV), reaching sensitivity limits of ~1.3 x 10^-14 erg/cm2/s (3-8 keV), ~3.4 x 10^-14 erg/cm2/s (8-24 keV) and ~3.0 x 10^-14 erg/cm2/s (3-24 keV). Fifty four (54) sources are detected over the full field, although five of these are found to lie below our significance threshold once contaminating flux from neighboring (i.e., blended) sources is taken into account. Of the remaining 49 that are significant, 19 are detected in the 8-24 keV band. The 8-24 keV to 3-8 keV band ratios of the twelve sources that are detected in both bands span the range 0.39-1.7, corresponding to a photon index range of Gamma ~ 0.5-2.3, with a median photon index of 1.70 +/- 0.52. The redshifts of the 49 sources in our main sample span the range z = 0.21-2.7, and their rest-frame 10-40 keV luminosities (derived from the observed 8-24 keV fluxes) span the range L(10-40 keV) ~ (0.7-300) x 10^43 erg/s, sampling below the "knee" of the X-ray luminosity function out to z ~ 0.8-1. Finally, we identify one NuSTAR source that has neither a Chandra nor an XMM-Newton counterpart, but that shows evidence of nuclear activity at infrared wavelengths, and thus may represent a genuine, new X-ray source detected by NuSTAR in the ECDFS.

NGC 7582 is a well-studied X-ray bright Seyfert 2 with moderately heavy ($N_{\text{H}}\sim10^{23}-10^{24}$~cm$^{-2}$), highly variable absorption and strong reflection spectral features. The spectral shape changed around the year 2000, dropping in observed flux and becoming much more highly absorbed. Two scenarios have been put forth to explain this spectral change: 1) the central X-ray source partially ``shut off'' around this time, decreasing in intrinsic luminosity, with a delayed decrease in reflection features due to the light-crossing time of the Compton-thick material or 2) the source became more heavily obscured, with only a portion of the power law continuum leaking through. NuSTAR observed NGC~7582 twice in 2012, two weeks apart, in order to quantify the reflection using high-quality data above 10 keV. We find that the most plausible scenario is that NGC 7582 has recently become more heavily absorbed by a patchy torus with a covering fraction of $\sim\,80-90\%$ and an equatorial column density of $\sim 3 \times10^{24}$ cm$^{-2}$. We find the need for an additional highly variable full-covering absorber with $N_{\text{H}}= 4-6 \times10^{23}$ cm$^{-2}$ in the line of sight, possibly associated with a hidden broad line region.

We present a broad-band (~0.3-70 keV) spectral and temporal analysis of NuSTAR observations of the luminous infrared galaxy NGC 6240, combined with archival Chandra, XMM-Newton and BeppoSAX data. NGC 6240 is a galaxy in a relatively early merger state with two distinct nuclei separated by ~1."5. Previous Chandra observations have resolved the two nuclei, showing that they are both active and obscured by Compton-thick material. Although they cannot be resolved by NuSTAR, thanks to the unprecedented quality of the NuSTAR data at energies >10 keV, we clearly detect, for the first time, both the primary and the reflection continuum components. The NuSTAR hard X-ray spectrum is dominated by the primary continuum piercing through an absorbing column density which is mildly optically thick to Compton scattering (tau ~ 1.2, N_H ~ 1.5 x 10^(24) cm^-2). We detect moderate hard X-ray (> 10 keV) flux variability up to 20% on short (15-20 ksec) timescales. The amplitude of the variability is maximum at ~30 keV and is likely to originate from the primary continuum of the southern nucleus. Nevertheless, the mean hard X-ray flux on longer timescales (years) is relatively constant. Moreover, the two nuclei remain Compton-thick, although we find evidence of variability of the material along the line of sight with column densities N_H <~ 2 x 10^(23) cm-2 over long (~3-15 years) timescales. The observed X-ray emission in the NuSTAR energy range is fully consistent with the sum of the best-fit models of the spatially resolved Chandra spectra of the two nuclei.

We examine the host morphologies of heavily obscured active galactic nuclei (AGN) at $z\sim1$ to test whether obscured supermassive black hole growth at this epoch is preferentially linked to galaxy mergers. Our sample consists of 154 obscured AGN with $N_{\rm H}>10^{23.5}$ cm$^{-2}$ and $z<1.5$. Using visual classifications, we compare the morphologies of these AGN to control samples of moderately obscured ($10^{22}$ cm$^{-2}$ $<N_{\rm H}< 10^{23.5}$ cm$^{-2}$) and unobscured ($N_{\rm H}<10^{22}$ cm$^{-2}$) AGN. These control AGN are matched in redshift and intrinsic X-ray luminosity to our heavily obscured AGN. We find that heavily obscured AGN at z~1 are twice as likely to be hosted by late-type galaxies relative to unobscured AGN ($65.3^{+4.1}_{-4.6}\%$ vs $34.5^{+2.9}_{-2.7}\%$) and three times as likely to exhibit merger or interaction signatures ($21.5^{+4.2}_{-3.3}\%$ vs $7.8^{+1.9}_{-1.3}\%$). The increased merger fraction is significant at the 3.8$\sigma$ level. We also find that the incidence of point-like morphologies is inversely proportional to obscuration. If we exclude all point sources and consider only extended hosts, we find the correlation between merger fraction and obscuration is still evident, however at a reduced statistical significance ($2.5\sigma$). The fact that we observe a different disk/spheroid fraction versus obscuration indicates that viewing angle cannot be the only thing differentiating our three AGN samples, as a simple unification model would suggest. The increased fraction of disturbed morphologies with obscuration supports an evolutionary scenario, in which Compton-thick AGN are a distinct phase of obscured SMBH growth following a merger/interaction event. Our findings also suggest that some of the merger-triggered SMBH growth predicted by recent AGN fueling models may be hidden among the heavily obscured, Compton-thick population.

We present two NuSTAR observations of the local Seyfert 2 active galactic nucleus (AGN) and an ultraluminous X-ray source (ULX) candidate in NGC 5643. Together with archival data from Chandra, XMM-Newton and Swift-BAT, we perform a high-quality broadband spectral analysis of the AGN over two decades in energy ($\sim$0.5-100 keV). Previous X-ray observations suggested that the AGN is obscured by a Compton-thick (CT) column of obscuring gas along our line-of-sight. However, the lack of high-quality $\gtrsim$ 10 keV observations, together with the presence of a nearby X-ray luminous source, NGC 5643 X-1, had left significant uncertainties in the characterization of the nuclear spectrum. NuSTAR now enables the AGN and NGC 5643 X-1 to be separately resolved above 10 keV for the first time and allows a direct measurement of the absorbing column density toward the nucleus. The new data show that the nucleus is indeed obscured by a CT column of $N_{\rm{H}}$ $\gtrsim$ 5 $\times$ 10$^{24}$ cm$^{-2}$. The range of 2-10 keV absorption-corrected luminosity inferred from the best fitting models is $L_{2-10,\rm{int}} =$ (0.8-1.7) $\times$ 10$^{42}$ erg s$^{-1}$, consistent with that predicted from multiwavelength intrinsic luminosity indicators. We also study the NuSTAR data for NGC 5643 X-1, and show that it exhibits evidence for a spectral cut-off at energy, $E$ $\sim$ 10 keV, similar to that seen in other ULXs observed by NuSTAR. Along with the evidence for significant X-ray luminosity variations in the 3-8 keV band from 2003-2014, our results further strengthen the ULX classification of NGC 5643 X-1.

The intrinsic column density (NH) distribution of quasars is poorly known. At the high obscuration end of the quasar population and for redshifts z<1, the X-ray spectra can only be reliably characterized using broad-band measurements which extend to energies above 10 keV. Using the hard X-ray observatory NuSTAR, along with archival Chandra and XMM-Newton data, we study the broad-band X-ray spectra of nine optically selected (from the SDSS), candidate Compton-thick (NH > 1.5e24 cm^-2) type 2 quasars (CTQSO2s); five new NuSTAR observations are reported herein, and four have been previously published. The candidate CTQSO2s lie at z<0.5, have observed [OIII] luminosities in the range 8.4 < log (L_[OIII]/L_solar) < 9.6, and show evidence for extreme, Compton-thick absorption when indirect absorption diagnostics are considered. Amongst the nine candidate CTQSO2s, five are detected by NuSTAR in the high energy (8-24 keV) band: two are weakly detected at the ~ 3 sigma confidence level and three are strongly detected with sufficient counts for spectral modeling (>~ 90 net source counts at 8-24 keV). For these NuSTAR-detected sources direct (i.e., X-ray spectral) constraints on the intrinsic AGN properties are feasible, and we measure column densities ~2.5-1600 times higher and intrinsic (unabsorbed) X-ray luminosities ~10-70 times higher than pre-NuSTAR constraints from Chandra and XMM-Newton. Assuming the NuSTAR-detected type 2 quasars are representative of other Compton-thick candidates, we make a correction to the NH distribution for optically selected type 2 quasars as measured by Chandra and XMM-Newton for 39 objects. With this approach, we predict a Compton-thick fraction of f_CT = 36^{+14}_{-12} %, although higher fractions (up to 76%) are possible if indirect absorption diagnostics are assumed to be reliable.

We present new NuSTAR and Chandra observations of NGC 3393, a galaxy reported to host the smallest separation dual AGN resolved in the X-rays. While past results suggested a 150 pc separation dual AGN, three times deeper Chandra imaging, combined with adaptive optics and radio imaging suggest a single, heavily obscured, radio-bright AGN. Using VLA and VLBA data, we find an AGN with a two-sided jet rather than a dual AGN and that the hard X-ray, UV, optical, NIR, and radio emission are all from a single point source with a radius <0.2". We find that the previously reported dual AGN is most likely a spurious detection resulting from the low number of X-ray counts (<160) at 6-7 keV and Gaussian smoothing of the data on scales much smaller than the PSF (0.25" vs. 0.80" FWHM). We show that statistical noise in a single Chandra PSF generates spurious dual peaks of the same separation (0.55$\pm$0.07" vs. 0.6") and flux ratio (39$\pm$9% vs. 32% of counts) as the purported dual AGN. With NuSTAR, we measure a Compton-thick source (NH=$2.2\pm0.4\times10^{24}$ cm$^{-2}$) with a large torus half-opening angle, {\theta}=79 which we postulate results from feedback from strong radio jets. This AGN shows a 2-10 keV intrinsic to observed flux ratio of 150. Using simulations, we find that even the deepest Chandra observations would severely underestimate the intrinsic luminosity of NGC 3393 above z>0.2, but would detect an unobscured AGN of this luminosity out to high redshift (z=5).

We aim to constrain the evolution of AGN as a function of obscuration using an X-ray selected sample of $\sim2000$ AGN from a multi-tiered survey including the CDFS, AEGIS-XD, COSMOS and XMM-XXL fields. The spectra of individual X-ray sources are analysed using a Bayesian methodology with a physically realistic model to infer the posterior distribution of the hydrogen column density and intrinsic X-ray luminosity. We develop a novel non-parametric method which allows us to robustly infer the distribution of the AGN population in X-ray luminosity, redshift and obscuring column density, relying only on minimal smoothness assumptions. Our analysis properly incorporates uncertainties from low count spectra, photometric redshift measurements, association incompleteness and the limited sample size. We find that obscured AGN with $N_{H}>{\rm 10^{22}\, cm^{-2}}$ account for ${77}^{+4}_{-5}\%$ of the number density and luminosity density of the accretion SMBH population with $L_{{\rm X}}>10^{43}\text{ erg/s}$, averaged over cosmic time. Compton-thick AGN account for approximately half the number and luminosity density of the obscured population, and ${38}^{+8}_{-7}\%$ of the total. We also find evidence that the evolution is obscuration-dependent, with the strongest evolution around $N_{H}\thickapprox10^{23}\text{ cm}^{-2}$. We highlight this by measuring the obscured fraction in Compton-thin AGN, which increases towards $z\sim3$, where it is $25\%$ higher than the local value. In contrast the fraction of Compton-thick AGN is consistent with being constant at $\approx35\%$, independent of redshift and accretion luminosity. We discuss our findings in the context of existing models and conclude that the observed evolution is to first order a side-effect of anti-hierarchical growth.

Between July 2012 and February 2013, NuSTAR and XMM-Newton performed four long-look joint obser- vations of the type 1.8 Seyfert, NGC 1365. We have analyzed the variable absorption seen in these observations in order to characterize the geometry of the absorbing material. Two of the observations caught NGC 1365 in an unusually low absorption state, revealing complexity in the multi-layer absorber which had previously been hidden. We find the need for three distinct zones of neutral absorption in addition to the two zones of ionized absorption and the Compton-thick torus previously seen in this source. The most prominent absorber is likely associated with broad line region clouds with column densities of around $\sim\,$10$^{23}$ cm$^{-2}$ and a highly clumpy nature as evidenced by an occultation event in February 2013. We also find evidence of a patchy absorber with a variable column around $\sim\,10^{22}$ cm$^{-2}$ and a line of sight covering fraction of 0.3$-$0.9 which responds directly to the intrinsic source flux, possibly due to a wind geometry. A full-covering, constant absorber with a low column density of $\sim\,1 \times$ 10$^{22}$ cm$^{-2}$ is also present, though the location of this low density haze is unknown.

The covering factor of Compton-thick obscuring material associated with the torus in active galactic nuclei (AGN) is at present best understood through the fraction of sources exhibiting Compton-thick absorption along the line of sight ($N_{H}>1.5\times10^{24}$ cm$^{-2}$) in the X-ray band, which reveals the average covering factor. Determining this Compton-thick fraction is difficult however, due to the extreme obscuration. With its spectral coverage at hard X-rays ($>$10 keV), NuSTAR is sensitive to the AGN covering factor since Compton scattering of X-rays off optically thick material dominates at these energies. We present a spectral analysis of 10 AGN observed with NuSTAR where the obscuring medium is optically thick to Compton scattering, so called Compton-thick (CT) AGN. We use the torus models of Brightman & Nandra which predict the X-ray spectrum from reprocessing in a torus and include the torus opening angle as a free parameter and aim to determine the covering factor of the Compton-thick gas in these sources individually. Across the sample we find mild to heavy Compton-thick columns, with $N_{H}$ measured from $10^{24}-10^{26}$ cm$^{-2}$, and a wide range of covering factors, where individual measurements range from 0.2-0.9. We find that the covering factor, $f_{c}$, is a strongly decreasing function of the intrinsic 2-10 keV luminosity, $L_{X}$, where $f_{c}=(-0.41\pm0.13)$log$_{10}$($L_{X}$/erg s$^{-1}$)$+18.31\pm5.33$, across more than two orders of magnitude in $L_{X}$ (10$^{41.5}-10^{44}$ erg s$^{-1}$). The covering factors measured here agree well with the obscured fraction as a function of $L_{X}$ as determined by studies of local AGN with $L_{X}>10^{42.5}$ erg s$^{-1}$.

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].

We present photometric redshifts and associated probability distributions for all detected sources in the Extended Chandra Deep Field South (ECDFS). The work makes use of the most up-to-date data from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) and the Taiwan ECDFS Near-Infrared Survey (TENIS) in addition to other data. We also revisit multi-wavelength counterparts for published X-ray sources from the 4Ms-CDFS and 250ks-ECDFS surveys, finding reliable counterparts for 1207 out of 1259 sources ($\sim 96\%$). Data used for photometric redshifts include intermediate-band photometry deblended using the TFIT method, which is used for the first time in this work. Photometric redshifts for X-ray source counterparts are based on a new library of AGN/galaxy hybrid templates appropriate for the faint X-ray population in the CDFS. Photometric redshift accuracy for normal galaxies is 0.010 and for X-ray sources is 0.014, and outlier fractions are $4\%$ and $5.4\%$ respectively. The results within the CANDELS coverage area are even better as demonstrated both by spectroscopic comparison and by galaxy-pair statistics. Intermediate-band photometry, even if shallow, is valuable when combined with deep broad-band photometry. For best accuracy, templates must include emission lines.

We present X-ray spectral analyses for three Seyfert 2 active galactic nuclei, NGC 424, NGC 1320, and IC 2560, observed by NuSTAR in the 3-79 keV band. The high quality hard X-ray spectra allow detailed modeling of the Compton reflection component for the first time in these sources. Using quasi-simultaneous NuSTAR and Swift/XRT data, as well as archival XMM-Newton data, we find that all three nuclei are obscured by Compton-thick material with column densities in excess of ~5 x $10^{24}$ cm$^{-2}$, and that their X-ray spectra above 3 keV are dominated by reflection of the intrinsic continuum on Compton-thick material. Due to the very high obscuration, absorbed intrinsic continuum components are not formally required by the data in any of the sources. We constrain the intrinsic photon indices and the column density of the reflecting medium through the shape of the reflection spectra. Using archival multi-wavelength data we recover the intrinsic X-ray luminosities consistent with the broadband spectral energy distributions. Our results are consistent with the reflecting medium being an edge-on clumpy torus with a relatively large global covering factor and overall reflection efficiency of the order of 1%. Given the unambiguous confirmation of the Compton-thick nature of the sources, we investigate whether similar sources are likely to be missed by commonly used selection criteria for Compton-thick AGN, and explore the possibility of finding their high-redshift counterparts.

We present Nustar 3-40 keV observations of the optically selected Type 2 quasar (QSO2) SDSS J1034+6001 or Mrk 34. The high-quality hard X-ray spectrum and archival XMM-Newton data can be fitted self-consistently with a reflection-dominated continuum and strong Fe Kalpha fluorescence line with equivalent-width >1 keV. Prior X-ray spectral fitting below 10 keV showed the source to be consistent with being obscured by Compton-thin column densities of gas along the line-of-sight, despite evidence for much higher columns from multiwavelength data. NuSTAR now enables a direct measurement of this column, and shows that Nh lies in the Compton-thick (CT) regime. The new data also show a high intrinsic 2-10 keV luminosity of L_{2-10}~10^{44} erg/s, in contrast to previous low-energy X-ray measurements for which L_{2-10}<~10^{43} erg/s (i.e. X-ray selection below 10 keV does not pick up this source as an intrinsically luminous obscured quasar). Both the obscuring column and the intrinsic power are about an order of magnitude (or more) larger than inferred from pre-NuSTAR X-ray spectral fitting. Mrk34 is thus a 'gold standard' CT QSO2 and is the nearest non-merging system in this class, in contrast to the other local CT quasar NGC6240 which is currently undergoing a major merger coupled with strong star-formation. For typical X-ray bolometric correction factors, the accretion luminosity of Mrk34 is high enough to potentially power the total infrared luminosity. X-ray spectral fitting also shows that thermal emission related to star-formation is unlikely to drive the observed bright soft component below ~3 keV, favoring photionization instead.

We present the results from an X-ray spectral analysis of active galactic nuclei (AGN) in the Chandra Deep Field-South, AEGIS-XD and Chandra-COSMOS surveys, focussing on the identification and characterisation of the most heavily obscured, Compton thick (CT, $N_{H}>10^{24} cm^{-2}$) sources. Our sample is comprised of 3184 X-ray selected extragalactic sources, which has a high rate of redshift completeness (97.6 %), and includes improved photometric redshifts over previous studies. We use spectral models designed for heavily obscured AGN which self consistently include all major spectral signatures of heavy absorption. We identify CT sources not selected through our spectral fitting method using X-ray colours, validate our spectral fitting method through simulations, and take considerations for the constraints on $N_{H}$ given the low count nature of many of our sources. After these considerations we identify a total of 100 CT AGN with best-fit $N_{H}>10^{24} cm^{-2}$ and $N_{H}$ constrained to be above $10^{23.5} cm^{-2}$ at 90% confidence. These sources cover an intrinsic $L_{X}$ range of $10^{42}-3x10^{45}$ ergs s$^{-1}$ and a redshift range of z=0.1-4. This sample will enable characterisation of these heavily obscured AGN across cosmic time and to ascertain their cosmological significance. Furthermore, by using the torus models to test different covering factor scenarios, and by investigating the inclusion of the soft scattered emission, we find evidence that the covering factor of the obscuring material decreases with $L_{X}$ for all redshifts, consistent with the receding torus model, and that this factor increases with redshift, consistent with an increase in the obscured fraction towards higher redshifts. We make $N_{H}$, $\Gamma$ (with uncertainties), observed X-ray fluxes and intrinsic $L_{X}$ for all sources analysed in this work publicly available in an online catalogue.

The Circinus galaxy is one of the nearest obscured AGN, making it an ideal target for detailed study. Combining archival Chandra and XMM-Newton data with new NuSTAR observations, we model the 2-79 keV spectrum to constrain the primary AGN continuum and to derive physical parameters for the obscuring material. Chandra's high angular resolution allows a separation of nuclear and off-nuclear galactic emission. In the off-nuclear diffuse emission we find signatures of strong cold reflection, including high equivalent-width neutral Fe lines. This Compton-scattered off-nuclear emission amounts to 18% of the nuclear flux in the Fe line region, but becomes comparable to the nuclear emission above 30 keV. The new analysis no longer supports a prominent transmitted AGN component in the observed band. We find that the nuclear spectrum is consistent with Compton-scattering by an optically-thick torus, where the intrinsic spectrum is a powerlaw of photon index $\Gamma = 2.2-2.4$, the torus has an equatorial column density of $N_{\rm H} = (6-10)\times10^{24}$cm$^{-2}$ and the intrinsic AGN $2-10$ keV luminosity is $(2.3-5.1)\times 10^{42}$ erg/s. These values place Circinus along the same relations as unobscured AGN in accretion rate-vs-$\Gamma$ and $L_X$-vs-$L_{IR}$ phase space. NuSTAR's high sensitivity and low background allow us to study the short time-scale variability of Circinus at X-ray energies above 10 keV for the first time. The lack of detected variability favors a Compton-thick absorber, in line with the the spectral fitting results.

We report on a NuSTAR and XMM-Newton program that has observed a sample of three extremely luminous, heavily obscured WISE-selected AGN at z~2 in a broad X-ray band (0.1 - 79 keV). The parent sample, selected to be faint or undetected in the WISE 3.4um (W1) and 4.6um (W2) bands but bright at 12um (W3) and 22um (W4), are extremely rare, with only ~1000 so-called W1W2-dropouts across the extragalactic sky. Optical spectroscopy reveals typical redshifts of z~2 for this population, implying rest-frame mid-IR luminosities of L(6um)~6e46 erg/s and bolometric luminosities that can exceed L(bol)~1e14 L(sun). The corresponding intrinsic, unobscured hard X-ray luminosities are L(2-10)~4e45 erg/s for typical quasar templates. These are amongst the most luminous AGN known, though the optical spectra rarely show evidence of a broad-line region and the selection criteria imply heavy obscuration even at rest-frame 1.5um. We designed our X-ray observations to obtain robust detections for gas column densities N(H)<1e24 /cm2. In fact, the sources prove to be fainter than these predictions. Two of the sources were observed by both NuSTAR and XMM-Newton, with neither being detected by NuSTAR and one being faintly detected by XMM-Newton. A third source was observed only with XMM-Newton, yielding a faint detection. The X-ray data require gas column densities N(H)>1e24 /cm2, implying the sources are extremely obscured, consistent with Compton-thick, luminous quasars. The discovery of a significant population of heavily obscured, extremely luminous AGN does not conform to the standard paradigm of a receding torus, in which more luminous quasars are less likely to be obscured. If a larger sample conforms with this finding, then this suggests an additional source of obscuration for these extreme sources.

AGN are known to have complex X-ray spectra that depend on both the properties of the accreting SMBH (e.g. mass, accretion rate) and the distribution of obscuring material in its vicinity ("torus"). Often however, simple and even unphysical models are adopted to represent the X-ray spectra of AGN. In the case of blank field surveys in particular, this should have an impact on e.g. the determination of the AGN luminosity function, the inferred accretion history of the Universe and also on our understanding of the relation between AGN and their host galaxies. We develop a Bayesian framework for model comparison and parameter estimation of X-ray spectra. We take into account uncertainties associated with X-ray data and photometric redshifts. We also demonstrate how Bayesian model comparison can be used to select among ten different physically motivated X-ray spectral models the one that provides a better representation of the observations. Despite the use of low-count spectra, our methodology is able to draw strong inferences on the geometry of the torus. For a sample of 350 AGN in the 4 Ms Chandra Deep Field South field, our analysis identifies four components needed to represent the diversity of the observed X-ray spectra: (abridged). Simpler models are ruled out with decisive evidence in favour of a geometrically extended structure with significant Compton scattering. Regarding the geometry of the obscurer, there is strong evidence against both a completely closed or entirely open toroidal geometry, in favour of an intermediate case. The additional Compton reflection required by data over that predicted by toroidal geometry models, may be a sign of a density gradient in the torus or reflection off the accretion disk. Finally, we release a catalogue with estimated parameters such as the accretion luminosity in the 2-10 keV band and the column density, $N_{H}$, of the obscurer.

Observations in the last decade have provided strong evidence that the growth of supermassive black holes at the centres of galaxies is among the most influential processes in galaxy evolution. Open questions that relate to our current understanding of black hole growth and its relation to the build-up of galaxies at redshifts z=1-4, when most black holes and stars we see in present-day galaxies were put in place, include: what is the nature of AGN feedback and whether it plays a significant role in the evolution of galaxies? what is the dominant population of accreting AGN at that critical epoch? is it dominated by obscured objects as required by many current observations and models? The Athena+ mission concept will provide the technological leap required for a breakthrough in our understanding of AGN and galaxy evolution at the heyday of the Universe. The high throughput of Athena+ will allow the systematic study of the incidence, nature and energetics of AGN feedback processes to z~4 via the identification and measurement of blue-shifted X-ray absorption lines. The excellent survey and spectral capabilities of the Athena+ Wide Field Imager will complete the census of black hole growth by yielding samples of up to 100 times larger than is currently possible of the most heavily obscured, including Compton thick, AGN to redshifts z~4. The demographics of this population relative to their hosts is fundamental for understanding how major black hole growth events relate to the build-up of galaxies.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

We present an investigation into how well the properties of the accretion flow onto a supermassive black hole may be coupled to those of the overlying hot corona. To do so, we specifically measure the characteristic spectral index, Gamma, of a power-law energy distribution, over an energy range of 2 to 10 keV, for X-ray selected, broad-lined radio-quiet AGN up to z~2 in COSMOS and E-CDF-S. We test the previously reported dependence between Gamma and black hole mass, FWHM and Eddington ratio using a sample of AGN covering a broad range in these parameters based on both the Mg ii and H-alpha emission lines with the later afforded by recent near infrared spectroscopic observations using Subaru/FMOS. We calculate the Eddington ratios, lambda_Edd, for sources where a bolometric luminosity (L_Bol) has been presented in the literature, based on SED fitting, or, for sources where these data do not exist, we calculate L_Bol using a bolometric correction to the LX, derived from a relationship between the bolometric correction, and LX/L3000. From a sample of 69 X-ray bright sources (>250 counts), where Gamma can be measured with greatest precision, with an estimate of L_Bol, we find a statistically significant correlation between Gamma and lambda_Edd, which is highly significant with a chance probability of 6.59x10^-8. A statistically significant correlation between Gamma and the FWHM of the optical lines is confirmed, but at lower significance than with lambda_Edd indicating that lambda_Edd is the key parameter driving conditions in the corona. Linear regression analysis reveals that Gamma=(0.32+/-0.05)log10 lambda_Edd+(2.27+/-0.06) and Gamma=(-0.69+/-0.11)log10(FWHM/km/s)+(4.44+/-0.42). Our results on Gamma-lambda_Edd are in very good agreement with previous results. (ABRIDGED)

The existence of an anti-correlation between the equivalent width (EW) of the narrow core of the iron Kalpha line and the luminosity of the continuum (i.e. the X-ray Baldwin effect) in type-I active galactic nuclei has been confirmed over the last years by several studies carried out with XMM-Newton, Chandra and Suzaku. However, so far no general consensus on the origin of this trend has been reached. Several works have proposed the decrease of the covering factor of the molecular torus with the luminosity (in the framework of the luminosity-dependent unification models) as a possible explanation for the X-ray Baldwin effect. Using the fraction of obscured sources measured by recent X-ray and IR surveys as a proxy of the half-opening angle of the torus, and the recent Monte-Carlo simulations of the X-ray radiation reprocessed by a structure with a spherical-toroidal geometry by Ikeda et al. (2009) and Brightman & Nandra (2011), we test the hypothesis that the X-ray Baldwin effect is related to the decrease of the half-opening angle of the torus with the luminosity. Simulating the spectra of an unabsorbed population with a luminosity-dependent covering factor of the torus as predicted by recent X-ray surveys, we find that this mechanism is able to explain the observed X-ray Baldwin effect. Fitting the simulated data with a log-linear L_{2-10keV}-EW relation, we found that in the Seyfert regime (L_{2-10keV}< 10^44.2 erg s^-1) luminosity-dependent unification produces a slope consistent with the observations for average values of the equatorial column densities of the torus of log N_H^T > 23.1. In the quasar regime (L_{2-10 keV}> 10^44.2 erg s^-1) a decrease of the covering factor of the torus with the luminosity slower than that observed in the Seyfert regime (as found by recent hard X-ray surveys) is able to reproduce the observations for 23.2 < log N_H^T < 24.2.

This paper investigates what constraints can be placed on the fraction of Compton-thick (CT) AGN in the Universe from the modeling of the spectrum of the diffuse X-ray background (XRB). We present a model for the synthesis of the XRB that uses as input a library of AGN X-ray spectra generated by the Monte Carlo simulations described by Brightman & Nandra. This is essential to account for the Compton scattering of X-ray photons in a dense medium and the impact of that process on the spectra of obscured AGN. We identify a small number of input parameters to the XRB synthesis code which encapsulate the minimum level of uncertainty in reconstructing the XRB spectrum. These are the power-law index and high energy cutoff of the intrinsic X-ray spectra of AGN, the level of the reflection component in AGN spectra and the fraction of CT AGN in the Universe. We then map the volume of the space allowed to these parameters by current observations of the XRB spectrum in the range 3-100 keV. One of the least constrained parameters is the fraction of CT AGN. Statistically acceptable fits to the XRB spectrum at the 68% confidence level can be obtained for CT fractions in the range 5-50%. This is because of degeneracies among input parameters to the XRB synthesis code and uncertainties in the modeling of AGN spectra (e.g. reflection). The most promising route for constraining the fraction of CT AGN in the Universe is via the direct detection of those sources in high energy (>10keV) surveys. It is shown that the observed fraction of CT sources identified in the SWIFT/BAT survey, limits the intrinsic fraction of CT AGN, at least at low redshift, to 10-20% (68% confidence level). We also make predictions on the number density of CT sources that current and future X-ray missions are expected to discover. Testing those predictions will constrain the intrinsic fraction of CT AGN as a function of redshift.

We present the results of a campaign of simultaneous X-ray and optical observations of `true' Type 2 Seyfert galaxies candidates, i.e. AGN without a Broad Line Region (BLR). Out of the initial sample composed by 8 sources, one object, IC1631, was found to be a misclassified starburst galaxy, another, Q2130-431, does show broad optical lines, while other two, IRAS01428-0404 and NGC4698, are very likely absorbed by Compton-thick gas along the line of sight. Therefore, these four sources are not unabsorbed Seyfert 2s as previously suggested in the literature. On the other hand, we confirm that NGC3147, NGC3660, and Q2131-427 belong to the class of true Type 2 Seyfert galaxies, since they do not show any evidence for a broad component of the optical lines nor for obscuration in their X-ray spectra. These three sources have low accretion rates ($\dot{m}= L_{bol}/L_{Edd} \la0.01$), in agreement with theoretical models which predict that the BLR disappears below a critical value of $L_{bol}/L_{Edd}$. The last source, Mrk273x, would represent an exception even of this accretion-dependent versions of the Unification Models, due to its high X-ray luminosity and accretion rate, and no evidence for obscuration. However, its optical classification as a Seyfert 2 is only based on the absence of a broad component of the H\beta, due to the lack of optical spectra encompassing the H\alpha\ band.

(Abridged) We present the results from the X-ray spectral analysis of high-z AGN in the CDFS, making use of the new 4Ms data set and new X-ray spectral models from Brightman & Nandra, which account for Compton scattering and the geometry of the circumnuclear material. Our goals are to ascertain to what extent the torus paradigm of local AGN is applicable at earlier epochs and to evaluate the evolution of the Compton thick fraction (f_CT) with z, important for XRB synthesis models and understanding the accretion history of the universe. In addition to the torus models, we measure the fraction of scattered nuclear light, f_scatt known to be dependant on covering factor of the circumnuclear materal, and use this to aid in our understanding of its geometry. We find that the covering factor of the circumnuclear material is correlated with NH, and as such the most heavily obscured AGN are in fact also the most geometrically buried. We come to these conclusions from the result that f_scatt decreases as NH increases and from the prevalence of the torus model with the smallest opening angle as best fit model in the fits to the most obscured AGN. We find that a significant fraction of sources (~ 20%) in the CDFS are likely to be buried in material with close to 4 pi coverage having been best fit by the torus model with a 0\degree opening angle. Furthermore, we find 41 CTAGN in the CDFS using the new torus models, 29 of which we report here for the first time. We bin our sample by z in order to investigate the evolution of f_CT. Once we have accounted for biases and incompleteness we find a significant increase in the intrinsic f_CT, normalised to LX= 10^43.5 erg/s, from \approx 20% in the local universe to \approx 40% at z=1-4.

We present a method for the identification of heavily absorbed AGN (NH>10^23 cm^-2) from X-ray photometric data. We do this using a set of XMM-Newton reference spectra of local galaxies for which we have accurate NH information, as described in Brightman & Nandra. The technique uses two rest-frame hardness ratios which are optimised for this task, which we designate HR1 (2-4/1-2 keV) and HR2 (4-16/2-4 keV). The selection method exploits the fact that while obscured AGN appear hard in HR2 due to absorption of the intrinsic source flux below ~4 keV, they appear soft in HR1 due to excess emission originating from scattered source light, thermal emission, or host galaxy emission. Such emission is ubiquitous in low redshift samples. The technique offers a very simple and straight forward way of estimating the fraction of obscured AGN in samples with relatively low signal-to-noise ratio in the X-ray band. We apply this technique to a moderate redshift (z~1) sample of AGN from the Chandra Deep Field North, finding that 61% of this sample has NH> 10^23 cm^-2. A clear and robust conclusion from our analysis, is that in deep surveys the vast majority of sources do not show hardness ratios consistent with a simple absorbed power-law. The ubiquity of complex spectra in turn shows that simple hardness ratio analysis will not yield reliable obscuration estimates, justifying the more complex colour-colour analysis described in this paper. While this method does very well at separating sources with NH> 10^23 cm^-2 from sources with lower NH, only X-ray spectroscopy can identify Compton thick sources, through the detection of the Fe Ka line. This will be made possible with the high throughput X-ray spectral capabilities of ATHENA.

We present a multi-waveband analysis of a 126 galaxy sub-sample of the 12 micron galaxy sample (12MGS), for which we have carried out a detailed X-ray spectral analysis of in a previous paper. We present optical activity classes and we characterise these by their X-ray, 12 um and [OIII] luminosities and by their X-ray spectral properties. Our most interesting results are: Seyfert (Sy) 1s and Sy 2s show significantly different LX distributions from each other; Sy 2s with a detection of an HBLR show a significantly higher LX than those without a detection; Sy 1s also present a significantly different L(12um) distribution from both Sy1.2-1.9s and Sy 2s; the Sy 2 fraction decreases towards high LX; X-ray indications of AGN power agree well with the optical classifications; There is X-ray evidence for the presence of an AGN in 17% of HII/AGN composites and 40% of LINERs; we advocate the use of a LX >10^41 erg/s in the X-ray selection of AGN, which gives a contamination rate of only 3% from star-forming galaxies; Sy 1s and Sy 2 have the same intrinsic {\Gamma}; in 24% of cases the absorption measured in X-rays does not correspond directly with that implied in the optical band from the visibility of the BLRs; the obscured fraction in AGN is a strong function of LX peaking at 10^42 erg/s; The average obscured and Compton thick fractions for this sample are 62\pm5% and 20\pm4% respectively, which are higher than hard X-ray and optically selected samples; we assess the use of the T ratio (FX/F[OIII]) for selecting Compton thick candidates. We conclude here that this quantity can often be unreliable due uncertainties in the extinction corrections to the [OIII] flux. These results have important impacts on AGN selection and unification and the results from the 12MGS are particularly useful as a local analogue to z=1 24 um samples.

We present an X-ray spectral analysis of 126 galaxies of the 12 micron galaxy sample. We pay particular attention to Compton thick AGN with the help of new spectral fitting models that we have produced, which are based on Monte-Carlo simulations of X-ray radiative transfer, using both a spherical and torus geometry, and taking into account Compton scattering and Fe fluorescence. We use this data to show that with a torus geometry, unobscured sight lines can achieve a maximum EW of the Fe K\alpha line of ~150 eV, originally shown by Ghisellini, Haardt & Matt (1994). In order for this to be exceeded, the line of sight must be obscured with N_H>10^23 cm^-2, as we show for one case, NGC 3690. We also calculate flux suppression factors from the simulated data, the main conclusion from which is that for N_H>10^25 cm^-2, the X-ray flux is suppressed by a factor of >10 in all X-ray bands and at all redshifts, revealing the biases present against these extremely heavily obscured systems inherent in all X-ray surveys. Furthermore, we confirm previous results from Murphy & Yaqoob (2009) that show that the reflection fraction determined from slab geometries is underestimated with respect to toroidal geometries. For the 12 micron selected galaxies, we investigate the distribution of X-ray power-law indices, finding that the mean $<\Gamma>=1.90_{-0.07}^{+0.05}$ and $\sigma_\Gamma = 0.31_{-0.05}^{+0.05}$) is consistent with previous works, and that the distribution of \Gamma for obscured and unobscured sources is consistent with the source populations being the same, in general support of unification schemes. We determine a Compton thick fraction for the X-ray AGN in our sample to be 18+/-5% which is higher than the hard X-ray (>10 keV) selected samples. Finally we find that the obscured fraction for our sample is a strong function of X-ray luminosity, peaking at L_X~10^42-43 ergs s^-1.

We present an analysis of six 12 um selected Seyfert 2 galaxies that have been reported to be unabsorbed in the X-ray. By comparing the luminosities of these galaxies in the mid-IR (12um), optical ([O III]) and hard X-ray (2-10 keV), we show that they are all under-luminous in the 2-10 keV X-ray band. Four of the objects exhibit X-ray spectra indicative of a hard excess, consistent with a heavily obscured X-ray component and hence a hidden nucleus. In these objects the softer X-rays may be dominated by a strong soft scattered continuum or contamination from the host galaxy, which is responsible for the unabsorbed X-ray spectra observed, and accounts for the anomalously low 2-10 keV X-ray luminosity. We confirm this assertion in NGC4501 with a Chandra observation, which shows hard X-ray emission coincident with the nucleus, consistent with heavy absorption, and a number of contaminating softer sources which account for the bulk of the softer emission. We point out that such 'Compton thick' sources need not necessarily present iron Ka emission of high EW. An example in our sample is IRASF01475-0740, which we know must host an obscured AGN as it hosts a HBLR seen in scattered light (Tran 2003). The X-ray spectrum is nonetheless relatively unobscured and the iron Ka line only moderate in strength. These observations can be reconciled if the hidden nuclear emission is dominated by transmitted, rather than reflected X-rays, which can then be weak compared to the soft scattered light or galactic emission. Despite these considerations, we conclude that NGC3147 & 3660 may intrinsically lack a BLR, confirming the recent results of Bianchi et al. (2008) for NGC3147. Neither X-ray spectrum shows signs of hidden hard emission and both sources exhibit X-ray variability leading us to believe we are viewing the nucleus directly.