Reprocessed emission from warped accretion discs with application to X-ray iron line profiles

Fluorescent iron line profiles currently provide the best diagnostic for engine geometries of active galactic nuclei (AGN). Here we construct a method for calculating the relativistic iron line profile from an arbitrarily warped accretion disc, illuminated from above and below by hard X-ray sources. This substantially generalizes previous calculations of reprocessing by accretion discs by including non-axisymmetric effects. We include a relativistic treatment of shadowing by ray-tracing photon paths along Schwarzschild geodesics. We apply this method to two classes of warped discs, and generate a selection of resulting line profiles. New profile features include a time-varying line profile if the warp precesses about the disc, profile differences between 'twisted' and 'twist-free' warps and the possibility of steeper red and softer blue fall-offs than for flat discs. We discuss some qualitative implications of the line profiles in the context of Type I and II Seyfert AGN and other sources.     

The effect of stellar feedback on the formation and evolution of gas and dust tori in AGN

Recently, the existence of geometrically thick dust structures in active galactic nuclei (AGN) has been directly proven with the help of interferometric methods in the mid-infrared. The observations are consistent with a two-component model made up of a geometrically thin and warm central disc, surrounded by a colder, fluffy torus component. Within the framework of an exploratory study, we investigate one possible physical mechanism, which could produce such a structure, namely the effect of stellar feedback from a young nuclear star cluster on the interstellar medium in centres of AGN. The model is realized by numerical simulations with the help of the hydrodynamics code tramp . We follow the evolution of the interstellar medium by taking discrete mass-loss and energy ejection due to stellar processes, as well as optically thin radiative cooling into account. In a post-processing step, we calculate observable quantities like spectral energy distributions (SEDs) and surface brightness distributions with the help of the radiative transfer code mc3d . The interplay between injection of mass, supernova explosions and radiative cooling leads to a two-component structure made up of a cold geometrically thin, but optically thick and very turbulent disc residing in the vicinity of the angular momentum barrier, surrounded by a filamentary structure. The latter consists of cold long radial filaments flowing towards the disc and a hot tenuous medium in between, which shows both inwards and outwards directed motions. With the help of this modelling, we are able to reproduce the range of observed neutral hydrogen column densities of a sample of Seyfert galaxies as well as the relation between them and the strength of the silicate 10 μm spectral feature. Despite being quite crude, our mean Seyfert galaxy model is even able to describe the SEDs of two intermediate type Seyfert galaxies observed with the Spitzer Space Telescope .     

On the hard X-ray spectra of radio-loud active galaxies

Over the last few years X-ray observations of broad-line radio galaxies (BLRGs) by ASCA , RXTE and BeppoSAX have shown that these objects seem to exhibit weaker X-ray reflection features (such as the iron K α line) than radio-quiet Seyferts. This has lead to speculation that the optically thick accretion disc in radio-loud active galactic nuclei (AGN) may be truncated to an optically thin flow in the inner regions of the source. Here, we propose that the weak reflection features are a result of reprocessing in an ionized accretion disc. This would alleviate the need for a change in accretion geometry in these sources. Calculations of reflection spectra from an ionized disc for situations expected in radio-loud AGN (high accretion rate, moderate-to-high black hole mass) predict weak reprocessing features. This idea was tested by fitting the ASCA spectrum of the bright BLRG 3C 120 with the constant density ionized disc models of Ross & Fabian. A good fit was found with an ionization parameter of   ξ ∼4000 erg cm s^-1  and the reflection fraction fixed at unity. If observations of BLRGs by XMM-Newton show evidence for ionized reflection then this would support the idea that a high accretion rate is likely required to launch powerful radio jets.     

Iron line profiles and self-shadowing from relativistic thick accretion discs

We present Fe Kα line profiles from and images of relativistic discs with finite thickness around a rotating black hole using a novel code. The line is thought to be produced by iron fluorescence of a relatively cold X-ray-illuminated material in the innermost parts of the accretion disc and provides an excellent diagnostic of accretion flows in the vicinity of black holes. Previous studies have concentrated on the case of a thin, Keplerian accretion disc. This disc must become thicker and sub-Keplerian with increasing accretion rates. These can affect the line profiles and in turn can influence the estimation of the accretion disc and black hole parameters from the observed line profiles. We here embark on, for the first time, a fully relativistic computation which offers key insights into the effects of geometrical thickness and the sub-Keplerian orbital velocity on the line profiles. We include all relativistic effects such as frame-dragging, Doppler boost, time dilation, gravitational redshift and light bending. We find that the separation and the relative height between the blue and red peaks of the line profile diminish as the thickness of the disc increases. This code is also well suited to produce accretion disc images. We calculate the redshift and flux images of the accretion disc and find that the observed image of the disc strongly depends on the inclination angle. The self-shadowing effect appears remarkable for a high inclination angle, and leads to the black hole shadow being in this case, completely hidden by the disc itself.     

A cloud model of active galactic nuclei: the iron Kα line diagnostics

The origin and the profile of the iron K α line from active galactic nuclei are studied. The model with a quasi-spherical distribution of clouds around a compact galactic nucleus is explored, in which the intrinsic spectral feature around 6–7 keV is substantially affected by Comptonization of primary X-rays. The predicted spectral profiles are influenced mainly by the motion of the clouds in the strong gravitational field of the centre, and by the orientation of the source with respect to the observer. The case of the Seyfert galaxy MCG–6-30-15 is examined in more detail and free parameters of the model are determined. A robust result is obtained: acceptable fits to the observed profile correspond to a spherical distribution of the clouds with strong self-obscuration effects and with a rather high value of ionization parameter. This is a different situation from that considered in the model of a cold illuminated disc. Nevertheless, geometrically flat (disc-type) configurations are not rejected.     

The evolution of misaligned accretion discs and spinning black holes

In this paper, we consider the process of alignment of a spinning black hole and a surrounding misaligned accretion disc. We use a simplified set of equations, that describe the evolution of the system in the case where the propagation of warping disturbances in the accretion disc occurs diffusively, a situation likely to be common in the thin discs in active galactic nuclei (AGN). We also allow the direction of the hole spin to move under the action of the disc torques. In such a way, the evolution of the hole–disc system is computed self-consistently. We consider a number of different situations and we explore the relevant parameter range, by varying the location of the warp radius R _w and the propagation speed of the warp. We find that the dissipation associated with the twisting of the disc results in a large increase in the accretion rate through the disc, so that AGN accreting from a misaligned disc are likely to be significantly more luminous than those accreting from a flat disc. We compute explicitly the time-scales for the warping of the disc and for the alignment process and compare our results with earlier estimates based on simplified steady-state solutions. We also confirm earlier predictions that, under appropriate circumstances, accretion can proceed in a counter-aligned fashion, so that the accreted material will spin-down the hole, rather than spinning it up. Our results have implication in a number of different observational features of AGN such as the orientation and shape of jets, the shape of X-ray iron lines and the possibility of obscuration and absorption of X-ray by the outer disc as well as the general issue of the spin history of black holes during their growth.     

RE J2248−511: not all variable, ultrasoft, X-ray AGN have narrow Balmer lines

We present ASCA data on RE J2248−511, extending existing optical and soft X-ray coverage to 10 keV, and monitoring the soft component. These data show that, despite a very strong ultrasoft X-ray excess below 0.3 keV and a soft 0.3–2 keV spectral index in earlier ROSAT data, the hard X-ray spectrum ( α ∼−0.8; 0.6–10 keV) is typical of type 1 active galactic nuclei (AGN), and the soft component has since disappeared. Optical data taken at two different epochs show that the big blue bump is also highly variable. The strength of the ultrasoft X-ray component and the extreme variability in RE J2248−511 are reminiscent of the behaviour observed in many narrow line Seyfert 1s (NLS1s). However, the high-energy end of the ROSAT spectrum, the ASCA spectrum and the Balmer line full widths at half maximum of ∼3000 km s⁻¹ in RE J2248−511 are typical of normal Seyfert 1 AGN.
The change in the soft X-ray spectrum as observed in the ROSAT and ASCA data is consistent with the behaviour of Galactic Black Hole Candidates (GBHCs) as they move from a high to a low state, i.e. a fall in the ultrasoft component and a hardening of the X-ray continuum. This GBHC analogy has also been proposed for NLS1s. Alternatively, the variability may be caused by opacity changes in a hot, optically thin corona which surrounds a cold, dense accretion disc; this was first suggested by Guainazzi et al. for 1H 0419−577, an object which shows remarkably similar properties to RE J2248−511.     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

Low-Luminosity Accretion in Black Hole X-Ray Binaries and Active Galactic Nuclei

At luminosities below a few percent of Eddington, accreting black holes switch to a hard spectral state which is very different from the soft blackbody-like spectral state that is found at higher luminosities. The hard state is well-described by a two-temperature, optically thin, geometrically thick, advection-dominated accretion flow (ADAF) in which the ions are extremely hot (up to 10¹² K near the black hole), the electrons are also hot (∼10^9−10.5 K), and thermal Comptonization dominates the X-ray emission. The radiative efficiency of an ADAF decreases rapidly with decreasing mass accretion rate, becoming extremely low when a source reaches quiescence. ADAFs are expected to have strong outflows, which may explain why relativistic jets are often inferred from the radio emission of these sources. It has been suggested that most of the X-ray emission also comes from a jet, but this is less well established.     

Effects of inflationary bubbles on the polarization and temperature anisotropies of the cosmic microwave background

The broad X-ray iron line, detected in many active galactic nuclei, is likely to be produced by fluorescence from the X-ray-illuminated central parts of an accretion disc close to a supermassive black hole. The time-averaged shape of the line can be explained most naturally by a combination of special and general relativistic effects. Such line profiles contain information about the black hole spin and the accretion disc, as well as the geometry of the emitting region, and may help to test general relativity in the strong gravity regime. In this paper we embark on the computation of the temporal response of the line to the illuminating flux. Previous studies concentrated on the calculation of reverberation signatures from static sources illuminating the disc. In this paper we focus on the more physically justified case of flares located above the accretion disc and corotating with it. We compute the time-dependent iron line, taking into account all general relativistic effects, and show that its shape is of a very complex nature, and we also present light curves accompanying the iron line variability. We suggest that present and future X-ray satellites like XMM or Constellation-X may be capable of detecting features present in the computed reverberation maps.     

The anticorrelation between the hard X-ray photon index and the Eddington ratio in low-luminosity active galactic nuclei

We find a significant anticorrelation between the hard X-ray photon index Γ and the Eddington ratio   L _bol/ L _Edd  for a sample of low-ionization nuclear emission-line regions and local Seyfert galaxies, compiled from literature with Chandra or XMM–Newton observations. This result is in contrast with the positive correlation found in luminous active galactic nuclei (AGN), while it is similar to that of X-ray binaries (XRBs) in the low/hard state. Our result is qualitatively consistent with the spectra produced from advection-dominated accretion flows (ADAFs). It implies that the X-ray emission of low-luminosity active galactic nuclei (LLAGN) may originate from the Comptonization process in ADAF, and the accretion process in LLAGN may be similar to that of XRBs in the low/hard state, which is different from that in luminous AGN.     

Broad reprocessed Balmer emission from warped accretion discs

Double-peaked broad emission lines in active galactic nuclei are generally considered to be formed in an accretion disc. In this paper, we compute the profiles of reprocessing emission lines from a relativistic, warped accretion disc around a black hole in order to explore the possibility that certain asymmetries in the double-peaked emission-line profile which cannot be explained by a circular Keplerian disc may be induced by disc warping. The disc warping also provides a solution for the energy budget in the emission-line region because it increases the solid angle of the outer disc portion subtended to the inner portion of the disc. We adopted a parametrized disc geometry and a central point-like source of ionizing radiation to capture the main characteristics of the emission-line profile from such discs. We find that the ratio between the blue and red peaks of the line profiles becoming less than unity can be naturally predicted by a twisted warped disc, and a third peak can be produced in some cases. We show that disc warping can reproduce the main features of multipeaked line profiles of four active galactic nuclei from the Sloan Digital Sky Survey.     

Testing a model of variability of X-ray reprocessing features in active galactic nuclei

A number of recent results from X-ray observations of active galactic nuclei involving the Fe K α line (reduction of line variability compared with the X-ray continuum variability, the X-ray 'Baldwin effect') were attributed to the presence of a hot, ionized skin of an accretion disc, suppressing emission of the line. The ionized skin appears as a result of the thermal instability of X-ray irradiated plasma. We test this hypothesis by computing the Thomson thickness of the hot skin on top of the αP _tot Shakura–Sunyaev disc, by simultaneously solving the vertical structure of both the hot skin and the disc. We then compute a number of relations between observable quantities, e.g. the hard X-ray flux, amplitude of the observed reprocessed component, relativistic smearing of the K α line and rms variability of the hard X-rays. These relations can be compared with present and future observations. We point out that this mechanism is unlikely to explain the behaviour of the X-ray source in MCG–6-30-15, where there are a number of arguments against the existence of a thick hot skin, but it can work for some other Seyfert 1 galaxies.     

X-ray spectra produced by a hot plasma containing cold clouds

We compute the hard X-ray spectra from a hot plasma pervaded by small cold dense clouds. The main cooling mechanism of the plasma is Compton cooling by the soft thermal emission from the clouds. We compute numerically the equilibrium temperature of the plasma together with the escaping spectrum. The spectrum depends mainly on the amount of cold clouds filling the hot phase. The clouds covering factor is constrained to be low in order to produce spectra similar to those observed in Seyfert galaxies and X-ray binaries, implying that an external reflector is required in order to reproduce the full range of observed reflection amplitudes. We also derive analytical estimates for the X-ray spectral slope and reflection amplitude using an escape probability formalism.     

Luminous hot accretion discs

We find a new two-temperature hot branch of equilibrium solutions for stationary accretion discs around black holes. In units of Eddington accretion rate defined as 10 L _Edd c ², the accretion rates to which these solutions correspond are within the range ̇ ₁≲ ̇ ≲1, where ̇ ₁ is the critical rate of advection-dominated accretion flow (ADAF). In these solutions, the energy loss rate of the ions by Coulomb energy transfer between the ions and electrons is larger than the viscously heating rate and it is the advective heating together with the viscous dissipation that balances the Coulomb cooling of ions. When ̇ ₁≲ ̇ ≲ ̇ ₂, where ̇ ₂∼5 ̇ ₁<1, the accretion flow remains hot throughout the disc. When ̇ ₂≲ ̇ ≲1, Coulomb interaction will cool the inner region of the disc within a certain radius ( r _tr∼several tens of Schwarzschild radii or larger depending on the accretion rate and the outer boundary condition) and the disc will collapse on to the equatorial plane and form an optically thick cold annulus. Compared with ADAF, these hot solutions are much more luminous because of the high accretion rate and efficiency; therefore, we call them luminous hot accretion discs.     

Three-dimensional smoothed particle hydrodynamics simulations of radiation-driven warped accretion discs

We present three-dimensional smoothed particle hydrodynamics calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc, resulting in a torque that acts on the disc to induce a twist or warp. Initially planar discs are unstable to warping driven by the radiation torque and, in general, the warps also precess in a retrograde direction relative to the orbital flow. We simulate a number of X-ray binary systems which have different mass ratios, using a number of different luminosities for each. Radiation-driven warping occurs for all systems simulated. For mass ratios   q ∼ 0.1  a moderate warp occurs in the inner disc while the outer disc remains in the orbital plane (cf. X 1916−053). For less extreme mass ratios, the entire disc tilts out of the orbital plane (cf. Her X–1). For discs that are tilted out of the orbital plane in which the outer edge material of the disc is precessing in a prograde direction, we obtain both positive and negative superhumps simultaneously in the dissipation light curve (cf. V603 Aql).     

How the X-ray spectrum of a narrow-line Seyfert 1 galaxy may be reflection-dominated

A model for the inner regions of accretion flows is presented where, owing to disc instabilities, cold and dense material is clumped into deep sheets or rings. Surrounding these density enhancements is hot, tenuous gas where coronal dissipation processes occur. We expect this situation to be most relevant when the accretion rate is close to Eddington and the disc is radiation-pressure-dominated, and so may apply to narrow-line Seyfert 1 (NLS1) galaxies. In this scenario, the hard X-ray source is obscured for most observers, and so the detected X-ray emission would be dominated by reflection off the walls of the sheets. A simple Comptonization calculation shows that the large photon-indices characteristic of NLS1s would be a natural outcome of two reprocessors closely surrounding the hard X-ray source. We test this model by fitting the XMM-Newton spectrum of the NLS1 1H  0707–495  between 0.5 and 11 keV with reflection-dominated ionized disc models. A very good fit is found with three different reflectors each subject to the same  Γ=2.35  power law. An iron overabundance is still required to fit the sharp drop in the spectrum at around 7 keV. We note that even a small corrugation of the accretion disc may result in  Γ>2  and a strong reflection component in the observed spectrum. Therefore, this model may also explain the strength and the variability characteristics of the MCG–6-30-15 Fe K α line. The idea needs to be tested with further broad-band XMM-Newton observations of NLS1s.     

On the complex disc–corona interactions in the soft spectral states of soft X-ray transients

Accreting black holes show a complex and diverse behaviour in their soft spectral states. Although these spectra are dominated by a soft, thermal component which almost certainly arises from an accretion disc, there is also a hard X-ray tail indicating that some fraction of the accretion power is instead dissipated in hot, optically thin coronal material. During such states, best observed in the early outburst of soft X-ray transients, the ratio of power dissipated in the hot corona to that in the disc can vary from ∼ 0 (pure disc accretion) to ∼ 1 (equal power in each). Here we present results of spectral analyses of a number of sources, demonstrating the presence of complex features in their energy spectra. Our main findings are: (1) the soft components are not properly described by a thermal emission from accretion discs: they are appreciably broader than can be described by disc blackbody models even including relativistic effects, and (2) the spectral features near     commonly seen in such spectra can be well described by reprocessing of hard X-rays by optically thick, highly ionized, relativistically moving plasma.     

Broad Iron Lines in AGN and X-Ray Binaries

Several AGN and black hole X-ray binaries show a clear very broad iron line, which is strong evidence that the black holes are rapidly spinning. Detailed analysis of these objects shows that the emission line is not significantly affected by absorption and that the source variability is principally due to variation in amplitude of a power-law. Underlying this is a much less variable, relativistically-smeared, reflection-dominated, component which carries the imprint of strong gravity at a few gravitational radii. The strong gravitational light bending in these regions then explains the power-law variability as due to changes in height of the primary X-ray source above the disc. The reflection component, in particular its variability and the profile of the iron line, enables us to study the innermost regions around an accreting, spinning, black hole.     

Long-term optical and X-ray variability of the Seyfert galaxy Markarian 79

We present the results of concurrent X-ray and optical monitoring of the Seyfert 1 galaxy Mrk 79 over a period of more than 5 yr. We find that on short to medium time-scales (days to a few tens of days) the 2–10 keV X-ray and optical u - and V -band fluxes are significantly correlated, with a delay between the bands consistent with 0 d. We show that most of these variations may be well reproduced by a model where the short-term optical variations originate from reprocessing of X-rays by an optically thick accretion disc. The optical light curves, however, also display long time-scale variations over thousands of days, which are not present in the X-ray light curve. These optical variations must originate from an independent variability mechanism and we show that they can be produced by variations in the (geometrically) thin disc accretion rate as well as by varying reprocessed fractions through changes in the location of the X-ray corona.