The dense environment of a galaxy cluster can radically transform the content of in-falling galaxies. Recent observations have found a significant population of active galactic nuclei (AGN) within "jellyfish galaxies," galaxies with trailing tails of gas and stars that indicate significant ram pressure stripping. The relationship between AGN and ram pressure stripping is not well understood. In this letter, we investigate the connection between AGN activity and ram pressure in a fully cosmological setting for the first time using the RomulusC simulation, one of the highest resolution simulations of a galaxy cluster to date. We find unambiguous morphological evidence for ram pressure stripping. For lower mass galaxies (with stellar masses < 10^9.5 solar masses) both star formation and black hole accretion are suppressed by ram pressure before they reach pericenter, whereas for more massive galaxies accretion onto the black hole is enhanced during pericentric passage. Our analysis also indicates that as long as the galaxy retains gas, AGN with higher Eddington ratios are more likely to be the found in galaxies experiencing higher ram pressure. We conclude that prior to quenching star formation, ram pressure triggers enhanced accretion onto the black hole, which then produces heating and outflows due to AGN feedback. AGN feedback may in turn serve to aid in the quenching of star formation in tandem with ram pressure.

We study the link between supermassive black hole growth and the stellar mass assembly of their host galaxies in the state-of-the-art Romulus suite of simulations. The cosmological simulations Romulus25 and RomulusC employ innovative recipes for the seeding, accretion, and dynamics of black holes in the field and cluster environments respectively. We find that the black hole accretion rate traces the star formation rate among star-forming galaxies. This result holds for stellar masses between 10^8 and 10^12 solar masses, with a very weak dependence on host halo mass or redshift. The inferred relation between accretion rate and star formation rate does not appear to depend on environment, as no difference is seen in the cluster/proto-cluster volume compared to the field. A model including the star formation rate, the black hole-to-stellar mass ratio, and the cold gas fraction can explain about 70 per cent of all variations in the black hole accretion rate among star forming galaxies. Finally, bearing in mind the limited volume and resolution of these cosmological simulations, we find no evidence for a connection between black hole growth and galaxy mergers, on any timescale and at any redshift. Black holes and their galaxies assemble in tandem in these simulations, regardless of the larger-scale intergalactic environment, suggesting that black hole growth simply follows star formation on galactic scales.

There exist hitherto unexplained fluctuations in the Cosmic Infrared Background (CIB) on arcminute scales and larger. These have been shown to cross-correlate with the Cosmic X-ray Background (CXB), leading several authors to attribute the excess to a high-redshift growing black hole population. In order to investigate potential sources that could explain this excess, in this paper, we develop a new framework to compute the power spectrum of undetected sources that do not have constant flux as a function of halo mass. In this formulation, we combine a semi-analytic model for black hole growth and their simulated spectra from hydrodynamical simulations. Revisiting the possible contribution of a high-redshift black hole population, we find that too much black hole growth is required at early epochs for z>6 accretion to explain these fluctuations. Examining a population of accreting black holes at more moderate redshifts, z\sim 2-3, we find that such models produce a poor fit to the observed fluctuations while simultaneously overproducing the local black hole mass density. Additionally, we rule out the hypothesis of a missing Galactic foreground of warm dust that produces coherent fluctuations in the X-ray via reflection of Galactic X-ray binary emission. Although we firmly rule out accreting massive black holes as the source of these missing fluctuations, additional studies will be required to determine their origin.

The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains challenging. Several new observables that encapsulate information about seeding have been proposed in recent years, and these offer exciting prospects for truly unraveling the nature of black hole seeds in the coming years. One of the key challenges for this task lies in the complexity of the problem, the required disentangling of the confounding effects of accretion physics and mergers, as mergers and accretion events over cosmic time stand to erase these initial conditions. Nevertheless, some unique signatures of seeding do survive and still exist in: local scaling relations between black holes and their galaxy hosts at low-masses; in high-redshift luminosity functions of accreting black holes; and in the total number and mass functions of gravitational wave coalescence events from merging binary black holes. One of the clearest discriminants for seed models are these high redshift gravitational wave detections of mergers from space detectable in the milliHertz range. These predicted event rates offer the most direct constraints on the properties of initial black hole seeds. Improving our theoretical understanding of black hole dynamics and accretion will also be pivotal in constraining seeding models in combination with the wide range of multi-messenger data.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that RomulusC results in an intracluster medium (ICM) consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster center. This enhanced quenched population extends beyond $R_{200}$ and is in place at high redshift. Similarly, we predict that SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest AGN in the most massive galaxies.

The origin and properties of the initial black hole seeds that grow to produce the observed population of accreting sources remain to be determined. It is a challenge to uniquely disentangle signatures of seeding from fueling and dynamics that shapes the assembly history of growing black holes. To address this, we use a semi-analytic model developed to track the growth of supermassive black holes adopting multiple prescriptions for accretion. In contrast with earlier treatments, we explore the interplay between seeding models and two accretion modes. We find that signatures of the initial seeding do survive in the following observational probes: the black hole occupation fraction; contribution to the unresolved X-ray background; low-luminosity and high-redshift luminosity functions; and in gravitational wave event signatures. We find that the behaviour of the low-mass end of the M-sigma relation is dominated by uncertainties in the adopted accretion prescriptions and does not offer clear discrimination between seeding models. We make concrete predictions for future surveys, particularly for the Lynx X-ray surveyor and LISA (The Laser Interferometer Space Antenna) mission, which will each provide different and yet strong constraints on the seed population. Black hole coalescences detected by LISA and high-redshift quasar luminosity functions observed by Lynx will offer the sharpest seeding discriminants. Although the signatures of the black hole seeding mechanism that persist remain linked to our understanding of black hole accretion and dynamics, we offer new insights on how these upcoming multi-wavelength observations could be leveraged to effectively disentangle them.

We develop a semi-analytic model to explore different prescriptions of supermassive black hole (SMBH) fuelling. This model utilises a merger-triggered burst mode in concert with two possible implementations of a long-lived steady mode for assembling the mass of the black hole in a galactic nucleus. We improve modelling of the galaxy-halo connection in order to more realistically determine the evolution of a halo's velocity dispersion. We use four model variants to explore a suite of observables: the M-sigma relation, mass functions of both the overall and broad-line quasar population, and luminosity functions as a function of redshift. We find that "downsizing" is a natural consequence of our improved velocity dispersion mappings, and that high-mass SMBHs assemble earlier than low-mass SMBHs. The burst mode of fuelling is sufficient to explain the assembly of SMBHs to z=2, but an additional steady mode is required to both assemble low-mass SMBHs and reproduce the low-redshift luminosity function. We discuss in detail the trade-offs in matching various observables and the interconnected modelling components that govern them. As a result, we demonstrate the utility as well as the limitations of these semi-analytic techniques.

Growing the supermassive black holes (~10^9 Msun) that power the detected luminous, highest redshift quasars (z > 6) from light seeds - the remnants of the first stars - within ~ 1 Gyr of the Big Bang poses a timing challenge for growth models. The formation of massive black hole seeds via direct collapse with initial masses ~ 10^4 - 10^5 Msun alleviates this problem. Physical conditions required to form these massive direct collapse black hole (DCBH) seeds are available in the early universe. These viable DCBH formation sites, satellite halos of star-forming galaxies, merge and acquire a stellar component. These produce a new, transient class of objects at high redshift, Obese Black hole Galaxies (OBGs), where the luminosity produced by accretion onto the black hole outshines the stellar component. Therefore, the OBG stage offers a unique way to discriminate between light and massive initial seeds. We predict the multi-wavelength energy output of OBGs and growing Pop III remnants at a fiducial redshift (z = 9), exploring both standard and slim disk accretion onto the growing central black hole for high and low metallicities of the associated stellar population. With our computed templates, we derive the selection criteria for OBGs, that comprise a pre-selection that eliminates blue sources; followed by color-color cuts ([F_{070W} - F_{220W}] > 0; -0.3 < [F_{200W} - F_{444W}] < 0.3) and when available, a high ratio of X-ray flux to rest-frame optical flux (F_X/F_{444W} >> 1) (Abridged).

Using {\em Chandra} observations in the 2.15 deg$^{2}$ COSMOS legacy field, we present one of the most accurate measurements of the Cosmic X-ray Background (CXB) spectrum to date in the [0.3-7] keV energy band. The CXB has three distinct components: contributions from two Galactic collisional thermal plasmas at kT$\sim$0.27 and 0.07 keV and an extragalactic power-law with photon spectral index $\Gamma$=1.45$\pm{0.02}$. The 1 keV normalization of the extragalactic component is 10.91$\pm{0.16}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. Removing all X-ray detected sources, the remaining unresolved CXB is best-fit by a power-law with normalization 4.18$\pm{0.26}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$ and photon spectral index $\Gamma$=1.57$\pm{0.10}$. Removing faint galaxies down to i$_{AB}\sim$27-28 leaves a hard spectrum with $\Gamma\sim$1.25 and a 1 keV normalization of $\sim$1.37 keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. This means that $\sim$91\% of the observed CXB is resolved into detected X-ray sources and undetected galaxies. Unresolved sources that contribute $\sim 8-9\%$ of the total CXB show a marginal evidence of being harder and possibly more obscured than resolved sources. Another $\sim$1\% of the CXB can be attributed to still undetected star forming galaxies and absorbed AGN. According to these limits, we investigate a scenario where early black holes totally account for non source CXB fraction and constrain some of their properties. In order to not exceed the remaining CXB and the $z\sim$6 accreted mass density, such a population of black holes must grow in Compton-thick envelopes with N$_{H}>$1.6$\times$10$^{25}$ cm$^{-2}$ and form in extremely low metallicity environments $(Z_\odot)\sim10^{-3}$.

Some reddened quasars appear to be transitional objects in the merger-induced black hole growth/galaxy evolution paradigm, where a heavily obscured nucleus starts to be unveiled by powerful quasar winds evacuating the surrounding cocoon of dust and gas. Hard X-ray observations are able to peer through this gas and dust, revealing the properties of circumnuclear obscuration. Here, we present NuSTAR and XMM-Newton/Chandra observations of FIRST-2MASS selected red quasars F2M 0830+3759 and F2M 1227+3214. We find that though F2M 0830+3759 is moderately obscured ($N_{\rm H,Z} = 2.1\pm0.2 \times10^{22}$ cm$^{-2}$) and F2M 1227+3214 is mildly absorbed ($N_{\rm H,Z} = 3.4^{+0.8}_{-0.7}\times10^{21}$ cm$^{-2}$) along the line-of-sight, heavier global obscuration may be present in both sources, with $N_{\rm H,S} = 3.7^{+4.1}_{-2.6} \times 10^{23}$ cm$^{-2}$ and $< 5.5\times10^{23}$ cm$^{-2}$, for F2M 0830+3759 and F2M 1227+3214, respectively. F2M 0830+3759 also has an excess of soft X-ray emission below 1 keV which is well accommodated by a model where 7% of the intrinsic AGN X-ray emission is scattered into the line-of-sight. While F2M 1227+3214 has a dust-to-gas ratio ($E(B-V)$/$N_{\rm H}$) consistent with the Galactic value, the $E(B-V)$/$N_{\rm H}$ value for F2M 0830+3759 is lower than the Galactic standard, consistent with the paradigm that the dust resides on galactic scales while the X-ray reprocessing gas originates within the dust-sublimation zone of the broad-line-region. The X-ray and 6.1$\mu$m luminosities of these red quasars are consistent with the empirical relations derived for high-luminosity, unobscured quasars, extending the parameter space of obscured AGN previously observed by NuSTAR to higher luminosities.

The presence of debris disks around young main sequence stars hints at the existence and structure of planetary systems. Millimeter-wavelength observations probe large grains that trace the location of planetesimal belts. The FEPS (Formation and Evolution of Planetary Systems) $Spitzer$ Legacy survey of nearby young solar analogues yielded a sample of five debris disk-hosting stars with millimeter flux suitable for interferometric follow-up. We present observations with the Submillimeter Array (SMA) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at ~2" resolution that spatially resolve the debris disks around these nearby ($d\sim$50 pc) stars. Two of the five disks (HD 377, HD 8907) are spatially resolved for the first time and one (HD 104860) is resolved at millimeter wavelengths for the first time. We combine our new observations with archival SMA and Atacama Large Millimeter/Submillimeter Array (ALMA) data to enable a uniform analysis of the full five-object sample. We simultaneously model the broad-band photometric data and resolved millimeter visibilities to constrain the dust temperatures and disk morphologies, and perform an MCMC analysis to fit for basic structural parameters. We find that the radii and widths of the cold outer belts exhibit properties consistent with scaled-up versions of the Solar System's Kuiper Belt. All the disks exhibit characteristic grain sizes comparable to the blowout size, and all the resolved observations of emission from large dust grains are consistent with an axisymmetric dust distribution to within the uncertainties. These results are consistent with comparable studies carried out at infrared wavelengths.

Massive black hole binaries (MBHBs) are a natural byproduct of galaxy mergers. Previous studies have shown that flares from stellar tidal disruption events (TDEs) are modified by the presence of a secondary perturber, causing interruptions in the light curve. We study the dynamics of TDE debris in the presence of a milliparsec-separated MBHB by integrating ballistic particle orbits in the time-varying potential of the binary. We find that gaps in the light curve appear when material misses the accretion radius on its first return to pericentre. Subsequent recurrences can be decomposed into "continuous" and "delayed" components, which exhibit different behaviour. We find that this potential can substantially alter the locations of stream self-intersections. When debris is confined to the plane, we find that close encounters with the secondary BH leave noticeable signatures on the fallback rate and can result in significant accretion onto the secondary BH. Tight, equal-mass MBHBs accrete equally, periodically trading the infalling stream.

The Event Horizon Telescope (EHT), a global sub-millimeter wavelength very long baseline interferometry array, is now resolving the innermost regions around the supermassive black holes Sgr A* and M87. Using black hole images from both simple geometric models and relativistic magnetohydrodynamical accretion flow simulations, we perform a variety of experiments to assess the promise of the EHT for studying strong gravity and accretion physics during the stages of its development. We find that (1) the addition of the LMT and ALMA along with upgraded instrumentation in the "Complete" stage of the EHT allow detection of the photon ring, a signature of Kerr strong gravity, for predicted values of its total flux; (2) the inclusion of coherently averaged closure phases in our analysis dramatically improves the precision of even the current array, allowing (3) significantly tighter constraints on plausible accretion models and (4) detections of structural variability at the levels predicted by the models. While observations at 345 GHz circumvent problems due to interstellar electron scattering in line-of-sight to the galactic center, short baselines provided by CARMA and/or the LMT could be required in order to constrain the overall shape of the accretion flow. Given the systematic uncertainties in the underlying models, using the full complement of two observing frequencies (230 and 345 GHz) and sources (Sgr A* and M87) may be critical for achieving transformative science with the EHT experiment.

HD 61005, also known as "The Moth," is one of only a handful of debris disks that exhibit swept-back "wings" thought to be caused by interaction with the ambient interstellar medium (ISM). We present 1.3 mm Submillimeter Array (SMA) observations of the debris disk around HD 61005 at a spatial resolution of 1.9 arcsec that resolve the emission from large grains for the first time. The disk exhibits a double-peaked morphology at millimeter wavelengths, consistent with an optically thin ring viewed close to edge-on. To investigate the disk structure and the properties of the dust grains we simultaneously model the spatially resolved 1.3 mm visibilities and the unresolved spectral energy distribution. The temperatures indicated by the SED are consistent with expected temperatures for grains close to the blowout size located at radii commensurate with the millimeter and scattered light data. We also perform a visibility-domain analysis of the spatial distribution of millimeter-wavelength flux, incorporating constraints on the disk geometry from scattered light imaging, and find suggestive evidence of wavelength-dependent structure. The millimeter-wavelength emission apparently originates predominantly from the thin ring component rather than tracing the "wings" observed in scattered light. The implied segregation of large dust grains in the ring is consistent with an ISM-driven origin for the scattered light wings.