We study outflows driven by Active Galactic Nuclei (AGNs) using high- resolution simulations of idealized z=2 isolated disk galaxies. Episodic accretion events lead to outflows with velocities >1000 km/s and mass outflow rates up to the star formatio
n rate (several tens of Msun/yr). Outflowing winds escape perpendicular to the disk with wide opening angles, and are typically asymmetric (i.e. unipolar) because dense gas above or below the AGN in the resolved disk inhibits outflow. Owing to rapid variability in the accretion rates, outflowing gas may be detectable even when the AGN is effectively off. The highest velocity outflows are sometimes, but not always, concentrated within 2-3 kpc of the galactic center during the peak accretion. With our purely thermal AGN feedback model -- standard in previous literature -- the outflowing material is mostly hot (10^6 K) and diffuse (nH<10^(-2) cm-3), but includes a cold component entrained in the hot wind. Despite the powerful bursts and high outflow rates, AGN feedback has little effect on the dense gas in the galaxy disk. Thus AGN-driven outflows in our simulations do not cause rapid quenching of star-formation, although they may remove significant amounts of gas of long (>Gyr) timescales.
We propose that star formation is delayed relative to the inflow rate in rapidly-accreting galaxies at very high redshift (z > 2) because of the energy conveyed by the accreting gas. Accreting gas streams provide fuel for star formation, but they sti
r the disk and increase turbulence above the usual levels compatible with gravitational instability, reducing the star formation efficiency in the available gas. After the specific inflow rate has sufficiently decreased - typically at z < 3 - galaxies settle in a self-regulated regime with efficient star formation. An analytic model shows that this interaction between infalling gas and young galaxies can significantly delay star formation and maintain high gas fractions (>40%) down to z = 2, in contrast to other galaxy formation models. Idealized hydrodynamic simulations of infalling gas streams onto primordial galaxies confirm the efficient energetic coupling at z > 2, and suggest that this effect is largely under-resolved in existing cosmological simulations.
Observations suggest that a large fraction of black hole growth occurs in normal star-forming disk galaxies. Here we describe simulations of black hole accretion in isolated disk galaxies with sufficient resolution (~5 pc) to track the formation of g
iant molecular clouds that feed the black hole. Black holes in z=2 gas-rich disks (fgas=50%) occasionally undergo ~10 Myr episodes of Eddington-limited accretion driven by stochastic collisions with massive, dense clouds. We predict that these gas-rich disks host weak AGNs 1/4 of the time, and moderate/strong AGNs 10% of the time. Averaged over 100 Myr timescales and the full distribution of accretion rates, the black holes grow at a few per cent of the Eddington limit -- sufficient to match observations and keep the galaxies on the MBH-Mbulge relation. This suggests that dense cloud accretion in isolated z=2 disks could dominate cosmic black hole growth. In z=0 disks with fgas=10%, Eddington-limited growth is extremely rare because typical gas clouds are smaller and more susceptible to disruption by AGN feedback. This results in an average black hole growth rate in high-fgas galaxies that is up to 1000 times higher than that in low-fgas galaxies. In all our simulations, accretion shows variability by factors of 10^4 on a variety of time scales, with variability at 1 Myr scales driven by the structure of the interstellar medium.
I highlight three results from cosmological hydrodynamic simulations that yield a realistic red sequence of galaxies: 1) Major galaxy mergers are not responsible for shutting off star-formation and forming the red sequence. Starvation in hot halos is
. 2) Massive galaxies grow substantially (about a factor of 2 in mass) after being quenched, primarily via minor (1:5) mergers. 3) Hot halo quenching naturally explains why galaxies are red when they either (a) are massive or (b) live in dense environments.
We examine the cosmic growth of the red sequence in a cosmological hydrodynamic simulation that includes a heuristic prescription for quenching star formation that yields a realistic passive galaxy population today. In this prescription, halos domina
ted by hot gas are continually heated to prevent their coronae from fueling new star formation. Hot coronae primarily form in halos above sim10^12 Modot, so that galaxies with stellar masses sim10^10.5 Modot are the first to be quenched and move onto the red sequence at z > 2. The red sequence is concurrently populated at low masses by satellite galaxies in large halos that are starved of new fuel, resulting in a dip in passive galaxy number densities around sim10^10 Modot. Stellar mass growth continues for galaxies even after joining the red sequence, primarily through minor mergers with a typical mass ratio sim1:5. For the most massive systems, the size growth implied by the distribution of merger mass ratios is typically sim2times the corresponding mass growth, consistent with observations. This model reproduces mass-density and colour-density trends in the local universe, with essentially no evolution to z = 1, with the hint that such relations may be washed out by z sim 2. Simulated galaxies are increasingly likely to be red at high masses or high local overdensities. In our model, the presence of surrounding hot gas drives the trends with both mass and environment.
Massive galaxies today typically are not forming stars despite being surrounded by hot gaseous halos with short central cooling times. This likely owes to some form of quenching feedback such as merger-driven quasar activity or radio jets emerging fr
om central black holes. Here we implement heuristic prescriptions for these phenomena on-the-fly within cosmological hydrodynamic simulations. We constrain them by comparing to observed luminosity functions and color-magnitude diagrams from SDSS. We find that quenching from mergers alone does not produce a realistic red sequence, because 1 - 2 Gyr after a merger the remnant accretes new fuel and star formation reignites. In contrast, quenching by continuously adding thermal energy to hot gaseous halos quantitatively matches the red galaxy luminosity function and produces a reasonable red sequence. Small discrepancies remain - a shallow red sequence slope suggests that our models underestimate metal production or retention in massive red galaxies, while a deficit of massive blue galaxies may reflect the fact that observed heating is intermittent rather than continuous. Overall, injection of energy into hot halo gas appears to be a necessary and sufficient condition to broadly produce red and dead massive galaxies as observed.
The bimodality in observed present-day galaxy colours has long been a challenge for hierarchical galaxy formation models, as it requires some physical process to quench (and keep quenched) star formation in massive galaxies. Here we examine phenomeno
logical models of quenching by post-processing the star formation histories of galaxies from cosmological hydrodynamic simulations that reproduce observations of star-forming galaxies reasonably well. We consider recipes for quenching based on major mergers, halo mass thresholds, gas temperature thresholds, and variants thereof. We compare the resulting simulated star formation histories to observed g-r colour-magnitude diagrams and red and blue luminosity functions from SDSS. The merger and halo mass quenching scenarios each yield a distinct red sequence and blue cloud of galaxies that are in broad agreement with data, albeit only under rather extreme assumptions. In detail, however, the simulated red sequence slope and amplitude in each scenario is somewhat discrepant, perhaps traceable to low metallicities in simulated galaxies. Merger quenching produces more massive blue galaxies, earlier quenching, and more frosting of young stars; comparing to relevant data tends to favor merger over halo mass quenching. Although physically-motivated quenching models can produce a red sequence, interesting generic discrepancies remain that indicate that additional physics is required to reproduce the star formation and enrichment histories of red and dead galaxies.
We post-process galaxy star formation histories in cosmological hydrodynamics simulations to test quenching mechanisms associated with AGN. By comparing simulation results to color-magnitude diagrams and luminosity functions of SDSS galaxies, we exam
ine whether quasar mode or radio mode AGN feedback can yield a realistic red sequence. Both cases yield red sequences distinct from the blue cloud, decent matches to the luminosity function, and galaxies that are too blue by about 0.1 magnitudes in g-r. Our merger-based prescription for quasar mode feedback, however, yields a red sequence build-up inconsistent with observations: the luminosity function lacks a characteristic knee, and the brightest galaxies include a small number of young stars.
We use HST/ACS images and a photometric catalog of the COSMOS field to analyze morphologies of the host galaxies of approximately 400 AGN candidates at redshifts 0.3 < z < 1.0. We compare the AGN hosts with a sample of non-active galaxies drawn from
the COSMOS field to match the magnitude and redshift distribution of the AGN hosts. We perform 2-D surface brightness modeling with GALFIT to yield host galaxy and nuclear point source magnitudes. X-ray selected AGN host galaxy morphologies span a substantial range that peaks between those of early-type, bulge-dominated and late-type, disk-dominated systems. We also measure the asymmetry and concentration of the host galaxies. Unaccounted for, the nuclear point source can significantly bias results of these measured structural parameters, so we subtract the best-fit point source component to obtain images of the underlying host galaxies. Our concentration measurements reinforce the findings of our 2-D morphology fits, placing X-ray AGN hosts between early- and late-type inactive galaxies. AGN host asymmetry distributions are consistent with those of control galaxies. Combined with a lack of excess companion galaxies around AGN, the asymmetry distributions indicate that strong interactions are no more prevalent among AGN than normal galaxies. In light of recent work, these results suggest that the host galaxies of AGN at these X-ray luminosities may be in a transition from disk-dominated to bulge-dominated, but that this transition is not typically triggered by major mergers.
