Black Holes In Galaxies Research Articles

Physical Properties of the First Quasars

AbstractSince the beginning of the new millennium, more than 100 z ~ 6 quasars have been discovered through several surveys and followed-up with multi-wavelength observations. These data provided a large amounts of information on the growth of supermassive black holes at the early epochs, the properties of quasar host galaxies and the joint formation and evolution of these massive systems. We review the properties of the highest z quasars known so far, especially focussing on some of the most recent results obtained in (sub-)millimetre bands. We discuss key observational challenges and open issues in theoretical models and highlight possible new strategies to improve our understanding of the galaxy black hole formation and evolution in the early Universe.

Tracking Galaxy Evolution Through Low-Frequency Radio Continuum Observations using SKA and Citizen-Science Research using Multi-Wavelength Data

We present a review on galaxy black hole co-evolution through merger, star formation and AGN-jet feedback. We highlight results on transitional galaxies (e.g. NGC1482, NGC6764, NGC3801, Speca, RAD-18 etc.) which has data from Giant Meterwave Radio Telescope (GMRT) and other sub-mm, IR, optical, UV and X-ray telescopes. The `smoking gun' relic-evidences of past AGN-jet feedback which is believed to have quenched star formation in transitional galaxies are still missing. Relic radio lobes, as old as a few hundred Myr, can be best detected at low radio frequencies with the GMRT, LOFAR and in future SKA. However, similar relic evidences of quasar activities, known as `Hanny's Voorwerp' discovered by Galaxy Zoo in optical data, are only around a few tens of thousand years old. More discoveries are needed to match these time-scales with time since the decline of star formation in transitional galaxies. Such faint fuzzy relic emissions in optical and angular-scale sensitive radio interferometric images can be discovered most efficiently by citizen-scientists but with a formal training. We describe RAD@home, the only Indian citizen-science research project in astronomy which takes such a modified approach. We present interesting objects, discovered from the TIFR GMRT Sky Survey (TGSS) by 69 trained citizen-scientists or e-astronomers, like relic radio lobes, episodic radio galaxies, jet-galaxy interaction, bent radio galaxy in filament etc.. This model can provide an equal opportunity of academic-growth to people even in the under-developed regions where we always need to establish our optical and radio telescopes. This can expand the research-activity of city-based research-institutes beyond their four brick walls, and alleviate various socio-economic and geo-political constraints on growth of citizens located in remote areas. #RADatHomeIndia #ABCDresearch

DISCOVERY OF A PSEUDOBULGE GALAXY LAUNCHING POWERFUL RELATIVISTIC JETS

ABSTRACT Supermassive black holes launching plasma jets at close to the speed of light, producing gamma-rays, have ubiquitously been found to be hosted by massive elliptical galaxies. Since elliptical galaxies are generally believed to be built through galaxy mergers, active galactic nuclei (AGN) launching relativistic jets are associated with the latest stages of galaxy evolution. We have discovered a pseudobulge morphology in the host galaxy of the gamma-ray AGN PKS 2004-447. This is the first gamma-ray emitter radio-loud AGN found to have been launched from a system where both the black hole and host galaxy have been actively growing via secular processes. This is evidence of an alternative black hole–galaxy co-evolutionary path to develop powerful relativistic jets, which is not merger driven.

STELLAR POPULATIONS ACROSS THE BLACK HOLE MASS–VELOCITY DISPERSION RELATION

ABSTRACT Coevolution between supermassive black holes (BH) and their host galaxies is universally adopted in models for galaxy formation. In the absence of feedback from active galactic nuclei (AGNs), simulated massive galaxies keep forming stars in the local universe. From an observational point of view, however, such coevolution remains unclear. We present a stellar population analysis of galaxies with direct BH mass measurements and the BH mass–σ relation as a working framework. We find that over-massive BH galaxies, i.e., galaxies lying above the best-fitting BH mass–σ line, tend to be older and more α-element-enhanced than under-massive BH galaxies. The scatter in the BH mass–σ–[α/Fe] plane is significantly lower than that in the standard BH mass–σ relation. We interpret this trend as an imprint of AGN feedback on the star formation histories of massive galaxies.

Unveiling early black holes with JWST

AbstractThe formation of direct collapse black hole seeds with masses ~104 − 105 ~M⊙ could help explain the assembly of supermassive black holes powering high redshift quasars. Conditions conducive to the formation of these massive initial seeds exist at high redshift. Halos hosting these massive seeds merge promptly with a nearby galaxy. These early stage mergers at high redshift produce a new class of transient galaxies that contain an accreting black hole that is over-massive compared to the newly acquired stellar component - Obese Black hole Galaxies (OBGs). During this phase, the accretion luminosity of the direct collapse black hole seed exceeds that of the acquired stellar component. Here we calculate the multi-wavelength spectrum of this short-lived OBG stage, and show that there exist unique observational signatures in long wavelengths spanning near, mid to far-infrared that should be detectable by instruments aboard the upcoming James Webb Space Telescope (JWST).

Can massive primordial black holes be produced in mild waterfall hybrid inflation?

We studied the possibility whether the massive primordial black holes (PBHs) surviving today can be produced in hybrid inflation. Though it is of great interest since such PBHs can be the candidate for dark matter or seeds of the supermassive black holes in galaxies, there have not been quantitatively complete works yet because of the non-perturbative behavior around the critical point of hybrid inflation. Therefore, combining the stochastic and δN formalism, we numerically calculated the curvature perturbations in a non-perturbative way and found, without any specific assumption of the types of hybrid inflation, PBHs are rather overproduced when the waterfall phase of hybrid inflation continues so long that the PBH scale is well enlarged and the corresponding PBH mass becomes sizable enough.

STELLAR AND BLACK HOLE MASS DENSITIES AS EMPIRICAL TRACERS OF CO-EVOLUTION SHOW LOCK-STEP GROWTH SINCE Z ∼ 3

ABSTRACT At redshifts beyond , measuring the black hole (BH) galaxy relations proves to be a difficult task. The bright light of the active galactic nuclei aggravates the deconvolution of BH and galaxy properties. However, high-redshift data on these relations are vital to understand the ways in which galaxies and BHs co-evolve and the ways in which they do not. In this work we use BH and stellar mass densities (BHMDs and SMDs) to constrain the possible co-evolution of BHs with their host galaxies since . The BHMDs are calculated from quasar luminosity functions using the Soltan argument, while we use integrals over stellar mass functions or the star-formation rate density to obtain values for the SMD. We find that both quantities grow in lock-step below redshifts of with a non-evolving BHMD to SMD ratio. A fit to the data assuming a power-law relation between the BHMD and the SMD yields exponents around unity (1.0–1.5). Up to the BHMD to SMD ratio does not show a strong evolution given the larger uncertainty in the completeness of high-redshift data sets. Our results, always applying the same analysis technique, seem to be consistent across all adopted data sets.

Bondi accretion in early-type galaxies

Accretion onto central massive black holes in galaxies is often modelled with the Bondi solution. In this paper we study a generalization of the classical Bondi accretion theory, considering the additional effects of the gravitational potential of the host galaxy, and of electron scattering in the optically thin limit. We provide a general analysis of the bias in the estimates of the Bondi radius and mass accretion rate, when adopting as fiducial values for the density and temperature at infinity the values of these quantities measured at finite distance from the central black hole. We also give general formulae to compute the correction terms of the critical accretion parameter in relevant asymptotic regimes. A full analytical discussion is presented in the case of an Hernquist galaxy, when the problem reduces to the discussion of a cubic equation, therefore allowing for more than one critical point in the accretion structure. The results are useful for observational works (especially in the case of low-luminosity systems), as well as for numerical simulations, where accretion rates are usually defined in terms of the gas properties near the

Selection bias in dynamically measured supermassive black hole samples: its consequences and the quest for the most fundamental relation

We compare the set of local galaxies having dynamically measured black holes with a large, unbiased sample of galaxies extracted from the Sloan Digital Sky Survey. We confirm earlier work showing that the majority of black hole hosts have significantly higher velocity dispersions sigma than local galaxies of similar stellar mass. We use Monte-Carlo simulations to illustrate the effect on black hole scaling relations if this bias arises from the requirement that the black hole sphere of influence must be resolved to measure black hole masses with spatially resolved kinematics. We find that this selection effect artificially increases the normalization of the Mbh-sigma relation by a factor of at least ~3; the bias for the Mbh-Mstar relation is even larger. Our Monte Carlo simulations and analysis of the residuals from scaling relations both indicate that sigma is more fundamental than Mstar or effective radius. In particular, the Mbh-Mstar relation is mostly a consequence of the Mbh-sigma and sigma-Mstar relations, and is heavily biased by up to a factor of 50 at small masses. This helps resolve the discrepancy between dynamically-based black hole-galaxy scaling relations versus those of active galaxies. Our simulations also disfavour broad distributions of black hole masses at fixed sigma. Correcting for this bias suggests that the calibration factor used to estimate black hole masses in active galaxies should be reduced to values of fvir~1. Black hole mass densities should also be proportionally smaller, perhaps implying significantly higher radiative efficiencies/black hole spins. Reducing black hole masses also reduces the gravitational wave signal expected from black hole mergers.

The 5–10 keV AGN luminosity function at 0.01 <z< 4.0

The active galactic nuclei X-ray luminosity function traces actively accreting supermassive black holes and is essential for the study of the properties of the active galactic nuclei (AGN) population, black hole evolution, and galaxy-black hole coevolution. Up to now, the AGN luminosity function has been estimated several times in soft (0.5-2 keV) and hard X-rays (2-10 keV). AGN selection in these energy ranges often suffers from identification and redshift incompleteness and, at the same time, photoelectric absorption can obscure a significant amount of the X-ray radiation. We estimate the evolution of the luminosity function in the 5-10 keV band, where we effectively avoid the absorbed part of the spectrum, rendering absorption corrections unnecessary up to NH=10^23 cm^-2. Our dataset is a compilation of six wide, and deep fields: MAXI, HBSS, XMM-COSMOS, Lockman Hole, XMM-CDFS, AEGIS-XD, Chandra-COSMOS, and Chandra-CDFS. This extensive sample of ~1110 AGN (0.01<z<4.0, 41<log L_x<46) is 98% redshift complete with 68% spectroscopic redshifts. We use Bayesian analysis to select the best parametric model from simple pure luminosity and pure density evolution to more complicated luminosity and density evolution and luminosity-dependent density evolution. We estimate the model parameters that describe best our dataset separately for each survey and for the combined sample. We show that, according to Bayesian model selection, the preferred model for our dataset is the luminosity-dependent density evolution (LDDE). Our estimation of the AGN luminosity function does not require any assumption on the AGN absorption and is in good agreement with previous works in the 2-10 keV energy band based on X-ray hardness ratios to model the absorption in AGN up to redshift three. Our sample does not show evidence of a rapid decline of the AGN luminosity function up to redshift four. [abridged]

THE CONTRIBUTION OF HOST GALAXIES TO THE INFRARED ENERGY OUTPUT OF z ≳ 5.0 QUASARS

ABSTRACT The infrared spectral energy distributions of z ≳ 5 quasars can be reproduced by combining a low-metallicity galaxy template with a standard active galactic nucleus (AGN) template. The host galaxy is represented by Haro 11, a compact, moderately low metallicity, starbursting galaxy that shares typical features of high-z galaxies. For the vast majority of z ≳ 5 quasars, the AGN contribution is well modeled by a standard empirical template with the contamination of star formation in the infrared subtracted. Together, these two templates can separate the contributions from the host galaxy and the AGN even in the case of limited data points, given that this model has only two free parameters. Using this method, we reanalyze 69 z ≳ 5 quasars with extensive Herschel observations and derive their AGN luminosities [L AGN = ( 0.78 &ndash; 27.4 ) &times; 10 13 L &CircleDot; ] , ?> infrared luminosities from star formation [ L SF,IR &lsim; ( 1.5 &ndash; 25.7 ) &times; 10 12 L &CircleDot; ] , ?> and corresponding star formation rates ( SFR &lsim; 290 &ndash; 2650 M &CircleDot; yr &minus; 1 ?> ). The average infrared luminosity from star formation and the average total AGN luminosity of the z ≳ 5 quasar sample follow the correlation defined by quasars at z &lt; 2.6. We assume that these quasar host galaxies maintain a constant average SFR (∼620 M ⊙ yr−1) during their mass assembly and estimate the stellar mass that could form prior to z ∼ 5−6 to be &lang; M &ast; &rang; &sim; ( 3 &ndash; 5 ) &times; 10 11 M &CircleDot; . ?> Combining with the black hole (BH) mass measurements, this stellar mass is adequate to establish a BH–galaxy mass ratio M BH / M &ast; ?> at 0.1%–1%, consistent with the local relation.

A sharper view of our galaxy's black hole

A sharper view of our galaxy's black hole

BLACK HOLE AND GALAXY COEVOLUTION FROM CONTINUITY EQUATION AND ABUNDANCE MATCHING

[abridged] We investigate the coevolution of galaxies and hosted supermassive black holes throughout the history of the Universe by a statistical approach based on the continuity equation and the abundance matching technique. Specifically, we present analytical solutions of the continuity equation without source term to reconstruct the supermassive black hole (BH) mass function from the AGN luminosity functions. Such an approach includes physically-motivated AGN lightcurves tested on independent datasets, which describe the evolution of the Eddington ratio and radiative efficiency from slim- to thin-disc conditions. We nicely reproduce the local estimates of the BH mass function, the AGN duty cycle as a function of mass and redshift, along with the Eddington ratio function and the fraction of galaxies with given stellar mass hosting an AGN with given Eddington ratio. We exploit the same approach to reconstruct the observed stellar mass function at different redshift from the UV and far-IR luminosity functions associated to star formation in galaxies. These results imply that the buildup of stars and BHs in galaxies occurs via in-situ processes, with dry mergers playing a marginal role at least for stellar masses < 3 10^11 M_sun and BH masses < 10^9 M_sun, where the statistical data are more secure and less biased by systematic errors. In addition, we develop an improved abundance matching technique to link the stellar and BH content of galaxies to the gravitationally dominant dark matter component. The resulting relationships constitute a testbed for galaxy evolution models, highlighting the complementary role of stellar and AGN feedback in the star formation process. Finally, the clustering properties of BHs and galaxies are found to be in full agreement with current observations, so further validating our results from the continuity equation.

Growth and activity of black holes in galaxy mergers with varying mass ratios

We study supermassive black holes (BHs) in merging galaxies, using a suite of hydrodynamical simulations with very high spatial (~10 pc) and temporal (~1 Myr) resolution, where we vary the initial mass ratio, the orbital configuration, and the gas fraction. (i) We address the question of when and why, during a merger, increased BH accretion occurs, quantifying gas inflows and BH accretion rates. (ii) We also quantify the relative effectiveness in inducing AGN activity of merger-related versus secular-related causes, by studying different stages of the encounter: the stochastic (or early) stage, the (proper) merger stage, and the remnant (or late) stage. (iii) We assess which galaxy mergers preferentially enhance BH accretion, finding that the initial mass ratio is the most important factor. (iv) We study the evolution of the BH masses, finding that the BH mass contrast tends to decrease in minor mergers and to increase in major mergers. This effect hints at the existence of a preferential range of mass ratios for BHs in the final pairing stages. (v) In both merging and dynamically quiescent galaxies, the gas accreted by the BH is not necessarily the gas with $low$ angular momentum, but the gas that $loses$ angular momentum.

Deep spectroscopy of the MV ∼ −14.8 host galaxy of a tidal disruption flare in A1795★

A likely tidal disruption of a star by the intermediate-mass black hole (IMBH) of a dwarf galaxy was recently identified in association with Abell 1795. Without deep spectroscopy for this very faint object, however, the possibility of a more massive background galaxy or even a disk-instability flare from a weak AGN could not be dismissed. We have now obtained 8 hours of Gemini spectroscopy which unambiguously demonstrate that the host galaxy is indeed an extremely low-mass $(M_\ast\sim 3\times 10^8\; {\rm M}_{\odot})$ galaxy in Abell 1795, comparable to the least-massive galaxies determined to host IMBHs via other studies. We find that the spectrum is consistent with the X-ray flare being due to a tidal disruption event rather than an AGN flare. We also set improved limits on the black hole mass $({\rm log}[M_{\bullet}/{\rm M}_{\odot}] \sim 5.3 - 5.7)$ and infer a 15-year X-ray variability of a factor of $> 10^4$. The confirmation of this galaxy-black hole system provides a glimpse into a population of galaxies that is otherwise difficult to study, due to the galaxies' low masses and intrinsic faintness, but which may be important contributors to the tidal disruption rate.

A fast and accurate method to compute the mass return from multiple stellar populations

The mass returned to the ambient medium by aging stellar populations over cosmological times sums up to a significant fraction (20% - 30% or more) of their initial mass. This continuous mass injection plays a fundamental role in phenomena such as galaxy formation and evolution, fueling of supermassive black holes in galaxies and the consequent (negative and positive) feedback phenomena, and the origin of multiple stellar populations in globular clusters. In numerical simulations the calculation of the mass return can be time consuming, since it requires at each time step the evaluation of a convolution integral over the whole star formation history, so the computational time increases quadratically with the number of time-steps. The situation can be especially critical in hydrodynamical simulations, where different grid points are characterized by different star formation histories, and the gas cooling and heating times are shorter by orders of magnitude than the characteristic stellar lifetimes. In this paper we present a fast and accurate method to compute the mass return from stellar populations undergoing arbitrarily complicated star formation histories. At each time-step the mass return is calculated from its value at the previous time, and the star formation rate over the last time-step only. Therefore in the new scheme there is no need to store the whole star formation history, and the computational time increases linearly with the number of time-steps.

Massive black holes in galaxies

Massive black holes in galaxies

Nuclear coups: dynamics of black holes in galaxy mergers

We study the dynamical evolution of supermassive black holes (BHs) in merging galaxies on scales of hundreds of kpc to 10 pc, to identify the physical processes that aid or hinder the orbital decay of BHs. We present hydrodynamical simulations of galaxy mergers with a resolution of $\leq$20 pc, chosen to accurately track the motion of the nuclei and provide a realistic environment for the evolution of the BHs. We find that, during the late stages of the merger, tidal shocks inject energy in the nuclei, causing one or both nuclei to be disrupted and leaving their BH `naked', without any bound gas or stars. In many cases, the nucleus that is ultimately disrupted is that of the larger galaxy (`nuclear coup'), as star formation grows a denser nuclear cusp in the smaller galaxy. We supplement our simulations with an analytical estimate of the orbital-decay time required for the BHs to form a binary at unresolved scales, due to dynamical friction. We find that, when a nuclear coup occurs, the time-scale is much shorter than when the secondary's nucleus is disrupted, as the infalling BH is more massive, and it also finds itself in a denser stellar environment.

Globular clusters and supermassive black holes in galaxies: further analysis and a larger sample

We explore several correlations between various large-scale galaxy properties, particularly total globular cluster population (N_GCS), the central black hole mass (M_BH), velocity dispersion (nominally sigma_e), and bulge mass (M_dyn). Our data sample of 49 galaxies, for which both N_GC and M_BH are known, is larger than used in previous discussions of these two parameters and we employ the same sample to explore all pairs of correlations. Further, within this galaxy sample we investigate the scatter in each quantity, with emphasis on the range of published values for sigma_e and effective radius (R_e). We find that these two quantities in particular are difficult to measure consistently and caution that precise intercomparison of galaxy properties involving R_e and sigma_e is particularly difficult. Using both chi^2 and Monte Carlo Markov Chain (MCMC) fitting techniques, we show that quoted observational uncertainties for all parameters are too small to represent the true scatter in the data. We find that the correlation between M_dyn and N_GC is stronger than either the M_BH-sigma_e or M_BH-N_GC relations. We suggest that this is because both the galaxy bulge population ans N_GC were fundamentally established at an early epoch during the same series of star-forming events. By contrast, although the seed for M_BH was likely formed at a similar epoch, its growth over time is less similar from galaxy to galaxy and thus less predictable.

THE UNIFICATION OF RELATIVISTIC JETS

I report about the unification of relativistic jets from compact objects. The mass range is between 1.4 and 10 billion solar masses (i.e. from neutron stars to supermassive black holes in galaxies).