Bondi or not Bondi: The impact of resolution on accretion and drag force modelling for Supermassive Black Holes
Monthly Notices of the Royal Astronomical Society Oxford University Press 478:1 (2018) 995-1016
Abstract:
Whilst in galaxy-size simulations, supermassive black holes (SMBH) are entirely handled by sub-grid algorithms, computational power now allows the accretion radius of such objects to be resolved in smaller scale simulations. In this paper, we investigate the impact of resolution on two commonly used SMBH sub-grid algorithms; the Bondi-Hoyle-Lyttleton (BHL) formula for accretion onto a point mass, and the related estimate of the drag force exerted onto a point mass by a gaseous medium. We find that when the accretion region around the black hole scales with resolution, and the BHL formula is evaluated using local mass-averaged quantities, the accretion algorithm smoothly transitions from the analytic BHL formula (at low resolution) to a supply limited accretion (SLA) scheme (at high resolution). However, when a similar procedure is employed to estimate the drag force it can lead to significant errors in its magnitude, and/or apply this force in the wrong direction in highly resolved simulations. At high Mach numbers and for small accretors, we also find evidence of the advective-acoustic instability operating in the adiabatic case, and of an instability developing around the wake's stagnation point in the quasi-isothermal case. Moreover, at very high resolution, and Mach numbers above $\mathcal{M}_\infty \geq 3$, the flow behind the accretion bow shock becomes entirely dominated by these instabilities. As a result, accretion rates onto the black hole drop by about an order of magnitude in the adiabatic case, compared to the analytic BHL formula.Gas flows in the circumgalactic medium around simulated high-redshift galaxies
Monthly Notices of the Royal Astronomical Society Oxford University Press 474:4 (2017) 4279-4301
Abstract:
We analyse the properties of circumgalactic gas around simulated galaxies in the redshift range z ≥ 3, utilizing a new sample of cosmological zoom simulations. These simulations are intended to be representative of the observed samples of Lyman α (Ly α) emitters recently obtained with the multi unit spectroscopic explorer (MUSE) instrument (halo masses ~ 10 10 - 10 11 M⊙). We show that supernova feedback has a significant impact on both the inflowing and outflowing circumgalactic medium (CGM) by driving outflows, reducing diffuse inflow rates, and by increasing the neutral fraction of inflowing gas. By temporally stacking simulation outputs, we find that significant net mass exchange occurs between inflowing and outflowing phases: none of the phases are mass-conserving. In particular, we find that the mass in neutral outflowing hydrogen declines exponentially with radius as gas flows outwards from the halo centre. This is likely caused by a combination of both fountain-like cycling processes and gradual photoionization/collisional ionization of outflowing gas. Our simulations do not predict the presence of fast-moving neutral outflows in the CGM. Neutral outflows instead move with modest radial velocities (~ 50 km s -1 ), and the majority of the kinetic energy is associated with tangential rather than radial motion.Cosmic evolution of stellar quenching by AGN feedback: clues from the Horizon-AGN simulation
Monthly Notices of the Royal Astronomical Society Oxford University Press 472:1 (2017) 949-965
Abstract:
The observed massive end of the local galaxy stellar mass function is steeper than its predicted dark matter (DM) halo counterpart in the standard $\Lambda $CDM paradigm. We investigate how active galactic nuclei (AGN) feedback can account for such a reduction in the stellar content of massive galaxies, through an influence on the gas content of their interstellar (ISM) and circum-galactic medium (CGM). We isolate the impact of AGNs by comparing two simulations from the HORIZON suite, which are identical except that one includes super massive black holes (SMBH) and related feedback. This allows us to cross-identify individual galaxies between these simulations and quantify the effect of AGN feedback on their properties, such as stellar mass and gas outflows. We find that the most massive galaxies ($ \rm M_{*} \geq 3 \times 10^{11} M_\odot $) are quenched to the extent that their stellar masses decrease by about 80% at $z=0$. More generally, SMBHs affect their host halo through a combination of outflows that reduce their baryonic mass, particularly for galaxies in the mass range $ \rm 10^9 M_\odot \leq M_{*} \leq 10^{11} M_\odot $, and a disruption of central gas inflows, which limits in-situ star formation, particularly massive galaxies with $ \rm M_{*} \approx10^{11} M_\odot $. As a result of these processes, net gas inflows onto massive galaxies drop by up to 70%. Finally, we measure a redshift evolution in the stellar mass ratio of twin galaxies with and without AGN feedback, with galaxies of a given stellar mass showing stronger signs of quenching earlier on. This evolution is driven by a progressive flattening of the $\rm M_{SMBH}-M_* $ relation for galaxies with $\rm M_{*} \leq 10^{10} M_\odot $ as redshift decreases, which translates into smaller SBMHs being harboured by galaxies of any fixed stellar mass, and indicates stronger AGN feedback at higher redshift.Fluctuating feedback-regulated escape fraction of ionizing radiation in low-mass, high-redshift galaxies
Monthly Notices of the Royal Astronomical Society Oxford University Press 470:1 (2017) 224-239
Abstract:
Low-mass galaxies are thought to provide the bulk of the ionizing radiation necessary to reionize the Universe. The amount of photons escaping the galaxies is poorly constrained theoretically, and difficult to measure observationally. Yet it is an essential parameter of reionization models.We study in detail how ionizing radiation can leak from high-redshift galaxies. For this purpose, we use a series of high-resolution radiation hydrodynamics simulations, zooming on three dwarf galaxies in a cosmological context. We find that the energy and momentum input from the supernova explosions has a pivotal role in regulating the escape fraction by disrupting dense star-forming clumps, and clearing sightlines in the halo. In the absence of supernovae, photons are absorbed very locally, within the birth clouds of massive stars. We follow the time evolution of the escape fraction and find that it can vary by more than six orders of magnitude. This explains the large scatter in the value of the escape fraction found by previous studies. This fast variability also impacts the observability of the sources of reionization: a survey even as deep as M 1500 = -14 would miss about half of the underlying population of Lyman-continuum emitters.Galaxy-halo alignments in the Horizon-AGN cosmological hydrodynamical simulation
Monthly Notices of the Royal Astronomical Society Oxford University Press (2017)