Search

'Bare' active galactic nuclei (AGN) are a subclass of Type 1 AGN that show little or no intrinsic absorption. They offer an unobscured view of the central regions of the AGN and therefore serve as ideal targets to study the relativistic reflection features originating from the innermost regions of the accretion disc. We present a detailed broadband spectral analysis ($0.3 - 70$ keV) of one of the most luminous bare AGN in the local universe, RBS 1124 ($z= 0.208$) using a new, co-ordinated high signal-to-noise observation obtained by $\textit{XMM-Newton}$ and $\textit{NuSTAR}$. The source exhibits a power-law continuum with $\Gamma \sim$ 1.8 along with a soft excess below 2 keV, a weak neutral iron line and curvature at high energies ($\sim 30$ keV). The broadband spectrum, including the soft excess and the high-energy continuum, is well fit by the relativistic reflection model when the accretion disc is allowed to have densities of log$(n_{\rm e}$/cm$^{-3}$) $\gtrsim 19.2$. Our analysis therefore suggests that when high-density effects are considered, relativistic reflection remains a viable explanation for the soft excess., Comment: 14 pages, 9 figures, accepted for publication in MNRAS

Galaxy clusters show large-scale azimuthal X-ray surface brightness fluctuations known as cold fronts. These are overdense (average density jumps $\sim 30\%$ or post-jump density $\sim 130\%$) and have milder discontinuity in pressure. Cold fronts are argued to originate due to sloshing driven by sub-halo passage at close proximity to the cluster center. While this is a viable source of large-scale perturbations, the physical mechanisms that can sustain such density structures (of specific geometry) are not clear. In this work, we explore whether long wavelength thermal instability is an explanation for cold front formation in a cluster core which is perturbed by sub-halos or AGN activity. Using global linear perturbation analysis, we show that internal gravity waves (thermally unstable) can form large-scale three-dimensional spiral structures, akin to observed cold fronts. We explore if the presence of magnetic field (along spherical $\hat{\phi}$) may support such structures (by suppressing small scale Kelvin-Helmholtz modes) or disrupt them (by promoting additional thermal instability). We find that latter happens at shorter wavelengths and only at frequencies above the characteristic buoyancy or Brunt V\"ais\"al\"a frequency ($>N_{\rm BV}$). Our work implies, firstly, that large-scale spirals may be formed and sustained over a long timescale ($>N^{-1}_{\rm BV}$) even in presence of aligned magnetic fields that is otherwise supportive against mixing at the interface. Secondly, short-wavelength (but relatively longer along the field) unstable compressive modes may form within or in the vicinity of such spirals. The instability is an overstable slow wave, and grows in 2D at timescales $\gtrsim 2-3$ times longer than the spiral growth timescale (via thermal instability). Thus we claim that this instability cannot destroy the large scale coherence., Comment: 14 pages, 8 figures in main content and 3 figures in Appendix, to be submitted to MNRAS. Comments are welcome

It has been recently discovered that a few super-Eddington sources undergoing black hole super-Eddington accretion exhibit X-ray reflection signatures. In such new systems, one expects that the coronal X-ray emissions are mainly reflected by optically thick super-Eddington winds instead of thin disks. In this paper, we conduct a series of general relativistic ray-tracing and Monte Carlo radiative transfer simulations to model the X-ray reflection signatures, especially the characteristic Fe K$\alpha$ line, produced from super-Eddington accretion flows. In particular, we allow the photons emitted by a lamppost corona to be reflected multiple times in a cone-like funnel surrounded by fast winds. We find that the Fe K$\alpha$ line profile most sensitively depends on the wind kinematics, while its exact shape also depends on the funnel open angle and corona height. Furthermore, very interestingly, we find that the Fe K$\alpha$ line can have a prominent double-peak profile in certain parameter spaces even with a face-on orientation. Moreover, we compare the Fe K$\alpha$ line profiles produced from super-Eddington and thin disks and show that such lines can provide important insights into the understanding of black hole systems undergoing super-Eddington accretion., Comment: 23 pages, 21 figures, 2 tables. Comments are welcome

The Advanced X-ray Imaging Satellite (AXIS) is a Probe-class concept that will build on the legacy of the Chandra X-ray Observatory by providing low-background, arcsecond-resolution imaging in the 0.3-10 keV band across a 450 arcminute$^2$ field of view, with an order of magnitude improvement in sensitivity. AXIS utilizes breakthroughs in the construction of lightweight segmented X-ray optics using single-crystal silicon, and developments in the fabrication of large-format, small-pixel, high readout rate CCD detectors with good spectral resolution, allowing a robust and cost-effective design. Further, AXIS will be responsive to target-of-opportunity alerts and, with onboard transient detection, will be a powerful facility for studying the time-varying X-ray universe, following on from the legacy of the Neil Gehrels (Swift) X-ray observatory that revolutionized studies of the transient X-ray Universe. In this paper, we present an overview of AXIS, highlighting the prime science objectives driving the AXIS concept and how the observatory design will achieve these objectives., Comment: Published in Proceedings of SPIE Optics & Photonics 2023, San Diego

We present results from a pair of high resolution, long timescale ($\sim10^5 GM/c^3$), global, three dimensional magnetohydrodynamical accretion disk simulations with differing initial magnetic plasma $\beta$ in order to study the effects of initial toroidal field strength on production of large-scale poloidal field. We initialize our disks in approximate equilibrium with purely toroidal magnetic fields of strength $\beta_0=5$ and $\beta_0=200$. We also perform a limited resolution study. We find that simulations of differing field strength diverge early in their evolution and remain distinct over the time studied, indicating that initial magnetic conditions leave a persistent imprint in our simulations. Neither simulation enters the Magnetically Arrested Disk (MAD) regime. Both simulations are able to produce poloidal fields from initially-toroidal fields, with the $\beta_0=5$ simulation evolving clear signs of a large-scale poloidal field. We make a cautionary note that computational artifacts in the form of large-scale vortices may be introduced in the combination of initially-weak field and disk-internal mesh refinement boundaries, as evidenced by the production of an $m=1$ mode overdensity in the weak field simulation. Our results demonstrate that the initial toroidal field strength plays a vital role in simulated disk evolution for the models studied., Comment: 17 pages, 11 figures, accepted to ApJ

Accretion occurs across a large range of scales and physical regimes. Despite this diversity in the physics, the observed properties show remarkably similarity. The theory of propagating fluctuations, in which broad-band variability within an accretion disc travel inwards and combine, has long been used to explain these phenomena. Recent numerical work has expanded on the extensive analytical literature but has been restricted to using the 1D diffusion equation for modelling the disc behaviour. In this work we present a novel numerical approach for 2D (vertically integrated), stochastically driven {\alpha}-disc simulations, generalising existing 1D models. We find that the theory of propagating fluctuations translates well to 2D. However, the presence of epicyclic motion in 2D (which cannot be captured within the diffusion equation) is shown to have an important impact on local disc dynamics. Additionally, there are suggestions that for sufficiently thin discs the log-normality of the light-curves changes. As in previous work, we find that the break frequency in the luminosity power spectrum is strongly dependent on the driving timescale of the stochastic perturbations within the disc, providing a possible observational signature for probing the magnetorotational instability (MRI) dynamo. We also find that thinner discs are significantly less variable than thicker ones, providing a compelling explanation for the greater variability seen in the hard state vs the soft state of X-ray binaries. Finally, we consider the wide-ranging applications of our numerical model for use in other simulations., Comment: 28 pages, 22 figures, 10 tables (including 3 appendices). Resubmitted to MNRAS following corrections

We report a detailed analysis of the XMM-Newton spectra of six Narrow-Line Seyfert 1 (NLS1) galaxies at redshift z = 0.35--0.92. Compared with the NLS1s at lower redshift in the previously most-studied sample, these NLS1s have larger black hole (BH) masses ($\log\,M_\text{BH}>7.5$) with similar or even lower Eddington ratios. Our extended XMM-Newton sample of NLS1s shows strong soft X-ray excess emission below 2 keV. The quantified soft excess strength does not show an obvious discrepancy from previous studies of the lower-redshift NLS1s. The systematic effect in the measurement of the Eddington ratio mainly lies in the bolometric correction factor. We also tentatively fit the spectra assuming two more physical models for the soft excess: warm Comptonization and relativistic reflection from the inner accretion disk. In the first scenario, we confirm the ubiquity of a warm and optically thick corona. The behavior of a single source can be better explained by relativistic reflection, although we cannot distinguish which model is a more favorable explanation for the soft excess based on the best-fit statistics., Comment: 13 pages, 10 figures

We discuss our recent constraints on the coupling of Very-Light Axion-Like Particles (of masses $<$$ 10^{-12} \ \mathrm{eV}$) to electromagnetism from $Chandra$ observations of the cluster-hosted Active Galactic Nuclei (AGN) H1821+643 and NGC1275. In both cases, the inferred high-quality AGN spectra excluded all photon-ALP couplings $g_\mathrm{a\gamma} > (6.3 - 8.0) \times 10^{-13} \ {\mathrm{GeV}}^{-1}$ at the $99.7\%$ level, respectively, based on the non-detection of spectral distortions attributed to photon-ALP inter-conversion along the cluster line-of-sight. Finally, we present the prospects of tightening current bounds on such ALPs by up to a factor of 10 with next-generation X-ray observatories such as $Athena$, $AXIS$ and $LEM$ given their improved spectral and spatial resolution and collecting area compared to current missions., Comment: 4 pages, 1 figure. Accepted for publication to Memorie della SAIt for the Proceedings of the European Astronomical Society 2022 (EAS 2022) Annual Meeting Symposium S3 "The Dark matter multi-messenger challenge"

We use numerical hydrodynamics simulations to study line driven winds launched from an accreting alpha-disc. Building on previous work where the driving radiation field is static, we compute a time-dependent radiation flux from the local, variable accretion rate of the disc. We find that prior to the establishment of a steady state in the disc, variations of ~ 15% in disc luminosity correlate with variations of ~ 2-3 in the mass flux of the wind. After a steady state is reached, when luminosity variations drop to ~ 3%, these correlations vanish as the variability in the mass flux is dominated by the intrinsic variability of the winds. This is especially evident in lower luminosity runs where intrinsic variability is higher due to a greater prevalence of failed winds. The changing mass flux occurs primarily due to the formation of clumps and voids near the disc atmosphere that propagate out into the low velocity part of the flow, a process that can be influenced by local variations in disc intensity. By computing the normalised standard deviation of the mass outflow, we show that the impact of luminosity variations on mass outflow is more visible at higher luminosity. However, the absolute change in mass outflow due to luminosity increases is larger for lower luminosity models due to the luminosity-mass flux scaling relation becoming steeper. We further discuss implications for CVs and AGN and observational prospects., Comment: 12 pages, 7 figures, accepted for publication in MNRAS

Axion-Like Particles (ALPs) are well-motivated extensions of the Standard Model of Particle Physics and a generic prediction of some string theories. X-ray observations of bright Active Galactic Nuclei (AGN) hosted by rich clusters of galaxies are excellent probes of very-light ALPs, with masses $\mathrm{log}(m_\mathrm{a}/\mathrm{eV}) < -12.0$. We evaluate the potential of future X-ray observatories, particularly $Athena$ and the proposed $AXIS$, to constrain ALPs via observations of cluster-hosted AGN, taking NGC 1275 in the Perseus cluster as our exemplar. Assuming perfect knowledge of instrument calibration, we show that a modest exposure (200-ks) of NGC 1275 by $Athena$ permits us to exclude all photon-ALP couplings $g_\mathrm{a\gamma} > 6.3 \times 10^{-14} \ {\mathrm{GeV}}^{-1}$ at the 95% level, as previously shown by $Conlon \ et \ al. \ (2018)$, representing a factor of 10 improvement over current limits. We then proceed to assess the impact of realistic calibration uncertainties on the $Athena$ projection by applying a standard $Cash$ likelihood procedure, showing the projected constraints on $g_\mathrm{a\gamma}$ weaken by a factor of 10 (back to the current most sensitive constraints). However, we show how the use of a deep neural network can disentangle the energy-dependent features induced by instrumental miscalibration and those induced by photon-ALP mixing, allowing us to recover most of the sensitivity to the ALP physics. In our explicit demonstration, the machine learning applied allows us to exclude $g_\mathrm{a\gamma} > 2.0 \times 10^{-13} \ {\mathrm{GeV}}^{-1}$, complementing the projected constraints of next-generation ALP dark matter birefringent cavity searches for very-light ALPs. Finally, we show that a 200-ks $AXIS$/on-axis observation of NGC 1275 will tighten the current best constraints on very-light ALPs by a factor of 3., Comment: 22 pages, 10 figures, 3 tables. Accepted for publication at the Astrophysical Journal. Comments welcome

Ark 564 is an extreme high-Eddington Narrow-line Seyfert 1 galaxy, known for being one of the brightest, most rapidly variable soft X-ray AGN, and for having one of the lowest temperature coronae. Here we present a 410-ks NuSTAR observation and two 115-ks XMM-Newton observations of this unique source, which reveal a very strong, relativistically broadened iron line. We compute the Fourier-resolved time lags by first using Gaussian processes to interpolate the NuSTAR gaps, implementing the first employment of multi-task learning for application in AGN timing. By fitting simultaneously the time lags and the flux spectra with the relativistic reverberation model RELTRANS, we constrain the mass at $2.3^{+2.6}_{-1.3} \times 10^6M_\odot$, although additional components are required to describe the prominent soft excess in this source. These results motivate future combinations of machine learning, Fourier-resolved timing, and the development of reverberation models., Comment: 19 pages, 9 figures. Accepted for publication in The Astrophysical Journal

We present the results from a coordinated XMM-Newton $+$ NuSTAR observation of the Seyfert 1 Galaxy ESO 511$-$G030. With this joint monitoring programme, we conduct a detailed variability and spectral analysis. The source remained in a low flux and very stable state throughout the observation period, although there are slight fluctuations of flux over long timescales. The broadband (0.3-78~keV) spectrum shows the presence of a power-law continuum with a soft excess below 2~keV, a relatively narrow iron K$\alpha$ emission ($\sim$6.4~keV), and an obvious cutoff at high energies. We find that the soft excess can be modeled by two different possible scenarios: a warm ($kT_{\rm e} \sim$ 0.19~keV) and optically thick ($\tau - 18\sim25$) Comptonizing corona or a relativistic reflection from a high-density ($\log [n_{\rm e}/{\rm cm}^{-3}]=17.1 \sim 18.5$) inner disc. All models require a low temperature ($kT_{\rm e} \sim$ 13~keV) for the hot corona., Comment: 18 pages, 15 figures. v2: refereed version

We present a study of the narrow Fe K$\alpha$ line in seven bright, nearby AGN that have been observed extensively with the Chandra High Energy Transmission Grating (HETG). The HETG data reveal a wider Fe K$\alpha$ line in the first order spectrum than in the second and third order spectra, which we interpret as the result of spatially extended Fe K$\alpha$ emission. We utilize these differences in narrow Fe K$\alpha$ line widths in the multi-order Chandra HETG spectra to determine the spatial extent and intrinsic velocity width of the emitting material in each object. We find that there is modest evidence for spatially extended emission in each object, corresponding to extension of $r\sim5-100$ pc. These distances are significantly larger than those inferred from velocity widths assuming gravitational motions, which give $r\sim0.01-1$ pc. This implies that either the gas is emitting at a range of radii, with smaller radii dominating the velocity width and larger radii dominating the spatial extent, or that the gas is exhibiting non-gravitational motions, which we suggest would be outflows due to slight excess redshift in the line and velocities that exceed the freefall velocity. We also use the spatial extent information to estimate the mass of the emitting gas by counting fluorescing iron atoms, finding masses on the order of $M_\mathrm{gas}\sim10^5-10^8\,M_\odot$. Future work with observatories like XRISM will be able to extend this study to a larger number of AGN and decrease uncertainties that arise due to the low signal-to-noise of the higher order HETG data., Comment: 22 pages, 4 figures, accepted for publication in ApJ

We study the effect of returning radiation on the shape of the X-ray reflection spectrum in the case of thin accretion disks. We show that the returning radiation mainly influences the observed reflection spectrum for a large black hole spin (a > 0.9) and a compact primary source of radiation close to the black hole at height h < 5 $r_\mathrm{g}$, and that it dominates the reflected flux for extreme values of spin and compactness. The main effect of the returning radiation is to increase the irradiating flux on to the outer parts of the accretion disk, leading to stronger reflection and a flatter overall emissivity profile. By analyzing simulated observations we show that neglecting returning radiation in existing studies of reflection dominated sources has likely resulted in overestimating the height of the corona above the black hole. An updated version of the publicly available relxill suite of relativistic reflection models which includes returning radiation is also presented., Comment: accepted for publication by MNRAS; the new RELXILL model can be downloaded at https://www.sternwarte.uni-erlangen.de/research/relxill/ including a new and detailed documentation (https://www.sternwarte.uni-erlangen.de/research/relxill/relxill_docu.pdf)

We present an analysis of deep $Chandra$ Low-Energy and High-Energy Transmission Grating archival observations of the extraordinarily luminous radio-quiet quasar H1821+643, hosted by a rich and massive cool-core cluster at redshift $z=0.3$. These datasets provide high-resolution spectra of the AGN at two epochs, free from contamination by the intracluster medium and from the effects of photon pile-up, providing a sensitive probe of the iron-$K$ band. At both epochs, the spectrum is well described by a power-law continuum plus X-ray reflection from both the inner accretion disc and cold, slowly-moving distant matter. Adopting this framework, we proceed to examine the properties of the inner disc and the black hole spin. Using Markov chain Monte Carlo (MCMC) methods, we combine constraints from the two epochs assuming that the black hole spin, inner disc inclination, and inner disc iron abundance are invariant. The black hole spin is found to be modest, with a 90$\%$ credible range of ${a}^{*}=0.62^{+0.22}_{-0.37}$; and, with a mass $M_\mathrm{BH}$ in the range $\log (M_\mathrm{BH}/M_\odot)\sim 9.2-10.5$, this is the most massive black hole candidate for which a well-defined spin constraint has yet been obtained. The modest spin of this black hole supports previous suggestions that the most massive black holes may grow via incoherent or chaotic accretion and/or SMBH-SMBH mergers., Comment: 13 pages, 8 figures, 3 tables. Accepted for publication in MNRAS. Amended to include two missing references

We present a detailed analysis of the XMM-Newton observations of five narrow-line Seyfert 1 galaxies (NLS1s). They all show very soft continuum emission in the X-ray band with a photon index of $\Gamma\gtrsim 2.5$. Therefore, they are referred to as `ultra-soft' NLS1s in this paper. By modeling their optical/UV-X-ray spectral energy distribution (SED) with a reflection-based model, we find indications that the disc surface in these ultra-soft NLS1s is in a higher ionisation state than other typical Seyfert 1 AGN. Our best-fit SED models suggest that these five ultra-soft NLS1s have an Eddington ratio of $\lambda_{\rm Edd}=1-20$ assuming available black hole mass measurements. In addition, our models infer that a significant fraction of the disc energy in these ultra-soft NLS1s is radiated away in the form of non-thermal emission instead of the thermal emission from the disc. Due to their extreme properties, X-ray observations of these sources in the iron band are particularly challenging. Future observations, e.g. from Athena, will enable us to have a clearer view of the spectral shape in the iron band and thus distinguish the reflection model from other interpretations of their broadband spectra., Comment: 16 pages, 11 figures; published in MNRAS, Volume 498, Issue 3, November 2020. The Reflionx model used in this work is available for download at https://www.michaelparker.space/reflionx-models

We present a re-analysis of the XMM-Newton and NuSTAR observing campaign for the well-studied, X-ray-bright AGN MCG-06-30-15. In particular, we consider a disc model with finite thickness. By fitting the disc reflection spectra in the data, we obtain a black hole spin of 0.87--0.99 (90\% confidence range) after taking the thickness of the disc into consideration. Spectral models with a grid of mass accretion rate from 0 to $30\%\dot{M}_{\rm Edd}$ are calculated for MCG-06-30-15. This result is obtained by considering a free disc reflection fraction parameter $f_{\rm refl}$ and is consistent with previous measurements based on razor-thin disc models. Besides, an isotropic, point-like geometry, i.e. the `lamppost' geometry, is assumed for the corona in our model. We find that such a geometry overestimates $f_{\rm refl}$ in the data. Therefore, thin disc models with consistent `lamppost' values of $f_{\rm refl}$ provide a worse fit than ones with a free $f_{\rm refl}$ parameter. We discuss possible reasons for the discrepancy between the observed and theoretical values of $f_{\rm refl}$ at the end of the paper. Modifications for the over-simplified lamppost model might be needed when the thickness of the thin disc is considered in future work., Comment: 15 pages, 14 figures, accepted for publication in MNRAS

Previously, we modelled the X-ray spectra of the narrow-line Seyfert 1 galaxy IRAS 13224$-$3809 using a disc reflection model with a fixed electron density of $10^{15}$ cm$^{-3}$. An additional blackbody component was required to fit the soft X-ray excess below 2 keV. In this work, we analyse simultaneously five flux-resolved XMM-Newton spectra of this source comprising data collected over 2 Ms. A disc reflection model with an electron density of $n_{\rm e}\approx10^{20}$ cm$^{-3}$ and an iron abundance of $Z_{\rm Fe}=3.2\pm0.5Z_{\odot}$ is used to fit the broad-band spectra of this source. No additional component is required to fit the soft excess. Our best-fit model provides consistent measurements of black hole spin and disc inclination angle as in previous models where a low disc density was assumed. In the end, we calculate the average illumination distance between the corona and the reflection region in the disc of IRAS 13224$-$3809 based on best-fit density and ionisation parameters, which changes from 0.43$\sqrt{f_{\rm AD}/f_{\rm INF}}$ $r_{\rm g}$ in the lowest flux state to 1.71$\sqrt{f_{\rm AD}/f_{\rm INF}}$ $r_{\rm g}$ in the highest flux state assuming a black hole mass of $2\times10^{6}M_{\odot}$. $f_{\rm AD}/f_{\rm INF}$ is the ratio between the flux of the coronal emission that reaches the accretion disc and infinity. This ratio depends on the geometry of the coronal region in IRAS 13224$-$3809. So we only discuss its value based on the simple `lamp-post' model, although detailed modelling of the disc emissivity profile of IRAS 13224$-$3809 is required in future to reveal the exact geometry of the corona., Comment: 17 pages, 15 figures including appendices; accepted for publication in MNRAS

Axion-like particles (ALPs) are a well-motivated extension to the standard model of particle physics, and X-ray observations of cluster-hosted AGN currently place the most stringent constraints on the ALP coupling to electromagnetism, $g_{a \gamma}$, for very light ALPs ($m_a\lesssim10^{-11}$ eV). We revisit limits obtained by Reynolds et al. (2020) using Chandra X-ray grating spectroscopy of NGC 1275, the central AGN in the Perseus cluster, examining the impact of the X-ray spectral model and magnetic field model. We also present a new publicly available code, ALPro, which we use to solve the ALP propagation problem. We discuss evidence for turbulent magnetic fields in Perseus and show that it can be important to resolve the magnetic field structure on scales below the coherence length. We re-analyse the NGC 1275 X-ray spectra using an improved data reduction and baseline spectral model. We find the limits are insensitive to whether a partially covering absorber is used in the fits. At low $m_a$ ($m_a\lesssim10^{-13}$ eV), we find marginally weaker limits on $g_{a \gamma}$ (by $0.1-0.3$ dex) with different magnetic field models, compared to Model B from Reynolds et al. (2020). A Gaussian random field (GRF) model designed to mimic $\sim50$ kpc scale coherent structures also results in only slightly weaker limits. We conclude that the existing Model B limits are robust assuming that $\beta_{\rm pl}\approx100$, and are insensitive to whether cell-based or GRF methods are used. However, astrophysical uncertainties regarding the strength and structure of cluster magnetic fields persist, motivating high sensitivity RM observations and tighter constraints on the radial profile of $\beta_{\rm pl}$., Comment: 17 pages, 9 figures, 2 tables. Accepted for publication in The Astrophysical Journal. ALPro code available at https://github.com/jhmatthews/alpro with accompanying documentation, as well as notebooks and data for reproducing the main results

Many recent works on the observed galaxy clusters in the X-rays highlight broadly two classes of exclusive energy carriers - sound waves and turbulence. In order to understand this dichotomy, we design an idealized three-dimensional hydrodynamic simulation of a cluster, to assess which of these carriers can dissipate energy in and around the core ($\gtrsim 100$ kpc) . Specifically, we explore how gentle (long-duration outbursts) and intermediate (shorter duration outbursts) feedback modes can function efficiently mediated by compressible (sound waves) and incompressible (g-modes/instabilities/turbulence) disturbances. Since g-modes are confined tightly to the central core, we attempt to maximise the flux of fast sound waves to distribute the feedback energy over a large distance. We find that the contribution to heat dissipation from sound and turbulence varies on the basis of the aforementioned feedback modes, namely: turbulence contributes relatively more than sound in the slow-piston regime and vice versa for the intermediate regime. For the first time in a 3D simulation, we show that up to $\lesssim 20\%$ (in some directions) of the injected power can be carried away by sound flux in the intermediate feedback but it reduces to $\lesssim 10 \%$ (in some directions) in the slow-piston regime. Lastly, we find that sound waves can be elusive if we deduce the equation-of-state (isobaric/isentropic) of the fluctuations from X-ray observations., Comment: 25 pages, 14 figures, 1 table (+2 figures in appendices). Accepted in MNRAS

Axion-Like Particles (ALPs) are predicted by several Beyond the Standard Model theories, in particular, string theory. In the presence of an external magnetic field perpendicular to the direction of propagation, ALPs can couple to photons. Therefore, if an X-ray source is viewed through a magnetised plasma, such as a luminous quasar in a galaxy cluster, we may expect spectral distortions that are well described by photon-ALP oscillations. We present a $571 \ \mathrm{ks}$ combined High and Low Energy Transmission Grating (HETG/LETG) Chandra observation of the powerful radio-quiet quasar H1821+643, hosted by a cool-core cluster at redshift $0.3$. The spectrum is well described by a double power-law continuum and broad$+$narrow iron line emission typical of type-1 Active Galactic Nuclei (AGN), with remaining spectral features $< 2.5\%$. Using a cell-based approach to describe the turbulent cluster magnetic field, we compare our spectrum with photon-ALP mixing curves for 500 field realisations assuming that the thermal-to-magnetic pressure ratio remains constant up to the virial radius. At $99.7\%$ credibility and taking $\beta = 100$, we exclude all couplings $g_\mathrm{a\gamma} > 6.3 \times 10^{-13} \ {\mathrm{GeV}}^{-1}$ for most ALP masses $< 10^{-12} \ \mathrm{eV}$. Our results are moderately more sensitive to constraining ALPs than the best previous result from Chandra observations of the Perseus cluster, albeit with a less constrained field model. We reflect on the promising future of ALP studies with bright AGN embedded in rich clusters, especially with the upcoming Athena mission., Comment: 14 pages, 7 figures. Accepted for publication in MNRAS on 24th November 2021. Comments welcome

NGC1275 is the Brightest Cluster Galaxy (BCG) in the Perseus cluster and hosts the active galactic nucleus (AGN) that is heating the central 100\,kpc of the intracluster medium (ICM) atmosphere via a regulated feedback loop. Here we use a deep 490ks Cycle-19 Chandra High-Energy Transmission Grating (HETG) observation of NGC1275 to study the anatomy of this AGN. The X-ray continuum is adequately described by an unabsorbed power-law with photon index $\Gamma\approx 1.9$, creating strong tension with the detected column of molecular gas seen via HCN and HCO$^+$ line absorption against the parsec-scale core/jet. This tension is resolved if we permit a composite X-ray source; allowing a column of $N_H\sim 8\times 10^{22}\,{\rm cm}^{-2}$ to cover $\sim 15$% of the X-ray emitter does produce a significant improvement in the statistical quality of the spectral fit. We suggest that the dominant unabsorbed component corresponds to the accretion disk corona, and the sub-dominant X-ray component is the jet working surface and/or jet cocoon that is expanding into clumpy molecular gas. We suggest that this may be a common occurence in BCG-AGN. We conduct a search for photoionized absorbers/winds and fail to detect such a component, ruling out columns and ionization parameters often seen in many other Seyfert galaxies. We detect the 6.4keV iron-K$\alpha$ fluorescence line seen previously by XMM-Newton and Hitomi. We describe an analysis methodology which combines dispersive HETG spectra, non-dispersive microcalorimeter spectra, and sensitive XMM-Newton/EPIC spectra in order to constrain (sub)arcsec-scale extensions of the iron-K$\alpha$ emission region., Comment: 13 pages, accepted for publication in MNRAS

Across a large range of scales, accreting sources show remarkably similar patterns of variability, most notably the log-normality of the luminosity distribution and the linear root-mean square (rms)-flux relationship. These results are often explained using the theory of propagating fluctuations in which fluctuations in the viscosity create perturbations in the accretion rate at all radii, propagate inwards and combine multiplicatively. While this idea has been extensively studied analytically in a linear regime, there has been relatively little numerical work investigating the non-linear behaviour. In this paper, we present a suite of stochastically driven 1-d $\alpha$-disc simulations, exploring the behaviour of these discs. We find that the eponymous propagating fluctuations are present in all simulations across a wide range of model parameters, in contradiction to previous work. Of the model parameters, we find by far the most important to be the timescale on which the viscosity fluctuations occur. Physically, this timescale will depend on the underlying physical mechanism, thought to be the magnetorotational instability (MRI). We find a close relationship between this fluctuation timescale and the break frequency in the power spectral density (PSD) of the luminosity, a fact which could allow observational probes of the behaviour of the MRI dynamo. We report a fitting formula for the break frequency as a function of the fluctuation timescale, the disc thickness and the mass of the central object., Comment: 18 pages, 14 figures, accepted for publication in MNRAS

We present a spectroscopic and imaging study of an abnormal active galactic nucleus (AGN), 2MASX J00423991+3017515. This AGN is newly identified in the hard X-rays by the Swift BAT All-Sky survey and found in an edge-on disk galaxy interacting with a nearby companion. Here, we analyze the first optical spectra obtained for this system (taken in 2011 and 2016), high-resolution imaging taken with the Hubble Space Telescope and Chandra X-ray Observatory, and 1" imaging with the Very Large Array. Two unique properties are revealed: the peaks of the broad Balmer emission lines (associated with gas orbiting very near the supermassive black hole) are blue shifted from the corresponding narrow line emission and host galaxy absorption by 1540 km/s, and the AGN is spatially displaced from the apparent center of its host galaxy by 3.8 kpc. We explore several scenarios to explain these features, along with other anomalies, and propose that 2MASX J00423991+3017515 may be an AGN with an unusually strong wind residing in a uniquely configured major merger, or that it is an AGN recoiling from either a gravitational "slingshot" in a three-body interaction or from a kick due to the asymmetric emission of gravitational waves following the coalescence of two progenitor supermassive black holes., Comment: 17 pages, 10 figures, 4 tables, accepted for publication in MNRAS

Outflows from active galactic nuclei (AGN) are one of the fundamental mechanisms by which the central supermassive black hole interacts with its host galaxy. Detected in $\ge 50\%$ of nearby AGN, these outflows have been found to carry kinetic energy that is a significant fraction of AGN power, and thereby give negative feedback to their host galaxies. To understand the physical processes that regulate them, it is important to have a robust estimate of their physical and dynamical parameters. In this review we summarize our current understanding on the physics of the ionized outflows detected in absorption in the UV and X-ray wavelength bands. We discuss the most relevant observations and our current knowledge and uncertainties in the measurements of the outflow parameters. We also discuss their origin and acceleration mechanisms. The commissioning and concept studies of large telescope missions with high resolution spectrographs in UV/optical and X-rays along with rapid advancements in simulations offer great promise for discoveries in this field over the next decade., Comment: An invited Review Article published in Nature Astronomy: https://www.nature.com/articles/s41550-020-01255-2

KUG 1141+371 is a Seyfert 1 galaxy that shows a simultaneous flux increase in the optical and UV bands in the past decade. For instance, the latest Swift observation in 2019 shows that the UVW2 flux of the AGN in KUG 1141+371 has increased by over one order of magnitude since 2009. Meanwhile, the soft X-ray flux of KUG 1141+371 also shows a steady increase by one order of magnitude since 2007. The significant multi-wavelength luminosity change is likely due to a boost of mass accretion rate from approximately 0.6% of the Eddington limit to 3.2%, assuming a black hole mass of $10^{8}M_{\odot}$. In this work, we conduct detailed multi-epoch X-ray spectral analysis focusing on the variability of the X-ray continuum emission and the puzzling soft excess emission. In addition, our SED models also suggest a simultaneous increase of disc temperature and a decreasing inner disc radius along with the increasing accretion rate. Finally, we discuss possible connection between KUG 1141+371 and black hole transients in outburst., Comment: 18 pages, 14 figures, accepted by MNRAS

The spin of a black hole is an important quantity to study, providing a window into the processes by which a black hole was born and grew. Further, spin can be a potent energy source for powering relativistic jets and energetic particle acceleration. In this review, I describe the techniques currently used to detect and measure the spins of black holes. It is shown that: (1) Two well understood techniques, X-ray reflection spectroscopy and thermal continuum fitting, can be used to measure the spins of black holes that are accreting at moderate rates. There is a rich set of other electromagnetic techniques allowing us to extend spin measurements to lower accretion rates. (2) Many accreting supermassive black holes are found to be rapidly-spinning, although a population of more slowly spinning black holes emerges at masses above $M>3\times 10^7\,M_\odot$ as expected from recent structure formation models. (3) Many accreting stellar-mass black holes in X-ray binary systems are rapidly spinning and must have been born in this state. (4) The advent of gravitational wave astronomy has enabled the detection of spin effects in merging binary black holes. Most of the pre-merger black holes are found to be slowly spinning, a notable exception being an object that may itself be a merger product. (5) The stark difference in spins between the black hole X-ray binary and the binary black hole populations shows that there is a diversity of formation mechanisms. Given the array of new electromagnetic and gravitational wave capabilities currently being planned, the future of black hole spin studies is bright., Comment: 42 pages, 6 figures. Invited review to be published in Annual Reviews of Astronomy and Astrophysics

We study Lense-Thirring precession of inviscid and viscous misaligned $\alpha-$discs around a black hole using a gravitomagnetic term in the momentum equation. For weak misalignments, $i \lesssim 10^{\circ}$, the discs behave like rigid bodies, undergoing the full suite of classical harmonic oscillator dynamics including, weak and critically damped motion (due to viscosity), precession (due to Lense-Thirring torque) and nutation (due to apsidal precession). For strong misalignments, $i \gtrsim 30^{\circ}$, we find sufficiently thin, $h/r \lesssim 0.05$ discs break, form a gap and the inner and outer sub-discs evolve quasi independently apart from slow mass transfer. Assuming the sound speed sets the communication speed of warps in the disc, we can estimate the breaking radius by requiring that the inner sub-disc precesses like a rigid body. We explicitly show for the first time using a grid code that an Einstein potential is needed to reproduce the analytic properties of the inner disc edge and find disc breaking. At large inclination angles we find multiple disc breaks, consistent with recent GRMHD simulations of highly inclined discs. Our results suggest that the inclusion of a gravitomagnetic term and appropriate pseudo-Newtonian potential captures the important quantitative features of misaligned discs., Comment: 13 pages, 10 figures

We present a comprehensive search for the 3.5 keV line, using $\sim$51 Ms of archival Chandra observations peering through the Milky Way's Dark Matter Halo from across the entirety of the sky, gathered via the Chandra Source Catalog Release 2.0. We consider the data's radial distribution, organizing observations into four data subsets based on angular distance from the Galactic Center. All data is modeled using both background-subtracted and background-modeled approaches to account for the particle instrument background, demonstrating statistical limitations of the currently-available $\sim$1 Ms of particle background data. A non-detection is reported in the total data set, allowing us to set an upper-limit on 3.5 keV line flux and constrain the sterile neutrino dark matter mixing angle. The upper-limit on sin$^2$(2$\theta$) is $2.58 \times 10^{-11}$ (though systematic uncertainty may increase this by a factor of $\sim$2), corresponding to the upper-limit on 3.5 keV line flux of $2.34 \times 10^{-7}$ ph s$^{-1}$ cm$^{-2}$. These limits show consistency with recent constraints and several prior detections. Non-detections are reported in all radial data subsets, allowing us to constrain the spatial profile of 3.5 keV line intensity, which does not conclusively differ from Navarro-Frenk-White predictions. Thus, while offering heavy constraints, we do not entirely rule out the sterile neutrino dark matter scenario or the more general decaying dark matter hypothesis for the 3.5 keV line. We have also used the non-detection of any unidentified emission lines across our continuum to further constrain the sterile neutrino parameter space., Comment: 37 pages, 27 figures, accepted by ApJ

We study how different opacity-temperature scalings affect the dynamical evolution of irradiated gas clouds using time-dependent, radiation-hydrodynamics (rad-HD) simulations. When clouds are optically thick, the bright side heats up and expands, accelerating the cloud via the rocket effect. Clouds that become more optically thick as they heat accelerate $\sim 35\%$ faster than clouds that become optically thin. An enhancement of $\sim 85\%$ in the acceleration can be achieved by having a broken powerlaw opacity profile, which allows the evaporating gas driving the cloud to become optically thin and not attenuate the driving radiation flux. We find that up to $\sim 2\%$ of incident radiation is re-emitted by accelerating clouds, which we estimate as the contribution of a single accelerating cloud to an emission or absorption line. Re-emission is suppressed by "bumps" in the opacity-temperature relation since these decrease the opacity of the hot, evaporating gas, primarily responsible for the re-radiation. If clouds are optically thin, they heat nearly uniformly, expand and form shocks. This triggers the Richtmyer-Meshkov instability, leading to cloud disruption and dissipation on thermal time-scales., Comment: 11 pages, 5 figures

We use multi-frequency radiation hydrodynamics (rad-HD) to simulate radiative acceleration of a spherically symmetric stellar wind. We demonstrate the rad-HD capabilities of Athena++ for a series of test problems with multi-group radiation transfer. We then model the radiative transfer of a single spectral line through a spherically symmetric, isothermal, "CAK"-type line driven wind. We find that correctly accounting for the Doppler shift of the absorbed radiation, the force is well described by the analytic Sobolev line transfer in the supersonic parts of the solution where the flow is stationary and the effects of Abbott waves is negligible. Unlike in the analytic, steady-state solution re-radiation is important and leads to non-trivial radiation energy density and fluxes in the outer parts of the wind. We discuss a variety of applications to these multi-group methods that are currently computationally tractable., Comment: 10 pages, 6 figures

We present a high density disc reflection spectral analysis of a sample of 17 Seyfert 1 galaxies to study the inner disc densities at different black hole mass scales and accretion rates. All the available XMM-Newton observations in the archive are used. OM observations in the optical/UV band are used to estimate their accretion rates. We find that 65% of sources in our sample show a disc density significantly higher than 1e15 cm^{-3}, which was assumed in previous reflection-based spectral analyses. The best-fit disc densities show an anti-correlation with black hole mass and mass accretion rate. High density disc reflection model can successfully explain the soft excess emission and significantly reduce inferred iron abundances. We also compare our black hole spin and disc inclination angle measurements with previous analyses., Comment: Accepted by MNRAS. The first paper of this series can be found here: arXiv:1901.01739

Axions/axion-like particles (ALPs) are a well motivated extension of the Standard Model and are generic within String Theory. The X-ray transparency of the intracluster medium (ICM) in galaxy clusters is a powerful probe of light ALPs (with mass $<10^{-11}\,{\rm eV}$); as X-ray photons from an embedded or background source propagate through the magnetized ICM, they may undergo energy-dependent quantum mechanical conversion into ALPs (and vice versa), imprinting distortions on the X-ray spectrum. We present Chandra data for the active galactic nucleus NGC1275 at the center of the Perseus cluster. Employing a 490ks High-Energy Transmission Gratings (HETG) exposure, we obtain a high-quality 1-9keV spectrum free from photon pileup and ICM contamination. Apart from iron-band features, the spectrum is described by a power-law continuum, with any spectral distortions at the $<3\%$ level. We compute photon survival probabilities as a function of ALP mass $m_a$ and ALP-photon coupling constant $g_{a\gamma}$ for an ensemble of ICM magnetic field models, and then use the NGC1275 spectrum to constraint the $(m_a, g_{a\gamma})$-plane. Marginalizing over magnetic field realizations, the 99.7% credible region limits the ALP-photon coupling to $g_{a\gamma}<6-8\times 10^{-13}\, {\rm GeV}^{-1}$ (depending upon magnetic field model) for masses $m_a<1\times 10^{-12}\,{\rm eV}$. These are the most stringent limit to date on $g_{a\gamma}$ for these light ALPs, and have already reached the sensitivity limits of next-generation helioscopes and light-shining-through-wall experiments. We highlight the potential of these studies with the next-generation X-ray observatories Athena and Lynx, but note the critical importance of advances in relative calibration of these future X-ray spectrometers., Comment: Published in The Astrophysical Journal, 2020, vol 890, pp59. This posting corrects typos in equation (2) and the likelihood function just prior to equation (3). Typos were introduced into the draft at late stage and the analysis is correct

We investigate the interaction between active galactic nuclei (AGN) jets and the intracluster medium (ICM) of galaxy clusters. Specifically, we study the efficiency with which jets can drive sound waves into the ICM. Previous works focused on this issue model the jet-ICM interaction as a spherically symmetric explosion, finding that <~ 12.5% of the blast energy is converted into sound waves, even for instantaneous energy injection. We develop a method for measuring sound wave energy in hydrodynamic simulations and measure the efficiency of sound wave driving by supersonic jets in a model ICM. Our axisymmetric fiducial simulations convert >~ 25% of the jet energy into strong, long-wavelength sound waves which can propagate to large distances. Vigorous instabilities driven by the jet-ICM interaction generate small-scale sound waves which constructively interfere, forming powerful large-scale waves. By scanning a parameter space of opening angles, velocities, and densities, we study how our results depend on jet properties. High velocity, wide angle jets produce sound waves most efficiently, yet the acoustic efficiency never exceeds 1/3 of the jet energy -- an indication that equipartition may limit the nonlinear energy conversion process. Our work argues that sound waves may comprise a significant fraction of the energy budget in cluster AGN feedback and underscores the importance of properly treating compressive wave dissipation in the weakly collisional, magnetized ICM., Comment: 17 pages, 14 figures, accepted for publication in ApJ

In 2015, a radio transient named Cygnus A-2 was discovered in Cygnus A with the Very Large Array. Because of its radio brightness ($\nu F_{\nu} \approx 6 \times 10^{39}$ erg s$^{-1}$), this transient likely represents a secondary black hole in orbit around the AGN. Using {\it Chandra} ACIS observations from 2015 to 2017, we have looked for an X-ray counterpart to Cygnus A-2. The separation of 0.42 arcsec means that Cygnus A-2 can not be spatially resolved, but by comparing the data with simulated \texttt{marx} data, we put an upper limit to the 2-10 keV X-ray luminosity of Cygnus A-2 of $1 \times 10^{43}$ erg s$^{-1}$. Using the Fundamental Plane for accreting black holes, we find that our upper limit to the X-ray flux of Cygnus A-2 in 2015-2017 disfavours the interpretation of Cygnus A-2 as a steadily accreting black hole. We suggest instead that Cygnus A-2 is the radio afterglow of a tidal disruption event (TDE), and that a peak in the 2-10 keV luminosity of the nuclear region in 2013, when it was observed by {\it Swift} and {\it NuSTAR}, is X-ray emission from the TDE. A TDE could naturally explain the X-ray light curve of the nuclear region, as well as the appearance of a short-lived, fast, and ionized outflow previously detected in the 2013 {\it NuSTAR} spectrum. Both the radio and X-ray luminosities fall in between typical luminosities for 'thermal' and 'jetted' TDE types, suggesting that Cygnus A-2 would be unlike previously seen TDE's., Comment: 14 pages, 7 figures, accepted for publication in MNRAS

The spin of a black hole retains the memory of how the black hole grew, and can be a potent source of energy for powering relativistic jets. To understand the diagnostic power and astrophysical significance of black hole spin, however, we must first devise observational methods for measuring spin. Here, I describe the current state of black hole spin measurements, highlighting the progress made by X-ray astronomers, as well as the current excitement of gravitational wave and radio astronomy based techniques. Today's spin measurements are already constraining models for the growth of supermassive black holes and giving new insights into the dynamics of stellar core-collapse, as well as hinting at the physics of relativistic jet production. Future X-ray, radio, and gravitational wave observatories will transform black hole spin into a precision tool for astrophysics and test fundamental theories of gravity., Comment: Review article published in Nature Astronomy, 2019, vol 3, p41 (pre-acceptance, author resubmitted version posted here as per Nature policy)

Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline., Comment: A Probe-class mission study commissioned by NASA for the NAS Astro2020 Decadal Survey. Cost section redacted. 66 pages, 41 figures. v2: minor fixes

Intermediate-mass black holes (IMBHs) with masses between $100 - 10^5M_{\odot}$ are crucial to our understanding of black hole seed formation and are the prime targets for LISA, yet black holes in this mass range have eluded detection by traditional optical spectroscopic surveys aimed at finding active galactic nuclei (AGNs). In this paper, we have modeled for the first time the dependence of the optical narrow emission line strengths on the black hole mass of accreting AGN over the range of $100-10^8M_{\odot}$. We show that as the black hole mass decreases, the hardening of the spectral energy distribution from the accretion disk changes the ionization structure of the nebula. The enhanced high energy emission from IMBHs results in a more extended partially ionized zone compared with models for higher mass black holes. This effect produces a net decrease in the predicted [OIII]/H$\beta$ and [NII]/H$\alpha$ emission line ratios. Based on this model, we demonstrate that the standard optical narrow emission line diagnostics used to identify massive black holes fail when black hole mass falls below $\approx10^4M_{\odot}$ for highly accreting IMBHs and for radiatively inefficient IMBHs with active star formation. Our models call into question the ability of common optical spectroscopic diagnostics to confirm AGN candidates in dwarf galaxies, and indicate that the low-mass black hole occupation fraction inferred from such diagnostics will be severely biased., Comment: 10 pages, 6 figures, accepted to ApJL

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$., Comment: 19 pages, 18 figures, 6 tables. Accepted for publication in PASJ

X-ray reverberation in AGN, believed to be the result of the reprocessing of corona photons by the underlying accretion disk, has allowed us to probe the properties of the inner-most regions of the accretion flow and the central black hole. This process is modeled via raytracing in the Kerr metric, with the disk thickness almost ubiquitously assumed to be negligible (razor-thin) and the corona commonly approximated as a point source located along the polar axis (a lamppost). In this work, we use the new raytracing suite, Fenrir, to explore the effect that accretion disk geometry has on reverberation signatures, assuming a lamppost configuration but allowing for a finite disk scale height. We characterize the signatures of finite disk thickness in the reverberation transfer-function and calculate how they might manifest in observed lag-frequency spectra. We also show that a disk-hugging corona (approximated by off-axis point-like flares) exhibits characteristics that are qualitatively different from observation, thus providing further evidence for a flaring corona that is separated from the underlying disk material., Comment: 21 pages, 14 figures. Accepted for publication in ApJ

Active galactic nuclei (AGN) release a huge amount of energy into the intracluster medium (ICM) with the consequence of offsetting cooling and star formation (AGN feedback) in the centers of cool core clusters. The Phoenix cluster is among the most massive clusters of galaxies known in the Universe. It hosts a powerful starburst of several hundreds of Solar masses per year and a large amount of molecular gas in the center. In this work we use the high-resolution Reflection Grating Spectrometer (RGS) on board XMM-Newton to study the X-ray emitting cool gas in the Phoenix cluster and heating-cooling balance. We detect for the first time evidence of O VIII and Fe XXI-XXII emission lines, the latter demonstrating the presence of gas below 2 keV. We find a cooling rate of 350 (-200,+250) Msun/year below 2 keV (at the 90% confidence level), which is consistent with the star formation rate in this object. This cooling rate is high enough to produce the molecular gas found in the filaments via instabilities during the buoyant rising time. The line broadening indicates that the turbulence (~ 300 km/s or less) is below the level required to produce and propagate the heat throughout the cool core. This provides a natural explanation to the coexistence of large amounts of cool gas, star formation and a powerful AGN in the core. The AGN activity may be either at a young stage or in a different feedback mode, due to a high accretion rate., Comment: Accepted for publication in MNRAS. 11 Pages, 8 figures, 3 tables

A fraction of active galactic nuclei (AGN) exhibit dramatic variability, which is observed on timescales down to minutes in the X-ray band. We introduce the case study of 1H 1934-063 (z = 0.0102), a Narrow-line Seyfert 1 (NLS1) among the brightest and most variable AGN ever observed with XMM-Newton. This work includes spectral and temporal analyses of a concurrent XMM-Newton and NuSTAR 2015 observation lasting 130 kiloseconds, during which the X-ray source exhibited a steep (factor of 6) plummet and subsequent full recovery of flux level, accompanied by deviation from a single log-normal flux distribution. We rule out Compton-thin obscuration as the cause for this dramatic variability observed even at NuSTAR energies. In order to constrain coronal geometry, dynamics, and emission/absorption processes, we compare detailed spectral fitting with Fourier-based timing analysis. Similar to other well-studied, highly variable Seyfert 1s, this AGN is X-ray bright and displays strong reflection features. We find a narrower broad iron line component compared to most Seyfert 1s, and constrain black hole spin to be < 0.1, one of the lowest yet discovered for such systems. Combined spectral and timing results are consistent with a dramatic change in the continuum on timescales as short as a few kiloseconds dictating the nature of this variability. We also discover a Fe-K time lag, measuring a delay of 20 seconds between relativistically-blurred reflection off the inner accretion flow and the hard X-ray continuum emission., Comment: 13 pages, 11 figures, Accepted to ApJ

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter., Comment: 16 pages, 8 figures, 4 tables. Accepted for publication in PASJ

Galaxy cluster cores are pervaded by hot gas which radiates at far too high a rate to maintain any semblance of a steady state; this is referred to as the cooling flow problem. Of the many heating mechanisms that have been proposed to balance radiative cooling, one of the most attractive is dissipation of acoustic waves generated by Active Galactic Nuclei (AGN). Fabian (2005) showed that if the waves are nearly adiabatic, wave damping due to heat conduction and viscosity must be well below standard Coulomb rates in order to allow the waves to propagate throughout the core. Because of the importance of this result, we have revisited wave dissipation under galaxy cluster conditions in a way that accounts for the self limiting nature of dissipation by electron thermal conduction, allows the electron and ion temperature perturbations in the waves to evolve separately, and estimates kinetic effects by comparing to a semi-collisionless theory. While these effects considerably enlarge the toolkit for analyzing observations of wavelike structures and developing a quantitative theory for wave heating, the drastic reduction of transport coefficients proposed in Fabian (2005) remains the most viable path to acoustic wave heating of galaxy cluster cores., Comment: Accepted for publication in Astrophys. J

Cluster cool cores possess networks of line-emitting filaments. These filaments are thought to originate via uplift of cold gas from cluster centers by buoyant active galactic nuclei (AGN) bubbles, or via local thermal instability in the hot intracluster medium (ICM). Therefore, the filaments are either the signatures of AGN feedback or feeding of supermassive black holes. Despite being characterized by very short cooling times, the filaments are significant H$\alpha$ emitters, which suggests that some process continuously powers these structures. Many cool cores host diffuse radio mini halos and AGN injecting radio plasma, suggesting that cosmic rays (CRs) and magnetic fields are present in the ICM. We argue that the excitation of Alfv\'en waves by CR streaming, and the replenishment of CR energy via accretion onto the filaments of high plasma-$\beta$ ICM characterized by low CR pressure support, can provide the adequate amount of heating to power and sustain the emission from these filaments. This mechanism does not require the CRs to penetrate the filaments even if the filaments are magnetically isolated from the ambient ICM and it may operate irrespectively of whether the filaments are dredged up from the center or form in situ in the ICM. This picture is qualitatively consistent with non-thermal line ratios seen in the cold filaments. Future X-ray observations of the iron line complex with XARM, Lynx, or Athena could help to test this model by providing constraints on the amount of CRs in the hot plasma that is cooling and accreting onto the filaments., Comment: ApJ, accepted

We investigate the role of AGN feedback in turbulent heating of galaxy clusters. Specifically, we analyze the production of turbulence by g-modes generated by the supersonic expansion and buoyant rise of AGN-driven bubbles. Previous work which neglects magnetic fields has shown that this process is inefficient, with less than 1% of the injected energy ending up in turbulence. This inefficiency is primarily due to the fact that the bubbles are shredded apart by hydrodynamic instabilities before they can excite sufficiently strong g-modes. Using a plane-parallel model of the ICM and 3D ideal MHD simulations, we examine the role of a large-scale magnetic field which is able to drape around these rising bubbles, preserving them from hydrodynamic instabilities. We find that, while magnetic draping appears better able to preserve AGN-driven bubbles, the driving of g-modes and the resulting production of turbulence is still inefficient. The magnetic tension force prevents g-modes from transitioning into the nonlinear regime, suppressing turbulence in our model ICM. Our work highlights the ways in which ideal MHD is an insufficient description for the cluster feedback process, and we discuss future work such as the inclusion of anisotropic viscosity as a means of simulating high $\beta$ plasma kinetic effects. These results suggest the hypothesis that other mechanisms of heating the ICM plasma such as sound waves or cosmic rays may be responsible for observed feedback in galaxy clusters., Comment: 14 pages, 10 figures, Accepted for publication in The Astrophysical Journal

We study a truncated accretion disk using a well-resolved, semi-global magnetohydrodynamic simulation that is evolved for many dynamical times (6096 inner disk orbits). The spectral properties of hard state black hole binary systems and low-luminosity active galactic nuclei are regularly attributed to truncated accretion disks, but a detailed understanding of the flow dynamics is lacking. In these systems the truncation is expected to arise through thermal instability driven by sharp changes in the radiative efficiency. We emulate this behavior using a simple bistable cooling function with efficient and inefficient branches. The accretion flow takes on an arrangement where a "transition zone" exists in between hot gas in the inner most regions and a cold, Shakura $\&$ Sunyaev thin disk at larger radii. The thin disk is embedded in an atmosphere of hot gas that is fed by a gentle outflow originating from the transition zone. Despite the presence of hot gas in the inner disk, accretion is efficient. Our analysis focuses on the details of the angular momentum transport, energetics, and magnetic field properties. We find that the magnetic dynamo is suppressed in the hot, truncated inner region of the disk which lowers the effective $\alpha$-parameter by $65\%$., Comment: 18 pgs, 16 figures, accepted for publication in ApJ

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band., Comment: 29 pages, 19 figures, 9 tables. Accepted for publication in PASJ

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission., Comment: 46 pages, 25 figures, 11 tables. Accepted for publication in PASJ

The relativistically-broadened reflection spectrum, observed in both AGN and X-ray binaries, has proven to be a powerful probe of the properties of black holes and the environments in which they reside. Being emitted from the inner-most regions of the accretion disk, this X-ray spectral component carries with it information not only about the plasma that resides in these extreme conditions, but also the black hole spin, a marker of the formation and accretion history of these objects. The models currently used to interpret the reflection spectrum are often simplistic, however, approximating the disk as an infinitely thin, optically thick plane of material orbiting in circular Keplerian orbits around the central object. Using a new relativistic ray tracing suite (Fenrir) that allows for more complex disk approximations, we examine the effects that disk thickness may have on the reflection spectrum. Assuming a lamp post corona, we find that finite disk thickness can have a variety of effects on the reflection spectrum, including a truncation of the blue wing (from self-shadowing of the accretion disk) and an enhancement of the red wing (from the irradiation of the central 'eye wall' of the inner disk). We make a first estimate on the systematic errors on black hole spin and height that may result from neglecting these effects., Comment: 20 pages, 10 figures, Accepted for publication in ApJ