Source author record

Christopher S. Reynolds

Christopher S. Reynolds appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph.CO astro-ph.GA astro-ph gr-qc hep-ph

Catalog footprint

What is connected

56works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A Disc Reflection Model for Ultra-Soft Narrow-Line Seyfert 1 Galaxies

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 $Γ\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 $λ_{\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.

preprint2022arXiv

Acoustic waves and g-mode turbulence as energy carriers in a viscous intracluster medium

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.

preprint2022arXiv

Black Hole Spin Measurements Based on a Thin Disc Model with Finite Thickness I. An example study of MCG-06-30-15

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.

preprint2022arXiv

Evidence for a moderate spin from X-ray reflection of the high-mass supermassive black hole in the cluster-hosted quasar H1821+643

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.

preprint2022arXiv

How do Magnetic Field Models Affect Astrophysical Limits on Light Axion-like Particles? An X-ray Case Study with NGC 1275

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 γ}$, 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 γ}$ (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 $β_{\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 $β_{\rm pl}$.

preprint2022arXiv

Probing the Extent of Fe K$α$ Emission in Nearby Active Galactic Nuclei using Multi-Order Analysis of Chandra High Energy Transmission Grating Data

We present a study of the narrow Fe K$α$ 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$α$ 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$α$ emission. We utilize these differences in narrow Fe K$α$ 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.

preprint2022arXiv

XMM-Newton Observations of the Narrow-Line Seyfert 1 Galaxy IRAS 13224$-$3809 : X-ray Spectral Analysis II

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.

preprint2021arXiv

New constraints on light Axion-Like Particles using Chandra Transmission Grating Spectroscopy of the powerful cluster-hosted quasar H1821+643

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 $β= 100$, we exclude all couplings $g_\mathrm{aγ} > 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.

preprint2020arXiv

Lense-Thirring Precession of Misaligned Discs I

We study Lense-Thirring precession of inviscid and viscous misaligned $α-$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.

preprint2020arXiv

Probing the Milky Way's Dark Matter Halo for the 3.5 keV Line

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$θ$) 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.

preprint2016arXiv

Revealing the broad iron Kalpha line in Cygnus X-1 through simultaneous XMM-Newton, RXTE, and INTEGRAL observations

We report on the analysis of the broad Fe Kalpha line feature of Cygnus X-1 in the spectra of four simultaneous hard intermediate state observations made with the X-ray Multiple Mirror mission (XMM-Newton), the Rossi X-ray Timing Explorer (RXTE), and the International Gamma-Ray Astrophysics Laboratory (INTEGRAL). The high quality of the XMM-Newton data taken in the Modified Timing Mode of the EPIC-pn camera provides a great opportunity to investigate the broadened Fe Kalpha reflection line at 6.4keV with a very high signal to noise ratio. The 4-500keV energy range is used to constrain the underlying continuum and the reflection at higher energies. We first investigate the data by applying a phenomenological model that consists of the sum of an exponentially cutoff power law and relativistically smeared reflection. Additionally, we apply a more physical approach and model the irradiation of the accretion disk directly from the lamp post geometry. All four observations show consistent values for the black hole parameters with a spin of $a\sim 0.9$, in agreement with recent measurements from reflection and disk continuum fitting. The inclination is found to be $i\sim30^\circ$, consistent with the orbital inclination and different from inclination measurements made during the soft state, which show a higher inclination. We speculate that the difference between the inclination measurements is due to changes in the inner region of the accretion disk.

preprint2016arXiv

The view of AGN-host alignment via reflection spectroscopy

The fuelling of active galactic nuclei (AGN) - via material propagated through the galactic disc or via minor mergers - is expected to leave an imprint on the alignment of the sub-pc disc relative to the host galaxy's stellar disc. Determining the inclination of the inner disc usually relies on the launching angle of the jet; here instead we use the inclination derived from reflection fits to a sample of AGN. We determine the distorting effect of unmodeled Fe XXV/XXVI features and, via extensive simulations, determine the difference in disc inclination resulting from the use of relxill compared to reflionx. We compare inner disc inclinations to those for the host galaxy stellar disc derived from the Hubble formula and, via Monte-Carlo simulations, find a strong lack of a correlation (at >> 5-sigma) implying either widespread feeding via mergers if we assume the sample to be homogeneous, or that radiative disc warps are distorting our view of the emission. However, we find that by removing a small (~1/5) subset of AGN, the remaining sample is consistent with random sampling of an underlying 1:1 correlation (at the 3-sigma level). A heterogenous sample would likely imply that our view is not dominated by radiative disc warps but instead by different feeding mechanisms with the majority consistent with coplanar accretion (although this may be the result of selection bias), whilst a smaller but not insignificant fraction may have been fuelled by minor mergers in the recent history of the host galaxy.

preprint2015arXiv

Accretion disk dynamo as the trigger for X-ray binary state transitions

Magnetohydrodynamic accretion disk simulations suggest that much of the energy liberated by the magnetorotational instability (MRI) can be channeled into large-scale toroidal magnetic fields through dynamo action. Under certain conditions, this field can dominate over gas and radiation pressure in providing vertical support against gravity, even close to the midplane. Using a simple model for the creation of this field, its buoyant rise, and its coupling to the gas, we show how disks could be driven into this magnetically dominated state and deduce the resulting vertical pressure and density profiles. Applying an established criterion for MRI to operate in the presence of a toroidal field, we show that magnetically supported disks can have two distinct MRI-active regions, separated by a "dead zone" where local MRI is suppressed, but where magnetic energy continues to flow upward from the dynamo region below. We suggest that the relative strengths of the MRI zones, and the local poloidal flux, determine the spectral states of X-ray binaries. Specifically, "intermediate" and "hard" accretion states occur when MRI is triggered in the hot, upper zone of the corona, while disks in "soft" states do not develop the upper MRI zone. We discuss the conditions under which various transitions should take place and speculate on the relationship of dynamo activity to the various types of quasi-periodic oscillations that sometimes appear in the hard spectral components. The model also explains why luminous accretion disks in the "soft" state show no signs of the thermal/viscous instability predicted by standard alpha models.

preprint2015arXiv

Inefficient Driving of Bulk Turbulence by Active Galactic Nuclei in a Hydrodynamic Model of the Intracluster Medium

Central jetted active galactic nuclei (AGN) appear to heat the core regions of the intracluster medium (ICM) in cooling-core galaxy clusters and groups, thereby preventing a cooling catastrophe. However, the physical mechanism(s) by which the directed flow of kinetic energy is thermalized throughout the ICM core remains unclear. We examine one widely discussed mechanism whereby the AGN induces subsonic turbulence in the ambient medium, the dissipation of which provides the ICM heat source. Through controlled inviscid 3-d hydrodynamic simulations, we verify that explosive AGN-like events can launch gravity waves (g-modes) into the ambient ICM which in turn decay to volume-filling turbulence. In our model, however, this process is found to be inefficient, with less than 1% of the energy injected by the AGN activity actually ending up in the turbulence of the ambient ICM. This efficiency is an order of magnitude or more too small to explain the observations of AGN-feedback in galaxy clusters and groups with short central cooling times. Atmospheres in which the g-modes are strongly trapped/confined have an even lower efficiency since, in these models, excitation of turbulence relies on the g-modes' ability to escape from the center of the cluster into the bulk ICM. Our results suggest that, if AGN-induced turbulence is indeed the mechanism by which the AGN heats the ICM core, its driving may rely on physics beyond that captured in our ideal hydrodynamic model.

preprint2015arXiv

NuSTAR observations of the powerful radio-galaxy Cygnus A

We present NuSTAR observations of the powerful radio galaxy Cygnus A, focusing on the central absorbed active galactic nucleus (AGN). Cygnus A is embedded in a cool-core galaxy cluster, and hence we also examine archival XMM-Newton data to facilitate the decomposition of the spectrum into the AGN and intracluster medium (ICM) components. NuSTAR gives a source-dominated spectrum of the AGN out to >70keV. In gross terms, the NuSTAR spectrum of the AGN has the form of a power law (Gamma~1.6-1.7) absorbed by a neutral column density of N_H~1.6x10^23 cm^-2. However, we also detect curvature in the hard (>10keV) spectrum resulting from reflection by Compton-thick matter out of our line-of-sight to the X-ray source. Compton reflection, possibly from the outer accretion disk or obscuring torus, is required even permitting a high-energy cutoff in the continuum source; the limit on the cutoff energy is E_cut>111keV (90% confidence). Interestingly, the absorbed power-law plus reflection model leaves residuals suggesting the absorption/emission from a fast (15,000-26,000km/s), high column-density (N_W>3x10^23 cm^-2), highly ionized (xi~2,500 erg cm/s) wind. A second, even faster ionized wind component is also suggested by these data. We show that the ionized wind likely carries a significant mass and momentum flux, and may carry sufficient kinetic energy to exercise feedback on the host galaxy. If confirmed, the simultaneous presence of a strong wind and powerful jets in Cygnus A demonstrates that feedback from radio-jets and sub-relativistic winds are not mutually exclusive phases of AGN activity but can occur simultaneously.

preprint2015arXiv

Unwrapping the X-ray Spectra of Active Galactic Nuclei

Active galactic nuclei (AGN) are complex phenomena. At the heart of an AGN is a relativistic accretion disk around a spinning supermassive black hole (SMBH) with an X-ray emitting corona and, sometimes, a relativistic jet. On larger scales, the outer accretion disk and molecular torus act as the reservoirs of gas for the continuing AGN activity. And on all scales from the black hole outwards, powerful winds are seen that probably affect the evolution of the host galaxy as well as regulate the feeding of the AGN itself. In this review article, we discuss how X-ray spectroscopy can be used to study each of these components. We highlight how recent measurements of the high-energy cutoff in the X-ray continuum by NuSTAR are pushing us to conclude that X-ray coronae are radiatively-compact and have electron temperatures regulated by electron-positron pair production. We show that the predominance of rapidly-rotating objects in current surveys of SMBH spin is entirely unsurprising once one accounts for the observational selection bias resulting from the spin-dependence of the radiative efficiency. We review recent progress in our understanding of fast (v~0.1-0.3c), highly-ionized (mainly visible in FeXXV and FeXXVI lines), high-column density winds that may dominate quasar-mode galactic feedback. Finally, we end with a brief look forward to the promise of Astro-H and future X-ray spectropolarimeters.

preprint2014arXiv

An XMM-Newton View of the Radio Galaxy 3C 411

We present the first high signal-to-noise XMM-Newton observations of the broad-line radio galaxy 3C 411. After fitting various spectral models, an absorbed double power-law continuum and a blurred relativistic disk reflection model (kdblur) are found to be equally plausible descriptions of the data. While the softer power-law component ($Γ$=2.11) of the double power-law model is entirely consistent with that found in Seyfert galaxies (and hence likely originates from a disk corona), the additional power law component is very hard ($Γ$=1.05); amongst the AGN zoo, only flat-spectrum radio quasars have such hard spectra. Together with the very flat radio-spectrum displayed by this source, we suggest that it should instead be classified as a FSRQ. This leads to potential discrepancies regarding the jet inclination angle, with the radio morphology suggesting a large jet inclination but the FSRQ classification suggesting small inclinations. The kdblur model predicts an inner disk radius of at most 20 r$_g$ and relativistic reflection.

preprint2014arXiv

Modeling Hot Gas Flow in the Low-Luminosity Active Galactic Nucleus of NGC3115

Based on the dynamical black hole (BH) mass estimates, NGC3115 hosts the closest billion solar mass BH. Deep studies of the center revealed a very underluminous active galactic nucleus (AGN) immersed in an old massive nuclear star cluster. Recent $1$~Ms \textit{Chandra} X-ray visionary project observations of the NGC3115 nucleus resolved hot tenuous gas, which fuels the AGN. In this paper we connect the processes in the nuclear star cluster with the feeding of the supermassive BH. We model the hot gas flow sustained by the injection of matter and energy from the stars and supernova explosions. We incorporate electron heat conduction as the small-scale feedback mechanism, the gravitational pull of the stellar mass, cooling, and Coulomb collisions. Fitting simulated X-ray emission to the spatially and spectrally resolved observed data, we find the best-fitting solutions with $χ^2/dof=1.00$ for $dof=236$ both with and without conduction. The radial modeling favors a low BH mass $<1.3\times10^{9}M_\odot$. The best-fitting supernova rate and the best-fitting mass injection rate are consistent with their expected values. The stagnation point is at $r_{\rm st}\lesssim1$arcsec, so that most of gas, including the gas at a Bondi radius $r_B=2-4$arcsec, outflows from the region. We put an upper limit on the accretion rate at $2\times10^{-3}M_\odot{\rm yr}^{-1}$. We find a shallow density profile $n\propto r^{-β}$ with $β\approx1$ over a large dynamic range. This density profile is determined in the feeding region $0.5-10$arcsec as an interplay of four processes and effects: (1) the radius-dependent mass injection, (2) the effect of the galactic gravitational potential, (3) the accretion flow onset at $r\lesssim1$arcsec, and (4) the outflow at $r\gtrsim1$arcsec. The gas temperature is close to the virial temperature $T_v$ at any radius.

preprint2014arXiv

Non-Linear Dynamics of Accretion Disks with Stochastic Viscosity

We present a non-linear numerical model for a geometrically thin accretion disk with the addition of stochastic non-linear fluctuations in the viscous parameter. These numerical realizations attempt to study the stochastic effects on the disk angular momentum transport. We show that this simple model is capable of reproducing several observed phenomenologies of accretion driven systems. The most notable of these is the observed linear rms-flux relationship in the disk luminosity. This feature is not formally captured by the linearized disk equations used in previous work. A Fourier analysis of the dissipation and mass accretion rates across disk radii show coherence for frequencies below the local viscous frequency. This is consistent with the coherence behavior observed in astrophysical sources such as Cygnus X-1.

preprint2014arXiv

The hard X-ray perspective on the soft X-ray excess

The X-ray spectra of many active galactic nuclei (AGN) exhibit a `soft excess' below 1keV, whose physical origin remains unclear. Diverse models have been suggested to account for it, including ionised reflection of X-rays from the inner part of the accretion disc, ionised winds/absorbers, and Comptonisation. The ionised reflection model suggests a natural link between the prominence of the soft excess and the Compton reflection hump strength above 10keV, but it has not been clear what hard X-ray signatures, if any, are expected from the other soft X-ray candidate models. Additionally, it has not been possible up until recently to obtain high-quality simultaneous measurements of both soft and hard X-ray emission necessary to distinguish these models, but upcoming joint XMM-NuSTAR programmes provide precisely this opportunity. In this paper, we present an extensive analysis of simulations of XMM+NuSTAR observations, using two candidate soft excess models as inputs, to determine whether such campaigns can disambiguate between them by using hard and soft X-ray observations in tandem. The simulated spectra are fit with the simplest "observer's model" of a black body and neutral reflection to characterise the strength of the soft and hard excesses. A plot of the strength of the hard excess against the soft excess strength provides a diagnostic plot which allows the soft excess production mechanism to be determined in individual sources and samples using current state-of-the-art and next generation hard X-ray enabled observatories. This approach can be straightforwardly extended to other candidate models for the soft excess.

preprint2014arXiv

The role of Compton heating in radiation-regulated accretion on to black holes

We investigate the role of Compton heating in radiation-regulated accretion on to black holes from a neutral dense medium using 1D radiation-hydrodynamic simulations. We focus on the relative effects of Compton-heating and photo-heating as a function of the spectral slope α, assuming a power-law spectrum in the energy range of 13.6 eV--100 keV. While Compton heating is dominant only close to the black hole, it can reduce the accretion rate to 0.1 % ($l \propto \dot{m}^2$ model)--0.01 % ($l \propto \dot{m}$ model) of the Bondi accretion rate when the BH radiation is hard (α ~ 1), where $l$ and $\dot{m}$ are the luminosity and accretion rate normalised by Eddington rates, respectively. The oscillatory behaviour otherwise typically seen in simulations with α > 1, become suppressed when α ~ 1 only for the $l \propto \dot{m}$ model. The relative importance of the Compton heating over photo-heating decreases and the oscillatory behaviour becomes stronger as the spectrum softens. When the spectrum is soft (α > 1.5), photo-heating prevails regardless of models making the effect of Compton heating negligible. On the scale of the ionization front, where the gas supply into the Strömgren sphere from large scale is regulated, photo-heating dominates. Our simulations show consistent results with the advection-dominated accretion flow ($l \propto \dot{m}^2$) where the accretion is inefficient and the spectrum is hard (α ~ 1).

preprint2014arXiv

The X-ray spectrum of the cooling-flow quasar H1821+643 : A massive black hole feeding off the intracluster medium

We present a deep Suzaku observation of H1821+643, an extremely rare example of a powerful quasar hosted by the central massive galaxy of a rich cooling-core cluster of galaxies. Informed by previous Chandra studies of the cluster, we achieve a spectral separation of emission from the active galactic nucleus (AGN) and the intracluster medium (ICM). With a high degree of confidence, we identify the signatures of X-ray reflection/reprocessing by cold and slowly moving material in the AGN's immediate environment. The iron abundance of this matter is found to be significantly sub-solar (Z~0.4Zsun), an unusual finding for powerful AGN but in line with the idea that this quasar is feeding from the ICM via a Compton-induced cooling flow. We also find a subtle soft excess that can be described phenomenologically (with an additional black body component) or as ionized X-ray reflection from the inner regions of a high inclination (i=57 degrees) accretion disk around a spinning (a>0.4) black hole. We describe how the ionization state of the accretion disk can be used to constrain the Eddington fraction of the source. Applying these arguments to our spectrum implies an Eddington fraction of 0.25-0.5, with an associated black hole mass of 3-6x10^9Msun.

preprint2013arXiv

An X-Ray View of the Jet-Cycle in the Radio Loud AGN 3C120

We present a study of the central engine in the broad-line radio galaxy 3C120 using a multi-epoch analysis of a deep XMM-Newton observation and two deep Suzaku pointings (in 2012). In order to place our spectral data into the context of the disk-disruption/jet-ejection cycles displayed by this object, we monitor the source in the UV/X-ray bands, and in the radio band. We find three statistically acceptable spectral models, a disk-reflection model, a jet-model and a jet+disk model. Despite being good descriptions of the data, the disk-reflection model violates the radio constraints on the inclination, and the jet-model has a fine-tuning problem, requiring a jet contribution exceeding that expected. Thus, we argue for a composite jet+disk model. Within the context of this model, we verify the basic predictions of the jet-cycle paradigm, finding a truncated/refilling disk during the Suzaku observations and a complete disk extending down to the innermost stable circular orbit (ISCO) during the XMM-Newton observation. The idea of a refilling disk is further supported by the detection of the ejection of a new jet knot approximately one month after the Suzaku pointings. We also discover a step-like event in one of the Suzaku pointings in which the soft band lags the hard band. We suggest that we are witnessing the propagation of a disturbance from the disk into the jet on a timescale set by the magnetic field.

preprint2013arXiv

GRB060218 as a Tidal Disruption of a White Dwarf by an Intermediate Mass Black Hole

A highly unusual pair of a gamma-ray burst (GRB) GRB060218 and an associated supernova SN2006aj has puzzled theorists for years. A supernova shock breakout and a jet from a newborn stellar mass compact object were put forward to explain its multiwavelength signature. We propose that the source is naturally explained by another channel, a tidal disruption of a white dwarf (WD) by an intermediate mass black hole (IMBH). The tidal disruption is accompanied by a tidal pinching, which leads to the ignition of a WD and a supernova. Some debris falls back onto the IMBH, forms a disk, which quickly amplifies the magnetic field, and launches a jet. We successfully fit soft X-ray spectrum with the Comptonized blackbody emission from a jet photosphere. The optical/UV emission is consistent with self-absorbed synchrotron from the expanding jet front. The accretion rate temporal dependence Mdot(t) in a tidal disruption provides a good fit to soft X-ray lightcurve. The IMBH mass is found to be about 10^4Msun in three independent estimates: (1) fitting tidal disruption Mdot(t) to soft X-ray lightcurve; (2) computing the jet base radius in a jet photospheric emission model; (3) inferring the central BH mass based on a host dwarf galaxy stellar mass. The supernova position is consistent with the center of the host galaxy, while low supernova ejecta mass is consistent with a WD mass. High expected rate of tidal disruptions in dwarf galaxies is consistent with one source observed by Swift satellite over several years at GRB060218 distance of 150Mpc. The encounters with the WDs provide a lot of fuel for IMBH growth.

preprint2013arXiv

Irradiation of an Accretion Disc by a Jet: General Properties and Implications for Spin Measurements of Black Holes

X-ray irradiation of the accretion disc leads to strong reflection features, which are then broadened and distorted by relativistic effects. We present a detailed, general relativistic approach to model this irradiation for different geometries of the primary X-ray source. These geometries include the standard point source on the rotational axis as well as more jet-like sources, which are radially elongated and accelerating. Incorporating this code in the relline model for relativistic line emission, the line shape for any configuration can be predicted. We study how different irradiation geometries affect the determination of the spin of the black hole. Broad emission lines are produced only for compact irradiating sources situated close to the black hole. This is the only case where the black hole spin can be unambiguously determined. In all other cases the line shape is narrower, which could either be explained by a low spin or an elongated source. We conclude that for all those cases and independent of the quality of the data, no unique solution for the spin exists and therefore only a lower limit of the spin value can be given.

preprint2013arXiv

Isotropic Heating of Galaxy Cluster Cores via Rapidly Reorienting AGN Jets

AGN jets carry more than sufficient energy to stave off catastrophic cooling of the intracluster medium (ICM) in the cores of cool-core clusters. However, in order to prevent catastrophic cooling, the ICM must be heated in a near-isotropic fashion and narrow bipolar jets with $P_{\rm jet}=10^{44-45}$ ergs/s, typical of radio AGNs at cluster centres, are inefficient at heating the gas in the transverse direction to the jets. We argue that due to existent conditions in cluster cores, the SMBHs will, in addition to accreting gas via radiatively inefficient flows, experience short stochastic episodes of enhanced accretion via thin discs. In general, the orientation of these accretion discs will be misaligned with the spin axis of the black holes and the ensuing torques will cause the black hole's spin axis (and therefore, the jet axis) to slew and rapidly change direction. This model not only explains recent observations showing successive generations of jet-lobes-bubbles in individual cool-core clusters that are offset from each other in the angular direction with respect to the cluster center, but also shows that AGN jets {\it can} heat the cluster core nearly isotropically on the gas cooling timescale. Our model {\it does} require that the SMBHs at the centers of cool-core clusters be spinning relatively slowly. Torques from individual misaligned discs are ineffective at tilting rapidly spinning black holes by more than a few degrees. Additionally, since SMBHs that host thin accretion discs will manifest as quasars, we predict that roughly 1--2 rich clusters within $z<0.5$ should have quasars at their centers.

preprint2013arXiv

Measuring Black Hole Spin using X-ray Reflection Spectroscopy

I review the current status of X-ray reflection (a.k.a. broad iron line) based black hole spin measurements. This is a powerful technique that allows us to measure robust black hole spins across the mass range, from the stellar-mass black holes in X-ray binaries to the supermassive black holes in active galactic nuclei. After describing the basic assumptions of this approach, I lay out the detailed methodology focusing on "best practices" that have been found necessary to obtain robust results. Reflecting my own biases, this review is slanted towards a discussion of supermassive black hole (SMBH) spin in active galactic nuclei (AGN). Pulling together all of the available XMM-Newton and Suzaku results from the literature that satisfy objective quality control criteria, it is clear that a large fraction of SMBHs are rapidly-spinning, although there are tentative hints of a more slowly spinning population at high (M>5*10^7Msun) and low (M<2*10^6Msun) mass. I also engage in a brief review of the spins of stellar-mass black holes in X-ray binaries. In general, reflection-based and continuum-fitting based spin measures are in agreement, although there remain two objects (GROJ1655-40 and 4U1543-475) for which that is not true. I end this review by discussing the exciting frontier of relativistic reverberation, particularly the discovery of broad iron line reverberation in XMM-Newton data for the Seyfert galaxies NGC4151, NGC7314 and MCG-5-23-16. As well as confirming the basic paradigm of relativistic disk reflection, this detection of reverberation demonstrates that future large-area X-ray observatories such as LOFT will make tremendous progress in studies of strong gravity using relativistic reverberation in AGN.

preprint2013arXiv

Rayleigh-Taylor instability of ionization front around black holes

We examine the role of ionizing radiation emitted from black holes (BHs) in suppressing the growth of the Rayleigh-Taylor instability (RTI) across the ionization front (I-front) that forms when the gas fueling the BH is neutral. We use radiation-hydrodynamic simulations to show that the RTI is suppressed for non-accelerating fronts on all scales resolved in our simulations. A necessary condition for the stability of the I-front is that the radius of the Strömgren sphere is larger than the Bondi radius. When this condition is violated the I-front collapses producing an accretion luminosity burst. Transient growth of the RTI occurs only during the accretion burst when the effective acceleration in the frame of reference of the I-front increases significantly due to the rapid expansion of the Strömgren sphere.

preprint2013arXiv

Tackling the soft X-ray excess in AGN with variability studies

The origin of the soft X-ray excess in AGN has been a mystery ever since its discovery. We present how the time variability of this spectral component can point towards its origin. Using the powerful technique of multi-epoch fitting, we study how the soft excess in a given object depends on other parameters of the continuum and the accretion disk possibly hinting at its nature. As an example, we present results from this technique applied to the Seyfert galaxy Mrk 841. We study all (3) XMM and some of the Suzaku pointings available and find that the source displays an impressive variability in the soft X-ray band on the timescale of years. We study several common soft excess models and their ability to physically consistently explain this spectral variability. Mrk 841 is found to show a distinct variability pattern that can be best explained by the soft excess originating mostly from a thermal Comptonization component. The variability timescale can be constrained to be on the order of a few days.

preprint2013arXiv

The Ionized Absorber and Nuclear Environment of IRAS 13349+2438: Multi-wavelength insights from coordinated Chandra HETGS, HST STIS, HET, and Spitzer IRS

We present results from a coordinated IR-to-X-ray spectral campaign of the QSO IRAS 13349+2438. Optical spectra reveal extreme Eigenvector-1 characteristics, but the H-beta line width argues against a NLS1 classification; we refine z=0.10853 based on [O III]. We estimate a BH mass=10^9 Msun using 2 independent methods (H-beta line width & SED fits). Blue-shifted absorption (-950km/s & -75km/s) is seen for the 1st time in STIS UV spectra from Ly-alpha, NV, & CIV. The higher velocity UV lines are coincident with the lower-ionisation (xi~1.6) X-ray warm absorber lines. A dusty multiple ionization absorber blueshifted by 700-900km/s is required to fit the X-ray data. Theoretical models comparing different ionising SEDs reveal that a UV-inclusive (i.e., the accretion disc) ionising continuum strongly impacts conclusions for the thermodynamic stability of the warm absorber. Specific to IRAS13349, an Xray-UV ionising SED favors a continuous distribution of ionisation states in a smooth flow (this paper), versus discrete clouds in pressure equilibrium (work by others where UV is omitted). Direct dust detections are seen in both the IR: PAH emission at (7.7 & 11.3)micron which may also be blended with forsterite, and (10 & 18)micron silicate emission, and X-rays: iron dust with a dust-to-gas ratio > 90%. We develop a geometrical model whereby the QSO nuclear region is viewed through the upper atmosphere of an obscuring torus. This sight line is obscured by dust that blocks a direct view of the UV/optical emission region but is largely transparent in X-rays since the gas is ionised. In our model, 20% of the intrinsic UV/optical continuum is scattered into our sight line by the far wall of an obscuring torus. An additional 2.4% of the direct light, which likely dominates the UV emission, is Thomson-scattered into our line-of-sight by another off-plane component of highly ionized gas.

preprint2013arXiv

The Spin of Supermassive Black Holes

Black hole spin is a quantity of great interest to both physicists and astrophysicists. We review the current status of spin measurements in supermassive black holes (SMBH). To date, every robust SMBH spin measurement uses X-ray reflection spectroscopy, so we focus almost exclusively on this technique as applied to moderately-luminous active galactic nuclei (AGN). After describing the foundations and uncertainties of the method, we summarize the current status of the field. At the time of writing, observations by XMM-Newton, Suzaku and NuSTAR have given robust spin constraints on 22 SMBHs. We find a significant number of rapidly-rotating SMBHs (with dimensionless spin parameters a>0.9) although, especially at the higher masses (M>4*10^7Msun), there are also some SMBHs with intermediate spin parameters. This may be giving us our first hint at a mass-dependent spin distribution which would, in turn, provide interesting constraints on models for SMBH growth. We also discuss the recent discovery of relativistic X-ray reverberation which we can use to "echo map" the innermost regions of the accretion disk. The ultimate development of these reverberation techniques, when applied to data from future high-throughput X-ray observatories such as LOFT, ATHENA+, and AXSIO, will permit the measurement of black hole spin by a characterization of strong-field Shapiro delays. We conclude with a brief discussion of other electromagnetic methods that have been attempted or are being developed to constrain SMBH spin.

preprint2013arXiv

What is on Tap? The Role of Spin in Compact Objects and Relativistic Jets

We examine the role of spin in launching jets from compact objects across the mass scale. Our work includes a total of 37 Seyferts, 11 stellar-mass black holes, and 13 neutron stars. We find that when the Seyfert reflection lines are modeled with Gaussian line features (a crude proxy for inner disk radius and therefore spin), only a slight inverse correlation is found between the Doppler-corrected radio luminosity at 5 GHz (a proxy for jet power) and line width. When the Seyfert reflection features are fit with relativistically-blurred disk reflection models that measure spin, there is a tentative positive correlation between the Doppler-corrected radio luminosity and the spin measurement. Further, when we include stellar-mass black holes in the sample, to examine the effects across the mass scale, we find a slightly stronger correlation with radio luminosity per unit mass and spin, at a marginal significance (2.3 sigma confidence level). Finally, when we include neutron stars, in order to probe lower spin values, we find a positive correlation (3.3 sigma confidence level) between radio luminosity per unit mass and spin. Although tentative, these results suggest that spin may have a role in determining the jet luminosity. In addition, we find a slightly more significant correlation (4.4 sigma confidence level) between radio luminosity per Bolometric luminosity and spin, using our entire sample of black holes and neutrons stars. Again, although tentative, these relations point to the possibility that the mass accretion rate, i.e. Bolometric luminosity, is also important in determining the jet luminosity, in addition to spin. Our analysis suggests that mass accretion rate and disk or coronal magnetic field strength may be the "throttle" in these compact systems, to which the Eddington limit and spin may set the maximum jet luminosity that can be achieved.

preprint2012arXiv

Buoyancy Instabilities in a Weakly Collisional Intracluster Medium

The intracluster medium of galaxy clusters is a weakly collisional, high-beta plasma in which the transport of heat and momentum occurs primarily along magnetic-field lines. Anisotropic heat conduction allows convective instabilities to be driven by temperature gradients of either sign, the magnetothermal instability (MTI) in the outskirts of non-isothermal clusters and the heat-flux buoyancy-driven instability (HBI) in their cooling cores. We employ the Athena MHD code to investigate the nonlinear evolution of these instabilities, self-consistently including the effects of anisotropic viscosity (i.e. Braginskii pressure anisotropy), anisotropic conduction, and radiative cooling. We highlight the importance of the microscale instabilities that inevitably accompany and regulate the pressure anisotropies generated by the HBI and MTI. We find that, in all but the innermost regions of cool-core clusters, anisotropic viscosity significantly impairs the ability of the HBI to reorient magnetic-field lines orthogonal to the temperature gradient. Thus, while radio-mode feedback appears necessary in the central few tens of kpc, conduction may be capable of offsetting radiative losses throughout most of a cool core over a significant fraction of the Hubble time. Magnetically-aligned cold filaments are then able to form by local thermal instability. Viscous dissipation during the formation of a cold filament produces accompanying hot filaments, which can be searched for in deep Chandra observations of nearby cool-core clusters. In the case of the MTI, anisotropic viscosity maintains the coherence of magnetic-field lines over larger distances than in the inviscid case, providing a natural lower limit for the scale on which the field can fluctuate freely. In the nonlinear state, the magnetic field exhibits a folded structure in which the field-line curvature and field strength are anti-correlated.

preprint2012arXiv

Constraints on Compton-thick winds from black hole accretion disks: can we see the inner disk?

Strong evidence is emerging that winds can be driven from the central regions of accretion disks in both active galactic nuclei (AGN) and Galactic black hole binaries (GBHBs). Direct evidence for highly-ionized, Compton-thin inner-disk winds comes from observations of blueshifted (v~0.05-0.1c) iron-K X-ray absorption lines. However, it has been suggested that the inner regions of black hole accretion disks can also drive Compton-thick winds --- such winds would enshroud the inner disk, preventing us from seeing direct signatures of the accretion disk (i.e. the photospheric thermal emission, or the Doppler/gravitationally broadened iron K-alpha line). Here, we show that, provided the source is sub-Eddington, the well-established wind driving mechanisms fail to launch a Compton-thick wind from the inner disk. For the accelerated region of the wind to be Compton-thick, the momentum carried in the wind must exceed the available photon momentum by a factor of at least 2/lambda, where lambda is the Eddington ratio of the source, ruling out radiative acceleration unless the source is very close to the Eddington limit. Compton-thick winds also carry large mass-fluxes, and a consideration of the connections between the wind and the disk show this to be incompatible with magneto-centrifugal driving. Finally, thermal driving of the wind is ruled out on the basis of the large Compton-radii that typify black hole systems. In the absence of some new acceleration mechanism, we conclude that the inner regions of sub-Eddington accretion disks around black holes are indeed naked.

preprint2012arXiv

General Relativistic Simulations of Magnetized Plasmas around Merging Supermassive Black Holes

Coalescing supermassive black hole binaries are produced by the mergers of galaxies and are the most powerful sources of gravitational waves accessible to space-based gravitational observatories. Some such mergers may occur in the presence of matter and magnetic fields and hence generate an electromagnetic counterpart. In this Letter, we present the first general relativistic simulations of magnetized plasma around merging supermassive black holes using the general relativistic magnetohydrodynamic code Whisky. By considering different magnetic field strengths, going from non-magnetically dominated to magnetically dominated regimes, we explore how magnetic fields affect the dynamics of the plasma and the possible emission of electromagnetic signals. In particular we observe a total amplification of the magnetic field of ~2 orders of magnitude which is driven by the accretion onto the binary and that leads to much stronger electromagnetic signals, more than a factor of 10^4 larger than comparable calculations done in the force-free regime where such amplifications are not possible.

preprint2012arXiv

Prompt Emission from Tidal Disruptions of White Dwarfs by Intermediate Mass Black Holes

We present a qualitative picture of prompt emission from tidal disruptions of white dwarfs (WD) by intermediate mass black holes (IMBH). The smaller size of an IMBH compared to a supermassive black hole and a smaller tidal radius of a WD disruption lead to a very fast event with high peak luminosity. Magnetic field is generated in situ following the tidal disruption, which leads to effective accretion. Since large-scale magnetic field is also produced, geometrically thick super-Eddington inflow leads to a relativistic jet. The dense jet possesses a photosphere, which emits quasi-thermal radiation in soft X-rays. The source can be classified as a long low-luminosity gamma-ray burst (ll-GRB). Tidal compression of a WD causes nuclear ignition, which is observable as an accompanying supernova. We suggest that GRB060218 and SN2006aj is such a pair of ll-GRB and supernova. We argue that in a flux-limited sample the disruptions of WDs by IMBHs are more frequent then the disruptions of other stars by IMBHs.

preprint2012arXiv

The Black Hole Spin and Soft X-ray Excess of the Luminous Seyfert Galaxy Fairall 9

We present an analysis of all XMM-Newton and Suzaku X-ray spectra of the nearby luminous Seyfert galaxy Fairall 9. Confirming previous analyses, we find robust evidence for a broad iron line associated with X-ray reflection from the innermost accretion disk. By fitting a spectral model that includes a relativistically ionized reflection component, we examine the constraints on the inclination of the inner accretion disk and the black hole spin, and the complications introduced by the presence of a photoionized emission line system. Employing multi-epoch fitting, we attempt to obtain robust and concordant measures of the accretion disk parameters. We also clearly see a soft X-ray excess in Fairall 9. During certain epochs, the soft excess can be described with same disk reflection component that produces the iron line. However, there are epochs where an additional soft component is required. This can be attributed to either an additional highly-ionized, strongly blurred disk reflection component, or a new X-ray continuum component.

preprint2012arXiv

The Central Engine Structure of 3C120: Evidence for a Retrograde Black Hole or a Refilling Accretion Disk

The broad-line radio galaxy 3C120 is a powerful source of both X-ray and radio emission including superluminal jet outflows. We report on our reanalysis of 160 ks of Suzaku data taken in 2006, previously examined by Kataoka et al. (2007). Spectral fits to the XIS and HXD/PIN data over a range of 0.7-45 keV reveal a well-defined iron K line complex with a narrow Ka core and relativistically broadened features consistent with emission from the inner regions of the accretion disk. Furthermore, the inner region of the disk appears to be truncated with an inner radius of r_in = 11.7^{+3.5}_{-5.2} r_g. If we assume that fluorescent iron line features terminate at the inner-most stable circular orbit (ISCO), we measure a black hole spin of a < -0.1 at a 90% confidence level. A rapidly spinning prograde black hole (a > 0.8) can be ruled out at the 99% confidence level. Alternatively, the disk may be truncated well outside of the ISCO of a rapid prograde hole. The most compelling scenario is the possibility that the inner regions of the disk were destroyed/ejected by catastrophic instabilities just prior to the time these observations were made.

preprint2012arXiv

X-ray dips in the Seyfert Galaxy Fairall 9: Compton-thick "comets" or a failed radio galaxy?

We investigate the spectral variability of the Seyfert galaxy Fairall 9 using almost 6 years of monitoring with the Rossi X-ray Timing Explorer (RXTE) with an approximate time resolution of 4 days. We discover the existence of pronounced and sharp dips in the X-ray flux, with a rapid decline of the 2--20 keV flux of a factor 2 or more followed by a recovery to pre-dip fluxes after ~10 days . These dips skew the flux distribution away from the commonly observed log-normal distribution. Dips may result from the eclipse of the central X-ray source by broad line region (BLR) clouds, as has recently been found in NGC 1365 and Mrk 766. Unlike these other examples, however, the clouds in Fairall 9 would need to be Compton-thick, and the non-dip state is remarkably free of any absorption features. A particularly intriguing alternative is that the accretion disk is undergoing the same cycle of disruption/ejection as seen in the accretion disks of broad line radio galaxies (BLRGs) such as 3C120 but, for some reason, fails to create a relativistic jet. This suggests that a detailed comparison of Fairall 9 and 3C120 with future high-quality data may hold the key to understanding the formation of relativistic jets in AGN.

preprint2011arXiv

A Chandra Study of the Radio Galaxy NGC 326: Wings, Outburst History, and AGN Feedback

NGC 326 is one of the most prominent X- or Z-shaped radio galaxies (XRGs/ZRGs) and has been the subject of several studies attempting to explain its morphology through either fluid motions or reorientation of the jet axis. We examine a 100 ks archival Chandra exposure and find several features associated with the radio galaxy: a high-temperature front that may indicate a shock, high-temperature knots around the rim of the radio emission, and a cavity associated with the eastern wing of the radio galaxy. A reasonable interpretation of these features in light of the radio data allows us to reconstruct the history of the AGN outbursts. The active outburst was likely once a powerful radio source which has since decayed, and circumstantial evidence favors reorientation as the means to produce the wings. Because of the obvious interaction between the radio galaxy and the ICM and the wide separation between the active lobes and wings, we conclude that XRGs are excellent sources in which to study AGN feedback in galaxy groups by measuring the heating rates associated with both active and passive heating mechanisms.

preprint2011arXiv

Fermi/LAT Observations of Swift/BAT Seyferts: on the Contribution of Radio-quiet AGN to the Extragalactic Gamma-ray Background

We present the analysis of 2.1 years of Fermi/LAT data on 491 Seyfert galaxies detected by the Swift/BAT survey. Only the two nearest objects, NGC 1068 and NGC 4945, which were identified in the Fermi First-year Catalog, are detected. Using the Swift/BAT and radio 20 cm fluxes, we define a new radio-loudness parameter $R_{X, BAT}$ where radio loud objects have $\log R_{X, BAT} > - 4.7$. Based on this parameter, only radio loud sources are detected by Fermi/LAT. An upper limit to the flux of the undetected sources is derived to be $\sim 2 \times 10^{-11}$ photons cm$^{-2}$ s$^{-1}$, approximately seven times lower than the observed flux of NGC 1068. Assuming a median redshift of 0.031, this implies an upper limit to the $γ$-ray (1--100 GeV) luminosity of $\lesssim 3 \times 10^{41}$ erg s$^{-1}$. In addition, we identified 120 new Fermi/LAT sources near the Swift/BAT Seyferts with significant Fermi/LAT detections. A majority of these objects do not have \bat counterparts, but their possible optical counterparts include blazars, FSRQs, and quasars.

preprint2011arXiv

Global Simulations of Accretion Disks I: Convergence and Comparisons with Local Models

Grid-based magnetohydrodynamic (MHD) simulations have proven invaluable for the study of astrophysical accretion disks. However, the fact that angular momentum transport in disks is mediated by MHD turbulence (with structure down to very small scales) raises the concern that the properties of the modeled accretion disks are affected by the finite numerical resolution of the simulation. By implementing an orbital advection algorithm into the Athena code in cylindrical geometry, we have performed a set of global (but unstratified) Newtonian disk simulations extending up to resolutions previously unattained. We study the convergence of these models as a function of spatial resolution and initial magnetic field geometry. The usual viscosity parameter ($α$) or the ratio of thermal-to-magnetic pressure ($β$) are found to be poor diagnostics of convergence, whereas the average tilt angle of the magnetic field in the $(r,ϕ)$-plane is a very good diagnostic of convergence. We suggest that this is related to the saturation of the MHD turbulence via parasitic modes of the magnetorotational instability. Even in the case of zero-net magnetic flux, we conclude that our highest resolution simulations (with 32-zones and 64-zones per vertical scale height) have achieved convergence. Our global simulations reach resolutions comparable to those used in local, shearing box models of MHD disk turbulence. We find that the saturation predictors derived from local simulations correspond well to the instantaneous correlations between local flux and stress found in our global simulations. However, the conservation of magnetic flux implicit in local models is not realized in our global disks. Thus, the magnetic connectivity of an accretion disk represents physics that is truly global and cannot be captured in any ab-initio local model.

preprint2011arXiv

Hydrodynamic Models of Radio Galaxy Morphology: Winged and X-shaped Sources

We present three-dimensional hydrodynamic models of radio galaxies interacting with initially relaxed hot atmospheres and explore the significant off-axis radio lobe structures which result under certain conditions. With a focus on the "winged" and "X-shaped" radio galaxy population, we confirm the importance of observed trends such as the connection of wing formation with jets co-aligned with the major axis of the surrounding atmosphere. These wings are formed substantially by the deflection of lobe plasma flowing back from the hot spots (backflow) and develop in two stages: supersonic expansion of an overpressured cocoon at early times followed by buoyant expansion at later times. We explore a limited parameter space of jet and atmosphere properties and find that the most prominent wings are produced when a decaying jet is injected into a small, dense, highly elliptical atmosphere. On the basis of this search, we argue that the deflection of backflow by gradients in the hot atmosphere is a strong candidate for forming observed wings but must work in tandem with some other mechanism for forming the initial wing channels. Our models indicate that lobe interaction with the hot atmosphere may play a dominant role in shaping the morphology of radio galaxies.

preprint2011arXiv

Low-Frequency Oscillations in Global Simulations of Black Hole Accretion

We have identified the presence of large-scale, low-frequency dynamo cycles in a long-duration, global, magnetohydrodynamic (MHD) simulation of black hole accretion. Such cycles had been seen previously in local shearing box simulations, but we discuss their evolution over 1,500 inner disk orbits of a global pi/4 disk wedge spanning two orders of magnitude in radius and seven scale heights in elevation above/below the disk midplane. The observed cycles manifest themselves as oscillations in azimuthal magnetic field occupying a region that extends into a low-density corona several scale heights above the disk. The cycle frequencies are ten to twenty times lower than the local orbital frequency, making them potentially interesting sources of low-frequency variability when scaled to real astrophysical systems. Furthermore, power spectra derived from the full time series reveal that the cycles manifest themselves at discrete, narrow-band frequencies that often share power across broad radial ranges. We explore possible connections between these simulated cycles and observed low-frequency quasi-periodic oscillations (LFQPOs) in galactic black hole binary systems, finding that dynamo cycles have the appropriate frequencies and are located in a spatial region associated with X-ray emission in real systems. Derived observational proxies, however, fail to feature peaks with RMS amplitudes comparable to LFQPO observations, suggesting that further theoretical work and more sophisticated simulations will be required to form a complete theory of dynamo-driven LFQPOs. Nonetheless, this work clearly illustrates that global MHD dynamos exhibit quasi-periodic behavior on timescales much longer than those derived from test particle considerations.

preprint2010arXiv

A Deep Chandra Observation of the X-shaped Radio Galaxy 4C +00.58: A Candidate for Merger-induced Reorientation?

Although rapid reorientation of a black hole spin axis (lasting less than a few Myr) has been suggested as a mechanism for the formation of wings in X-shaped radio galaxies (XRGs), to date no convincing case of reorientation has been found in any XRG. Alternative wing formation models such as the hydrodynamic backflow models are supported by observed trends indicating that XRGs form preferentially with jets aligned along the major axis of the surrounding medium and wings along the minor axis. In this Letter, we present a deep Chandra observation of 4C +00.58, an oddball XRG with its jet oriented along the minor axis. By using the X-ray data in tandem with available radio and optical data, we estimate relevant timescales with which to evaluate wing formation models. The hydrodynamic models have difficulty explaining the long wings, whereas the presence of X-ray cavities (suggesting jet activity along a prior axis) and a potential stellar shell (indicating a recent merger) favor a merger-induced reorientation model.

preprint2010arXiv

Connections Between Local and Global Turbulence in Accretion Disks

We analyze a suite of global magnetohydrodynamic (MHD) accretion disk simulations in order to determine whether scaling laws for turbulence driven by the magnetorotational instability, discovered via local shearing box studies, are globally robust. The simulations model geometrically-thin disks with zero net magnetic flux and no explicit resistivity or viscosity. We show that the local Maxwell stress is correlated with the self-generated local vertical magnetic field in a manner that is similar to that found in local simulations. Moreover, local patches of vertical field are strong enough to stimulate and control the strength of angular momentum transport across much of the disk. We demonstrate the importance of magnetic linkages (through the low-density corona) between different regions of the disk in determining the local field, and suggest a new convergence requirement for global simulations -- the vertical extent of the corona must be fully captured and resolved. Finally, we examine the temporal convergence of the average stress, and show that an initial long-term secular drift in the local flux-stress relation dies away on a time scale that is consistent with turbulent mixing of the initial magnetic field.

preprint2010arXiv

Radiative and dynamic stability of a dilute plasma

We analyze the linear stability of a dilute, hot plasma, taking into account the effects of stratification and anisotropic thermal conduction. The work is motivated by attempts to understand the dynamics of the intracluster medium in galaxy clusters. We show that magnetic field configurations that nominally stabilize either the heat-flux driven buoyancy instability (associated with a positive thermal gradient) or the magnetothermal instability (negative thermal gradient) can lead to previously unrecognized g-mode overstabilities. The driving source of the overstability is either radiative cooling (positive temperature gradient) or the heat flux itself (negative temperature gradient). While the implications of these overstabilities have yet to be explored, we speculate that the cold fronts observed in many relaxed galaxy clusters may be related to their non-linear evolution.

preprint2009arXiv

Constraints on the absorption-dominated model for the X-ray spectrum of MCG-6-30-15

Complexities in the X-ray spectrum of the nearby Seyfert 1.2 galaxy MCG-6-30-15 are commonly interpreted in terms of a broad iron line and the associated Compton reflection hump from the innermost relativistic regions of an accretion disk around a rapidly spinning black hole. However, an alternative model has recently been proposed in which these spectral features are caused entirely by complex (ionized and partial-covering) absorption. By considering the fluorescent emission that must accompany photoelectric absorption, we show that the absorption-dominated model over-predicts the 6.4keV iron line flux unless the marginally Compton-thick absorber responsible for the hard X-ray hump satisfies very restrictive geometric constraints. In the absence of a specific model that both obeys these geometrical constraints and is physically-plausible, the relativistic-reflection model is favoured.

preprint2009arXiv

Modeling Flows Around Merging Black Hole Binaries

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal.

preprint2009arXiv

Relativistic Broadening of Iron Emission Lines in a Sample of AGN

We present a uniform X-ray spectral analysis of eight type-1 active galactic nuclei (AGN) that have been previously observed with relativistically broadened iron emission lines. Utilizing data from the XMM-Newton European Photon Imaging Camera (EPIC-pn) we carefully model the spectral continuum, taking complex intrinsic absorption and emission into account. We then proceed to model the broad Fe K feature in each source with two different accretion disk emission line codes, as well as a self-consistent, ionized accretion disk spectrum convolved with relativistic smearing from the inner disk. Comparing the results, we show that relativistic blurring of the disk emission is required to explain the spectrum in most sources, even when one models the full reflection spectrum from the photoionized disk.

preprint2009arXiv

The Chandra View of Nearby X-shaped Radio Galaxies

We present new and archival Chandra X-ray Observatory observations of X-shaped radio galaxies within z < 0.1 alongside a comparison sample of normal double-lobed FR I and II radio galaxies. By fitting elliptical distributions to the observed diffuse hot X-ray emitting atmospheres, we find that the ellipticity and the position angle of the hot gas follows that of the stellar light distribution for radio galaxy hosts in general. Moreover, compared to the control sample, we find a strong tendency for X-shaped morphology to be associated with wings directed along the minor axis of the hot gas distribution. Taken at face value, this result favors the hydrodynamic backflow models for the formation of X-shaped radio galaxies which naturally explain the geometry; the merger-induced rapid reorientation models make no obvious prediction about orientation.

preprint2008arXiv

On the time variability of geometrically-thin black hole accretion disks I : the search for modes in simulated disks

We present a detailed temporal analysis of a set of hydrodynamic and magnetohydrodynamic (MHD) simulations of geometrically-thin (h/r~0.05) black hole accretion disks. The black hole potential is approximated by the Paczynski-Wiita pseudo-Newtonian potential. In particular, we use our simulations to critically assess two widely discussed models for high-frequency quasi-periodic oscillations, global oscillation modes (diskoseismology) and parametric resonance instabilities. We find that initially disturbed hydrodynamic disks clearly display the trapped global g-mode oscillation predicted by linear perturbation theory. In contrast, the sustained turbulence produced in the simulated MHD disks by the magneto-rotational instability does not excite these trapped g-modes. We cannot say at present whether the MHD turbulence actively damps the hydrodynamic g-mode. Our simulated MHD disks also fail to display any indications of a parametric resonance instability between the vertical and radial epicyclic frequencies. On the other hand, we do see characteristic frequencies at any given radius in the disk corresponding to local acoustic waves. We also conduct a blind search for any quasi-periodic oscillation in a proxy lightcurve based on the instantaneous mass accretion rate of the black hole, and place an upper limit of 2% on the total power in any such feature. We highlight the importance of correcting for secular changes in the simulated accretion disk when performing temporal analyses.

preprint2008arXiv

On the Time Variability of Geometrically-Thin Black Hole Accretion Disks II: Viscosity-Induced Global Oscillation Modes in Simulated Disks

We examine the evolution and influence of viscosity-induced diskoseismic modes in simulated black hole accretion disks. Understanding the origin and behavior of such oscillations will help us to evaluate their potential role in producing astronomically observed high-frequency quasi-periodic oscillations in accreting black hole binary systems. Our simulated disks are geometrically-thin with a constant half-thickness of five percent the radius of the innermost stable circular orbit. A pseudo-Newtonian potential reproduces the relevant effects of general relativity, and an alpha-model viscosity achieves angular momentum transport and the coupling of orthogonal velocity components in an otherwise ideal hydrodynamic numerical treatment. We find that our simulated viscous disks characteristically develop and maintain trapped global mode oscillations with properties similar to those expected of trapped g-modes and inner p-modes in a narrow range of frequencies just below the maximum radial epicyclic frequency. Although the modes are driven in the inner portion of the disk, they generate waves that propagate at the trapped mode frequencies out to larger disk radii. This finding is contrasted with the results of global magnetohydrodynamic disk simulations, in which such oscillations are not easily identified. Such examples underscore fundamental physical differences between accretion systems driven by the magneto-rotational instability and those for which alpha viscosity serves as a proxy for the physical processes that drive accretion, and we explore potential approaches to the search for diskoseismic modes in full magnetohydrodynamic disks.

preprint2008arXiv

Probing the accretion disk and central engine structure of NGC4258 with Suzaku and XMM-Newton observations

[abridged] We present an X-ray study of the low-luminosity active galactic nucleus (AGN) in NGC4258 using data from Suzaku, XMM-Newton, and the Swift/BAT survey. We find that signatures of X-ray reprocessing by cold gas are very weak in the spectrum of this Seyfert-2 galaxy; a weak, narrow fluorescent-Kalpha emission line of cold iron is robustly detected in both the Suzaku and XMM-Newton spectra but at a level much below that of most other Seyfert-2 galaxies. We conclude that the circumnuclear environment of this AGN is very "clean" and lacks the Compton-thick obscuring torus of unified Seyfert schemes. From the narrowness of the iron line, together with evidence for line flux variability between the Suzaku and XMM-Newton observations, we constrain the line emitting region to be between $3\times 10^3r_g$ and $4\times 10^4r_g$ from the black hole. We show that the observed properties of the iron line can be explained if the line originates from the surface layers of a warped accretion disk. In particular, we present explicit calculations of the expected iron line from a disk warped by Lens-Thirring precession from a misaligned central black hole. Finally, the Suzaku data reveal clear evidence for large amplitude 2-10keV variability on timescales of 50ksec as well as smaller amplitude flares on timescales as short as 5-10ksec. If associated with accretion disk processes, such rapid variability requires an origin in the innermost regions of the disk ($r\approx 10r_g$ or less).

preprint2006arXiv

AGN Heating of Cooling Flow Clusters: Problems with 3D Hydrodynamic Models

We have performed a set of high resolution, three dimensional simulations of a jetted AGN embedded in a relaxed cooling cluster. To the best of our knowledge, these are the first simulation to include both full jet dynamics and a feedback model. These ideal hydrodynamic simulations show that, although there is enough energy present to offset cooling on average, the jet heating is not spatially deposited in a way that can prevent catastrophic cooling of the cluster.

preprint2000arXiv

Simulations of accretion flows crossing the last stable orbit

We use three dimensional magnetohydrodynamic simulations, in a pseudo-Newtonian potential, to study geometrically thin accretion disc flows crossing the marginally stable circular orbit around black holes. We concentrate on vertically unstratified and isothermal disk models, but also consider a model that includes stratification. In all cases, we find that the sonic point lies just inside the last stable orbit, with a modest increase in the importance of magnetic field energy, relative to the thermal energy, observed interior to this radius. The time-averaged gradient of the specific angular momentum of the flow, (dl/dr), is close to zero within the marginally stable orbit, despite the presence of large fluctuations and continuing magnetic stress in the plunging region. The result that the specific angular momentum is constant within the last stable orbit is in general agreement with traditional disk models computed using a zero-torque boundary condition at the last stable orbit.

Christopher S. Reynolds

What is connected

Connect this record

See the researcher in context

Building this map preview

56 published item(s)

A Disc Reflection Model for Ultra-Soft Narrow-Line Seyfert 1 Galaxies

Acoustic waves and g-mode turbulence as energy carriers in a viscous intracluster medium

Black Hole Spin Measurements Based on a Thin Disc Model with Finite Thickness I. An example study of MCG-06-30-15

Evidence for a moderate spin from X-ray reflection of the high-mass supermassive black hole in the cluster-hosted quasar H1821+643

How do Magnetic Field Models Affect Astrophysical Limits on Light Axion-like Particles? An X-ray Case Study with NGC 1275

Probing the Extent of Fe K$α$ Emission in Nearby Active Galactic Nuclei using Multi-Order Analysis of Chandra High Energy Transmission Grating Data

XMM-Newton Observations of the Narrow-Line Seyfert 1 Galaxy IRAS 13224$-$3809 : X-ray Spectral Analysis II

New constraints on light Axion-Like Particles using Chandra Transmission Grating Spectroscopy of the powerful cluster-hosted quasar H1821+643

Lense-Thirring Precession of Misaligned Discs I

Probing the Milky Way's Dark Matter Halo for the 3.5 keV Line

Revealing the broad iron Kalpha line in Cygnus X-1 through simultaneous XMM-Newton, RXTE, and INTEGRAL observations

The view of AGN-host alignment via reflection spectroscopy

Accretion disk dynamo as the trigger for X-ray binary state transitions

Inefficient Driving of Bulk Turbulence by Active Galactic Nuclei in a Hydrodynamic Model of the Intracluster Medium

NuSTAR observations of the powerful radio-galaxy Cygnus A

Unwrapping the X-ray Spectra of Active Galactic Nuclei

An XMM-Newton View of the Radio Galaxy 3C 411

Modeling Hot Gas Flow in the Low-Luminosity Active Galactic Nucleus of NGC3115

Non-Linear Dynamics of Accretion Disks with Stochastic Viscosity

The hard X-ray perspective on the soft X-ray excess

The role of Compton heating in radiation-regulated accretion on to black holes

The X-ray spectrum of the cooling-flow quasar H1821+643 : A massive black hole feeding off the intracluster medium

An X-Ray View of the Jet-Cycle in the Radio Loud AGN 3C120

GRB060218 as a Tidal Disruption of a White Dwarf by an Intermediate Mass Black Hole

Irradiation of an Accretion Disc by a Jet: General Properties and Implications for Spin Measurements of Black Holes

Isotropic Heating of Galaxy Cluster Cores via Rapidly Reorienting AGN Jets

Measuring Black Hole Spin using X-ray Reflection Spectroscopy

Rayleigh-Taylor instability of ionization front around black holes

Tackling the soft X-ray excess in AGN with variability studies

The Ionized Absorber and Nuclear Environment of IRAS 13349+2438: Multi-wavelength insights from coordinated Chandra HETGS, HST STIS, HET, and Spitzer IRS

The Spin of Supermassive Black Holes

What is on Tap? The Role of Spin in Compact Objects and Relativistic Jets

Buoyancy Instabilities in a Weakly Collisional Intracluster Medium

Constraints on Compton-thick winds from black hole accretion disks: can we see the inner disk?

General Relativistic Simulations of Magnetized Plasmas around Merging Supermassive Black Holes

Prompt Emission from Tidal Disruptions of White Dwarfs by Intermediate Mass Black Holes

The Black Hole Spin and Soft X-ray Excess of the Luminous Seyfert Galaxy Fairall 9

The Central Engine Structure of 3C120: Evidence for a Retrograde Black Hole or a Refilling Accretion Disk

X-ray dips in the Seyfert Galaxy Fairall 9: Compton-thick "comets" or a failed radio galaxy?

A Chandra Study of the Radio Galaxy NGC 326: Wings, Outburst History, and AGN Feedback

Fermi/LAT Observations of Swift/BAT Seyferts: on the Contribution of Radio-quiet AGN to the Extragalactic Gamma-ray Background

Global Simulations of Accretion Disks I: Convergence and Comparisons with Local Models

Hydrodynamic Models of Radio Galaxy Morphology: Winged and X-shaped Sources

Low-Frequency Oscillations in Global Simulations of Black Hole Accretion

A Deep Chandra Observation of the X-shaped Radio Galaxy 4C +00.58: A Candidate for Merger-induced Reorientation?

Connections Between Local and Global Turbulence in Accretion Disks

Radiative and dynamic stability of a dilute plasma

Constraints on the absorption-dominated model for the X-ray spectrum of MCG-6-30-15

Modeling Flows Around Merging Black Hole Binaries

Relativistic Broadening of Iron Emission Lines in a Sample of AGN

The Chandra View of Nearby X-shaped Radio Galaxies

On the time variability of geometrically-thin black hole accretion disks I : the search for modes in simulated disks

On the Time Variability of Geometrically-Thin Black Hole Accretion Disks II: Viscosity-Induced Global Oscillation Modes in Simulated Disks

Probing the accretion disk and central engine structure of NGC4258 with Suzaku and XMM-Newton observations

AGN Heating of Cooling Flow Clusters: Problems with 3D Hydrodynamic Models

Simulations of accretion flows crossing the last stable orbit