Search

We present preliminary results of a Chandra Large Program to monitor the ultraluminous X-ray source (ULX) populations of three nearby, ULX-rich galaxies over the course of a year, finding the ULX population to show a variety of long-term variability behaviours. Of a sample of 36 ULXs, some show persistent or moderately variable flux, often with a significant relationship between hardness and luminosity, consistent with a supercritically accreting source with varying accretion rates. Six show very high-amplitude variability with no strong relationship between luminosity and hardness, though not all of them show evidence of any long-term periodicity, nor of the bimodal distribution indicative of the propeller effect. We find evidence of additional eclipses for two previously-identified eclipsing ULXs. Additionally, many sources that were previously identified as ULXs in previous studies were not detected at ULX luminosities during our monitoring campaign, indicating a large number of transient ULXs., Comment: 7 pages, 4 figures, to be published in Proceedings of the XMM-Newton Workshop 2024 "The X-ray Mysteries of Neutron Stars and White Dwarfs", Astronomische Nachrichten, in press

Recent models suggest approximately half of all accreting supermassive black holes (SMBHs; $M_{\rm BH}$ $\gtrsim$ 10$^{5}$ M$_{\odot}$) are expected to undergo intense growth phases behind Compton-thick ($N_{\rm H}$ $>$ 1.5 $\times$ 10$^{24}$ cm$^{-2}$) veils of obscuring gas. However, despite being a viable source for the seeding of SMBHs, there are currently no examples known of a Compton-thick accreting intermediate mass black hole (IMBH; $M_{\rm BH}$ $\sim$ 10$^{2}$ $-$ 10$^{5}$ M$_{\odot}$). We present a detailed X-ray spectral analysis of IC 750 $-$ the only AGN to-date with a precise megamaser-based intermediate mass $<$ 10$^{5}$ M$_{\odot}$. We find the equivalent width of neutral 6.4 keV Fe K$\alpha$ to be 1.9$^{+2.2}_{-1.0}$ keV via phenomenological modelling of the co-added 177 ks Chandra spectrum. Such large equivalent widths are seldom produced by processes other than fluorescence from dense obscuration. We fit three physically-motivated X-ray spectral models to infer a range of possible intrinsic 2$-$10 keV luminosity posteriors that encompass the systematic uncertainties associated with a choice of model. Despite a wide range of predicted intrinsic 2$-$10 keV luminosities between $\sim$ 10$^{41}$ and 10$^{43}$ erg s$^{-1}$, all three models agree that IC 750 has a Compton-thick line-of-sight column density to $>$ 99\% confidence. Compton-thick obscuration is well-documented to impinge substantial bias on the pursuit of SMBH AGN. Our results thus provide the first indication that Compton-thick obscuration should also be properly considered to uncover and understand the IMBH population in an unbiased manner., Comment: Accepted for publication in ApJ. 15 pages, 6 figures

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

The discovery of pulsations in (at least) six ultraluminous X-ray sources (ULXs) has shown that neutron stars can accrete at (highly) super-Eddington rates, challenging the standard accretion theories. M51 ULX-7, with a spin signal of $P\simeq2.8$ s, is the pulsating ULX (PULX) with the shortest known orbital period ($P_\mathrm{orb}\simeq2$ d) and has been observed multiple times by XMM-Newton, Chandra, and NuSTAR. We report on the timing and spectral analyses of three XMM-Newton observations of M51 ULX-7 performed between the end of 2021 and the beginning of 2022, together with a timing re-analysis of XMM-Newton, Chandra, and NuSTAR archival observations. We investigated the spin signal by applying accelerated search techniques and studied the power spectrum through the fast Fourier transform, looking for (a)periodic variability in the source flux. We analysed the energy spectra of the 2021-2022 observations and compared them to the older ones. We report the discovery of a recurrent, significant ($>$3$\sigma$) broad complex at mHz frequencies in the power spectra of M51 ULX-7. We did not detect the spin signal, setting a 3$\sigma$ upper limit on the pulsed fraction of $\lesssim10\%$ for the single observation. The complex is significantly detected also in five Chandra observations performed in 2012. M51 ULX-7 represents the second PULX for which we have a significant detection of mHz-QPOs at super-Eddington luminosities. These findings suggest that one should avoid using the observed QPO frequency to infer the mass of the accretor in a ULX. The absence of spin pulsations when the broad complex is detected suggests that the mechanism responsible for the aperiodic modulation also dampens the spin signal's pulsed fraction. If true, this represents an additional obstacle in the detection of new PULXs, suggesting an even larger occurrence of PULXs among ULXs., Comment: 14 pages (12 main text + Appendix), 5 figures. Accepted for publication on A&A

We observed the nearby and relatively understudied ultraluminous X-ray source (ULX) NGC 4190 ULX-1 jointly with NICER and NuSTAR to investigate its broadband spectrum, timing properties, and spectral variation over time. We found NGC 4190 ULX-1 to have a hard spectrum characterized by two thermal components (with temperatures ~0.25keV and ~1.6keV) and a high-energy excess typical of the ULX population, although the spectrum turns over at an unusually low energy. While no pulsations were detected, (with pulsed fraction 3-sigma upper limits of 16% for NICER and 35% for NuSTAR), the source shows significant stochastic variability and the covariance spectrum indicates the presence of a high-energy cut-off power-law component, potentially indicative of an accretion column. Additionally, when fitting archival XMM-Newton data with a similar model, we find that the luminosity-temperature evolution of the hot thermal component follows the behavior of a super-Eddington slim disk though the expected spectral broadening for such a disk is not seen, suggesting that the inner accretion disk may be truncated by a magnetic field. Therefore, despite the lack of detected pulsations, there is tantalizing evidence for NGC 4190 ULX-1 being a candidate neutron star accretor, although further broadband observations will be required to confirm this behavior., Comment: 14 pages, 9 figures, 2 tables, accepted for publication in The Astrophysical Journal

HEX-P will launch at a time when the sky is being routinely scanned for transient gravitational wave, electromagnetic and neutrino phenomena that will require the capabilities of a sensitive, broadband X-ray telescope for follow up studies. These include the merger of compact objects such as neutron stars and black holes, stellar explosions, and the birth of new compact objects. \hexp\ will probe the accretion and ejecta from these transient phenomena through the study of relativistic outflows and reprocessed emission, provide unique capabilities for understanding jet physics, and potentially revealing the nature of the central engine., Comment: Submitted to Frontiers in Astronomy and Space Sciences

Ultraluminous X-ray sources (ULXs) represent an extreme class of accreting compact objects: from the identification of some of the accretors as neutron stars to the detection of powerful winds travelling at 0.1-0.2 c, the increasing evidence points towards ULXs harbouring stellar-mass compact objects undergoing highly super-Eddington accretion. Measuring their intrinsic properties, such as the accretion rate onto the compact object, the outflow rate, the masses of accretor/companion -- hence their progenitors, lifetimes, and future evolution -- is challenging due to ULXs being mostly extragalactic and in crowded fields. Yet ULXs represent our best opportunity to understand super-Eddington accretion physics and the paths through binary evolution to eventual double compact object binaries and gravitational wave sources. Through a combination of end-to-end and single-source simulations, we investigate the ability of HEX-P to study ULXs in the context of their host galaxies and compare it to XMM-Newton and NuSTAR, the current instruments with the most similar capabilities. HEX-P's higher sensitivity, which is driven by its narrow point-spread function and low background, allows it to detect pulsations and broad spectral features from ULXs better than XMM-Newton and NuSTAR. We describe the value of HEX-P in understanding ULXs and their associated key physics, through a combination of broadband sensitivity, timing resolution, and angular resolution, which make the mission ideal for pulsation detection and low-background, broadband spectral studies., Comment: 24 pages, 9 figures. Accepted for publication in Frontiers of Astronomy and Space Science. Minor corrections included

We present NuSTAR observations of the nearby SN 2023ixf in M101 (d=6.9 Mpc) which provide the earliest hard X-ray detection of a non-relativistic stellar explosion to date at $\delta t\approx$4-d and $\delta t\approx$11-d. The spectra are well described by a hot thermal bremsstrahlung continuum with $T>25 \rm{keV}$ shining through a thick neutral medium with a neutral hydrogen column that decreases with time (initial $N_{\rm{Hint}}=2.6\times 10^{23} \rm{cm^{-2}}$). A prominent neutral Fe K$\alpha$ emission line is clearly detected, similar to other strongly interacting SNe such as SN2020jl. The rapidly decreasing intrinsic absorption with time suggests the presence of a dense but confined circumstellar medium (CSM). The absorbed broadband X-ray luminosity (0.3--79 keV) is $L_{X} \approx 2.5 \times 10^{40}$ erg s$^{-1}$ during both epochs, with the increase in overall X-ray flux related to the decrease in the absorbing column. Interpreting these observations in the context of thermal bremsstrahlung radiation originating from the interaction of the SN shock with a dense medium we infer large particle densities in excess of $n_{\rm{CSM}}\approx 4\times 10^{8} \rm{cm^{-3}}$ at $r<10^{15} \rm{cm}$, corresponding to an enhanced progenitor mass-loss rate of $\dot M \approx 3\times 10^{-4}$ M$_{\odot}$ yr$^{-1}$ for an assumed wind velocity of $v_w=50$ km s$^{-1}$., Comment: 8 pages, 2 figures. Submitted

Ultraluminous X-ray sources (ULXs) are our best laboratories for studying extreme super-Eddington accretion. Most studies of these objects are of relatively persistent sources, however there is growing evidence to suggest a large fraction of these sources are transient. Here we present a sample of five newly reported transient ULXs in the galaxies NGC 4945, NGC 7793 and M81 serendipitously discovered in Swift/XRT observations. Swift monitoring of these sources have provided well sampled lightcurves, allowing for us to model the lightcurves with the disk instability model of Hameury & Lasota (2020) which implies durations of 60-400 days and that the mass accretion rate through the disk is close to or greater than the Eddington rate. Of the three source regions with prior HST imaging, color magnitude diagrams of the potential stellar counterparts show varying ages of the possible stellar counterparts. Our estimation of the rates of these sources in these three galaxies is 0.4-1.3 year$^{-1}$. We find that while persistent ULXs dominate the high end of galaxy luminosity functions, the number of systems that produce ULX luminosities are likely dominated by transient sources., Comment: Accepted for publication in ApJ

We present analysis of NuSTAR X-ray observations of three AGN that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs. We previously observed these AGN with Chandra and found no differences between their low energy X-ray properties and the larger AGN population. However some models predict differences to be more prominent at energies higher than probed by Chandra. We find that even at the higher energies probed by NuSTAR, the spectra of these AGN are indistinguishable from the larger AGN population. This could rule out models predicting large differences in the X-ray spectra in the NuSTAR bands. Alternatively, it might mean that these three AGN are not binary SMBHs.

We present optical, radio and X-ray observations of a rapidly-evolving transient AT2019wxt (PS19hgw), discovered during the search for an electromagnetic (EM) counterpart to the gravitational-wave (GW) trigger S191213g (LIGO Scientific Collaboration & Virgo Collaboration 2019a). Although S191213g was not confirmed as a significant GW event in the off-line analysis of LIGO-Virgo data, AT2019wxt remained an interesting transient due its peculiar nature. The optical/NIR light curve of AT2019wxt displayed a double-peaked structure evolving rapidly in a manner analogous to currently know ultra-stripped supernovae (USSNe) candidates. This double-peaked structure suggests presence of an extended envelope around the progenitor, best modelled with two-components: i) early-time shock-cooling emission and ii) late-time radioactive $^{56}$Ni decay. We constrain the ejecta mass of AT2019wxt at $M_{ej} \approx{0.20 M_{\odot}}$ which indicates a significantly stripped progenitor that was possibly in a binary system. We also followed-up AT2019wxt with long-term Chandra and Jansky Very Large Array observations spanning $\sim$260 days. We detected no definitive counterparts at the location of AT2019wxt in these long-term X-ray and radio observational campaigns. We establish the X-ray upper limit at $9.93\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ and detect an excess radio emission from the region of AT2019wxt. However, there is little evidence for SN1993J- or GW170817-like variability of the radio flux over the course of our observations. A substantial host galaxy contribution to the measured radio flux is likely. The discovery and early-time peak capture of AT2019wxt in optical/NIR observation during EMGW follow-up observations highlights the need of dedicated early, multi-band photometric observations to identify USSNe., Comment: 20 pages, 14 figures, Submitted to ApJ

We report the results of a 2019-2021 monitoring campaign with Swift and associated target-of-opportunity observations with XMM-Newton and NuSTAR, examining the spectral and timing behavior of the highly variable ultraluminous X-ray source (ULX) NGC 925 ULX-3. We find that the source exhibits a 127-128 day periodicity, with fluxes typically ranging from 1e-13 to 8e-13 ergs/s/cm2. We do not find strong evidence for a change in period over the time that NGC 925 ULX-3 has been observed, although the source may have been in a much lower flux state when first observed with Chandra in 2005. We do not detect pulsations, and we place an upper limit on the pulsed fraction of ~40% in the XMM-Newton band, consistent with some previous pulsation detections at low energies in other ULXs. The source exhibits a typical ULX spectrum that turns over in the NuSTAR band and can be fitted using two thermal components. These components have a high temperature ratio that may indicate the lack of extreme inner disk truncation by a magnetar-level magnetic field. We examine the implications for a number of different models for superorbital periods in ULXs, finding that a neutron star with a magnetic field of ~10^12 G may be plausible for this source. The future detection of pulsations from this source would allow for the further testing and constraining of such models., Comment: 13 pages, 6 figures, 2 tables. Accepted for publication in ApJ

We report on VLT/FORS2 imaging polarimetry observations in the $R_{\rm special}$ band of WISE J011601.41-050504.0 (W0116-0505), a heavily obscured hyper-luminous quasar at $z=3.173$ classified as a Hot, Dust-Obscured Galaxy (Hot DOG) based on its mid-IR colors. Recently, Assef et al. (2020) identified W0116-0505 as having excess rest-frame optical/UV emission, and concluded this excess emission is most likely scattered light from the heavily obscured AGN. We find that the broad-band rest-frame UV flux is strongly linearly polarized (10.8$\pm$1.9\%, with a polarization angle of 74$\pm$9~deg), confirming this conclusion. We analyze these observations in the context of a simple model based on scattering either by free electrons or by optically thin dust, assuming a classical dust torus with polar openings. Both can replicate the degree of polarization and the luminosity of the scattered component for a range of geometries and column densities, but we argue that optically thin dust in the ISM is the more likely scenario. We also explore the possibility that the scattering medium corresponds to an outflow recently identified for W0116-0505. This is a feasible option if the outflow component is bi-conical with most of the scattering occurring at the base of the receding outflow. In this scenario the quasar would still be obscured even if viewed face on, but might appear as a reddened type 1 quasar once the outflow has expanded. We discuss a possible connection between blue-excess Hot DOGs, extremely red quasars (ERQs), reddened type 1 quasars, and unreddened quasars that depends on a combination of evolution and viewing geometry., Comment: 19 pages, 10 figures, 2 tables. Resubmitted to ApJ after first round of referee comments

The memo describes the methods used to track the long-term gain variations in the NuSTAR detectors. It builds on the analysis presented in Madsen et al. (2015) using the deployable calibration source to measure the gain drift in the NuSTAR CdZnTe detectors. This is intended to be a live document that is periodically updated as new entries are required in the NuSTAR gain CALDB files. This document covers analysis up through early-2022 and the gain v011 CALDB file released in version 20240226., Comment: 11 page, 7 figures. Intended as a living, easy-to-find document. No intention of submitting this to a journal

Type 2 active galactic nuclei (AGN) show signatures of accretion onto a supermassive black hole through strong, high-ionization, narrow emission lines extended on scales of 100s to 1000s of parsecs, but they lack the broad emission lines from close in to the black hole that characterize type 1 AGN. The lack of broad emission could indicate obscuration of the innermost nuclear regions, or could indicate that the black hole is no longer strongly accreting. Since high-energy X-rays can penetrate thick obscuring columns, they have the power to distinguish these two scenarios. We present high-energy NuSTAR observations of 9 Seyfert 2 AGN from the IRAS 12 micron survey, supplemented with low-energy X-ray observations from Chandra, XMM-Newton, and Swift. The galaxies were selected to have anomalously low observed 2-10 keV luminosities compared to their [O III] optical luminosities, a traditional diagnostic of heavily obscured AGN, reaching into the Compton-thick regime for the highest hydrogen column densities ($N_{\rm H} > 1.5 \times 10^{24}\, {\rm cm}^{-2}$). Based on updated [O III] luminosities and intrinsic X-ray luminosities based on physical modeling of the hard X-ray spectra, we find that one galaxy was misclassified as type 2 (NGC 5005) and most of the remaining AGN are obscured, including three confirmed as Compton-thick (IC 3639, NGC 1386, and NGC 3982). One galaxy, NGC 3627, appears to be recently deactivated. Compared to the original sample the 9 AGN were selected from, this is a rate of approximately 1%. We also find a new X-ray changing-look AGN in NGC 6890., Comment: 42 pages, Accepted for publication in ApJ

We report on the discovery of an absorption line at $E=8.56^{+0.05}_{-0.11}$ keV detected with a significance of $>3.3\sigma$ in the NuSTAR and XMM-Newton spectra of a newly discovered hyperluminous X-ray source (HLX, $L_{\rm X}>10^{41}$ ergs$^{-1}$) in the galaxy NGC 4045 at a distance of 32 Mpc. The source was first discovered serendipitously in a Swift/XRT observation of the galaxy, and Swift monitoring reveals a highly variable source changing by over an order of magnitude from maximum to minimum. The origin of the absorption line appears likely to be by highly ionized iron with a blue shift of 0.19$c$, indicating an ultrafast outflow (UFO). However, the large equivalent width of the line (EW$=-0.22^{+0.08}_{-0.09}$ keV) paired with the lack of other absorption lines detected are difficult to reconcile with models. An alternative explanation is that the line is due to a cyclotron resonance scattering feature (CRSF), produced by the interaction of X-ray photons with the powerful magnetic field of a neutron star., Comment: Accepted for publication in ApJ

This paper describes a method by which the metrology system of the Nuclear Spectroscopic Telescope Array (NuSTAR) X-ray space observatory, which uses two lasers to characterize the relative motion of the optics and focal plane benches, can be approximated should one laser fail. The two benches are separated by a ten-meter-long rigid mast that undergoes small amounts of thermal flexing which need to be compensated for in order to produce a non-blurred image. We analyze the trends of mast motion by archival observation parameters in order to discover whether the mast motion in future observations can be predicted. We find that, by using the solar aspect angle (SAA), observation date, and orbital phase, we can simulate the motion of one laser by translating the track produced by the other and applying modifications to the resulting mast aspect solution, allowing the reconstruction of a minimally distorted point spread function in most cases. We will implement the generation of simulated mast files alongside the usual NuSTAR data reduction pipeline for contingency purposes. This work has implications for reducing the risk of implementing laser metrology systems on future missions that use deployable masts to achieve the long focal lengths required in high-energy astronomy by mitigating the impact of a metrology laser failure in the extended phase of a mission., Comment: 36 pages, 20 figures, 1 table. Accepted for publication in JATIS

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

M82 X-2 is the first pulsating ultraluminous X-ray source (PULX) discovered. The luminosity of these extreme pulsars, if isotropic, implies an extreme mass transfer rate. An alternative is to assume a much lower mass transfer rate, but with an apparent luminosity boosted by geometrical beaming. Only an independent measurement of the mass transfer rate can help discriminate between these two scenarios. In this Paper, we follow the orbit of the neutron star for seven years, measure the decay of the orbit ($\dot{P}_{orb}/{P}_{orb}\approx-8\cdot10^{-6}\mathrm{yr}^{-1}$), and argue that this orbital decay is driven by extreme mass transfer of more than 150 times the mass transfer limit set by the Eddington luminosity. If this is true, the mass available to the accretor is more than enough to justify its luminosity, with no need for beaming. This also strongly favors models where the accretor is a highly-magnetized neutron star., Comment: 18 pages, 6 figures, accepted for publication in ApJ

M51 ULX7 is among a small group of known ultraluminous X-ray pulsars (ULXP). The neutron star powering the source has a spin period of 2.8s, orbits its companion star with a period of 2 days, and a super-orbital period of 38 days is evident in its X-ray lightcurve. Here we present NuSTAR and XMM-Newton data on the source from 2019 obtained when the source was near its peak brightness. We detect the pulsations, having spun up at a rate of 3$\pm0.5\times10^{-10}$ s s$^{-1}$ since they were previously detected in 2018. The data also provide the first high-quality broadband spectrum of the source. We find it to be very similar to that of other ULXPs, with two disk-like components, and a high energy tail. When combined with XMM-Newton data obtained in 2018, we explore the evolution of the spectral components with super-orbital phase, finding that the luminosity of the hotter component drives the super-orbital flux modulation. The inclination the disk components appear to change with phase, which may support the idea that these super-orbital periods are caused by disk precession. We also reexamine the super-orbital period with 3 years of Swift/XRT monitoring, finding that the period is variable, increasing from 38.2$\pm0.5$ days in 2018--2019 to 44.2$\pm0.9$ days in 2020--2021, which rules out alternative explanations for the super-orbital period., Comment: Accepted for publication in ApJ

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

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

We present new radio observations of the binary neutron star merger GW170817 carried out with the Karl G. Jansky Very large Array (VLA) more than 3\,yrs after the merger. Our combined dataset is derived by co-adding more than $\approx32$\,hours of VLA time on-source, and as such provides the deepest combined observation (rms sensitivity $\approx 0.99\,\mu$Jy) of the GW170817 field obtained to date at 3\,GHz. We find no evidence for a late-time radio re-brightening at a mean epoch of $t\approx 1200$\,d since merger, in contrast to a $\approx 2.1\,\sigma$ excess observed at X-ray wavelengths at the same mean epoch. Our measurements agree with expectations from the post-peak decay of the radio afterglow of the GW170817 structured jet. Using these results, we constrain the parameter space of models that predict a late-time radio re-brightening possibly arising from the high-velocity tail of the GW170817 kilonova ejecta, which would dominate the radio and X-ray emission years after the merger (once the structured jet afterglow fades below detection level). Our results point to a steep energy-speed distribution of the kilonova ejecta (with energy-velocity power law index $\alpha \gtrsim 5$). We suggest possible implications of our radio analysis, when combined with the recent tentative evidence for a late-time re-brightening in the X-rays, and highlight the need for continued radio-to-X-ray monitoring to test different scenarios., Comment: 6 pages, 2 figures. Version published in ApJL (includes minor updates to figures, data and refs)

Here, we present MAXI, SWIFT, NICER, NuSTAR and Chandra observations of the X-ray transient AT2019wey (SRGA J043520.9+552226, SRGE J043523.3+552234). From spectral and timing analyses we classify it as a Galactic low-mass X-ray binary (LMXB) with a black hole (BH) or neutron star (NS) accretor. AT2019wey stayed in the low/hard state (LHS) from 2019 December to 2020 August 21, and the hard-intermediate state (HIMS) from 2020 August 21 to 2020 November. For the first six months of the LHS, AT2019wey had a flux of $\sim 1$ mCrab, and displayed a power-law X-ray spectrum with photon index $\Gamma = 1.8$. From 2020 June to August, it brightened to $\sim 20$ mCrab. Spectral features characteristic of relativistic reflection became prominent. On 2020 August 21, the source left the "hard line" on the rms--intensity diagram, and transitioned from LHS to HIMS. The thermal disk component became comparable to the power-law component. A low-frequency quasi-periodic oscillation (QPO) was observed. The QPO central frequency increased as the spectrum softened. No evidence of pulsation was detected. We are not able to decisively determine the nature of the accretor (BH or NS). However, the BH option is favored by the position of this source on the $\Gamma$--$L_{\rm X}$, $L_{\rm radio}$--$L_{\rm X}$, and $L_{\rm opt}$--$L_{\rm X}$ diagrams. We find the BH candidate XTE J1752-223 to be an analog of AT2019wey. Both systems display outbursts with long plateau phases in the hard states. We conclude by noting the potential of SRG in finding new members of this emerging class of low luminosity and long-duration LMXB outbursts., Comment: Accepted by ApJ

We present the discovery of a luminous X-ray transient, serendipitously detected by Swift's X-ray Telescope (XRT) on 2020 February 5, located in the nucleus of the galaxy SDSS J143359.16+400636.0 at z=0.099 (luminosity distance $D_{\rm L}=456$ Mpc). The transient was observed to reach a peak luminosity of $\sim10^{44}$ erg s$^{-1}$ in the 0.3--10 keV X-ray band, which was $\sim20$ times more than the peak optical/UV luminosity. Optical, UV, and X-ray lightcurves from the Zwicky Transient Facility (ZTF) and Swift show a decline in flux from the source consistent with $t^{-5/3}$, and observations with NuSTAR and Chandra show a soft X-ray spectrum with photon index $\Gamma=2.9\pm0.1$. The X-ray/UV properties are inconsistent with well known AGN properties and have more in common with known X-ray tidal disruption events (TDE), leading us to conclude that it was likely a TDE. The broadband spectral energy distribution (SED) can be described well by a disk blackbody model with an inner disk temperature of $7.3^{+0.3}_{-0.8}\times10^{5}$ K, with a large fraction ($>40$%) of the disk emission up-scattered into the X-ray band. An optical spectrum taken with Keck/LRIS after the X-ray detection reveals LINER line ratios in the host galaxy, suggesting low-level accretion on to the supermassive black hole prior to the event, but no broad lines or other indications of a TDE were seen. The stellar velocity dispersion implies the mass of the supermassive black hole powering the event is log($M_{\rm BH}$/$M_{\odot}$)$=7.41\pm0.41$, and we estimate that at peak the Eddington fraction of this event was $\sim$50%. This likely TDE was not identified by wide-field optical surveys, nor optical spectroscopy, indicating that more events like this would be missed without wide-field UV or X-ray surveys., Comment: Accepted for publication by ApJ. Accepted version now replaces initial submission

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 evaluated the impact of differences in data processing or flux determination methods used by different groups, and attempted to mitigate these differences to provide a more uniform dataset. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.92\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.584\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.168\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., Comment: Accepted by ApJ. Revised X-ray and up-to-date afterglow measurements are available in ASCII format on GitHub (https://github.com/kmooley/GW170817/) and mirrored on tauceti.caltech.edu (http://www.tauceti.caltech.edu/kunal/gw170817/). The online dataset will be continuously updated as new measurements become available

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., Comment: The paper has been updated for the NuSTARDAS release v. 2.0.0 and CALDB version 20200813

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., Comment: In press at ApJ

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., Comment: 9 pages, 3 figures, 1 table. Accepted for publication in ApJ; updated to accepted version with minor revisions

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

We present results from modeling the broadband X-ray spectrum of the Type 1 AGN 2MASX J193013.80+341049.5 using NuSTAR, Swift and archival XMM-Newton observations. We find this source to be highly X-ray obscured, with column densities exceeding 10$^{23}$ cm$^{-2}$ across all epochs of X-ray observations, spanning an 8 year period. However, the source exhibits prominent broad optical emission lines, consistent with an unobscured Type 1 AGN classification. We fit the X-ray spectra with both phenomenological reflection models and physically-motivated torus models to model the X-ray absorption. We examine the spectral energy distribution of this source and investigate some possible scenarios to explain the mismatch between X-ray and optical classifications. We compare the ratio of reddening to X-ray absorbing column density ($E_{B-V}/N_{H}$) and find that 2MASX J193013.80+341049.5 likely has a much lower dust-to-gas ratio relative to the Galactic ISM, suggesting that the Broad Line Region (BLR) itself could provide the source of extra X-ray obscuration, being composed of low-ionization, dust-free gas., Comment: 13 pages, 9 figures. Accepted for publication in ApJ

We present a 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., Comment: Published in A&A. Xspec tables and VR videos of the UXClumpy model can be found at https://github.com/JohannesBuchner/xars/blob/master/doc/README.rst , among other geometries. The X-ray ray tracing code XARS is open source and available at https://github.com/JohannesBuchner/xars . Contributions and collaborations are welcome

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

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., Comment: 15 pages, 7 figures, 5 tables. Revised version

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., Comment: 19 pages, 11 figures, 6 tables. Accepted for publication by ApJ

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., Comment: Astro 2020 Decadal science White Paper

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., Comment: 9 pages, 4 figures, accepted to MNRAS, updated following corrections

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., Comment: Accepted for publication in the Astrophysical Journal

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

We present the 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., Comment: 15 pages, 11 figures, published in the Astrophysics Journal Supplement Series

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

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

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

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., Comment: Author version of the paper published in Nature Astronomy on 26 February 2018. See free-to-view published version here: http://rdcu.be/HQpR

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

We study the geometry of the AGN obscurer in IRAS 09104+4109, an IR-luminous, radio-intermediate FR-I source at $z=0.442$, using infrared data from Spitzer and Herschel, X-ray data from Nustar, Swift, Suzaku, and Chandra, and an optical spectrum from Palomar. The infrared data imply a total rest-frame 1-1000$\mu$m luminosity of $5.5\times10^{46}$erg s$^{-1}$ and require both an AGN torus and starburst model. The AGN torus has an anisotropy-corrected IR luminosity of $4.9\times10^{46}$erg s$^{-1}$, and a viewing angle and half opening angle both of approximately $36$ degrees from pole-on. The starburst has a star formation rate of $(110\pm34)$M$_{\odot}$ yr$^{-1}$ and an age of $<50$Myr. These results are consistent with two epochs of luminous activity in 09104: one approximately $150$Myr ago, and one ongoing. The X-ray data suggest a photon index of $\Gamma \simeq 1.8$ and a line-of-sight column of $N_{\rm H} \simeq 5\times10^{23}$cm$^{-2}$. This argues against a reflection-dominated hard X-ray spectrum, which would have implied a much higher $N_{\rm H}$ and luminosity. The X-ray and infrared data are consistent with a bolometric AGN luminosity of $L_{\rm bol}\sim(0.5-2.5)\times10^{47}$erg s$^{-1}$. The X-ray and infrared data are further consistent with coaligned AGN obscurers in which the line of sight "skims" the torus. This is also consistent with the optical spectra, which show both coronal iron lines and broad lines in polarized but not direct light. Combining constraints from the X-ray, optical, and infrared data suggests that the AGN obscurer is within a vertical height of $20$pc, and a radius of $125$pc, of the nucleus., Comment: Accepted to ApJ

We 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}$ $

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., Comment: Accepted for publication by ApJ