Philip F. Hopkins

A Unified, Merger-driven Model of the Origin of Starbursts, Quasars, the Cosmic X-Ray Background, Supermassive Black Holes, and Galaxy Spheroids

We present an evolutionary model for starbursts, quasars, and spheroidal galaxies in which mergers between gas-rich galaxies drive nuclear inflows of gas, producing 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 quasar luminosity, and we determine this dependence using a large set of galaxy merger simulations varying galaxy properties, orbital geometry, and gas physics. We use these fits to deconvolve observed quasar luminosity functions and obtain the evolution of the formation rate of quasars with peak luminosity, (Lpeak, z). Quasars spend extended periods at luminosities well below peak, so (Lpeak) has a maximum corresponding to the break in the observed luminosity function. From (Lpeak) and our simulations, we obtain self-consistent hard and soft X-ray and optical luminosity functions and predict many observables at multiple redshifts, including column density distributions of optical and X-ray samples, the luminosity function of broad-line quasars in X-ray samples and broad-line fraction versus luminosity, active BH mass functions, the distribution of Eddington ratios, the mass function of relic BHs and total BH mass density, and the cosmic X-ray background. In every case, our predictions agree well with observed estimates, without invoking ad hoc assumptions about source properties or distributions. We provide a library of Monte Carlo realizations of our models for comparison with observations.

An Observational Determination of the Bolometric Quasar Luminosity Function

We combine a large set of quasar luminosity function (QLF) measurements from the rest-frame optical, soft and hard X-ray, and near- and mid-IR bands to determine the bolometric QLF in the redshift interval z = 0-6. Accounting for the observed distributions of quasar column densities and variation of SED shapes, as well as their dependence on luminosity, makes it possible to integrate the observations in a reliable manner and provides a baseline in redshift and luminosity larger than that of any individual survey. We infer the QLF break luminosity and faint-end slope out to z ~ 4.5 and confirm at high significance (10 σ) previous claims of a flattening in both the faint- and bright-end slopes with redshift. With the best-fit estimates of the column density distribution and quasar SED, which both depend on luminosity, a single bolometric QLF self-consistently reproduces the observed QLFs in all bands and at all redshifts for which we compile measurements. Ignoring this luminosity dependence does not yield a self-consistent bolometric QLF and there is no evidence for any additional dependence on redshift. We calculate the expected relic black hole mass function and mass density, cosmic X-ray background, and ionization rate as a function of redshift and find that they are consistent with existing measurements. The peak in the total quasar luminosity density is well constrained at z = 2.15 ± 0.05. We provide a number of fitting functions to the bolometric QLF and its manifestations in various bands, as well as a script to return the QLF at arbitrary frequency and redshift from these fits.

A semi-analytic model for the co-evolution of galaxies, black holes and active galactic nuclei

We present a new semi-analytic model that self-consistently traces the growth of supermassive black holes (BH) and their host galaxies within the context of the Lambda cold dark matter (� CDM) cosmological framework. In our model, the energy emitted by accreting black holes regulates the growth of the black holes themselves, drives galactic scale winds that can remove cold gas from galaxies, and produces powerful jets that heat the hot gas atmospheres surrounding groups and clusters. We present a comprehensive comparison of our model predictions with observational measurements of key physical properties of low-redshift galaxies, such as cold gas fractions, stellar metallicities and ages, and specific star formation rates. We find that our new models successfully reproduce the exponential cut-off in the stellar mass function and the stellar and cold gas mass densities at z ∼ 0, and predict that star formation should be largely, but not entirely, quenched in massive galaxies at the present day. We also find that our model of self-regulated BH growth naturally reproduces the observed relation between BH mass and bulge mass. We explore the global formation history of galaxies and black holes in our models, presenting predictions for the cosmic histories of star formation, stellar mass assembly, cold gas and metals. We find that models assuming the ‘concordance’ � CDM cosmology overproduce star formation and stellar mass at high redshift (z 2). A model with less small-scale power predicts less star formation at high redshift, and excellent agreement with the observed stellar mass assembly history, but may have difficulty accounting for the cold gas in quasar absorption systems at high redshift (z ∼ 3–4).

A Cosmological Framework for the Co-Evolution of Quasars, Supermassive Black Holes, and Elliptical Galaxies. I. Galaxy Mergers and Quasar Activity

We develop a model for the cosmological role of mergers in the evolution of starbursts, quasars, and spheroidal galaxies. By combining theoretically well-constrained halo and subhalo mass functions as a function of redshift and environment with empirical halo occupation models, we can estimate where galaxies of given properties live at a particular epoch. This allows us to calculate, in an a priori cosmological manner, where major galaxy-galaxy mergers occur and what kinds of galaxies merge, at all redshifts. We compare this with the observed mass functions, clustering, fractions as a function of halo and galaxy mass, and small-scale environments of mergers, and we show that this approach yields robust estimates in good agreement with observations and can be extended to predict detailed properties of mergers. Making the simple Ansatz that major, gas-rich mergers cause quasar activity (but not strictly assuming they are the only triggering mechanism), we demonstrate that this model naturally reproduces the observed rise and fall of the quasar luminosity density at -->z = 0–6, as well as quasar luminosity functions, fractions, host galaxy colors, and clustering as a function of redshift and luminosity. The recent observed excess of quasar clustering on small scales at -->z ~ 0.2–2.5 is a natural prediction of our model, as mergers will preferentially occur in regions with excess small-scale galaxy overdensities. In fact, we demonstrate that quasar environments at all observed redshifts correspond closely to the empirically determined small group scale, where major mergers of ~L* gas-rich galaxies will be most efficient. We contrast this with a secular model in which quasar activity is driven by bars or other disk instabilities, and we show that, while these modes of fueling probably dominate the high Eddington ratio population at Seyfert luminosities (significant at -->z = 0), the constraints from quasar clustering, observed pseudobulge populations, and disk mass functions suggest that they are a small contributor to the -->z 1 quasar luminosity density, which is dominated by massive BHs in predominantly classical spheroids formed in mergers. Similarly, low-luminosity Seyferts do not show a clustering excess on small scales, in agreement with the natural prediction of secular models, but bright quasars at all redshifts do so. We also compare recent observations of the colors of quasar host galaxies and show that these correspond to the colors of recent merger remnants, in the transition region between the blue cloud and the red sequence, and are distinct from the colors of systems with observed bars or strong disk instabilities. Even the most extreme secular models, in which all bulge (and therefore BH) formation proceeds via disk instability, are forced to assume that this instability acts before the (dynamically inevitable) mergers, and therefore predict a history for the quasar luminosity density that is shifted to earlier times, in disagreement with observations. Our model provides a powerful means to predict the abundance and nature of mergers and to contrast cosmologically motivated predictions of merger products such as starbursts and active galactic nuclei.

Galaxies on FIRE (Feedback In Realistic Environments): Stellar feedback explains cosmologically inefficient star formation

We present a series of high-resolution cosmological simulations of galaxy formation to z = 0, spanning halo masses ∼10^8–10^(13) M⊙, and stellar masses ∼10^4–10^(11) M⊙. Our simulations include fully explicit treatment of the multiphase interstellar medium and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M_(halo) ∼ 10^(12) M⊙. We predict weak redshift evolution in the M*–M_(halo) relation, consistent with current constraints to z > 6. We find that the M*–M_(halo) relation is insensitive to numerical details, but is sensitive to feedback physics. Simulations with only supernova feedback fail to reproduce observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photoheating and radiation pressure) is necessary to destroy giant molecular clouds and enable efficient coupling of later supernovae to the gas. Star formation rates (SFRs) agree well with the observed Kennicutt relation at all redshifts. The galaxy-averaged Kennicutt relation is very different from the numerically imposed law for converting gas into stars, and is determined by self-regulation via stellar feedback. Feedback reduces SFRs and produces reservoirs of gas that lead to rising late-time star formation histories, significantly different from halo accretion histories. Feedback also produces large short-time-scale variability in galactic SFRs, especially in dwarfs. These properties are not captured by common ‘sub-grid’ wind models.

Black Holes in Galaxy Mergers: Evolution of Quasars

Basedonnumericalsimulationsofgas-richgalaxymergers,wediscussamodelinwhichquasaractivityistiedto the self-regulated growth of supermassive black holes in galaxies. The nuclear inflow of gas attending a galaxy collisiontriggersastarburstandfeedsblackholegrowth,butformostofthedurationofthestarburst,theblackhole is ‘‘buried,’’ being heavily obscured by surrounding gas and dust, limiting the visibility of the quasar, especially at optical and ultraviolet wavelengths. As the black hole grows, feedback energy from accretion heats the gas and eventuallyexpelsitinapowerfulwind, leaving behinda‘‘deadquasar.’’Betweenthe buried anddeadphases, there is a window in time 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 our estimated black hole growth timescales � 100 Myr, naturally producing a substantial population of buried quasars.However,theobservedandintrinsicenergyoutputsconvergeintheIRandhardX-raybandsasattenuation becomes weaker and chances of observation greatly increase. We also obtain the distribution of column densities along sight lines 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 the appropriate luminosity thresholds.Ourmodelreproducesawiderangeofquasarphenomena,includingobservedquasarlifetimes,intrinsic lifetimes, column density distributions, and differences between optical and X-ray samples, having properties consistent with observations across more than 5 orders of magnitude in bolometric luminosity from 10 9 to 10 14 L� (� 17PMB P� 30). Subject headingg cosmology: theory — galaxies: active — galaxies: evolution — galaxies: nuclei — quasars: general

HOW DO DISKS SURVIVE MERGERS

We develop a general physical model for how galactic disks survive and/or are destroyed in mergers and interactions. Based on simple dynamical arguments, we show that gas primarily loses angular momentum to internal torques in a merger, induced by the gravity of the secondary. Gas within some characteristic radius, determined by the efficiency of this angular momentum loss (itself a function of the orbital parameters, mass ratio, and gas fraction of the merging galaxies), will quickly lose angular momentum to the stars sharing the perturbed host disk, fall to the center, and be consumed in a starburst. We use a similar analysis to determine where violent relaxation of the premerger stellar disks is efficient on final coalescence. Our model describes both the dissipational and dissipationless components of the merger, and allows us to predict, for a given arbitrary encounter, the stellar and gas content of the material that will survive (without significant angular momentum loss or violent relaxation) to re-form a disk in the merger remnant, versus being dissipationlessly and violently relaxed or dissipationally losing angular momentum and forming a compact central starburst. We test these predictions with a large library of hydrodynamic merger simulations, and show that they agree well (with small scatter) with the properties of simulated merger remnants as a function of merger mass ratio, orbital parameters, and gas distributions in simulations which span a wide range of parameter space in these properties as well as prescriptions for gas physics, stellar and active galactic nucleus (AGN) feedback, halo and initial disk structural properties, redshift, and galaxy masses. We show that, in an immediate (short-term) sense, the amount of stellar or gaseous disk that survives or re-forms following a given interaction can be understood purely in terms of simple, well understood gravitational physics, independent of the details of the interstellar medium gas physics or stellar and AGN feedback. This allows us to demonstrate and quantify how these physics are in fact important, in an indirect sense, to enable disks to survive mergers by lowering star formation efficiencies in low-mass systems (allowing them to retain large gas fractions) and distributing the gas to large radii. The efficiency of disk destruction in mergers is a strong function of gas content—our model allows us to explicitly predict and demonstrate how, in sufficiently gas-rich mergers (with quite general orbital parameters), even 1:1 mass-ratio mergers can yield disk-dominated remnants, and more realistic 1:3-1:4 mass-ratio major mergers can yield systems with <20% of their mass in bulges. We discuss a number of implications of this modeling for the abundance and morphology of bulges as a function of mass and redshift, and provide simple prescriptions for the implementation of our results in analytic or semianalytic models of galaxy formation.

The Kinematic Structure of Merger Remnants

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 remnants, and all of the dissipationless ones, are slowly rotating and anisotropic. The remnants of gas-rich major mergers can well reproduce the observed distribution of projected ellipticities, rotation parameter (V/?)*, kinematic misalignments, ?, and isophotal shapes. The dissipationless remnants are a poor match to this data. We also investigate the properties of merger remnants as a function of initial disk gas fraction, orbital angular momentum, and the mass of the progenitor galaxies. 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.

The fundamental scaling relations of elliptical galaxies

We examine the fundamental scaling relations of elliptical galaxies formed through mergers. Using hundreds of simulations to judge the impact of progenitor galaxy properties on the properties of merger remnants, we find that gas dissipation provides an important contribution to tilt in the Fundamental Plane relation. Dissipationless mergers of disks produce remnants that occupy a plane similar to that delineated by the virial relation. As the gas content of progenitor disk galaxies is increased, the ti lt of the resulting Fundamental Plane relation increases and the slope of the Re - M⋆ relation steepens. For gas fractions fgas > 30%, the simulated Fundamental Plane scalings (Re ∝ σ 1.55 I -0.82 e ) approach those observed in the K-band (Re ∝ σ 1.53 I -0.79 e ). The dissipationless merging of spheroidal galaxies and the re-merging of disk galaxy remnants roughly maintain the tilt of the Fundamental Plane occupied by the progenitor ellipticals, approximately independent of the orbital energy or angular momentum. Dry merging of spheroidal systems at redshifts z < 1 is then expected to maintain the stellar-mass Fundamental Plane relations imprinted by gas-rich merging during the epoch of rapid spheroid and supermassive black hole growth at redshifts z ≈ 1 - 3. In our simulations, feedback from supermassive black hole growth has only a minor influence on the stellar-mass sca ling relations of spheroidal galaxies, but may play a role in maintaining the observed Fundamental Plane tilt at optical wavelengths by suppressing residual star formation in merger remnants. We estimate that ≈ 40 - 100% of the Fundamental Plane tilt induced by structural properties, as opposed to stellar population effects, owes to trends in the central total-to-stellar mass ratio Mtotal/M⋆ produced by dissipation. Gas cooling allows for an increase in central s tellar phase-space density relative to dissipationless mergers, thereby decreasing the central Mtotal/M⋆. Lower mass systems obtain greater phase-space densities than higher mass systems, producing a galaxy mass-dependent central Mtotal/M⋆ and a corresponding tilt in the Fundamental Plane. We account for these trends in the importance of dissipation with galaxy mass in terms of the inefficient cooling of collisionally heated gas in massi ve halos and dynamically varying gas consumption timescales in smaller systems.

Stellar feedback in galaxies and the origin of galaxy‐scale winds

Feedback from massive stars is believed to play a critical role in driving galactic superwinds that enrich the intergalactic medium and shape the galaxy mass function, massmetallicity relation, and other global galaxy properties. In previous papers, we have introduced new numerical methods for implementing stellar feedback on sub-GMC through galactic scales in numerical simulations of galaxies; the key physical processes include radiation pressure in the UV through IR, supernovae (Type-I & II), stellar winds (“fast” O star through “slow” AGB winds), and HII photoionization. Here, we show that these feedback mechanisms drive galactic winds with outflow rates as high as 10 20 times the galaxy star formation rate. The mass-loading efficiency (wind mass loss rate divided by the star formation rate) scales roughly as _ Mwind= _ M / V 1 c (where Vc is the galaxy circular velocity), consistent with simple momentum-conservation expectations. We use our suite of simulations to study the relative contribution of each feedback mechanism to the generation of galactic winds in a range of galaxy models, from SMC-like dwarfs and Milky-way analogues to z 2 clumpy disks. In massive, gas-rich systems (local starbursts and high-z galaxies), radiation pressure dominates the wind generation. By contrast, for MW-like spirals and dwarf galaxies the gas densities are much lower and sources of shock-heated gas such as supernovae and stellar winds dominate the production of large-scale outflows. In all of our models, however, the winds have a complex multi-phase structure that depends on the interaction between multiple feedback mechanisms operating on different spatial and time scales: any single feedback mechanism fails to reproduce the winds observed. We use our simulations to provide fitting functions to the wind mass-loading and velocities as a function of galaxy properties, for use in cosmological simulations and semi-analytic models. These differ from typically-adopted formulae with an explicit dependence on the gas surface density that can be very important in both low-density dwarf galaxies and high-density gas-rich galaxies.