Constraining the number of compact remnants near Sgr A*

Due to dynamical friction stellar mass black holes and neutron stars are expected to form high-density cusps in the inner parsec of our Galaxy. These compact remnants, expected to number around 20 000, may be accreting cold dense gas present there, and give rise to potentially observable X-ray emission. Here we build a simple but detailed time-dependent model of such emission. The possibility that these accretion flows are radiatively inefficient is taken into account and brings in some uncertainty in the conclusions. Despite this uncertainty, we find that at least several X-ray sources of this nature should be detectable with Chandra at any one time. Turning this issue around, we also ask a question of what current observational constraints might be telling us about the total number of compact remnants. In our 'best guess' model, a cusp of ∼40 000 remnants overpredicts the number of discrete sources and the total X-ray luminosity of the inner parsec, and is hence ruled out. In the most radiatively inefficient scenario that we consider, the radiative efficiency is set to be as small as  ɛ= 10⁻⁵  . In this rather unlikely scenario, a cusp of ∼40 000 black holes would be allowed by the data, but several individual sources should still be visible. Future observations of the distribution and orbits of the cold ionized gas in the inner parsec of our Galaxy will put tighter constraints on the cusp of compact remnants.     

Possible evidence for the disc origin for the powering of jets in Sgr A* and nearby elliptical galaxies

Recent VLBA observation indicates the existence of an elongated (jet) structure in the compact radio source Sgr A*. This is hard to explain in the context of advection-dominated accretion flow (ADAF) model for this source. On the other hand, the mass accretion rate favoured by ADAF is 10–20 times smaller than that favoured by the hydrodynamical simulation based on Bondi capture. If the latter were adopted, the predicted radio flux would significantly exceed the observation. A similar situation exists in the case of nearby giant ellipticals, where the canonical ADAF model – the widely assumed standard model for these sources – also significantly overpredicts the radio flux. Based on these facts, in this paper we propose a truncated ADAF model for Sgr A* and three ellipticals M87, NGC 4649 and NGC 4636. We assume that the accretion disc is truncated at a certain radius R _tr within which the jet forms by extracting the energy of the disc. The radio flux is greatly suppressed owing to the radiative truncation of the disc and the fits to the observational data are excellent. For example, for Sgr A*, the model fits the observational spectrum very well from radio including the 'excess' below the break frequency to hard X-ray under a high accretion rate near the simulation value, and the predicted size-frequency relationship is also in excellent agreement with the observation; for M87, the predicted upper limit of the jet location is 24 R _g, in excellent agreement with the observational result that the jet is formed on scales smaller than 30 R _g, and the ≈20 per cent variability at ∼1 keV – which is hard to explain in another model that succeeded in explaining the low radio flux of M87 – is also marginally interpreted. The success of the model supplies possible evidence for the disc rather than the hole origin for the powering of jets.     

The circumnuclear material in the Galactic Centre: a clue to the accretion process

On the basis of 'sticky particle' calculations, it is argued that the gas features observed within 10 pc of the Galactic Centre — the circumnuclear disc (CND) and the ionized gas filaments — as well as the newly formed stars in the inner 1 pc can be understood in terms of tidal capture and disruption of gas clouds on low angular momentum orbits in a potential containing a point mass. The calculations demonstrate that a dissipative component forms a 'dispersion ring', an asymmetric elliptical torus precessing counter to the direction of rotation, and that this shape can be maintained for many orbital periods. For a range of plausible initial conditions, such a structure can explain the morphology and kinematics of the CND and of the most conspicuous ionized filament. While forming the dispersion ring, a small cloud with low specific angular momentum is drawn into a long filament which repeatedly collides with itself at high velocity. The compression in strong shocks is likely to lead to star formation even in the near tidal field of the point mass. This process may have general relevance to accretion on to massive black holes in normal and active galactic nuclei.     

Variable accretion and emission from the stellar winds in the Galactic Centre

We present numerical simulations of stellar wind dynamics in the central parsec of the Galactic Centre, studying in particular the accretion of gas on to Sgr A*, the supermassive black hole. Unlike our previous work, here we use state-of-the-art observational data on orbits and wind properties of individual wind-producing stars. Since wind velocities were revised upwards and non-zero eccentricities were considered, our new simulations show fewer clumps of cold gas and no conspicuous disc-like structure. The accretion rate is dominated by a few close 'slow-wind stars' ( v _w≤ 750 km s⁻¹), and is consistent with the Bondi estimate, but variable on time-scales of tens to hundreds of years. This variability is due to the stochastic infall of cold clumps of gas, as in earlier simulations, and to the eccentric orbits of stars. The present models fail to explain the high luminosity of Sgr A* a few hundred years ago implied by Integral observations, but we argue that the accretion of a cold clump with a small impact parameter could have caused it. Finally, we show the possibility of constraining the total mass-loss rate of the 'slow-wind stars' using near infrared observations of gas in the central few arcseconds.     

X-rays from cusps of compact remnants near galactic centres

Compact remnants – stellar mass black holes and neutron stars formed in the inner few parsec of galactic centres are predicted to sink into the central parsec due to dynamical friction on low-mass stars, forming a high concentration cusp. Same physical region may also contain very high-density molecular clouds and accretion discs that are needed to fuel supermassive black hole (SMBH) activity. Here we estimate gas capture rates on to the cusp of stellar remnants, and the resulting X-ray luminosity, as a function of the accretion disc mass. At low disc masses, most compact objects are too dim to be observable, whereas in the high disc case most of them are accreting at their Eddington rates. We find that for low accretion disc masses, compact remnant cusps may be more luminous than the central SMBHs. This 'diffuse' emission may be of importance for local moderately bright active galactic nuclei (AGNs), especially low-luminosity AGNs. We also briefly discuss how this expected emission can be used to put constraints on the black hole cusp near our Galactic Centre.     

Component thermalization time-scale estimate for the advection dominated accretion flow around Sgr A

We report here on a calculation of thermalization time-scale of the two temperature advection dominated accretion flow (ADAF) model. It is established that time required to equalize the electron and ion temperatures via electron-ion collisions in the ADAF with plausible physical parameters greatly exceeds the accretion time, which corroborates validity one of the crucial assumptions of the ADAF model, namely the existence of a hot two temperature plasma. This work is motivated by the recent success of the ADAF model [Nature 394 (1998) 651; MNRAS 304 (1999) 501] in explaining the emitted spectrum of Sgr A*, and it is complementary to the similar analysis of Mahadevan and Quataert [ApJ 490 (1997) 605].     

Variability associated with α accretion disc theory for standard and advection-dominated discs

The α turbulent viscosity formalism for accretion discs must be interpreted as a mean field theory, modelling a steady state only on spatial or time-scales greater than those of the turbulence. The extent of the scale separation determines the relative precision error (RPE) of the predicted luminosity L _ν. Turbulence and the use of α implies that (1) field line stretching gives a magnetic pressure  α²/6 of the total pressure generally, and a one-to-one relation between α and the pressure ratio for thin discs, and (2) large turbulent scales in advection-dominated accretion flows (ADAFs) predict a lower L _ν precision than thin discs for a given observation duration and central mass. The allowed variability (or RPE) at frequency ν increases with the size of the contributing region. For X-ray binary ADAFs, the RPE ∼ 5 per cent at R  ≤ 1000 Schwarzchild radii ( R _s) for averages over  1000 s. However, current data for galaxies like NGC 4258 and M87 give RPEs in L _ν of 50–100 per cent even at R  ≤ 100  R _S. More data are required, but systematic deviations from ADAF predictions are more significant than random deviations, and may constrain properties of the turbulence, the accretion mode, the assumption of a steady state or the accretion rate.     

The effect of radiation stress on gas motions in AGN accretion discs

In this paper we look at one of the effects of irradiation on a warped accretion disc in the context of active galactic nuclei (AGN). A warp will catch a substantial amount of the radiation emitted by the central object. We consider the fluid motions that may arise inside a warped disc when the surface is subject to a radiation stress, and also the net mass flows that result. We find that, to first order, we have a balance of the viscous and Coriolis-type forces. The radial radiation stress causes outward motion of the surface layer, but only the azimuthal Poynting–Robertson drag leads to an increase in the net accretion rate. We investigate the distribution of the velocity perturbations and find them to be significant in determining the local structure of the disc.
An unexpected result is that the picture changes significantly when we take into account the periodic illumination of the warped disc. A type of resonance at the local Keplerian rotation frequency causes a flow that penetrates the whole thickness of the disc; these flows are faster than the flows due to unchanging illumination. They totally dominate the induced flows in terms of sheer mass, but significant impact on disc structure still occurs only near the surface, where velocity perturbations typically go up to some kilometres per second.     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

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.     

The fraction of galaxies that contain active nuclei and their accretion rates

We investigate the relationship between the present-day optical luminosity function of galaxies and the X-ray luminosity function of Seyfert 1s to determine the fraction of galaxies that host Seyfert 1 nuclei and their Eddington ratios. The local type 1 active galactic nuclei (AGN) X-ray luminosity function is well reproduced if ∼1 per cent of all galaxies are type 1 Seyferts which have Eddington ratios of ∼10⁻³. However, in such a model the X-ray luminosity function is completely dominated by AGN in E and S0 galaxies, contrary to the observed mix of Seyfert host galaxies. To obtain a plausible mix of AGN host galaxy morphologies requires that the most massive black holes in E and S0 galaxies accrete with lower Eddington ratios, or have a lower incidence of Seyfert activity, than the central black holes of later-type galaxies.