The coalescence of a binary black hole can be accompanied by a large gravitational recoil due to anisotropic emission of gravitational waves. A recoiling supermassive black hole (SBH) can subsequently undergo long-lived oscillations in the potential
well of its host galaxy, suggesting that offset SBHs may be common in the cores of massive ellipticals. We have analyzed HST archival images of 14 nearby core ellipticals, finding evidence for small ($lesssim 10$ pc) displacements between the AGN (locating the SBH) and the center of the galaxy (the mean photocenter) in 10 of them. Excluding objects that may be affected by large-scale isophotal asymmetries, we consider six galaxies to have detected displacements, including M87, where a displacement was previously reported by Batcheldor et al. 2010. In individual objects, these displacements can be attributed to residual gravitational recoil oscillations following a major or minor merger within the last few Gyr. For plausible merger rates, however, there is a high probability of larger displacements than those observed, if SBH coalescence took place in these galaxies. Remarkably, the AGN-photocenter displacements are approximately aligned with the radio source axis in four of the six galaxies with displacements, including three of the four having relatively powerful kpc-scale jets. This suggests intrinsic asymmetries in radio jet power as a possible displacement mechanism, although approximate alignments are also expected for gravitational recoil. Orbital motion in SBH binaries and interactions with massive perturbers can produce the observed displacement amplitudes but do not offer a ready explanation for the alignments.
We present new CRIRES spectroscopic observations of BrGamma in the nuclear region of the Circinus galaxy, obtained with the aim of measuring the black hole (BH) mass with the spectroastrometric technique. The Circinus galaxy is an ideal benchmark for
the spectroastrometric technique given its proximity and secure BH measurement obtained with the observation of its nuclear H2O maser disk. The kinematical data have been analyzed both with the classical method based on the analysis of the rotation curves and with the new method developed by us and based on spectroastrometry. The classical method indicates that the gas disk rotates in the gravitational potential of an extended stellar mass distribution and a spatially unresolved mass of (1.7 +- 0.2) 10^7 Msun, concentrated within r < 7 pc. The new method is capable of probing gas rotation at scales which are a factor ~3.5 smaller than those probed by the rotation curve analysis. The dynamical mass spatially unresolved with the spectroastrometric method is a factor ~2 smaller, 7.9 (+1.4 -1.1) 10^6 Msun indicating that spectroastrometry has been able to spatially resolve the nuclear mass distribution down to 2 pc scales. This unresolved mass is still a factor ~4.5 larger than the BH mass measurement obtained with the H2O maser emission indicating that it has not been possible to resolve the sphere of influence of the BH. Based on literature data, this spatially unresolved dynamical mass distribution is likely dominated by molecular gas and it has been tentatively identified with the circum-nuclear torus which prevents a direct view of the central BH in Circinus. This mass distribution, with a size of ~2pc, is similar in shape to that of the star cluster of the Milky Way suggesting that a molecular torus, forming stars at a high rate, might be the earlier evolutionary stage of the nuclear star clusters which are common in late type spirals.
Most galaxy evolutionary models require quasar feedback to regulate star formation in their host galaxies. In particular, at high redshift, models expect that feedback associated with quasar-driven outflows is so efficient that the gas in the host ga
laxy is largely swept away or heated up, hence suppressing star formation in massive galaxies. We observationally investigate this phenomenon by using VLT-SINFONI integral field spectroscopy of the luminous quasar 2QZJ002830.4-281706 at z=2.4. The spectra sample the optical emission lines redshifted into the near-IR. The [OIII]5007 emission-line kinematics map reveals a massive outflow on scales of several kpc. The detection of narrow Halpha emission reveals star formation in the quasar host galaxy, with SFR=100 Msun/yr. However, the star formation is not distributed uniformly, but is strongly suppressed in the region with the highest outflow velocity and highest velocity dispersion. This result indicates that star formation in this region is strongly quenched by the quasar outflow, which is cleaning the galaxy disk of its molecular gas. This is one of the first direct observational proofs of quasar feedback quenching the star formation at high redshift.
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 f
or 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/Modot)=7.5 pm 0.1 which corresponds to MBH = 9.6(+2.5-1.8) times 10^7 Modot 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.
Galaxy dynamical masses are important physical quantities to constrain galaxy evolutionary models, especially at high redshifts. However, at z~2 the limited signal to noise ratio and spatial resolution of the data usually do not allow spatially resol
ved kinematical modeling and very often only virial masses can be estimated from line widths. But even such estimates require a good knowledge of galaxy size, which may be smaller than the spatial resolution. Spectroastrometry is a technique which combines spatial and spectral resolution to probe spatial scales significantly smaller than the spatial resolution of the observations. Here we apply it to the case of high-z galaxies and present a method based on spectroastrometry to estimate dynamical masses of high z galaxies, which overcomes the problem of size determination with poor spatial resolution. We construct and calibrate a spectroastrometric virial mass estimator, modifying the classical virial mass formula. We apply our method to the [O III] or H{alpha} emission line detected in z~2-3 galaxies from AMAZE, LSD and SINS samples and we compare the spectroastrometric estimator with dynamical mass values resulting from full spatially resolved kinematical modeling. The spectroastrometric estimator is found to be a good approximation of dynamical masses, presenting a linear relation with a residual dispersion of only 0.15 dex. This is a big improvement compared to the classical virial mass estimator which has a non linear relation and much larger dispersion (0.47 dex) compared to dynamical masses. By applying our calibrated estimator to 16 galaxies from the AMAZE and LSD samples, we obtain masses in the ~10^7-10^10 Modot range extending the mass range attainable with dynamical modeling.
We present a SINFONI integral field kinematical study of 33 galaxies at z~3 from the AMAZE and LSD projects which are aimed at studying metallicity and dynamics of high-redshift galaxies. The number of galaxies analyzed in this paper constitutes a si
gnificant improvement compared to existing data in the literature and this is the first time that a dynamical analysis is obtained for a relatively large sample of galaxies at z~3. 11 galaxies show ordered rotational motions (~30% of the sample), in these cases we estimate dynamical masses by modeling the gas kinematics with rotating disks and exponential mass distributions. We find dynamical masses in the range 2 times 10^9 Modot - 2 times 10^11 Modot with a mean value of ~ 2 times 10^10 Modot. By comparing observed gas velocity dispersion with that expected from models, we find that most rotating objects are dynamically hot, with intrinsic velocity dispersions of the order of ~90 km s-1. The median value of the ratio between the maximum disk rotational velocity and the intrinsic velocity dispersion for the rotating objects is 1.6, much lower than observed in local galaxies value (~10) and slightly lower than the z~2 value (2 - 4). Finally we use the maximum rotational velocity from our modeling to build a baryonic Tully-Fisher relation at z~3. Our measurements indicate that z~3 galaxies have lower stellar masses (by a factor of ten on average) compared to local galaxies with the same dynamical mass. However, the large observed scatter suggests that the Tully-Fisher relation is not yet in place at these early cosmic ages, possibly due to the young age of galaxies. A smaller dispersion of the Tuly-Fisher relation is obtained by taking into account the velocity dispersion with the use of the S_0.5 indicator, suggesting that turbulent motions might have an important dynamical role.
This is the first in a series of papers in which we study the application of spectroastrometry in the context of gas kinematical studies aimed at measuring the mass of supermassive black holes. The spectroastrometrical method consists in measuring th
e photocenter of light emission in different wavelength or velocity channels. In particular we explore the potential of spectroastrometry of gas emission lines in galaxy nuclei to constrain the kinematics of rotating gas disks and to measure the mass of putative supermassive black holes. By means of detailed simulations and test cases, we show that the fundamental advantage of spectroastrometry is that it can provide information on the gravitational potential of a galaxy on scales significantly smaller (~ 1/10) than the limit imposed by the spatial resolution of the observations. We then describe a simple method to infer detailed kinematical informations from spectroastrometry in longslit spectra and to measure the mass of nuclear mass concentrations. Such method can be applied straightforwardly to integral field spectra, which do not have the complexities due to a partial spatial covering of the source in the case of longslit spectra.
The application of the virial theorem to the Broad Line Region of Active Galactic Nuclei allows Black Hole mass estimates for large samples of objects at all redshifts. In a recent paper we showed that ionizing radiation pressure onto BLR clouds affe
cts virial BH mass estimates and we provided empirically calibrated corrections. More recently, a new test of the importance of radiation forces has been proposed: the MBH-sigma relation has been used to estimate MBH for a sample of type-2 AGN and virial relations (with and without radiation pressure) for a sample of type-1 AGN extracted from the same parent population. The observed L/LEdd distribution based on virial BH masses is in good agreement with that based on MBH-sigma only if radiation pressure effects are negligible, otherwise significant discrepancies are observed. In this paper we investigate the effects of intrinsic dispersions associated to the virial relations providing MBH, and we show that they explain the discrepancies between the observed L/LEdd distributions of type-1 and type-2 AGN. These discrepancies in the L/LEdd distributions are present regardless of the general importance of radiation forces, which must be negligible only for a small fraction of sources with large L/LEdd. Average radiation pressure corrections should then be applied in virial MBH estimators until their dependence on observed source physical properties has been fully calibrated. Finally, the comparison between MBH and L/LEdd distributions derived from sigma-based and virial estimators can constrain the variance of BLR physical properties in AGN.
The application of the virial theorem provides a tool to estimate supermassive black hole (BH) masses in large samples of active galactic nuclei (AGN) with broad emission lines at all redshifts and luminosities, if the broad line region (BLR) is grav
itationally bound. In this paper we discuss the importance of radiation forces on BLR clouds arising from the deposition of momentum by ionizing photons. Such radiation forces counteract gravitational ones and, if not taken into account, BH masses can be severely underestimated. We provide virial relations corrected for the effect of radiation pressure and we discuss their physical meaning and application. If these corrections to virial masses, calibrated with low luminosity objects, are extrapolated to high luminosities then the BLRs of most quasars might be gravitationally unbound. The importance of radiation forces in high luminosity objects must be thoroughly investigated to assess the reliability of quasar BH masses.
