Person:

Snyder, Gregory

In extragalactic astronomy, a central challenge is that we cannot directly watch what happens to galaxies before and after they are observed. This dissertation focuses on linking predictions of galaxy time-evolution directly with observations, evaluating how interactions, mergers, and other processes affect the appearance of elliptical galaxies. The primary approach is to combine hydrodynamical simulations of galaxy formation, including all major components, with dust radiative transfer to predict their observational signatures. The current paradigm implies that a quiescent elliptical emerges following a formative starburst event. These trigger accretion onto the central supermassive black hole (SMBH), which then radiates as an active galactic nucleus (AGN). However, it is not clear the extent to which SMBH growth is fueled by these events nor how important is their energy input at setting the appearance of the remnant. This thesis presents results drawing from three phases in the formation of a typical elliptical: 1) I evaluate how to disentangle AGN from star formation signatures in mid-infrared spectra during a dust-enshrouded starburst, making testable predictions for robustly tracing SMBH growth with the James Webb Space Telescope ; 2) I develop a model for the rate of merger-induced post-starburst galaxies selected from optical spectra, resolving tension between their observed rarity and merger rates from other estimates; and 3) I present results from Hubble Space Telescope imaging of elliptical galaxies in galaxy clusters at 1 < z < 2, the precursors of present-day massive clusters with (M \sim10^{15}M_{\odot}), demonstrating that their stars formed over an extended period and ruling out the simplest model for their formation history. These results lend support to a stochastic formation history for ellipticals driven by mergers or interactions. However, significant uncertainties remain in how to evaluate the implications of galaxy appearance, in particular their morphologies across cosmic time. In the final chapter, I outline an approach to build a "mock observatory" from cosmological hydrodynamical simulations, with which observations of all types, including at high spatial and spectral resolutions, can be brought to bear in directly constraining the physics of galaxy formation and evolution.

We analyze the link between active galactic nuclei (AGNs) 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-4. This highly obscured burst begins star-formation-dominated with significant polycyclic aromatic hydrocarbon (PAH) emission, and ends with a ~109 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 PAHs, and enhances the 9.7 μm 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 μm, 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.

Supermassive black hole (BH) mergers produce powerful gravitational wave emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding their observability. By extracting merging BHs and host galaxy properties from the Illustris cosmological simulations, we have developed a comprehensive model for recoiling AGN. Here, for the first time, we model the effects of BH spin alignment and recoil dynamics based on the gas richness of host galaxies. We predict that if BH spins are not highly aligned, seeing-limited observations could resolve offset AGN, making them promising targets for all-sky surveys. For randomly oriented spins, ≲ 10 spatially offset AGN may be detectable in Hubble Space Telescope-Cosmological Evolution Survey, and >103 could be found with the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS), the Large Synoptic Survey Telescope (LSST), Euclid, and the Wide-Field Infrared Survey Telescope (WFIRST). Nearly a thousand velocity offset AGN are predicted within the Sloan Digital Sky Survey (SDSS) footprint; the rarity of large broad-line offsets among SDSS quasars is likely due in part to selection effects but suggests that spin alignment plays a role in suppressing recoils. None the less, in our most physically motivated model where alignment occurs only in gas-rich mergers, hundreds of offset AGN should be found in all-sky surveys. Our findings strongly motivate a dedicated search for recoiling AGN.