Spherical Accretion in Nearby Weakly Active Galaxies

We consider the sample of weakly active galaxies situated in 'Local Universe' collected in the paper of Pellegrini (2005) with inferred accretion efficiencies from $10^{-2}$ to $10^{-7}$. We apply a model of spherically symmetrical Bondi accretion for given parameters ($M_{BH}$,$T_{\infty}$,$\rho_{\infty}$,) taken from observation. We calculate spectra emitted by the gas accreting onto its central objects using Monte Carlo method including synchrotron and bremsstrahlung photons as seed photons. We compare our results with observed nuclear X-ray luminosities $L_{X,nuc}$ (0.3-10 keV) of the sample. Model is also tested for different external medium parameters ($\rho_{\infty}$ and $T_{\infty}$) and different free parameters of the model. Our model is able to explain most of the observed nuclear luminosities $L_X$ under an assumption that half of the compresion energy is transfered directly to the electrons.Comment: 19 pages, 11 figures, accepted to A&

GMOS Integral Field Spectroscopy of a Merging System with Enhanced Balmer Absorption

In this paper we present the three dimensional dynamics of the galaxy SDSS J101345.39+011613.66, selected for its unusually strong Balmer absorption lines (Wo(H-delta)=7.5A). Using the GMOS-South IFU in Nod & Shuffle mode we have mapped the continuum and optical absorption lines of this z=0.1055 field galaxy. This galaxy has a disturbed morphology, with a halo of diffuse material distributed asymmetrically toward the north. Using the [OII] emission line (Wo([OII])=4.1A) we find that the gas and hot OB stars are offset from the older stars in the system. The gas also has a spatially extended and elongated morphology with a velocity gradient of 100+/-20km/s across 6kpc in projection. Using the strong H-gamma and H-delta absorption lines we find that the A- stars are widely distributed across the system and are not centrally concentrated arguing that the A-star population has formed in molecular clouds outside the nucleus. By cross correlating the spectra from the datacube with an A-star template we find evidence that the A-star population has a 40km/s shear in the same direction as the gas. The disturbed morphology, strong colour gradients and strong H-delta and H-gamma absorption lines in SDSS J101345.39 argue that this is a recent tidal interaction/merger between a passive elliptical and star-forming galaxy. Although based on a single object, these results show that we can spatially resolve and constrain the dynamics of this short lived (yet important) phase of galaxy formation in which the evolutionary process take galaxies from star-forming to their quiescent end products.Comment: 7 pages, 7 figures. Accepted for publication in Ap

Time Dependent Models of Flares from Sagittarius A*

The emission from Sgr A*, the supermassive black hole in the Galactic Center, shows order of magnitude variability ("flares") a few times a day that is particularly prominent in the near-infrared (NIR) and X-rays. We present a time-dependent model for these flares motivated by the hypothesis that dissipation of magnetic energy powers the flares. We show that episodic magnetic reconnection can occur near the last stable circular orbit in time-dependent magnetohydrodynamic simulations of black hole accretion - the timescales and energetics of these events are broadly consistent with the flares from Sgr A*. Motivated by these results, we present a spatially one-zone time-dependent model for the electron distribution function in flares, including energy loss due to synchrotron cooling and adiabatic expansion. Synchrotron emission from transiently accelerated particles can explain the NIR/X-ray lightcurves and spectra of a luminous flare observed 4 April 2007. A significant decrease in the magnetic field strength during the flare (coincident with the electron acceleration) is required to explain the simultaneity and symmetry of the simultaneous lightcurves. Our models predict that the NIR and X-ray spectral indices differ by 0.5 and that there is only modest variation in the spectral index during flares. We also explore implications of this model for longer wavelength (radio-submm) emission seemingly associated with X-ray and NIR flares; we argue that a few hour decrease in the submm emission is a more generic consequence of large-scale magnetic reconnection than delayed radio emission from adiabatic expansion.Comment: 18 pages, 10 figures, ApJ accepte

Dark matter haloes determine the masses of supermassive black holes

The energy and momentum deposited by the radiation from accretion onto the supermassive black holes (BHs) that reside at the centres of virtually all galaxies can halt or even reverse gas inflow, providing a natural mechanism for supermassive BHs to regulate their growth and to couple their properties to those of their host galaxies. However, it remains unclear whether this self-regulation occurs on the scale at which the BH is gravitationally dominant, on that of the stellar bulge, the galaxy, or that of the entire dark matter halo. To answer this question, we use self-consistent simulations of the co-evolution of the BH and galaxy populations that reproduce the observed correlations between the masses of the BHs and the properties of their host galaxies. We first confirm unambiguously that the BHs regulate their growth: the amount of energy that the BHs inject into their surroundings remains unchanged when the fraction of the accreted rest mass energy that is injected, is varied by four orders of magnitude. The BHs simply adjust their masses so as to inject the same amount of energy. We then use simulations with artificially reduced star formation rates to demonstrate explicitly that BH mass is not set by the stellar mass. Instead, we find that it is determined by the mass of the dark matter halo with a secondary dependence on the halo concentration, of the form that would be expected if the halo binding energy were the fundamental property that controls the mass of the BH. We predict that the logarithmic slope of the relation between dark matter halo mass and black hole mass is 1.55+/-0.05 and that the scatter around the mean relation in part reflects the scatter in the halo concentration-mass relation.Comment: MNRAS accepted. 6 pages, 3 figures. v2: Minor changes in response to referee comment

Highly ionized gas on galaxy scales: mapping the interacting Seyfert galaxy LEDA 135736

We have used the VIMOS IFU to map the properties of the Seyfert 1.9 galaxy LEDA 135736. These maps reveal a number of interesting features including: an Extended Narrow Line Region detectable out to 9 kpc, an area of intense star formation located at a projected distance of 12 kpc from the centre, an elliptical companion galaxy, and kinematic features, aligned along the long-axis of the ENLR, that are consistent with radio jet-driven mass outflow. We propose that the ENLR results from extra-planar gas ionized by the AGN, and that the AGN in turn might be triggered by interaction with the companion galaxy, which can also explain the burst of star formation and morphological features. Only about two percent of the ENLR's kinetic energy is in the mass outflow. We infer from this that the bulk of mechanical energy imparted by the jet is used to heat this gas.Comment: Accepted for publication in MNRAS Letter

Direct cosmological simulations of the growth of black holes and galaxies

We investigate the coupled formation and evolution of galaxies and their embedded supermassive black holes using state-of-the-art hydrodynamic simulations of cosmological structure formation. For the first time, we self-consistently follow the dark matter dynamics, radiative gas cooling, star formation, as well as black hole growth and associated feedback processes, starting directly from initial conditions appropriate for the LambdaCDM cosmology. Our modeling of the black hole physics is based on an approach we have developed in simulations of isolated galaxy mergers. Here we examine: (i) the predicted global history of black hole mass assembly (ii) the evolution of the local black hole-host mass correlations and (iii) the conditions that allow rapid growth of the first quasars, and the properties of their hosts and descendants today. We find a total black hole mass density in good agreement with observational estimates. The black hole accretion rate density peaks at lower redshift and evolves more strongly at high redshift than the star formation rate density, but the ratio of black hole to stellar mass densities shows only a moderate evolution at low redshifts. We find strong correlations between black hole masses and properties of the stellar systems, agreeing well with the measured local M_BH-sigma and M_BH -M_* relationships, but also suggesting (dependent on the mass range) a weak evolution with redshift in the normalization and the slope. Our simulations also produce massive black holes at high redshift, due to extended periods of exponential growth in regions that collapse early and exhibit strong gas inflows. These first supermassive BH systems however are not necessarily the most massive ones today, since they are often overtaken in growth by quasars that form later. (abridged)Comment: 22 pages, 17 figures, submitted to Ap

The Accretion Disc Particle Method for Simulations of Black Hole Feeding and Feedback

Black holes grow by accreting matter from their surroundings. However, angular momentum provides an efficient natural barrier to accretion and so only the lowest angular momentum material will be available to feed the black holes. The standard sub-grid model for black hole accretion in galaxy formation simulations - based on the Bondi-Hoyle method - does not account for the angular momentum of accreting material, and so it is unclear how representative the black hole accretion rate estimated in this way is likely to be. In this paper we introduce a new sub-grid model for black hole accretion that naturally accounts for the angular momentum of accreting material. Both the black hole and its accretion disc are modelled as a composite accretion disc particle. Gas particles are captured by the accretion disc particle if and only if their orbits bring them within its accretion radius R_acc, at which point their mass is added to the accretion disc and feeds the black hole on a viscous timescale t_visc. The resulting black hole accretion rate (dM/dt)_BH powers the accretion luminosity L_acc ~ (dM/dt)_BH, which drives black hole feedback. Using a series of controlled numerical experiments, we demonstrate that our new accretion disc particle method is more physically self-consistent than the Bondi-Hoyle method. We also discuss the physical implications of the accretion disc particle method for systems with a high degree of rotational support, and we argue that the M_BH-sigma relation in these systems should be offset from the relation for classical bulges and ellipticals, as appears to be observed.Comment: Accepted for publication in MNRAS; 9 pages, 5 figure

Measuring the kinetic power of AGN in the radio mode

(Abridged) We have studied the relationship among nuclear radio and X-ray power, Bondi rate and the kinetic luminosity of sub-Eddington active galactic nuclear (AGN) jets. Besides the recently discovered correlation between jet kinetic and Bondi power, we show that a clear correlation exists also between Eddington-scaled kinetic power and bolometric luminosity, given by: Log(L_kin/L_Edd)=0.49*Log(L_bol/L_Edd)-0.78. The measured slope suggests that these objects are in a radiatively inefficient accretion mode, and has been used to put stringent constraints on the properties of the accretion flow. We found no statistically significant correlations between Bondi power and bolometric AGN luminosity, apart from that induced by their common dependence on L_kin. Analyzing the relation between kinetic power and radio core luminosity, we are then able to determine, statistically, both the probability distribution of the mean jets Lorentz factor, peaking at \Gamma~7, and the intrinsic relation between kinetic and radio core luminosity, that we estimate as: Log(L_kin)=0.81*Log(L_R)+11.9, in good agreement with theoretical predictions of synchrotron jet models. With the aid of these findings, quantitative assessments of kinetic feedback from supermassive black holes in the radio mode will be possible based on accurate determinations of the central engine properties alone. As an example, Sgr A* may follow the correlations of radio mode AGN, based on its observed radiative output and on estimates of the accretion rate both at the Bondi radius and in the inner flow. If this is the case, the SMBH in the Galactic center is the source of ~ 5 times 10^38 ergs/s of mechanical power, equivalent to about 1.5 supernovae every 10^5 years.Comment: 13 pages, 6 figures. Accepted for publication in MNRA

Spectroastrometry of rotating gas disks for the detection of supermassive black holes in galactic nuclei. II. Application to the galaxy Centaurus A (NGC 5128)

We measure the black hole mass in the nearby active galaxy Centaurus A (NGC 5128) using a new method based on spectroastrometry of a rotating gas disk. The spectroastrometric approach consists in measuring the photocenter position of emission lines for different velocity channels. In a previous paper we focused on the basic methodology and the advantages of the spectroastrometric approach with a detailed set of simulations demonstrating the possibilities for black hole mass measurements going below the conventional spatial resolution. In this paper we apply the spectroastrometric method to multiple longslit and integral field near infrared spectroscopic observations of Centaurus A. We find that the application of the spectroastrometric method provides results perfectly consistent with the more complex classical method based on rotation curves: the measured BH mass is nearly independent of the observational setup and spatial resolution and the spectroastrometric method allows the gas dynamics to be probed down to spatial scales of ~0.02", i.e. 1/10 of the spatial resolution and ~1/50 of BH sphere of influence radius. The best estimate for the BH mass based on kinematics of the ionized gas is then log(MBH (sin i)^2/M\odot)=7.5 \pm 0.1 which corresponds to MBH = 9.6(+2.5-1.8) \times 10^7 M\odot for an assumed disk inclination of i = 35deg. The complementarity of this method with the classic rotation curve method will allow us to put constraints on the disk inclination which cannot be otherwise derived from spectroastrometry. With the application to Centaurus A, we have shown that spectroastrometry opens up the possibility of probing spatial scales smaller than the spatial resolution, extending the measured MBH range to new domains which are currently not accessible: smaller BHs in the local universe and similar BHs in more distant galaxies