The paper examines the early growth of supermassive black holes (SMBHs) in cosmological hydrodynamic simulations with different BH seeding scenarios. Employing the constrained Gaussian realization, we reconstruct the initial conditions in the large-volume BlueTides simulation and run them to $z=6$ to cross-validate that the method reproduces the first quasars and their environments. Our constrained simulations in a volume of $(15\, h^{-1}{\rm Mpc})^3$ successfully recover the evolution of large-scale structure and the stellar and BH masses in the vicinity of a $\sim10^{12}\, M_{\odot}$ halo which we identified in BlueTides at $z\sim7$ hosting a $\sim10^9\, M_{\odot}$ SMBH. Among our constrained simulations, only the ones with a low-tidal field and high-density peak in the initial conditions induce the fastest BH growth required to explain the $z>6$ quasars. We run two sets of simulations with different BH seed masses of $5\times10^3$, $5\times10^4$, and $5\times10^5\, h^{-1}M_{\odot}$, (a) with the same ratio of halo to BH seed mass and (b) with the same halo threshold mass. At $z=6$, all the SMBHs converge in mass to $\sim10^9\, M_{\odot}$ except for the one with the smallest seed in (b) undergoing critical BH growth and reaching $10^8$ -- $10^9\, M_{\odot}$, albeit with most of the growth in (b) delayed compared to set (a). The finding of eight BH mergers in the small-seed scenario (four with masses $10^4$ -- $10^6\, M_{\odot}$ at $z>12$), six in the intermediate-seed scenario, and zero in the large-seed scenario suggests that the vast BHs in the small-seed scenario merge frequently during the early phases of the growth of SMBHs. The increased BH merger rate for the low-mass BH seed and halo threshold scenario provides an exciting prospect for discriminating BH formation mechanisms with the advent of multi-messenger astrophysics and next-generation gravitational wave facilities.

We examine the scaling relations between the mass of a supermassive black hole (SMBH) and its host galaxy properties at $1.2<z<1.7$ using both observational data and simulations. Recent measurements of 32 X-ray-selected broad-line Active Galactic Nucleus (AGNs) are compared with two independent state-of-the-art efforts, including the hydrodynamic simulation MassiveBlackII (MBII) and a semi-analytic model (SAM). After applying an observational selection function to the simulations, we find that both MBII and SAM agree well with the data, in terms of the central distribution. However, the dispersion in the mass ratio between black hole mass and stellar mass is significantly more consistent with the MBII prediction ($\sim0.3~$dex), than with the SAM ($\sim0.7~$dex), even when accounting for observational uncertainties. Hence, our observations can distinguish between the different recipes adopted in the models. The mass relations in the MBII are highly dependent on AGN feedback while the relations in the SAM are more sensitive to galaxy merger events triggering nuclear activity. Moreover, the intrinsic scatter in the mass ratio of our high-$z$ sample is comparable to that observed in the local sample, all but ruling out the proposed scenario the correlations are purely stochastic in nature arising from some sort of cosmic central limit theorem. Our results support the hypothesis of AGN feedback being responsible for a causal link between the SMBH and its host galaxy, resulting in a tight correlation between their respective masses.

The Quasar Luminosity Function (QLF) offers insight into the early co-evolution of black holes and galaxies. It has been characterized observationally up to redshift $z\sim6$ with clear evidence of a double power-law shape, in contrast to the Schechter-like form of the underlying dark-matter halo mass function. We investigate a physical origin for the difference in these distributions by considering the impact of stochasticity induced by the processes that determine the quasar luminosity for a given host halo and redshift. We employ a conditional luminosity function and construct the relation between median quasar magnitude versus halo mass $M_{UV,\rm{c}}(M_{\rm{h}})$ with log-normal in luminosity scatter $\Sigma$, and duty-cycle $\epsilon_{\rm{DC}}$, and focus on high redshift $z\gtrsim4$. We show that, in order to reproduce the observed QLF, the $\Sigma=0$ abundance matching requires all of the brightest quasars to be hosted in the rarest most massive dark-matter halos (with an increasing $M_{UV,\rm{c}}/M_{\rm{h}}$ in halo mass). Conversely, for $\Sigma>0$ the brightest quasars can be over-luminous outliers hosted in relatively common dark-matter halos. In this case, the median quasar magnitude versus halo mass relation, $M_{UV,\rm{c}}$, flattens at the high-end, as expected in self-regulated growth due to feedback. We sample the parameter space of $\Sigma$ and $\epsilon_{\rm{DC}}$ and show that $M_{UV,\rm{c}}$ flattens above $M_{\rm{h}}\sim 10^{12}M_{\odot}$ for $\epsilon_{\rm{DC}}<10^{-2}$. Models with $\epsilon_{\rm{DC}}\sim1$ instead require a high mass threshold close to $M_{\rm{h}}\gtrsim10^{13}M_{\odot}$. We investigate the impact of $\epsilon_{\rm{DC}}$ and $\Sigma$ on measurements of clustering and find there is no luminosity dependence on clustering for $\Sigma>0.3$, consistent with recent observations from Subaru HSC.

What triggers collimated relativistic outflows or jets, from the centers of galaxies remains a fundamental question in astrophysics. The merging of two galaxies has been proposed to realize the conditions to successfully launch and drive such jets into the intergalactic medium. However, evidences for the operation of this mechanism are scarce. Here we report the first unambiguous detection of an ongoing merger of a narrow-line Seyfert 1 galaxy, TXS 2116$-$077, hosting a closely aligned, $\gamma$-ray emitting relativistic jet with a Seyfert 2 galaxy at a separation of $\sim$12 kpc, using the observations taken with 8.2 m Subaru telescope. Our subsequent followup observations with 10.4 m Gran Telescopio Canarias, 4.2 m William Herschel Telescope, and Chandra X-ray observatory have provided what is likely to be the first glimpse of the merging environment hosting a closely aligned relativistic jet. Our finding that the jet is considerably younger than the merger demonstrates that jet activity can be triggered by galaxy mergers and that $\gamma$-ray detected narrow-line Seyfert 1 galaxies represent the beginning phase of that activity. These results also highlight the crucial role of mergers in shaping the fate of galaxies in their cosmological evolution and are consistent with recent studies focused on the host galaxy imaging of this enigmatic class of active galactic nuclei.

Nebular emission lines associated with galactic HII regions carry information about both physical properties of the ionised gas and the source of ionising photons as well as providing the opportunity of measuring accurate redshifts and thus distances once a cosmological model is assumed. While nebular line emission has been extensively studied at lower redshift there are currently only few constraints within the epoch of reionisation (EoR, $z>6$), chiefly due to the lack of sensitive near-IR spectrographs. However, this will soon change with the arrival of the Webb Telescope providing sensitive near-IR spectroscopy covering the rest-frame UV and optical emission of galaxies in the EoR. In anticipation of Webb we combine the large cosmological hydrodynamical simulation Bluetides with photoionisation modelling to predict the nebular emission line properties of galaxies at $z=8\to 13$. We find good agreement with the, albeit limited, existing direct and indirect observational constraints on equivalent widths though poorer agreement with luminosity function constraints.

We study the alignments of satellite galaxies, and their anisotropic distribution, with respect to location and orientation of their host central galaxy in MassiveBlack-II and IllustrisTNG simulations. We find that: the shape of the satellite system in halos of mass ($> 10^{13}h^{-1}M_{\odot}$) is well aligned with the shape of the central galaxy at $z=0.06$ with the mean alignment between the major axes being $\sim \Delta \theta = 12^{\circ}$ when compared to a uniform random distribution; that satellite galaxies tend to be anisotropically distributed along the major axis of the central galaxy with a stronger alignment in halos of higher mass or luminosity; and that the satellite distribution is more anisotropic for central galaxies with lower star formation rate, which are spheroidal, and for red central galaxies.Radially we find that satellites tend to be distributed along the major axis of the shape of the stellar component of central galaxies at smaller scales and the dark matter component on larger scales. We find that the dependence of satellite anisotropy on central galaxy properties and the radial distance is similar in both the simulations with a larger amplitude in MassiveBlack-II. The orientation of satellite galaxies tends to point toward the location of the central galaxy at small scales and this correlation decreases with increasing distance, and the amplitude of satellite alignment is higher in high mass halos. However, the projected ellipticities do not exhibit a scale-dependent radial alignment, as has been seen in some observational measurements.

We examine multiple AGN systems (triples and quadruples, in particular) in the \texttt{MassiveBlackII} simulation over a redshift range of $0.06\lesssim z \lesssim 4$. We identify AGN systems (with bolometric luminosity $L_{\mathrm{bol}}>10^{42}~\mathrm{ergs/sec}$) at different scales~(defined by the maximum distance between member AGNs) to determine the AGN multiplicity functions. This is defined as the volume/ surface density of AGN systems per unit \textit{richness} $R$, the number of AGNs in a system. We find that gravitationally bound multiple AGN systems tend to populate scales of $\lesssim0.7~\mathrm{cMpc}/h$; this corresponds to angular separations of $\lesssim100~\mathrm{arcsec}$ and a line of sight velocity difference $\lesssim200~\mathrm{km/sec}$. The simulation contains $\sim 10$ and $\sim100$ triples/quadruples per $\mathrm{deg}^2$ up to depths of DESI ($g\lesssim24$) and LSST ($g\lesssim26$) imaging respectively; at least $20\%$ of these should be detectable in spectroscopic surveys. The simulated quasar ($L_{\mathrm{bol}}>10^{44}~\mathrm{ergs/sec}$) triples and quadruples predominantly exist at $1.5\lesssim z \lesssim 3$. Their members have black hole masses $10^{6.5}\lesssim M_{bh}\lesssim 10^{9}~M_{\odot}/h$ and live in separate (one central and multiple satellite) galaxies with stellar masses $10^{10}\lesssim M_{*}\lesssim 10^{12}~M_{\odot}/h$. They live in the most massive haloes (for e.g. $\sim 10^{13}~M_{\odot}/h$ at $z=2.5$; $\sim 10^{14}~M_{\odot}/h$ at $z=1$) in the simulation. Their detections provide an exciting prospect for understanding massive black hole growth and their merger rates in galaxies in the era of multi-messenger astronomy.

The Wide Field Infrared Survey Telescope (WFIRST), with its field of view and high sensitivity will make surveys of cosmological large-scale structure possible at high redshifts. We investigate the possibility of detecting Baryon Acoustic Oscillations (BAO) at redshifts $z>7.6$ for use as a standard ruler. We use data from the hydrodynamic simulation \textsc{BlueTides} in conjunction with the gigaparsec-scale Outer Rim simulation and a model for patchy reionization to create mock WFIRST High Latitude Survey grism data for \lya\ emission line selected galaxies at redshifts $z=7.4$ to $z=10$, covering 2280 square degrees. We measure the monopoles of galaxies in the mock catalogues and fit the BAO features. We find that for a line flux of $L = 7\times 10^{-17} \ {\rm erg/s/cm}^{2}$, the $5 \sigma$ detection limit for the current design, the BAO feature is partially detectable (measured in three out of four survey quadrants analysed independently). The resulting root mean square error on the angular diameter distance to $z=7.7$ is 7.9$\%$. If we improve the detection sensitivity by a factor of two (i.e. $L = 3.5\times 10^{-17} \ {\rm erg/s/cm}^{2}$), the distance error reduces to $1.4\%$. We caution that many more factors are yet to be modelled, including dust obscuration, the damping wing due to the intergalactic medium, and low redshift interlopers. If these issues do not strongly affect the results, or different observational techniques (such as use of multiple lines) can mitigate them, WFIRST or similar instruments may be able to constrain the angular diameter distance to the high redshift Universe.

Galaxy intrinsic alignments (IAs) have long been recognised as a significant contaminant to weak lensing-based cosmological inference. In this paper we seek to quantify the impact of a common modelling assumption in analytic descriptions of intrinsic alignments: that of spherically symmetric dark matter halos. Understanding such effects is important as the current generation of intrinsic alignment models are known to be limited, particularly on small scales, and building an accurate theoretical description will be essential for fully exploiting the information in future lensing data. Our analysis is based on a catalogue of 113,560 galaxies between $z=0.06-1.00$ from MassiveBlack-II, a hydrodynamical simulation of box length $100 h^{-1}$ Mpc. We find satellite anisotropy contributes at the level of $\geq 30-40\%$ to the small scale alignment correlation functions. At separations larger than $1 h^{-1}$ Mpc the impact is roughly scale-independent, inducing a shift in the amplitude of the IA power spectra of $\sim20\%$. These conclusions are consistent across the redshift range and between the MassiveBlack-II and Illustris simulations. The cosmological implications of these results are tested using a simulated likelihood analysis. Synthetic cosmic shear data is constructed with the expected characteristics (depth, area and number density) of a future LSST-like survey. Our results suggest that modelling alignments using a halo model based upon spherical symmetry could potentially induce cosmological parameter biases at the $\sim 1.5\sigma$ level for $S_8$ and $w$.

We examine the properties of the host galaxies of $z=7$ quasars using the large volume, cosmological hydrodynamical simulation BlueTides. We find that the most massive black holes and quasars are hosted by massive galaxies with stellar masses $\log(M_\ast/M_\odot)=10.8\pm0.2$, and $10.2\pm0.4$, which have a wide range of star formation rates, of $50\substack{+119 \\ -35}M_\odot/\rm{yr}$ and $19\substack{+28 \\ -12}M_\odot/\rm{yr}$, respectively. The hosts of the most massive black holes and quasars in BlueTides are generally bulge-dominated, with bulge-to-total mass ratio $B/T\simeq0.85\pm0.1$, however their morphologies are not biased relative to the overall $z=7$ galaxy sample. We find that the hosts of the most massive black holes and quasars are significantly more compact, with half-mass radii $R_{0.5}=0.29\substack{+0.15 \\ -0.10}$ kpc and $0.28\substack{+0.08 \\ -0.06}$ kpc respectively, relative to galaxies with similar masses, which have $R_{0.5}=0.5\substack{+0.3 \\ -0.2}$ kpc. We make mock James Webb Space Telescope (JWST) images of these quasars and their host galaxies. We find that distinguishing the host from the quasar emission will be possible but still challenging with JWST, due to the small sizes of quasar hosts. We find that currently observable quasar samples are biased tracers of the intrinsic black hole--stellar mass relations, following a relation that is 0.27 dex higher than that of the full galaxy sample. Finally, we find that black hole hosts are generally quite isolated. However, the most massive black holes are more likely to be found in denser environments than the typical $M_{\textrm{BH}}>10^{6.5}M_\odot$ black hole, indicating that minor mergers at least play some role in growing black holes in the early Universe.

High-$z$ AGNs hosted in gas rich galaxies are expected to grow through significantly obscured accretion phases. This may limit or bias their observability. In this work, we use \textsc{BlueTides}, a large volume cosmological simulation of galaxy formation to examine quasar obscuration for the highest-redshift ($z \geq 7$) supermassive black holes residing in the center of galaxies. We find that for the bright quasars, most of the high column density gas ($>90\%$) resides in the innermost regions of the host galaxy, (typically within $< 10$ ckpc), while the gas in the outskirts is a minor contributor to the $N_\mathrm H$. The brightest quasars can have large angular variations in galactic obscuration, over 2 orders of magnitude, where the lines of sight with the lowest obscuration are those formed via strong gas outflows driven by AGN feedback. We find that for the overall AGN population, the mean $N_\mathrm H$ is generally larger for high luminosity and BH mass, while the $N_\mathrm H$ distribution is significantly broadened, developing a low $N_\mathrm H $ wing due to the angular variations driven by the AGN outflows/feedback. The obscured fraction P($N_{\rm H} > 10^{23} {\rm cm}^{-2}$) typically range from 0.6 to 1.0 for increasing $L_{X}$ (with $L_X > 10^{43} \rm{ergs/s}$), with no clear trend of redshift evolution. With respect to the galaxy host property, we find a linear relation between $N_{\rm H}$, $M_*$ and $M_{\rm H_2}$ with $\log N_{\rm H} = (0.24 \pm 0.03) \log M_{*} + (20.7 \pm 0.3)$ and $\log N_{\rm H} = (0.47 \pm 0.03) \log M_{\rm H_2} + (18.4 \pm 0.3)$. The dust optical depth in the UV band $\tau_{\mathrm UV}$ has tight positive correlation with $N_{\rm H}$. Our dust extincted UVLF is about 1.5 dex lower than the intrinsic UVLF, implying that more than 99\% of the $z \sim 7$ AGNs are heavily dust extincted and therefore would be missed by the UV band observation.

We investigate the redshift evolution of the intrinsic alignments (IA) of galaxies in the \texttt{MassiveBlackII} (MBII) simulation. We select galaxy samples above fixed subhalo mass cuts ($M_h>10^{11,12,13}~M_{\odot}/h$) at $z=0.6$ and trace their progenitors to $z=3$ along their merger trees. Dark matter components of $z=0.6$ galaxies are more spherical than their progenitors while stellar matter components tend to be less spherical than their progenitors. The distribution of the galaxy-subhalo misalignment angle peaks at $\sim10~\mathrm{deg}$ with a mild increase with time. The evolution of the ellipticity-direction~(ED) correlation amplitude $\omega(r)$ of galaxies (which quantifies the tendency of galaxies to preferentially point towards surrounding matter overdensities) is governed by the evolution in the alignment of underlying dark matter~(DM) subhaloes to the matter density of field, as well as the alignment between galaxies and their DM subhaloes. At scales $\sim1~\mathrm{cMpc}/h$, the alignment between DM subhaloes and matter overdensity gets suppressed with time, whereas the alignment between galaxies and DM subhaloes is enhanced. These competing tendencies lead to a complex redshift evolution of $\omega(r)$ for galaxies at $\sim1~\mathrm{cMpc}/h$. At scales $>1~\mathrm{cMpc}/h$, alignment between DM subhaloes and matter overdensity does not evolve significantly; the evolution of the galaxy-subhalo misalignment therefore leads to an increase in $\omega(r)$ for galaxies by a factor of $\sim4$ from $z=3$ to $0.6$ at scales $>1~\mathrm{cMpc}/h$. The balance between competing physical effects is scale dependant, leading to different conclusions at much smaller scales($\sim0.1~\mathrm{Mpc}/h$).

During the last decades, rapid progress has been made in measurements of the rest-frame ultraviolet (UV) luminosity function (LF) for high-redshift galaxies ($z \geq 6$). The faint-end of the galaxy LF at these redshifts provides powerful constraints on different dark matter models that suppress small-scale structure formation. In this work we perform full hydrodynamical cosmological simulations of galaxy formation using an alternative DM model composed of extremely light bosonic particles ($m \sim 10^{-22}$ eV), also known as fuzzy dark matter (FDM), and examine the predictions for the galaxy stellar mass function and luminosity function at $z \geq 6$ for a range of FDM masses. We find that for FDM models with bosonic mass $m = 5\times10^{-22}$ eV, the number density of galaxies with stellar mass $\rm M_* \sim 10^7 M_{\odot}$ is suppressed by $\sim 40\%$ at z = 9, $\sim 20\%$ at z = 5, and the UV LFs within magnitude range of -16 < $M_{\rm UV}$ < -14 is suppressed by $\sim 60\%$ at $z = 9$, $\sim 20\%$ at $z = 5$ comparing to the CDM counterpart simulation. Comparing our predictions with current measurements of the faint-end LFs ($-18 \leqslant M_{\rm UV} \leqslant -14$), we find that FDM models with $m_{22} < 5\times10^{-22}$ are ruled out at $3\sigma$ confidence level. We expect that future LF measurements by James Webb Space Telescope (JWST), which will extend down to $M_{\rm UV} \sim -13$ for $z \lesssim 10$, with a survey volume that is comparable to the Hubble Ultra Deep Field (HUDF) would have the capability to constrain FDM models to $m\; \gtrsim 10^{-21}$ eV.

We present an analysis of the fundamental plane of black hole accretion, an empirical correlation of the mass of a black hole ($M$), its 5 GHz radio continuum luminosity ($\nu L_{\nu}$), and its 2-10 keV X-ray power-law continuum luminosity ($L_X$). We compile a sample of black holes with primary, direct black hole-mass measurements that also have sensitive, high-spatial-resolution radio and X-ray data. Taking into account a number of systematic sources of uncertainty and their correlations with the measurements, we use Markov chain Monte Carlo methods to fit a mass-predictor function of the form $\log(M/10^{8}\,M_{\scriptscriptstyle \odot}) = \mu_0 + \xi_{\mu R} \log(L_R / 10^{38}\,\mathrm{erg\,s^{-1}}) + \xi_{\mu X} \log(L_X / 10^{40}\,\mathrm{erg\,s^{-1}})$. Our best-fit results are $\mu_0 = 0.55 \pm 0.22$, $\xi_{\mu R} = 1.09 \pm 0.10$, and $\xi_{\mu X} = -0.59^{+0.16}_{-0.15}$ with the natural logarithm of the Gaussian intrinsic scatter in the log-mass direction $\ln\epsilon_\mu = -0.04^{+0.14}_{-0.13}$. This result is a significant improvement over our earlier mass scaling result because of the increase in active galactic nuclei sample size (from 18 to 30), improvement in our X-ray binary sample selection, better identification of Seyferts, and improvements in our analysis that takes into account systematic uncertainties and correlated uncertainties. Because of these significant improvements, we are able to consider potential influences on our sample by including all sources with compact radio and X-ray emission but ultimately conclude that the fundamental plane can empirically describe all such sources. We end with advice for how to use this as a tool for estimating black hole masses.

The most distant known quasar recently discovered by Ba\~nados et al. (2018) is at $z=7.5$ (690 Myr after the Big Bang), at the dawn of galaxy formation. We explore the host galaxy of the brightest quasar in the large volume cosmological hydrodynamic simulation BlueTides, which in Phase II has reached these redshifts. The brightest quasar in BlueTides has a luminosity of a $\sim$ few $10^{13} L_{\odot}$ and a black hole mass of $6.4 \times 10^{8} M_{\odot}$ at $z \sim 7.5$, comparable to the observed quasar (the only one in this large volume). The quasar resides in a rare halo of mass $M_{H} \sim 10^{12} M_{\odot}$ and has a host galaxy of stellar mass of $4 \times 10^{10}M_{\odot}$ with an ongoing (intrinsic) star formation rate of $\sim 80 M_{\odot} yr^{-1}$. The corresponding intrinsic UV magnitude of the galaxy is $-23.1$, which is roughly $2.7$ magnitudes fainter than the quasar's magnitude of $-25.9$. We find that the galaxy is highly metal enriched with a mean metallicity equal to the solar value. We derive quasar and galaxy spectral energy distribution (SED) in the mid and near infrared JWST bands. We predict a significant amount of dust attenuation in the rest-frame UV corresponding to $A_{1500} \sim 1.7$ giving an UV based SFR of $\sim 14 M_{\odot} yr^{-1}$. We present mock JWST images of the galaxy with and without central point source, in different MIRI and NIRCam filters. The host galaxy is detectable in NIRCam filters, but it is extremely compact ($R_{E}=0.35$ kpc). It will require JWST's exquisite sensitivity and resolution to separate the galaxy from the central point source. Finally within the FOV of the quasar in BlueTides there are two more sources that would be detectable by JWST.

Many theoretical models predict that quasar driven outflows account for the observed quenching of star formation in massive galaxies. There is growing observational evidence for quasar-launched massive outflows at z>6, while the details of outflow-host galaxy interaction remain obscure. In this paper, we study the feedback around the highest redshift quasar in the BlueTides simulation, the largest volume cosmological hydrodynamic simulation so far carried out. We present predictions for gas outflows around the brightest $z = 7.54$ quasar which hosts the most massive black hole in the simulation volume, which has grown to black hole mass $6.7\times 10^{8}{\rm M}_\odot$ consistent with the current record holder for high-z quasars. We introduce a method to identify and trace the gas outflowing from the halo. We find that the total mass of the outflow gas is about $3.6\times 10^{9}{\rm M}_\odot$, constituting 6\% of the total gas in the halo. The outflow gas contains a cold, dense molecular component with mass about $2.6\times 10^{8}{\rm M}_\odot$, that originates from the inner region of the halo, within a few kpc of the central black hole. The velocities of the outflow gas reach thousands of km/s, within which the molecular component has mass averaged outward radial velocity of $1300$ km/s, consistent with observations. The averaged outflow rate is about $200-300 {\rm M}_\odot/yr$, with the outflowing gas mainly in a hotter ($T \sim 10^7$ K) and lower density state than the average of the host halo gas.

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs.

Observations of the 21 cm line from neutral hydrogen indicate that an Epoch of Heating (EoH) might have preceded the later Epoch of Reionization (EoR). Here we study the effects on the ionization state and the thermal history of the Intergalactic Medium (IGM) during the EoH induced by different assumptions on ionizing sources in the high redshift Universe: (i) stars, (ii) X-ray binaries (XRBs), (iii) thermal bremsstrahlung of the hot Interstellar Medium (ISM), and (iv) accreting nuclear black holes (BHs). To this aim, we post-process outputs from the ($100 h^{-1}$ cMpc)$^3$ hydrodynamical simulation MassiveBlack-II with the cosmological 3D radiative transfer code CRASH, which follows the propagation of UV and X-ray photons, computing the thermal and ionization state of hydrogen and helium through the EoH. We find that stars determine the fully ionized morphology of the IGM, while the spectrally hard XRBs pave way for efficient subsequent heating and ionization by the spectrally softer ISM. With the seeding prescription in MassiveBlack-II, BHs do not contribute significantly to either ionization or heating. With only stars, most of the IGM remains in a cold state (with a median $T=11$ K at $z=10$), however, the presence of more energetic sources raises the temperature of regions around the brightest and more clustered sources above that of the CMB, opening the possibility to observing the 21 cm signal in emission.

The scaling relations between the mass of supermassive black holes ($M_{\bullet}$) and host galaxy properties (stellar mass, $M_{\star}$, and velocity dispersion, $\sigma$), provide a link between the growth of black holes (BHs) and that of their hosts. Here we investigate if and how the BH-galaxy relations are established in the high-$z$ universe using \textsc{BlueTides}, a high-resolution large volume cosmological hydrodynamic simulation. We find the $M_{\bullet}-M_{\star}$ and $M_{\bullet}-\sigma$ relations at $z=8$: $\log_{10}(M_{\bullet}) = 8.25 + 1.10 \ \log_{10}(M_{\star}/10^{11}M_{\odot})$ and $\log_{10}(M_{\bullet}) = 8.35 + 5.31 \ \log_{10}(\sigma/200kms^{-1})$ at $z=8$, both fully consistent with the local measurements. The slope of the $M_{\bullet}-\sigma$ relation is slightly steeper for high star formation rate and $M_{\star}$ galaxies while it remains unchanged as a function of Eddington accretion rate onto the BH. The intrinsic scatter in $M_{\bullet}-\sigma$ relation in all cases ($\epsilon \sim 0.4$) is larger at these redshifts than inferred from observations and larger than in $M_{\bullet}-M_{\star}$ relation ($\epsilon \sim 0.14$). We find the gas-to-stellar ratio $f=M_{\rm gas}/M_{\star}$ in the host (which can be very high at these redshifts) to have the most significant impact setting the intrinsic scatter of $M_{\bullet}-\sigma$. The scatter is significantly reduced when galaxies with high gas fractions ($\epsilon = 0.28$ as $f<10$) are excluded (making the sample more comparable to low-$z$ galaxies); these systems have the largest star formation rates and black hole accretion rates, indicating that these fast-growing systems are still moving toward the relation at these high redshifts. Examining the evolution (from $z=10$ to 8) of high mass black holes in $M_{\bullet}-\sigma$ plane confirms this trend.

We study the clustering of the highest-z galaxies (from ~ $0.1$ to a few tens Mpc scales) using the BLUETIDES simulation and compare it to current observational constraints from Hubble legacy and Hyper Suprime Cam (HSC) fields (at $z=6-7.2$). With a box length of $400$ $Mpc/h$ on each side and $0.7$ trillion particles, BLUETIDES is the largest high resolution cosmological hydrodynamic simulation to date ideally suited for studies of high-z galaxies. We find that galaxies with magnitude $m_{UV}<27.7$ have a bias ($b_g$) of $8.1\pm 1.2$ at $z=8$, and typical halo masses $M_H \gtrsim 6\times10^{10} M_{\odot}$. Given the redshift evolution between $z=8$ to $z=10$ ($b_g\propto(1+z)^{1.6}$), our inferred values of the bias and halo masses are consistent with measured angular clustering at $z \sim 6.8$ from these brighter samples. The bias of fainter galaxies (in the Hubble legacy field at $H_{160} \lesssim29.5$) is $5.9\pm0.9$ at $z=8$ corresponding to halo masses $M_H \gtrsim 10^{10} M_{\odot}$. We investigate directly the 1-halo term inthe clustering and show that it dominates on scales $r \lesssim 0.1$ Mpc/$h$ ($\Theta \lesssim 3"$) with non-linear effect at transition scales between the 1-halo and 2-halo term affecting scales 0.1 $\lesssim r \lesssim $ 20 Mpc/$h$ ($3"\lesssim \Theta \lesssim 90"$). Current clustering measurements probe down to the scales in the transition between 1-halo to 2-halo regime where non-linear effects are important. The amplitude of the 1-halo term implies that occupation numbers for satellites in \texttt{BLUETIDES} are somewhat higher than standard HODs adopted in these analyses (which predict amplitudes in the 1-halo regime suppressed by a factor 2-3). That possibly implies a higher number of galaxies detected by JWST (at small scales and even fainter magnitudes) observing these fields.

Motivated by recent observational constraints on dust reprocessed emission in star forming galaxies at $z\sim 6$ and above we use the very-large cosmological hydrodynamical simulation \bluetides\ to explore predictions for the amount of dust obscured star formation in the early Universe ($z>8$). \bluetides\ matches current observational constraints on both the UV luminosity function and galaxy stellar mass function and predicts that approximately $90\%$ of the star formation in high-mass ($M_{*}>10^{10}\,{\rm M_{\odot}}$) galaxies at $z=8$ is already obscured by dust. The relationship between dust attenuation and stellar mass predicted by \bluetides\ is consistent with that observed at lower redshift. However, observations of several individual objects at $z>6$ are discrepant with the predictions, though it is possible their uncertainties may have been underestimated. We find that the predicted surface density of $z\ge 8$ sub-mm sources is below that accessible to current {\em Herschel}, SCUBA-2, and ALMA sub-mm surveys. However, as ALMA continues to accrue additional surface area the population of $z>8$ dust-obscured galaxies may become accessible in the near future.

Supermassive blackholes with masses of a billion solar masses or more are known to exist up to $z=7$. However, the present-day environments of the descendants of first quasars is not well understood and it is not known if they live in massive galaxy clusters or more isolated galaxies at $z=0$. We use a dark matter-only realization (BTMassTracer) of the BlueTides cosmological hydrodynamic simulation to study the halo properties of the descendants of the most massive black holes at $z=8$. We find that the descendants of the quasars with most massive black holes are not amongst the most massive halos. They reside in halos of with group-like ($\sim 10^{14}M_{\odot}$) masses, while the most massive halos in the simulations are rich clusters with masses $\sim 10^{15} M_{\odot}$. The distribution of halo masses at low redshift is similar to that of the descendants of least massive black holes, for a similar range of halo masses at $z=8$, which indicates that they are likely to exist in similar environments. By tracing back to the $z = 8$ progenitors of the most massive (cluster sized) halos at $z=0$; we find that their most likely black hole mass is less than $10^7 M_{\odot}$; they are clearly not amongst the most massive black holes. We also provide estimates for the likelihood of finding a high redshift quasar hosting a black hole with masses above $10^{7} M_{\odot}$ for a given halo mass at $z=0$. For halos above $10^{15} M_{\odot}$, there is only $20 \%$ probability that their $z=8$ progenitors hosted a black hole with mass above $10^{7} M_{\odot}$.

A strong correlation has been measured between the observed centripetal accelerations in galaxies and the accelerations implied by the baryonic components of galaxies. This empirical radial acceleration relation must be accounted for in any viable model of galaxy formation. We measure and compare the radial accelerations contributed by baryons and by dark matter in disk galaxies in the MassiveBlack-II hydrodynamic galaxy formation simulation. The sample of 1594 galaxies spans three orders of magnitude in luminosity and four in surface brightness, comparable to the observed sample from the Spitzer Photometry & Accurate Rotation Curves (SPARC) dataset used by McGaugh et al. (2016). We find that radial accelerations contributed by baryonic matter only and by total matter are highly correlated, with only small scatter around their mean or median relation, despite the wide ranges of galaxy luminosity and surface brightness. We further find that the radial acceleration relation in this simulation differs from that of the SPARC sample, and can be described by a simple power law in the acceleration range we are probing.

We study the shapes and intrinsic alignments of disks and elliptical galaxies in the MassiveBlack-II (MBII) and Illustris cosmological hydrodynamic simulations, with volumes of $(100h^{-1}Mpc)^{3}$ and $(75h^{-1}Mpc)^{3}$ respectively. We find that simulated disk galaxies are more oblate in shape and more misaligned with the shape of their host dark matter subhalo when compared with ellipticals. The disk major axis is found to be oriented towards the location of nearby elliptical galaxies. We also find that the disks are thinner in MBII and misalignments with dark matter halo orientations are smaller in both disks and ellipticals when compared with Illustris. As a result, the intrinsic alignment correlation functions at fixed mass have a higher amplitude in MBII than in Illustris. Finally, at scales above $\sim 0.1$\hmpc, the intrinsic alignment two-point correlation functions for disk galaxies in both simulations are consistent with a null detection, unlike those for ellipticals. Despite significant differences in the treatments of hydrodynamics and baryonic physics in the simulations, we find that the $w_{\delta +}$ correlation function scales similarly with transverse separation. However, the less massive galaxies show different scale dependence in the ED correlation. This result indicates that, while hydrodynamic simulations are a promising tool to study intrinsic alignments, further study is needed to understand the impact of differences in the implementations of hydrodynamics and baryonic feedback.

The growth of the most massive black holes in the early universe, consistent with the detection of highly luminous quasars at $z> 6$ implies sustained, critical accretion of material to grow and power them. Given a black hole seed scenario, it is still uncertain which conditions in the early Universe allow the fastest black hole growth. Large scale hydrodynamical cosmological simulations of structure formation allow us to explore the conditions conducive to the growth of the earliest supermassive black holes. We use the cosmological hydrodynamic simulation BlueTides, which incorporates a variety of baryon physics in a (400 Mpc/h)^3 volume with 0.7 trillion particles to follow the earliest phases of black hole critical growth. At z=8 the most massive black holes (a handful) approach masses of 10^8 Msun with the most massive (with M_BH = 4 x 10^8 Msun ) being found in an extremely compact spheroid-dominated host galaxy. Examining the large-scale environment of hosts, we find that the initial tidal field is more important than overdensity in setting the conditions for early BH growth. In regions of low tidal fields gas accretes 'cold' onto the black hole and falls along thin, radial filaments straight into the center forming the most compact galaxies and most massive black holes at earliest times. Regions of high tidal fields instead induce larger coherent angular momenta and influence the formation of the first population of massive compact disks. The extreme early growth depends on the early interplay of high gas densities and the tidal field that shapes the mode of accretion. Mergers play a minor role in the formation of the first generation, rare massive BHs.

We use the BlueTides simulation to predict the properties of the high-$z$ galaxy and active galactic nuclei (AGN) populations for the planned 2200deg$^2$ Wide-Field Infrared Survey Telescope's (WFIRST)-AFTA High Latitude Survey (HLS). BlueTides is a cosmological hydrodynamic simulation, which incorporates a variety of baryon physics in a $(400h^{-1} \mathrm{Mpc})^3$ volume evolved to $z=8$ with 0.7 trillion particles. The galaxy luminosity functions in the simulation show good agreement with all the current observational constraints (up to $z=11$) and predicts an enhanced number of UV bright galaxies. At the proposed depth of the HLS ($m < 26.75$), BlueTides predicts $10^6$ galaxies at $z=8$ with a few up to $z\sim 15$ due to the enhanced bright end of the galaxy luminosity function. At $z=8$, galaxies in the mock HLS have specific star formation rates of $\sim 10 {\rm Gyr}^{-1}$ and ages of $\sim 80 {\rm Myr}$ (both evolving linearly with redshift) and a non-evolving mass-metallicity relation. BlueTides also predicts $\sim 10^4$ AGN in WFIRST HLS from $z=8$ out to $z\sim 14$. These AGN host black holes of $M\sim 10^6-10^8 M_\odot$ accreting close to their Eddington luminosity. Galaxies and AGN have host halo masses of $M_{halo}\sim 10^{11-12} M_\odot$ and a linear bias $b\approx 13-20$. Given the expected galaxy space densities, their high bias and large volume probed we speculate that it may be feasible for WFIRST HLS detect the Baryon Acoustic Oscillation peak in the galaxy power spectrum out to $z=8-9$.

Using deep Hubble and Spitzer observations Oesch et al. (2016) have identified a bright ($M_{\rm UV}\approx -22$) star forming galaxy candidate at $z \approx 11$. The presence of GN-$z11$ implies a number density $\sim 10^{-6}\,{\rm Mpc^{-3}}$, roughly an order of magnitude higher than the expected value based on extrapolations from lower redshift. Using the unprecedented volume and high resolution of the BlueTides cosmological hydrodynamical simulation, we study the population of luminous rare objects at $z > 10$. The luminosity function in BlueTides implies an enhanced number of massive galaxies, consistent with the observation of GN-$z11$. We find about 30 galaxies at $M_{\rm UV}\approx -22$ at $z = 11$ in the BlueTides volume, including a few objects about 1.5 magnitudes brighter. The probability of observing GN-$z11$ in the volume probed by Oesch et al. (2016) is $\sim 13$ per cent. The predicted properties of the rare bright galaxies at $z = 11$ in BlueTides closely match those inferred from the observations of GN-$z11$. BlueTides predicts a negligible contribution from faint AGN in the observed SED. The enormous increase in volume surveyed by WFIRST will provide observations of $\sim1000$ galaxies with $M_{\rm UV} < -22$ beyond $z = 11$ out to $z = 13.5$.

The Lyman Continuum photon production efficiency ($\xi_{\rm ion}$) is a critical ingredient for inferring the number of photons available to reionise the intergalactic medium. To estimate the theoretical production efficiency in the high-redshift Universe we couple the BlueTides cosmological hydrodynamical simulation with a range of stellar population synthesis models. We find Lyman Continuum photon production efficiencies of $\log_{10}(\xi_{\rm ion}/{\rm erg^{-1}\, Hz})\approx 25.1-25.5$ depending on the choice of stellar population synthesis model. These results are broadly consistent with recent observational constraints at high-redshift though favour a model incorporating the effects of binary evolution

We examine the luminosity function (LF) of [OII] emission-line galaxies in the high-resolution cosmological simulation MassiveBlack-II (MBII). From the spectral energy distribution of each galaxy, we select a sub-sample of star-forming galaxies at $0.06 \le z \le 3.0$ using the [OII] emission line luminosity L([OII]). We confirm that the specific star formation rate matches that in the GAMA survey. We show that the [OII] LF at z=1.0 from the MBII shows a good agreement with the LFs from several surveys below L([OII])=$10^{43.0}$ erg/s while the low redshifts ($z \le 0.3$) show an excess in the prediction of bright [OII] galaxies, but still displaying a good match with observations below L([OII])=$10^{41.6}$ erg/s. Based on the validity in reproducing the properties of [OII] galaxies at low redshift ($z \le 1$), we forecast the evolution of the [OII] LF at high redshift ($z \le 3$), which can be tested by upcoming surveys such as the HETDEX and DESI. The slopes of the LFs at bright and faint ends range from -3 to -2 showing minima at z=2. The slope of the bright end evolves approximately as 1/(z+1) at z=2 while the faint end evolves as ~3/(z+1) at $0.6 \le z \le 2$. In addition, a similar analysis is applied for the evolution of [OIII] LFs, which is to be explored in the forthcoming survey WFIRST-AFTA. Finally, we show that the auto-correlation function of [OII] and [OIII] emitting galaxies shows a rapid evolution from z=2 to 1.

Whether among the myriad tiny proto-galaxies there exists a population with similarities to present day galaxies is an open question. We show, using BlueTides, the first hydrodynamic simulation large enough to resolve the relevant scales, that the first massive galaxies to form are %in fact predicted to have extensive rotationally-supported disks. Although their morphology resembles in some ways Milky-way types seen at much lower redshifts, these high-redshift galaxies are smaller, denser, and richer in gas than their low redshift counterparts. From a kinematic analysis of a statistical sample of 216 galaxies at redshift $z=8-10$ we have found that disk galaxies make up 70\% of the population of galaxies with stellar mass $10^{10} M_\odot$ or greater. Cold Dark Matter cosmology therefore makes specific predictions for the population of large galaxies 500 million years after the Big Bang. We argue that wide-field satellite telescopes (e.g. WFIRST) will in the near future discover these first massive disk galaxies. The simplicity of their structure and formation history should make possible new tests of cosmology.

We compare the shapes and intrinsic alignments of galaxies in the MassiveBlack-II cosmological hydrodynamic simulation (MBII) to those in a dark matter-only (DMO) simulation performed with the same volume (100$h^{-1}$Mpc)$^{3}$, cosmological parameters, and initial conditions. Understanding the impact of baryonic physics on galaxy shapes and alignments and their relation to the dark matter distribution should prove useful to map the intrinsic alignments of galaxies from hydrodynamic to dark matter-only simulations. We find that dark matter subhalos are typically rounder in MBII, and the shapes of stellar matter in low mass galaxies are more misaligned with the shapes of the dark matter of the corresponding subhalos in the DMO simulation. At $z=0.06$, the fractional difference in the mean misalignment angle between MBII and DMO simulations varies from $\sim 28 \% - 12 \%$ in the mass range $10^{10.8} - 6.0 \times 10^{14} h^{-1}M_{\odot}$. We study the dark matter halo shapes and alignments as a function of radius, and find that while galaxies in MBII are more aligned with the inner parts of their dark matter subhalos, there is no radial trend in their alignments with the corresponding subhalo in the DMO simulation. This result highlights the importance of baryonic physics in determining the alignment of the galaxy with respect to the inner parts of the halo. Finally, we compare the ellipticity-direction (ED) correlation for galaxies to that for dark matter halos, finding that it is suppressed on all scales by stellar-dark matter misalignment. In the projected shape-density correlation ($w_{\delta+}$), which includes ellipticity weighting, this effect is partially canceled by the higher mean ellipticities of the stellar component, but differences of order $30-40\%$ remain on scales $> 1$ Mpc over a range of subhalo masses, with scale-dependent effects below $1$ Mpc.

The intrinsic alignment of galaxies with the large-scale density field is an important astrophysical contaminant in upcoming weak lensing surveys. We present detailed measurements of the galaxy intrinsic alignments and associated ellipticity-direction (ED) and projected shape ($w_{g+}$) correlation functions for galaxies in the cosmological hydrodynamic MassiveBlack-II (MB-II) simulation. We carefully assess the effects on galaxy shapes, misalignment of the stellar component with the dark matter shape and two-point statistics of iterative weighted (by mass and luminosity) definitions of the (reduced and unreduced) inertia tensor. We find that iterative procedures must be adopted for a reliable measurement of the reduced tensor but that luminosity versus mass weighting has only negligible effects. Both ED and $w_{g+}$ correlations increase in amplitude with subhalo mass (in the range of $10^{10} - 6.0\times 10^{14}h^{-1}M_{\odot}$), with a weak redshift dependence (from $z=1$ to $z=0.06$) at fixed mass. At $z \sim 0.3$, we predict a $w_{g+}$ that is in reasonable agreement with SDSS LRG measurements and that decreases in amplitude by a factor of $\sim 5$--18 for galaxies in the LSST survey. We also compared the intrinsic alignments of centrals and satellites, with clear detection of satellite radial alignments within their host halos. Finally, we show that $w_{g+}$ (using subhalos as tracers of density) and $w_{\delta+}$ (using dark matter density) predictions from the simulations agree with that of non-linear alignment models (NLA) at scales where the 2-halo term dominates in the correlations (and tabulate associated NLA fitting parameters). The 1-halo term induces a scale dependent bias at small scales which is not modeled in the NLA model.

We use the high-resolution simulation MassiveBlackII to examine scaling relations between black hole mass (MBH) and host galaxy properties (sigma, M*, and LV), finding good agreement with observational data, especially at the high-mass end. The simulations have less intrinsic scatter than observations, and the MBH-LV correlation has the largest scatter, suggesting it may the the least fundamental of the three relations. We find Gaussian scatter about all three relations, except among the highest mass galaxies, which host more massive black holes. Below z~2 the slopes for the full population remain roughly z-independent, and only steepen by 50% by z~4. The normalization of the sigma, LV relations evolve by 0.3, 0.43 dex, while the MBH correlation does not evolve to at least z~2. Testing for selection biases, we find samples selected by MBH or M* have steeper slopes than randomly selected samples. If unaccounted for, such a selection function would find faster evolution than inferred from a randomly selected sample, as objects at the highend of the relation tend to evolve more rapidly. We find a potential bias among high-LBH subsamples (tending to reside in higher mass galaxies), but these bright-AGN exhibit no intrinsic bias relative to fainter ones in equivalent-mass hosts, nor is there a significant difference between active- and inactive-samples. Finally we characterize the evolution of individual black holes along the scaling planes. Below the local relation, black holes grow faster than their host (72% of black holes 0.3 dex below the mean relation have a MBH-M* trajectory steeper than the local relation), while those above have shallower trajectories (only 14% are steeper than local). Black holes tend to grow faster than their hosts until surpassing the local relation, at which point their growth is suppressed while their hosts continue to grow, returning them to the mean relation.

We investigate the role of Compton heating in radiation-regulated accretion on to black holes from a neutral dense medium using 1D radiation-hydrodynamic simulations. We focus on the relative effects of Compton-heating and photo-heating as a function of the spectral slope {\alpha}, assuming a power-law spectrum in the energy range of 13.6 eV--100 keV. While Compton heating is dominant only close to the black hole, it can reduce the accretion rate to 0.1 % ($l \propto \dot{m}^2$ model)--0.01 % ($l \propto \dot{m}$ model) of the Bondi accretion rate when the BH radiation is hard ({\alpha} ~ 1), where $l$ and $\dot{m}$ are the luminosity and accretion rate normalised by Eddington rates, respectively. The oscillatory behaviour otherwise typically seen in simulations with {\alpha} > 1, become suppressed when {\alpha} ~ 1 only for the $l \propto \dot{m}$ model. The relative importance of the Compton heating over photo-heating decreases and the oscillatory behaviour becomes stronger as the spectrum softens. When the spectrum is soft ({\alpha} > 1.5), photo-heating prevails regardless of models making the effect of Compton heating negligible. On the scale of the ionization front, where the gas supply into the Str\"omgren sphere from large scale is regulated, photo-heating dominates. Our simulations show consistent results with the advection-dominated accretion flow ($l \propto \dot{m}^2$) where the accretion is inefficient and the spectrum is hard ({\alpha} ~ 1).

The intrinsic alignment of galaxy shapes with the large-scale density field is a contaminant to weak lensing measurements, as well as being an interesting signature of galaxy formation and evolution (albeit one that is difficult to predict theoretically). Here we investigate the shapes and relative orientations of the stars and dark matter of halos and subhalos (central and satellite) extracted from the MassiveBlack-II simulation, a state-of-the-art high resolution hydrodynamical cosmological simulation which includes stellar and AGN feedback in a volume of $(100{h^{-1}\mathrm{Mpc}})^3$. We consider redshift evolution from $z=1$ to $0.06$ and mass evolution within the range of subhalo masses, $10^{10} -6.0 \times 10^{14.0}{h^{-1}M_{\odot}}$. The shapes of the dark matter distributions are generally more round than the shapes defined by stellar matter. The projected root-mean-square (RMS) ellipticity per component for stellar matter is measured to be $e_\text{rms} = 0.28$ at $z=0.3$ for $M_{subhalo}> 10^{12.0}{h^{-1}M_{\odot}}$, which compares favourably with observational measurements. We find that the shapes of stellar and dark matter are more round for less massive subhalos and at lower redshifts. By directly measuring the relative orientation of the stellar matter and dark matter of subgroups, we find that, on average, the misalignment between the two components is larger for less massive subhalos. The mean misalignment angle varies from $\sim 30^{\circ}-10^{\circ}$ for $M \sim 10^{10} - 10^{14} {h^{-1}M_{\odot}}$ and shows a weak dependence on redshift. We also compare the misalignment angles in central and satellite subhalos at fixed subhalo mass, and find that centrals are more misaligned than satellites. We present fitting formulae for the shapes of dark and stellar matter in subhalos and also the probability distributions of misalignment angles.

(Abridged for arXiv)We investigate the properties of halos, galaxies and blackholes to z=0 in the high resolution hydrodynamical simulation MassiveBlack-II (MBII) which evolves a LCDM cosmology in a comoving volume Vbox=100(Mpc/h)^3. MBII is the highest resolution simulation of this size which includes a self-consistent model for star formation, black hole accretion and associated feedback. We provide a simulation browser web application which enables interactive search and tagging of halos, subhalos and their properties and publicly release our galaxy catalogs. Our analysis of the halo mass function (MF) in MBII reveals that baryons have strong effects, with changes in the halo abundance of 20-35% below the knee of the MF (Mhalo < 10^13.2 Msun/h at z=0) when compared to fits based on dark matter only simulations. We provide a fitting function for the halo MF out to redshift z=11 and discuss how the onset of non-universality in the MF limits the accuracy of our fit. We study the halo occupation distribution and clustering of galaxies, in particular the evolution and scale dependence of stochasticity and bias finding reasonable agreement with observational data. The shape of the cosmic spectral energy distribution predicted by MBII is consistent with observations, but lower in amplitude. The Galaxy Stellar Mass Function (GSMF) function is broadly consistent with observations at z>=2. At z<2, the population of passive low mass (for M*<10^9 Msun) galaxies in MBII makes the GSMF too steep compared to observations whereas at the high mass end (M*>10^11 Msun) galaxies hosting bright AGNs make significant contributions to the GSMF. The quasar bolometric luminosity function is also largely consistent with observations. We note however that more efficient AGN feedback (beyond simple thermal coupling used here) is likely necessary for the largest, rarest objects/clusters at low redshifts.

We use zoom-in techniques to re-simulate three high-redshift (z > 5.5) halos which host 10^9 solar mass blackholes from the ~ Gpc volume, MassiveBlack cosmological hydrodynamic simulation. We examine a number of factors potentially affecting supermassive blackhole growth at high redshift in cosmological simulations. These include numerical resolution, feedback prescriptions and formulation of smoothed particle hydrodynamics. We find that varying the size of the region over which feedback energy is deposited directly, either for fixed number of neighbours or fixed volume makes very little difference to the accretion history of blackholes. Changing mass resolution by factors of up to 64 also does not change the blackhole growth history significantly. We find that switching from the density-entropy formulation to the pressure-entropy formulation of smoothed particle hydrodynamics slightly increases the accretion rate onto blackholes. In general numerical details appear to have small effects on the main fueling mechanism for blackholes at these high redshifts. We examine the fashion by which this occurs, finding that the insensitivity to simulation technique seems to be a hallmark of the cold flow feeding picture of these high-z supermassive blackholes. We show that the gas that participates in critical accretion phases, in these massive objects at z > 6~7 is in all cases colder, denser, and forms more coherent streams than the average gas in the halo. This is also mostly the case when the blackhole accretion is feedback regulated (z < 6), however the distinction is less prominent. For our resimulated halos, cold flows appear to be a viable mechanism for forming the most massive blackholes in the early universe, occurring naturally in LambdaCDM models of structure formation. Not requiring fine tuning of numerical parameters, they seem to be physically inevitable in these objects.

We examine the role of ionizing radiation emitted from black holes (BHs) in suppressing the growth of the Rayleigh-Taylor instability (RTI) across the ionization front (I-front) that forms when the gas fueling the BH is neutral. We use radiation-hydrodynamic simulations to show that the RTI is suppressed for non-accelerating fronts on all scales resolved in our simulations. A necessary condition for the stability of the I-front is that the radius of the Str\"omgren sphere is larger than the Bondi radius. When this condition is violated the I-front collapses producing an accretion luminosity burst. Transient growth of the RTI occurs only during the accretion burst when the effective acceleration in the frame of reference of the I-front increases significantly due to the rapid expansion of the Str\"omgren sphere.

By combining optical and near-IR observations from the Hubble Space Telescope with NIR photometry from the Spitzer Space Telescope it is possible to measure the rest-frame UV-optical colours of galaxies at z=4-8. The UV-optical spectral energy distribution of star formation dominated galaxies is the result of several different factors. These include the joint distribution of stellar masses, ages, and metallicities, and the subsequent reprocessing by dust and gas in the ISM. Using a large cosmological hydrodynamical simulation we investigate the predicted spectral energy distributions of galaxies at high-redshift with a particular emphasis on assessing the potential contribution of nebular emission. We find that the average pure stellar UV-optical colour correlates with both luminosity and redshift such that galaxies at lower-redshift and higher-luminosity are typically redder. Assuming the escape fraction of ionising photons is close to zero, the effect of nebular emission is to redden the UV-optical 1500-V_w colour by, on average, 0.4 mag at z=8 declining to 0.25 mag at z=4. Young and low-metallicity stellar populations, which typically have bluer pure stellar UV-optical colours, produce larger ionising luminosities and are thus more strongly affected by the reddening effects of nebular emission. This causes the distribution of 1500-V_w colours to narrow and the trends with luminosity and redshift to weaken. The strong effect of nebular emission leaves observed-frame colours critically sensitive to the source redshift. For example, increasing the redshift by 0.1 can result in observed frame colours changing by up to ~0.6. These predictions reinforce the need to include nebular emission when modelling the spectral energy distributions of galaxies at high-redshift and also highlight the difficultly in interpreting the observed colours of individual galaxies without precise redshifts.

We calculate the distribution of HI within 750 proper kpc/h of a quasar, Lbol = 1.62e13 Lsun, powered by an SMBH, Mbh = 4.47e8 Msun, at z = 3. Our numerical model includes a cosmological hydrodynamic simulation that tracks the self consistent growth and thermal feedback of black holes calculated using GADGET-3 as well as a detailed post-processing ray tracing treatment of the non-uniform ionizing radiation field calculated using SPHRAY, which naturally accounts for the self shielding of optically thick systems. We show that the correct treatment of self shielding introduces a flattening feature into the neutral column density distribution around Log NHI = 20 and that regions with the lowest neutral fractions are not those with the highest density gas. For comparison, we solve a Ricatti equation which determines the equilibrium Hydrogen ionization fractions in the presence of a radiation field that falls off as 1/r^2 with regions above a given gas density threshold completely shielded from ionizing radiation. We demonstrate that these semi analytic models cannot reproduce the HI field calculated using SPHRAY. We conclude by comparing our models of this single proximity zone to observations by Hennawi and Prochaska of the absorption spectra of background quasars which are coincident on the sky with foreground quasars in their Quasars Probing Quasars (QPQ) series of papers. Compared to the QPQ sample, we find a factor of 3 fewer optically thick (Log NHI > 17.2) systems around our quasar, however the dark matter halo that hosts our simulated quasar, Mhalo = 5.25e12 Msun, is less massive than the typical QPQ host halo by a factor of four. Allowing for a linear scaling between halo mass, baryonic overdensity and number of absorbers, we estimate the typical host halo mass in the QPQ sample as 1.92e13 Msun.

The observed UV continuum slope of star forming galaxies is strongly affected by the presence of dust. Its observation is then a potentially valuable diagnostic of dust attenuation, particularly at high-redshift where other diagnostics are currently inaccesible. Interpreting the observed UV continuum slope in the context of dust attenuation is often achieved assuming the empirically calibrated Meurer et al. (1999) relation. Implicit in this relation is the assumption of an intrinsic UV continuum slope ($\beta=-2.23$). However, results from numerical simulations suggest that the intrinsic UV continuum slopes of high-redshift star forming galaxies are bluer than this, and moreover vary with redshift. Using values of the intrinsic slope predicted by numerical models of galaxy formation combined with a Calzetti et al. (2000) reddening law we infer UV attenuations ($A_{1500}$) $0.35-0.5\,{\rm mag}$ ($A_{V}$: $0.14-0.2\,{\rm mag}$ assuming Calzetti et al. 2000) greater than simply assuming the Meurer relation. This has significant implications for the inferred amount of dust attenuation at very-high ($z\approx 7$) redshift given current observational constraints on $\beta$, combined with the Meurer relation, suggest dust attenuation to be virtually zero in all but the most luminous systems.

We investigate the evolution of the galaxy stellar mass function at high-redshift ($z\ge 5$) using a pair of large cosmological hydrodynamical simulations: {\em MassiveBlack} and {\em MassiveBlack-II}. By combining these simulations we can study the properties of galaxies with stellar masses greater than $10^{8}\,{\rm M_{\odot}}\,h^{-1}$ and (co-moving) number densities of $\log_{10}(\phi\, [{\rm Mpc^{-3}\,dex^{-1}}\,h^{3}])>-8$. Observational determinations of the galaxy stellar mass function at very-high redshift typically assume a relation between the observed UV luminosity and stellar mass-to-light ratio which is applied to high-redshift samples in order to estimate stellar masses. This relation can also be measured from the simulations. We do this, finding two significant differences with the usual observational assumption: it evolves strongly with redshift and has a different shape. Using this relation to make a consistent comparison between galaxy stellar mass functions we find that at $z=6$ and above the simulation predictions are in good agreement with observed data over the whole mass range. Without using the correct UV luminosity and stellar mass-to-light ratio, the discrepancy would be up to two orders of magnitude for large galaxies $>10^{10}\,{\rm M_{\odot}}\,h^{-1}$. At $z=5$, however the stellar mass function for low mass $<10^{9}\,{\rm M_{\odot}}\,h^{-1}$ galaxies is overpredicted by factors of a few, consistent with the behaviour of the UV luminosity function, and perhaps a sign that feedback in the simulation is not efficient enough for these galaxies.

We use radiative transfer to study the growth of ionized regions around the brightest, z=8 quasars in a large cosmological hydrodynamic simulation that includes black hole growth and feedback (the MassiveBlack simulation). We find that in the presence of the quasar s the comoving HII bubble radii reach 10 Mpc/h after 20 My while with the stellar component alone the HII bubbles are smaller by at least an order of magnitude. Our calculations show that several features are not captured within an analytical growth model of Stromgren spheres. The X-ray photons from hard quasar spectra drive a smooth transition from fully neutral to partially neutral in the ionization front. However the transition from partially neutral to fully ionized is significantly more complex. We measure the distance to the edge of bubbles as a function of angle and use the standard deviation of these distances as a diagnostic of the isotropy of ionized regions. We find that the overlapping of nearby ionized regions from clustered halos not only increases the anisotropy, but also is the main mechanism which allows the outer radius to grow. We therefore predict that quasar ionized bubbles at this early stage in the reionization process should be both significantly larger and more irregularly shaped than bubbles around star forming galaxies. Before the star formation rate increases and the Universe fully reionizes, quasar bubbles will form the most striking and recognizable features in 21cm maps.

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

We present Chandra observations of 12 galaxies that contain supermassive black holes with dynamical mass measurements. Each galaxy was observed for 30 ksec and resulted in a total of 68 point source detections in the target galaxies including supermassive black hole sources, ultraluminous X-ray sources, and extragalactic X-ray binaries. Based on our fits of the X-ray spectra, we report fluxes, luminosities, Eddington ratios, and slope of the power-law spectrum. Normalized to the Eddington luminosity, the 2--10 keV band X-ray luminosities of the SMBH sources range from $10^{-8}$ to $10^{-6}$, and the power-law slopes are centered at $\sim2$ with a slight trend towards steeper (softer) slopes at smaller Eddington fractions, implying a change in the physical processes responsible for their emission at low accretion rates. We find 20 ULX candidates, of which six are likely ($>90%$ chance) to be true ULXs. The most promising ULX candidate has an isotropic luminosity in the 0.3--10 keV band of $1.0_{-0.3}^{+0.6} \times 10^{40}$ erg/s.

Using a new large-scale (~ 0.75 Gpc)^3 hydrodynamic cosmological simulation we investigate the growth rate of supermassive black holes in the early universe (z > 4.75). Remarkably, we find a clear peak in the typical Eddington ratio at black hole masses of 4-8 * 10^7 solar masses (typically found in halos of ~7 * 10^11 to 10^12 solar masses), independent of redshift and indicative that most of BH growth occurs in the cold-flow dominated regime. Black hole growth is by and large regulated by the evolution of gas density. The typical Eddington ratio at a given mass scales simply as cosmological density (1+z)^3 and the peak is caused by the competition between increased gas density available in more massive hosts, and a decrease due to strong AGN feedback that deprives the black hole of sufficient gas to fuel further rapid growth in the high mass end. In addition to evolution in the mean Eddington ratio, we show that the distribution of Eddington ratio among both mass-selected and luminosity-selected samples is approximately log-normal. We combine these findings into a single log-normal fitting formula for the distribution of Eddington ratios as a function of (M_BH,z). This formula can be used in analytic and semi analytic models for evolving black hole populations, predicting black hole masses of observed quasars, and, in conjunction with the observed distribution of Eddington ratios, can be used to constrain the black hole mass function.

We investigate the formation and properties of galaxies hosting z~6 quasars, in the gigaparsec scale cosmological hydrodynamical simulation: MassiveBlack, which includes a self-consistent model for star formation, black hole accretion and associated feedback. We show that the MassiveBlack reproduces current estimates of the galaxy stellar mass function z=5, 6. We find that quasar hosts in the simulation are compact gas rich systems with high star formations rates of SFR ~ 100-1000 Msun/yr consistent with observed properties of Sloan quasar hosts in the redshift range 5.5 < z < 6.5. We show that the star-forming gas in these galaxies predominantly originates from high density cold streams which are able to penetrate the halo and grow the galaxy at the center. MassiveBlack predicts a deviation from the local Mbh-sigma and Mbh-Mstar relation implying that black holes are relatively more massive for a given stellar host at these redshifts.

Using a fully cosmological hydrodynamic simulation that self-consistently incorporates the growth and feedback of supermassive black holes and the physics of galaxy formation, we examine the effects of environmental factors (e.g., local gas density, black hole feedback) on the halo occupation distribution of low luminosity active galactic nuclei (AGN). We decompose the mean occupation function into central and satellite contribution and compute the conditional luminosity functions (CLF). The CLF of the central AGN follows a log-normal distribution with the mean increasing and scatter decreasing with increasing redshifts. We analyze the light curves of individual AGN and show that the peak luminosity of the AGN has a tighter correlation with halo mass compared to instantaneous luminosity. We also compute the CLF of satellite AGN at a given central AGN luminosity. We do not see any significant correlation between the number of satellites with the luminosity of the central AGN at a fixed halo mass. We also show that for a sample of AGN with luminosity above 10^42 ergs/s the mean occupation function can be modeled as a softened step function for central AGN and a power law for the satellite population. The radial distribution of AGN inside halos follows a power law at all redshifts with a mean index of -2.33 +/- 0.08. Incorporating the environmental dependence of supermassive black hole accretion and feedback, our formalism provides a theoretical tool for interpreting current and future measurements of AGN clustering.

We use a large dark matter simulation of a LambdaCDM model to investigate the clustering and environmental dependence of the number of substructures in a halo. Focusing on redshift z=1, we find that the halo occupation distribution is sensitive at the tens of percent level to the surrounding density and to a lesser extent to asymmetry of the surrounding density distribution. We compute the autocorrelation function of halos as a function of occupation, building on the finding of Wechsler et al. (2006) and Gao and White (2007) that halos (at fixed mass) with more substructure are more clustered. We compute the relative bias as a function of occupation number at fixed mass, finding a strong relationship. At fixed mass, halos in the top 5% of occupation can have an autocorrelation function ~ 1.5-2 times higher than the mean. We also compute the bias as a function of halo mass, for fixed halo occupation. We find that for group and cluster sized halos, when the number of subhalos is held fixed, there is a strong anticorrelation between bias and halo mass. Such a relationship represents an additional challenge to the halo model.

We examine predictions for the quasar luminosity functions (QLF) and quasar clustering at high redshift (z > 4.75) using MassiveBlack, our new hydrodynamic cosmological simulation which includes a self-consistent model for black hole growth and feedback. We show that the model reproduces the Sloan QLF within observational constraints at z >= 5. We find that the high-z QLF is consistent with a redshift-independent occupation distribution of BHs among dark matter halos (which we provide) such that the evolution of the QLF follows that of the halo mass function. The sole exception is the bright-end at z=6 and 7, where BHs in high-mass halos tend to be unusually bright due to extended periods of Eddington growth caused by high density cold flows into the halo center. We further use these luminosity functions to make predictions for the number density of quasars in upcoming surveys, predicting there should be ~119+-28 (~87+-28) quasars detectable in the F125W band of the WIDE (DEEP) fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) from z=5-6, ~19+-7 (~18+-9) from z=6-7, and ~1.7+-1.5 (~1.5+-1.5) from z=7-8. We also investigate quasar clustering, finding that the correlation length is fully consistent with current constraints for Sloan quasars (r_0~17 h^{-1} Mpc at z=4 for quasars above m_i = 20.2), and grows slowly with redshift up to z=6 (r_0~22 h^{-1} Mpc). Finally, we note that the quasar clustering strength depends weakly on luminosity for low L_BH, but gets stronger at higher L_BH as the BHs are found in higher mass halos.

The increasing size of cosmological simulations has led to the need for new visualization techniques. We focus on Smoothed Particle Hydrodynamical (SPH) simulations run with the GADGET code and describe methods for visually accessing the entire simulation at full resolution. The simulation snapshots are rastered and processed on supercomputers into images that are ready to be accessed through a web interface (GigaPan). This allows any scientist with a web-browser to interactively explore simulation datasets in both in spatial and temporal dimensions, datasets which in their native format can be hundreds of terabytes in size or more. We present two examples, the first a static terapixel image of the MassiveBlack simulation, a P-GADGET SPH simulation with 65 billion particles, and the second an interactively zoomable animation of a different simulation with more than one thousand frames, each a gigapixel in size. Both are available for public access through the GigaPan web interface. We also make our imaging software publicly available.

Observations of the most distant bright quasars imply that billion solar mass supermassive black holes (SMBH) have to be assembled within the first eight hundred million years. Under our standard galaxy formation scenario such fast growth implies large gas densities providing sustained accretion at critical or supercritical rates onto an initial black hole seed. It has been a long standing question whether and how such high black hole accretion rates can be achieved and sustained at the centers of early galaxies. Here we use our new cosmological hydrodynamic simulation (MassiveBlack) covering a volume (0.75 \Gpc)^3 appropriate for studying the rare first quasars to show that steady high density cold gas flows responsible for assembling the first galaxies produce the high gas densities that lead to sustained critical accretion rates and hence rapid growth commensurate with the existence of ~10^9 solar mass black holes as early as z~7. We find that under these conditions quasar feedback is not effective at stopping the cold gas from penetrating the central regions and hence cannot quench the accretion until the host galaxy reaches M_halo > 10^{12} solar masses. This cold-flow driven scenario for the formation of quasars implies that they should be ubiquitous in galaxies in the early universe and that major (proto)galaxy mergers are not a requirement for efficient fuel supply and growth, particularly for the earliest SMBHs.

We investigate the halo occupation distribution (HOD) of black holes within a hydrodynamic cosmological simulation that directly follows black hole growth. Similar to the HOD of galaxies/subhalos, we find that the black hole occupation number can be described by the form N_BH proportional to 1+ (M_Host)^alpha where alpha evolves mildly with redshift indicating that a given mass halo (M_Host) at low redshift tends to host fewer BHs than at high redshift (as expected as a result of galaxy and BH mergers). We further divide the occupation number into contributions from black holes residing in central and satellite galaxies within a halo. The distribution of M_BH within halos tends to consist of a single massive BH (distributed about a peak mass strongly correlated with M_Host), and a collection of relatively low-mass secondary BHs, with weaker correlation with M_Host. We also examine the spatial distribution of BHs within their host halos, and find they typically follow a power-law radial distribution (i.e. much more centrally concentrated than the subhalo distribution). Finally, we characterize the host mass for which BH growth is feedback dominated (e.g. star formation quenched). We show that halos with M_Host > 3 * 10^12 M_sun have primary BHs that are feedback dominated by z~3 with lower mass halos becoming increasingly more affected at lower redshift.

We use a large N-body simulation to examine the detectability of HI in emission at redshift z ~ 1, and the constraints imposed by current observations on the neutral hydrogen mass function of galaxies at this epoch. We consider three different models for populating dark matter halos with HI, designed to encompass uncertainties at this redshift. These models are consistent with recent observations of the detection of HI in emission at z ~ 0.8. Whilst detection of 21 cm emission from individual halos requires extremely long integrations with existing radio interferometers, such as the Giant Meter Radio Telescope (GMRT), we show that the stacked 21 cm signal from a large number of halos can be easily detected. However, the stacking procedure requires accurate redshifts of galaxies. We show that radio observations of the field of the DEEP2 spectroscopic galaxy redshift survey should allow detection of the HI mass function at the 5-12 sigma level in the mass range 10^(11.4) M_sun/h < M_halo < 10^(12.5)M_sun/h, with a moderate amount of observation time. Assuming a larger noise level that corresponds to an upper bound for the expected noise for the GMRT, the detection significance for the HI mass function is still at the 1.7-3 sigma level. We find that optically undetected satellite galaxies enhance the HI emission profile of the parent halo, leading to broader wings as well as a higher peak signal in the stacked profile of a large number of halos. We show that it is in principle possible to discern the contribution of undetected satellites to the total HI signal, even though cosmic variance limitation make this challenging for some of our models.

Both types of long and short gamma ray bursts involve a stage of a hyper-Eddington accretion of hot and dense plasma torus onto a newly born black hole. The prompt gamma ray emission originates in jets at some distance from this 'central engine' and in most events is rapidly variable, having a form of spikes and subpulses. This indicates at the variable nature of the engine itself, for which a plausible mechanism is an internal instability in the accreting flow. We solve numerically the structure and evolution of the neutrino-cooled torus. We take into account the detailed treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. The models are calculated for both Schwarzschild and Kerr black holes. We find that for sufficiently large accretion rates (> 10 Msun/s for non-rotating black hole, and >1 Msun/s for rotating black hole, depending on its spin), the inner regions of the disk become opaque, while the helium nuclei are being photodissociated. The sudden change of pressure in this region leads to the development of a viscous and thermal instability, and the neutrino pressure acts similarly to the radiation pressure in sub-Eddington disks. In the case of rapidly rotating black holes, the instability is enhanced and appears for much lower accretion rates. We also find the important and possibly further destabilizing role of the energy transfer from the rotating black hole to the torus via the magnetic coupling.

We examine the clustering properties of a population of quasars drawn from fully hydrodynamic cosmological simulations that directly follow black hole growth. We find that the black hole correlation function is best described by two distinct components: contributions from BH pairs occupying the same dark matter halo ('1-halo term') which dominate at scales below 300 kpc/h, and contributions from BHs occupying separate halos ('2-halo term') which dominate at larger scales. From the 2-halo BH term we find a typical host halo mass for faint-end quasars (those probed in our simulation volumes) ranging from 10^11 to a few 10^12 solar masses from z=5 to z=1 respectively (consistent with the mean halo host mass). The BH correlation function shows a luminosity dependence as a function of redshift, though weak enough to be consistent with observational constraints. At small scales, the high resolution of our simulations allows us to probe the 1-halo clustering in detail, finding that the 1-halo term follows an approximate power law, lacking the characteristic decrease in slope at small scales found in 1-halo terms for galaxies and dark matter. We show that this difference is a direct result of a boost in the small-scale quasar bias caused by galaxies hosting multiple quasars (1-subhalo term) following a merger event, typically between a large central subgroup and a smaller, satellite subgroup hosting a relatively small black hole. We show that our predicted small-scale excess caused by such mergers is in good agreement with both the slope and amplitude indicated by recent small-scale measurements. Finally, we note the excess to be a strong function of halo mass, such that the observed excess is well matched by the multiple black holes of intermediate mass (10^7-10^8 solar masses) found in hosts of 4-8*10^11 solar masses, a range well probed by our simulations.

We investigate the predictions for the faint-end quasar luminosity function (QLF) and its evolution using fully cosmological hydrodynamic simulations which self-consistently follow star formation, black hole growth and associated feedback processes. We find remarkably good agreement between predicted and observed faint end of the optical and X-ray QLFs (the bright end is not accessible in our simulated volumes) at z < 2. At higher redshifts our simulations tend to overestimate the QLF at the faintest luminosities. We show that although the low (high) luminosity ranges of the faint-end QLF are dominated by low (high) mass black holes, a wide range of black hole masses still contributes to any given luminosity range. This is consistent with the complex lightcurves of black holes resulting from the detailed hydrodynamics followed in the simulations. Consistent with the results on the QLFs, we find good agreement for the evolution of the comoving number density (in optical, soft and hard X-ray bands) of AGN for luminosities above 10^43 erg/s. However, the luminosity density evolution from the simulation appears to imply a peak at higher redshift than constrained from hard X-ray data (but not in optical). Our predicted excess at the faintest fluxes at z >= 2 does not lead to an overestimate to the total X-ray background and its contribution is at most a factor of two larger than the unresolved fraction of the 2-8 keV background. Even though this could be explained by some yet undetected, perhaps heavily obscured faint quasar population, we show that our predictions for the faint sources at high redshifts (which are dominated by the low mass black holes) in the simulations are likely affected by resolution effects.

Black hole accretion and jet production are areas of intensive study in astrophysics. Recent work has found a relation between radio luminosity, X-ray luminosity, and black hole mass. With the assumption that radio and X-ray luminosity are suitable proxies for jet power and accretion power, respectively, a broad fundamental connection between accretion and jet production is implied. In an effort to refine these links and enhance their power, we have explored the above relations exclusively among black holes with direct, dynamical mass-measurements. This approach not only eliminates systematic errors incurred through the use of secondary mass measurements, but also effectively restricts the range of distances considered to a volume-limited sample. Further, we have exclusively used archival data from the Chandra X-ray Observatory to best isolate nuclear sources. We find log(L_R) = (4.80 +/- 0.24) + (0.78 +/- 0.27) log(M_BH) + (0.67 +/- 0.12) log(L_X), in broad agreement with prior efforts. Owing to the nature of our sample, the plane can be turned into an effective mass predictor. When the full sample is considered, masses are predicted less accurately than with the well-known M-sigma relation. If obscured AGN are excluded, the plane is potentially a better predictor than other scaling measures.

In this paper we examine the effect of X-ray and Lyalpha photons on the intergalactic medium temperature. We calculate the photon production from a population of stars and micro-quasars in a set of cosmological hydrodynamic simulations which self-consistently follow the dark matter dynamics, radiative processes as well as star formation, black hole growth and associated feedback processes. We find that, (i) IGM heating is always dominated by X-rays unless the Lyalpha photon contribution from stars in objects with mass M<10^8 Msun becomes significantly enhanced with respect to the X-ray contribution from BHs in the same halo (which we do not directly model). (ii) Without overproducing the unresolved X-ray background, the gas temperature becomes larger than the CMB temperature, and thus an associated 21 cm signal should be expected in emission, at z<11.5. We discuss how in such a scenario the transition redshift between a 21 cm signal in absorption and in emission could be used to constraint BHs accretion and associated feedback processes.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee outlines some of the outstanding questions regarding the assembly history of Massive Black Holes in the nuclei of galaxies and the revolutionary contributions anticipated in this field from low-frequency gravitational wave astronomy.

Quasar feedback has most likely a substantial but only partially understood impact on the formation of structure in the universe. A potential direct probe of this feedback mechanism is the Sunyaev-Zeldovich effect: energy emitted from quasar heats the surrounding intergalactic medium and induce a distortion in the microwave background radiation passing through the region. Here we examine the formation of such hot quasar bubbles using a cosmological hydrodynamic simulation which includes a self-consistent treatment of black hole growth and associated feedback, along with radiative gas cooling and star formation. From this simulation, we construct microwave maps of the resulting Sunyaev-Zeldovich effect around black holes with a range of masses and redshifts. The size of the temperature distortion scales approximately with black hole mass and accretion rate, with a typical amplitude up to a few micro-Kelvin on angular scales around 10 arcseconds. We discuss prospects for the direct detection of this signal with current and future single-dish and interferometric observations, including ALMA and CCAT. These measurements will be challenging, but will allow us to characterize the evolution and growth of supermassive black holes and the role of their energy feedback on galaxy formation.

We study the effect of quasar feedback on distributions of baryons and properties of intracluster medium in galaxy groups using high-resolution numerical simulations. We use the entropy-conserving Gadget code that includes gas cooling and star formation, modified to include a physically-based model of quasar feedback. For a sample of ten galaxy group-sized dark matter halos with masses in the range of 1 to $5\times 10^{13} M_{\odot}/h$, star formation is suppressed by more than 50% in the inner regions due to the additional pressure support by quasar feedback, while gas is driven from the inner region towards the outer region of the halos. As a result, the average gas density is 50% lower in the inner region and 10% higher in the outer region in the simulation, compared to a similar simulation with no quasar feedback. Gas pressure is lowered by about 40% in the inner region and higher in the outer region, while temperature and entropy are enhanced in the inner region by about 20-40%. The total group gas fraction in the two simulations generally differs by less than 10%. We also find a small change of the total thermal Sunyaev-Zeldovich distortion, leading to 10% changes in the microwave angular power spectrum at angular scales below two arcminutes.

We make use of the first high--resolution hydrodynamic simulations of structure formation which self-consistently follows the build up of supermassive black holes introduced in Di Matteo et al. (2007) to investigate the relation between black holes (BH), host halo and large--scale environment. There are well--defined relations between halo and black hole masses and between the activities of galactic nuclei and halo masses at low redshifts. A large fraction of black holes forms anti--hierarchically, with a higher ratio of black hole to halo mass at high than at low redshifts. At $z=1$, we predict group environments (regions of enhanced local density) to contain the highest mass and most active (albeit with a large scatter) BHs while the rest of the BH population to be spread over all densities from groups to filaments and voids. Density dependencies are more pronounced at high rather than low redshift. These results are consistent with the idea that gas rich mergers are likely the main regulator of quasar activity. We find star formation to be a somewhat stronger and tighter function of local density than BH activity, indicating some difference in the triggering of the latter versus the former. There exists a large number of low--mass black holes, growing slowly predominantly through accretion, which extends all the way into the most underdense regions, i.e. in voids.

We explore galaxy properties and their link with environment and clustering using a population of ~1000 galaxies formed in a high resolution hydrodynamic simulation of the Lambda CDM cosmology. At the redshift we concentrate on, z=1, the spatial resolution is 1.4 proper kpc/h and Milky-way sized disk galaxies contain ~10^5 particles within their virial radii. We include supermassive black hole accretion and feedback as well as a multiphase model for star formation. We find that a number of familiar qualitative relationships hold approximately between galaxy properties, for example, galaxies lie between two broad extremes of type, where ``late'' types tend to be smaller in size, have lower circular velocities, younger stars, higher star formation rates, larger disk to bulge ratios and lower Sersic indices than ``early types''. As in previous studies the stellar component of disk galaxies is not as rotationally supported as in observations. Bulges contain too much of the stellar mass, although disks do have scale lengths compatible with observations. The addition of black hole physics to the simulations does not appear to have an impact on the angular momentum results, nor do we find that it is affected in an identical simulation with significantly lower mass resolution. Despite this, we can profitably use the rank order of either disk to total ratio, Sersic index, or galaxy age to separate galaxies into morphological classes and examine the density-morphology relation and morphology dependence of clustering. We find that while at redshift z=0, the well known preponderance of early types in dense environments is seen, at z=1 the density-morphology relation becomes flatter and late type galaxies are even seen to have a higher clustering amplitude than early types (abridged).

We assess the impact of starburst and AGN feedback-driven winds on the CO emission from galaxy mergers, and, in particular, search for signatures of these winds in the simulated CO morphologies and emission line profiles. We do so by combining a 3D non-LTE molecular line radiative transfer code with smoothed particle hydrodynamics (SPH) simulations of galaxy mergers that include prescriptions for star formation, black hole growth, a multiphase interstellar medium (ISM), and the winds associated with star formation and black hole growth. Our main results are: (1) Galactic winds can drive outflows of masses ~10^8-10^9 Msun which may be imaged via CO emission line mapping. (2) AGN feedback-driven winds are able to drive imageable CO outflows for longer periods of time than starburst-driven winds owing to the greater amount of energy imparted to the ISM by AGN feedback compared to star formation. (3) Galactic winds can control the spatial extent of the CO emission in post-merger galaxies, and may serve as a physical motivation for the sub-kiloparsec scale CO emission radii observed in local advanced mergers. (4) Secondary emission peaks at velocities greater than the circular velocity are seen in the CO emission lines in all models. In models with winds, these high velocity peaks are seen to preferentially correspond to outflowing gas entrained in winds, which is not the case in the model without winds. The high velocity peaks seen in models without winds are typically confined to velocity offsets (from the systemic) < 1.7 times the circular velocity, whereas the models with AGN feedback-driven winds can drive high velocity peaks to ~2.5 times the circular velocity.

Notwithstanding the advent of the Gamma-ray Large Area Telescope, theoretical models predict that a significant fraction of the cosmic gamma-ray background (CGB), at the level of 20% of the currently measured value, will remain unresolved. The angular power spectrum of intensity fluctuations of the CGB contains information on its origin. We show that probing the latter from a few tens of arcmin to several degree scales, together with complementary GLAST observations of gamma-ray emission from galaxy clusters and the blazars luminosity function, can discriminate between a background that originates from unresolved blazars or cosmic rays accelerated at structure formation shocks.

We discuss a numerical model for black hole growth and its associated feedback processes that for the first time allows cosmological simulations of structure formation to self-consistently follow the build up of the cosmic population of galaxies and active galactic nuclei. Our model assumes that seed black holes are present at early cosmic epochs at the centres of forming halos. We then track their growth from gas accretion and mergers with other black holes in the course of cosmic time. For black holes that are active, we distinguish between two distinct modes of feedback, depending on the black hole accretion rate itself. Black holes that accrete at high rates are assumed to be in a `quasar regime', where we model their feedback by thermally coupling a small fraction of their bolometric luminosity to the surrounding gas. For black holes with low accretion rates, we conjecture that most of their feedback occurs in mechanical form, where AGN-driven bubbles are injected into a gaseous environment. Using our new model, we carry out TreeSPH cosmological simulations on the scales of individual galaxies to those of massive galaxy clusters, both for isolated systems and for cosmological boxes. We demonstrate that our model produces results for the black hole and stellar mass densities in broad agreement with observational constraints. We find that the black holes significantly influence the evolution of their host galaxies, changing their star formation history, their amount of cold gas, and their colours. Also, the properties of intracluster gas are affected strongly by the presence of massive black holes in the cores of galaxy clusters, leading to shallower metallicity and entropy profiles, and to a suppression of strong cooling flows. [Abridged]

We investigate the nature of CO emission from z~6 quasars by combining non-LTE radiative transfer calculations with merger-driven models of z~6 quasar formation that arise naturally in LCDM cosmological simulations. We consider four model quasars formed in 10^12-10^13 M_sun halos from different merging histories. Our main results follow. Owing to massive starbursts and funneling of dense gas into the nuclear regions of merging galaxies, the CO is highly excited and the flux density peaks between J=5-8. The CO morphology of z~6 quasars often exhibits multiple emission peaks which arise from H2 concentrations which have not yet fully coalesced. Quasars at z~6 display a large range of sightline dependent line widths such that the lines are narrowest when the rotating H2 gas associated with the quasar is viewed face-on (when L_B is largest), and broadest when the gas is seen edge-on (when L_B is lowest). Thus for all models selection effects exist such that quasars selected for optical luminosity are preferentially face-on which may result in detected CO line widths narrower than the median. The sightline averaged line width is reflective of the circular velocity (V_c) of the host halo, and ranges from sigma~300-650 km/s. For optically selected QSOs, 10-25% (halo-mass dependant) of sightlines have narrow line widths compatible with the sole CO detection at z>6, J1148+5251. When accounting for both the temporal evolution of CO line widths, as well as the redshift evolution of halo V_c, these models self-consistently account for the CO line widths of both z~2 sub-mm galaxies and QSO's. Finally, the dynamical mass derived from the sightline averaged line widths provides a good estimate of the total mass, and allows for a stellar bulge and SMBH consistent with the local M_BH-M_bulge relation. [abridged]

The discovery of luminous quasars at redshift z ~ 6 indicates the presence of supermassive black holes (SMBHs) of mass ~10^9 Msun when the Universe was less than one billion years old. This finding presents several challenges for theoretical models. Here, we present the first multi-scale simulations that, together with a self-regulated model for the SMBH growth, produce a luminous quasar at z ~ 6.5 in the LCDM paradigm. We follow the hierarchical assembly history of the most massive halo in a ~ 3 Gpc^3 volume, and find that this halo of ~ 8x 10^{12} Msun forming at z ~ 6.5 after several major mergers is able to reproduce a number of observed properties of SDSS J1148+5251, the most distant quasar detected at z =6.42 (Fan et al. 2003). Moreover, the SMBHs grow through gas accretion below the Eddington limit in a self-regulated manner owing to feedback. We find that the progenitors experience significant star formation (up to 10^4 Msun/yr) preceding the major quasar phase such that the stellar mass of the quasar host reaches 10^{12} Msun at z ~ 6.5, consistent with observations of significant metal enrichment in SDSS J1148+5251. Our results provide a viable formation mechanism for z ~ 6 quasars in the standard LCDM cosmology, and demonstrate a common, merger-driven origin for the rarest quasars and the fundamental SMBH-host correlation in a hierarchical Universe.(Abridged)

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)

We investigate the properties and evolution of accretion tori formed after the coalescence of two compact objects. At these extreme densities and temperatures, the accreting torus is cooled mainly by neutrino emission produced primarily by electron and positron capture on nucleons (beta reactions). We solve for the disc structure and its time evolution by introducing a detailed treatment of the equation of state which includes photodisintegration of helium, the condition of beta-equilibrium, and neutrino opacities. We self-consistently calculate the chemical equilibrium in the gas consisting of helium, free protons, neutrons and electron-positron pairs and compute the chemical potentials of the species, as well as the electron fraction throughout the disc. We find that, for sufficiently large accretion rates (> 10 solar masses per second), the inner regions of the disk become opaque and develop a viscous and thermal instability. The identification of this instability might be relevant for GRB observations.

We investigate the physical conditions for the growth of intermediate mass seed black holes assumed to have formed from remnants of the first generation of massive stars. We follow the collapse of high-sigma halos with Tvir > 1e4 K using cosmological, smooth-particle hydrodynamic (SPH) simulations in the standard LCDM model. During collapse of the parent halo the seed holes are incorporated through mergers into larger systems and accrete mass from the surrounding gas. We include a self-consistent treatment of star formation, black hole accretion and associated feedback processes. Even under optimistic assumptions for the seed black hole mass and for efficient merger rates, we find that seed holes in halos M<1e10 Msun never reach the conditions for critical Eddington growth. Most of the black hole growth in this regime is determined by the initial mass and the merger rates. Critical accretion rates are reached, albeit only after a significant delay, at the time of collapse z~7) for 3-4 sigma halos of M~1e11 Msun. Our results imply M_BH = 5e6 Msun (M_halo/1e11 Msun)^0.78 at the time of collapse. The required conditions of Eddington growth to explain the build-up of supermassive black holes (~1e9 Msun), as implied by Sloan quasars at z>6, are therefore hard to meet in such a scenario. Without a 'jump-start' these conditions may be only achieved in extremely rare halos with M_halo > 1e13 Msun that collapsed before z~6. The sub-Eddington regime in which black holes holes accrete at early time implies a small contribution to the reionization by miniquasar but still sufficient to cause appreciable heating of the IGM at z<15-18.

We calculate infrared spectral energy distributions (SEDs) from simulations of major galaxy mergers and study the effect of AGN and starburst driven feedback on the evolution of the SED as a function of time. We use a self-consistent three-dimensional radiative equilibrium code to calculate the emergent SEDs and to make images. To facilitate a simple description of our findings, we describe our results in reference to an approximate analytic solution for the far-IR SED. We focus mainly on the luminous infrared galaxy (LIRG) and ultraluminous infrared galaxy (ULIRG) phases of evolution. We contrast the SEDs of simulations performed with AGN feedback to simulations performed with starburst driven wind feedback. We find that the feedback processes critically determine the evolution of the SED. Changing the source of illumination (whether stellar or AGN) has virtually no impact on the reprocessed far-infrared SED. We find that AGN feedback is particularly effective at dispersing gas and rapidly injecting energy into the ISM. The observational signature of such powerful feedback is a warm SED. In general, simulations performed with starburst driven winds have colder spectra and reprocess more of their emission into the infrared, resulting in higher infrared to bolometric luminosities compared to (otherwise equivalent) simulations performed with AGN feedback. We depict our results in IRAS bands, as well as in Spitzer's MIPS bands, and in Herschel's PACS bands.

We use numerical simulations to study the kinematic structure of remnants formed from mergers of equal-mass disk galaxies. In particular, we show that remnants of dissipational mergers, which include the radiative cooling of gas, star formation, feedback from supernovae, and the growth of supermassive black holes, are smaller, rounder, have, on average, a larger central velocity dispersion, and show significant rotation compared to remnants of dissipationless mergers. The increased rotation speed of dissipational remnants owes its origin to star formation that occurs in the central regions during the galaxy merger. We have further quantified the anisotropy, three-dimensional shape, minor axis rotation, and isophotal shape of each merger remnant, finding that dissipational remnants are more isotropic, closer to oblate, have the majority of their rotation along their major axis, and are more disky than dissipationless remnants. Individual remnants display a wide variety of kinematic properties. A large fraction of the dissipational remnants are oblate isotropic rotators. Many dissipational, and all of the dissipationless, are slowly rotating and anisotropic. The remnants of gas-rich major mergers can well-reproduce the observed distribution of projected ellipticities, rotation parameter (V/\sigma)*, kinematic misalignments, Psi, and isophotal shapes. The dissipationless remnants are a poor match to this data. Our results support the merger hypothesis for the origin of low-luminosity elliptical galaxies provided that the progenitor disks are sufficiently gas-rich, however our remnants are a poor match to the bright ellipticals that are slowly rotating and uniformly boxy.

We examine X-ray emission produced from hot gas during collisions and mergers of disk galaxies. To study this process, we employ simulations that incorporate cosmologically motivated disk-galaxy models and include the effects of radiative cooling, star formation, supernova feedback, and accreting supermassive black holes. We find that during a merger, the colliding gas in the disks is shock-heated to X-ray-emitting temperatures. The X-ray luminosity is spatially extended, rises during the initial stages of the merger, and peaks when the galactic centers coalesce. When a physical model for accreting black holes is included, the resulting feedback can drive powerful winds that contribute significantly to the amount and metallicity of hot gas, both of which increase the X-ray luminosity. In terms of their stellar kinematics and structural properties, the merger remnants in our simulations resemble elliptical galaxies. We find that the X-ray luminosities of the remnants with B-band luminosities in the range L_B ~ 10^10 - 10^11 Lsun are consistent with observations, while remnants with smaller or larger masses are underluminous in X-rays. Moreover, because the majority of the merger remnants are broadly consistent with the observed scaling relations between temperature, B-band luminosity and X-ray luminosity we conclude that major mergers are a viable mechanism for producing the X-ray halos of large, luminous elliptical galaxies.

We study the effects of feedback from active galactic nuclei (AGN) on emission from molecular gas in galaxy mergers by combining hydrodynamic simulations which include black holes with a three-dimensional, non-local thermodynamic equilibrium (LTE) radiative transfer code. We find that molecular clouds entrained in AGN winds produce an extended CO morphology with significant off-nuclear emission, which may be detectable via contour mapping. Furthermore, kinematic signatures of these molecular outflows are visible in emission line profiles when the outflow has a large line of sight velocity. Our results can help interpret current and upcoming observations of luminous infrared galaxies, as well as provide a detailed test of subresolution prescriptions for supermassive black hole growth in galaxy-scale hydrodynamic simulations.

(Abridged) The violent hierarchical nature of the LCDM cosmology poses serious difficulties for the formation of disk galaxies. To help resolve these issues, we describe a new, merger-driven scenario for the cosmological formation of disk galaxies at high redshifts that supplements the standard model based on dissipational collapse.In this picture, large gaseous disks may be produced from high-angular momentum mergers of systems that are gas-dominated, i.e. M_gas/(M_gas +M_star > 0.5 at the height of the merger. Pressurization from the multiphase structure of the interstellar medium prevents the complete conversion of gas into stars during the merger, and if enough gas remains to form a disk, the remnant eventually resembles a disk galaxy. We perform numerical simulations of galaxy mergers to study how supernovae feedback strength, supermassive black hole growth and feedback, progenitor gas fraction, merger mass-ratio, and orbital geometry impact the formation of remnant disks. We find that disks can build angular momentum through mergers and the degree of rotational support of the baryons in the merger remnant is primarily related to feedback processes associated with star formation. Disk-dominated remnants are restricted to form in mergers that are gas-dominated at the time of final coalescence and gas-dominated mergers typically require extreme progenitor gas fractions (>80%). We also show that the formation of rotationally-supported stellar systems in mergers is not restricted to idealized orbits, or major or minor mergers. We suggest that the hierarchical nature of the LCDM cosmology and the physics of the interstellar gas may act together to form spiral galaxies by building the angular momentum of disks through early, gas-dominated mergers.

We investigate the effect of embedded active galactic nuclei (AGN) in galaxy mergers on the CO molecular line emission by combining non-local thermodynamic equilibrium (LTE) radiative transfer calculations with hydrodynamic simulations. We find that AGN feedback energy in gas rich galaxy mergers can contribute to large molecular outflows which may be detectable via velocity-integrated emission contour maps, as well as through kinematic features in the emission line profiles.

We compare all the available observational data on the redshift evolution of the total stellar mass and star formation rate density in the Universe with the mass and accretion rate density evolution of supermassive black holes, estimated from the hard X-ray selected luminosity function of quasars and active galactic nuclei. We find that on average black hole mass must have been higher at higher redshift for given spheroid stellar mass. Moreover, we find negative redshift evolution of the disk/irregulars to spheroid mass ratio. The total accretion efficiency is constrained to be between 0.06 and 0.12, depending on the exact value of the local SMBH mass density, and on the critical accretion rate below which radiatively inefficient accretion may take place.

The fundamental plane of black hole activity is a non-linear correlation among radio core luminosity, X-ray luminosity and mass of all accreting black holes, both of stellar mass and supermassive, found by Merloni, Heinz and Di Matteo (2003) and, independently, by Falcke, Koerding and Markoff (2004). Here we further examine a number of statistical issues related to this correlation. In particular, we discuss the issue of sample selection and quantify the bias introduced by the effect of distance in two of the correlated quantities. We demonstrate that the fundamental plane relation cannot be a distance artifact, and that its non-linearity must represent an intrinsic characteristic of accreting black holes. We also discuss possible future observational strategies to improve our understanding of this correlation.

(Abridged) We examine the evolution of the black hole mass - stellar velocity dispersion (M-sigma) relation over cosmic time using simulations of galaxy mergers that include feedback from supermassive black hole growth. We consider mergers of galaxies varying the properties of the progenitors to match those expected at redshifts z=0-6. We find that the slope of the resulting M-sigma relation is the same at all redshifts considered. For the same feedback efficiency that reproduces the observed amplitude of the M-sigma relation at z=0, there is a weak redshift-dependence to the normalization that results from an increasing velocity dispersion for a given galactic stellar mass. We develop a formalism to connect redshift evolution in the M-sigma relation to the scatter in the local relation at z=0. We show that the scatter in the local relation places severe constraints on the redshift evolution of both the normalization and slope of the M-sigma relation. Furthermore, we demonstrate that cosmic downsizing introduces a black hole mass-dependent dispersion in the M-sigma relation and that the skewness of the distribution about the locally observed M-sigma relation is sensitive to redshift evolution in the normalization and slope. In principle, these various diagnostics provide a method for differentiating between theories for producing the M-sigma relation. In agreement with existing constraints, our simulations imply that hierarchical structure formation should produce the relation with small intrinsic scatter.

(ABRIDGED) We present in detail our new 3D numerical models for the accretion of stellar winds on to Sgr A*. In our most sophisticated models, we put stars on realistic orbits around Sgr A*, include `slow' winds (300 km/s), and account for radiative cooling. We first model only one phase `fast' stellar winds (1000 km/s). For wind sources fixed in space, the accretion rate is Mdot ~ 1e-5 Msun/yr, fluctuates by < 10%, and is in a good agreement with previous models. In contrast, Mdot decreases by an order of magnitude for stars following circular orbits, and fluctuates by ~ 50%. Then we allow a fraction of stars to produce slow winds. Much of these winds cool radiatively, forming cold clumps immersed into the X-ray emitting gas. We test two orbital configurations for the stars in this scenario, an isotropic distribution and two rotating discs with perpendicular orientation. The morphology of cold gas is quite sensitive to the orbits. In both cases, however, most of the accreted gas is hot, with an almost constant Mdot ~ 3e-6 Msun/yr, consistent with Chandra observations. The cold gas accretes in intermittent, short but powerful episodes which may give rise to large amplitude variability in the luminosity of Sgr A* on time scales of 10s to 100s of years. The circularisation radii for the flows are ~ 1e3 and 1e4 Rsch, for the one and two-phase wind simulations, respectively, never forming the quasi-spherical accretion flows suggested in some previous work. Our work suggests that, averaged over time scales of 100s to 1000s of years, the radiative and mechanical luminosity of Sgr A* may be substantially higher than it is in its current state. Further improvements of the wind accretion modelling of Sgr A* will rely on improved observational constraints for the wind properties and stellar orbits.

We present an evolutionary model for starbursts, quasars, and spheroidal galaxies in which mergers between gas-rich galaxies drive nuclear inflows of gas, producing intense starbursts and feeding the buried growth of supermassive black holes (BHs) until feedback expels gas and renders a briefly visible optical quasar. The quasar lifetime and obscuring column density depend on both the instantaneous and peak luminosity of the quasar, and we determine this dependence using a large set of simulations of galaxy mergers varying host galaxy properties, orbital geometry, and gas physics. We use these fits to deconvolve observed quasar luminosity functions (LFs) and obtain the evolution of the formation rate of quasars with a certain peak luminosity, n(L_peak,z). Quasars spend extended periods of time at luminosities well below peak, and so n(L_peak) has a maximum corresponding to the 'break' in the observed LF, falling off at both brighter and fainter luminosities. From n(L_peak) and our simulation results, we obtain self-consistent fits to hard and soft X-ray and optical quasar LFs and predict many observables, including: column density distributions of optical and X-ray samples, the LF of broad-line quasars in X-ray samples and the broad-line fraction as a function of luminosity, active BH mass functions, the distribution of Eddington ratios at z~0-2, the z=0 mass function of relic BHs and total mass density of BHs, and the cosmic X-ray background. In every case, our predictions agree well with observed estimates, and unlike previous modeling attempts, we are able to reproduce them without invoking any ad hoc assumptions about source properties or distributions. We provide a library of Monte Carlo realizations of our models for comparison with observations. (Abridged)

(Abridged) Based on numerical simulations of galaxy mergers that incorporate black hole (BH) growth, we predict the faint end slope of the quasar luminosity function (QLF) and its evolution with redshift. Our simulations have yielded a new model for quasar lifetimes where the lifetime depends on both the instantaneous and peak quasar luminosities. This motivates a new interpretation of the QLF in which the bright end consists of quasars radiating at nearly their peak luminosities, but the faint end is mostly made up of quasars in less luminous phases of evolution. The faint-end QLF slope is then determined by the faint-end slope of the quasar lifetime for quasars with peak luminosities near the observed break. We determine this slope from the quasar lifetime as a function of peak luminosity, based on a large set of simulations spanning a wide variety of host galaxy, merger, BH, and ISM gas properties. Brighter peak luminosity (higher BH mass) systems undergo more violent evolution, and expel and heat gas more rapidly in the final stages of quasar evolution, resulting in a flatter faint-end slope (as these objects fall below the observed break in the QLF more rapidly). Therefore, as the QLF break luminosity moves to higher luminosities with increasing redshift, implying a larger typical quasar peak luminosity, the faint-end QLF slope flattens. From the quasar lifetime as a function of peak luminosity and this interpretation of the QLF, we predict the faint-end QLF slope and its evolution with redshift in good agreement with observations. Although BHs grow anti-hierarchically (with lower-mass BHs formed primarily at lower redshifts), the observed change in slope and differential or luminosity dependent density evolution in the QLF is completely determined by the luminosity-dependent quasar lifetime and physics of quasar feedback.

We consider implications of our new model of quasar lifetimes and light curves for the quasar luminosity function (LF) at different frequencies and redshifts. In our picture, quasars evolve rapidly and the lifetime depends on both their instantaneous and peak luminosities. The bright end of the LF traces the peak intrinsic quasar activity, but the faint end consists of quasars which are either undergoing exponential growth to much larger masses and luminosities, or are in sub-Eddington quiescent states going into or coming out of a period of peak activity. The 'break' in the observed LF corresponds directly to the maximum in the intrinsic distribution of peak luminosities, which falls off at both brighter and fainter luminosities. We study this model using simulations of galaxy mergers which successfully reproduce a wide range of observed quasar phenomena, including the observed column density distribution. By combining quasar lifetimes and the distribution of maximum quasar luminosities determined from the observed hard X-ray LF with the corresponding luminosity and host-system dependent column densities, we produce the expected soft X-ray and B-band LFs. Our predictions agree exceptionally well with the observed LFs at all observed luminosities, over the redshift range considered (z < 1), without invoking any ad hoc assumptions about an obscured population of sources. Our results also suggest that observed correlations in hard X-ray samples between the obscured fraction of quasars and luminosity can be explained in the context of our model by the expulsion of surrounding gas due to heating from accretion feedback energy as a quasar nears its peak luminosity and final black hole mass.

Based on numerical simulations of gas-rich galaxy mergers, we discuss a model in which quasar activity is tied to the self-regulated growth of supermassive black holes in galaxies. Nuclear inflow of gas attending a galaxy collision triggers a starburst and feeds black hole growth, but for most of the duration of the starburst, the black hole is heavily obscured by surrounding gas and dust which limits the visibility of the quasar, especially at optical and UV wavelengths. Eventually, feedback energy from accretion heats the gas and expels it in a powerful wind, leaving a 'dead quasar'. Between buried and dead phases there is a window during which the galaxy would be seen as a luminous quasar. Because the black hole mass, radiative output, and distribution of obscuring gas and dust all evolve strongly with time, the duration of this phase of observable quasar activity depends on both the waveband and imposed luminosity threshold. We determine the observed and intrinsic lifetimes as a function of luminosity and frequency, and calculate observable lifetimes ~10 Myr for bright quasars in the optical B-band, in good agreement with empirical estimates and much smaller than the black hole growth timescales ~100 Myr, naturally producing a substantial population of 'buried' quasars. However, observed and intrinsic energy outputs converge in the IR and hard X-ray bands as attenuation becomes weaker and chances of observation greatly increase. We obtain the distribution of column densities along sightlines in which the quasar is seen above a given luminosity, and find that our result agrees remarkably well with observed estimates of the column density distribution from the SDSS for appropriate luminosity thresholds. (Abridged)

We propose a new interpretation of the quasar luminosity function (LF), derived from physically motivated models of quasar lifetimes and light curves. In our picture, quasars evolve rapidly and their lifetime depends on both their instantaneous and peak luminosities. We study this model using simulations of galaxy mergers that successfully reproduce a wide range of observed quasar phenomena. With lifetimes inferred from the simulations, we deconvolve the observed quasar LF from the distribution of peak luminosities, and show that they differ qualitatively, unlike for the simple models of quasar lifetimes used previously. We find that the bright end of the LF traces the intrinsic peak quasar activity, but that the faint end consists of quasars which are either undergoing exponential growth to much larger masses and higher luminosities, or are in sub-Eddington quiescent states going into or coming out of a period of peak activity. The 'break' in the LF corresponds directly to the maximum in the intrinsic distribution of peak luminosities, which falls off at both brighter and fainter luminosities. Our interpretation of the quasar LF provides a physical basis for the nature and slope of the faint-end distribution, as well as the location of the break luminosity.

Sgr A* is currently being fed by winds from a cluster of gravitationally bound young mass-loosing stars. Using observational constraints on the orbits, mass loss rates and wind velocities of these stars, we numerically model the distribution of gas in the ~ 0.1--10'' region around Sgr A*. We find that radiative cooling of recently discovered slow winds leads to the formation of many cool filaments and blobs, and to a thin and rather light accretion disc of about an arcsecond scale. The disc however does not extend all the way to our inner boundary. Instead, hot X-ray emitting gas dominates the inner arcsecond. In our simulations, cool streams of gas frequently enter this region on low angular momentum orbits, and are then disrupted and heated up to the ambient hot gas temperature. The accreting gas around Sgr A* is thus two-phase, with a hot component, observable at X-ray wavelengths, and a cool component, which may be responsible for the majority of time variability of Sgr A* emission on hundred and thousand years time-scales. We obtain an accretion rate of a few times 10^-6 msun/year, consistent with Chandra estimates, but variable on time-scales even shorter than hundred years. These results strongly depend on the chosen stellar orbits and wind parameters. Further observational input is thus key to a better modelling of Sgr A* wind accretion.

We propose a model of quasar lifetimes in which observational quasar lifetimes and an intrinsic lifetime of rapid accretion are strongly distinguished by the physics of obscuration by surrounding gas and dust. Quasars are powered by gas funneled to galaxy centers, but for a large part of the accretion lifetime are heavily obscured by the large gas densities powering accretion. In this phase, starbursts and black hole growth are fueled but the quasar is buried. Eventually, feedback from accretion energy disperses surrounding gas, creating a window in which the black hole is observable optically as a quasar, until accretion rates drop below those required to maintain a quasar luminosity. We model this process and measure the unobscured and intrinsic quasar lifetimes in a hydrodynamical simulation of a major galaxy merger. The source luminosity is determined from the black hole accretion rate, calculated from local gas properties. We calculate the column density of hydrogen to the source along multiple lines of sight and use these column densities and gas metallicities to determine B-band attenuation of the source. Defining the observable quasar lifetime as the total time with an observed B-band luminosity above some limit L_B,min, we find lifetimes ~10-20 Myr for L_B,min=10^11 L_sun (M_B=-23), in good agreement with observationally determined quasar lifetimes. This is significantly smaller than the intrinsic lifetime ~100 Myr obtained if attenuation is neglected. The ratio of observed to intrinsic lifetime is also strong function of both the limiting luminosity and the observed frequency.

In the early Universe, while galaxies were still forming, black holes as massive as a billion solar masses powered quasars. Supermassive black holes are found at the centers of most galaxies today, where their masses are related to the velocity dispersions of stars in their host galaxies and hence to the mass of the central bulge of the galaxy. This suggests a link between the growth of the black holes and the host galaxies, which has indeed been assumed for a number of years. But the origin of the observed relation between black hole mass and stellar velocity dispersion, and its connection with the evolution of galaxies have remained unclear. Here we report simulations that simultaneously follow star formation and the growth of black holes during galaxy-galaxy collisions. We find that in addition to generating a burst of star formation, a merger leads to strong inflows that feed gas to the supermassive black hole and thereby power the quasar. The energy released by the quasar expels enough gas to quench both star formation and further black hole growth. This determines the lifetime of the quasar phase (approaching 100 million years) and explains the relationship between the black hole mass and the stellar velocity dispersion.

We describe techniques for incorporating feedback from star formation and black hole accretion into simulations of isolated and merging galaxies. At present, the details of these processes cannot be resolved in simulations on galactic scales. Our basic approach therefore involves forming coarse-grained representations of the properties of the interstellar medium and black hole accretion starting from basic physical assumptions, so that the impact of these effects can be included on resolved scales. We illustrate our method using a multiphase description of star-forming gas. Feedback from star formation pressurises highly overdense gas, altering its effective equation of state. We show that this allows the construction of stable galaxy models with much larger gas fractions than possible in earlier numerical work. We extend the model by including a treatment of gas accretion onto central supermassive black holes in galaxies. Assuming thermal coupling of a small fraction of the bolometric luminosity of accreting black holes to the surrounding gas, we show how this feedback regulates the growth of black holes. In gas-rich mergers of galaxies, we observe a complex interplay between starbursts and central AGN activity when the tidal interaction triggers intense nuclear inflows of gas. Once an accreting supermassive black hole has grown to a critical size, feedback terminates its further growth, and expels gas from the central region in a powerful quasar-driven wind. Our simulation methodology is therefore able to address the coupled processes of gas dynamics, star formation, and black hole accretion during the formation of galaxies.

We examine the disc--jet connection in stellar mass and supermassive black holes by investigating the properties of their compact emission in the hard X-ray and radio bands. We compile a sample of ~100 active galactic nuclei with measured mass, 5 GHz core emission, and 2-10 keV luminosity, together with 8 galactic black holes with a total of about 50 simultaneous observations in the radio and X-ray bands. Using this sample, we study the correlations between the radio (L_R) and the X-ray (L_X) luminosity and the black hole mass (M). We find that the radio luminosity is correlated with both M and L_X, at a highly significant level. We show how this result can be used to extend the standard unification by orientation scheme to encompass unification by mass and accretion rate.

We use hydrodynamical simulations to study the color transformations induced by star formation and active galactic nuclei (AGN) during major mergers of spiral galaxies. Our modeling accounts for radiative cooling, star formation, and supernova feedback. Moreover, we include a treatment of accretion onto supermassive black holes embedded in the nuclei of the merging galaxies. We assume that a small fraction of the bolometric luminosity of an accreting black hole couples thermally to surrounding gas, providing a feedback mechanism that regulates its growth. The encounter and coalescence of the galaxies triggers nuclear gas inflow which fuels both a powerful starburst and strong black hole accretion. Comparing simulations with and without black holes, we show that AGN feedback can quench star formation and accretion on a short timescale,particularly in large galaxies where the black holes can drive powerful winds once they become sufficiently massive. The color evolution of the remnant differs markedly between mergers with and without central black holes. Without AGN, gas-rich mergers lead to ellipticals which remain blue owing to residual star formation, even after more than 7 Gyrs have elapsed. In contrast, mergers with black holes produce ellipticals that redden much faster, an effect that is more pronounced in massive remnants where a nearly complete termination of star formation occurs, allowing them to redden to u-r ~ 2.3 in less than one Gyr. AGN feedback may thus be required to explain the population of extremely red massive early type-galaxies, and it appears to be an important driver in generating the observed bimodal color distribution of galaxies in the Local Universe.

We examine possible phenomenological constraints for the joint evolution of supermassive black holes (SMBH) and their host spheroids. We compare all the available observational data on the redshift evolution of the total stellar mass and star formation rate density in the Universe with the mass and accretion rate density evolution of supermassive black holes, estimated from the hard X-ray selected luminosity function of quasars and active galactic nuclei (AGN) for a given radiative efficiency, \epsilon. We assume that the ratio of the stellar mass in spheroids to the black hole mass density evolves as (1+z)^{-\alpha}, while the ratio of the stellar mass in disks+irregulars to that in spheroids evolves as (1+z)^{-\beta}, and we derive constraints on \alpha, \beta and \epsilon. We find that \alpha>0 at the more than 4-sigma level, implying a larger black hole mass at higher redshift for a given spheroid stellar mass. The favored values for \beta are typically negative, suggesting that the fraction of stellar mass in spheroids decreases with increasing redshift. This is consistent with recent determinations that show that the mass density at high redshift is dominated by galaxies with irregular morphology. In agreement with earlier work, we constrain \epsilon to be between 0.04 and 0.11, depending on the exact value of the local SMBH mass density, but almost independently of \alpha and \beta.

Rapid, hyper-Eddington accretion is likely to power the central engines of gamma-ray bursts (GRBs). In the extreme conditions of densities and temperatures the accreting torus is cooled by neutrino emission rather than by radiation. Another important cooling mechanism is the advection of energy into the central black hole. We compute the time evolution of a neutrino-dominated disc that proceeds during the burst and investigate the changes in its density and temperature. The discrimination between short and long bursts is made on the basis of the different rates of material inflow to the outer parts of the disc, thus favoring the binary merger scenario for the short GRBs and the collapsar scenario for the long ones. Within the context of the collapsar model, we also study the evolution of the photon luminosity of the remnant disc up to times of the order of 1 day, and discuss its implications for the production of emission lines in GRB spectra.

We present the two--point correlation function (2PCF) of narrow-line active galactic nuclei (AGN) selected within the First Data Release of the Sloan Digital Sky Survey. Using a sample of 13605 AGN in the redshift range 0.055 < z < 0.2, we find that the AGN auto--correlation function is consistent with the observed galaxy auto--correlation function on scales 0.2h^{-1}Mpc to >100h^{-1}Mpc. The AGN hosts trace an intermediate population of galaxies and are not detected in either the bluest (youngest) disk--dominated galaxies or many of the reddest (oldest) galaxies. We show that the AGN 2PCF is dependent on the luminosity of the narrow [OIII] emission line (L_{[OIII]}), with low L_{[OIII]} AGN having a higher clustering amplitude than high L_{[OIII]} AGN. This is consistent with lower activity AGN residing in more massive galaxies than higher activity AGN, and L_{[OIII]} providing a good indicator of the fueling rate. Using a model relating halo mass to black hole mass in cosmological simulations, we show that AGN hosted by ~ 10^{12} M_{odot} dark matter halos have a 2PCF that matches that of the observed sample. This mass scale implies a mean black hole mass for the sample of M_{BH} ~ 10^8 M_{odot}.

Fluctuations in the redshifted 21 centimeter emission from neutral hydrogen probe the epoch of reionization. We examine the observability of this signal and the impact of extragalactic foreground radio sources. We use cosmological simulations to predict the angular correlation functions of intensity fluctuations due to unresolved radio galaxies, cluster radio halos and relics and free-free emission from the interstellar and intergalactic medium at the frequencies and angular scales relevant for the proposed 21cm tomography. In accord with previous findings, the brightness temperature fluctuations due to foreground sources are much larger than those from the primary 21cm signal at all scales. In particular, diffuse cluster radio emission, which has been previously neglected, provides the most significant foreground contamination. However, we show that the contribution to the angular fluctuations at scales \theta > 1' is dominated by the spatial clustering of bright foreground sources. This excess can be removed if sources above flux levels S > 0.1 mJy are detected and removed. Hence, efficient source removal may be sufficient to allow the detection of angular fluctuations in the 21cm emission free of extragalactic foregrounds at \theta > 1 arcmin. In addition, the removal of sources above S=0.1 mJy also reduces the foreground fluctuations to roughly the same level as the 21cm signal at scales \theta < 1 arcmin. This should allow the substraction of the foreground components in frequency space, making it possible to observe in detail the topology and history of reionization.

We report on XMM-Newton spectroscopy of the low-luminosity active galaxies (LLAGN) M81 and NGC4579 both of which have known black hole masses and well-sampled spectral energy distributions (SED). The iron Kalpha line profiles from both the LLAGN can be described in terms of two components - a narrow line at 6.4 keV and a moderately broad line (FWHM = 20000 km/s) arising from highly ionized, He-like or H-like species (E \sim 6.8 keV). We interpret the broad lines arising from an accretion disk the inner edge of which is restricted to large radii (r_{in} \sim 100 r_g). However, the Eddington ratio, {L}/{L_{Edd}}, of these sources, is 3 - 4 orders of magnitude lower than that required to photo-ionize a cold disk to He-like iron. We suggest that the lines can be explained as collisionally ionized X-ray lines arising from the transition region between a hot (radiatively inefficient) flow in the inner regions and a cold disk outside r \sim 100r_g. The accretion flow geometry probed by our XMM-Newton observations is consistent with the truncated disk models proposed to explain the SED of LLAGNs.

We study the gas metallicity of quasar hosts using cosmological hydrodynamic simulations of the Lambda-cold dark matter model. Galaxy formation in the simulations is coupled with a prescription for black hole activity enabling us to study the evolution of the metal enrichment in quasar hosts and hence explore the relationship between star/spheroid formation and black hole growth/activity. We find a steep radial metallicity gradient in quasar host galaxies, with gas metallicities close to solar values in the outer parts but becoming supersolar in the center. The hosts of the rare bright quasars at z~5-6 have star formation rates of several hundred solar masses per year and halo masses of order ~10^12 solar masses. Already at these redshifts they have supersolar (Z ~2-3 solar) central metallicities, with a mild dependence of metallicity on luminosity, consistent with observed trends. The mean value of metallicity is sensitive to the assumed quasar lifetime, providing a useful new probe of this parameter. We find that lifetimes from 10^7-4x10^7yr are favored by comparison to observational data. In both the models and observations, the rate of evolution of the mean quasar metallicity as a function of redshift is generally flat out to z ~ 4-5. Beyond the observed redshift range and out to redshift z ~ 6-8, we predict a slow decline of the mean central metallicity towards solar and slightly subsolar values as we approach the epoch of the first significant star formation activity

We examine the disc--jet connection in stellar mass and supermassive black holes by investigating the properties of their compact emission in the X-ray and radio bands. We compile a sample of ~100 active galactic nuclei with measured mass, 5 GHz core emission, and 2-10 keV luminosity, together with 8 galactic black holes with a total of ~50 simultaneous observations in the radio and X-ray bands. Using this sample, we study the correlations between the radio (L_{R}) and the X-ray (L_{X}) luminosity and the black hole mass (M). We find that the radio luminosity is correlated with {\em both} M and L_{X}, at a highly significant level. In particular, we show that the sources define a ``fundamental plane'' in the three-dimensional (log L_{R},log L_{X},log M) space, given by log L_{R}=(0.60^{+0.11}_{-0.11}) log L_{X} +(0.78^{+0.11}_{-0.09}) log M + 7.33^{+4.05}_{-4.07}, with a substantial scatter of \sigma_{R}=0.88. We compare our results to the theoretical relations between radio flux, black hole mass, and accretion rate derived by Heinz and Sunyaev (2003). Such relations depend only on the assumed accretion model and on the observed radio spectral index. Therefore, we are able to show that the X-ray emission from black holes accreting at less than a few per cent of the Eddington rate is unlikely to be produced by radiatively efficient accretion, and is marginally consistent with optically thin synchrotron emission from the jet. On the other hand, models for radiatively inefficient accretion flows seem to agree well with the data.

(abridged) The observed properties of supermassive black holes suggest a fundamental link between their assembly and the formation of their host spheroids. We model the growth and activity of black holes in galaxies using LambdaCDM cosmological hydrodynamic simulations by following the evolution of the baryonic mass component in galaxy potential wells. We find that the observed steep relation between black hole mass and spheroid velocity dispersion, M_BH propto \sigma^4, is reproduced if the gas mass in bulges is linearly proportional to the black hole mass. In this model, black hole growth saturates because of the competition with star-formation and feedback, both of which determine the gas fraction available for accretion. Unless other processes also operate, we predict that the M_BH-sigma$ relation is not set in primordial structures but is fully established at low redshifts, $z \approxlt 2$, and is shallower at earlier times. We find that that central black hole masses are related to their dark matter halos simply via M_BH ~ M_DM^4/3. We assume that galaxies undergo a quasar phase with a typical lifetime, t_Q ~ 2\times 10^7 yr and show that star-formation regulated depletion of gas in spheroids can explain the decrease of the quasar population at redshift z<3 in the optical blue band. However, with the simplest assumption of a redshift independent quasar lifetime, the model overpredicts optical quasar numbers at high redshifts although it yields the observed evolution of number density of X-ray quasars over the redshift range 1 < z< 6. Finally, we find that the majority of black hole mass is assembled in galaxies by z ~ 3 and that the black hole accretion rate density peaks in rough correspondence to the star formation rate density at z ~ 4-5.

Chandra X-ray observations of the giant elliptical galaxy M87 resolve the thermal state of the hot interstellar medium into the accretion (Bondi) radius of its central 3 10^9 Msun black hole. We measure the X-ray gas temperature and density profiles and calculate the Bondi accretion rate, Mdot_Bondi \sim 0.1 Msun/yr. The X-ray luminosity of the active nucleus of M87 observed with Chandra is L_{x, 0.5-7 \keV} \sim 7 \times 10^{40}erg/s. This value is much less than the predicted nuclear luminosity, L_{Bondi} \sim 5 \times 10^{44} erg/s, for accretion at the Bondi rate with a canonical accretion radiative efficiency of 10%. If the black hole in M87 accretes at this rate it must do so at a much lower radiative efficiency than the canonical value. The multiwavelength spectrum of the nucleus is consistent with that predicted by an advection-dominated flow. However, as is likely, the X-ray nucleus is dominated by jet emission then the properties of flow must be modified, possibly by outflows. We show that the overall energetics of the system are just consistent with the predicted Bondi nuclear power. This suggests that either most of the accretion energy is released in the relativistic jet or that the central engine of M87 undergoes on-off activity cycles. We show that, at present, the energy dumped into the ISM by the jet may reduce the accretion rate onto the black hole by a factor \propto (v_j/c_s)^{-2}, where v_j is the jet velocity and c_s the ISM sound speed, and that this is sufficient to account for the low nuclear luminosity.

On two occasions, we obtained nearly simultaneous ~ 4 kilosecond snapshot observations of the Galactic black hole and microquasar XTE J1550-564 with Chandra and RXTE near the peak of its May, 2000 outburst. The low-energy sensitivity of Chandra and the resolution of the High Energy Transmission Grating Spectrometer (HETGS), coupled with the broad energy range and large collecting area of RXTE, have allowed us to place constraints on the outburst accretion flow geometry of this source in the ``intermediate'' X-ray state. The 0.65-25.0 keV continuum spectra are well-described by relatively hot (kT ~ 0.8 keV) accretion disk and hard (Gamma ~ 2.3) coronal power-law components. Broad, relatively strong Fe K-alpha emission line (EW ~170 eV) and smeared absorption edge components consistent with Fe XXV are strongly required in joint spectral fits. The resolution of the Chandra/HETGS reveals that the broad Fe K-alpha emission lines seen clearly in the individual RXTE spectra are not due to an intrinsically narrow line.

Many models of Gamma Ray Bursts invoke a central engine consisting of a black hole of a few solar masses accreting matter from a disk at a rate of a fraction to a few solar masses per second. Popham et al. and Narayan et al. have shown that, for Mdot >~ 0.1 Msun/s, accretion proceeds via neutrino cooling and neutrinos can carry away a significant amount of energy from the inner regions of the disks. We improve on these calculations by including a simple prescription for neutrino transfer and neutrino opacities in such regions. We find that the flows become optically thick to neutrinos inside a radius R~6-40R_s for Mdot in the range of 0.1 -10 Msun/s, where R_s is the black hole Schwarzchild radius. Most of the neutrino emission comes from outside this region and, the neutrino luminosity stays roughly constant at a value L_{\nu} ~ 10^{53} erg/s. We show that, for Mdot > 1 Msun/s, neutrinos are sufficiently trapped that energy advection becomes the dominant cooling mechanism in the flow. These results imply that neutrino annihilation in hyperaccreting black holes is an inefficient mechanism for liberating large amounts of energy. Extraction of rotational energy by magnetic processes remains the most viable mechanism.

We observed the Galactic black hole Cygnus X-1 with the Chandra High Energy Transmission Grating Spectrometer for 30 kiloseconds on 4 January, 2001. The source was in an intermediate state, with a flux that was approximately twice that commonly observed in its persistent low/hard state. Our best-fit model for the X-ray spectrum includes narrow Gaussian emission line (E = 6.415 +/- 0.007 keV, FWHM = 80 (+28, -19) eV, W = 16 (+3, -2) eV) and broad line (E = 5.82 (+0.06, -0.07) keV, FWHM = 1.9 (+0.5, -0.3) keV, W = 140 (+70, -40) eV) components, and a smeared edge at 7.3 +/- 0.2 keV (tau ~ 1.0). The broad line profile is not as strongly skewed as those observed in some Seyfert galaxies. We interpret these features in terms of an accretion disk with irradiation of the inner disk producing a broad Fe K-alpha emission line and edge, and irradiation of the outer disk producing a narrow Fe K-alpha emission line. The broad line is likely shaped predominantly by Doppler shifts and gravitational effects, and to a lesser degree by Compton scattering due to reflection. We discuss the underlying continuum X-ray spectrum and these line features in the context of diagnosing the accretion flow geometry in Cygnus X-1 and other Galactic black holes.

Absorption or emission against the cosmic microwave background radiation (CMB) may be observed in the redshifted 21cm line if the spin temperature of the neutral intergalactic medium prior to reionization differs from the CMB temperature. This so-called 21cm tomography should reveal important information on the physical state of the intergalactic medium at high redshifts. The fluctuations in the redshifted 21 cm, due to gas density inhomogeneities at early times, should be observed at meter wavelengths by the next generation radio telescopes such as the proposed {\it Square Kilometer Array (SKA)}. Here we show that the extra-galactic radio sources provide a serious contamination to the brightness temperature fluctuations expected in the redshifted 21 cm emission from the IGM at high redshifts. Unless the radio source population cuts off at flux levels above the planned sensitivity of SKA, its clustering noise component will dominate the angular fluctuations in the 21 cm signal. The integrated foreground signal is smooth in frequency space and it should nonetheless be possible to identify the sharp spectral feature arising from the non-uniformities in the neutral hydrogen density during the epoch when the first UV sources reionize the intergalactic medium.

We consider the X-ray emission due to bremsstrahlung processes from hot, low radiative-efficiency accretion flows around supermassive and galactic black holes. We calculate surface brightness profiles and Michelson visibility functions for a range of density profiles, rho ~ r^(-3/2+p), with 0 < p < 1, to allow for the presence of outflows. We find that although the 1 keV emitting region in these flows can always extend up to 10^6 Schwarzschild radii (R_S), their surface brightness profiles and visibility functions are strongly affected by the specific density profile. The advection-dominated solutions with no outflows (p=0) lead to centrally peaked profiles with characteristic sizes of only a few tens of R_S. Solutions with strong outflows (p~1) lead to flat intensity profiles with significantly larger characteristic sizes of up to 10^6 R_S. This implies that low luminosity galactic nuclei, such as M87, may appear as extended X-ray sources when observed with current X-ray imaging instruments. We show that X-ray brightness profiles and their associated visibility functions may be powerful probes for determining the relevant mode of accretion and, in turn, the properties of hot accretion flows. We discuss the implications of our results for observations with the Chandra X-ray Observatory and the planned X-ray interferometer MAXIM.

Chandra observations of low-luminosity supermassive black holes in nearby elliptical galaxies provide tight limits on both their nuclear luminosities and on their Bondi accretion rates. We examine Chandra constraints on NGC 6166, the dominant galaxy in the cluster Abell 2199, which hosts a \sim 10^9 \Msun black hole. We measure a nuclear X-ray luminosity L_{x} \sim 10^{40} erg/s and show that the density and temperature profiles of the hot interstellar medium imply a Bondi accretion rate of Mdot_{Bondi} \sim 3 10^{-2} Msun/yr. This accretion rate predicts a nuclear luminosity of \sim 10^{44} erg/s for a canonical radiative efficiency of 10%. Unless the Bondi estimate is inappropriate and/or the accretion rate onto the black hole is significantly reduced, the observed nuclear flux constrains the radiative efficiency of the accretion flow to be \eta \sim 10^{-5}. We show that low-radiative efficiency accretion flows with radial density profiles \rho \approxpropto r^{-3/2} (and not significantly flatter) can explain the observed nuclear X-ray luminosity but that the power output from the jets in NGC 6166 is also important to the energetics of the system.

We show how the simultaneous presence of a strong quasi periodic oscillation (QPO) of period of about 10 seconds in the optical and X-ray lightcurves of the X-ray transient XTE J1118+480 can be used to obtain information about the nature of the accretion flow around the source. The unusually high optical-to-X-ray flux ratio and the QPO observed simultaneously in both energy bands suggest that a significant fraction of the optical flux might originate close to the central source, where most of the X-rays are produced and be indicative of a magnetically dominated corona. We also show how the temporal evolution of the QPO can provide us with information on both the inner radius and the viscous properties of the optically thick accretion disc.

(Abridged) We use a hydrodynamic simulation of a LambdaCDM model to predict the extragalactic X-ray background (XRB), focussing on emission from the intergalactic medium (IGM). We also include X-rays from point sources associated with galaxies in the simulation, and make maps of the angular distribution of the emission. We find that filaments in the maps are not evident, being diluted by projection. In the soft (0.5-2 keV) band, the mean intensity of radiation from intergalactic and cluster gas is 2.3*10^-12 ergdeg^-2cm^-2s^-1, 35% of the total soft band emission. This is compatible at the ~1 sigma level with estimates of the unresolved soft background from ROSAT and {\it Chandra}. Only 4% of the hard (2-10 keV) emission is associated with the IGM. Relative to AGN flux, the IGM component peaks at a lower redshift (median z~0.45) so its clustering makes an important contribution to that of the total XRB. The angular correlations on 0.1-10 arcmin scales are significant, with an amplitude roughly consistent with an extrapolation of recent ROSAT results to small scales. A cross-correlation of the XRB against nearby galaxies taken from a simulated redshift survey also yields a strong signal from the IGM. Although some recent papers have argued that the expected soft band intensity from gas in galaxy, group, and cluster halos would exceed XRB limits unless much of the gas is expelled by supernova feedback, we obtain reasonable compatibility with current observations in a simulation that incorporates cooling, star formation, and only modest feedback. A prediction of our model is that the unresolved portion of the soft XRB will remain mostly unresolved.

The simultaneous presence of a strong quasi periodic oscillation of period of about 10 seconds in the optical and X-ray lightcurves of the X-ray transient XTE J1118+480 suggests that a significant fraction of the optical flux originates from the inner part of the accretion flow, where most of the X-rays are produced. We present a model of magnetic flares in an accretion disc corona where thermal cyclo-synchrotron emission contributes significantly to the optical emission, while the X-rays are produced by inverse Compton scattering of the soft photons produced by dissipation in the underlying disc and by the synchrotron process itself. Given the observational constraints, we estimate the values for the coronal temperature, optical depth and magnetic field intensity, as well as the accretion rate for the source. Within our model we predict a correlation between optical and hard X-ray variability and an anticorrelation between optical and soft X-rays. We also expect optical variability on flaring timescales (about tens of milliseconds), with a power density spectrum similar to the one observed in the X-ray band. Finally we use both the available optical/EUV/X-ray spectral energy distribution and the low frequency time variability to discuss limits on the inner radius of the optically thick disc.

The radio emission from supermassive black holes in nearby early-type galaxies can be used to test possible explanations for their low luminosities. We calculate the expected contribution from thermal synchrotron emission from hot accretion flows to the high radio frequency observations of NGC 2300, NGC 1399, NGC 4278 and NGC 4594. We find that, in all cases, and in accordance with our previous findings, hot flows accreting close to their Bondi rates overestimate significantly the observed fluxes. This implies that simply assuming a low radiative efficiency for the accreting gas is not enough to explain their low luminosities. Smaller central densities and accretion rates, as expected in the presence of strong mass loss or convection in the flows, can help reconcile the models with observations. We also show that a significant contribution to the high-frequency radio spectra can be due to non-thermal synchrotron emission from the small scale radio jets observed in these systems, allowing for even lower accretion rates in the inflows. We suggest that these outflows or jets may dump significant energy into the surronding medium close to the accretion radius and so reduce the accretion rates onto these systems. We discuss the relationship between the radio flux and black hole mass for the observed sample and its potential importance for probing accretion models.

Current estimates of number counts of radio sources in the frequency range where the most sensitive Cosmic Microwave Background (CMB) experiments are carried out significantly under-represent sources with strongly inverted spectra. Hot accretion flows around supermassive black holes in the nuclei of nearby galaxies are expected to produce inverted radio spectra by thermal synchrotron emission. We calculate the temperature fluctuations and power spectra of these sources in the Planck Surveyor 30 GHz energy channel, where their emission is expected to peak. We find that their potential contribution is generally comparable to the instrumental noise, and approaches the CMB anisotropy level at small angular scales. Forthcoming CMB missions, which will provide a large statistical sample of inverted-spectra sources, will be crucial for determining the distribution of hot accretion flows in nearby quiescent galactic nuclei. Detection of these sources in different frequency channels will help constrain their spectral characteristics, hence their physical properties.

Popular models for the formation of X-ray spectra in AGN assume that a large fraction of the disk's angular momentum dissipation takes place in a hot corona that carries a small amount of the accreting mass. Here I discuss the formation of a magnetically-structured accretion disk corona, generated by buoyancy instability in the disk and the heating of localized flare regions up to the canonical X-ray emitting temperatures. I also examine the analogy between accreting disk-coronae and ADAFs and discuss the relevant emission mechanism in these two accretion models and how observational constraints can allow us to discriminate between them.

We explore the implications of the discovery of hard, power-law X-ray sources in the spectra of nearby elliptical galaxies for the origin of the X-ray background. The spectra of these sources are consistent with models of thermal bremsstrahlung emission from low radiative efficiency accretion flows around central supermassive black holes and are unique in that they approximately match that of the hard XRB. If such sources, with luminosities consistent with those observed in nearby ellipticals, are present in most early-type galaxies, then their integrated emission may contribute significantly to the XRB. These sources may also contribute to the hard source counts detected in deep X-ray surveys.

We calculate upper bounds on the inner radii of geometrically thin accretion disks in galactic black-hole systems by relating their rapid variability properties to those of neutron stars. We infer that the inner disk radii do not exhibit large excursions between different spectral states, in contrast with the concept that the disk retreats significantly during the soft-to-hard state transition. We find that, in the hard state, the accretion disks extend down to radii less than 6-25 GM/c^2 and discuss the implications of our results for models of black-hole X-ray spectra.

We report the detection of hard X-ray emission components in the spectra of six nearby, giant elliptical galaxies observed with the ASCA satellite. The systems studied, which exhibit strong dynamical evidence for supermassive black holes in their nuclei, are M87, NGC 1399 and NGC 4696 (the dominant galaxies of the Virgo, Fornax and Centaurus clusters, respectively) and NGC 4472, 4636 and 4649 (three further giant ellipticals in the Virgo cluster). The ASCA data for all six sources provide clear evidence for hard emission components, which can be parameterized by power-law models with photon indices in the range \Gamma = 0.6-1.5 (mean value 1.2) and intrinsic 1-10 keV luminosities of 2 \times 10^{40} - 2 \times 10^{42} erg/s. Our results imply the identification of a new class of accreting X-ray source, with X-ray spectra significantly harder than those of binary X-ray sources, Seyfert nuclei or low luminosity AGN, and bolometric luminosities relatively dominated by their X-ray emission. We discuss various possible origins for the hard X-ray emission and argue that it is most likely to be due to accretion onto the central supermassive black holes, via low-radiative efficiency accretion coupled with strong outflows. In the case of M87, our detected power-law flux is in good agreement with a previously-reported measurement from ROSAT High Resolution Imager observations, which were able to resolve the jet from the nuclear X-ray emission components. We confirm previous results showing that the use of multiphase models in the analysis of the ASCA data leads to determinations of approximately solar emission-weighted metallicities for the X-ray gas in the galaxies. We also present results on the individual element abundances in NGC 4636.

M81 and NGC 4579 are two of the few low-luminosity active galactic nuclei which have an estimated mass for the central black hole, detected hard X-ray emission, and detected optical/UV emission. In contrast to the canonical ``big blue bump,'' both have optical/UV spectra which decrease with increasing frequency in a $\nu L_\nu$ plot. Barring significant reddening by dust and/or large errors in the black hole mass estimates, the optical/UV spectra of these systems require that the inner edge of a geometrically thin, optically thick, accretion disk lies at roughly 100 Schwarzschild radii. The observed X-ray radiation can be explained by an optically thin, two temperature, advection-dominated accretion flow at smaller radii.

The discovery of hard, power-law X-ray emission from a sample of six nearby elliptical galaxies, including the dominant galaxies of the Virgo,Fornax and Centaurus clusters (M87, NGC 1399 and NGC 4696, respectively), and NGC 4472, 4636 and 4649 in the Virgo cluster, has important implications for the study of quiescent supermassive black holes. We describe how the broad band spectral energy distributions of these galaxies, which accrete from their hot gaseous halos at rates comparable to their Bondi rates, can be explained by low-radiative efficiency accretion flows in which a significant fraction of the mass, angular momentum and energy is removed from the flows by winds. The observed suppression of the synchrotron component in the radio band and the systematically hard X-ray spectra, which are interpreted as thermal bremsstrahlung emission, support the conjecture that significant mass outflow is a natural consequence of systems accreting at low-radiative efficiencies. We briefly discuss an alternative model for the X-ray emission, namely that it is due to nonthermal synchrotron-self Compton processes in the accretion flow, or wind. This appears to require implausibly weak magnetic fields. Emission from a collimated jet viewed off axis should be distinguishable from the bremsstrahlung model by variability and thermal line emission studies. We argue that the difference in radiative efficiency between the nuclei of spiral and elliptical galaxies arises from the different manner in which interstellar gas is fed into the nuclei. In ellipticals, matter fed from the hot (slowly cooling) ISM is likely to be highly magnetized and with low specific angular momentum, both of which favor low-radiative efficiency accretion solutions and possibly the formation of the observed jets.

We discuss a series of observations of the black hole candidate GX 339-4 in low luminosity, spectrally hard states. We present spectral analysis of three separate archival Advanced Satellite for Cosmology and Astrophysics (ASCA) data sets and eight separate Rossi X-ray Timing Explorer (RXTE) data sets. Three of the RXTE observations were strictly simultaneous with 843 MHz and 8.3-9.1 GHz radio observations. All of these observations have (3-9 keV) flux approximately < 10^{-9} ergs s^{-1} cm^{-2}. The ASCA data show evidence for an 6.4 keV Fe line with equivalent width 40 eV, as well as evidence for a soft excess that is well-modeled by a power law plus a multicolor blackbody spectrum with peak temperature 150-200 eV. The RXTE data sets also show evidence of an Fe line with equivalent widths 20-140 eV. Reflection models show a hardening of the RXTE spectra with decreasing X-ray flux; however, these models do not exhibit evidence of a correlation between the photon index of the incident power law flux and the solid angle subtended by the reflector. `Sphere+disk' Comptonization models and Advection Dominated Accretion Flow (ADAF) models also provide reasonable descriptions of the RXTE data. The former models yield coronal temperatures in the range 20-50 keV and optical depths of \tau ~ 3. The model fits to the X-ray data, however, do not simultaneously explain the observed radio properties. The most likely source of the radio flux is synchrotron emission from an extended outflow of size greater than O(10^7 GM/c^2).

Several models for the hard X-ray Background (XRB) suggest that it is due to the emission from heavily obscured AGN. Recent studies have revealed the presence of a new population of hard X-ray sources which must contribute significantly to the XRB. However, whether the majority of these sources are obscured AGN or some other type of object at present remains unclear. Here, we further examine the possibility that a significant fraction of the XRB comes from a population of galaxies undergoing advection-dominated accretion in the high-\mdot regime and thus produce intrinsically hard spectra. When the accretion rate is close to its critical value, above which an advection-dominated accretion flow (ADAF) no longer exists, the major contribution to X-ray emission is due to inverse Compton scattering of the soft seed photons produced by cyclo-synchrotron radiation. In this regime, the resulting ADAF spectra are relatively hard with a fairly constant X-ray spectral index \alpha \sim 0.2-0.4 and a spectral cut-off at \sim 200 keV. We show that the integrated emission from such sources can provide a good fit to the hard (> 2 keV) X-ray background, provided that the spectrum is dominated by the contribution from objects located at redshifts z \sim 2-3. The model requires most of the contribution to the XRB to be due to objects accreting at the critical rate.

The growing evidence for supermassive black holes in the centres of nearby galaxies has brought into sharper focus the question of why elliptical galaxies, rich in hot gas, do not possess quasar-like luminosities. Recent studies suggest that advection-dominated accretion flows (ADAFs), with associated low radiative efficiency, might provide an explanation for the observed quiescence. Although ADAF models have been applied to a number of low-luminosity systems, compelling observational evidence for their existence is still required. Here, we examine new high-frequency radio observations of the giant, low-luminosity elliptical galaxies NGC 4649, NGC 4472 and NGC 4636 made with the Very Large Array (VLA) and the sub-millimetre common-user bolometer array (SCUBA) on the JCMT. At these frequencies the predictions are very precise and an ADAF is unequivocally characterised by a slowly rising spectrum with a sharp spectral cut-off produced by thermal synchrotron radiation. The new radio limits severely disagree with the canonical ADAF predictions which significantly overestimate the observed flux. While our observations do not rule out the presence of an ADAF in these objects, they place strong constraints on the model. If the accretion in these objects occurs in an advection-dominated mode then our radio limits imply that the emission from their central regions must be suppressed. We examine the possibility that the magnetic field in the flow is extremely low or that synchrotron emission is free-free absorbed by cold material in the accretion flow. We also discuss whether slow non-radiating accretion flows may drive winds/outflows to remove energy, angular momentum and mass so that the central densities, pressures and emissivities are much smaller than in a standard ADAF.

We present a model for the different X-ray spectral states displayed by Galactic Black Hole Candidates (GBHC). We discuss the physical and spectral implications for a magnetically structured corona in which magnetic flares result from reconnection of flux tubes rising from the accretion disk by the magnetic buoyancy instability. Using observations of one of the best studied examples, GX339-4, we identify the geometry and the physical conditions characterizing each of these states. We find that, in the Soft state, flaring occurs at small scale heights above the accretion disk. The soft thermal-like spectrum is the result of heating and consequent re-radiation of the hard X-rays produced by such flares. The hard tail is produced by Comptonization of the soft field radiation. Conversely, the hard state is the result of flares triggered high above the underlying accretion disk which produce X-rays via Comptonization of either internal synchrotron radiation or soft disk photons. The spectral characteristics of the different states are naturally accounted for by the choice of geometry: when flares are triggered high above the disk the system is photon-starved, hence the hard Comptonized spectrum of the hard state. Intense flaring close to the disk greatly enhances the soft-photon field with the result that the spectrum softens. We interpret the two states as being related to two different phases of magnetic energy dissipation. In the Soft state, Parker instability in the disk favours the emergence of large numbers of relatively low magnetic field flux tubes. In the hard state, only intense magnetic fields become buoyant. The model can also qualitatively account for the observed short timescale variability and the characteristics of the X-ray reflected component of the hard state.

A magnetically-structured accretion disk corona, generated by buoyancy instability in the disk, can account for observations of flare--like events in Active Galactic Nuclei. We examine how Petschek magnetic reconnection, associated with MHD turbulence, can result in a violent release of energy and heat the magnetically closed regions of the corona up to canonical X-ray emitting temperatures. X-ray magnetic flares, the after effect of the energy released in slow shocks, can account for the bulk of the X-ray luminosity from Seyfert galaxies and consistently explain the observed short-timescale variability.

We discuss the role of thermal cyclo-synchrotron emission in a magnetically dominated corona above an accretion disk in an active galaxy or a galactic black hole candidate. The dissipation occurs in localized active regions around the central black hole. Cyclo-synchrotron radiation is found to be an important process. In the case of galactic black hole candidates, its emission can dominate the inverse Compton scattering of the soft photon field produced by the disc by thermal electrons. We discuss observational predictions and the detectability of cyclo-synchrotron radiation both for these sources and radio-quiet active galactic nuclei.

We show that coronal magnetic dissipation in thin active sheets that sandwich standard thin accretion disks in active galactic nuclei may account for canonical electron temperatures of a few $\times 10^9$K if protons acquire most of the dissipated energy. Coulomb collisions transfer energy from the ions to the electrons, which subsequently cool rapidly by inverse-Compton scattering. In equilibrium, the proton energy density likely exceeds that of the magnetic field and both well exceed the electron and photon energy densities. The Coulomb energy transfer from protons to electrons is slow enough to maintain a high proton temperature, but fast enough to explain observed rapid X-ray variabilities in Seyferts. The $\sim 10^9$K electron temperature is insensitive to the proton temperature when the latter is $\ge 10^{12}$K.

It has been suggested that the final stages of accretion in present--day giant elliptical galaxies occur in an advection-dominated mode. The poor radiative efficiency of this accretion solution accounts for the fact that massive black holes, which are believed to reside in the centres of these galaxies, do not have the luminosities expected if accretion from the hot interstellar medium occurs at the Bondi rate. We discuss the advection-dominated solution for the nucleus of M60 (NGC 4649) and show that accretion at the Bondi rate is consistent with the core flux from radio to X-ray. This solution allows for a black hole mass of ~ 10^9 solar masses as required by independent arguments. The successful application of this model to both M60, and previously to the well known nucleus of M87, suggests that accretion of hot gas in an elliptical galaxy creates the ideal circumstances for advection-dominated accretion flows to operate.

It is believed that most giant elliptical galaxies possess nuclear black holes with masses in excess of $10^8\Msun$. Bondi accretion from the interstellar medium might then be expected to produce quasar-like luminosities from the nuclei of even quiescent elliptical galaxies. It is a puzzle that such luminosities are not observed. Motivated by this problem, Fabian & Rees have recently suggested that the final stages of accretion in these objects occurs in an advection-dominated mode with a correspondingly small radiative efficiency. Despite possessing a long-known active nucleus and dynamical evidence for a black hole, the low radiative and kinetic luminosities of the core of M87 provide the best illustration of this problem. We examine an advection-dominated model for the nucleus of M87 and show that accretion at the Bondi rate is compatible with the best known estimates for the core flux from radio through to X-ray wavelengths. The success of this model prompts us to propose that FR-I radio galaxies and quiescent elliptical galaxies accrete in an advection dominated mode whereas FR-II type radio-loud nuclei possess radiatively efficient thin accretion disks.

We discuss the possibility that a significant contribution of the hard X-ray Background is the integrated emission from a population of galaxies undergoing advection-dominated accretion in their nuclei. Owing to poor coupling between ions and electrons and to efficient radiative cooling of the electrons, the accreting plasma is two-temperature, with the ions being generally much hotter than the electrons and forming an ion-supported torus. We show that the electron te mperature then saturates at approximately 100keV independent of model parameters. At this temperature the hard X-ray emission is dominated by bremsstrahlung radiation. We find that this physical model gives an excellent fit to the spectrum of the XRB in the 3-60 keV range, provided that there is some evolution associated with the spectral emissivity which must peak at a redshift of about 2. We estimate that such galaxies contribute only to a small fraction of the local X-ray volume emissivity. The model implies a higher mean black hole mass than is obtained from the evolution of quasars alone.