Source author record

Philip F. Hopkins

Philip F. Hopkins appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.HE astro-ph.SR astro-ph.EP astro-ph physics.flu-dyn astro-ph.IM physics.plasm-ph hep-ph physics.comp-ph Machine Learning physics.ao-ph physics.space-ph

Catalog footprint

What is connected

140works

14topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

Early Stages of Dusty Tori: The First Infrared Spectra from a Highly Multiscale Quasar Simulation

We present the first infrared spectral predictions from a self-consistent simulation of the formation of a quasar in a starburst galaxy, spanning the cosmological environment to scales well below the dust sublimation region. The infrared (IR) emission is dominated by a torus-like dust structure composed of the highly magnetized, turbulence-supported outer accretion disk and of accreting gas tidally torn from the interstellar medium (ISM). At these early stages, the active galactic nuclei (AGN) is buried and Compton-thick. The near- to mid-IR escaping luminosity varies by almost an order of magnitude across sightlines, largely due to extinction from the inflowing stream of cold dust. Self-absorption within the torus suppresses silicate emission features, and further reprocessing by the ambient ISM leads to prominent silicate absorption and colder IR emission. The sublimation structure is stratified by composition and size, producing sightline-dependent extinction curves that intrinsically vary in shape. However, after repeated scattering in the optically thick dusty medium, these curves emerge substantially grayed. We also demonstrate that bipolar outflows from the central black hole that carve biconical cavities and reveal the central engine in later stages can preserve IR anisotropy and silicate features. These results suggest that dusty starburst quasars can undergo a buried, IR-bright phase early in their evolution.

preprint2022arXiv

A Consistent Reduced-Speed-of-Light Formulation of Cosmic Ray Transport Valid in Weak and Strong-Scattering Regimes

We derive a consistent set of moments equations for CR-magnetohydrodynamics, assuming a gyrotropic distribution function (DF). Unlike previous efforts we derive a closure, akin to the M1 closure in radiation hydrodynamics (RHD), that is valid in both the nearly-isotropic-DF and/or strong-scattering regimes, and the arbitrarily-anisotropic DF or free-streaming regimes, as well as allowing for anisotropic scattering and transport/magnetic field structure. We present the appropriate two-moment closure and equations for various choices of evolved variables, including the CR phase space distribution function, number density, total energy, kinetic energy, and their fluxes or higher moments, and the appropriate coupling terms to the gas. We show that this naturally includes and generalizes a variety of terms including convection/fluid motion, anisotropic CR pressure, streaming, diffusion, gyro-resonant/streaming losses, and re-acceleration. We discuss how this extends previous treatments of CR transport including diffusion and moments methods and popular forms of the Fokker-Planck equation, as well as how this differs from the analogous M1-RHD equations. We also present two different methods for incorporating a reduced speed of light (RSOL) to reduce timestep limitations: in both we carefully address where the RSOL (versus true c) must appear for the correct behavior to be recovered in all interesting limits, and show how current implementations of CRs with a RSOL neglect some additional terms.

preprint2022arXiv

Amplified J-factors in the Galactic Center for velocity-dependent darkmatter annihilation in FIRE simulations

We use FIRE-2 zoom cosmological simulations of Milky Way size galaxy halos to calculate astrophysical J-factors for dark matter annihilation and indirect detection studies. In addition to velocity-independent (s-wave) annihilation cross sections $σ_v$, we also calculate effective J-factors for velocity-dependent models, where the annihilation cross section is either either p-wave ($\propto v^2/c^2$) or d-wave ($\propto v^4/c^4$). We use 12 pairs of simulations, each run with dark-matter-only (DMO) physics and FIRE-2 physics. We observe FIRE runs produce central dark matter velocity dispersions that are systematically larger than in DMO runs by factors of $\sim 2.5-4$. They also have a larger range of central ($\sim 400$ pc) dark matter densities than the DMO runs ($ρ_{\rm FIRE}/ρ_{\rm DMO} \simeq 0.5 - 3$) owing to the competing effects of baryonic contraction and feedback. At 3 degrees from the Galactic Center, FIRE J-factors are $5-50$ (p-wave) and $15-500$ (d-wave) times higher than in the DMO runs. The change in s-wave signal at 3 degrees is more modest and can be higher or lower ($\sim 0.3-6$), though the shape of the emission profile is flatter (less peaked towards the Galactic Center) and more circular on the sky in FIRE runs. Our results for s-wave are broadly consistent with the range of assumptions in most indirect detection studies. We observe p-wave J-factors that are significantly enhanced compared to most past estimates. We find that thermal models with p-wave annihilation may be within range of detection in the near future.

preprint2022arXiv

Black hole -- galaxy scaling relations in FIRE: the importance of black hole location and mergers

The concurrent growth of supermassive black holes (SMBHs) and their host galaxies remains to be fully explored, especially at high redshift. While often understood as a consequence of self-regulation via AGN feedback, it can also be explained by alternative SMBH accretion models. Here, we expand on previous work by studying the growth of SMBHs with the help of a large suite of cosmological zoom-in simulations (MassiveFIRE) that are part of the Feedback in Realistic Environments (FIRE) project. The growth of SMBHs is modelled in post-processing with different black hole accretion models, placements, and merger treatments, and validated by comparing to on-the-fly calculations. Scaling relations predicted by the gravitational torque driven accretion (GTDA) model agree with observations at low redshift without the need for AGN feedback, in contrast to models in which the accretion rate depends strongly on SMBH mass. At high redshift, we find deviations from the local scaling relations in line with previous theoretical results. In particular, SMBHs are under-massive, presumably due to stellar feedback, but start to grow efficiently once their host galaxies reach $M_* \sim 10^{10} M_{\odot}$. We analyse and explain these findings in the context of a simple analytic model. Finally, we show that the predicted scaling relations depend sensitively on the SMBH location and the efficiency of SMBH merging, particularly in low-mass systems. These findings highlight the relevance of understanding the evolution of SMBH-galaxy scaling relations to predict the rate of gravitational wave signals from SMBH mergers across cosmic history.

preprint2022arXiv

Cluster assembly and the origin of mass segregation in the STARFORGE simulations

Stars form in dense, clustered environments, where feedback from newly formed stars eventually ejects the gas, terminating star formation and leaving behind one or more star clusters. Using the STARFORGE simulations, it is possible to simulate this process in its entirety within a molecular cloud, while explicitly evolving the gas radiation and magnetic fields and following the formation of individual, low-mass stars. We find that individual star-formation sites merge to form ever larger structures, while still accreting gas. Thus clusters are assembled through a series of mergers. During the cluster assembly process a small fraction of stars are ejected from their clusters; we find no significant difference between the mass distribution of the ejected stellar population and that of stars inside clusters. The star-formation sites that are the building blocks of clusters start out mass segregated with one or a few massive stars at their center. As they merge the newly formed clusters maintain this feature, causing them to have mass-segregated substructures without themselves being centrally condensed. The merged clusters relax to a centrally condensed mass segregated configuration through dynamical interactions between their members, but this process does not finish before feedback expels the remaining gas from the cluster. In the simulated runs the gas-free clusters then become unbound and break up. We find that turbulent driving and a periodic cloud geometry can significantly reduce clustering and prevent gas expulsion. Meanwhile, the initial surface density and level of turbulence have little qualitative effect on cluster evolution, despite the significantly different star formation histories.

preprint2022arXiv

Cosmic-Ray Driven Outflows to Mpc Scales from $L_{\ast}$ Galaxies

We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and inter-galactic medium (CGM/IGM), in high-resolution, fully-cosmological FIRE-2 simulations (accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive ($M_{\rm halo}\gtrsim 10^{11}\,M_{\odot}$), low-redshift ($z\lesssim 1-2$) halos can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive halos "traps" galactic outflows near the disk, so they recycle. With CRs injected in supernovae as modeled here, the low-pressure halo allows "escape" and CR pressure gradients continuously accelerate this material well into the IGM in "fast" outflows, while lower-density gas at large radii is accelerated in-situ into "slow" outflows that extend to $>$Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disk to $>$Mpc. The CR-driven outflows are primarily cool ($T\sim10^{5}\,$K) and low-velocity. All of these effects weaken and eventually vanish at lower halo masses ($\lesssim 10^{11}\,M_{\odot}$) or higher redshifts ($z\gtrsim 1-2$), reflecting the ratio of CR to thermal+gravitational pressure in the outer halo. We present a simple analytic model which explains all of the above phenomena.

preprint2022arXiv

Dissipative Dark Matter on FIRE: II. Observational signatures and constraints from local dwarf galaxies

We analyze the first set of cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modelled as a constant fractional energy dissipation ($f_{\rm diss}=0.5$). In this paper, we examine the properties of dwarf galaxies with $M_{\ast} \sim 10^{5}\operatorname{-}10^{9}\,{\rm M}_{\odot}$ in both isolation and within Milky Way-mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of stellar/neutral gas disk formation in dSIDM with $(σ/m)\leq 1\,{\rm cm^2\,g^{-1}}$ but they are still consistent with observed galaxy sizes and masses. In addition, as a result of the steeper central density profiles developed in dSIDM, the sub-kpc circular velocities of isolated dwarfs in models with $(σ/m)\geq 0.1\,{\rm cm^2\,g^{-1}}$ are enhanced by about a factor of two, which are still consistent with the measured stellar velocity dispersions of Local Group dwarfs but in tension with the HI rotation curves of more massive field dwarfs. Meanwhile, for satellites of the simulated Milky Way-mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help address the missing compact dwarf satellites in CDM. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with constant cross-section $(σ/m) \gtrsim 0.1\,{\rm cm^2\,g^{-1}}$ (assuming $f_{\rm diss}=0.5$) are effectively ruled out in bright dwarfs ($M_{\rm halo}\sim 10^{11}\,{\rm M}_{\odot}$) by circular velocity constraints. However, models with lower effective cross-sections (at this halo mass/velocity scale) are still viable and can give rise to non-trivial observable signatures.

preprint2022arXiv

Effects of Different Cosmic Ray Transport Models on Galaxy Formation

Cosmic rays (CRs) with ~GeV energies can contribute significantly to the energy and pressure budget in the interstellar, circumgalactic, and intergalactic medium (ISM, CGM, IGM). Recent cosmological simulations have begun to explore these effects, but almost all studies have been restricted to simplified models with constant CR diffusivity and/or streaming speeds. Physical models of CR propagation/scattering via extrinsic turbulence and self-excited waves predict transport coefficients which are complicated functions of local plasma properties. In a companion paper, we consider a wide range of observational constraints to identify proposed physically-motivated cosmic-ray propagation scalings which satisfy both detailed Milky Way (MW) and extra-galactic $γ$-ray constraints. Here, we compare the effects of these models relative to simpler 'diffusion+streaming' models on galaxy and CGM properties at dwarf through MW mass scales. The physical models predict large local variations in CR diffusivity, with median diffusivity increasing with galacto-centric radii and decreasing with galaxy mass and redshift. These effects lead to a more rapid dropoff of CR energy density in the CGM (compared to simpler models), in turn producing weaker effects of CRs on galaxy star formation rates (SFRs), CGM absorption profiles and galactic outflows. The predictions of the more physical CR models tend to lie 'in between' models which ignore CRs entirely and models which treat CRs with constant diffusivity.

preprint2022arXiv

Effects of the environment and feedback physics on the initial mass function of stars in the STARFORGE simulations

One of the key mysteries of star formation is the origin of the stellar initial mass function (IMF). The IMF is observed to be nearly universal in the Milky Way and its satellites, and significant variations are only inferred in extreme environments, such as the cores of massive elliptical galaxies. In this work we present simulations from the STARFORGE project that are the first cloud-scale RMHD simulations that follow individual stars and include all relevant physical processes. The simulations include detailed gas thermodynamics, as well as stellar feedback in the form of protostellar jets, stellar radiation, winds and supernovae. In this work we focus on how stellar radiation, winds and supernovae impact star-forming clouds. Radiative feedback plays a major role in quenching star formation and disrupting the cloud, however the IMF peak is predominantly set by protostellar jet physics. We find the effect of stellar winds is minor, and supernovae occur too late}to affect the IMF or quench star formation. We also investigate the effects of initial conditions on the IMF. The IMF is insensitive to the initial turbulence, cloud mass and cloud surface density, even though these parameters significantly shape the star formation history of the cloud, including the final star formation efficiency. The characteristic stellar mass depends weakly on metallicity and the interstellar radiation field. Finally, while turbulent driving and the level of magnetization strongly influences the star formation history, they only influence the high-mass slope of the IMF.

preprint2022arXiv

Evidence for a Vast Prograde Stellar Stream in the Solar Vicinity

Massive dwarf galaxies that merge with the Milky Way on prograde orbits can be dragged into the disk plane before being completely disrupted. Such mergers can contribute to an accreted stellar disk and a dark matter disk. We present evidence for Nyx, a vast new stellar stream in the vicinity of the Sun, that may provide the first indication that such an event occurred in the Milky Way. We identify about 500 stars that have coherent radial and prograde motion in this stream using a catalog of accreted stars built by applying deep learning methods to the second Gaia data release. Nyx is concentrated within $\pm 2$ kpc of the Galactic midplane and spans the full radial range studied (6.5-9.5 kpc). The kinematics of Nyx stars are distinct from those of both the thin and thick disk. In particular, its rotational speed lags the disk by $\sim 80$ km/s and its stars follow more eccentric orbits. A small number of Nyx stars have chemical abundances or inferred ages; from these, we deduce that Nyx stars have a peak metallicity of [Fe/H] $\sim -0.5$ and ages $\sim $10-13 Gyr. Taken together with the kinematic observations, these results strongly favor the interpretation that Nyx is the remnant of a disrupted dwarf galaxy. To further justify this interpretation, we explicitly demonstrate that metal-rich, prograde streams like Nyx can be found in the disk plane of Milky Way-like galaxies using the FIRE hydrodynamic simulations. Future spectroscopic studies will be able to validate whether Nyx stars originate from a single progenitor.

preprint2022arXiv

Exploring metallicity-dependent rates of Type Ia supernovae and their impact on galaxy formation

Type Ia supernovae play a critical role in stellar feedback and elemental enrichment in galaxies. Recent transient surveys like the All-Sky Automated Survey for Supernova (ASAS-SN) and the Dark Energy Survey (DES) find that the specific Ia rate at z ~ 0 may be ~ 15-50 times higher in lower-mass galaxies than at Milky Way-mass. Independently, Milky Way observations show that the close-binary fraction of solar-type stars is higher at lower metallicity. Motivated by these observations, we use the FIRE-2 cosmological zoom-in simulations to explore the impact of varying Ia rate models, including metallicity dependence, on galaxies across a range of stellar masses: 10^7 Msun - 10^{11} Msun. First, we benchmark our simulated star-formation histories (SFHs) against observations. We show that assumed SFHs and stellar mass functions play a major role in determining the degree of tension between observations and metallicity-independent Ia rate models, and potentially cause ASAS-SN and DES observations to be much more consistent with each other than might naively appear. Models in which the Ia rate increases with decreasing metallicity (as ~ Z^{-0.5} to Z^{-1}) provide significantly better agreement with observations. Encouragingly, these increases in Ia rate (> 10 times in low-mass galaxies) do not significantly impact galaxy stellar masses and morphologies: effective radii, axis ratios, and v/sigma remain largely unaffected except for our most extreme rate models. We explore implications for both [Fe/H] and [alpha/Fe] enrichment: metallicity-dependent Ia rate models can improve agreement with observed stellar mass-metallicity relations in low-mass galaxies. Our results demonstrate that a wide range of metallicity-dependent Ia models are viable for galaxy formation and motivate future work in this area.

preprint2022arXiv

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario.

preprint2022arXiv

Gas infall and radial transport in cosmological simulations of Milky Way-mass disks

Observations indicate that a continuous supply of gas is needed to maintain observed star formation rates in large, disky galaxies. To fuel star formation, gas must reach the inner regions of such galaxies. Despite its crucial importance for galaxy evolution, how and where gas joins galaxies is poorly constrained observationally and is rarely explored in fully cosmological simulations. To investigate gas accretion in the vicinity of galaxies, we analyze the FIRE-2 cosmological zoom-in simulations for 4 Milky Way mass galaxies (M_halo ~ 10E12 solar masses), focusing on simulations with cosmic ray physics. We find that at z~0, gas approaches the disk with angular momentum similar to the gaseous disk edge and low radial velocities, piling-up near the edge and settling into full rotational support. Accreting gas moves predominantly parallel to the disk with small but nonzero vertical velocity components, and joins the disk largely in the outskirts as opposed to "raining" down onto the disk. Once in the disk, gas trajectories are complex, being dominated by spiral arm induced oscillations and feedback. However, time and azimuthal averages show clear but slow net radial infall with transport speeds of 1-3 km/s and net mass fluxes through the disk on the order of one solar mass per year, comparable to the star formation rates of the galaxies and decreasing towards galactic center as gas is sunk into star formation. These rates are slightly higher in simulations without cosmic rays (1-7 km/s, ~4-5 solar masses per year). We find overall consistency of our results with observational constraints and discuss prospects of future observations of gas flows in and around galaxies.

preprint2022arXiv

Hot-mode accretion and the physics of thin-disk galaxy formation

We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (>~20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-like configuration. We show that the existence of this "rotating cooling flow" accretion mode is strongly correlated with the fraction of stars forming in a thin disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe.

preprint2022arXiv

Numerical Study of Cosmic Ray Confinement through Dust Resonant Drag Instabilities

We investigate the possibility of cosmic ray (CR) confinement by charged dust grains through resonant drag instabilities (RDIs). We perform magnetohydrodynamic particle-in-cell simulations of magnetized gas mixed with charged dust and cosmic rays, with the gyro-radii of dust and GeV CRs on $\sim\mathrm{AU}$ scales fully resolved. As a first study, we focus on one type of RDI wherein charged grains drift super-Alfv{é}nically, with Lorentz forces strongly dominating over drag forces. Dust grains are unstable to the RDIs and form concentrated columns and sheets, whose scale grows until saturating at the simulation box size. Initially perfectly-streaming CRs are strongly scattered by RDI-excited Alfv{é}n waves, with the growth rate of the CR perpendicular velocity components equaling the growth rate of magnetic field perturbations. These rates are well-predicted by analytic linear theory. CRs finally become isotropized and drift at least at $\sim v_\mathrm{A}$ by unidirectional Alfvén waves excited by the RDIs, with a uniform distribution of the pitch angle cosine $μ$ and a flat profile of the CR pitch angle diffusion coefficient $D_{μμ}$ around $μ= 0$, without the "$90$ degree pitch angle problem." With CR feedback on the gas included, $D_{μμ}$ decreases by a factor of a few, indicating a lower CR scattering rate, because the backreaction on the RDI from the CR pressure adds extra wave damping, leading to lower quasi-steady-state scattering rates. Our study demonstrates that the dust-induced CR confinement can be very important under certain conditions, e.g., the dusty circumgalactic medium around quasars or superluminous galaxies.

preprint2022arXiv

Rapid disc settling and the transition from bursty to steady star formation in Milky Way-mass galaxies

Recent observations and simulations indicate substantial evolution in the properties of galaxies with time, wherein rotationally-supported and steady thin discs (like those frequently observed in the local universe) emerge from galaxies that are clumpy, irregular, and have bursty star formation rates (SFRs). To better understand the progenitors of local disc galaxies we carry out an analysis of three FIRE-2 simulated galaxies with a mass similar to the Milky Way at redshift z=0. We show that all three galaxies transition from bursty to steady SFRs at a redshift between z=0.5 and z=0.8, and that this transition coincides with a rapid (< ~1 Gyr) emergence of a rotationally-supported interstellar medium (ISM).In the late phase with steady SFR, the rotational energy comprises > ~90% of the total kinetic + thermal energy in the ISM, and is roughly half the gravitational energy. By contrast, during the early phase with bursty star formation, the ISM has a quasi-spheroidal morphology and its energy budget is dominated by quasi-isotropic flows including turbulence and coherent inflows/outflows. This result, that rotational support is subdominant at early times, challenges the common application of equilibrium disc models to the high-redshift progenitors of Milky Way-like galaxies. We further find that the formation of a rotation-supported ISM coincides with the formation of a thermal energy-supported inner circumgalactic medium (CGM). Before this transition, the inner CGM is also supported by turbulence and coherent flows, indicating that at early times there is no clear boundary between the ISM and inner CGM.

preprint2022arXiv

Shapes of Milky-Way-Mass Galaxies with Self-Interacting Dark Matter

Self-interacting dark matter (SIDM) models offer one way to reconcile inconsistencies between observations and predictions from collisionless cold dark matter (CDM) models on dwarf-galaxy scales. In order to incorporate the effects of both baryonic and SIDM interactions, we study a suite of cosmological-baryonic simulations of Milky-Way (MW)-mass galaxies from the Feedback in Realistic Environments (FIRE-2) project where we vary the SIDM self-interaction cross-section $σ/m$. We compare the shape of the main dark matter (DM) halo at redshift $z=0$ predicted by SIDM simulations (at $σ/m=0.1$, $1$, and $10$ cm$^2$ g$^{-1}$) with CDM simulations using the same initial conditions. In the presence of baryonic feedback effects, we find that SIDM models do not produce the large differences in the inner structure of MW-mass galaxies predicted by SIDM-only models. However, we do find that the radius where the shape of the total mass distribution begins to differ from that of the stellar mass distribution is dependent on $σ/m$. This transition could potentially be used to set limits on the SIDM cross-section in the MW.

preprint2022arXiv

Simulating the Diverse Instabilities of Dust in Magnetized Gas

Recently Squire & Hopkins showed that charged dust grains moving through magnetized gas under the influence of any external force (e.g. radiation pressure, gravity) are subject to a spectrum of instabilities. Qualitatively distinct instability families are associated with different Alfvenic or magnetosonic waves and drift or gyro motion. We present a suite of simulations exploring these instabilities, for grains in a homogeneous medium subject to an external acceleration. We vary parameters such as the ratio of Lorentz-to-drag forces on dust, plasma $β$, size scale, and acceleration. All regimes studied drive turbulent motions and dust-to-gas fluctuations in the saturated state, can rapidly amplify magnetic fields into equipartition with velocity fluctuations, and produce instabilities that persist indefinitely (despite random grain motions). Different parameters produce diverse morphologies and qualitatively different features in dust, but the saturated gas state can be broadly characterized as anisotropic magnetosonic or Alfvenic turbulence. Quasi-linear theory can qualitatively predict the gas turbulent properties. Turbulence grows from small to large scales, and larger-scale modes usually drive more vigorous gas turbulence, but dust velocity and density fluctuations are more complicated. In many regimes, dust forms structures (clumps, filaments, sheets) that reach extreme over-densities (up to $\gg 10^{9}$ times mean), and exhibit substantial sub-structure even in nearly-incompressible gas. These can be even more prominent at lower dust-to-gas ratios. In other regimes, dust self-excites scattering via magnetic fluctuations that isotropize and amplify dust velocities, producing fast, diffusive dust motions.

preprint2022arXiv

Spatially Resolved Gas-phase Metallicity in FIRE-2 Dwarfs: Late-Time Evolution of Metallicity Relations in Simulations with Feedback and Mergers

We present an analysis of spatially resolved gas-phase metallicity relations in five dwarf galaxies ($M_{halo} \approx 10^{11} M_\odot$, $M_\star \approx 10^{8.8}-10^{9.6} M_\odot$) from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, which include an explicit model for sub-grid turbulent mixing of metals in gas, near $z\approx 0$, over a period of 1.4 Gyrs, and compare our findings with observations. While these dwarf galaxies represent a diverse sample, we find that all simulated galaxies match the observed mass-metallicity (MZR) and mass-metallicity gradient (MZGR) relations. We note that in all five galaxies, the metallicities are effectively identical between phases of the interstellar medium (ISM), with 95$\%$ being within $\pm$0.1 dex between various ISM phases, including the cold and dense gas ($T < 500$ K and $n_{\rm H} > 1$ cm$^{-3}$), ionized gas (near the H$α$ $T \approx 10^4$ K ridge-line), and nebular regions (ionized gas where the 10 Myr-averaged star formation rate is non-zero). We find that most of the scatter in relative metallicity between cold and dense gas and ionized gas/nebular regions can be attributed to either local starburst events or metal-poor inflows. We also note the presence of a major merger in one of our galaxies, m11e, with a substantial impact on the metallicity distribution in the spatially resolved map, showing two strong metallicity peaks and triggering a starburst in the main galaxy.

preprint2022arXiv

Testing Physical Models for Cosmic Ray Transport Coefficients on Galactic Scales: Self-Confinement and Extrinsic Turbulence at GeV Energies

The microphysics of ~GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant-diffusivity). We explore different physically-motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological FIRE simulations of dwarf and ~$L_{\ast}$ galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfven-wave damping, etc. We self-consistently predict observables including $γ$-rays ($L_γ$), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multi-phase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or $L_γ$), as transport is rate-limited by the ionized WIM and 'inner CGM' gaseous halo ($10^{4}-10^{6}$ K gas within 10-30 kpc), and $L_γ$ can be dominated by trapping in small 'patches.' Most physical ET models contribute negligible scattering of ~1-10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely-assumed scalings for SC models produce excessive confinement by factors >100 in the WIM and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.

preprint2022arXiv

The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert

We analyze the first giant molecular cloud (GMC) simulation to follow the formation of individual stars and their feedback from jets, radiation, winds, and supernovae, using the STARFORGE framework in the GIZMO code. We evolve the GMC for $\sim 9 \rm Myr$, from initial turbulent collapse to dispersal by feedback. Protostellar jets dominate feedback momentum initially, but radiation and winds cause cloud disruption at $\sim 8\%$ star formation efficiency (SFE), and the first supernova at $8.3 \rm Myr$ comes too late to influence star formation significantly. The per-freefall SFE is dynamic, accelerating from 0 to $\sim 18\%$ before dropping quickly to <1%, but the estimate from YSO counts compresses it to a narrower range. The primary cluster forms hierarchically and condenses to a brief ($\sim 1\,\mathrm{Myr}$) compact ($\sim 1 \rm pc$) phase, but does not virialize before the cloud disperses, and the stars end as an unbound expanding association. The initial mass function resembles the Chabrier (2005) form with a high-mass slope $α=-2$ and a maximum mass of $55 M_\odot$. Stellar accretion takes $\sim 400 \rm kyr$ on average, but $\gtrsim 1\rm Myr$ for $>10 M_\odot$ stars, so massive stars finish growing latest. The fraction of stars in multiples increases as a function of primary mass, as observed. Overall, the simulation much more closely resembles reality, compared to variations which neglect different feedback physics entirely. But more detailed comparison with synthetic observations is necessary to constrain the theoretical uncertainties.

preprint2022arXiv

The Galactic Dust-Up: Modeling Dust Evolution in FIRE

Recent strides have been made developing dust evolution models for galaxy formation simulations but these approaches vary in their assumptions and degree of complexity. Here we introduce and compare two separate dust evolution models (labelled 'Elemental' and 'Species'), based on recent approaches, incorporated into the GIZMO code and coupled with FIRE-2 stellar feedback and ISM physics. Both models account for turbulent dust diffusion, stellar production of dust, dust growth via gas-dust accretion, and dust destruction from time-resolved supernovae, thermal sputtering in hot gas, and astration. The "Elemental" model tracks the evolution of generalized dust species and utilizes a simple, 'tunable' dust growth routine, while the "Species" model tracks the evolution of specific dust species with set chemical compositions and incorporates a physically motivated, two-phase dust growth routine. We test and compare these models in an idealized Milky Way-mass galaxy and find that while both produce reasonable galaxy-integrated dust-to-metals (D/Z) ratios and predict gas-dust accretion as the main dust growth mechanism, a chemically motivated model is needed to reproduce the observed scaling relation between individual element depletions and D/Z with column density and local gas density. We also find the inclusion of theoretical metallic iron and O-bearing dust species are needed in the case of specific dust species in order to match observations of O and Fe depletions, and the integration of a sub-resolution dense molecular gas/CO scheme is needed to both match observed C depletions and ensure carbonaceous dust is not overproduced in dense environments.

preprint2022arXiv

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of "non-disky" dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy's intrinsic shape -- which is a direct quantification of the distribution of the galaxy's stellar content -- to interrogate the feedback implementation in simulated galaxies.

preprint2022arXiv

The observability of galaxy merger signatures in nearby gas-rich spirals

Galaxy mergers are crucial to understanding galaxy evolution, therefore we must determine their observational signatures to select them from large IFU galaxy samples such as MUSE and SAMI. We employ 24 high-resolution idealised hydrodynamical galaxy merger simulations based on the "Feedback In Realistic Environment" (FIRE-2) model to determine the observability of mergers to various configurations and stages using synthetic images and velocity maps. Our mergers cover a range of orbital configurations at fixed 1:2.5 stellar mass ratio for two gas rich spirals at low redshift. Morphological and kinematic asymmetries are computed for synthetic images and velocity maps spanning each interaction. We divide the interaction sequence into three: (1) the pair phase; (2) the merging phase; and (3) the post-coalescence phase. We correctly identify mergers between first pericentre passage and 500 Myr after coalescence using kinematic asymmetry with 66% completeness, depending upon merger phase and the field-of-view of the observation. We detect fewer mergers in the pair phase (40%) and many more in the merging and post-coalescence phases (97%). We find that merger detectability decreases with field-of-view, except in retrograde mergers, where centrally concentrated asymmetric kinematic features enhances their detectability. Using a cut-off derived from a combination of photometric and kinematic asymmetry, we increase these detections to 89% overall, 79% in pairs, and close to 100% in the merging and post-coalescent phases. By using this combined asymmetry cut-off we mitigate some of the effects caused by smaller fields-of-view subtended by massively multiplexed integral field spectroscopy programmes.

preprint2022arXiv

Why do Black Holes Trace Bulges (& Central Surface Densities), Instead of Galaxies as a Whole?

Previous studies of fueling black holes (BHs) in galactic nuclei have argued (on scales ~0.01-1000pc) accretion is dynamical with inflow rates $\dot{M}\simη\,M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass $M_{\rm gas}$, dynamical time $t_{\rm dyn}$, and some $η$. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency $f_{\rm wind}=M_{\rm ejected}/M_{\rm total}$ scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim\bar{a}_{\rm grav}/\langle\dot{p}/m_{\ast}\rangle\sim Σ_{\rm eff}/Σ_{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\,M_{\rm tot}(<r)/r^{2}$ and $\langle\dot{p}/m_{\ast}\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $η\rightarrow η\,(1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting e.g. the slope and normalization of the $M-σ$ and $M-M_{\rm bulge}$ relation, $L_{\rm AGN}-$SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient ($f_{\rm wind}\sim1$) below total-mass surface density $M_{\rm tot}/π\,r^{2}<Σ_{\rm crit}\sim3\times10^{9}\,M_{\odot}\,{\rm kpc^{-2}}$ (where $Σ_{\rm crit}=\langle\dot{p}/m_{\ast}\rangle/(π\,G)$), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness $Σ_{\rm crit}$ (B-band $μ_{\rm B}^{\rm crit}\sim 19\,{\rm mag\,arcsec^{-2}}$). This naturally explains why BH masses preferentially reflect bulge properties or central surface-densities ($Σ_{1\,{\rm kpc}}$), not 'total' galaxy properties.

preprint2021arXiv

Fiery Cores: Bursty and Smooth Star Formation Distributions across Galaxy Centers in Cosmological Zoom-in Simulations

We present an analysis of the $R\lesssim 1.5$ kpc core regions of seven simulated Milky Way mass galaxies, from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, for a finely sampled period ($Δt = 2.2$ Myr) of 22 Myr at $z \approx 0$, and compare them with star formation rate (SFR) and gas surface density observations of the Milky Way's Central Molecular Zone (CMZ). Despite not being tuned to reproduce the detailed structure of the CMZ, we find that four of these galaxies are consistent with CMZ observations at some point during this 22 Myr period. The galaxies presented here are not homogeneous in their central structures, roughly dividing into two morphological classes; (a) several of the galaxies have very asymmetric gas and SFR distributions, with intense (compact) starbursts occurring over a period of roughly 10 Myr, and structures on highly eccentric orbits through the CMZ, whereas (b) others have smoother gas and SFR distributions, with only slowly varying SFRs over the period analyzed. In class (a) centers, the orbital motion of gas and star-forming complexes across small apertures ($R \lesssim 150$pc, analogously $|l|<1^\circ$ in the CMZ observations) contributes as much to tracers of star formation/dense gas appearing in those apertures, as the internal evolution of those structures does. These asymmetric/bursty galactic centers can simultaneously match CMZ gas and SFR observations, demonstrating that time-varying star formation can explain the CMZ's low star formation efficiency.

preprint2021arXiv

On the optical properties of Resonant Drag Instabilities: Variability of Asymptotic Giant Branch and R Coronae Borealis stars

In dusty cool-star outflow or ejection events around AGB or RCB-like stars, dust is accelerated by radiation from the star and coupled to the gas via collisional drag forces. But it has recently been shown that such dust-gas mixtures are unstable to a super-class of instabilities called the resonant drag instabilities (RDIs), which promote dust clustering. We therefore consider idealized simulations of the RDIs operating on a spectrum of dust grain sizes subject to radiative acceleration (allowing for different grain optical properties), coupled to the gas with a realistic drag law, including or excluding the effects of magnetic fields and charged grains, and calculate for the first time how the RDIs could contribute to observed variability. We show that the RDIs naturally produce significant variations ($\sim 10-20\%$ $1σ$-level) in the extinction, corresponding to $\sim 0.1-1\,$mag level in the stellar types above, on timescales of order months to a year. The fluctuations are surprisingly robust to finite source-size effects as they are dominated by large-scale modes, which also means their spatial structure could be resolved in some nearby systems. We also quantify how this produces variations in the line-of-sight grain size-distribution. All of these variations are similar to those observed, suggesting that the RDIs may play a key role driving observed variability in dust extinction within dusty outflow/ejection events around cool stars. We further propose that the measured variations in grain sizes could directly be used to identify the presence of the RDIs in close by systems with observations.

preprint2021arXiv

Probing Hot Gas Components of Circumgalactic Medium in Cosmological Simulations with the Thermal Sunyaev-Zel'dovich Effect

The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful tool with the potential for constraining directly the properties of the hot gas that dominates dark matter halos because it measures pressure and thus thermal energy density. Studying this hot component of the circumgalactic medium (CGM) is important because it is strongly impacted by star-formation and active galactic nuclei (AGN) activity in galaxies, participating in the feedback loop that regulates star and black hole mass growth in galaxies. We study the tSZ effect across a wide halo mass range using three cosmological hydrodynamical simulations: Illustris-TNG, EAGLE, and FIRE-2. Specifically, we present the scaling relation between tSZ signal and halo mass and radial profiles of gas density, temperature, and pressure for all three simulations. The analysis includes comparisons to Planck tSZ observations and to the thermal pressure profile inferred from the Atacama Cosmology Telescope (ACT) measurements. We compare these tSZ data to the simulations to interpret the measurements in terms of feedback and accretion processes in the CGM. We also identify as-yet unobserved potential signatures of these processes that may be visible in future measurements, which will have the capability of measuring tSZ signals to even lower masses. We also perform internal comparisons between runs with different physical assumptions. We conclude: (1) there is strong evidence for the impact of feedback at $R_{500}$ but that this impact decreases by $5R_{500}$, and (2) the thermodynamic profiles of the CGM are highly dependent on the implemented model, such as cosmic-ray or AGN feedback prescriptions.

preprint2021arXiv

STARFORGE: The effects of protostellar outflows on the IMF

The initial mass function (IMF) of stars is a key quantity affecting almost every field of astrophysics, yet it remains unclear what physical mechanisms determine it. We present the first runs of the STARFORGE project, using a new numerical framework to follow the formation of individual stars in giant molecular clouds (GMCs) using the GIZMO code. Our suite include runs with increasingly complex physics, starting with isothermal ideal magnetohydrodynamic (MHD) and then adding non-isothermal thermodynamics and protostellar outflows. We show that without protostellar outflows the resulting stellar masses are an order of magnitude too high, similar to the result in the base isothermal MHD run. Outflows disrupt the accretion flow around the protostar, allowing gas to fragment and additional stars to form, thereby lowering the mean stellar mass to a value similar to that observed. The effect of jets upon global cloud evolution is most pronounced for lower-mass GMCs and dense clumps, so while jets can disrupt low-mass clouds, they are unable to regulate star formation in massive GMCs, as they would turn an order unity fraction of the mass into stars before unbinding the cloud. Jets are also unable to stop the runaway accretion of massive stars, which could ultimately lead to the formation of stars with masses $\mathrm{>500\,M_\odot}$. Although we find that the mass scale set by jets is insensitive to most cloud parameters (i.e., surface density, virial parameter), it is strongly dependent on the momentum loading of the jets (which is poorly constrained by observations) as well the the temperature of the parent cloud, which predicts slightly larger IMF variations than observed. We conclude that protostellar jets play a vital role in setting the mass scale of stars, but additional physics are necessary to reproduce the observed IMF.

preprint2020arXiv

A Census of Sub-kiloparsec Resolution Metallicity Gradients in Star-forming Galaxies at Cosmic Noon from HST Slitless Spectroscopy

We present hitherto the largest sample of gas-phase metallicity radial gradients measured at sub-kiloparsec resolution in star-forming galaxies in the redshift range of $z\in[1.2, 2.3]$. These measurements are enabled by the synergy of slitless spectroscopy from the Hubble Space Telescope near-infrared channels and the lensing magnification from foreground galaxy clusters. Our sample consists of 76 galaxies with stellar mass ranging from 10$^7$ to 10$^{10}$ $M_\odot$, instantaneous star-formation rate in the range of [1, 100] $M_\odot$/yr, and global metallicity [$\frac{1}{12}$, 2] solar. At 2-$σ$ confidence level, 15/76 galaxies in our sample show negative radial gradients, whereas 7/76 show inverted gradients. Combining ours and all other metallicity gradients obtained at similar resolution currently available in the literature, we measure a negative mass dependence of $Δ\log({\rm O/H})/Δr~ [\mathrm{dex~kpc^{-1}}] = \left(-0.020\pm0.007\right) + \left(-0.016\pm0.008\right) \log(M_\ast/10^{9.4} M_\odot)$ with the intrinsic scatter being $σ=0.060\pm0.006$ over four orders of magnitude in stellar mass. Our result is consistent with strong feedback, not secular processes, being the primary governor of the chemo-structural evolution of star-forming galaxies during the disk mass assembly at cosmic noon. We also find that the intrinsic scatter of metallicity gradients increases with decreasing stellar mass and increasing specific star-formation rate. This increase in the intrinsic scatter is likely caused by the combined effect of cold-mode gas accretion and merger-induced starbursts, with the latter more predominant in the dwarf mass regime of $M_\ast\lesssim10^9 M_\odot$.

preprint2020arXiv

A dark matter profile to model diverse feedback-induced core sizes of $Λ$CDM haloes

We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cuspy. We characterize our results using a new analytic density profile that extends the standard Einasto form to allow for a pronounced constant-density core in the resolved innermost radius. With one additional core-radius parameter, $r_{c}$, this "core-Einasto" profile is able to characterize the shape and normalization of our feedback-impacted dark matter halos. In order to enable comparisons with observations, we provide fitting functions for $r_{c}$ and other profile parameters as a function of both $M_\star$ and $M_{\star}/M_{\rm halo}$. In agreement with similar studies done in the literature, we find that dark matter core formation is most efficient at the characteristic stellar-mass to halo-mass ratio $M_\star/M_{\rm halo} \simeq 5 \times 10^{-3}$, or $M_{\star} \sim 10^9 \, M_{\odot}$, with cores that are roughly the size of the galaxy half-light radius, $r_{c} \simeq 1-5$ kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with $M_{\star}/M_{\rm halo} < 5\times 10^{-4}$ or $M_\star \lesssim 10^6 \, M_{\odot}$. For Milky Way-size galaxies, baryonic contraction often makes halos significantly more concentrated and dense at the stellar half-light radius than dark matter only runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of $\simeq 0.5-2$ kpc in size. Recent evidence for a ${\sim} 2$ kpc core in the Milky Way's dark matter halo is consistent with this expectation.

preprint2020arXiv

A general-purpose timestep criterion for simulations with gravity

We describe a new adaptive timestep criterion for integrating gravitational motion, which uses the tidal tensor to estimate the local dynamical timescale and scales the timestep proportionally. This provides a better candidate for a truly general-purpose gravitational timestep criterion than the usual prescription derived from the gravitational acceleration, which does not respect the equivalence principle, breaks down when $\mathbf{a}=0$, and does not obey the same dimensional scaling as the true timescale of orbital motion. We implement the tidal timestep criterion in the simulation code GIZMO, and examine controlled tests of collisionless galaxy and star cluster models, as well as fully-dynamic galaxy merger and cosmological dark matter simulations. The tidal criterion estimates the dynamical time faithfully, and generally provides a more efficient timestepping scheme compared to an acceleration criterion. Specifically, the tidal criterion achieves order-of-magnitude smaller energy errors for the same number of force evaluations in potentials with inner profiles shallower than $ρ\propto r^{-1}$ (ie. where $\mathbf{a}\rightarrow 0$), such as star clusters and cored galaxies. For a given problem these advantages must be weighed against the additional overhead of computing the tidal tensor on-the-fly, but in many cases this overhead is small.

preprint2020arXiv

A profile in FIRE: resolving the radial distributions of satellite galaxies in the Local Group with simulations

While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances <~ 100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.

preprint2020arXiv

A Relationship Between Stellar Metallicity Gradients and Galaxy Age in Dwarf Galaxies

We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between $10^{5.5}$ and $10^{8.6} \msun$. Whilst gas-phase metallicty gradients are generally weak in our simulated galaxies, we find that stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is correlated with galaxy-wide median stellar age, such that galaxies with younger stellar populations have flatter gradients. Stellar metallicty gradients are set by two competing processes: (1) the steady "puffing" of old, metal-poor stars by feedback-driven potential fluctuations, and (2) the accretion of extended, metal-rich gas at late times, which fuels late-time metal-rich star formation. If recent star formation dominates, then extended, metal-rich star formation washes out pre-existing gradients from the "puffing" process. We use published results from ten Local Group dwarf galaxies to show that a similar relationship between age and stellar metallicity-gradient strength exists among real dwarfs. This suggests that observed stellar metallicity gradients may be driven largely by the baryon/feedback cycle rather than by external environmental effects.

preprint2020arXiv

A Stable Finite-Volume Method for Scalar-Field Dark Matter

We describe and test a family of new numerical methods to solve the Schrodinger equation in self-gravitating systems, e.g. Bose-Einstein condensates or 'fuzzy'/ultra-light scalar field dark matter. The methods are finite-volume Godunov schemes with stable, higher-order accurate gradient estimation, based on a generalization of recent mesh-free finite-mass Godunov methods. They couple easily to particle-based N-body gravity solvers (with or without other fluids, e.g. baryons), are numerically stable, and computationally efficient. Different sub-methods allow for manifest conservation of mass, momentum, and energy. We consider a variety of test problems and demonstrate that these can accurately recover solutions and remain stable even in noisy, poorly-resolved systems, with dramatically reduced noise compared to some other proposed implementations (though certain types of discontinuities remain challenging). This is non-trivial because the "quantum pressure" is neither isotropic nor positive-definite and depends on higher-order gradients of the density field. We implement and test the method in the code GIZMO.

preprint2020arXiv

But What About... Cosmic Rays, Magnetic Fields, Conduction, & Viscosity in Galaxy Formation

We present a suite of high-resolution cosmological simulations, using the FIRE-2 feedback physics together with explicit treatment of magnetic fields, anisotropic conduction and viscosity, and cosmic rays (CRs) injected by supernovae (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultra-faint dwarf ($M_{\ast}\sim 10^{4}\,M_{\odot}$, $M_{\rm halo}\sim 10^{9}\,M_{\odot}$) through Milky Way masses, systematically vary CR parameters (e.g. the diffusion coefficient $κ$ and streaming velocity), and study an ensemble of galaxy properties (masses, star formation histories, mass profiles, phase structure, morphologies). We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$pc) scales have small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast} \ll 10^{10}\,M_{\odot}$, $M_{\rm halo} \lesssim 10^{11}\,M_{\odot}$), or at high redshifts ($z\gtrsim 1-2$), for any physically-reasonable parameters. However at higher masses ($M_{\rm halo} \gtrsim 10^{11}\,M_{\odot}$) and $z\lesssim 1-2$, CRs can suppress star formation by factors $\sim 2-4$, given relatively high effective diffusion coefficients $κ\gtrsim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$. At lower $κ$, CRs take too long to escape dense star-forming gas and lose energy to hadronic collisions, producing negligible effects on galaxies and violating empirical constraints from $γ$-ray emission. But around $κ\sim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which supports dense, cool ($T\ll 10^{6}$ K) gas that would otherwise rain onto the galaxy. CR heating (from collisional and streaming losses) is never dominant.

preprint2020arXiv

Cataloging Accreted Stars within Gaia DR2 using Deep Learning

The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significantly reduces the number of classifiable stars. The approach advocated here is applicable to a much larger portion of Gaia DR2. A method known as "transfer learning" is shown to be effective through extensive testing on a set of mock Gaia catalogs that are based on the FIRE cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies. The machine is first trained on simulated data using only 5D kinematics as inputs and is then further trained on a cross-matched Gaia/RAVE data set, which improves sensitivity to properties of the real Milky Way. The result is a catalog that identifies around 767,000 accreted stars within Gaia DR2. This catalog can yield empirical insights into the merger history of the Milky Way and could be used to infer properties of the dark matter distribution.

preprint2020arXiv

Measuring dynamical masses from gas kinematics in simulated high-redshift galaxies

Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disk galaxies, in which gas rotation accurately traces the gravitational potential, the interstellar medium (ISM) of z>1 galaxies is typically more dynamic and exhibits elevated turbulence. If not properly modeled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyze the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a range of galaxy properties from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M_* > 10^11 M_sun at z=1). Selecting only snapshots where a disk is present, we calculate the rotational profile v_phi(r) of the cool (10^3.5 K < T < 10^4.5 K) gas and compare it to the circular velocity v_c=sqrt(GM/r). In the simulated galaxies, the gas rotation traces the circular velocity at intermediate radii, but the two quantities diverge significantly in the center and in the outer disk. Our simulations appear to over-predict observed rotational velocities in the centers of massive galaxies (likely from a lack of black hole feedback), so we focus on larger radii. Gradients in the turbulent pressure at these radii can provide additional radial support and bias dynamical mass measurements low by up to 40%. In both the interior and exterior, the gas' motion can be significantly non-circular due to e.g. bars, satellites, and inflows/outflows. We discuss the accuracy of commonly-used analytic models for pressure gradients (or "asymmetric drift") in the ISM of high-redshift galaxies.

preprint2020arXiv

Most Stars (and Planets?) Are Born in Intense Radiation Fields

Protostars and young stars are strongly spatially "clustered" or "correlated" within their natal giant molecular clouds (GMCs). We demonstrate that such clustering leads to the conclusion that the incident bolometric radiative flux upon a random young star/disc is enhanced (relative to volume-averaged fluxes) by a factor which increases with the total stellar mass of the complex. Because the Galactic cloud mass function is top-heavy, the typical star in our Galaxy experienced a much stronger radiative environment than those forming in well-observed nearby (but relatively small) clouds, exceeding fluxes in the Orion Nebular Cluster by factors of $\gtrsim$30. Heating of the circumstellar disc around a median young star is dominated by this external radiation beyond $\sim 50\,$AU. And if discs are not well-shielded by ambient dust, external UV irradiation can dominate over the host star down to sub-AU scales. Another consequence of stellar clustering is an extremely broad Galaxy-wide distribution of incident flux (spanning $>10$ decades), with half the Galactic star formation in a substantial "tail" towards even more intense background radiation. We also show that the strength of external irradiation is amplified super-linearly in high-density environments such as the Galactic centre, starbursts, or high-redshift galaxies.

preprint2020arXiv

No missing photons for reionization: moderate ionizing photon escape fractions from the FIRE-2 simulations

We present the escape fraction of hydrogen ionizing photons (f_esc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z>5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 halos in M_vir~10^8--10^{12} M_sun (M_star~10^4--10^{10} M_sun) at z=5--12. We find the sample average <f_esc> increases with halo mass for M_vir~10^8--10^{9.5} M_sun, becomes nearly constant for M_vir~10^{9.5}--10^{11} M_sun, and decreases at M_vir>10^{11} M_sun. Equivalently, <f_esc> increases with stellar mass up to M_star~10^8 M_sun and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate <f_esc>~0.2 for galaxies at M_star~10^8 M_sun. Nearly half of the escaped ionizing photons come from stars 1--3 Myr old and the rest from stars 3--10 Myr old. Binaries only have a modest effect, boosting <f_esc> by ~25--35% and the number of escaped photons by 60--80%. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of <f_esc> at the high-mass end is due to dust attenuation, while at the low-mass end, <f_esc> decreases owing to inefficient star formation (and hence feedback). At fixed mass, <f_esc> tends to increase with redshift. Our simulations produce sufficient ionizing photons for cosmic reionization.

preprint2020arXiv

Physical models of streaming instabilities in protoplanetary disks

We develop simple, physically motivated models for drag-induced dust-gas streaming instabilities, which are thought to be crucial for clumping grains to form planetesimals in protoplanetary disks. The models explain, based on the physics of gaseous epicyclic motion and dust-gas drag forces, the most important features of the streaming instability and its simple generalisation, the disk settling instability. Some of the key properties explained by our models include the sudden change in the growth rate of the streaming instability when the dust-to-gas-mass ratio surpasses one, the slow growth rate of the streaming instability compared to the settling instability for smaller grains, and the main physical processes underlying the growth of the most unstable modes in different regimes. As well as providing helpful simplified pictures for understanding the operation of an interesting and fundamental astrophysical fluid instability, our models may prove useful for analysing simulations and developing nonlinear theories of planetesimal growth in disks.

preprint2020arXiv

Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $κ\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $γ$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these "CR-dominated" halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the "low" and "mid" ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.

preprint2020arXiv

Radiative Stellar Feedback in Galaxy Formation: Methods and Physics

Radiative feedback (RFB) from stars plays a key role in galaxies, but remains poorly-understood. We explore this using high-resolution, multi-frequency radiation-hydrodynamics (RHD) simulations from the Feedback In Realistic Environments (FIRE) project. We study ultra-faint dwarf through Milky Way mass scales, including H+He photo-ionization; photo-electric, Lyman Werner, Compton, and dust heating; and single+multiple scattering radiation pressure (RP). We compare distinct numerical algorithms: ray-based LEBRON (exact when optically-thin) and moments-based M1 (exact when optically-thick). The most important RFB channels on galaxy scales are photo-ionization heating and single-scattering RP: in all galaxies, most ionizing/far-UV luminosity (~1/2 of lifetime-integrated bolometric) is absorbed. In dwarfs, the most important effect is photo-ionization heating from the UV background suppressing accretion. In MW-mass galaxies, meta-galactic backgrounds have negligible effects; but local photo-ionization and single-scattering RP contribute to regulating the galactic star formation efficiency and lowering central densities. Without some RFB (or other 'rapid' FB), resolved GMCs convert too-efficiently into stars, making galaxies dominated by hyper-dense, bound star clusters. This makes star formation more violent and 'bursty' when SNe explode in these hyper-clustered objects: thus, including RFB 'smoothes' SFHs. These conclusions are robust to RHD methods, but M1 produces somewhat stronger effects. Like in previous FIRE simulations, IR multiple-scattering is rare (negligible in dwarfs, ~10% of RP in massive galaxies): absorption occurs primarily in 'normal' GMCs with A_v~1.

preprint2020arXiv

Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.

preprint2020arXiv

Synthetic Gaia surveys from the FIRE cosmological simulations of Milky Way-mass galaxies

With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky-Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments (FIRE) project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 $M_{\odot}$) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at $z = 0$ for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.

preprint2020arXiv

The Bolometric Quasar Luminosity Function at z = 0-7

In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from $z=0$ to $z=7$. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft and hard X-ray in past decades. Our method follows Hopkins et al. 2007 with an updated quasar SED model and bolometric and extinction corrections. The new best-fit bolometric quasar luminosity function behaves qualitatively different from the Hopkins et al. 2007 model at high redshift. Compared with the old model, the number density normalization decreases towards higher redshift and the bright-end slope is steeper at $z\gtrsim 2$. Due to the paucity of measurements at the faint end, the faint end slope at $z\gtrsim 5$ is quite uncertain. We present two models, one featuring a progressively steeper faint-end slope at higher redshift and the other featuring a shallow faint-end slope at $z\gtrsim 5$. Further multi-band observations of the faint-end QLF are needed to distinguish between these models. The evolutionary pattern of the bolometric QLF can be interpreted as an early phase likely dominated by the hierarchical assembly of structures and a late phase likely dominated by the quenching of galaxies. We explore the implications of this model on the ionizing photon production by quasars, the CXB spectrum, the SMBH mass density and mass functions. The predicted hydrogen photoionization rate contributed by quasars is subdominant during the epoch of reionization and only becomes important at $z\lesssim 3$. The predicted CXB spectrum, cosmic SMBH mass density and SMBH mass function are generally consistent with existing observations.

preprint2020arXiv

The IRX-$β$ relation of high-redshift galaxies

The relation between infrared excess (IRX) and UV spectral slope ($β_{\rm UV}$) is an empirical probe of dust properties of galaxies. The shape, scatter, and redshift evolution of this relation are not well understood, however, leading to uncertainties in estimating the dust content and star formation rates (SFRs) of galaxies at high redshift. In this study, we explore the nature and properties of the IRX-$β_{\rm UV}$ relation with a sample of $z=2-6$ galaxies ($M_*\approx 10^9-10^{12}\,M_\odot$) extracted from high-resolution cosmological simulations (MassiveFIRE) of the Feedback in Realistic Environments (FIRE) project. The galaxies in our sample show an IRX-$β_{\rm UV}$ relation that is in good agreement with the observed relation in nearby galaxies. IRX is tightly coupled to the UV optical depth, and is mainly determined by the dust-to-star geometry instead of total dust mass, while $β_{\rm UV}$ is set both by stellar properties, UV optical depth, and the dust extinction law. Overall, much of the scatter in the IRX-$β_{\rm UV}$ relation of our sample is found to be driven by variations of the intrinsic UV spectral slope. We further assess how the IRX-$β_{\rm UV}$ relation depends on viewing direction, dust-to-metal ratio, birth-cloud structures, and the dust extinction law and we present a simple model that encapsulates most of the found dependencies. Consequently, we argue that the reported `deficit' of the infrared/sub-millimetre bright objects at $z>5$ does not necessarily imply a non-standard dust extinction law at those epochs.

preprint2020arXiv

The Robustness of Cosmological Hydrodynamic Simulation Predictions to Changes in Numerics and Cooling Physics

We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$α$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

preprint2020arXiv

The universal acceleration scale from stellar feedback

It has been established for decades that rotation curves deviate from the Newtonian gravity expectation given baryons alone below a characteristic acceleration scale $g_{\dagger}\sim 10^{-8}\,\rm{cm\,s^{-2}}$, a scale promoted to a new fundamental constant in MOND. In recent years, theoretical and observational studies have shown that the star formation efficiency (SFE) of dense gas scales with surface density, SFE $\sim Σ/Σ_{\rm crit}$ with $Σ_{\rm crit} \sim \langle\dot{p}/m_{\ast}\rangle/(π\,G)\sim 1000\,\rm{M_{\odot}\,pc^{-2}}$ (where $\langle \dot{p}/m_{\ast}\rangle$ is the momentum flux output by stellar feedback per unit stellar mass in a young stellar population). We argue that the SFE, more generally, should scale with the local gravitational acceleration, i.e. that SFE $\sim g_{\rm tot}g_\mathrm{crit} \equiv (G\,M_{\rm tot}/R^{2}) / \langle\dot{p}/m_{\ast}\rangle$, where $M_{\rm tot}$ is the total gravitating mass and $g_\mathrm{crit}=\langle\dot{p}/m_{\ast}\rangle = π\,G\,Σ_{\rm crit} \approx 10^{-8}\,\rm{cm\,s^{-2}} \approx g_{\dagger}$. Hence the observed $g_\dagger$ may correspond to the characteristic acceleration scale above which stellar feedback cannot prevent efficient star formation, and baryons will eventually come to dominate. We further show how this may give rise to the observed acceleration scaling $g_{\rm obs}\sim(g_{\rm baryon}\,g_{\dagger})^{1/2}$ (where $g_{\rm baryon}$ is the acceleration due to baryons alone) and flat rotation curves. The derived characteristic acceleration $g_{\dagger}$ can be expressed in terms of fundamental constants (gravitational constant, proton mass, and Thomson cross section): $g_{\dagger}\sim 0.1\,G\,m_{p}/σ_{\rm T}$.

preprint2019arXiv

High-redshift Galaxy Formation with Self-consistently Modeled Stars and Massive Black Holes: Stellar Feedback and Quasar Growth

As computational resolution of modern cosmological simulations reach ever so close to resolving individual star-forming clumps in a galaxy, a need for "resolution-appropriate" physics for a galaxy-scale simulation has never been greater. To this end, we introduce a self-consistent numerical framework that includes explicit treatments of feedback from star-forming molecular clouds (SFMCs) and massive black holes (MBHs). In addition to the thermal supernovae feedback from SFMC particles, photoionizing radiation from both SFMCs and MBHs is tracked through full 3-dimensional ray tracing. A mechanical feedback channel from MBHs is also considered. Using our framework, we perform a state-of-the-art cosmological simulation of a quasar-host galaxy at z~7.5 for ~25 Myrs with all relevant galactic components such as dark matter, gas, SFMCs, and an embedded MBH seed of ~> 1e6 Ms. We find that feedback from SFMCs and an accreting MBH suppresses runaway star formation locally in the galactic core region. Newly included radiation feedback from SFMCs, combined with feedback from the MBH, helps the MBH grow faster by retaining gas that eventually accretes on to the MBH. Our experiment demonstrates that previously undiscussed types of interplay between gas, SFMCs, and a MBH may hold important clues about the growth and feedback of quasars and their host galaxies in the high-redshift Universe.

preprint2019arXiv

On the Survival of Cool Clouds in the Circum-Galactic Medium

We explore the survival of cool clouds in multi-phase circum-galactic media. We revisit the "cloud crushing problem" in a large survey of simulations including radiative cooling, self-shielding, self-gravity, magnetic fields, and anisotropic Braginskii conduction and viscosity (with saturation). We explore a wide range of parameters including cloud size, velocity, ambient temperature and density, as well as a variety of magnetic field configurations and cloud turbulence. We find that realistic magnetic fields and turbulence have weaker effects on cloud survival; the most important physics is radiative cooling and conduction. Self-gravity and self-shielding are important for clouds which are initially Jeans-unstable, but largely irrelevant otherwise. Non-self-gravitating, realistically magnetized clouds separate into four regimes: (1) At low column densities, clouds evaporate rapidly via conduction. (2) A "failed pressure confinement" regime, where the ambient hot gas cools too rapidly to provide pressure confinement for the cloud. (3) An "infinitely long-lived" regime, in which the cloud lifetime becomes longer than the cooling time of gas swept up in the leading bow shock, so the cloud begins to accrete and grow. (4) A "classical cloud destruction" regime, where clouds are eventually destroyed by instabilities. In the final regime, the cloud lifetime can exceed the naive cloud-crushing time owing to conduction-induced compression. However, small and/or slow-moving clouds can also evaporate more rapidly than the cloud-crushing time. We develop simple analytic models that explain the simulated cloud destruction times in this regime.

preprint2019arXiv

Stars made in outflows may populate the stellar halo of the Milky Way

We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at $z=0$.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to $\sim 10 \%$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $\sim 40\%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell $\&$ Sandange, who suggested that halo stars formed in rapidly $\textit{infalling}$ gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in the Milky Way.

preprint2019arXiv

Swirls of FIRE: Spatially Resolved Gas Velocity Dispersions and Star Formation Rates in FIRE-2 Disk Environments

We study the spatially resolved (sub-kpc) gas velocity dispersion ($σ$)--star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way mass disk galaxies at late times. In agreement with observations, we find a relatively flat relationship, with $σ\approx 15-30$ km/s in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant $σ$. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher $σ$ at constant SFR. The limits of the $σ$-$Σ_{\rm SFR}$ relation correspond to the onset of strong outflows. We see evidence of "on-off" cycles of star formation in the simulations, corresponding to feedback injection timescales of 10-100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disk (Toomre's $Q =1$) model predictions, and the data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e., ${\rm 1\%}$ per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of $σ$, it appears to be strongly subdominant to stellar feedback in the simulated galaxy disks.

preprint2018arXiv

Cosmic Rays or Turbulence can Suppress Cooling Flows (Where Thermal Heating or Momentum Injection Fail)

The quenching `maintenance' and `cooling flow' problems are important from the Milky Way through massive cluster elliptical galaxies. Previous work has shown that some source of energy beyond that from stars and pure magnetohydrodynamic processes is required, perhaps from AGN, but even the qualitative form of this energetic input remains uncertain. Different scenarios include thermal `heating,' direct wind or momentum injection, cosmic ray heating or pressure support, or turbulent `stirring' of the intra-cluster medium (ICM). We investigate these in $10^{12}-10^{14}\,{\rm M}_{\odot}$ halos using high-resolution non-cosmological simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, including simplified toy energy-injection models, where we arbitrarily vary the strength, injection scale, and physical form of the energy. We explore which scenarios can quench without violating observational constraints on energetics or ICM gas. We show that turbulent stirring in the central $\sim100\,$kpc, or cosmic-ray injection, can both maintain a stable low-SFR halo for $>$Gyr timescales with modest energy input, by providing a non-thermal pressure which stably lowers the core density and cooling rates. In both cases, associated thermal-heating processes are negligible. Turbulent stirring preserves cool-core features while mixing condensed core gas into the hotter halo and is by far the most energy efficient model. Pure thermal heating or nuclear isotropic momentum injection require vastly larger energy, are less efficient in lower-mass halos, easily over-heat cores, and require fine-tuning to avoid driving unphysical temperature gradients or gas expulsion from the halo center.

preprint2016arXiv

Binary Stars Can Provide the "Missing Photons" Needed for Reionization

Empirical constraints on reionization require galactic ionizing photon escape fractions fesc>20%, but recent high-resolution radiation-hydrodynamic calculations have consistently found much lower values ~1-5%. While these models include strong stellar feedback and additional processes such as runaway stars, they almost exclusively consider stellar evolution models based on single (isolated) stars, despite the fact that most massive stars are in binaries. We re-visit these calculations, combining radiative transfer and high-resolution cosmological simulations from the Feedback in Realistic Environments (FIRE) project. For the first time, we use a stellar evolution model that includes a physically and observationally motivated treatment of binaries (the BPASS model). Binary mass transfer and mergers enhance the population of massive stars at late times (>3 Myr) after star formation, which in turn strongly enhances the late-time ionizing photon production (especially at low metallicities). These photons are produced after feedback from massive stars has carved escape channels in the ISM, and so efficiently leak out of galaxies. As a result, the time-averaged "effective" escape fraction (ratio of escaped ionizing photons to observed 1500 A photons) increases by factors 4-10, sufficient to explain reionization. While important uncertainties remain, we conclude that binary evolution may be critical for understanding the ionization of the Universe.

preprint2016arXiv

Breathing FIRE: How Stellar Feedback Drives Radial Migration, Rapid Size Fluctuations, and Population Gradients in Low-Mass Galaxies

We examine the effects of stellar feedback and bursty star formation on low-mass galaxies ($M_{\rm star}=2\times10^6-5\times10^{10}{\rm M_{\odot}}$) using the FIRE (Feedback in Realistic Environments) simulations. While previous studies emphasized the impact of feedback on dark matter profiles, we investigate the impact on the stellar component: kinematics, radial migration, size evolution, and population gradients. Feedback-driven outflows/inflows drive significant radial stellar migration over both short and long timescales via two processes: (1) outflowing/infalling gas can remain star-forming, producing young stars that migrate $\sim1{\rm\,kpc}$ within their first $100 {\rm\,Myr}$, and (2) gas outflows/inflows drive strong fluctuations in the global potential, transferring energy to all stars. These processes produce several dramatic effects. First, galaxies' effective radii can fluctuate by factors of $>2$ over $\sim200 {\rm\,Myr}$, and these rapid size fluctuations can account for much of the observed scatter in radius at fixed $M_{\rm star}.$ Second, the cumulative effects of many outflow/infall episodes steadily heat stellar orbits, causing old stars to migrate outward most strongly. This age-dependent radial migration mixes---and even inverts---intrinsic age and metallicity gradients. Thus, the galactic-archaeology approach of calculating radial star-formation histories from stellar populations at $z=0$ can be severely biased. These effects are strongest at $M_{\rm star}\approx10^{7-9.6}{\rm M_{\odot}}$, the same regime where feedback most efficiently cores galaxies. Thus, detailed measurements of stellar kinematics in low-mass galaxies can strongly constrain feedback models and test baryonic solutions to small-scale problems in $Λ$CDM.

preprint2016arXiv

Evolution of Intrinsic Scatter in the SFR-Stellar Mass Correlation at 0.5<z<3

We present estimates of intrinsic scatter in the Star Formation Rate (SFR) - Stellar Mass (M*) correlation in the redshift range 0.5 < z < 3.0 and in the mass range 10^7 < M* < 10^11 Msun. We utilize photometry in the Hubble Ultradeep Field (HUDF12), Ultraviolet Ultra Deep Field (UVUDF) campaigns and CANDELS/GOODS-S. We estimate SFR, M* from broadband Spectral Energy Distributions (SEDs) and the best available redshifts. The maximum depth of the HUDF photometry (F160W 29.9 AB, 5 sigma depth) probes the SFR-M* correlation down to M* ~ 10 ^7 Msun, a factor of 10-100X lower in M* than previous studies, and comparable to dwarf galaxies in the local universe. We find the slope of the SFR-M* relationship to be near unity at all redshifts and the normalization to decrease with cosmic time. We find a moderate increase in intrinsic scatter with cosmic time from 0.2 to 0.4 dex across the epoch of peak cosmic star formation. None of our redshift bins show a statistically significant increase in intrinsic scatter at low mass. However, it remains possible that intrinsic scatter increases at low mass on timescales shorter than ~ 100 Myr. Our results are consistent with a picture of gradual and self-similar assembly of galaxies across more than three orders of magnitude in stellar mass from as low as 10^7 Msun.

preprint2016arXiv

Giant clumps in the FIRE simulations: a case study of a massive high-redshift galaxy

The morphology of massive star-forming galaxies at high redshift is often dominated by giant clumps of mass ~10^8-10^9 Msun and size ~100-1000 pc. Previous studies have proposed that giant clumps might have an important role in the evolution of their host galaxy, particularly in building the central bulge. However, this depends on whether clumps live long enough to migrate from their original location in the disc or whether they get disrupted by their own stellar feedback before reaching the centre of the galaxy. We use cosmological hydrodynamical simulations from the FIRE (Feedback in Realistic Environments) project that implement explicit treatments of stellar feedback and ISM physics to study the properties of these clumps. We follow the evolution of giant clumps in a massive (stellar mass ~10^10.8 Msun at z=1), discy, gas-rich galaxy from redshift z>2 to z=1. Even though the clumpy phase of this galaxy lasts over a gigayear, individual gas clumps are short-lived, with mean lifetime of massive clumps of ~20 Myr. During that time, they turn between 0.1% and 20% of their gas into stars before being disrupted, similar to local GMCs. Clumps with M>10^7 Msun account for ~20% of the total star formation in the galaxy during the clumpy phase, producing ~10^10 Msun of stars. We do not find evidence for net inward migration of clumps within the galaxy. The number of giant clumps and their mass decrease at lower redshifts, following the decrease in the overall gas fraction and star-formation rate.

preprint2016arXiv

Gravitational torque-driven black hole growth and feedback in cosmological simulations

We investigate black hole-host galaxy scaling relations in cosmological simulations with a self-consistent black hole growth and feedback model. The sub-grid accretion model captures the key scalings governing angular momentum transport from galactic scales down to parsec scales, while our kinetic feedback implementation enables the injection of outflows with properties chosen to match observed nuclear outflows. We show that "quasar mode" feedback can have a large impact on the thermal properties of the intergalactic medium and the growth of galaxies and massive black holes for kinetic feedback efficiencies as low as 0.1% relative to the bolometric luminosity. Nonetheless, our simulations suggest that the black hole-host scaling relations are only weakly dependent on the effects of black hole feedback on galactic scales, owing to feedback suppressing the growth of galaxies and massive black holes by a similar amount. In contrast, the rate at which gravitational torques feed the central black hole relative to the host galaxy star formation rate governs the slope and normalization of the black hole-host correlations. Our results suggest that a common gas supply regulated by gravitational torques is the primary driver of the observed co-evolution of black holes and galaxies.

preprint2016arXiv

MUFASA: Galaxy Formation Simulations With Meshless Hydrodynamics

We present the MUFASA suite of cosmological hydrodynamic simulations, which employs the GIZMO meshless finite mass (MFM) code including H2-based star formation, nine-element chemical evolution, two-phase kinetic outflows following scalings from the Feedback in Realistic Environments zoom simulations, and evolving halo mass-based quenching. Our fiducial (50 Mpc/h)^3 volume is evolved to z=0 with a quarter billion particles, The predicted galaxy stellar mass functions (GSMF) reproduce observations from z=4-0 to <1.2sigma in cosmic variance, providing an unprecedented match to this key diagnostic. The cosmic star formation history and stellar mass growth show general agreement with data, with a strong archaeological downsizing trend such that dwarf galaxies form the majority of their stars after z~1. We run 25 Mpc/h and 12.5 Mpc/h volumes to z=2 with identical feedback prescriptions, the latter resolving all hydrogen-cooling halos, and the three runs display fair resolution convergence. The specific star formation rates broadly agree with data at z=0, but are underpredicted at z~2 by a factor of three, re-emphasizing a longstanding puzzle in galaxy evolution models. We compare runs using MFM and two flavours of Smoothed Particle Hydrodynamics, and show that the GSMF is sensitive to hydrodynamics methodology at the ~x2 level, which is sub-dominant to choices for parameterising feedback.

preprint2016arXiv

Reconciling dwarf galaxies with LCDM cosmology: Simulating a realistic population of satellites around a Milky Way-mass galaxy

Low-mass "dwarf" galaxies represent the most significant challenges to the cold dark matter (CDM) model of cosmological structure formation. Because these faint galaxies are (best) observed within the Local Group (LG) of the Milky Way (MW) and Andromeda (M31), understanding their formation in such an environment is critical. We present first results from the Latte Project: the Milky Way on FIRE (Feedback in Realistic Environments). This simulation models the formation of a MW-mass galaxy to z = 0 within LCDM cosmology, including dark matter, gas, and stars at unprecedented resolution: baryon particle mass of 7070 Msun with gas kernel/softening that adapts down to 1 pc (with a median of 25 - 60 pc at z = 0). Latte was simulated using the GIZMO code with a mesh-free method for accurate hydrodynamics and the FIRE-2 model for star formation and explicit feedback within a multi-phase interstellar medium. For the first time, Latte self-consistently resolves the spatial scales corresponding to half-light radii of dwarf galaxies that form around a MW-mass host down to Mstar > 10^5 Msun. Latte's population of dwarf galaxies agrees with the LG across a broad range of properties: (1) distributions of stellar masses and stellar velocity dispersions (dynamical masses), including their joint relation; (2) the mass-metallicity relation; and (3) a diverse range of star-formation histories, including their mass dependence. Thus, Latte produces a realistic population of dwarf galaxies at Mstar > 10^5 Msun that does not suffer from the "missing satellites" or "too big to fail" problems of small-scale structure formation. We conclude that baryonic physics can reconcile observed dwarf galaxies with standard LCDM cosmology.

preprint2016arXiv

Star Formation in a Turbulent Framework: From Giant Molecular Clouds to Protostars

Turbulence is thought to be a primary driving force behind the early stages of star formation. In this framework large, self gravitating, turbulent clouds fragment into smaller clouds which in turn fragment into even smaller ones. At the end of this cascade we find the clouds which collapse into protostars. Following this process is extremely challenging numerically due to the large dynamical range so in this paper we propose a semi analytic framework which is able to follow star formation from the largest, giant molecular cloud (GMC) scale, to the final protostellar size scale. Due to the simplicity of the framework it is ideal for theoretical experimentation to explore the principal processes behind different aspects of star formation, at the cost of strong assumptions. The basic version of the model discussed in this paper only contains turbulence, gravity and crude assumptions about feedback, nevertheless it can reproduce the observed core mass function (CMF) and provide the protostellar system mass function (PSMF), which shows a striking resemblance to the observed IMF. Furthermore we find that to produce a universal IMF protostellar feedback must be taken into account otherwise the PSMF peak shows a strong dependence on the background temperature.

preprint2016arXiv

Stellar and Quasar Feedback in Concert: Effects on AGN Accretion, Obscuration, and Outflows

We study the interaction of feedback from active galactic nuclei (AGN) and a multi-phase interstellar medium (ISM), in simulations including explicit stellar feedback, multi-phase cooling, accretion-disk winds, and Compton heating. We examine radii ~0.1-100 pc around a black hole (BH), where the accretion rate onto the BH is determined and where AGN-powered winds and radiation couple to the ISM. We conclude: (1) The BH accretion rate is determined by exchange of angular momentum between gas and stars in gravitational instabilities. This produces accretion rates ~0.03-1 Msun/yr, sufficient to power luminous AGN. (2) The gas disk in the galactic nucleus undergoes an initial burst of star formation followed by several Myrs where stellar feedback suppresses the star formation rate (SFR). (3) AGN winds injected at small radii with momentum fluxes ~L/c couple efficiently to the ISM and have dramatic effects on ISM properties within ~100 pc. AGN winds suppress the nuclear SFR by factors ~10-30 and BH accretion rate by factors ~3-30. They increase the outflow rate from the nucleus by factors ~10, consistent with observational evidence for galaxy-scale AGN-driven outflows. (4) With AGN feedback, the predicted column density distribution to the BH is consistent with observations. Absent AGN feedback, the BH is isotropically obscured and there are not enough optically-thin sightlines to explain Type-I AGN. A 'torus-like' geometry arises self-consistently as AGN feedback evacuates gas in polar regions.

preprint2016arXiv

The formation of massive, quiescent galaxies at cosmic noon

The cosmic noon (z~1.5-3) marked a period of vigorous star formation for most galaxies. However, about a third of the more massive galaxies at those times were quiescent in the sense that their observed stellar populations are inconsistent with rapid star formation. The reduced star formation activity is often attributed to gaseous outflows driven by feedback from supermassive black holes, but the impact of black hole feedback on galaxies in the young Universe is not yet definitively established. We analyze the origin of quiescent galaxies with the help of ultra-high resolution, cosmological simulations that include feedback from stars but do not model the uncertain consequences of black hole feedback. We show that dark matter halos with specific accretion rates below ~0.25-0.4 per Gyr preferentially host galaxies with reduced star formation rates and red broad-band colors. The fraction of such halos in large dark matter only simulations matches the observed fraction of massive quiescent galaxies (~10^10-10^11 Msun). This strongly suggests that halo accretion rate is the key parameter determining which massive galaxies at z~1.5-3 become quiescent. Empirical models that connect galaxy and halo evolution, such as halo occupation distribution or abundance matching models, assume a tight link between galaxy properties and the masses of their parent halos. These models will benefit from adding the specific accretion rate of halos as a second model parameter.

preprint2016arXiv

The impact of stellar feedback on hot gas in galaxy haloes: the Sunyaev-Zel'dovich effect and soft X-ray emission

The thermal Sunyaev-Zel'dovich (SZ) effect and soft X-ray emission are routinely observed around massive galaxies and in galaxy groups and clusters. We study these observational diagnostics of galaxy haloes for a suite of cosmological `zoom-in' simulations from the `Feedback In Realistic Environments' project, which spans a large range in halo mass 10^10-10^13 Msun). We explore the effect of stellar feedback on the hot gas observables. The properties of our simulated groups, such as baryon fractions, SZ flux, and X-ray luminosities (L_X), are broadly consistent with existing observations, even though feedback from active galactic nuclei is not included. We make predictions for future observations of lower-mass objects for both SZ and diffuse X-ray measurements, finding that they are not just scaled-down versions of massive galaxies, but more strongly affected by galactic winds driven by star formation. Low-mass haloes (<~10^11 Msun) retain a low fraction of their baryons, which results in a strong suppression of the SZ signal. Our simulations therefore predict a scaling with halo mass that is steeper than self-similar for haloes less massive than 10^13 Msun. For halo masses <~10^12 Msun, L_X is time-variable and correlated primarily with the star formation rate (SFR). For these objects, the diffuse X-ray emission is powered mostly by galactic winds and the gas dominating the X-ray emission is flowing out with radial velocities close to the halo's circular velocity. For halo masses >~10^13 Msun, on the other hand, L_X is much less variable and not correlated with the SFR, because the emission originates from the quasi-hydrostatic, virialized halo gas.

preprint2015arXiv

A Simple Phenomenological Model for Grain Clustering in Turbulence

We propose a simple model for density fluctuations of aerodynamic grains, embedded in a turbulent, gravitating gas disk. The model combines a calculation for the behavior of a group of grains encountering a single turbulent eddy, with a hierarchical approximation of the eddy statistics. This makes analytic predictions for a range of quantities including: distributions of grain densities, power spectra and correlation functions of fluctuations, and maximum grain densities reached. We predict how these scale as a function of grain drag time t_stop, spatial scale, grain-to-gas mass ratio, strength of turbulence (alpha), and detailed disk properties. We test these against numerical simulations with various turbulence-driving mechanisms. The simulations agree well with the predictions, spanning t_stop*Omega ~ 1e-4 - 10, alpha ~ 1e-10 - 1e-2, and grain-to-gas mass ratio ~0-3. Results from 'turbulent concentration' simulations and laboratory experiments are also predicted as a special case. Vortices on a wide range of scales disperse and concentrate grains hierarchically. For small grains this is most efficient in eddies with turnover time comparable to the stopping time, but fluctuations are also damped by local gas-grain drift. For large grains, shear and gravity lead to a much broader range of eddy scales driving fluctuations, with most power on the largest scales. The grain density distribution has a log-Poisson shape, with fluctuations for large grains up to factors >1000. We provide simple analytic expressions for the predictions, and discuss implications for planetesimal formation, grain growth, and the structure of turbulence.

preprint2015arXiv

Accurate, Meshless Methods for Magneto-Hydrodynamics

Recently, we developed a pair of meshless finite-volume Lagrangian methods for hydrodynamics: the 'meshless finite mass' (MFM) and 'meshless finite volume' (MFV) methods. These capture advantages of both smoothed-particle hydrodynamics (SPH) and adaptive mesh-refinement (AMR) schemes. Here, we extend these to include ideal magneto-hydrodynamics (MHD). The MHD equations are second-order consistent and conservative. We augment these with a divergence-cleaning scheme, which maintains div*B~0 to high accuracy. We implement these in the code GIZMO, together with a state-of-the-art implementation of SPH MHD. In every one of a large suite of test problems, the new methods are competitive with moving-mesh and AMR schemes using constrained transport (CT) to ensure div*B=0. They are able to correctly capture the growth and structure of the magneto-rotational instability (MRI), MHD turbulence, and the launching of magnetic jets, in some cases converging more rapidly than AMR codes. Compared to SPH, the MFM/MFV methods exhibit proper convergence at fixed neighbor number, sharper shock capturing, and dramatically reduced noise, div*B errors, and diffusion. Still, 'modern' SPH is able to handle most of our tests, at the cost of much larger kernels and 'by hand' adjustment of artificial diffusion parameters. Compared to AMR, the new meshless methods exhibit enhanced 'grid noise' but reduced advection errors and numerical diffusion, velocity-independent errors, and superior angular momentum conservation and coupling to N-body gravity solvers. As a result they converge more slowly on some problems (involving smooth, slowly-moving flows) but more rapidly on others (involving advection or rotation). In all cases, divergence-control beyond the popular Powell 8-wave approach is necessary, or else all methods we consider will systematically converge to unphysical solutions.

preprint2015arXiv

Do We Expect Most AGN to Live in Disks?

Recent observations have indicated that a large fraction of the low to intermediate luminosity AGN population lives in disk-dominated hosts, while the more luminous quasars live in bulge-dominated hosts, in conflict with some previous model predictions. We therefore build and compare a semi-empirical model for AGN fueling which accounts for both merger and non-merger 'triggering.' In particular, we show that the 'stochastic accretion' model - in which fueling in disk galaxies is essentially a random process arising whenever dense gas clouds reach the nucleus - provides a good match to the present observations at low/intermediate luminosities. However it falls short of the high-luminosity population. We combine this with models for major merger-induced AGN fueling, which lead to rarer but more luminous events, and predict the resulting abundance of disk-dominated and bulge-dominated AGN host galaxies as a function of luminosity and redshift. We compile and compare observational constraints from z~0-2. The models and observations generically show a transition from disk to bulge dominance in hosts near the Seyfert-quasar transition, at all redshifts. 'Stochastic' fueling dominates AGN by number (dominant at low luminosity), and dominates BH growth below the knee in the present-day BH mass function (<10^7 M_sun). However it accounts for just ~10% of BH mass growth at masses >10^8 M_sun. In total, fueling in disky hosts accounts for ~30% of the total AGN luminosity density/BH mass density. The combined model also accurately predicts the AGN luminosity function and clustering/bias as a function of luminosity and redshift; however, we argue that these are not sensitive probes of BH fueling mechanisms.

preprint2015arXiv

Forged in FIRE: cusps, cores, and baryons in low-mass dwarf galaxies

We present ultra-high resolution cosmological hydrodynamic simulations of $M_*\simeq10^{4-6.3}M_{\odot}$ dwarf galaxies that form within $M_{v}=10^{9.5-10}M_{\odot}$ dark matter halos. Our simulations rely on the FIRE implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve stellar and dark matter structure on the ~200 pc scales of interest for classical and ultra-faint dwarfs in the Local Group. The resultant galaxies sit on the $M_*$ vs. $M_{v}$ relation required to match the Local Group stellar mass function. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. We demonstrate that it is possible to create a large (~1 kpc) dark matter core in a cosmological simulation of an $M_*\simeq10^{6.5}M_{\odot}$ dwarf galaxy that resides within an $M_{v}=10^{10}M_{\odot}$ halo -- precisely the scale of interest for resolving the Too Big to Fail problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late ($z\lesssim2$), after the early epoch of cusp building mergers has ended. Our $M_*\simeq10^4M_{\odot}$ dwarf retains a cuspy dark matter halo density profile that matches almost identically that of a dark-matter only run of the same system. Despite forming in a field environment, this very low mass dwarf has observable properties that match closely to those of ultra-faint satellite galaxies of the Milky Way, including a uniformly old stellar population (>10 Gyr). Though ancient, most of the stars in our ultra-faint form after reionization; the UV field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the proto-dwarf.

preprint2015arXiv

GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods

We present two new Lagrangian methods for hydrodynamics, in a systematic comparison with moving-mesh, SPH, and stationary (non-moving) grid methods. The new methods are designed to simultaneously capture advantages of both smoothed-particle hydrodynamics (SPH) and grid-based/adaptive mesh refinement (AMR) schemes. They are based on a kernel discretization of the volume coupled to a high-order matrix gradient estimator and a Riemann solver acting over the volume 'overlap.' We implement and test a parallel, second-order version of the method with self-gravity & cosmological integration, in the code GIZMO: this maintains exact mass, energy and momentum conservation; exhibits superior angular momentum conservation compared to all other methods we study; does not require 'artificial diffusion' terms; and allows the fluid elements to move with the flow so resolution is automatically adaptive. We consider a large suite of test problems, and find that on all problems the new methods appear competitive with moving-mesh schemes, with some advantages (particularly in angular momentum conservation), at the cost of enhanced noise. The new methods have many advantages vs. SPH: proper convergence, good capturing of fluid-mixing instabilities, dramatically reduced 'particle noise' & numerical viscosity, more accurate sub-sonic flow evolution, & sharp shock-capturing. Advantages vs. non-moving meshes include: automatic adaptivity, dramatically reduced advection errors & numerical overmixing, velocity-independent errors, accurate coupling to gravity, good angular momentum conservation and elimination of 'grid alignment' effects. We can, for example, follow hundreds of orbits of gaseous disks, while AMR and SPH methods break down in a few orbits. However, fixed meshes minimize 'grid noise.' These differences are important for a range of astrophysical problems.

preprint2015arXiv

Gusty, gaseous flows of FIRE: galactic winds in cosmological simulations with explicit stellar feedback

We present an analysis of the galaxy-scale gaseous outflows from the FIRE (Feedback in Realistic Environments) simulations. This suite of hydrodynamic cosmological zoom simulations resolves formation of star-forming giant molecular clouds to $z=0$, and features an explicit stellar feedback model on small scales. Our simulations reveal that high redshift galaxies undergo bursts of star formation followed by powerful gusts of galactic outflows that eject much of the ISM and temporarily suppress star formation. At low redshift, however, sufficiently massive galaxies corresponding to L*-progenitors develop stable disks and switch into a continuous and quiescent mode of star formation that does not drive outflows far into the halo. Mass-loading factors for winds in L*-progenitors are $η\approx 10$ at high redshift, but decrease to $η\ll 1$ at low redshift. Although lower values of $η$ are expected as halos grow in mass over time, we show that the strong suppression of outflows with decreasing redshift cannot be explained by mass evolution alone. Circumgalactic outflow velocities are variable and broadly distributed, but typically range between one and three times the circular velocity of the halo. Much of the ejected material builds a reservoir of enriched gas within the circumgalactic medium, some of which could be later recycled to fuel further star formation. However, a fraction of the gas that leaves the virial radius through galactic winds is never regained, causing most halos with mass $M_h \le 10^{12} M_{\odot}$ to be deficient in baryons compared to the cosmic mean by $z=0$.

preprint2015arXiv

Jumping the Gap: The Formation Conditions and Mass Function of Pebble-Pile Planetesimals

In a turbulent proto-planetary disk, dust grains undergo large density fluctuations and under the right circumstances, these grain overdensities can overcome shear, turbulent, and gas pressure support to collapse under self-gravity (forming a "pebble pile" planetesimal). Using a simple analytic model for the fluctuations predicted in simulations, we estimate the rate-of-formation and mass function of self-gravitating, rapidly-collapsing planetesimal-mass bodies formed by this mechanism. The statistics of this process depend sensitively on the size/stopping time of the grains, disk surface density, and turbulent Mach numbers. However, when it occurs, we predict that the resulting planetesimal mass function is broad and quasi-universal, with a slope dN/dM~1/M^(1-2), spanning a size/mass range ~10-1e4 km (~1e-9-5.0 M_Earth). Collapse to planetesimal through super-Earth masses is possible. The key condition is that grain density fluctuations reach large amplitudes on large scales, where gravitational instability proceeds most easily (collapse of small grains is strongly suppressed by turbulent vorticity). We show this leads to a new criterion for 'pebble-pile' formation in terms of the dimensionless particle stopping time (tau_stop > f(Q,Z,alpha)). In a MMSN, this requires grains larger than a=(50,1,0.1)cm at r=(1,30,100)au. At small radii, it would depend on the existence of large boulders. However, because density fluctuations depend super-exponentially on tau_stop (inversely proportional to disk surface density), lower-density disks are more unstable. In fact, we predict that cm-sized grains at ~1au will form pebble piles in a disk with ~10% the MMSN density, so planet formation at ~au may generically occur late, as disks are evaporating. We also predict that conditions become progressively more favorable for pebble-pile formation around lower-mass, cooler stars.

preprint2015arXiv

Sweating the small stuff: simulating dwarf galaxies, ultra-faint dwarf galaxies, and their own tiny satellites

We present FIRE/Gizmo hydrodynamic zoom-in simulations of isolated dark matter halos, two each at the mass of classical dwarf galaxies ($M_{\rm vir} \simeq 10^{10} M_{\odot}$) and ultra-faint galaxies ($M_{\rm vir} \simeq 10^9 M_{\odot}$), and with two feedback implementations. The resultant central galaxies lie on an extrapolated abundance matching relation from $M_{\star} \simeq 10^6$ to $10^4 M_{\odot}$ without a break. Every host is filled with subhalos, many of which form stars. Our dwarfs with $M_{\star} \simeq 10^6 M_{\odot}$ each have 1-2 well-resolved satellites with $M_{\star} = 3-200 \times 10^3 M_{\odot}$. Even our isolated ultra-faint galaxies have star-forming subhalos. If this is representative, dwarf galaxies throughout the universe should commonly host tiny satellite galaxies of their own. We combine our results with the ELVIS simulations to show that targeting $\sim 50~ \rm kpc$ regions around nearby isolated dwarfs could increase the chances of discovering ultra-faint galaxies by $\sim 35\%$ compared to random halo pointings, and specifically identify the region around the Phoenix dwarf galaxy as a good potential target. The well-resolved ultra-faint galaxies in our simulations ($M_{\star} \simeq 3 - 30 \times 10^3 M_{\odot}$) form within $M_{\rm peak} \simeq 0.5 - 3 \times 10^9 M_{\odot}$ halos. Each has a uniformly ancient stellar population ($ > 10~ \rm Gyr$) owing to reionization-related quenching. More massive systems, in contrast, all have late-time star formation. Our results suggest that $M_{\rm halo} \simeq 5 \times 10^9 M_{\odot}$ is a probable dividing line between halos hosting reionization "fossils" and those hosting dwarfs that can continue to form stars in isolation after reionization.

preprint2015arXiv

The creation and persistence of a misaligned gas disc in a simulated early-type galaxy

Massive early-type galaxies commonly have gas discs which are kinematically misaligned with the stellar component. These discs feel a torque from the stars and the angular momentum vectors are expected to align quickly. We present results on the evolution of a misaligned gas disc in a cosmological simulation of a massive early-type galaxy from the Feedback In Realistic Environments project. This galaxy experiences a merger which, together with a strong galactic wind, removes most of the original gas disc. The galaxy subsequently reforms a gas disc through accretion of cold gas, but it is initially 120 degrees misaligned with the stellar rotation axis. This misalignment persists for about 2 Gyr before the gas-star misalignment angle drops below 20 degrees. The time it takes for the gaseous and stellar components to align is much longer than previously thought, because the gas disc is accreting a significant amount of mass for about 1.5 Gyr after the merger, during which the angular momentum change induced by accreted gas dominates over that induced by stellar torques. Once the gas accretion rate has decreased sufficiently, the gas disc decouples from the surrounding halo gas and realigns with the stellar component in about 6 dynamical times. During the late evolution of the misaligned gas disc, the centre aligns faster than the outskirts, resulting in a warped disc. We discuss the observational consequences of the long survival of our misaligned gas disc and how our results can be used to calibrate merger rate estimates from observed gas misalignments.

preprint2015arXiv

The Difficulty of Getting High Escape Fractions of Ionizing Photons from High-redshift Galaxies: a View from the FIRE Cosmological Simulations

We present a series of high-resolution (20-2000 Msun, 0.1-4 pc) cosmological zoom-in simulations at z~6 from the Feedback In Realistic Environment (FIRE) project. These simulations cover halo masses 10^9-10^11 Msun and rest-frame ultraviolet magnitude Muv = -9 to -19. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback, which produce reasonable galaxy properties at z = 0-6. We post-process the snapshots with a radiative transfer code to evaluate the escape fraction (fesc) of hydrogen ionizing photons. We find that the instantaneous fesc has large time variability (0.01%-20%), while the time-averaged fesc over long time-scales generally remains ~5%, considerably lower than the estimate in many reionization models. We find no strong dependence of fesc on galaxy mass or redshift. In our simulations, the intrinsic ionizing photon budgets are dominated by stellar populations younger than 3 Myr, which tend to be buried in dense birth clouds. The escaping photons mostly come from populations between 3-10 Myr, whose birth clouds have been largely cleared by stellar feedback. However, these populations only contribute a small fraction of intrinsic ionizing photon budgets according to standard stellar population models. We show that fesc can be boosted to high values, if stellar populations older than 3 Myr produce more ionizing photons than standard stellar population models (as motivated by, e.g., models including binaries). By contrast, runaway stars with velocities suggested by observations can enhance fesc by only a small fraction. We show that "sub-grid" star formation models, which do not explicitly resolve star formation in dense clouds with n >> 1 cm^-3, will dramatically over-predict fesc.

preprint2015arXiv

The Fundamentally Different Dynamics of Dust and Gas in Molecular Clouds

We study the behavior of large dust grains in turbulent molecular clouds (MCs). In primarily neutral regions, dust grains move as aerodynamic particles, not necessarily with the gas. We therefore directly simulate, for the first time, the behavior of aerodynamic grains in highly supersonic, magnetohydrodynamic turbulence typical of MCs. We show that, under these conditions, grains with sizes a >0.01 micron exhibit dramatic (exceeding factor ~1000) fluctuations in the local dust-to-gas ratio (implying large small-scale variations in abundances, dust cooling rates, and dynamics). The dust can form highly filamentary structures (which would be observed in both dust emission and extinction), which can be much thinner than the characteristic width of gas filaments. Sometimes, the dust and gas filaments are not even in the same location. The 'clumping factor' of the dust (critical for dust growth/coagulation/shattering) can reach ~100, for grains in the ideal size range. The dust clustering is maximized around scales ~0.2pc*(a/micron)*(100 cm^-3/n_gas), and is 'averaged out' on larger scales. However, because the density varies widely in supersonic turbulence, the dynamic range of scales (and interesting grain sizes) for these fluctuations is much broader than in the subsonic case. Our results are applicable to MCs of essentially all sizes and densities, but we note how Lorentz forces and other physics (neglected here) may change them in some regimes. We discuss the potentially dramatic consequences for star formation, dust growth and destruction, and dust-based observations of MCs.

preprint2015arXiv

The Origin and Evolution of the Galaxy Mass-Metallicity Relation

We use high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass-metallicity relations (MZR) from z=0-6. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback. The simulations cover halo masses Mhalo=10^9-10^13 Msun and stellar mass Mstar=10^4-10^11 Msun at z=0 and have been shown to produce many observed galaxy properties from z=0-6. For the first time, our simulations agree reasonably well with the observed mass-metallicity relations at z=0-3 for a broad range of galaxy masses. We predict the evolution of the MZR from z=0-6 as log(Zgas/Zsun)=12+log(O/H)-9.0=0.35[log(Mstar/Msun)-10]+0.93 exp(-0.43 z)-1.05 and log(Zstar/Zsun)=[Fe/H]-0.2=0.40[log(Mstar/Msun)-10]+0.67 exp(-0.50 z)-1.04, for gas-phase and stellar metallicity, respectively. Our simulations suggest that the evolution of MZR is associated with the evolution of stellar/gas mass fractions at different redshifts, indicating the existence of a universal metallicity relation between stellar mass, gas mass, and metallicities. In our simulations, galaxies above Mstar=10^6 Msun are able to retain a large fraction of their metals inside the halo, because metal-rich winds fail to escape completely and are recycled into the galaxy. This resolves a long-standing discrepancy between "sub-grid" wind models (and semi-analytic models) and observations, where common sub-grid models cannot simultaneously reproduce the MZR and the stellar mass functions.

preprint2014arXiv

Galactic r-process enrichment by neutron star mergers in cosmological simulations of a Milky Way-mass galaxy

We quantify the stellar abundances of neutron-rich r-process nuclei in cosmological zoom-in simulations of a Milky Way-mass galaxy from the Feedback In Realistic Environments project. The galaxy is enriched with r-process elements by binary neutron star (NS) mergers and with iron and other metals by supernovae. These calculations include key hydrodynamic mixing processes not present in standard semi-analytic chemical evolution models, such as galactic winds and hydrodynamic flows associated with structure formation. We explore a range of models for the rate and delay time of NS mergers, intended to roughly bracket the wide range of models consistent with current observational constraints. We show that NS mergers can produce [r-process/Fe] abundance ratios and scatter that appear reasonably consistent with observational constraints. At low metallicity, [Fe/H]<-2, we predict there is a wide range of stellar r-process abundance ratios, with both supersolar and subsolar abundances. Low-metallicity stars or stars that are outliers in their r-process abundance ratios are, on average, formed at high redshift and located at large galactocentric radius. Because NS mergers are rare, our results are not fully converged with respect to resolution, particularly at low metallicity. However, the uncertain rate and delay time distribution of NS mergers introduces an uncertainty in the r-process abundances comparable to that due to finite numerical resolution. Overall, our results are consistent with NS mergers being the source of most of the r-process nuclei in the Universe.

preprint2014arXiv

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_sun, and stellar masses ~10^4-10^11. Our simulations include fully explicit treatment of both the multi-phase ISM (molecular through hot) and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of stellar feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M_halo~10^12 M_sun (including dwarfs, satellites, MW-mass disks, and small groups). By extension, this leads to reasonable agreement with the stellar mass function for M_star<10^11 M_sun. We predict weak redshift evolution in the M_star-M_halo relation, consistent with current constraints to z>6. We find that the M_star-M_halo relation is insensitive to numerical details, but is sensitive to the feedback physics. Simulations with only supernova feedback fail to reproduce the observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photo-heating and radiation pressure) is necessary to disrupt GMCs and enable efficient coupling of later supernovae to the gas. Star formation rates 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 in the simulation, and is instead determined by self-regulation via stellar feedback. Feedback reduces star formation rates considerably and produces a reservoir of gas that leads to rising late-time star formation histories significantly different from the halo accretion history. Feedback also produces large short-timescale variability in galactic SFRs, especially in dwarfs. Many of these properties are not captured by common 'sub-grid' galactic wind models.

preprint2014arXiv

No More Active Galactic Nuclei in Clumpy Disks Than in Smooth Galaxies at z~2 in CANDELS / 3D-HST

We use CANDELS imaging, 3D-HST spectroscopy, and Chandra X-ray data to investigate if active galactic nuclei (AGNs) are preferentially fueled by violent disk instabilities funneling gas into galaxy centers at 1.3<z<2.4. We select galaxies undergoing gravitational instabilities using the number of clumps and degree of patchiness as proxies. The CANDELS visual classification system is used to identify 44 clumpy disk galaxies, along with mass-matched comparison samples of smooth and intermediate morphology galaxies. We note that, despite being being mass-matched and having similar star formation rates, the smoother galaxies tend to be smaller disks with more prominent bulges compared to the clumpy galaxies. The lack of smooth extended disks is probably a general feature of the z~2 galaxy population, and means we cannot directly compare with the clumpy and smooth extended disks observed at lower redshift. We find that z~2 clumpy galaxies have slightly enhanced AGN fractions selected by integrated line ratios (in the mass-excitation method), but the spatially resolved line ratios indicate this is likely due to extended phenomena rather than nuclear AGNs. Meanwhile the X-ray data show that clumpy, smooth, and intermediate galaxies have nearly indistinguishable AGN fractions derived from both individual detections and stacked non-detections. The data demonstrate that AGN fueling modes at z~1.85 - whether violent disk instabilities or secular processes - are as efficient in smooth galaxies as they are in clumpy galaxies.

preprint2014arXiv

Some Stars are Totally Metal: A New Mechanism Driving Dust Across Star-Forming Clouds, and Consequences for Planets, Stars, and Galaxies

Dust grains in neutral gas behave as aerodynamic particles, so they can develop large density fluctuations independent of gas density fluctuations. Specifically, gas turbulence can drive order-of-magnitude 'resonant' fluctuations in the dust on scales where the gas stopping/drag timescale is comparable to the turbulent eddy turnover time. Here we show that for large grains (size >0.1 micron, containing most grain mass) in sufficiently large molecular clouds (radii >1-10 pc, masses >10^4 M_sun), this scale becomes larger than the characteristic sizes of pre-stellar cores (the sonic length), so large fluctuations in the dust-to-gas ratio are imprinted on cores. As a result, star clusters and protostellar disks formed in large clouds should exhibit significant abundance spreads in the elements preferentially found in large grains. This naturally predicts populations of carbon-enhanced stars, certain highly unusual stellar populations observed in nearby open clusters, and may explain the 'UV upturn' in early-type galaxies. It will also dramatically change planet formation in the resulting protostellar disks, by preferentially 'seeding' disks with an enhancement in large carbonaceous or silicate grains. The relevant threshold for this behavior scales simply with cloud densities and temperatures, making straightforward predictions for clusters in starbursts and high-redshift galaxies. Because of the selective sorting by size, this process is not necessarily visible in extinction mapping. We also predict the shape of the abundance distribution -- when these fluctuations occur, a small fraction of the cores may actually be seeded with abundances ~100 times the mean, such that they are almost 'totally metal' (Z~1)! Assuming the cores collapse, these totally metal stars would be rare (1 in 10^4 in clusters where this occurs), but represent a fundamentally new stellar evolution channel.

preprint2014arXiv

The Origin and Universality of the Stellar Initial Mass Function

We review current theories for the origin of the Stellar Initial Mass Function (IMF) with particular focus on the extent to which the IMF can be considered universal across various environments. To place the issue in an observational context, we summarize the techniques used to determine the IMF for different stellar populations, the uncertainties affecting the results, and the evidence for systematic departures from universality under extreme circumstances. We next consider theories for the formation of prestellar cores by turbulent fragmentation and the possible impact of various thermal, hydrodynamic and magneto-hydrodynamic instabilities. We address the conversion of prestellar cores into stars and evaluate the roles played by different processes: competitive accretion, dynamical fragmentation, ejection and starvation, filament fragmentation and filamentary accretion flows, disk formation and fragmentation, critical scales imposed by thermodynamics, and magnetic braking. We present explanations for the characteristic shapes of the Present-Day Prestellar Core Mass Function and the IMF and consider what significance can be attached to their apparent similarity. Substantial computational advances have occurred in recent years, and we review the numerical simulations that have been performed to predict the IMF directly and discuss the influence of dynamics, time-dependent phenomena, and initial conditions.

preprint2013arXiv

A General Theory of Turbulent Fragmentation

We develop an analytic framework to understand fragmentation in turbulent, self-gravitating media. Previously, we showed some properties of turbulence can be predicted with the excursion-set formalism. Here, we generalize to fully time-dependent gravo-turbulent fragmentation & collapse. We show that turbulent systems are always gravitationally unstable (in a probabilistic sense). The fragmentation mass spectra, size/mass relations, correlation functions, range of scales over which fragmentation occurs, & time-dependent rates of fragmentation are predictable. We show how this depends on bulk turbulent properties (Mach numbers & power spectra). We also generalize to include rotation, complicated equations of state, collapsing/expanding backgrounds, magnetic fields, intermittency, & non-normal statistics. We derive how fragmentation is suppressed with 'stiffer' equations of state or different driving mechanisms. Suppression appears at an 'effective sonic scale' where Mach(R,rho)~1. Gas becomes stable below this scale for polytropic gamma>4/3, but fragmentation still occurs on larger scales. The scale-free nature of turbulence and gravity generically drives mass spectra and correlation functions towards universal shapes, with weak dependence on many properties of the media. Correlated fluctuation structures, non-Gaussian density distributions, & intermittency have surprisingly small effects on the fragmentation process. This is because fragmentation cascades on small scales are 'frozen in' when large-scale modes push the 'parent' region above the collapse threshold; though they collapse, their statistics are only weakly modified by the collapse process. With thermal support, structure develops 'top-down' in time via fragmentation cascades; but strong rotational support reverses this to 'bottom-up' growth via mergers & introduces a maximal instability scale distinct from the Toomre scale.

preprint2013arXiv

A Model for (Non-Lognormal) Density Distributions in Isothermal Turbulence

We propose a new, physically motivated fitting function for density PDFs in turbulent gas. Although it is known that when gas is isothermal, the PDF is approximately lognormal in the core, high-resolution simulations show large deviations from exact log-normality. The proposed function provides an extraordinarily accurate description of the density PDFs in simulations with Mach numbers ~0.1-15 and dispersion in log(rho) from ~0.01-4 dex. Compared to a lognormal or lognormal-skew-kurtosis model, the fits are improved by orders of magnitude in the wings of the distribution (with fewer free parameters). This is true in simulations using distinct numerical methods, including or excluding magnetic fields. Deviations from lognormality are represented by a parameter T that increases with the compressive Mach number. The proposed distribution can be derived from intermittent cascade models of the longitudinal (compressive) velocity differences, which should be directly related to density fluctuations, and we provide a simple interpretation of the density PDF as the product of a continuous-time relaxation process. As such the parameter T is consistent with the same parameter needed to explain the (intermittent) velocity structure functions; its behavior is consistent with turbulence becoming more intermittent as it becomes more dominated by strong shocks. It provides a new and unique probe of the role of intermittency in the density (not just velocity) structure of turbulence. We show that this naturally explains some apparent contradictory results in the literature based on use of different moments of the density PDF, as well as differences based on whether volume-weighted or mass-weighted quantities are measured. We show how these are fundamentally related to the fact that mass conservation requires violations of log-normal statistics.

preprint2013arXiv

Accretion Does Not Drive the Turbulence in Galactic Disks

Rapid accretion of cold gas plays a crucial role in getting gas into galaxies. It has been suggested that this accretion proceeds along narrow streams that might also directly drive the turbulence in galactic gas, dynamical disturbances, and bulge formation. In cosmological simulations, however, it is impossible to isolate and hence disentangle the effect of accretion from internal instabilities and mergers. Moreover, in most cosmological simulations, the phase structure and turbulence in the ISM arising from stellar feedback are treated in a sub-grid manner, so that feedback cannot generate ISM turbulence. In this paper we therefore test the effects of cold streams in extremely high-resolution simulations of otherwise isolated galaxy disks using detailed models for star formation and feedback; we then include or exclude mock cold flows falling onto the galaxies with accretion rates, velocities and geometry set to maximize their effect on the disk. We find: (1) Turbulent velocity dispersions in gas disks are identical with or without the cold flow; the energy injected by the flow is dissipated where it meets the disk. (2) In runs without stellar feedback, the presence of a cold flow has essentially no effect on runaway local collapse, resulting in star formation rates (SFRs) that are far too large. (3) Disks in runs with feedback and cold flows have higher SFRs, but only insofar as they have more gas. (4) Because flows are extended relative to the disk, they do not trigger strong resonant responses and so induce weak morphological perturbation (bulge formation via instabilities is not accelerated). (5) However, flows can thicken the disk by direct contribution of out-of-plane streams. We conclude that while inflows are critical over cosmological timescales to determine the supply and angular momentum of gas disks, they have weak instantaneous dynamical effects on galaxies.

preprint2013arXiv

Dense Molecular Gas: A Sensitive Probe of Stellar Feedback Models

We show that the mass fraction of GMC gas (n>100 cm^-3) in dense (n>>10^4 cm^-3) star-forming clumps, observable in dense molecular tracers (L_HCN/L_CO(1-0)), is a sensitive probe of the strength and mechanism(s) of stellar feedback. Using high-resolution galaxy-scale simulations with pc-scale resolution and explicit models for feedback from radiation pressure, photoionization heating, stellar winds, and supernovae (SNe), we make predictions for the dense molecular gas tracers as a function of GMC and galaxy properties and the efficiency of stellar feedback. In models with weak/no feedback, much of the mass in GMCs collapses into dense sub-units, predicting L_HCN/L_CO(1-0) ratios order-of-magnitude larger than observed. By contrast, models with feedback properties taken directly from stellar evolution calculations predict dense gas tracers in good agreement with observations. Changing the strength or timing of SNe tends to move systems along, rather than off, the L_HCN-L_CO relation (because SNe heat lower-density material, not the high-density gas). Changing the strength of radiation pressure (which acts efficiently in the highest density gas), however, has a much stronger effect on L_HCN than on L_CO. We predict that the fraction of dense gas (L_HCN/L_CO(1-0)) increases with increasing GMC surface density; this drives a trend in L_HCN/L_CO(1-0) with SFR and luminosity which has tentatively been observed. Our results make specific predictions for enhancements in the dense gas tracers in unusually dense environments such as ULIRGs and galactic nuclei (including the galactic center).

preprint2013arXiv

Modeling Mid-Infrared Diagnostics of Obscured Quasars and Starbursts

We analyze the link between active galactic nuclei (AGN) and mid-infrared flux using dust radiative transfer calculations of starbursts realized in hydrodynamical simulations. Focusing on the effects of galaxy dust, we evaluate diagnostics commonly used to disentangle AGN and star formation in ultraluminous infrared galaxies (ULIRGs). We examine these quantities as a function of time, viewing angle, dust model, AGN spectrum, and AGN strength in merger simulations representing two possible extremes of the ULIRG population: one is a typical gas-rich merger at z ~ 0, and the other is characteristic of extremely obscured starbursts at z ~ 2 to 4. This highly obscured burst begins star-formation-dominated with significant PAH emission, and ends with a ~10^9 yr period of red near-IR colors. At coalescence, when the AGN is most luminous, dust obscures the near-infrared AGN signature, reduces the relative emission from polycyclic aromatic hydrocarbons (PAHs), and enhances the 9.7 micron absorption by silicate grains. Although generally consistent with previous interpretations, our results imply none of these indicators can unambiguously estimate the AGN luminosity fraction in all cases. Motivated by the simulations, we show that a combination of the extinction feature at 9.7 micron, the PAH strength, and a near-infrared slope can simultaneously constrain the AGN fraction and dust grain distribution for a wide range of obscuration. We find that this indicator, accessible to the James Webb Space Telescope, may estimate the AGN power as tightly as the hard X-ray flux alone, thereby providing a valuable future cross-check and constraint for large samples of distant ULIRGs.

preprint2013arXiv

Resolving The Generation of Starburst Winds in Galaxy Mergers

We study galaxy super-winds driven in major mergers, using pc-resolution simulations with detailed models for stellar feedback that can self-consistently follow the formation/destruction of GMCs and generation of winds. The models include molecular cooling, star formation at high densities in GMCs, and gas recycling and feedback from SNe (I&II), stellar winds, and radiation pressure. We study mergers of systems from SMC-like dwarfs and Milky Way analogues to z~2 starburst disks. Multi-phase super-winds are generated in all passages, with outflow rates up to ~1000 M_sun/yr. However, the wind mass-loading efficiency (outflow rate divided by SFR) is similar to that in isolated galaxy counterparts of each merger: it depends more on global galaxy properties (mass, size, escape velocity) than on the dynamical state of the merger. Winds tend to be bi- or uni-polar, but multiple 'events' build up complex morphologies with overlapping, differently-oriented bubbles/shells at a range of radii. The winds have complex velocity and phase structure, with material at a range of speeds up to ~1000 km/s, and a mix of molecular, ionized, and hot gas that depends on galaxy properties and different feedback mechanisms. These simulations resolve a problem in some 'sub-grid' models, where simple wind prescriptions can dramatically suppress merger-induced starbursts. But despite large mass-loading factors (>~10) in the winds, the peak SFRs are comparable to those in 'no wind' simulations. Wind acceleration does not act equally, so cold dense gas can still lose angular momentum and form stars, while blowing out gas that would not have participated in the starburst in the first place. Considerable wind material is not unbound, and falls back on the disk at later times post-merger, leading to higher post-starburst SFRs in the presence of stellar feedback. This may require AGN feedback to explain galaxy quenching.

preprint2013arXiv

Star Formation in Galaxy Mergers with Realistic Models of Stellar Feedback & the Interstellar Medium

We use simulations with realistic models for stellar feedback to study galaxy mergers. These high resolution (1 pc) simulations follow formation and destruction of individual GMCs and star clusters. The final starburst is dominated by in situ star formation, fueled by gas which flows inwards due to global torques. The resulting high gas density results in rapid star formation. The gas is self gravitating, and forms massive (~10^10 M_sun) GMCs and subsequent super-starclusters (masses up to 10^8 M_sun). However, in contrast to some recent simulations, the bulk of new stars which eventually form the central bulge are not born in superclusters which then sink to the center of the galaxy, because feedback efficiently disperses GMCs after they turn several percent of their mass into stars. Most of the mass that reaches the nucleus does so in the form of gas. The Kennicutt-Schmidt law emerges naturally as a consequence of feedback balancing gravitational collapse, independent of the small-scale star formation microphysics. The same mechanisms that drive this relation in isolated galaxies, in particular radiation pressure from IR photons, extend over seven decades in SFR to regulate star formation in the most extreme starbursts (densities >10^4 M_sun/pc^2). Feedback also drives super-winds with large mass loss rates; but a significant fraction of the wind material falls back onto the disks at later times, leading to higher post-starburst SFRs in the presence of stellar feedback. Strong AGN feedback is required to explain sharp cutoffs in star formation rate. We compare the predicted relic structure, mass profile, morphology, and efficiency of disk survival to simulations which do not explicitly resolve GMCs or feedback. Global galaxy properties are similar, but sub-galactic properties and star formation rates can differ significantly.

preprint2013arXiv

The Meaning and Consequences of Star Formation Criteria in Galaxy Models with Resolved Stellar Feedback

We consider the effects of different star formation criteria on galactic scales, in high-resolution simulations with explicitly resolved GMCs and stellar feedback. We compare: (1) a self-gravity criterion (based on the local virial parameter and the assumption that self-gravitating gas collapses to high density in a free-fall time), (2) a fixed density threshold, (3) a molecular-gas law, (4) a temperature threshold, (5) a Jeans-instability requirement, (6) a criteria that cooling times be shorter than dynamical times, and (7) a convergent-flow criterion. We consider these both MW-like and high-density (starburst) galaxies. With feedback present, all models produce identical integrated star formation rates (SFRs), in agreement with the Kennicutt relation. Without feedback all produce orders-of-magnitude excessive SFRs. This is totally dependent on feedback and independent of the SF law. However, the spatial and density distribution of SF depend strongly on the SF criteria. Because cooling rates are generally fast and gas is turbulent, criteria (4)-(7) are 'weak' and spread SF uniformly over the disk (above densities n~0.01-0.1 cm^-3). A molecular criterion (3) localizes to higher densities, but still a wide range; for Z Z_solar, it is similar to a density threshold at n~1 cm^-3 (well below mean densities in the MW center or starbursts). Fixed density thresholds (2) can always select the highest densities, but must be adjusted for simulation resolution and galaxy properties; the same threshold that works in a MW-like simulation will select nearly all gas in a starburst. Binding criteria (1) tend to adaptively select the largest over-densities, independent of galaxy model or resolution, and automatically predict clustered SF. We argue that this SF model is most physically-motivated and presents significant numerical advantages in large-dynamic range simulations.

preprint2013arXiv

Turbulent Disks are Never Stable: Fragmentation and Turbulence-Promoted Planet Formation

A fundamental assumption in our understanding of disks is that when the Toomre Q>>1, the disk is stable against fragmentation into self-gravitating objects (and so cannot form planets via direct collapse). But if disks are turbulent, this neglects a spectrum of stochastic density fluctuations that can produce rare, high-density mass concentrations. Here, we use a recently-developed analytic framework to predict the statistics of these fluctuations, i.e. the rate of fragmentation and mass spectrum of fragments formed in a turbulent Keplerian disk. Turbulent disks are never completely stable: we calculate the (always finite) probability of forming self-gravitating structures via stochastic turbulent density fluctuations in such disks. Modest sub-sonic turbulence above Mach number ~0.1 can produce a few stochastic fragmentation or 'direct collapse' events over ~Myr timescales, even if Q>>1 and cooling is slow (t_cool>>t_orbit). In trans-sonic turbulence this extends to Q~100. We derive the true Q-criterion needed to suppress such events, which scales exponentially with Mach number. We specify to turbulence driven by MRI, convection, or spiral waves, and derive equivalent criteria in terms of Q and the cooling time. Cooling times >~50*t_dyn may be required to completely suppress fragmentation. These gravoturbulent events produce mass spectra peaked near ~M_disk*(Q*M_disk/M_star)^2 (rocky-to-giant planet masses, increasing with distance from the star). We apply this to protoplanetary disk models and show that even minimum mass solar nebulae could experience stochastic collapse events, provided a source of turbulence.

preprint2013arXiv

Variations in the Stellar IMF: from Bottom to Top

We use a recently-developed analytic model for the ISM structure from scales of GMCs through star-forming cores to explore how the pre-stellar core mass function (CMF) and, by extrapolation, stellar initial mass function (IMF) should depend on both local and galactic properties. If the ISM is supersonically turbulent, the statistical properties of the density field follow from the turbulent velocity spectrum, and the excursion set formalism can be applied to analytically calculate the mass function of collapsing cores on the smallest scales on which they are self-gravitating (non-fragmenting). Two parameters determine the model: the disk-scale Mach number M_h (which sets the shape of the CMF), and the absolute velocity (to assign an absolute scale). For 'normal' variation in disk properties and core gas temperatures in the MW and local galaxies, there is almost no variation in the predicted high-mass behavior of the CMF/IMF. The slope is always close to Salpeter down to <1 M_sun. We predict modest variation in the sub-solar regime, mostly from variation in M_h, but within the observed scatter in sub-solar IMFs in local regions. For fixed galaxy properties, there is little variation in shape or 'upper mass limit' with parent GMC mass. However, in extreme starbursts (e.g. ULIRGs) we predict a bottom-heavy CMF. This agrees with the IMF inferred for the centers of Virgo ellipticals, believed to form in such a nuclear starburst. The CMF is bottom heavy despite the gas temperature being an order of magnitude larger, because M_h is also much larger. Larger M_h values make the 'parent' cloud mass (turbulent Jeans mass) larger, but promote fragmentation to smaller scales; this steepens the slope of the low-mass CMF and shifts the turnover mass. The model may predict a top-heavy CMF for the sub-pc disks around Sgr A*, but the relevant input parameters are uncertain.

preprint2012arXiv

3-D Radiative Transfer Calculations of Radiation Feedback from Massive Black Holes: Outflow of Mass from the Dusty "Torus"

Observational and theoretical arguments suggest that the momentum carried in mass outflows from AGN can reach several times L / c, corresponding to outflow rates of hundreds of solar masses per year. Radiation pressure on lines alone may not be sufficient to provide this momentum deposition, and the transfer of reprocessed IR radiation in dusty nuclear gas has been postulated to provide the extra enhancement. The efficacy of this mechanism, however, will be sensitive to multi-dimensional effects such as the tendency for the reprocessed radiation to preferentially escape along sight-lines of lower column density. We use Monte Carlo radiative transfer calculations to determine the radiation force on dusty gas residing within approximately 10 parsecs from an accreting super-massive black hole. We calculate the net rate of momentum deposition in the surrounding gas and estimate the mass-loss rate in the resulting outflow as a function of solid angle for different black hole luminosities, sightline-averaged column densities, clumping parameters, and opening angles of the dusty gas. We find that these dust-driven winds carry momentum fluxes of 1-5 times L / c and correspond to mass-loss rates of 10-100 solar masses per year for a 10^8 solar mass black hole radiating at or near its Eddington limit. These results help to explain the origin of high velocity molecular and atomic outflows in local ULIRGs, and can inform numerical simulations of galaxy evolution including AGN feedback.

preprint2012arXiv

An Analytic Model of Angular Momentum Transport by Gravitational Torques: From Galaxies to Massive Black Holes

We present analytic calculations of angular momentum transport and gas inflow in galaxies, from scales of ~kpc to deep in the potential of a central black hole (BH). We compare these analytic calculations to numerical simulations and use them to develop a sub-grid model of BH growth that can be incorporated into semi-analytic models or cosmological simulations. Both analytic calculations and simulations argue that the strongest torque on gas arises when non-axisymmetric perturbations to the stellar gravitational potential produces orbit crossings and shocks in the gas. This is true both at large radii ~0.01-1 kpc, where bar-like modes dominate the non-axisymmetric potential, and at smaller radii <10 pc, where a lopsided/eccentric disk dominates. The traditional orbit crossing criterion is not always adequate to predict the locations of, and inflow due to, shocks in gas+stellar disks with finite sound speeds. We derive a modified criterion that predicts the presence of shocks in stellar dominated systems even absent formal orbit crossing. We then derive analytic expressions for the loss of angular momentum and the resulting gas inflow rates in the presence of such shocks. We test our analytic predictions using hydrodynamic simulations at a range of galactic scales, and show that they successfully predict the mass inflow rates and quasi-steady gas surface densities with small scatter (0.3 dex). We use our analytic results to construct a new estimate of the BH accretion rate given galaxy properties at larger radii. This captures the key scalings in the numerical simulations. Alternate estimates such as the local viscous accretion rate or the spherical Bondi rate fail systematically to reproduce the simulations.

preprint2012arXiv

An Excursion-Set Model for the Structure of GMCs and the ISM

The ISM is governed by supersonic turbulence on a range of scales. We use this to develop a rigorous excursion-set model for the formation and time evolution of dense gas structures (GMCs, massive clumps, and cores). Supersonic turbulence drives the density distribution to a lognormal with dispersion increasing with Mach number; we generalize this to include scales >h (the disk scale height), and use it to construct the statistical properties of the density field smoothed on a scale R. We then compare conditions for self-gravitating collapse including thermal, turbulent, and rotational support. We show this becomes a well-defined barrier crossing problem. As such, an exact 'bound object mass function' can be derived, from scales of the sonic length to above the disk Jeans mass. This agrees remarkably well with observed GMC mass functions in the MW and other galaxies; the only inputs are the mass and size of the galaxies (to normalize the model). This explains the mass function cutoff and its power-law slope (close to, but shallower than, -2). The model also predicts the linewidth-size and size-mass relations of clouds and the dependence of their residuals on surface density/pressure. We use this to predict the spatial correlation function/clustering of clouds and star clusters; these also agree well with observations. We predict the size/mass function of ISM 'bubbles' or 'holes', and show this can account for observed HI hole distributions without any local feedback. We generalize the model to construct time-dependent 'merger/fragmentation trees' which can be used to follow cloud evolution and construct semi-analytic models for the ISM. We provide explicit recipes to construct the trees. We use a simple example to show that, if clouds are not destroyed in ~1-5 crossing times, then all ISM mass would be trapped in collapsing objects even if the large-scale turbulence were maintained.

preprint2012arXiv

An Explanation for the Slopes of Stellar Cusps in Galaxy Spheroids

The stellar surface mass density profiles at the centers of typical ~L* and lower-mass spheroids exhibit power law 'cusps' with $Σ\propto R^(-n)$, where 0.5<n<1 for radii ~1-100 pc. Observations and theory support models in which these cusps are formed by dissipative gas inflows and nuclear starbursts in gas-rich mergers. At these comparatively large radii, stellar relaxation is unlikely to account for or strongly modify the cuspy stellar profiles. We argue that the power-law surface density profiles observed are a natural consequence of the gravitational instabilities that dominate angular momentum transport in the gravitational potential of a central massive black hole. The dominant mode at these radii is an m=1 lopsided/eccentric disk instability, in which stars torquing the gas can drive rapid inflow and accretion. Such a mode first generically appears at large radii and propagates inwards by exciting eccentricities at smaller and smaller radii, where M*(<R)<<M_BH. When the stellar surface density profile is comparatively shallow with n<1/2, the modes cannot efficiently propagate to R=0 and so gas piles up and star formation steepens the profile. But if the profile is steeper than n=1, the inwards propagation of eccentricity is strongly damped, suppressing inflow and bringing n down again. Together these results produce an equilibrium slope of 1/2 < n < 1 in the potential of the central black hole. These physical arguments are supported by nonlinear numerical simulations of gas inflow in galactic nuclei. Together, these results naturally explain the observed stellar density profiles of 'cusp' elliptical galaxies.

preprint2012arXiv

CANDELS: The progenitors of compact quiescent galaxies at z~2

We combine high-resolution HST/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (log(M*) > 10) galaxies at redshifts z = 1.4 - 3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5 - 3. At z > 2 most cSFGs have specific star-formation rates (sSFR = 10^-9 yr^-1) half that of typical, massive SFGs at the same epoch, and host X-ray luminous AGN 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2 - 3 and fade to cQGs by z = 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary scenarios of QG formation: an early (z > 2), fast-formation path of rapidly-quenched cSFGs that evolve into cQGs that later enlarge within the quiescent phase, and a slow, late-arrival (z < 2) path for SFGs to form QGs without passing through a compact state.

preprint2012arXiv

Stellar Feedback & Bulge Formation in Clumpy Disks

We use numerical simulations of isolated galaxies to study the effects of stellar feedback on the formation and evolution of giant star-forming gas 'clumps' in high-redshift, gas-rich galaxies. Such galactic disks are unstable to the formation of bound gas-rich clumps whose properties initially depend only on global disk properties, not the microphysics of feedback. In simulations without stellar feedback, clumps turn an order-unity fraction of their mass into stars and sink to the center, forming a large bulge and kicking most of the stars out into a much more extended stellar envelope. By contrast, strong radiative stellar feedback disrupts even the most massive clumps after they turn ~10-20% of their mass into stars, in a timescale of ~10-100 Myr, ejecting some material into a super-wind and recycling the rest of the gas into the diffuse ISM. This suppresses the bulge formation rate by direct 'clump coalescence' by a factor of several. However, the galactic disks do undergo significant internal evolution in the absence of mergers: clumps form and disrupt continuously and torque gas to the galactic center. The resulting evolution is qualitatively similar to bar/spiral evolution in simulations with a more homogeneous ISM.

preprint2012arXiv

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 super-winds that enrich the IGM and shape the galaxy mass function and mass-metallicity relation. In previous papers, we introduced new numerical methods for implementing stellar feedback on sub-GMC through galactic scales in galaxy simulations. This includes radiation pressure (UV through IR), SNe (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 SFR. The mass-loading efficiency (wind mass loss rate divided by SFR) scales inversely with circular velocity, consistent with momentum-conservation expectations. We study the contributions of each feedback mechanism to galactic winds in a range of galaxy models, from SMC-like dwarfs & MW-analogues to z~2 clumpy disks. In massive, gas-rich systems (local starbursts and high-z galaxies), radiation pressure dominates the wind generation. For MW-like spirals and dwarf galaxies the gas densities are much lower, and shock-heated gas from SNe and stellar winds dominates production of large-scale outflows. In all models, however, winds have a multi-phase structure that depends on interactions between multiple feedback mechanisms operating on different spatial & time scales: any single mechanism fails to reproduce the winds observed. We provide fitting functions for 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 explicit dependence on gas surface density that can be very important in both low-density dwarf galaxies and high-density gas-rich galaxies.

preprint2012arXiv

Submillimetre galaxies in a hierarchical universe: number counts, redshift distribution, and implications for the IMF

High-redshift submillimetre galaxies (SMGs) are some of the most rapidly star-forming galaxies in the Universe. Historically, galaxy formation models have had difficulty explaining the observed number counts of SMGs. We combine a semi-empirical model with 3-D hydrodynamical simulations and 3-D dust radiative transfer to predict the number counts of unlensed SMGs. Because the stellar mass functions, gas and dust masses, and sizes of our galaxies are constrained to match observations, we can isolate uncertainties related to the dynamical evolution of galaxy mergers and the dust radiative transfer. The number counts and redshift distributions predicted by our model agree well with observations. Isolated disc galaxies dominate the faint (S_1.1 < ~1 mJy, or S_850 < ~2 mJy) population. The brighter sources are a mix of merger-induced starbursts and galaxy-pair SMGs; the latter subpopulation accounts for ~30 - 50 per cent of all SMGs at all S_1.1 > ~0.5 mJy (S_850 > ~1 mJy). The mean redshifts are ~3.0 - 3.5, depending on the flux cut, and the brightest sources tend to be at higher redshifts. Because the galaxy-pair SMGs will be resolved into multiple fainter sources by ALMA, the bright ALMA counts should be as much as 2 times less than those observed using single-dish telescopes. The agreement between our model, which uses a Kroupa IMF, and observations suggests that the IMF in high-redshifts starbursts need not be top-heavy; if the IMF were top-heavy, our model would over-predict the number counts. We conclude that the difficulty some models have reproducing the observed SMG counts is likely indicative of more general problems -- such as an under-prediction of the abundance of massive galaxies or a star formation rate--stellar mass relation normalisation lower than that observed -- rather than a problem specific to the SMG population.

preprint2012arXiv

The Stellar IMF, Core Mass Function, & The Last-Crossing Distribution

Hennebelle & Chabrier 2008 (HC08) attempted to derive the stellar IMF as a consequence of turbulent density fluctuations, using an argument similar to Press & Schechter 1974 for Gaussian random fields. Like that example, however, this solution does not resolve the 'cloud in cloud' problem; it also does not extend to large scales that dominate the velocity/density fluctuations. In principle, these can change the results at the order-of-magnitude level. Here, we use the results from Hopkins 2011 (H11) to generalize the excursion set formalism and derive the exact solution in this regime. We argue that the stellar IMF and core mass function (CMF) should be associated with the last-crossing distribution, i.e. the mass spectrum of bound objects defined on the smallest scale on which they are self-gravitating. This differs from the first-crossing distribution (mass function on the largest self-gravitating scale) which is defined cosmologically and which H11 show corresponds to the GMC mass function in disks. We derive an analytic equation for the last-crossing distribution that can be used for an arbitrary collapse threshold in ISM and cosmological studies. With this, we show that the same model that predicts the GMC mass function and large-scale structure of galaxy disks also predicts the CMF (and by extrapolation IMF) in good agreement with observations. The only adjustable parameter in the model is the turbulent velocity power spectrum, which in the range p~5/3-2 gives similar results. We also use this to justify why the approximate solution in HC08 is reasonable (up to a normalization) over the CMF/IMF mass range; however there are significant corrections at intermediate and high masses. We discuss how the exact solutions here can be embedded into time-dependent models that follow density fluctuations, fragmentation, successive generations of star formation.

preprint2012arXiv

The Structure of the Interstellar Medium of Star Forming Galaxies

We present numerical methods for including stellar feedback in galaxy-scale simulations. We include heating by SNe (I & II), gas recycling and shock-heating from O-star & AGB winds, HII photoionization, and radiation pressure from stellar photons. The energetics and time-dependence are taken directly from stellar evolution models. We implement these in simulations with pc-scale resolution, modeling galaxies from SMC-like dwarfs and MW analogues to massive z~2 starburst disks. Absent feedback, gas cools and collapses without limit. With feedback, the ISM reaches a multi-phase steady state in which GMCs continuously form, disperse, and re-form. Our primary results include: (1) Star forming galaxies generically self-regulate at Toomre Q~1. Most of the volume is in diffuse hot gas with most of the mass in dense GMC complexes. The phase structure and gas mass at high densities are much more sensitive probes of stellar feedback physics than integrated quantities (Toomre Q or gas velocity dispersion). (2) Different feedback mechanisms act on different scales: radiation & HII pressure are critical to prevent runaway collapse of dense gas in GMCs. SNe and stellar winds dominate the dynamics of volume-filling hot gas; however this primarily vents out of the disk. (3) The galaxy-averaged SFR is determined by feedback. For given feedback efficiency, restricting star formation to molecular gas or modifying the cooling function has little effect; but changing feedback mechanisms directly translates to shifts off the Kennicutt-Schmidt relation. (4) Self-gravity leads to marginally-bound GMCs with an ~M^-2 mass function with a cutoff at the Jeans mass; they live a few dynamical times before being disrupted by stellar feedback and turn ~1-10% of their mass into stars (increasing from dwarfs through starburst galaxies). Low-mass GMCs are preferentially unbound.

preprint2012arXiv

Why Are AGN and Host Galaxies Misaligned?

It is well-established observationally that the characteristic angular momentum axis on small scales around AGN, traced by radio jets and the putative torus, is not well-correlated with the large-scale angular momentum axis of the host galaxy. In this paper, we show that such misalignments arise naturally in high-resolution simulations in which we follow angular momentum transport and inflows from galaxy to sub-pc scales near AGN, triggered either during galaxy mergers or by instabilities in isolated disks. Sudden misalignments can sometimes be caused by single massive clumps falling into the center slightly off-axis, but more generally, they arise even when the gas inflows are smooth and trace only global gravitational instabilities. When several nested, self-gravitating modes are present, the inner ones can precess and tumble in the potential of the outer modes. Resonant angular momentum exchange can flip or re-align the spin of an inner mode on a short timescale, even without the presence of massive clumps. We therefore do not expect that AGN and their host galaxies will be preferentially aligned, nor should the relative alignment be an indicator of the AGN fueling mechanism. We discuss implications of this conclusion for AGN feedback and BH spin evolution. The misalignments may mean that even BHs accreting from smooth large-scale disks will not be spun up to maximal rotation, and so have more modest radiative efficiencies and inefficient jet formation. Even more random orientations are possible if there is further, un-resolved clumpiness in the gas, and more ordered accretion may occur if the inflow is slower and not self-gravitating.

preprint2012arXiv

Why Do Stars Form In Clusters? An Analytic Model for Stellar Correlation Functions

Recently, we have shown that if the ISM is governed by super-sonic turbulent flows, the excursion-set formalism can be used to calculate the statistics of self-gravitating objects over a wide range of scales. On the largest self-gravitating scales ('first crossing'), these correspond to GMCs, and on the smallest non-fragmenting self-gravitating scales ('last crossing'), to protostellar cores. Here, we extend this formalism to rigorously calculate the auto and cross-correlation functions of cores (and by extension, young stars) as a function of spatial separation and mass, in analogy to the cosmological calculation of halo clustering. We show that this generically predicts that star formation is very strongly clustered on small scales: stars form in clusters, themselves inside GMCs. Outside the binary-star regime, the projected correlation function declines as a weak power-law, until a characteristic scale which corresponds to the characteristic mass scale of GMCs. On much larger scales the clustering declines such that star formation is not strongly biased on galactic scales, relative to the actual dense gas distribution. The precise correlation function shape depends on properties of the turbulent spectrum, but its qualitative behavior is quite general. The predictions agree well with observations of young star and core autocorrelation functions over ~4 dex in radius. Clustered star formation is a generic consequence of supersonic turbulence if most of the power in the velocity field, hence the contribution to density fluctuations, comes from large scales. The distribution of self-gravitating masses near the sonic length is then imprinted by fluctuations on larger scales. We similarly show that the fraction of stars formed in 'isolated' modes should be small (\lesssim10%).

preprint2011arXiv

A Multiwavelength Study of Binary Quasars and Their Environments

We present Chandra X-ray imaging and spectroscopy for 14 quasars in spatially resolved pairs, part of a complete sample of binary quasars with small transverse separations drawn from Sloan Digital Sky Survey (DR6) photometry. We find no significant difference in X-ray properties when compared with large control samples of isolated quasars. We present infrared photometry from our observations with SWIRC at the MMT, and from the WISE Preliminary Data Release, and fit simple spectral energy distributions to all 14 QSOs. We find preliminary evidence that substantial contributions from star formation are required, but possibly no more so than for isolated X-ray-detected QSOs. Sensitive searches of the X-ray images for extended emission, and the optical images for optical galaxy excess show that these binary QSOs are not preferentially found in rich cluster environments. While larger binary QSO samples with richer far-IR and sub-millimeter multiwavelength data might better reveal signatures of merging and triggering, optical color-selection of QSO pairs may be biased against such signatures. X-ray and/or variability selection of QSO pairs, while challenging, should be attempted. We present in our Appendix a primer on X-ray flux and luminosity calculations.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5σ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5σ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

preprint2011arXiv

Dynamical Delays Between Starburst and AGN Activity in Galaxy Nuclei

Observations of AGN have suggested a possible delay between the peak of star formation (on some scale) and AGN activity. Inefficient fueling (and/or feedback) from fast stellar winds has been invoked to explain this, but we argue this is unlikely in bright systems accreting primarily cold dense gas. We show that such a delay can arise even in bright quasars for purely dynamical reasons. If some large-scale process produces rapid inflow, smaller scales will quickly become gas-dominated. As the gas density peaks, so does the SFR. However, gravitational torques which govern further inflow are relatively inefficient in gas-dominated systems; as more gas is turned into stars, the stars provide an efficient angular momentum sink allowing more rapid inflow. Moreover, the gas provided to the central regions in mergers or strong disk instabilities will typically be ~100 times larger than that needed to fuel the BH; the system is effectively in the 'infinite gas supply' limit. BH growth can therefore continue for some time while the gas supply exhausts, until it has significantly depleted to the point where the BH is finally 'starved.' Both of these effects act together with comparable magnitude, and mean that the peak of BH growth can lag the peak in the SFR measured at a given scale by a timescale corresponding to the gas exhaustion time on that scale (~ 10-100 local dynamical times). This predicts that the inferred delay will vary in a specific manner with the radius over which the star formation rate is measured. We discuss possible implications for the role of AGN feedback in suppressing star formation activity.

preprint2011arXiv

Relation Between Globular Clusters and Supermassive Black Holes in Ellipticals as a Manifestation of the Black Hole Fundamental Plane

We analyze the relation between the mass of the central supermassive black hole (Mbh) and the number of globular clusters (Ngc) in elliptical galaxies and bulges as a ramification of the black hole fundamental plane, the theoretically predicted and observed multi-variable correlation between Mbh and bulge binding energy. Although the tightness of the Mbh-Ngc correlation suggests an unlikely causal link between supermassive black holes and globular clusters, such a correspondence can exhibit small scatter even if the physical relationship is indirect. We show that the relatively small scatter of the Mbh-Ngc relation owes to the mutual residual correlation of Mbh and Ngc with stellar mass when the velocity dispersion is held fixed. Thus, present observations lend evidence for feedback-regulated models in which the bulge binding energy is most important; they do not necessarily imply any `special' connection between globular clusters and Mbh. This raises the question of why Ngc traces the formation of ellipticals and bulges sufficiently well to be correlated with binding energy.

preprint2011arXiv

Self-Regulated Star Formation in Galaxies via Momentum Input from Massive Stars

Feedback from massive stars is believed to play a critical role in shaping the galaxy mass function, the structure of the interstellar medium (ISM), and the low efficiency of star formation, but the exact form of the feedback is uncertain. In this paper, the first in a series, we present and test a novel numerical implementation of stellar feedback resulting from momentum imparted to the ISM by radiation, supernovae, and stellar winds. We employ a realistic cooling function, and find that a large fraction of the gas cools to <100K, so that the ISM becomes highly inhomogeneous. Despite this, our simulated galaxies reach an approximate steady state, in which gas gravitationally collapses to form giant molecular clouds (GMCs), dense clumps, and stars; subsequently, stellar feedback disperses the GMCs, repopulating the diffuse ISM. This collapse and dispersal cycle is seen in models of SMC-like dwarfs, the Milky-Way, and z~2 clumpy disk analogues. The simulated global star formation efficiencies are consistent with the observed Kennicutt-Schmidt relation. Moreover, the star formation rates are nearly independent of the numerically imposed high-density star formation efficiency, density threshold, and density scaling. This is a consequence of the fact that, in our simulations, star formation is regulated by stellar feedback limiting the amount of very dense gas available for forming stars. In contrast, in simulations without stellar feedback, i.e. under the action of only gravity and gravitationally-induced turbulence, the ISM experiences runaway collapse to very high densities. In these simulations without feedback, the global star formation rates exceed observed galactic star formation rates by 1-2 orders of magnitude, demonstrating that stellar feedback is crucial to the regulation of star formation in galaxies.

preprint2011arXiv

The Origins of AGN Obscuration: The 'Torus' as a Dynamical, Unstable Driver of Accretion

Multi-scale simulations have made it possible to follow gas inflows onto massive black holes (BHs) from galactic scales to the accretion disk. When sufficient gas is driven towards the BH, gravitational instabilities generically form lopsided, eccentric disks that propagate inwards. The lopsided stellar disk exerts a strong torque on the gas disk, driving inflows that fuel rapid BH growth. Here, we investigate whether the same gas disk is the 'torus' invoked to explain obscured AGN. The disk is generically thick and has characteristic ~1-10 pc sizes and masses resembling those required of the torus. The scale heights and obscured fractions of the predicted torii are substantial even in the absence of strong stellar feedback providing the vertical support. Rather, they can be maintained by strong bending modes and warps excited by the inflow-generating instabilities. Other properties commonly attributed to feedback processes may be explained by dynamical effects: misalignment between torus and host galaxy, correlations between local SFR and turbulent gas velocities, and dependence of obscured fractions on AGN luminosity or SFR. We compare the predicted torus properties with observations of gas surface density profiles, kinematics, scale heights, and SFR densities in AGN nuclei, and find that they are consistent. We argue that it is not possible to reproduce these observations and the observed column density (N_H) distribution without a clumpy gas distribution, but allowing for clumping on small scales the predicted N_H distribution is in good agreement with observations from 10^20-27 cm^-2. We examine how N_H scales with galaxy and AGN properties, and find that AGN feedback may be necessary to explain some trends with luminosity and/or redshift. The torus is not merely a bystander or passive fuel source for accretion, but is itself the mechanism driving accretion.

preprint2010arXiv

A Non-Parametric Estimate of Mass 'Scoured' in Galaxy Cores

We present a simple estimate of the mass 'deficits' in cored spheroids, as a function of galaxy mass and radius within the galaxy. Previous attempts to measure such deficits depended on fitting some functional form to the profile at large radii and extrapolating inwards; this is sensitive to the assumed functional form and does not allow for variation in nuclear profile shapes. We take advantage of larger data sets to directly construct stellar mass profiles of observed systems and measure the stellar mass enclosed in a series of physical radii (M(<R)), for samples of cusp and core spheroids at the same stellar mass. There is a significant bimodality in this distribution at small radii, and we non-parametrically measure the median offset between core and cusp populations (the deficit Delta_M(<R)). We construct the scoured mass profile as a function of radius, without reference to any assumed functional form. The mass deficit rises in power-law fashion (Delta_M(<R) R^{1.3-1.8}) from a significant but small mass at R<10pc, to asymptote to a maximum ~0.5-2 M_BH at ~100pc. At larger radii there is no statistically significant separation between populations; the upper limit to the cumulative scoured mass at ~kpc is ~2-4 M_BH. This does not depend strongly on stellar mass. The dispersion in M(<R) appears larger in the core population, possibly reflecting the fact that scouring increases the scatter in profile shapes. These results are in good agreement with models of scouring from BH binary systems.

preprint2010arXiv

A physical model for the origin of the diffuse cosmic infrared background

We present a physical model for origin of the cosmic diffuse infrared background (CDIRB). By utilizing the observed stellar mass function and its evolution as input to a semi-empirical model of galaxy formation, we isolate the physics driving diffuse IR emission. The model includes contributions from three primary sources of IR emission: steady-state star formation owing to isolated disk galaxies, interaction-driven bursts of star formation owing to close encounters and mergers, and obscured active galactic nuclei (AGN). We find that most of the CDIRB is produced by equal contributions from objects at z=0.5-1 and z>1, as suggested by recent observations. Of those sources, the vast majority of the emission originates in systems with low to moderate IR luminosities (L_{IR}<10^{12} $L_sun); the most luminous objects contribute significant flux only at high-redshifts (z>2). All star formation in ongoing mergers accounts for <10% of the total at all wavelengths and redshifts, while emission directly attributable to the interaction-driven burst itself accounts for <5%. We furthermore find that obscured AGN contribute <1-2% of the CDIRB at all wavelengths and redshifts, with a strong upper limit of less than 4% of the total emission. Finally, since electron-positron pair production interactions with the CDIRB represent the primary source of opacity to very high energy (VHE: E_γ> 1 TeV) γ-rays, the model provides predictions for the optical depth of the Universe to the most energetic photons. We find that these predictions agree with observations of high-energy cutoffs at TeV energies in nearby blazars, and suggest that while the Universe is extremely optically thick at >10 TeV, the next generation of VHE γ-ray telescopes can reasonably expect detections from out to 50-150 Mpc.

preprint2010arXiv

Do sub-millimeter galaxy number counts provide evidence for a top-heavy IMF?

Theoretical models have had difficulty matching the observed number density of sub-millimeter galaxies (SMGs), causing some authors (e.g., Baugh et al. 2005) to suggest that SMGs provide evidence for a top-heavy initial mass function (IMF). To test this claim, we have, for the first time, combined high-resolution 3-D hydrodynamic simulations of isolated and merging massive, gas-rich galaxies, radiative transfer, and a semi-empirical merger rate model to predict the number density of SMGs. Our model can reproduce the observed SMG number density even when using a standard (Kroupa) IMF. The agreement is due to a combination of relatively long sub-mm duty cycles for mergers (a few times 10^8 years for our most massive models), which owe to our combination of high-resolution 3-D hydrodynamic simulations and dust radiative transfer; sufficient number densities of massive, gas-rich mergers; and the decrease in sub-mm counts observed by recent deep/wide surveys (e.g., Austermann et al. 2010) relative to previous surveys. Our results suggest that the observed SMG number counts do not provide evidence for a top-heavy IMF at high redshift.

preprint2010arXiv

How Do Massive Black Holes Get Their Gas?

We use multi-scale SPH simulations to follow the inflow of gas from galactic scales to <0.1pc, where the gas begins to resemble a traditional Keplerian accretion disk. The key ingredients are gas, stars, black holes (BHs), self-gravity, star formation, and stellar feedback. We use ~100 simulations to survey a large parameter space of galaxy properties and subgrid models for the ISM physics. We generate initial conditions for our simulations of galactic nuclei (<~300pc) using galaxy scale simulations, including both major mergers and isolated bar-(un)stable disk galaxies. For sufficiently gas-rich, disk-dominated systems, a series of gravitational instabilities generates large accretion rates of up to 1-10 M_sun/yr onto the BH (at <<0.1pc); sufficient to fuel the most luminous quasars. The BH accretion rate is highly time variable, given fixed conditions at ~kpc. At >~10pc, our simulations resemble the 'bars within bars' model, but the gas exhibits diverse morphologies, including spirals, rings, clumps, and bars; their duty cycle is modest, complicating attempts to correlate BH accretion with nuclear morphology. At ~1-10pc, the gravitational potential becomes dominated by the BH and bar-like modes are no longer present. However, the gas becomes unstable to a standing, eccentric disk or a single-armed spiral mode (m=1), driving the gas to sub-pc scales. Proper treatment of this mode requires including star formation and the self-gravity of both the stars and gas. We predict correlations between BHAR and SFR at different galactic nuclei: nuclear SF is more tightly coupled to AGN activity, but correlations exist at all scales.

preprint2010arXiv

Mergers and Bulge Formation in Lambda-CDM: Which Mergers Matter?

We use a suite of semi-empirical models to predict galaxy merger rates and contributions to bulge growth as functions of merger mass, redshift, and mass ratio. The models use empirical halo occupation constraints to identify mergers, together with high-resolution simulations to quantify how mergers with different properties contribute to the bulge population. We find good agreement with a variety of observational constraints, and provide fitting functions for merger rates and contributions to bulge growth. We identify several robust conclusions. (1) Major mergers dominate formation and assembly of L* bulges and the spheroid mass density, minor mergers contribute ~30%. (2) This is mass-dependent: bulge formation is dominated by more minor mergers in lower-mass systems. At higher masses, bulges form in major mergers near L*, but subsequently assemble in minor mergers. (3) The minor/major contribution is also morphology-dependent: higher B/T systems form in more major mergers, lower B/T systems form in situ from minor mergers. (4) Low-mass galaxies, being gas-rich, require more major mergers to reach the same B/T as high-mass systems. (5) Absolute merger rates increase with galaxy mass. (6) Predicted rates agree well with observations, but suggest that some morphology-selected samples include contamination from minor mergers. (7) Predicted rates agree with integrated growth in bulge mass with cosmic time, but with factor ~2 uncertainty - half the bulge mass density could come from non-mergers. We consider ~1000 model variations and quantify resulting uncertainties. Conclusions regarding the major/minor contribution to bulge growth are very robust, absolute merger rates have systematic factor ~2 uncertainties.

preprint2010arXiv

Mergers in Lambda-CDM: Uncertainties in Theoretical Predictions and Interpretations of the Merger Rate

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

preprint2010arXiv

On Sizes, Kinematics, M/L Gradients, and Light Profiles of Massive Compact Galaxies at z~2

We present a detailed analysis of the structure and resolved stellar populations of simulated merger remnants, and compare them to observations of compact quiescent galaxies at z ~ 2. We find that major merging is a viable mechanism to produce systems of ~ 10^11 Msun and ~ 1 kpc size, provided the gas fraction at the time of final coalescence is high (~ 40%), and provided that the progenitors are compact star-forming galaxies, as expected at high redshift. Their integrated spectral energy distributions and velocity dispersions are in good agreement with the observations, and their position in the (v_{maj}/sigma, ellipticity) diagram traces the upper envelope of the distribution of lower redshift early-type galaxies. The simulated merger remnants show time- and sightline-dependent M/L ratio gradients that result from a superposition of radially dependent stellar age, stellar metallicity, and extinction. The median ratio of effective radius in rest-frame V-band light to that in mass surface density is ~ 2 during the quiescent remnant phase. This is typically expressed by a negative color gradient (i.e., red core), which we expect to correlate with the integrated color of the system. Finally, the simulations differ from the observations in their surface brightness profile shape. The simulated remnants are typically best fit by high (n >> 4) Sersic indices, whereas observed quiescent galaxies at z ~ 2 tend to be less cuspy (median n ~ 2.3). Limiting early star formation in the progenitors may be required to prevent the simulated merger remnants from having extended wings.

preprint2010arXiv

Origins, Structure, and Inflows of m=1 Modes in Quasi-Keplerian Disks

Simulations show eccentric disks (m=1 modes) forming around quasi-Keplerian potentials, a topic of interest for fueling quasars, forming super-massive BHs, planet formation and migration, explaining the origin and properties of nuclear eccentric stellar disks like that in M31, and driving the formation of the obscuring AGN torus. We consider the global, linear normal m=1 modes in collisionless disks, without the restriction that the disk mass be negligible relative to the central (Keplerian) mass. We derive their structure and key resonance features, and show how they arise, propagate inwards, and drive both inflow/outflow and eccentricities in the disk. We compare with hydrodynamic simulations of such disks around a super-massive BH, with star formation, gas cooling, and feedback. We derive the dependence of the normal mode structure on disk structure, mass profiles, and thickness, and mode pattern speeds and growth rates. We show that, if the disk at some radii has mass of >~10% the central point mass, the modes are linearly unstable and are self-generating. They arise as 'fast modes' with pattern speed of order the local angular velocity at these radii. The characteristic global normal modes have pattern speeds comparable to the linear growth rate, of order (G*M_0*R_0^{-3})^{1/2}, where M_0 is the central mass and R_{0} is the radius where the enclosed disk mass ~M_{0}. They propagate inwards by exciting eccentricities towards smaller and smaller radii, until at small radii these are 'slow modes.' With moderate amplitude, the global normal modes can lead to shocks and significant gas inflows at near-Eddington rates at all radii inside several ~R_0.

preprint2010arXiv

The Nuclear Stellar Disk in Andromeda: A Fossil from the Era of Black Hole Growth

The physics of angular momentum transport from galactic scales (~10-100 pc) to much smaller radii is one of the oustanding problems in our understanding of the formation and evolution of super-massive black holes (BHs). Seemingly unrelated observations have discovered that there is a lopsided stellar disk of unknown origin orbiting the BH in M31, and possibly many other systems. We show that these nominally independent puzzles are in fact closely related. Multi-scale simulations of gas inflow from galactic to BH scales show that when sufficient gas is driven towards a BH, gravitational instabilities form a lopsided, eccentric disk that propagates inwards from larger radii. The lopsided stellar disk exerts a strong torque on the remaining gas, driving inflows that fuel the growth of the BH and produce quasar-level luminosities. The same disk can produce significant obscuration along many sightlines and thus may be the putative 'torus' invoked to explain obscured active galactic nuclei and the cosmic X-ray background. The stellar relic of this disk is long lived and retains the eccentric pattern. Simulations that yield quasar-level accretion rates produce relic stellar disks with kinematics, eccentric patterns, precession rates, and surface density profiles in reasonable agreement with observations of M31. The observed properties of nuclear stellar disks can thus be used to constrain the formation history of super-massive BHs.

preprint2009arXiv

A Maximum Stellar Surface Density in Dense Stellar Systems

We compile observations of the surface mass density profiles of dense stellar systems, including globular clusters in the Milky Way and nearby galaxies, massive star clusters in nearby starbursts, nuclear star clusters in dwarf spheroidals and late-type disks, ultra-compact dwarfs, and galaxy spheroids spanning the range from low-mass cusp bulges and ellipticals to massive core ellipticals. We show that in all cases the maximum stellar surface density attained in the central regions of these systems is similar, Sigma_max ~ 10^11 M_sun/kpc^2 (~20 g/cm^2), despite the fact that the systems span 7 orders of magnitude in total stellar mass M_star, 5 in effective radius R_e, and have a wide range in effective surface density M_star/R_e^2. The surface density limit is reached on a wide variety of physical scales in different systems and is thus not a limit on three-dimensional stellar density. Given the very different formation mechanisms involved in these different classes of objects, we argue that a single piece of physics likely determines Sigma_max. The radiation fields and winds produced by massive stars can have a significant influence on the formation of both star clusters and galaxies, while neither supernovae nor black hole accretion are important in star cluster formation. We thus conclude that feedback from massive stars likely accounts for the observed Sigma_max, plausibly because star formation reaches an Eddington-like flux that regulates the growth of these diverse systems. This suggests that current models of galaxy formation, which focus on feedback from supernovae and active galactic nuclei, are missing a crucial ingredient.

preprint2009arXiv

A New Empirical Method to Infer the Starburst History of the Universe from Local Galaxy Properties

The centers of bulges are formed dissipationally via gas inflows over short timescales: the 'starburst' mode of star formation (SF). Recent work has shown that detailed observations can be used to separate the stellar mass profile of these 'burst relic' components in local systems. Together with the assumption that some Kennicutt-Schmidt law holds, and that the burst was indeed a dissipational gas-rich event, we show that the observed profiles can be inverted to obtain the time and space-dependent SF history of each burst. Performing this with a large sample of well-studied spheroids, we show that the implied bursts scale in magnitude, mass, peak SFR, and spatial extent with galaxy mass in simple manner, and provide fits to these correlations. Burst masses are ~10% the total spheroid mass; timescales a mass-independent ~10^8 yr; peak SFR ~M_burst/t_burst; and they decay in power-law fashion ~t^-2.4. Sizes are ~0.1 R_e(spheroid), but grow with time. Combined with measurements of the nuclear stellar population ages of these systems (i.e. burst times), it is possible to re-construct the burst contribution to the distribution of SFRs and IR luminosity functions at all redshifts. The burst LF agrees well with observed IR LFs at the brightest luminosities, at z=0-2. At low-L, bursts are always unimportant; the transition to their importance increases from ULIRG luminosities at z~0 to HyLIRG luminosities at z~2. At all redshifts, bursts are a small fraction (~5-10%) of the total SFR density. We discuss possible implications of tension between maximum relic stellar mass densities in massive systems, and estimated number counts of brightest sub-millimeter galaxies.

preprint2009arXiv

Are Most Low-Luminosity AGN Really Obscured?

At low Eddington ratio (mdot), two effects make it harder to detect AGN given some selection criteria. First, even with fixed accretion physics, AGN are diluted/less luminous relative to their hosts; the magnitude of this depends on host properties and so on luminosity and redshift. Second, they may transition to a radiatively inefficient state, changing SED shape and dramatically decreasing in optical/IR luminosity. These effects lead to differences in observed AGN samples, even at fixed bolometric luminosity and after correction for obscuration. The true Eddington ratio distribution may depend strongly on luminosity, but this will be seen only in surveys robust to dilution and radiative inefficiency (X-ray or narrow-line samples); selection effects imply that AGN in optical samples will have uniformly high mdot. This also implies that different selection methods yield systems with different hosts: the clustering of faint optical/IR sources will be weaker than that of X-ray sources, and optical/IR Seyferts will reside in more disk-dominated galaxies while X-ray selected Seyferts will preferentially occupy early-type systems. If observed mdot distributions are correct, a large fraction of low-luminosity AGN currently classified as 'obscured' are in fact diluted and/or radiatively inefficient, not obscured by gas or dust. This is equally true if X-ray hardness is used as a proxy for obscuration, since radiatively inefficient SEDs near mdot~0.01 are X-ray hard. These effects can explain most of the claimed luminosity/redshift dependence in the 'obscured' AGN population, with the true obscured fraction as low as 20%.

preprint2009arXiv

Compact High-Redshift Galaxies Are the Cores of the Most Massive Present-Day Spheroids

Observations suggest that effective radii of high-z massive spheroids are as much as a factor ~6 smaller than low-z galaxies of comparable mass. Given the apparent absence of low-z counterparts, this has often been interpreted as indicating that the high density, compact red galaxies must be 'puffed up' by some mechanism. We compare the ensemble of high-z observations with large samples of well-observed low-z ellipticals. At the same physical radii, the stellar surface mass densities of low and high-z systems are comparable. Moreover, the abundance of high surface density material at low redshift is comparable to or larger than that observed at z>1-2, consistent with the continuous buildup of spheroids over this time. The entire population of compact, high-z red galaxies may be the progenitors of the high-density cores of present-day ellipticals, with no need for a decrease in stellar density from z=2 to z=0. The primary difference between low and high-z systems is thus the observed low-density material at large radii in low-z spheroids (rather than the high-density material in high-z spheroids). Such low-density material may either (1) assemble at z<2 or (2) be present, but not yet detected, at z>2. Mock observations of low-z massive systems show that the high-z observations do not yet probe sufficiently low surface brightness material to detect the low surface density 'wings' (if present). Thus, if the high-z galaxies resemble the most massive systems today, their inferred effective radii could be under-estimated by factors ~2-4. This difference arises because massive systems at low redshift are not well-fit by single Sersic profiles. We discuss implications of our results for physical models of galaxy evolution.

preprint2009arXiv

Discriminating Between the Physical Processes that Drive Spheroid Size Evolution

Massive galaxies at high-z have smaller effective radii than those today, but similar central densities. Their size growth therefore relates primarily to the evolving abundance of low-density material. Various models have been proposed to explain this evolution, which have different implications for galaxy, star, and BH formation. We compile observations of spheroid properties as a function of redshift and use them to test proposed models. Evolution in progenitor gas-richness with redshift gives rise to initial formation of smaller spheroids at high-z. These systems can then evolve in apparent or physical size via several channels: (1) equal-density 'dry' mergers, (2) later major or minor 'dry' mergers with less-dense galaxies, (3) adiabatic expansion, (4) evolution in stellar populations & mass-to-light-ratio gradients, (5) age-dependent bias in stellar mass estimators, (6) observational fitting/selection effects. If any one of these is tuned to explain observed size evolution, they make distinct predictions for evolution in other galaxy properties. Only model (2) is consistent with observations as a dominant effect. It is the only model which allows for an increase in M_BH/M_bulge with redshift. Still, the amount of merging needed is larger than that observed or predicted. We therefore compare cosmologically motivated simulations, in which all these effects occur, & show they are consistent with all the observational constraints. Effect (2), which builds up an extended low-density envelope, does dominate the evolution, but effects 1,3,4, & 6 each contribute ~20% to the size evolution (a net factor ~2). This naturally also predicts evolution in M_BH-sigma similar to that observed.

preprint2009arXiv

Mergers, AGN, and 'Normal' Galaxies: Contributions to the Distribution of Star Formation Rates and Infrared Luminosity Functions

We use a novel method to predict the contribution of normal star-forming galaxies, merger-induced bursts, and obscured AGN, to IR luminosity functions (LFs) and global SFR densities. We use empirical halo occupation constraints to populate halos with galaxies and determine the distribution of normal and merging galaxies. Each system can then be associated with high-resolution hydrodynamic simulations. We predict the distribution of observed luminosities and SFRs, from different galaxy classes, as a function of redshift from z=0-6. We provide fitting functions for the predicted LFs, quantify the uncertainties, and compare with observations. At all redshifts, 'normal' galaxies dominate the LF at moderate luminosities ~L* (the 'knee'). Merger-induced bursts increasingly dominate at L>>L*; at the most extreme luminosities, AGN are important. However, all populations increase in luminosity at higher redshifts, owing to increasing gas fractions. Thus the 'transition' between normal and merger-dominated sources increases from the LIRG-ULIRG threshold at z~0 to bright Hyper-LIRG thresholds at z~2. The transition to dominance by obscured AGN evolves similarly, at factor of several higher L_IR. At all redshifts, non-merging systems dominate the total luminosity/SFR density, with merger-induced bursts constituting ~5-10% and AGN ~1-5%. Bursts contribute little to scatter in the SFR-stellar mass relation. In fact, many systems identified as 'ongoing' mergers will be forming stars in their 'normal' (non-burst) mode. Counting this as 'merger-induced' star formation leads to a stronger apparent redshift evolution in the contribution of mergers to the SFR density.

preprint2009arXiv

Metallicity Gradients at Large Galactocentric Radii Using the Near-infrared Calcium Triplet

We describe a new spectroscopic technique for measuring radial metallicity gradients out to large galactocentric radii. We use the DEIMOS multi-object spectrograph on the Keck telescope and the galaxy spectrum extraction technique of Proctor et al. (2009). We also make use of the metallicity sensitive near-infrared (NIR) Calcium triplet (CaT) features together with single stellar population models to obtain metallicities. Our technique is applied as a pilot study to a sample of three relatively nearby (<30 Mpc) intermediate-mass to massive early-type galaxies. Results are compared with previous literature inner region values and generally show good agreement. We also include a comparison with profiles from dissipational disk-disk major merger simulations. Based on our new extended metallicity gradients combined with other observational evidence and theoretical predictions, we discuss possible formation scenarios for the galaxies in our sample. The limitations of our new technique are also discussed.

preprint2009arXiv

Quasar Feedback: More Bang for Your Buck

We propose a two-stage model for the effects of feedback from a bright quasar on the cold gas in a galaxy. It is difficult for feedback from near the accretion disk to directly impact dense molecular clouds at ~kpc. But if such feedback can drive a weak wind or outflow in the hot, diffuse ISM (a relatively 'easy' task), then in the wake of such an outflow passing over a cold cloud, a combination of instabilities will drive the cloud material to effectively expand in the direction perpendicular to the outflow. Such expansion dramatically increases the effective cross section of the cloud material and makes it more susceptible to ionization and momentum coupling from absorption of the incident quasar radiation field. Even a moderate effect of this nature can dramatically alter the ability of clouds at large radii to be fully ionized and driven into a secondary outflow by radiation pressure. Since the amount of momentum and volume which can be ionized by observed quasar radiation field is more than sufficient to affect the entire cold gas supply once it has been altered in this manner (and the 'initial' feedback need only initiate a moderate wind in the low-density hot gas), this reduces by an order of magnitude the required energy budget for feedback to affect a host galaxy. Instead of ~5% of the radiated energy (~100% momentum) needed if the initial feedback must directly heat or blow out the galactic gas, if only ~0.5% of the luminosity (~10% momentum) can couple to drive the initial hot outflow, this mechanism could be efficient. This amounts to hot gas outflow rates from near the accretion disk of only 5-10% of the BH accretion rate.

preprint2009arXiv

Recovering Stellar Population Properties and Redshifts from Broad-Band Photometry of Simulated Galaxies: Lessons for SED Modeling

We present a detailed analysis of our ability to determine stellar masses, ages, reddening and extinction values, and star formation rates of high-redshift galaxies by modeling broad-band SEDs with stellar population synthesis. In order to do so, we computed synthetic optical-to-NIR SEDs for model galaxies taken from hydrodynamical merger simulations placed at redshifts 1.5 < z < 3. Viewed under different angles and during different evolutionary phases, the simulations represent a wide variety of galaxy types (disks, mergers, spheroids). We show that simulated galaxies span a wide range in SEDs and color, comparable to these of observed galaxies. In all star-forming phases, dust attenuation has a large effect on colors, SEDs, and fluxes. The broad-band SEDs were then fed to a standard SED modeling procedure and resulting stellar population parameters were compared to their true values. Disk galaxies generally show a decent median correspondence between the true and estimated mass and age, but suffer from large uncertainties. During the merger itself, we find larger offsets (e.g., log M_recovered - log M_true = -0.13^{+0.10}_{-0.14}). E(B-V) values are generally recovered well, but the estimated total visual absorption Av is consistently too low, increasingly so for larger optical depths. Since the largest optical depths occur during the phases of most intense star formation, it is for the highest SFRs that we find the largest underestimates. The masses, ages, E(B-V), Av, and SFR of merger remnants (spheroids) are very well reproduced. We discuss possible biases in SED modeling results caused by mismatch between the true and template star formation history, dust distribution, metallicity variations and AGN contribution.

preprint2009arXiv

The Effects of Gas on Morphological Transformation in Mergers: Implications for Bulge and Disk Demographics

Transformation of disks into spheroids via mergers is a well-accepted element of galaxy formation models. However, recent simulations have shown that bulge formation is suppressed in increasingly gas-rich mergers. We investigate the global implications of these results in a cosmological framework, using independent approaches: empirical halo-occupation models (where galaxies are populated in halos according to observations) and semi-analytic models. In both, ignoring the effects of gas in mergers leads to the over-production of spheroids: low and intermediate-mass galaxies are predicted to be bulge-dominated (B/T~0.5 at <10^10 M_sun), with almost no bulgeless systems), even if they have avoided major mergers. Including the different physical behavior of gas in mergers immediately leads to a dramatic change: bulge formation is suppressed in low-mass galaxies, observed to be gas-rich (giving B/T~0.1 at <10^10 M_sun, with a number of bulgeless galaxies in good agreement with observations). Simulations and analytic models which neglect the similarity-breaking behavior of gas have difficulty reproducing the strong observed morphology-mass relation. However, the observed dependence of gas fractions on mass, combined with suppression of bulge formation in gas-rich mergers, naturally leads to the observed trends. Discrepancies between observations and models that ignore the role of gas increase with redshift; in models that treat gas properly, galaxies are predicted to be less bulge-dominated at high redshifts, in agreement with the observations. We discuss implications for the global bulge mass density and future observational tests.

preprint2009arXiv

The Small Scatter in BH-Host Correlations & The Case for Self-Regulated BH Growth

Supermassive black holes (BHs) obey tight scaling relations between their mass and their host galaxy properties such as total stellar mass, velocity dispersion, and potential well depth. This has led to the development of self-regulated models for BH growth, in which feedback from the central BH halts its own growth upon reaching a critical threshold. However, models have also been proposed in which feedback plays no role: so long as a fixed fraction of the host gas supply is accreted, relations like those observed can be reproduced. Here, we argue that the scatter in the observed BH-host correlations, and its run with scale, presents a demanding constraint on any model for these correlations, and that it favors self-regulated models of BH growth. We show that the scatter in the stellar mass fraction within a radius R in observed ellipticals and spheroids increases strongly at small R. At fixed total stellar mass (or host velocity dispersion), on very small scales near the BH radius of influence, there is an order-of-magnitude scatter in the amount of gas that must have entered and formed stars. In short, the BH appears to 'know more' about the global host galaxy potential on large scales than the stars and gas supply on small scales. This is predicted in self-regulated models; however, models where there is no feedback would generically predict order-of-magnitude scatter in the BH-host correlations. Likewise, models in which the BH feedback in the 'bright' mode does not regulate the growth of the BH itself, but sets the stellar mass of the galaxy by inducing star formation or blowing out a mass in gas much larger than the galaxy stellar mass, are difficult to reconcile with the scatter on small scales.

preprint2009arXiv

When Should We Treat Galaxies as Isolated?

Traditionally, secular evolution is defined as evolution of systems where the internal growth of structure and instabilities dominates the growth via external drivers (e.g. accretion/mergers). Most study has focused on 'isolated' galaxies, where seed asymmetries may represent realistic cosmological substructure, but subsequent evolution ignores galaxy growth. Large-scale modes in the disk then grow on a timescale of order a disk rotation period (0.1-1 Gyr). If, however, galaxies evolve cosmologically on a shorter timescale, then it may not be appropriate to consider them 'isolated.' We outline simple scalings to ask whether the timescale for secular evolution is shorter than the timescale for cosmological accretion and mergers. This is the case in a narrow, but important range of perturbation amplitudes corresponding to substructure or mode/bar fractional amplitudes 0.01-0.1, a range of interest for observed strong bars and pseudobulges. At smaller amplitudes <<0.1, systems are not isolated: typical disks will grow by accretion at a comparable level over even a single dynamical time. At larger amplitudes >>0.1, the evolution is no longer secular; direct gravitational evolution of the seed swamps the internal disk response. We derive criteria for when disks can be well-approximated as 'isolated' as a function of mass, redshift, and disk stability. The relevant parameter space shrinks at higher mass, higher disk stability, and higher-z as accretion rates increase. Cosmological rates of galaxy evolution also define a maximum bar/mode lifetime of practical interest, of ~0.1/H(z). Longer-lived modes will de-couple from their drivers (if driven) and encounter cosmological effects.

preprint2008arXiv

A Characteristic Division Between the Fueling of Quasars and Seyferts: Five Simple Tests

Given the existence of the M_BH-sigma relation, models of self-regulated black hole (BH) growth require both a fuel supply and growth of the host bulge to deepen the potential, or else the system will either starve or self-regulate without sustained activity. This suggests that bright quasars must be triggered in major mergers: a large fraction of the galaxy must be converted to new bulge mass in a dynamical time or less. Low-luminosity AGN, in contrast, require little bulge growth and small gas supplies, and could be triggered in more common non-merger events. This predicts a transition to merger-induced fueling around the traditional quasar-Seyfert luminosity divide (growth of BH masses above/below 10^7 M_sun). We compile observations to test several predictions of such a division, including: (1) A transition to bulge-dominated hosts. (2) A transition between 'pseudobulges' and 'classical' bulges hosting the remnant BHs: pseudobulges are formed in secular processes and minor mergers, whereas classical bulges are relics of major mergers. (3) An increase in the amplitude of small-scale clustering where mergers are more efficient. (4) Different redshift evolution, with gas-rich merger rates rising to redshifts z>2 while secular processes are relatively constant in time. (5) An increasing prominence of post-starburst features in more luminous systems. Our compilation of observations provides tentative evidence for the predicted division around the Seyfert-quasar threshold. We discuss how future observations can improve these constraints and break degeneracies between different fueling models.

preprint2008arXiv

Dissipation and Extra Light in Galactic Nuclei: IV. Evolution in the Scaling Relations of Spheroids

We develop a model for the origins and redshift evolution of spheroid scaling relations. We consider spheroid sizes, velocity dispersions, masses, profile shapes (Sersic indices), and black hole (BH) masses, and their related scalings. Our approach combines advantages of observational constraints in halo occupation models and hydrodynamic merger simulations. This allows us to separate the relative roles of dissipation, dry mergers, formation time, and progenitor evolution, and identify their effects on scalings at each redshift. Dissipation is the most important factor determining spheroid sizes and fundamental plane (FP) scalings, and can account for the FP tilt and differences between disk and spheroid scalings. Because disks at high-z have higher gas fractions, mergers are more gas-rich, yielding more compact spheroids. This predicts mass-dependent evolution in spheroid sizes, in agreement with observations. This relates to subtle evolution in the FP, important to studies that assume a fixed intrinsic FP. This also predicts mild evolution in BH-host correlations, towards larger BHs at higher z. Dry mergers are significant, but only for massive systems which form early: they form compact, but undergo dry mergers (consistent with observations) such that their sizes at later times are similar to spheroids of similar mass formed more recently. We model descendants of observed compact high-z spheroids: most will become cores of BCGs, with sizes, velocity dispersions, and BH masses consistent with observations, but we identify a fraction that might survive to z=0 intact.

preprint2008arXiv

How Do Disks Survive Mergers?

We develop a physical model for how galactic disks survive and/or are destroyed in interactions. Based on dynamical arguments, we show gas primarily loses angular momentum to internal torques in a merger. Gas within some characteristic radius (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 disk, fall to the center and be consumed in a starburst. A similar analysis predicts where violent relaxation of the stellar disks is efficient. Our model allows us to predict the stellar and gas content that will survive to re-form a disk in the remnant, versus being violently relaxed or contributing to a starburst. We test this in hydrodynamic simulations and find good agreement as a function of mass ratio, orbital parameters, and gas fraction, in simulations spanning a wide range in these properties and others, including different prescriptions for gas physics and feedback. In an immediate sense, the amount of disk that re-forms can be understood in terms of well-understood gravitational physics, independent of details of ISM gas physics or feedback. This allows us to explicitly quantify the requirements for such feedback to (indirectly) enable disk survival, by changing the pre-merger gas content and distribution. The efficiency of disk destruction is a strong function of gas content: we show how and why sufficiently gas-rich major mergers can, under general conditions, yield systems with small bulges (B/T<0.2). We provide prescriptions for inclusion of our results in semi-analytic models.

preprint2007arXiv

Observational Evidence for the Co-evolution of Galaxy Mergers, Quasars, and the Blue/Red Galaxy Transition

We compile a number of observations to estimate the time-averaged rate of formation or buildup of red sequence galaxies, as a function of mass and redshift. Comparing this with the mass functions of mergers and quasar hosts, and independently comparing their clustering properties as a function of redshift, we find that these populations trace the same mass distribution, with similar evolution, at redshifts 0<z<~1.5. Knowing one of the quasar, merger, or elliptical mass/luminosity functions, it is possible to predict the others. Allowing for greater model dependence, we compare the rate of early-type buildup with the implied merger and quasar triggering rates as a function of mass and redshift and find agreement. Over this redshift range, observed merger fractions can account for the entire bright quasar luminosity function and buildup of the red sequence at all but the highest masses at low redshift (>~10^11 M_solar at z<~0.3) where 'dry' mergers appear to dominate. This supports a necessary prediction of theories where mergers between gas-rich galaxies produce ellipticals with an associated phase of quasar activity, after which the remnant becomes red. These populations trace a similar characteristic transition mass, possibly reflecting the mass above which the elliptical population is mostly (>~50%) assembled at a given redshift, which increases with redshift over the observed range in a manner consistent with suggestions that cosmic downsizing may apply to red galaxy assembly as well as star formation. These mass distributions as a function of redshift do not uniformly trace the all/red/blue galaxy population, ruling out models in which quasar activity is generically associated with star formation or is long lived in 'old' systems.

preprint2007arXiv

Photometric Properties of the Most Massive High-Redshift Galaxies

We calculate the observable properties of the most massive high-redshift galaxies in the hierarchical formation scenario where stellar spheroid and supermassive black hole growth are fueled by gas-rich mergers. Combining high-resolution hydrodynamical simulations of the hierarchical formation of a z~6 quasar, stellar population synthesis models, template AGN spectra, prescriptions for interstellar and intergalactic absorption, and the response of modern telescopes, the photometric evolution of galaxies destined to host z~6 quasars are modeled at redshifts z~4-14. These massive galaxies, with enormous stellar masses of M_star ~10^11.5-10^12 M_sun. and star formation rates of SFR~10^3-10^4 M_sun yr^-1 at z>~7, satisfy a variety of photometric selection criteria based on Lyman-break techniques including V-band dropouts at z>~5, i-band dropouts at z>~6, and z-band dropouts at z>~7. The observability of the most massive high-redshift galaxies is assessed and compared with a wide range of existing and future photometric surveys including SDSS, GOODS/HUDF, NOAO WDFS, UKIDSS, the IRAC Shallow Survey, Pan-STARRS, LSST, and SNAP. Massive stellar spheroids descended from z~6 quasars will likely be detected at z~4 by existing surveys, but owing to their low number densities the discovery of quasar progenitor galaxies at z>7 will likely require future surveys of large portions of the sky (>~0.5%) at wavelengths lambda>1 micron. The detection of rare, star-bursting, massive galaxies at redshifts z>~6 would provide support for the hierarchical formation of the earliest quasars and characterize the primitive star-formation histories of the most luminous elliptical galaxies.

preprint2006arXiv

An Upper Limit to the Degree of Evolution Between Supermassive Black Holes and their Host Galaxies

We describe a model-independent integral constraint which defines an upper limit to the allowed degree of evolution in the ratio of black hole (BH) mass to host galaxy luminosity or mass, as a function of redshift. Essentially, if the BH/host ratio is excessive at redshift z, then it would imply that the total mass density in BHs above some M_BH(min) is larger at that redshift than at z=0, which is impossible. This argument requires no knowledge of host or BH properties, only a lower limit to the observed luminosity density in the brightest galaxies at some z. We calculate this upper limit from a variety of luminosity and mass functions in different bands from redshifts z=0-2. We show that it is consistent with passive evolution of spheroid populations (with a fixed M_BH/M_host relation) in all cases, and provides tighter constraints than have generally been obtained previously, ruling out at >6 sigma observational and theoretical estimates suggesting that M_BH/M_host was significantly larger at high redshifts than locally, although relatively weak (factor ~2 by z=2) evolution is still allowed. We discuss a variety of possible 'loopholes' or changes in the BH/host populations and correlations, and show that they typically lower the upper limits and strengthen our conclusions.

preprint2006arXiv

The Relation Between Quasar and Merging Galaxy Luminosity Functions and the Merger-Induced Star Formation Rate of the Universe

Using a model for self-regulated growth of black holes (BHs) in mergers involving gas-rich galaxies, we study the relationship between quasars and the population of merging galaxies and predict the merger-induced star formation rate density of the Universe. Mergers drive nuclear gas inflows, fueling starbursts and 'buried quasars' until accretion feedback expels the gas, rendering a briefly visible optical quasar. Star formation is shut down and accretion declines, leaving a passively evolving remnant with properties typical of red, elliptical galaxies. Based on evolution of these events in our simulations, we demonstrate that the observed statistics of merger rates, luminosity functions (LFs) and mass functions, SFR distributions, specific SFRs, quasar and quasar host galaxy LFs, and elliptical/red galaxy LFs are self-consistent and follow from one another as predicted by the merger hypothesis. We use our simulations to de-convolve both quasar and merging galaxy LFs to determine the birthrate of black holes of a given final mass and merger rates as a function of stellar mass. We use this to predict the merging galaxy LF in several observed wavebands, color-magnitude relations, mass functions, absolute and specific SFR distributions and SFR density, and quasar host galaxy LFs, as a function of redshift from z=0-6. We invert this and predict e.g. quasar LFs from observed merger LFs or SFR distributions. Our results agree well with observations, but idealized models of quasar lightcurves are ruled out by comparison of merger and quasar observations at >99.9% confidence. Using only observations of quasars, we estimate the contribution of mergers to the SFR density of the Universe even to high redshifts z~4.

preprint2005arXiv

The Evolution of the M-sigma Relation

(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.

Philip F. Hopkins

What is connected

Connect this record

See the researcher in context

Building this map preview

140 published item(s)

Early Stages of Dusty Tori: The First Infrared Spectra from a Highly Multiscale Quasar Simulation

A Consistent Reduced-Speed-of-Light Formulation of Cosmic Ray Transport Valid in Weak and Strong-Scattering Regimes

Amplified J-factors in the Galactic Center for velocity-dependent darkmatter annihilation in FIRE simulations

Black hole -- galaxy scaling relations in FIRE: the importance of black hole location and mergers

Cluster assembly and the origin of mass segregation in the STARFORGE simulations

Cosmic-Ray Driven Outflows to Mpc Scales from $L_{\ast}$ Galaxies

Dissipative Dark Matter on FIRE: II. Observational signatures and constraints from local dwarf galaxies

Effects of Different Cosmic Ray Transport Models on Galaxy Formation

Effects of the environment and feedback physics on the initial mass function of stars in the STARFORGE simulations

Evidence for a Vast Prograde Stellar Stream in the Solar Vicinity

Exploring metallicity-dependent rates of Type Ia supernovae and their impact on galaxy formation

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

Gas infall and radial transport in cosmological simulations of Milky Way-mass disks

Hot-mode accretion and the physics of thin-disk galaxy formation

Numerical Study of Cosmic Ray Confinement through Dust Resonant Drag Instabilities

Rapid disc settling and the transition from bursty to steady star formation in Milky Way-mass galaxies

Shapes of Milky-Way-Mass Galaxies with Self-Interacting Dark Matter

Simulating the Diverse Instabilities of Dust in Magnetized Gas

Spatially Resolved Gas-phase Metallicity in FIRE-2 Dwarfs: Late-Time Evolution of Metallicity Relations in Simulations with Feedback and Mergers

Testing Physical Models for Cosmic Ray Transport Coefficients on Galactic Scales: Self-Confinement and Extrinsic Turbulence at GeV Energies

The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert

The Galactic Dust-Up: Modeling Dust Evolution in FIRE

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

The observability of galaxy merger signatures in nearby gas-rich spirals

Why do Black Holes Trace Bulges (& Central Surface Densities), Instead of Galaxies as a Whole?

Fiery Cores: Bursty and Smooth Star Formation Distributions across Galaxy Centers in Cosmological Zoom-in Simulations

On the optical properties of Resonant Drag Instabilities: Variability of Asymptotic Giant Branch and R Coronae Borealis stars

Probing Hot Gas Components of Circumgalactic Medium in Cosmological Simulations with the Thermal Sunyaev-Zel'dovich Effect

STARFORGE: The effects of protostellar outflows on the IMF

A Census of Sub-kiloparsec Resolution Metallicity Gradients in Star-forming Galaxies at Cosmic Noon from HST Slitless Spectroscopy

A dark matter profile to model diverse feedback-induced core sizes of $Λ$CDM haloes

A general-purpose timestep criterion for simulations with gravity

A profile in FIRE: resolving the radial distributions of satellite galaxies in the Local Group with simulations

A Relationship Between Stellar Metallicity Gradients and Galaxy Age in Dwarf Galaxies

A Stable Finite-Volume Method for Scalar-Field Dark Matter

But What About... Cosmic Rays, Magnetic Fields, Conduction, & Viscosity in Galaxy Formation

Cataloging Accreted Stars within Gaia DR2 using Deep Learning

Measuring dynamical masses from gas kinematics in simulated high-redshift galaxies

Most Stars (and Planets?) Are Born in Intense Radiation Fields

No missing photons for reionization: moderate ionizing photon escape fractions from the FIRE-2 simulations

Physical models of streaming instabilities in protoplanetary disks

Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos

Radiative Stellar Feedback in Galaxy Formation: Methods and Physics

Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies

Synthetic Gaia surveys from the FIRE cosmological simulations of Milky Way-mass galaxies

The Bolometric Quasar Luminosity Function at z = 0-7

The IRX-$β$ relation of high-redshift galaxies

The Robustness of Cosmological Hydrodynamic Simulation Predictions to Changes in Numerics and Cooling Physics

The universal acceleration scale from stellar feedback

High-redshift Galaxy Formation with Self-consistently Modeled Stars and Massive Black Holes: Stellar Feedback and Quasar Growth

On the Survival of Cool Clouds in the Circum-Galactic Medium

Stars made in outflows may populate the stellar halo of the Milky Way

Swirls of FIRE: Spatially Resolved Gas Velocity Dispersions and Star Formation Rates in FIRE-2 Disk Environments

Cosmic Rays or Turbulence can Suppress Cooling Flows (Where Thermal Heating or Momentum Injection Fail)

Binary Stars Can Provide the "Missing Photons" Needed for Reionization

Breathing FIRE: How Stellar Feedback Drives Radial Migration, Rapid Size Fluctuations, and Population Gradients in Low-Mass Galaxies

Evolution of Intrinsic Scatter in the SFR-Stellar Mass Correlation at 0.5<z<3

Giant clumps in the FIRE simulations: a case study of a massive high-redshift galaxy

Gravitational torque-driven black hole growth and feedback in cosmological simulations

MUFASA: Galaxy Formation Simulations With Meshless Hydrodynamics

Reconciling dwarf galaxies with LCDM cosmology: Simulating a realistic population of satellites around a Milky Way-mass galaxy

Star Formation in a Turbulent Framework: From Giant Molecular Clouds to Protostars

Stellar and Quasar Feedback in Concert: Effects on AGN Accretion, Obscuration, and Outflows

The formation of massive, quiescent galaxies at cosmic noon

The impact of stellar feedback on hot gas in galaxy haloes: the Sunyaev-Zel'dovich effect and soft X-ray emission

A Simple Phenomenological Model for Grain Clustering in Turbulence

Accurate, Meshless Methods for Magneto-Hydrodynamics

Do We Expect Most AGN to Live in Disks?

Forged in FIRE: cusps, cores, and baryons in low-mass dwarf galaxies

GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods

Gusty, gaseous flows of FIRE: galactic winds in cosmological simulations with explicit stellar feedback

Jumping the Gap: The Formation Conditions and Mass Function of Pebble-Pile Planetesimals

Sweating the small stuff: simulating dwarf galaxies, ultra-faint dwarf galaxies, and their own tiny satellites

The creation and persistence of a misaligned gas disc in a simulated early-type galaxy